Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Refining”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

Asynchronous-many-task systems: Challenges and opportunities - Scaling an AMR astrophysics code on exascale machines using Kokkos and HPX

Dynamic and adaptive mesh refinement is pivotal in high-resolution, multi-physics, multi-model simulations, necessitating precise physics resolution in localized areas across expansive domains. Today’s supercomputers’ extreme heterogeneity presents a significant challenge for dynamically adaptive codes, highlighting the importance of achieving performance portability at scale. Our research focuses on astrophysical simulations, particularly stellar mergers, to elucidate early universe dynamics. Here, we present Octo-Tiger, leveraging Kokkos, HPX, and SIMD for portable performance at scale in complex, massively parallel adaptive multi-physics simulations. Octo-Tiger supports diverse processors, accelerators, and network backends. Experiments demonstrate exceptional scalability across several heterogeneous supercomputers including Perlmutter, Frontier, and Fugaku, encompassing major GPU architectures and x86, ARM, and RISC-V CPUs. Parallel efficiency of 47.59% (110,080 cores and 6880 hybrid A100 GPUs) on a full-system run on Perlmutter (26% HPCG peak performance) and 51.37% (using 32,768 cores and 2048 MI250X) on Frontier are achieved.

97 MATHEMATICS AND COMPUTING↗

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗

mpfit: a robust method for fitting atomic resolution images with multiple Gaussian peaks

The standard technique for sub-pixel estimation of atom positions from atomic resolution scanning transmission electron microscopy images relies on fitting intensity maxima or minima with a two-dimensional Gaussian function. While this is a widespread method of measurement, it can be error prone in images with non-zero aberrations, strong intensity differences between adjacent atoms or in situations where the neighboring atom positions approach the resolution limit of the microscope. Here we demonstrate mpfit, an atom finding algorithm that iteratively calculates a series of overlapping two-dimensional Gaussian functions to fit the experimental dataset and then subsequently uses a subset of the calculated Gaussian functions to perform sub-pixel refinement of atom positions. Based on both simulated and experimental datasets presented in this work, this approach gives lower errors when compared to the commonly used single Gaussian peak fitting approach and demonstrates increased robustness over a wider range of experimental conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗

Compositional Reasoning for Hierarchical State Machines

Harel statecharts and its derivatives are popular graphical languages for specifying discrete control systems via hierarchical state machines. Separately, there has been a long line of work on specifying concurrent systems with process calculi which come equipped with an algebraic theory, the ability reason compositionally about various temporal properties, and strong type systems. While these two approaches to modeling systems are tantalizingly similar, the integrated reasoning principles that exist for process calculi have not been demonstrated in hierarchical state machines. A key issue is that operational theories for process calculi do not behave like control systems, and thus, there is virtually no tool support for modeling control systems with such languages. For a control system designer, bringing the integrated, more scalable reasoning from the process calculi to state-machine languages would enable the specification of more complex systems and a more modular systems development process. Our insight is that we can recover many important results from the process calculi in hierarchical state machines with local scope. We employ a structural operational semantics, which is ubiquitous in process and 𝜆-calculi but uncommon in hierarchical statemachine formalizations, to enable inductive reasoning about behavior. Taking inspiration from the structure of process calculi metatheories, we define a calculus of refinement and equivalence that we prove sound with respect to local notion of (bi)simulation. Furthermore, we prove that the calculus preserves the behavioral properties of reactivity, observational determinism, traces, and linear temporal properties. Our results are mechanized in the Rocq proof assistant.

97 MATHEMATICS AND COMPUTING↗

Applying Time-Parallelization to Turbulent Flows

Parallelization of the temporal domain is explored for the solution of turbulent flows. Multigrid reduction-in-time (MGRIT) is used to advance the large-scale fluid dynamics in time sequentially on the coarsest space-time grid but propagate the information in time parallel on all other levels. The goal of this process is to accurately and efficiently resolve the coarse-scale turbulence structure and use that to drive the fine-scales of the turbulent flow. The extra forcing from nonlinear multigrid facilitates the coupling and interaction between fine and coarse scales, through which the multiscale nonlinear physics is properly captured. Adaptive mesh refinement is employed to finely resolve only the regions with strong gradients, which provides further computational efficiency. The underlying computational fluid dynamics solver is a fourth-order finite-volume scheme with the standard 4-stage Runge-Kutta method. An advanced approach is devised and implemented to enable MGRIT to solve highly turbulent flows successfully. Furthermore, the method is applied to solve a Taylor-Green vortex problem and a doubleshear-layer turbulent mixing flow. Results are promising, validating that MGRIT with the filtering approach has the potential to efficiently solve general turbulent flows.

