Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Timescale”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Stability of plutonium oxide nanoparticles in the presence of montmorillonite and implications for colloid facilitated transport

The physical and chemical stability of PuO 2 nanoparticles (intrinsic colloids) in groundwater will control their transport and may affect the performance of high-level radioactive waste repositories. In this study, we examined the chemical stability of two types of PuO 2 nanoparticles at both 25 and 80 °C. The “alkaline” PuO 2 nanoparticles were prepared by neutralizing an acidic Pu(IV) solution with dilute NaOH to pH ~9.5. The “acidic” PuO 2 nanoparticles were precipitated from 0.1 M nitric acid by heating a Pu(IV) solution at 60–80 °C for 30 min. The chemical stability of these PuO 2 nanoparticles was tested in the presence of montmorillonite, a common mineral in the environment and potentially relevant backfill material in some nuclear waste repository designs. The “alkaline” PuO 2 nanoparticles were found to be unstable over a timescale of months at both 25 and 80 °C, with elevated temperature enhancing their dissolution rates and sorption to montmorillonite. PuO 2 nanoparticle dissolution rates decreased with increasing Pu concentration, consistent with solution saturation. The “acidic” PuO 2 nanoparticles appeared to remain stable for much longer than the “alkaline” PuO2 nanoparticles. The “alkaline” PuO 2 nanoparticle dissolution rate constants were as high as 10 -11.9±0.4 mol m -2 s -1 and 10 -11.2±0.2 mol m -2 s -1 at 25 and 80 °C, respectively, while the “acidic” PuO 2 nanoparticle dissolution rate at 80 °C was 10 -13.5 mol m -2 s -1 . Based on transmission electron microscopy, the “alkaline” and “acidic” PuO 2 nanoparticles were of similar size (2.5–4.5 and 2–3 nm nanoparticles, respectively). However, the “acidic” PuO 2 nanoparticles formed more ordered nanoparticle aggregates. Our results suggest that the specific conditions experienced during PuO 2 nanoparticle formation could significantly affect the stability of PuO 2 in the presence of competing sorption processes and, in turn, the relative importance of intrinsic versus pseudocolloid transport in the environment.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗
DOE OSTI2020

Optical follow-up of gravitational wave triggers with DECam during the first two LIGO/VIRGO observing runs

Gravitational wave (GW) events detectable by LIGO and Virgo have several possible progenitors, including black hole mergers, neutron star mergers, black hole–neutron star mergers, supernovae, and cosmic string cusps. A subset of GW events is expected to produce electromagnetic (EM) emission that, once detected, will provide complementary information about their astrophysical context. To that end, the LIGO–Virgo Collaboration (LVC) sends GW candidate alerts to the astronomical community so that searches for their EM counterparts can be pursued. The DESGW group, consisting of members of the Dark Energy Survey (DES), the LVC, and other members of the astronomical community, uses the Dark Energy Camera (DECam) to perform a search and discovery program for optical signatures of LVC GW events. DESGW aims to use a sample of GW events as standard sirens for cosmology. Due to the short decay timescale of the expected EM counterparts and the need to quickly eliminate survey areas with no counterpart candidates, it is critical to complete the initial analysis of each night’s images as quickly as possible. We discuss our search area determination, imaging pipeline, and candidate selection processes. We review results from the DESGW program during the first two LIGO–Virgo observing campaigns and introduce other science applications that our pipeline enables.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Triboelectric backgrounds to radio-based polar ultra-high energy neutrino (UHEN) experiments

Here, in the hopes of observing the highest-energy neutrinos (E EeV) populating the Universe, both past (RICE, AURA, ANITA) and current (RNO-G, ARIANNA, ARA and TAROGE-M) polar-sited experiments exploit the impulsive radio emission produced by neutrino interactions. In such experiments, rare single event candidates must be unambiguously identified above backgrounds. Background rejection strategies to date primarily target thermal noise fluctuations and also impulsive radio-frequency signals of anthropogenic origin. In this paper, we consider the possibility that ‘fake’ neutrino signals may also be generated naturally via the ‘triboelectric effect.’ This broadly describes any process in which force applied at a boundary layer results in displacement of surface charge, leading to the production of an electrostatic potential difference V. Wind blowing over granular surfaces such as snow can induce such a potential difference, with subsequent coronal discharge. Discharges over timescales as short as nanoseconds can then lead to radio-frequency emissions at characteristic MHz–GHz frequencies. Using data from various past (RICE, AURA, SATRA, ANITA) and current (RNO-G, ARIANNA and ARA) neutrino experiments, we find evidence for such backgrounds, which are generally characterized by: (a) a threshold wind velocity which likely depends on the experimental trigger criteria and layout; for the experiments considered herein, this value is typically $\mathscr{O}$(10 m/s), (b) frequency spectra generally shifted to the low-end of the frequency regime to which current radio experiments are typically sensitive (100–200 MHz), (c) for the strongest background signals, an apparent preference for discharges from above-surface structures, although the presence of more isotropic, lower amplitude triboelectric discharges cannot be excluded.

