This work presents measurements of liquid drop deformation and breakup time behind approximately conical shock waves and evaluates the predictive capabilities of low-order models and correlations developed using planar shock experiments. A conical shock was approximated by firing a bullet at Mach 4.5 past a vertical column of water drops with a mean initial diameter of 192 µm. The time-resolved drop position and maximum transverse dimension were characterized using backlit stereo images taken at 500 kHz. The gas density and velocity fields experienced by the drops were estimated using a Reynolds-averaged Navier-Stokes simulation of the bullet. Classical correlations predict drop breakup times and deformation in error by a factor of 3 or more. The Taylor analogy breakup (TAB) model predicts deformed drop diameters that agree within the confidence bounds of the ensemble-averaged experimental values using a dimensionless constant C2 = 2 compared to the accepted value C2 = 2/3. Results demonstrate existing correlations are inadequate for predicting the drop response to the three-dimensional relaxation of the flowfield downstream of a conical-like shock and suggest the TAB model results represent a path toward improved predictions.
Villa, Daniel L.; Schostek, Tyler; Govertsen, Krissy; Macmillan, Madeline
Applying extreme temperature events for future conditions is not straightforward for infrastructure resilience analyses. This work introduces a stochastic model that fills this gap. The model uses at least 50 years of daily extreme temperature records, climate normals with 10%–90% confidence intervals, and shifts/offsets for increased frequency and intensity of heat wave events. Intensity and frequency are shifted based on surface temperature anomaly from 1850–1900 for 32 models from CMIP6. A case study for Worcester, Massachusetts passed 85% of cases using the two-sided Kolmogorov–Smirnov p-value test with 95% confidence for both temperature and duration. Future shifts for several climate scenarios to 2020, 2040, 2060, and 2080 had acceptable errors between the shifted model and 10- and 50-year extreme temperature event thresholds with the largest error being 2.67°C. The model is likely to be flexible enough for other patterns of extreme weather such as extreme precipitation and hurricanes.
The purpose of this protocol is to define procedures and practices to be used by the PACT center for field testing of metal halide perovskite (MHP) photovoltaic (PV) modules. The protocol defines the physical, electrical, and analytical configuration of the tests and applies equally to mounting systems at a fixed orientation or sun tracking systems. While standards exist for outdoor testing of conventional PV modules, these do not anticipate the unique electrical behavior of perovskite cells. Further, the existing standards are oriented toward mature, relatively stable products with lifetimes that can be measured on the scale of years to decades. The state of the art for MHP modules is still immature with considerable sample to sample variation among nominally identical modules. Version 0.0 of this protocol does not define a minimum test duration, although the intent is for modules to be fielded for periods ranging for weeks to months. This protocol draws from relevant parts of existing standards, and where necessary includes modifications specific to the behavior of perovskites.
The Sandia National Laboratories, in California (Sandia/CA) is a research and development facility, owned by the U.S. Department of Energy’s National Nuclear Security Administration agency (DOE/NNSA). The laboratory is located in the City of Livermore (the City) and is comprised of approximately 410 acres. The Sandia/CA facility is operated by National Technology and Engineering Solutions of Sandia, LLC (NTESS) under a contract with the DOE/NNSA. The DOE/ NNSA’s Sandia Field Office (SFO) oversees the operations of the site. North of the Sandia/CA facility is the Lawrence Livermore National Laboratory (LLNL), in which Sandia/CA’s sewer system combines with before discharging to the City’s Publicly Owned Treatment Works (POTW) for final treatment and processing. The City’s POTW authorizes the wastewater discharge from Sandia/CA via the assigned Wastewater Discharge Permit #1251 (the Permit), which is issued to the DOE/NNSA’s main office for Sandia National Laboratories, located in New Mexico (Sandia/NM). The Permit requires the submittal of this Monthly Sewer Monitoring Report to the City by the twenty-fifth day of each month.
This document provides the instructions for participating in the 2021 blind photovoltaic (PV) modeling intercomparison organized by the PV Performance Modeling Collaborative (PVPMC). It describes the system configurations, metadata, and other information necessary for the modeling exercise. The practical details of the validation datasets are also described. The datasets were published online in open access in April 2023, after completing the analysis of the results.
Diesel piston-bowl shape is a key design parameter that affects spray-wall interactions and turbulent flow development, and in turn affects the engine’s thermal efficiency and emissions. It is hypothesized that thermal efficiency can be improved by enhancing squish-region vortices as they are hypothesized to promote fuel-air mixing, leading to faster heat-release rates. However, the strength and longevity of these vortices decrease with advanced injection timings for typical stepped-lip (SL) piston geometries. Dimple stepped-lip (DSL) pistons enhance vortex formation at early injection timings. Previous engine experiments with such a bowl show 1.4% thermal efficiency gains over an SL piston. However, soot was increased dramatically [SAE 2022-01-0400]. In a previous study, a new DSL bowl was designed using non-combusting computational fluid dynamic simulations. This improved DSL bowl is predicted to promote stronger, more rotationally energetic vortices than the baseline DSL piston: it employs shallower, narrower, and steeper-curved dimples that are placed further out into the squish region. In the current experimental study, this improved bowl is tested in a medium-duty diesel engine and compared against the SL piston over an injection timing sweep at low-load and part-load operating conditions. No substantial thermal efficiency gains are achieved at the early injection timing with the improved DSL design, but soot emissions are lowered by 45% relative to the production SL piston, likely due to improved air utilization and soot oxidation. However, these benefits are lost at late injection timings, where the DSL piston renders a lower thermal efficiency than that of the SL piston. Energy balance analyses show higher wall heat transfer with the DSL piston than with the SL piston despite a 1.3% reduction in the piston surface area. Vortex enhancement may not necessarily lead to improved efficiency as more energetic squish-region vortices can lead to higher convective heat transfer losses.
Solving large number of small linear systems is increasingly becoming a bottleneck in computational science applications. While dense linear solvers for such systems have been studied before, batched sparse linear solvers are just starting to emerge. In this paper, we discuss algorithms for solving batched sparse linear systems and their implementation in the Kokkos Kernels library. The new algorithms are performance portable and map well to the hierarchical parallelism available in modern accelerator architectures. The sparse matrix vector product (SPMV) kernel is the main performance bottleneck of the Krylov solvers we implement in this work. The implementation of the batched SPMV and its performance are therefore discussed thoroughly in this paper. The implemented kernels are tested on different Central Processing Unit (CPU) and Graphic Processing Unit (GPU) architectures. We also develop batched Conjugate Gradient (CG) and batched Generalized Minimum Residual (GMRES) solvers using the batched SPMV. Our proposed solver was able to solve 20,000 sparse linear systems on V100 GPUs with a mean speedup of 76x and 924x compared to using a parallel sparse solver with a block diagonal system with all the small linear systems, and compared to solving the small systems one at a time, respectively. We see mean speedup of 0.51 compared to dense batched solver of cuSOLVER on V100, while using lot less memory. Thorough performance evaluation on three different architectures and analysis of the performance are presented.
In magneto-inertial fusion, the ratio of the characteristic fuel length perpendicular to the applied magnetic field R to the α-particle Larmor radius Q α is a critical parameter setting the scale of electron thermal-conduction loss and charged burn-product confinement. Using a previously developed deep-learning-based Bayesian inference tool, we obtain the magnetic-field fuel-radius product B R ∝ R / Q α from an ensemble of 16 magnetized liner inertial fusion (MagLIF) experiments. Observations of the trends in BR are consistent with relative trade-offs between compression and flux loss as well as the impact of mix from 1D resistive radiation magneto-hydrodynamics simulations in all but two experiments, for which 3D effects are hypothesized to play a significant role. Finally, we explain the relationship between BR and the generalized Lawson parameter χ. Our results indicate the ability to improve performance in MagLIF through careful tuning of experimental inputs, while also highlighting key risks from mix and 3D effects that must be mitigated in scaling MagLIF to higher currents with a next-generation driver.
Sandia National Laboratories (SNL) has developed a novel reduced order modeling approach. Prioritization of inputs is accomplished using Sobo' indices obtained through a more efficient variance-based global sensitivity analysis. To determine the Sobo' functions, simulated input values are aligned to collocation points to permit the use of Gauss-Lobatto integration, thereby reducing the number of simulation trials needed by more than an order of magnitude compared to standard Monte Carlo approaches. Furthermore, by leveraging the orthogonality of Legendre polynomials in conjunction with those same simulations at the collocation nodes, an efficient fitting method is developed to represent the Sobo' functions from which a reduced order model (ROM) is constructed. The developed method is both more efficient computationally, and the resulting ROM is more accurate. The efficacy of this technique is demonstrated on a nonlinear polynomial test function as well as the nonlinear Ishigami and Sobo' g functions.
Induced seismicity is an inherent risk associated with geologic carbon storage (GCS) in deep rock formations that could contain undetected faults prone to failure. Modeling-based risk assessment has been implemented to quantify the potential of injection-induced seismicity, but typically simplified multiscale geologic features or neglected multiphysics coupled mechanisms because of the uncertainty in field data and computational cost of field-scale simulations, which may limit the reliable prediction of seismic hazard caused by industrial-scale CO2 storage. The degree of lateral continuity of the stratigraphic interbedding below the reservoir and depth-dependent fault permeability can enhance or inhibit pore-pressure diffusion and corresponding poroelastic stressing along a basement fault. This study presents a rigorous modeling scheme with optimal geological and operational parameters needed to be considered in seismic monitoring and mitigation strategies for safe GCS.
Deploying Tidal Energy Converters for electricity generation requires prior-knowledge of the potential Annual Energy Production (AEP) at the site, Ideally using a year-long tidal current record at the proposed site to minimize uncertainty. However, such records are often unavailable. Fortunately, using the periodic nature of tidal variability, the International Electrotechnical Commission Technical Specification for tidal energy resource assessment requires AEP calculation using at least 90 days of tidal current records at each turbine location. The sensitivity of AEP to different record durations has not been fully assessed. This is the goal of our study. The study utilized the U.S. tidal energy geodatabase to simulate tidal currents with various lengths, during 100 years of the 21st century. We then consider two frameworks for evaluating AEP: (a) The long-term (months) fixed instrument (FI) measurement at each proposed tidal turbine location, and (b) one FI measurement and short-term (hours) boat-based moving vessel measurements. Under the two scenarios, we examine the AEP assessed from short tidal current records, including how the AEP uncertainties vary spatially and temporally, and how they are associated with various astronomical factors. This helps provide guidance on choosing the appropriate assessment methodologies to reduce the AEP uncertainties and project cost.
We propose primal–dual mesh optimization algorithms that overcome shortcomings of the standard algorithm while retaining some of its desirable features. “Hodge-Optimized Triangulations” defines the “HOT energy” as a bound on the discretization error of the diagonalized Delaunay Hodge star operator. HOT energy is a natural choice for an objective function, but unstable for both mathematical and algorithmic reasons: it has minima for collapsed edges, and its extrapolation to non-regular triangulations is inaccurate and has unbounded minima. We propose a different extrapolation with a stronger theoretical foundation, and avoid extrapolation by recalculating the objective just beyond the flip threshold. We propose new objectives, based on normalizations of the HOT energy, with barriers to edge collapses and other undesirable configurations. We propose mesh improvement algorithms coupling these. When HOT optimization nearly collapses an edge, we actually collapse the edge. Otherwise, we use the barrier objective to update positions and weights and remove vertices. By combining discrete connectivity changes with continuous optimization, we more fully explore the space of possible meshes and obtain higher quality solutions.
HyRAM+ is a toolkit that includes fast-running models for the unconstrained (i.e., no wall interactions) dispersion and flames for non-premixed fuels. The models were developed for use with hydrogen, but the toolkit was expanded to include propane and methane in a recent release. In this work we validate the dispersion and flame models for these additional fuels, based on reported literature data. The validation efforts spanned a range of release conditions, from subsonic to underexpanded jets and flames for a range of mass flow rates. In general, the dispersion model works well for both propane and methane although the width of the jet/plume is predicted to be wider than observed in some cases. The flame model tends to over-predict the induced buoyancy for low-momentum flames, while the radiative heat flux agrees with the experimental data reasonably well, for both fuels. The models could be improved but give acceptable predictions for propane and methane behavior for the purposes of risk assessment.