Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MIXTURES”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effect of particle shape on stratification in drying films of binary colloidal mixtures

The role of particle shape in evaporation-induced auto-stratification in polydisperse colloidal suspensions is explored with molecular dynamics simulations of mixtures of spheres and aspherical particles. A unified framework based on the competition between diffusion and diffusiophoresis is proposed to understand the effects of shape and size dispersity. In general, particles diffusing more slowly (e.g., larger particles) tend to accumulate more strongly at the evaporation front. However, larger particles have larger surface areas and therefore greater diffusiophoretic mobility. Hence, they are more likely to be driven away from the evaporation front via diffusiophoresis. For a rapidly dried bidisperse suspension containing small and large spheres, the competition leads to “small-on-top” stratification. Here, we employ a computational model in which the diffusion coefficient is inversely proportional to particle mass. For a mixture of spheres and aspherical particles with similar mass, the diffusion contrast is reduced, and the spheres are always enriched at the evaporation front as they have the smallest surface area for a given mass and, therefore, the lowest diffusiophoretic mobility. Furthermore, for a mixture of solid and hollow spheres that have the same outer radius and thus the same surface area, the diffusiophoretic contrast is suppressed, and the system is dominated by diffusion. Consequently, the solid spheres, which have a larger mass and diffuse more slowly, accumulate on top of the hollow spheres. Finally, for a mixture of thin disks and long rods that differ significantly in shape but have similar mass and surface area, both diffusion and diffusiophoresis contrasts are suppressed, and the mixture does not stratify.

Classical molecular dynamic simulations↗
DOE OSTI2024

A mixture of grass–legume cover crop species may ameliorate water stress in a changing climate

Climate change models predict increasing precipitation variability in the mid-latitude regions of Earth, generating a need to reduce the negative impacts of these changes on crop production. Despite considerable research on how cover crops support agriculture in a changing climate, understanding is limited of how climate change influences the growth of cover crops. We investigated the early development of two common cover crop species—crimson clover (Trifolium incarnatum) and rye (Secale cereale)—and hypothesized that growing them in the mixture would ameliorate stress from drought or waterlogging. This hypothesis was tested in a 25-day greenhouse experiment, where the two factors (species number and water stress) were fully crossed in randomized blocks, and plant responses were quantified through survival, growth rate, biomass production and root morphology. Water stress negatively influenced the early growth of these two species in contrasting ways: crimson clover was susceptible to drought while rye performed poorly under waterlogging. Per-plant biomass in rye was always greater in mixture than in monoculture, while per-plant biomass of crimson clover was greater in mixture under drought. Both species grew longer roots in mixture than in monoculture under drought, and total biomass of mixtures did not differ significantly from the more-productive monoculture (rye) in any water condition. In the face of increasingly variable precipitation, growing crimson clover and rye together has potential to ameliorate water stress, a possibility that should be further investigated in field experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Extreme compression of planetary gases: High-accuracy pressure-density measurements of hydrogen-helium mixtures above fourfold compression

Hydrogen (H 2 ) and helium (He), the most abundant elements in the universe, pose a unique challenge in measuring the equation of state of the mixture, owing to their differing physical properties. There remains a need for data with high enough precision to discriminate between existing equation of state (EOS) mix models in order to understand the internal structure of gas-giant planets. Here, we have measured the EOS of precompressed H 2 - He mixtures at conditions directly relevant to the planetary interiors using hypervelocity gas guns and Sandia’s Z machine with less than 10% uncertainty in density, enabling validation of mixture models. We precompressed 50:50 molar mixtures of H 2 -He to 0.1–0.2 GPa and directly measured particle velocity (in gas-gun experiments) and shock velocities (in Z-machine experiments). To complement the experimental efforts, we also computed the Hugoniots of precompressed H 2 -He mixtures using density-functional-theory-based molecular dynamics. Furthermore, we observe approximately 3- to 4.3-fold density compression at pressures up to 44 GPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Learning Molecular Mixture Property Using Chemistry-Aware Graph Neural Network