Background rejection↗
DOE OSTI2022

Lifetime of actin-dependent protein nanoclusters

Protein nanoclusters (PNCs) are dynamic collections of a few proteins that spatially organize in nanometer-length clusters. PNCs are one of the principal forms of spatial organization of membrane proteins, and they have been shown or hypothesized to be important in various cellular processes, including cell signaling. PNCs show remarkable diversity in size, shape, and lifetime. In particular, the lifetime of PNCs can vary over a wide range of timescales. The diversity in size and shape can be explained by the interaction of the clustering proteins with the actin cytoskeleton or the lipid membrane, but very little is known about the processes that determine the lifetime of the nanoclusters. In this paper, using mathematical modeling of the cluster dynamics, we model the biophysical processes that determine the lifetime of actin-dependent PNCs. In particular, we investigated the role of actin aster fragmentation, which had been suggested to be a key determinant of the PNC lifetime, and we found that it is important only for a small class of PNCs. A simple extension of our model allowed us to investigate the kinetics of protein-ligand interaction near PNCs. We found an anomalous increase in the lifetime of ligands near PNCs, which agrees remarkably well with experimental data on RAS-RAF kinetics. In particular, analysis of the RAS-RAF data through our model provides falsifiable predictions and novel hypotheses that will not only shed light on the role of RAS-RAF kinetics in various cancers, but also will be useful in studying membrane protein clustering in general.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Spatiotemporal 4D Whole-cell Modeling of a Minimal Autotroph Reveals Central Carbon Metabolism Regulated Locally by Protein Megacomplexes via Post-translational Modifications under Light Disturbance

Photosynthetic microorganisms rely on multiple pathways in central carbon metabolism to adapt to fluctuating light and energy availability across diel cycles. Mechanistic insight into the regulatory dynamics of this adaptation requires integrating processes spanning disparate timescales, from rapid redox-dependent post-translational modifications (PTMs) to slower changes in protein expression and metabolic pathway usage. To address this complexity beyond genome-based inference and traditional modeling, we develop a whole-cell four-dimensional (3D + time) model of the marine cyanobacterium Prochlorococcus marinus MED4 that explicitly represents the spatial organization of enzymatic and molecular processes in central carbon metabolism under light perturbation. We employ a perturbation-based research design to experimentally generate time-series, multi-omics measurements that provide molecular descriptors and cryo-ET derived 3D segmented volumes as constraints for this dynamic 4D framework. The integration of experiments and modeling across defined light regimes enables quantitative validation of system-level responses and forecasting under distinct light disturbances. We test the hypothesis that light-dependent redox PTMs regulating the structural assembly of a protein megacomplex, the “dark complex,” modulate metabolic flux at a conserved regulatory node of the Calvin–Benson cycle (CBC) in cyanobacteria. Our model shows that subcellular spatial organization buffers rapid light-induced changes in thylakoid reaction rates, which are followed by redox-PTM-mediated sequestration or release of CBC enzymes in the dark complex, ultimately impacting carbon fixation dynamics within carboxysomes. Comparison with an equivalently parameterized well-mixed stochastic model demonstrates that post-translational regulation not only buffers transcriptional noise and diffusion-driven fluctuations but also stabilizes phenotypic outcomes, underscoring the importance of spatial heterogeneity in phenotypic robustness. This ability to probe adaptive, spatiotemporally resolved mechanisms in photosynthetic machinery and central carbon metabolism addresses a critical gap in genotype-to-phenotype inference and expands modeling and design capabilities for understudied or genetically intractable autotrophs such as P. marinus MED4.

Johnson, Connah G.↗
DOE OSTI2024

NO formation by N 2 /O 2 plasma catalysis: The impact of surface reactions, gas-phase reactions, and mass transport

