SEARCH · Search NASA

Results for “mixture models”

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.

SourceResearch connectionFrom yearThrough year

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

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Coupling Remote Sensing With a Process Model for the Simulation of Rangeland Carbon Dynamics

Rangelands provide significant environmental benefits through many ecosystem services, which may include soil organic carbon (SOC) sequestration. However, quantifying SOC stocks and monitoring carbon (C) fluxes in rangelands are challenging due to the considerable spatial and temporal variability tied to rangeland C dynamics as well as limited data availability. We developed the Rangeland Carbon Tracking and Management (RCTM) system to track long-term changes in SOC and ecosystem C fluxes by leveraging remote sensing inputs and environmental variable data sets with algorithms representing terrestrial C-cycle processes. Bayesian calibration was conducted using quality-controlled C flux data sets obtained from 61 Ameriflux and NEON flux tower sites from Western and Midwestern US rangelands to parameterize the model according to dominant vegetation classes (perennial and/or annual grass, grass-shrub mixture, and grass-tree mixture). The resulting RCTM system produced higher model accuracy for estimating annual cumulative gross primary productivity (GPP) (R 2 > 0.6, RMSE <390 g C m -2 ) relative to net ecosystem exchange of CO 2 (NEE) (R 2 > 0.4, RMSE <180 g C m -2 ). Model performance in estimating rangeland C fluxes varied by season and vegetation type. The RCTM captured the spatial variability of SOC stocks with R 2 = 0.6 when validated against SOC measurements across 13 NEON sites. Model simulations indicated slightly enhanced SOC stocks for the flux tower sites during the past decade, which is mainly driven by an increase in precipitation. Future efforts to refine the RCTM system will benefit from long-term network-based monitoring of vegetation biomass, C fluxes, and SOC stocks.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Thermochemical Modeling in Molten Fluoride Salts for Radionuclide Speciation

An important aspect of the licensing process for nuclear reactors is providing a reasonable assurance of safety to the general public. This includes modeling potential radionuclide releases from the reactor during normal operations and accident scenarios, which is known as the reactor’s source term. A new class of advanced (non-LWR) reactors are being developed which use molten salts as the coolant fluid. Because the molten salt coolant represents a credited barrier for radionuclide transport between the fuel and the environment, a necessary aspect of mechanistic source term (MST) analysis for the KP-FHR is modeling the thermochemistry of molten salts. Provided here is a review of the theory of the thermodynamic principles governing multicomponent phase equilibria, the background of molten salt thermochemistry research, and a summary of the thermochemical data relevant to the KP-FHR coolant salt, Li 2 BeF 4 , commonly referred to as “FLiBe”. A review of literature is followed by a brief introduction to methods that can be used to model the thermochemical behavior of molten salt mixtures. The methodology outlined is based on the use of a commercial thermodynamic modeling software called FactSage, which is one of only a few available softwares based on the modified quasichemical model(MQM), which is the recommended solution model for molten salts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

A unified non-equilibrium phase change model for injection flow modeling

The homogenous relaxation model (HRM) is one of the most widely used models to describe the liquid- gas phase transition. However, in its original formulation, it is unable to handle multispecies vapor-liquid equilibrium (VLE), which limits its applicability to single-component fluids. In this work, a unified non-equilibrium phase change model that considers the VLE of multicomponent mixtures is proposed building upon the HRM's structure. A time factor is introduced to mimic the effect of different phase change timescales due to different mechanisms, e.g., cavitation, flash-boiling, and evaporation. Here to assess the model's performance, computational fluid dynamics simulations of the internal and near-nozzle injection flow of the Engine Combustion Network's Spray G injector were performed using the nine-component PACE-20 fuel with both the unified model and the original HRM. The predicted fuel density in the near-nozzle region matched well with X-ray tomography measurements. The simulation results indicated that, whereas the HRM failed to capture the vaporization due to convective mixing between the fuel and ambient gas, the unified model performed well in predicting the mixing-driven vaporization and the corresponding evaporative cooling. Further comparisons using the nine-component fuel formula and a single-component fuel surrogate demonstrated the unified model's ability to predict preferential vaporization, which affects the predictions of local mixture composition and rate of vaporization. Finally, it is shown that the unified model is capable of representing multiple phase change mechanisms, and the relaxation time factor plays an important role in determining the degree of phase change due to the different mechanisms.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2026

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Unveiling the Dynamics of KRAS4b on Lipid Model Membranes