Recent advances in machine learning (ML) are expediting materials discovery and design. One significant challenge facing ML for materials is the expansive combinatorial space of potential materials formed by diverse constituents and their flexible configurations. This complexity is particularly evident in molecular mixtures, a frequently explored space for materials, such as battery electrolytes. Owing to the complex structures of molecules and the sequence-independent nature of mixtures, conventional ML methods have difficulties in modeling such systems. Here, we present MolSets, a specialized ML model for molecular mixtures, to overcome the difficulties. Representing individual molecules as graphs and their mixture as a set, MolSets leverages a graph neural network and the deep sets architecture to extract information at the molecular level and aggregate it at the mixture level, thus addressing local complexity while retaining global flexibility. We demonstrate the efficacy of MolSets in predicting the conductivity of lithium battery electrolytes and highlight its benefits in the virtual screening of the combinatorial chemical space. Published by the American Physical Society 2024

Zhang, Hengrui (ORCID:0000000231831654)↗
DOE OSTI2023

Development of Binder Jet Compatible Silicon Carbide Mixtures for Increased Sintered Density

Binder jet additive manufacturing can process a variety of ceramic powder feedstocks to manufacture near net shape green parts that are subsequently sintered to achieve densification. In particular, silicon carbide (SiC) is a promising ceramic material for heat exchangers in energy storage applications and binder jet is a well-suited process for manufacturing these heat exchangers due to the often geometrically complex designs. However, the need for good powder flowability in binder jet printing often limits the feedstock powder to be larger particles rather than that of ultra fine submicron powders used in conventional ceramic processing, which are more readily densified during sintering. The present work explores the use of various powder mixtures and sintering settings to identify strategies for achieving high density sintered parts in binder jet compatible systems without the use of infiltration. Liquid-phase sintering additives and bimodal powder mixtures are explored alongside different sintering temperatures, rates, atmospheres, and hold times. Two mixtures were found to be good candidates for further exploration: 1) 90 wt% 0.6 µm SiC + 6.25 wt% Al 2 O 3 + 3.75 wt% Y 2 O 3 mixture achieved the highest sintered relative density of 91.2%; and 2) a bimodal mixture comprised of 45 wt% 10 µm SiC + 45 wt% 0.6 µm SiC + 6.25 wt% Al 2 O 3 + 3.75 wt% Y 2 O 3 had the best powder flowability of the powders measured and a 75.1% sintered relative density.

25 ENERGY STORAGE↗
DOE OSTI2025

Finch: Toxicity Dose Response Curve Prediction of Chemical Compounds and Mixtures

A paradigm shift in chemical risk assessment is emphasizing mixture testing over single compound analysis, eliminating animal testing, and adopting advanced modeling approaches to understand mixture activity profiles. However, existing computational models largely focus on single chemicals, with few effective solutions for modeling complex mixtures that account for synergistic or antagonistic effects and multiple Modes of Action (MoA). Conventional methods like concentration addition (CA) and independent action (IA) are insufficient for this task as they are designed for simplistic interactions and struggle to account for the dynamic and multifaceted nature of chemical mixtures, such as overlapping MoA and non-linear interactions. Finch offers a novel approach utilizing deep learning (DL) embeddings and multi-task quantitative structure-activity relationship (QSAR) models to improve chemical exposure prediction. By leveraging molecular descriptors, physiochemical properties, and large language model (LLM) embeddings from SMILES inputs, Finch preserves critical information in a latent space thereby enhancing predictive accuracy. The multi-task learning aspect of Finch is highly advantageous, as it simultaneously optimizes multiple loss functions, leveraging all available data across tasks to develop generalized representations that effectively capture complex ingredient interactions within mixtures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Increasing Compressed Gas Energy Storage Density Using CO2–N2 Gas Mixture

This paper demonstrates a new method by which the energy storage density of compressed air systems is increased by 56.8% by changing the composition of the compressed gas to include a condensable component. A higher storage density of 7.33 MJ/m3 is possible using a mixture of 88% CO2 and 12% N2 compared to 4.67 MJ/m3 using pure N2. This ratio of gases representing an optimum mixture was determined through computer simulations that considered a variety of different proportions from pure CO2 to pure N2. The computer simulations are based on a thermodynamic equilibrium model that predicts the mixture composition as a function of volume and pressure under progressive compression to ultimately identify the optimal mixture composition (88% CO2 + 12% N2). The model and simulations predict that the optimal gas mixture attains a higher energy storage density than using either of the pure gases.