Pathways and timescales relevant to facilitate plasma-assisted N 2 -O 2 reactions are assessed by measuring the consumption of plasma-derived N and the formation of NO in the gas phase and over Ag catalytic surfaces. These measurements are enabled by a setup that enables N 2 activation in an atmospheric pressure RF plasma jet, enables O 2 addition in the plasma afterglow, facilitates reactions over an Ag wire catalyst, and allows species density quantification by molecular beam mass spectrometry. Gas-phase reactions consume N but do not form NO with high selectivity. The presence of the non-porous Ag wire catalyst increases the rate of N conversion to NO, though mass transfer processes, not surface reactions, dictate the rate of N consumption. When O 2 concentrations and the ratio of the surface area of the catalyst to the void volume of the reactor are high (3–5 mol% O 2 , 10900 m –1 ), N conversion to NO reaches 100 % selectivity. When both N 2 and O 2 are fed through the plasma jet, gas-phase NO production increases 10×, although plasma and gas-phase processes do not exclusively produce NO. Above a threshold NO density, N cannot diffuse to the catalyst surface faster than it is consumed in the gas phase by reactions with NO. Furthermore, the use of heterogeneous catalysts to enhance plasma-driven N x O y formation and control N x O y product selectivity is limited to cases where diffusive transport of N from the gas phase to the catalyst surface is faster than consumption of N from gas-phase reactions with NO.

Engineering↗
DOE OSTI2024

Validating the rhenium proxy for rock organic carbon oxidation using weathering profiles

Chemical weathering over geological timescales acts as a source or sink of atmospheric carbon dioxide (CO 2 ), while influencing long-term redox cycling and atmospheric oxygen (O 2 ) at Earth's surface. There is a growing recognition that the oxidative weathering of rock organic carbon (OC petro ) can release more CO 2 than is locally drawn down by silicate weathering, and may vary due to changes in erosion and climate. The element rhenium (Re) has emerged as a proxy to track the oxidative weathering of OC petro , yet uncertainties in its application remain namely that we lack a systematic assessment of the comparative mobility of Re and OC petro during sedimentary rock weathering. Here we measure Re and OC petro loss across gradients in rock weathering at 9 global sites, spanning a range of initial OC petro values from ~0.2 % to >10 %. We use titanium to account for volume changes during weathering and assess Re and OC petro loss alongside major elements that reflect silicate (Na, Mg), carbonate (Ca, Mg) and sulfide (S) weathering. Across the dataset, Re loss is correlated with OC petro loss but not with loss of any other major element. Further, across the weathering profiles, the average molar ratio of OC petro to Re loss was 0.84 ± 0.15, with 8 out of 9 sites having a ratio >0.74. At one site (Marcellus Shale), the average ratio was lower at 0.58 ± 0.11. The excess loss of Re matches expectations that, typically, between ~0 and 20 % of the Re liberated by sedimentary rock weathering derives from silicate or sulfide phases, while some OC petro may be physically or chemically protected from weathering. Overall, our measurements provide validation for the Re proxy of OC petro oxidation and allow future work to further improve our knowledge of regional and global-scale rates of this important source of CO 2 in the geochemical carbon cycle.

58 GEOSCIENCES↗
DOE OSTI2022

Plutonium mobilization from contaminated estuarine sediments, Esk Estuary (UK)

Since 1952, liquid radioactive effluent containing 238-242 Pu, 241 Am, 237 Np, 137 Cs, and 99 Tc has been released with authorization from the Sellafield nuclear complex (UK) into the Irish Sea. This represents the largest source of plutonium (Pu) discharged in all western Europe, with 276 kg having been released. In the Eastern Irish Sea, the majority of the transuranic activity has settled into an area of sediments (Mudpatch) located off the Cumbrian coast. Radionuclides from the Mudpatch have been re-dispersed via particulate transport in fine-grained estuarine and intertidal sediments to the North-East Irish Sea, including the intertidal saltmarsh located at the mouth of the Esk Estuary. Saltmarshes are highly dynamic systems which are vulnerable to external agents (sea level change, erosion, sediment supply, and freshwater inputs), and their stability remains uncertain under current sea level rise projections and possible increases in storm activity. In this work, we examined factors affecting Pu mobility in contaminated sediments collected from the Esk Estuary by conducting leaching experiments under both anoxic and oxic conditions. Here, leaching experiments were conducted over a 9-month period and were periodically sampled to determine solution phase Pu via multicollector-inductively coupled plasma-mass spectrometry (MC-ICP-MS), and to measure redox indicators (Eh, pH and extractable Fe(II)). Microbial community composition was also characterized in the sediments, and at the beginning and end of the anoxic/oxic experiments. Results show that: 1) Pu leaching is about three times greater in solutions leached under anoxic conditions compared to oxic conditions, 2) the sediment slurry microbial communities shift as conditions change from anoxic to oxic, 3) Pu leaching is enhanced in the shallow sediments (0–10 cm depth), and 4) the magnitude of Pu leached from sediments is not correlated with total Pu, indicating that the biogeochemistry of sediment-associated Pu is spatially heterogeneous. These findings provide constraints on the stability of redox sensitive Pu in biogeochemically dynamic/transient environments on a timescale of months and suggests that anoxic conditions can enhance Pu mobility in estuarine systems.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Phonon-mediated lipid raft formation in biological membranes