Abstract Small GTPase proteins are ubiquitous and responsible for regulating several processes related to cell growth and differentiation. Mutations that stabilize their active state can lead to uncontrolled cell proliferation and cancer. Although these proteins are well characterized at the cellular scale, the molecular mechanisms governing their functions are still poorly understood. In addition, there is limited information about the regulatory function of the cell membrane which supports their activity. Thus, we have studied the dynamics and conformations of the farnesylated KRAS4b in various membrane model systems, ranging from binary fluid mixtures to heterogeneous raft mimics. Our approach combines long time-scale coarse-grained (CG) simulations and Markov state models to dissect the membrane-supported dynamics of KRAS4b. Our simulations reveal that protein dynamics is mainly modulated by the presence of anionic lipids and to some extent by the nucleotide state (activation) of the protein. In addition, our results suggest that both the farnesyl and the polybasic hypervariable region (HVR) are responsible for its preferential partitioning within the liquid-disordered (Ld) domains in membranes, potentially enhancing the formation of membrane-driven signaling platforms. Graphic Abstract

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Dynamics and chemical mode analysis of plasma thermal-chemical instability

The stability of the weakly ionized plasma and the transition from a stable homogeneous discharge to unstable filaments play an important role in gas laser physics, plasma-assisted combustion, chemical reforming, and material synthesis. Here, theoretical stability analysis and thermal-chemical mode analysis were performed to understand the mechanism of plasma thermal-chemical instability by using a zero-dimensional plasma system with both simplified and detailed chemical kinetics of H 2 /O 2 /N 2 mixtures. The plasma dynamic and kinetic models accounted for multiple physical mechanisms in the chemically-reactive weakly ionized plasma, including ionization, attachment/detachment, recombination, vibrational and electronic energy relaxation, convective and diffusive species/heat removal, Joule heating, and detailed chemical kinetics. An analytical criterion and the explosive mode species/temperature pointers were formulated while the representative active species were identified for different thermal-chemical modes. The results showed that in addition to the classical thermal-ionization mechanism, various chemical modes from chemical heat imbalance and elementary kinetics significantly modified the time dynamics and the stability of the weakly ionized plasma. The present analysis provides insights and guidance to control plasma instability using chemical kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Implications of real-gas behavior on refractive index calculations for optical diagnostics of fuel–air mixing at high pressures

Three models to compute the refractive index of gaseous mixtures at real gas conditions are presented with the purpose to improve the accuracy of state relationships between refractive index and thermodynamic properties. Models are compared with experimental data to determine one that is applicable to high pressure mixtures with non-ideal thermodynamic behavior near or above a fluids’ critical point. The optimal model is applied to analyze adiabatic thermal mixing of fuel and air at typical diesel engine conditions. Our results show that the ideal gas mixture law is appropriate in the vapor region, assuming it is known, for example, from experimental measurements. Finally, the model is applied for binary fuel-air mixing at supercritical conditions to demonstrate its full potential.

02 PETROLEUM↗

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

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↗

DOE OSTI2021

Surfactants and cloud droplet activation: A systematic extension of Köhler theory based on analysis of droplet stability

The activation of aerosol particles to form cloud droplets, a necessary first step in cloud formation, controls much of the impact that aerosols have on clouds and climate. Recently, there has been a surge of interest in extending the Köhler theory of cloud droplet activation to include surface active (typically organic) as well as water-soluble (typically inorganic) aerosol components, but a systematic framework for doing this has yet to be developed. Here, we apply a droplet stability analysis to this end. Ideal and Szyszkowski–Langmuir surfactant models are analyzed to demonstrate the new approach, but the underlying theoretical framework is fundamental and model free. A key finding is that superficial densities at the cloud activation threshold (Köhler maximum) are significantly sub-monolayer, with fractional coverage ranging from 69% to 85% for the organic compounds and mixtures studied. The result, significant for model inventories of cloud condensation nuclei, is a weakening of the surfactant effect relative to expectations based on bulk sample measurements. Analytical results are obtained for the loci of Köhler maxima and applied to aerosol mixtures containing an arbitrary number of water-soluble and surfactant components.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Measurement of diamond nucleation rates from hydrocarbons at conditions comparable to the interiors of icy giant planets

We present measurements of the nucleation rate into a diamond lattice in dynamically compressed polystyrene obtained in a pump-probe experiment using a high-energy laser system and in situ femtosecond x-ray diffraction. Different temperature-pressure conditions that occur in planetary interiors were probed. For a single shock reaching 70 GPa and 3000 K no diamond formation was observed, while with a double shock driving polystyrene to pressures around 150 GPa and temperatures around 5000 K nucleation rates between 10 29 and 10 34 m –3 s –1 were recorded. These nucleation rates do not agree with predictions of the state-of-the-art theoretical models for carbon-hydrogen mixtures by many orders of magnitude. Our data suggest that there is significant diamond formation to be expected inside icy giant planets like Neptune and Uranus.