Computational Fluid Dynamics↗

High-Fidelity Simulations of Plasma-Assisted Oxidation of Hydrocarbon Fuels Using Nanosecond Pulsed Discharges

Next, a robust and efficient framework for simulating NSPD in multiple dimensions is developed. The reactive Navier-Stokes equations are extended to include a drift-diffusion plasma-fluid model with a local field approximation (LFA) in a finite-volume solver, which uses an adaptive mesh refinement (AMR) strategy to address the wide separation of length scales in the problem. A two-way coupling strategy is used whereby the plasma-fluid model and reactive Navier-Stokes equations are integrated simultaneously. The oxidation of ethylene/air mixtures mediated by NSPD is simulated in a pin-to-pin configuration. All phases of the plasma discharge are simulated explicitly (including streamer ignition, propagation, and connection, as well as the subsequent spark phase), along with the evolution of the plasma during the inter-pulse period. Temporally and spatially-resolved results are presented, with an emphasis on the analysis of heating and energy deposition, as well as of the evolution of the concentration of active particles generated during the NSPD and their influence on ignition.

adaptive mesh refinement↗

DNS of Flame Stabilization Dynamics of a Swirl-Stabilized Spray Burner Using Sustainable Aviation Fuels

Concern for emission reduction has motivated the development of new cost-effective alternative sustainable aviation fuels (SAFs). Drop-in SAFs with blends of the certified and the alternative fuels are beneficial because they do not require engine modifications for use in current aviation engines. Combustion characteristics are of particular concern when comparing a new fuel with Jet-A. Lean blow-off and flame stabilization have been identified to be governed by the fuel propensity to autoignition, i.e. derived cetane number. However, fundamental investigation is required to identify the impact of derived cetane number on the combustion mode of flame stabilization in a realistic combustor. The combustion mode has also a direct impact on turbulent combustion modeling closures. DNS simulations are performed in the low-Mach solver of the Pele Suite called PeleLMeX. Lagrangian multi-phase modeling is used to capture the liquid spray injection of Jet-A (reference fuel) and C1 as a representative of a low cetane number SAF. Adaptive Mesh Refinement (AMR) is used to enable a more efficient simulation of a more realistic domain size and embedded boundary treatment is used to model a bluff-body geometry. Local extinction and edge flame propagation were observed for both fuels. The edge flame propagation mode was quantified in terms of a Damkohler number defined as the ratio between progress variable reaction rate and its diffusive flux. Initial analysis suggests that a mixed mode combustion occurs for the edge flame propagation, with flame propagation assisted by ignition and autoignition co-existing for both Jet-A and C1 flames. The analysis of the extinction region shows a larger progress variable for Jet-A which can further increase the local displacement speed and present a faster reignition of the stoichiometric mixture fraction.

adaptive mesh refinement↗

Effect of Feedstock Variability, Feedstock Blends, and Pretreatment Conditions on Sugar Yield and Production Costs

Mass production of sustainable biofuels replacing fossil fuels will significantly reduce global greenhouse gas emissions. However, biomass feedstock variability, availability, and affordability pose significant risks for the sustainable production of biofuels in biorefinery plants. To derisk the sustainable production of biofuels, this work investigates the sugar yields and production costs of four different feedstocks as well as their blends, using different pretreatment methods: deacetylation and dilute acid pretreatment (DDA) and deacetylation and mechanical refining pretreatment (DMR). We found that not only feedstock species variability, but also feedstock harvesting techniques have large impacts on both sugar yields and production costs. Single-pass harvested corn stover (SPCS) showed higher sugar yields and lower production costs than multipass harvested corn stover, followed by sorghum and switchgrass. The highest glucose yields of 91% and 95.3%, which corresponds to the lowest sugar production costs of $0.2286/lb. and $0.2490/lb., were observed for SPCS using DDA and DMR, respectively. The difference in sugar yields and production costs due to feedstock variability can be largely mitigated by increasing deacetylation severity, especially in DMR. Higher deacetylation severity resulted in decreased sugar production costs in the case of DMR in contrast to increased costs in the case of DDA. Blending switchgrass and sorghum into corn stover, which can significantly increase the feedstock availability and sustainability, resulted in sugar yields and production costs of roughly the weighted average of those derived from constituent individual feedstocks, respectively. Thus, blending feedstocks of different biomass species is a potential successful strategy for sustainable production of biofuels to replace fossil-derived fuels.