25 ENERGY STORAGE↗
DOE OSTI2020

A new approach for measuring the carbon and oxygen content of atmospherically relevant compounds and mixtures

Abstract. Due to its complexity, gas- and particle-phase organic carbon in the atmosphere is often classified by its bulk physicochemical properties. However, there is a dearth of robust, moderate-cost approaches to measure the bulk chemical composition of organic carbon in the atmosphere. This is particularly true for the degree of oxygenation, which critically affects the properties and impacts of organic carbon but for which routine measurement approaches are lacking. This gap has limited the understanding of a wide range of atmospheric components, including particulate matter, the mass of which is monitored worldwide due to its health and environmental effects but the chemical characterization of which requires relatively high capital costs and complex operation by highly trained technical personnel. In this work, we demonstrate a new approach to estimate the mass of carbon and oxygen in analytes and mixtures that relies only on robust, moderate-cost detectors designed for use with gas chromatography. Organic compounds entering a flame ionization detector were found to be converted with approximately complete efficiency to CO2, which was analyzed downstream using an infrared detector to measure the mass of carbon analyzed. The ratio of the flame ionization detector (FID) signal generated to CO2 formed (FID∕CO2) was shown to be strongly correlated (R2=0.89) to the oxygen-to-carbon ratio (O∕C) of the analyte. Furthermore, simple mixtures of analytes behaved as the weighted average of their components, indicating that this correlation extends to mixtures. These properties were also observed to correlate well with the sensitivity of the FID estimated by structure activity relationships (quantified as the relative effective carbon number). The relationships between measured FID∕CO2, analyte O∕C, and FID sensitivity allow the estimation of one property from another with <15 % error for mixtures and <20 % error for most individual analytes. The approach opens the possibility of field-deployable, autonomous measurement of the carbon and oxygen content of particulate matter using time-tested, low-maintenance detectors, though such an application would require some additional testing on complex mixtures. With some instrumental modifications, similar measurements on gas-phase species may be feasible. Moreover, the potential expansion to additional gas chromatography detectors may provide concurrent measurement of other elements (e.g., sulfur, nitrogen).

42 ENGINEERING↗
DOE OSTI2022

Impact of Variable Gas Mixtures on Bubble Size Distribution and Mass Transfer in Gas Fermentation Reactors

Gas fermentation has emerged as a promising new technology for the generation of fuels and chemicals from mixtures of greenhouse and energy rich gas streams (CO2/CH4/H2/CO) via microbial bioreaction. Example pathways include biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Effective Gas-liquid mass-transfer is an important physical phenomenon that determines the design and scale-up of these systems. There is currently a knowledge-gap regarding bubble-size distributions when using a mixture of gases with vastly different properties, which can have a significant impact on overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a detailed multiphase computational-fluid-dynamics (CFD) model to study the impact of gas-mixtures on overall mass-transfer in bubble column and air-lift reactors. The CFD tool previously developed by the authors (1) for simulating aerobic fermentation reactors at scale is used in this study. The Reynolds-averaged mass, momentum, energy, and species transport equations are solved for interpenetrating gas and liquid phase in this model. We use a population balance-based bubble-size-distribution model that is validated against small-scale experiments in our solver. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient is analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗
DOE OSTI2025

In situ deformation of antigorite-olivine two-phase mixtures: Implications for dynamics and seismic anisotropy in the mantle wedge