Short-wavelength collective molecular motions, also known as phonons, have recently attracted much interest in revealing dynamic properties of biological membranes through the use of neutron and x-ray scattering, infrared and Raman spectroscopies, and molecular dynamics simulations. Experimentally detecting unique vibrational patterns such as, shear phonon excitations, viscoelastic crossovers, transverse acoustic phonon gaps, and continuous and truncated optical phonon modes in cellular membranes, to name a few, has proven non-trivial. Here, we review recent advances in liquid thermodynamics that have resulted in the development of the phonon theory of 1 liquids. The theory has important predictions regarding the shear vibrational spectra of fluids, namely the emergence of viscoelastic crossovers and transverse acoustic phonon gaps. Furthermore, we show that these vibrational patterns are common in soft (non-crystalline) materials, including, but not limited to liquids, colloids, liquid crystals (mesogens), block copolymers, and biological membranes. The existence of viscoelastic crossovers and acoustic phonon gaps define the self-diffusion properties of cellular membranes and provide a molecular picture of the transient nature of lipid rafts.1 Importantly, the timescales (picoseconds) for the formation and dissolution of transient lipid rafts match the lifetime of the formation and breakdown of interfacial water hydrogen bonds. Apart from acoustic propagating phonon modes, biological membranes can also support more energetic non-propagating optical phonon excitations, also known as standing waves or breathing modes. Importantly, optical phonons can be truncated due to the existence of finite size nanodomains made up of strongly correlated lipid-cholesterol molecular pairs. These strongly coupled molecular pairs can serve as nucleation centers for the formation of stable rafts at larger length scales, due to correlations of spontaneous fluctuations (Onsager’s regression hypothesis). Finally and importantly, molecular level viscoelastic crossovers, acoustic phonon gaps, and continuous and truncated optical phonon modes may offer insights as to how lipid-lipid and lipid-protein interactions enable biological function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
DOE OSTI2020

Slow manifold reduction for plasma science

The classical Chapman–Enskog procedure admits a substantial geometrical generalization known as slow manifold reduction. This generalization provides a paradigm for deriving and understanding most reduced models in plasma physics that are based on controlled approximations applied to problems with multiple timescales. In this Review we develop the theory of slow manifold reduction with a plasma physics audience in mind. In particular we illustrate (a) how the slow manifold concept may be used to understand breakdown of a reduced model over sufficiently-long time intervals, and (b) how a discrete-time analogue of slow manifold theory provides a useful framework for developing implicit integrators for temporally-stiff plasma models. For readers with more advanced mathematical training we also use slow manifold reduction to explain the phenomenon of inheritance of Hamiltonian structure in dissipation-free reduced plasma models. Various facets of the theory are illustrated in the context of the Abraham–Lorentz model of a single charged particle experiencing its own radiation drag. As a culminating example we derive the slow manifold underlying kinetic quasineutral plasma dynamics up to first-order in perturbation theory. This first-order result incorporates several physical effects associated with small deviations from exact charge neutrality that lead to slow drift away from predictions based on the leading-order approximation n e = Z i n i .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Simulated production of OH, HO 2 , CH 2 O, and CO 2 during dilute fuel oxidation can predict 1st-stage ignition delays

Chemical kinetics simulations are used to explore whether detailed measurements of relevant chemical species during the oxidation of very dilute fuels (less than 1 Torr partial pressure) in a high-pressure plug flow reactor (PFR) can predict autoignition propensity. We find that for many fuels the timescale for the onset of spontaneous oxidation in dilute fuel/air mixtures in a simple PFR is similar to the 1st-stage ignition delay time (IDT) at stoichiometric engine-relevant conditions. For those fuels that deviate from this simple trend, the deviation is closely related to the peak rate of production of OH, HO 2 , CH 2 O, and CO 2 formed during oxidation. We use these insights to show that an accurate correlation between simulated profiles of these species in a PFR and 1st-stage IDT can be developed using convolutional neural networks. Our simulations suggest that the accuracy of such a correlation is 10–50%, which is appropriate for rapid fuel screening and may be sufficient for predictive fuel performance modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Evaluation of semi-empirical soot models for nonpremixed flames with increased stoichiometric mixture fraction and strain