09 BIOMASS FUELS↗

A Possible Radiation-Induced Transition from Monazite-(Ce) to Xenotime-(Y)

This study examines two pegmatitic monazite samples (2a and 4b, these numbers are related to a previous study) to determine their crystal chemistry and effects of internal radiation damage using synchrotron high-resolution powder X-ray diffraction and electron-probe micro-analysis. Both the huttonite and cheralite substitutions are discussed. Rietveld structure refinement of sample 2a shows three different phases [2a = monazite-(Ce), 2b = monazite-(Ce), and 2c = xenotime-(Y)] with distinct structural parameters. The changes among the unit-cell parameters between the two monazite-(Ce) phases is more pronounced in the a followed by the b and c unit-cell parameters. Sample 4a is a single-phase monazite-(Sm) that contains 0.164 apfu Th. Phase 2c with space group I41/amd arises from redistribution of La, Ce, Pr, Nd, Sm, Gd, Dy, Si, and Y atoms from those in monazite (space group P21/n). A possible cause for the phase transition from monazite-(Ce) to xenotime-(Y) is α-radiation events over a long geological time. However, other chemical processes cannot be ruled out as a cause for the transition.

structural variations↗

Crystal Structure of Moganite and Its Anisotropic Atomic Displacement Parameters Determined by Synchrotron X-ray Diffraction and X-ray/Neutron Pair Distribution Function Analyses

The crystal structure of moganite from the Mogán formation on Gran Canaria has been re-investigated using high-resolution synchrotron X-ray diffraction (XRD) and X-ray/neutron pair distribution function (PDF) analyses. Our study for the first time reports the anisotropic atomic displacement parameters (ADPs) of a natural moganite. Rietveld analysis of synchrotron XRD data determined the crystal structure of moganite with the space group I2/a. The refined unit-cell parameters are a = 8.7363(8), b = 4.8688(5), c = 10.7203(9) Å, and β = 90.212(4)°. The ADPs of Si and O in moganite were obtained from X-ray and neutron PDF analyses. The shapes and orientations of the anisotropic ellipsoids determined from X-ray and neutron measurements are similar. The anisotropic ellipsoids for O extend along planes perpendicular to the Si-Si axis of corner-sharing SiO4 tetrahedra, suggesting precession-like movement. Neutron PDF result confirms the occurrence of OH over some of the tetrahedral sites. We postulate that moganite nanomineral is stable with respect to quartz in hypersaline water. The ADPs of moganite show a similar trend as those of quartz determined by single-crystal XRD. In short, the combined methods can provide high-quality structural parameters of moganite nanomineral, including its ADPs and extra OH position at the surface. This approach can be used as an alternative means for solving the structures of crystals that are not large enough for single-crystal XRD measurements, such as fine-grained and nanocrystalline minerals formed in various geological environments.

58 GEOSCIENCES↗

Crystal Structure of an Anisotropic Pyrope Garnet That Contains Two Cubic Phases

The crystal structure of two different samples of pyrope garnet, ideally Mg 3 Al 2 Si 3 O 12 , from South Africa was refined using the Rietveld method, space group Ia3¯d, and monochromatic synchrotron high-resolution powder X-ray diffraction (HRPXRD) data. Sample 1 from Wesselton Mine is a single cubic phase and is optically isotropic. Electron-probe microanalysis (EPMA) provided an average composition {Mg 2.30 Fe 2+ 0.26 Ca 0.42 Mn 2+ 0.02 } Σ3 [Al 1.53 Fe 3+ 0.06 Cr 3+ 0.40 Ti 4+ 0.01 Fe 2+ 0.01 ] Σ2 Si 3 O 12 , which contains a significant amount of Cr cations. The unit-cell parameter (Å) and bond distances (Å) are a = 11.56197(1) Å, average = 2.2985, Al-O = 1.9101(4), and Si-O = 1.6343(3) Å. Sample 2 from De Beers Diamond Mine has an average composition {Mg 2.33 Fe 2+ 0.33 Ca 0.33 Mn 2+ 0.01 } Σ3 [Al 1.73 Fe 3+ 0.12 Cr 3+ 0.06 Ti 4+ 0.05 Fe 2+ 0.05 ] Σ2 Si 3 O 12 and is a fine-scale intergrowth of two cubic phases. The weight percentage, unit-cell parameter (Å), and bond distances (Å) for phase 2a are 62.2(1)%, a = 11.56185(1) Å, average = 2.3006, Al-O = 1.9080(4), Si-O = 1.6334(4) Å. The corresponding values for phase 2b are 37.8(1)%, a = 11.53896(1) Å, average = 2.2954, Al-O = 1.9020(6), Si-O = 1.6334(6) Å. The two cubic phases in sample 2 cause the crystal to be optically anisotropic because of strain induce birefringence. The unit-cell parameter and bond distances for sample 1 are similar to those in phase 2a.