Water released from hydrous minerals in subducting slabs reacts with the overlying plate, resulting in widespread serpentinization in the mantle wedge. Deformation of serpentinized peridotites has been invoked to explain forearc seismic anisotropy, yet studies of the mechanical properties and deformation behaviors of serpentine-bearing multiphase aggregates remain limited. Here we deformed olivine-antigorite mixtures containing 70, 50, and 20 vol.% of antigorite at 2.5 – 7.6 GPa, 673 K and strain rates of ∼10 –5 –10 –4 s –1 . Elasto-viscoplastic self-consistent simulations, constrained by synchrotron X-ray diffraction (XRD) data, were used to estimate lattice strain, stress–strain partitioning, crystallographic preferred orientations (CPO) development, and aggregate strength. Selected run products were also analyzed by electron backscatter diffraction for comparison with the CPO results obtained from XRD experiments. We found olivine transitions from A- or B-type to C-type when antigorite fraction drops to 20 vol.%, coinciding with a microstructural change from interconnected weak layers to a load-bearing framework (LBF). An additional run on a sample Atg50/Ol50 with preexisting microstructures suggested the formation of LBF was promoted by these microstructures, although the preexisting antigorite CPO has been overprinted at 20.8 % strain and could be erased completely by subsequent deformation in nature. Estimated viscosity of the two-phase mixtures suggests that low-degree serpentinization (≤20 %) in the mantle wedge may increase the strength of olivine-rich peridotite and hinder slab-mantle decoupling, whereas high-degree serpentinization (≥50–70 %) weakens the peridotite and favors decoupling if sufficient viscosity contrast (>10) develops. Seismic anisotropy shows a nonlinear dependence on antigorite fraction: antigorite CPO governs the anisotropy of the mixtures with ≥50 vol.% antigorite, whereas olivine CPO dominates at low fractions (∼20 vol.%). The presence of pre-existing microstructures reduces seismic anisotropy of the deformed mixtures, however the persistence of pre-existing CPO in actively subducting slabs remains uncertain, making their significance over geological timescales questionable.

Crystallographic preferred orientation↗
DOE OSTI2021

Multi-objective optimization of low-GWP mixture composition and heat exchanger circuitry configuration for improved system performance and reduced refrigerant flammability

The recently introduced hydrofluoroolefin (HFO) refrigerants, including R1234yf and R1234ze(E), have significantly lower global warming potentials (GWPs) than traditional hydrofluorocarbon (HFC) refrigerants like R410A. However, prior tests show that direct drop-in of pure R1234yf or R1234ze(E) into equipment designed for R410A results in a decrease in heat exchanger capacity and the system coefficient of performance. The primary reason is the lower in-tube heat transfer performance of R1234yf and R1234ze(E) compared with that of R410A. To address this issue, previous studies have mixed the mildly flammable HFC R32 with HFOs to improve system performance, with HFC R125 also added to suppress flammability. Previous studies selected compositions based on simple cycle analyses and did not consider modifications of the heat exchanger circuitry configuration to adapt to the new refrigerants. This study presents a novel multi-objective optimization approach to design a refrigerant composition that maximizes energy efficiency within flammability and GWP limits. The approach in this work simultaneously optimizes mixture composition and heat exchanger circuitry configuration. A case study on a rooftop unit indicates that, compared with mixture-only optimization, simultaneous optimization of mixture and heat exchanger circuitry yields a 5.9% improvement in cycle efficiency and a 48.6% reduction in refrigerant flammability with a GWP of 268. Circuitry optimization using refrigerants with different temperature glides shows that the larger the temperature glide is, the larger EER improvement is obtained. The results show that zeotropic blends with a large temperature glide are more sensitive to the refrigerant circuitry than pure refrigerants and may suffer significant performance degradation with subpar heat exchanger circuitry design. The proposed optimization approach is generally applicable to mixtures with any number of components. Using this approach to design a HVAC system can yield higher system efficiency within flammability and GWP constraints.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Volatility Basis Set Distributions and Viscosity of Organic Aerosol Mixtures: Insights from Chemical Characterization Using Temperature-Programmed Desorption–Direct Analysis in Real-Time High-Resolution Mass Spectrometry