Due to the size and complexity of industrial-scale systems and computing limitations, semi-empirical soot models are often employed in CFD simulations rather than computationally-expensive detailed models. Many of these models were developed for specific applications with unique characteristic timescales and/or validated only under fuel-air combustion conditions. Hence, their use in different contexts, such as oxy-combustion, could lead to inaccurate predictions. In this study, twelve semi-empirical models (1-step or 2-step) are evaluated on their ability to respond to changes in stoichiometric mixture fraction (Z st ) and strain in a series of ethylene counterflow flames which span the experimental sooting-to-non-sooting (yellow to blue) transition. A unique approach of plotting soot formation rate in normalized, local equivalence ratio space is introduced to aid in the analysis. Results show that no existing model is able to predict a blue flame when Z st is increased beyond the experimentally-measured sooting limit. Many models give the counter-intuitive result of increasing peak and/or integrated soot volume fraction (svf) as Z st increased, contrary to experimental observations. Two models are able to successfully predict a blue flame upon increase of strain to the sooting limit value. For two-step models, the strong dependence of growth rate on surface area results in high sensitivity of svf to changing flame boundary conditions, but with mixed accuracy. There is still a significant need for a robust semi-empirical model which can accurately predict soot fraction in systems with oxygen-enrichment or variable strain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Eulerian thermo-mechanical simulations of heterogeneous solid propellants using an approximate projection method

We present a mathematical model that describes thermo-mechanical deformations and thermal gradients on the unsteady burning of a heterogeneous solid propellant. A scaling study shows that the deformations in solid at combustion timescales can be treated as quasi-static. The resulting thermo-mechanical formulation is formulated on a Cartesian grid and makes use of a weak form of Chorin-type projection method to deal with large difference in shear modulus of constituent materials. Further, a one-dimensional verification study is carried out by comparing numerical simulations with those of an analytical model. In addition, convergence studies for a two-dimensional propellant sandwich configuration are presented for the stress, velocity, and reference map components. Finally, simulations are carried out for a two-dimensional random propellant pack and the time-averaged burn rate and strains are reported.

42 ENGINEERING↗
DOE OSTI2020

Cloud based tool for analysis of chemical kinetic mechanisms

The size and complexity of chemical kinetic models used to simulate combustion of hydrocarbon fuels continue to increase with our improved understanding of the underlying physical mechanisms, along with the improved ability to utilize such models as the result of increased computational power and efficiency. As mechanisms grow beyond thousands of species and tens of thousands of reactions, it becomes increasingly difficult to manually check for errors and physical inconsistencies. In this work, we present several automated methods to check for such issues in the specification of chemical reaction models. First, we demonstrate how discontinuities in thermodynamic data can cause simulation difficulties manifested as long wall-clock times and convergence failures. To correct this type of problem, we describe an automated method for refitting thermodynamic parameters. We also outline several methods to check the timescales of reaction rate coefficients to ensure physical consistency. All the methods are made available through a web-based tool (https://combustiontools.llnl.gov) intended to aid mechanism developers and users to improve the accuracy and performance of fuel models used by the combustion community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Analysis of direct and indirect noise in a next-generation aviation gas turbine combustor

A large-eddy simulation (LES) of a next-generation combustor is performed to examine effects of this combustor concepts on direct and indirect combustion noise characteristics. Direct noise is computed considering the unsteady heat release predictions while entropy fluctuations in the downstream part of the combustor are used to estimate indirect noise through the transfer function of the outlet nozzle. A low-order acoustic reconstruction technique, which utilizes the Green’s function of the configuration, is developed to compute the acoustics inside the combustor. Comparisons to experimental results are reported, showing the reasonable accuracy of the LES and combustion noise computations. The direct noise spectrum peaks at frequencies around 3 kHz because of the fast timescales associated with the chemical reactions inside the combustor. In addition, the indirect noise is dominant at low frequencies, i.e., less than 400 Hz, because of the characteristics of the entropy spectrum at the outlet nozzle. Further, compared to a conventional rich-quench-lean (RQL) combustor configuration, significant differences in the acoustics are observed, and are explained by two specific technological novelties. Firstly, the direct noise spectrum peaks at significantly higher frequencies due to the combustor’s compactness. Secondly, the entropy inhomogeneties downstream of the combustion region are an order of magnitude smaller in amplitude due to the lack of dilution with cold air. This leads to a significant reduction in the indirect noise amplitude at low frequencies, where it is the dominant source of noise. These results suggest opportunities for advanced combustion technologies to achieve substantial reductions in combustion-related noise.

42 ENGINEERING↗