36 MATERIALS SCIENCE↗

Transforming Energy through Computational Excellence. Exascale Computing: Combustion; Simulating Effects of Fuel Injection Location in Supersonic Jet Engines

Computational tractable simulations using an adaptive-mesh-refinement solverfor compressible reacting flows help researchers understand how variations in fuel injection location within the supersonic flow cavity impacts combustion efficiency. By identifying the important physical determinants of the combustion processes, this study shows a promising pathway to improving flame stability and combustion efficiency, as well as reducing emissions.

adaptive mesh refinement↗

Computational Study of Variable Fuel-Air Ratio and Hydrogen Doping in a Rotating Detonation Engine

Rotating detonation engines (RDEs) hold great promise as a technology for powering aviation engines. The pressure gain that results from continuously traversing detonation waves enables high theoretical performance compared to deflagration-based combustion. However, the harsh operating conditions inside an RDE restrict the availability of suitable experimental diagnostics to study these systems in detail and the range of spatial and temporal scales render numerical simulations of these systems difficult. In this study, we have combined adaptive mesh refinement and a robust shock-capturing approach to study RDE systems computationally using the compressible, reactive solver PeleC. Multi-species transport along with compressible Navier-Stokes equations are solved in the model along with finite-rate chemistry. Complex geometries are represented using an embedded boundary method with second-order spatial accuracy and the system is evolved in time using a second-order Runge-Kutta method. We present studies of a methane-air RDE system that is doped with varying levels of hydrogen. Our analysis focusses on how fuel mixture and doping levels effect the detonation flame structure and the generation of multiple wave modes.

ADVANCED PROPULSION SYSTEMS↗

Numerical Study on the Effect of Methane Doping in Hydrogen-Air Rotating Detonation Engines for Various Temperatures and Pressures

Rotating detonation engines (RDEs) have gained attention as a promising technology for future aviation engines. However, the numerical studies of these systems pose severe challenges due to the broad range of spatial and temporal scales. In this study, we use an adaptive mesh refinement based compressible, reactive solver PeleC to resolve the broad range of scales and accurately capture shock and detonation waves using high-resolution numerical schemes. Multi- species transport along with compressible Navier-Stokes equations are solved in the model along with a finite-rate based chemistry model. Embedded boundary method is used to model the complex geometry consisting of discrete fuel nozzles and the combustion chamber. The fuel consisting predominantly of hydrogen is doped with varying levels of methane while air is used as the oxidizer. For a specified total pressure and temperature, the number of stable detonation waves is found to decrease with increasing methane concentration in the fuel mixture. Additionally, no stable detonation solutions are observed for methane composition higher than 20% by volume for the range of operating conditions studied (300-900K, 10 Atm). The increased presence of high temperature zones is also indicative of higher thermal NOx emissions at low methane concentrations. The effect of fuel-air mixture composition and temperature on the detonability, detonation wave structure, mode transitions and their stability are analyzed in this study in addition to their implication on NOx emission.

automatic mesh refinement↗

Simulations of Fuel-Air Mixing in a 7 Element Lean Direct Injection (LDI) Aviation Combustor