Quantitative assessment of gas-particle partitioning of individual components within complex atmospheric organic aerosols (OA) mixtures is critical for predicting and comprehending the formation and evolution of OA particles in the atmosphere. This investigation leverages previously documented data obtained through a temperature programmed desorption - direct analysis in real time – high resolution mass spectrometry (TPD-DART-HRMS) platform. This methodology facilitates the bottom-up construction of volatility basis set (VBS) distributions for constituents found in three biogenic secondary organic aerosol (SOA) mixtures produced through the ozonolysis of a-pinene, limonene, and ocimene. The apparent enthalpies (ΔH*, kJ mol -1 ) and saturated vapor mass concentrations (C T *, µg∙m -3 ) of individual SOA components, determined as a function of temperature (T, K), facilitated an assessment of changes in VBS distributions and gas-particle partitioning with respect to T and atmospheric total organic mass loadings (tOM, µg∙m -3 ). Further, the VBS distributions reveal distinct differences in volatilities among monomers, dimers, and trimers, enabling their categorization into separate volatility bins. At the ambient temperature of T = 298 K, only monomers efficiently partition between gas and particle phases across a broad range of atmospherically relevant total organic mass loadings (tOM) values of 1–100 µg∙m -3 . Partitioning of dimers and trimers becomes notable only at T > 360 K and T > 420 K, respectively. The viscosity of SOA mixtures is assessed using a bottom-up calculation approach, incorporating the input of elemental formulas, ΔH*, C T *, and particle-phase mass fractions of the SOA components. Through this approach, we are able to accurately estimate the variations in SOA viscosity that result from the evaporation of its components. These variations are, in turn, influenced by atmospherically relevant changes in tOM and T. Comparison of the calculated SOA viscosity and diffusivity values with literature reported experimental results shows close agreement, thereby validating the employed calculation approach. These findings underscore the significant potential for TPD-DART-HRMS measurements in enabling the untargeted analysis of organic molecules within OA mixtures. This approach facilitates quantitative assessment of their gas-particle partitioning and allows for the estimation of their viscosity and condensed-phase diffusion, thereby contributing valuable insights to atmospheric models.

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

Modeled Emission of Hydroxyl and Ozone Reactivity from Evaporation of Fragrance Mixtures

Volatile chemical products (VCPs) account for increasing fractions of organic carbon emitted to the atmosphere, particularly in urban areas. Fragrances are potentially reactive components that are added to many VCPs. To better constrain these emissions, 11 commercially available liquid fragrance mixtures were characterized for their composition and their evaporation modeled. Here, emissions of mass, hydroxyl reactivity, and ozone reactivity were estimated by modeling under four different scenarios. Fragrance compounds were generally less than one-half the mass of fragrance mixtures, with the balance comprised of solvents and plasticizers and unresolved mass thought to be dominated by plasticizers. The results showed that terpenes and terpenoids account for nearly all of the emitted mass and reactivity while only comprising ~10% w/w on average of the liquid fragrance mixtures. Most of the reactivity is emitted within hours, with ozone reactivity evolving more rapidly than OH reactivity and comprised almost entirely of terpenes. Limonene, a common fragrance constituent, dominates the reactivity of emitted carbon. Generally, 20–40% of the potential hydroxyl reactivity contained in the fragrance mixture does not evaporate on time scales sufficient to have an impact on local or regional air quality.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Prediction of Self-Diffusion in Binary Fluid Mixtures Using Artificial Neural Networks

Artificial neural networks (ANNs) were developed to accurately predict the self-diffusion constants for individual components in binary fluid mixtures. The ANNs were tested on an experimental database of 4328 self-diffusion constants from 131 mixtures containing 75 unique compounds. The presence of strong hydrogen bonding molecules may lead to clustering or dimerization resulting in non-linear diffusive behavior. To address this, self- and binary association energies were calculated for each molecule and mixture to provide information on intermolecular interaction strength and were used as input features to the ANN. An accurate, generalized ANN model was developed with an overall average absolute deviation of 4.1%. Forward input feature selection reveals the importance of critical properties and self-association energies along with other fluid properties. Additional ANNs were developed with subsets of the full input feature set to further investigate the impact of various properties on model performance. The results from two specific mixtures are discussed in additional detail: one providing an example of strong hydrogen bonding and the other an example of extreme pressure changes, with the ANN models predicting self-diffusion well in both cases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Filling the Gaps: A Bayesian Mixture Model for Imputing Missing Soil Water Content Data