43 PARTICLE ACCELERATORS↗

DOE OSTI2022

Thermal properties of HfN-MoW surrogate cermet fuel for nuclear thermal propulsion

Here, the thermal properties of MoW-HfN, a surrogate cermet for MoW-UN nuclear thermal propulsion fuel, were characterized over a wide range of elevated temperatures. Thermal diffusivity, coefficient of thermal expansion (CTE), and heat capacity were measured. Optical and scanning electron microscopy were performed to characterize the microstructure and draw structure–property correlations. The thermal diffusivity was obtained using the laser flash method. Diffusivity values ranged from about 0.18 cm 2 /s at 200°C to 0.15 cm 2 /s at 1800°C. The CTE was measured using push-rod dilatometry up to 1600°C, giving values between 6.0 and 9.0 μm/m. A scientific rationalization of the effective material properties is made using the rule-of-mixtures and other effective properties models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Hydrogenous Content in Residues from Aqueous Chloride Pu and Am Recovery Operations

Byproduct residues generated from Aqueous Chloride Pu and Am recovery operations at Los Alamos National Laboratory have been investigated and determined to contain high weight % chloride salts and minimal hydrogen content. Hydrogen content is the result of limited surface area adsorption of ambient moisture. The assumption used in technical analysis NCS-TECH-19-029 to model potentially hydrogenous items as mixtures of polyethylene is shown to be overly conservative. Our results and literature precedent indicate that water adsorption of 1-3 % by weight is expected for hydroxide cakes and Pu oxides, and a conservative bounding assumption for these materials is 5% water by weight. LOH results for a dissolution heel stored long term showed 6.25% weight loss, which was conservatively assumed to be water. A conservative bounding assumption for dissolution heels is 8 wt.% water content. This equates to 0.55 % hydrogen by weight for hydroxide cakes, and 0.88% hydrogen by weight for dissolution heels. The bulk composition of dissolution heels reflects the composition of the original byproduct feed material, chloride salts of potassium, sodium, magnesium and calcium. The bulk composition of calcined hydroxide precipitation residues is composed of potassium chloride as a result of potassium hydroxide neutralization of hydrochloric acid solutions. Metal hydroxides generated during the hydroxide precipitation process are decomposed under the conditions of calcination to metal oxide forms.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Reduced-order condensed-phase kinetic models for polyethylene, polypropylene and polystyrene thermochemical recycling

Thermochemical recycling of plastic waste (PW) into chemicals and energy vectors requires coupling particle and reactor-scale simulations to accurate condensed phase pyrolysis mechanisms for each constituent. This work proposes a methodology to derive reduced-order condensed-phase kinetic models from validated semi-detailed kinetic mechanisms. Two types of kinetic models are obtained for polyethylene (PE), polypropylene (PP) and polystyrene (PS): reduced semi-detailed models and multi-step fully lumped ones. These families offer different compromises between accuracy and computational cost. The former employ 50–100 gas + liquid species and describe both the radical degradation and the detailed carbon distribution of the products. Conversely, the latter involves 5–10 species per polymer tracking only the main petroleum cuts. The kinetic mechanisms are complemented by the definition of thermochemical properties of gas, liquid, and solid-phase species, accounting for phase-transitions through pseudo-chemical reactions. Model validations are performed by comparison with experimental data and the original semi-detailed mechanisms in terms of mass loss, heat fluxes and product distribution profiles. The resulting CHEMKIN-like condensed-phase models are attached as Supplementary Material and as a GitHub repository. Extending the proposed approach to other polymers and coupling it with existing subsets in the CRECK kinetic framework (e.g., biomass, PVC, PET) offers a powerful tool to model thermochemical recycling of PW and biomass/PW mixtures.

kinetics↗

DOE OSTI2026

Application of an Empirical Density Law via Python for Aqueous Plutonium Chloride Systems for MCNP6