The increased accessibility of commercial aviation to the general population has in-creased the conventional jet fuel consumption, thereby causing increased CO2 emissions. One of the strategies to reduce the environmental impact is to use sustainable aviation fuels (SAFs). SAFs are jet fuels obtained from renewable resources such as biomass and waste streams such as plastics/municipal solid wastes. Even though the use of SAFs has been gaining importance over the last two decades, only a few numerical studies on their usage have been reported so far. In this study, two SAFs namely, Alcohol to Jet (ATJ) and Hydrotreated Esters and Fatty Acids (HEFA), are chosen and their use in a simple, Lead Direct Injection (LDI) combustor is studied. The focus of this study is limited to the thermo-chemical properties of the fuels and their impact on fuel distribution and evaporation. Hence, only non-reactive simulations are performed. A high accuracy, adaptive mesh refinement based low-Mach solver, PeleLMeX, is used for performing the simulations. The time-averaged flow results indicate faster evaporation of HEFA when compared to ATJ and Jet-A. The difference in fuel evaporation rates is attributed to the differences in fuel viscosity and surface tension properties and the concomitant variations in the fuel droplet size distributions.

adaptive mesh refinement↗

Simulations of Fuel-Air Mixing in a 7 Element Lean Direct Injection (LDI) Aviation Combustor: Preprint

The increased accessibility of commercial aviation to the general population has in-creased the conventional jet fuel consumption, thereby causing increased CO2emissions. One of the strategies to reduce the environmental impact is to use sustainable aviation fuels (SAFs). SAFs are jet fuels obtained from renewable resources such as biomass and waste streams such as plastics/municipal solid wastes. Even though the use of SAFs has been gaining importance over the last two decades, only a few numerical studies on their usage have been reported so far. In this study, two SAFs namely, alcohol to jet (ATJ) and hydrotreated esters and fatty acids (HEFA), are chosen and their use in a simple, lead direct injection (LDI) combustor is studied. The focus of this study is limited to the thermo-chemical properties of the fuels and their impact on fuel distribution and evaporation. Hence, only non-reactive simulations are performed. A high accuracy, adaptive mesh refinement based low-Mach solver, PeleLMeX, is used for performing the simulations. The time-averaged flow results indicate faster evaporation of HEFA when compared to ATJ and Jet-A. The difference in fuel evaporation rates is attributed to the differences in fuel viscosity and surface tension properties and the concomitant variations in the fuel droplet size distributions.

adaptive mesh refinement↗

A decade of curtailment studies demonstrates a consistent and effective strategy to reduce bat fatalities at wind turbines in North America

Abstract There is a rapid, global push for wind energy installation. However, large numbers of bats are killed by turbines each year, raising concerns about the impacts of wind energy expansion on bat populations. Preventing turbine blades from spinning at low wind speeds, referred to as curtailment, is a method to reduce bat fatalities, but drawing consistent inference across studies has been challenging. We compiled publicly available studies that evaluated curtailment at six wind energy facilities in North America across 10 years. We used meta‐regression of 29 implemented treatments to determine fatality reduction efficacy as well as sources of variation influencing efficacy. We also estimated species‐specific fatality reduction for three species that comprise most fatalities in North America: hoary bat ( Lasiurus cinereus ), eastern red bat ( Lasiurus borealis ) and silver‐haired bat ( Lasionycteris noctivagans ). We found that curtailment reduced total bat fatalities by 33% with every 1.0 ms −1 increase in curtailment wind speed. Estimates of the efficacy for the three target species were similar (hoary bats: 28% per ms −1 , 95% CI: 0.4%–48%, eastern red bats: 32% per ms −1 , 95% CI: 13%–47% and silver‐haired bats: 32% per ms −1 , 95% CI: 3%–53%). Across multiple facilities and years, a 5.0 ms −1 cut‐in speed was estimated to reduce total bat fatalities by an average of 62% (95% CI: 54%–69%). Mortality reductions at individual facilities in any given year were estimated to fall between 33%–79% (95% prediction interval). Inter‐annual differences rather than inter‐site or turbine characteristics accounted for most of the variation in efficacy rates. Species‐specific average mortality reduction at 5.0 ms −1 curtailment wind speed was 48% (95% CI: 24%–64%) for hoary bats, 61% (95% CI: 42%–74%) for eastern red bats and 52% (95% CI: 30%–66%) for silver‐haired bats. Practical implication . curtailment reduced bat mortality at wind turbines in this North American study. Efficacy increased proportionally as curtailment speed is raised, and patterns and rates of efficacy were similar across species. This indicates that curtailment is an effective strategy to reduce bat fatalities at wind energy facilities, but exploration of further refinements could both minimize bat mortality and maximize energy production.

17 WIND ENERGY↗