ABSTRACT Soil water content (SWC) data are central to evaluating how soil moisture varies over time and space and influences critical plant and ecosystem functions, especially in water‐limited drylands. However, sensors that record SWC at high frequencies often malfunction, leading to incomplete timeseries and limiting our understanding of dryland ecosystem dynamics. We developed an analytical approach to impute missing SWC data, which we tested at six eddy flux tower sites along an elevation gradient in the southwestern United States. We impute missing data as a mixture of linearly interpolated SWC between the observed endpoints of a missing data gap and SWC simulated by an ecosystem water balance model (SOILWAT2). Within a Bayesian framework, we allowed the relative utility (mixture weight) of each component (linearly interpolated vs. SOILWAT2) to vary by depth, site and gap characteristics. We explored “fixed” weights versus “dynamic” weights that vary as a function of cumulative precipitation, average temperature, and time since the start of the gap. Both models estimated missing SWC data well ( R 2 = 0.70–0.88 vs. 0.75–0.91 for fixed vs. dynamic weights, respectively), but the utility of linearly interpolated versus SOILWAT2 values depended on site and depth. SOILWAT2 was more useful for more arid sites, shallower depths, longer and warmer gaps and gaps that received greater precipitation. Overall, the mixture model reliably gap‐fills SWC, while lending insight into processes governing SWC dynamics. This approach to impute missing data could be adapted to accommodate more than two mixture components and other types of environmental timeseries.

Ogle, Kiona [School of Informatics, Computing, and↗
DOE OSTI2020

Investigating the effect of different shielding gas mixtures on microstructure and mechanical properties of 410 stainless steel fabricated via large scale additive manufacturing

Metal Big Area Additive Manufacturing (mBAAM) offers the potential to fabricate large scale tools at high deposition rates (15 lb/h+). 410 martensitic steel is a potential tooling material, owing to its low cost, good machinability and reasonable printability. During the mBAAM process, the shielding gas can have a significant impact on the material properties as well as the process cost. Therefore, the current study aims to understand the effect of different shielding gas mixtures on large-scale additive manufacturing of 410 martensitic stainless steel. We show that an argon mixture with 3% nitrogen gas produced the best performance in terms of maximum hardness and tensile strength, with much less scatter in tensile strength. He-Ar-CO 2 or tri-mix shielded samples showed a low tensile strength with wide scatter, due to stabilized delta ferrite in microstructure during printing. Both tri-mix and Ar-CO 2 shielded samples showed slightly higher porosity. Thus, we recommend the use of argon-3% nitrogen as a shielding gas mixture for processing 410 steel for tool applications, based on the relatively low cost of this gas mixture and the resulting higher hardness, higher dimensional stability, and lower porosity.

36 MATERIALS SCIENCE↗
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 OSTI2020

Combustion of Silane-Hydrogen-Nitrous oxide Mixtures: Analysis of Laminar Flame Propagation and Condensed Products

The laminar burning rate, the explosion pressure, and the pressure rise coefficient have been measured for silane-nitrous oxide-argon mixtures using the spherically expanding flame technique in a constant volume combustion chamber. For these three parameters, the values obtained were higher than for hydrogen-nitrous oxide-argon and typical hydrocarbon-based mixtures. As an example, a maximum burning rate of 1800 g/m 2 s was measured at 101 kPa, whereas under similar conditions, a maximum burning rate around 950 g/m 2 s has been reported for hydrogen-nitrous oxide-argon mixtures. Here, the condensed combustion products formed in the silane-nitrous oxide-argon flames were sampled and characterized using electron micrograph, electronic diffraction, energy-dispersive spectroscopy, and X-ray powder diffraction. For all equivalence ratios, silica spherical particles with a mean diameter in the range 200-300 nm were observed. In addition, for mixtures with Φ ≥2.2, silicon nanowires were formed. X-ray diffraction experiments showed that the silicon nanowires are composed of a material referred to as “metal silicon” characterized by a cubic structure (lattice parameter: a = 5.425 Å) with the Fm-3m space group.

42 ENGINEERING↗