Criticality safety models for aqueous plutonium chloride systems often contain a significant bias due to assumptions in material compositions. These systems are currently modeled as a fictitious metal-water mixture because little is known about the true solution density. Furthermore, no predictive density tools or capabilities for modeling aqueous plutonium chloride systems are approved for use at Los Alamos National Laboratory. Recent density measurements of this ternary system (PuCl 3 -HCl-H 2 O) have allowed for the development of a more realistic density law, which is applied in this work via an empirical method based in Python. This tool, entitled PuCS (Plutonium Chloride Solution tool) may be used to determine solution density and composition based on the plutonium content, acid content, and temperature for MCNP6 inputs. PuCS has been found to predict density within 2% of experimental data. In conclusion, MCNP6 calculations have found that crediting minimal amounts of free acid (0.5 M) may correspond to a ~12% decrease in peak reactivity in comparison to current modeling methods.

42 ENGINEERING↗

DOE OSTI2021

Investigation of post-injection strategies for diesel engine Catalyst Heating Operation using a vapor-liquid-equilibrium-based spray model

Most multidimensional engine simulations spend much time solving for non-equilibrium spray dynamics (atomization, collision, vaporization). However, their accuracy is limited by significant grid dependency, and the need for extensive calibration. This is critical for modeling cold-start diesel fuel post injections, which occur at low temperatures and pressures, far from typical model validation ranges. At the same time, resolving micron-scale spray phenomena would render full Eulerian multiphase calculations prohibitive. In this study, an improved phase equilibrium based approach was implemented and assessed for simulating diesel catalyst heating operation strategies. A phase equilibrium solver based on the model by Yue and Reitz [1] was implemented: a fully multiphase CFD solver is employed with an engineering-size engine grid, and fuel injection is modeled using the standard Lagrangian parcels approach. Mass and energy from the liquid parcels are released to the Eulerian multiphase mixture according to an equilibrium-based liquid jet model. An improved phase equilibrium solver was developed to handle large real-gas mixtures such as those from accurate chemical kinetics mechanisms. The liquid-jet model was improved such that momentum transfer to the Eulerian solver better reproduces the physical spray jet structure. Validation of liquid/vapor penetration predictions showed that the model yields accurate results with very limited tuning and low sensitivity to the few calibration constants. In-cylinder simulations of diesel catalyst heating operation strategies showed that capturing spray structure is paramount when short, transient injection pulses and low temperatures are present. Furthermore, the EP model provides improved predictions of post-injection spray structure and ignitability, while conventional spray modeling does not capture the increase of liquid penetration during the expansion stroke. Finally, the only important EP model calibration constant, Cliq, does not affect momentum transfer, but it changes the local charge cooling distribution through the local energy transfer, which makes it candidate to additional research. The results confirm that non-equilibrium spray processes do not need to be resolved in engineering simulations of high-pressure diesel sprays.

42 ENGINEERING↗

DOE OSTI2023

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid–liquid mixtures

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

An experimental and kinetic modeling study of cyclopentane and dimethyl ether blends

Cyclopentane is a suitable naphthene, or cycloalkane, in a palette for multi-component gasoline surrogate fuels due to its presence in market fuels and its relevance to alkyl substituted cyclopentanes also present. However, the previous oxidation studies of cyclopentane have primarily focused on neat mixtures. Blending cyclopentane with dimethyl ether in this work therefore serves to inform our understanding of, and improve predictive models for, multi-component mixtures. Here, the auto-ignition of cyclopentane/dimethyl ether blends was studied in a high-pressure shock tube and in a rapid compression machine. A wide range of temperatures (650 – 1350 K) and elevated pressures of 20 and 40 bar were studied at equivalence ratios of 0.5, 1.0 and 2.0 in air for two blending ratios (30/70 and 70/30 mole% cyclopentane/di-methyl ether mixtures). A detailed kinetic model for cyclopentane was revised to capture the measured ignition delay times and apparent heat release rates in this study. Literature ignition delay time, jet-stirred reactor, and laminar burning velocity measurements of neat cyclopentane were used as additional validation. Improvements to the kinetic model were based on recent literature studies related to sub-models including cyclopentene and cyclopentadiene which allowed the removal of previous local rate-constant optimizations. Low temperature reactivity of cyclopentane was found to be controlled by the branching ratio between concerted elimination of HO͘ 2 and the strained formation of Q&#775OOH radicals in agreement with previous studies. In this study, the low branching ratio of Q&#775OOH formation increases the influence of a competing consumption pathway for cyclopentyl-peroxy (CPTO͘ 2 J) radicals. The sensitivity of the simulated ignition delay times to the formation of cyclopentyl hydroperoxide (CPTO 2 H), from CPTO͘ 2 J and HO͘ 2 , is discussed. The current model is used to analyze the influence of dimethyl ether on the reactivity of cyclopentane in the context of previous literature studies of dimethyl ether binary blends with ethanol and toluene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