Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixture modeling”

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 379 records · Page 21

Designing Molten Salt Eutectics: A Combined Thermodynamic Modeling and Machine Learning Approach

Designing stable electrolytes with target properties is an important challenge in realizing next generation energy storage devices. Molten salt eutectics-based electrolytes are known for their stability with minimal parasitic reactions when compared to traditional organic electrolytes and are an attractive option for different battery chemistries. The operating temperature of the molten salt batteries depends on the melting temperature of the eutectic and hence there is a necessity to discover novel low melting temperature molten salt eutectic mixtures for energy storage applications. In this work we develop a high throughput computational screening approach for molten salt mixtures using thermodynamic modeling and machine learning (ML). COSMO-SAC model and ML approaches were independently developed based on the existing experimental data and these models were further used to predict the eutectic melting temperature and composition of several new binary, ternary, and quaternary mixtures. We show that combining ML and thermodynamic modeling strategies is effective in exploring the vast design space of molten salt mixtures.

Thermodynamics↗

Robust and optimal alignment of high-dimensional data using maximum likelihood estimation through a random sample consensus framework

Abstract Correcting spatial orientations of groups of high-dimensional data sets such that they are all in a consistent coordinate system is often a time-consuming and error-prone process. Automation of this process can be accomplished by using Generalized Procrustes Analysis to estimate the relative orientations among a population of high-dimensional data sets. A least squares Procrustes solution is applied through a maximum likelihood estimation and random sample consensus framework for robustness. The likelihood model is comprised of a mixture distribution where inliers are modeled using t -distribution and outliers from a uniform distribution. Applications will focus on a synthetic data set that emulates triaxial acceleration data and also real shock data from a population of triaxial accelerometers. Outliers represent either non-rigid body responses, environmental noise, and/or sensor and data acquisition issues. The intended application for the methodology is to robustly automate the rotation of populations of experimentally collected triaxial accelerometer data sets to a single global coordinate system.

LOSAC↗

Bayesian model mixing with multireference energy density functional

Reliably predicting nuclear properties across the entire chart of isotopes is important for applications ranging from nuclear astrophysics to superheavy science to nuclear technology. To this day, however, all the theoretical models that can scale at the level of the chart of isotopes remain semiphenomenological. Because they are fitted locally, their predictive power can vary significantly; different versions of the same theory provide different predictions. Bayesian model mixing takes advantage of such imperfect models to build a local mixture of a set of models to make improved predictions. Earlier attempts to use Bayesian model mixing for mass table calculations relied on models treated at single-reference energy density functional level, which fail to capture some of the correlations caused by configuration mixing or the restoration of broken symmetries. In this study we have applied Bayesian model mixing techniques within a multireference energy density functional (MR-EDF) framework. We considered predictions of two-particle separation energies from particle number projection or angular momentum projection with four different energy density functionals—a total of eight different MR-EDF models. We used a hierarchical Bayesian stacking framework with a Dirichlet prior distribution over weights together with an inverse log-ratio transform to enable positive correlations between different models. We found that Bayesian model mixing provides significantly improved predictions compared to the participating models. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A robust and efficient generalized real cascade model

This paper describes a new separation cascade model that allows for arbitrary cascade connectivity and includes the effects of mixing, material losses and generalized separation element performance. The model also accommodates mixtures with an arbitrary number of component species. The mathematical and numerical framework of this model includes both time-dependent and steady state cascade simulation. Robust numerical approximations of the continuous mathematical model are accomplished using the Sundials suite of nonlinear algebraic and differential-algebraic equation solvers. Verification of the model is accomplished using manufactured solutions while validation is demonstrated by comparison to published experimental results and benchmarks.

36 MATERIALS SCIENCE↗

Gesture Based Control and EMG Decomposition

This paper presents two probabilistic developments for use with Electromyograms (EMG). First described is a new-electric interface for virtual device control based on gesture recognition. The second development is a Bayesian method for decomposing EMG into individual motor unit action potentials. This more complex technique will then allow for higher resolution in separating muscle groups for gesture recognition. All examples presented rely upon sampling EMG data from a subject's forearm. The gesture based recognition uses pattern recognition software that has been trained to identify gestures from among a given set of gestures. The pattern recognition software consists of hidden Markov models which are used to recognize the gestures as they are being performed in real-time from moving averages of EMG. Two experiments were conducted to examine the feasibility of this interface technology. The first replicated a virtual joystick interface, and the second replicated a keyboard. Moving averages of EMG do not provide easy distinction between fine muscle groups. To better distinguish between different fine motor skill muscle groups we present a Bayesian algorithm to separate surface EMG into representative motor unit action potentials. The algorithm is based upon differential Variable Component Analysis (dVCA) [l], [2] which was originally developed for Electroencephalograms. The algorithm uses a simple forward model representing a mixture of motor unit action potentials as seen across multiple channels. The parameters of this model are iteratively optimized for each component. Results are presented on both synthetic and experimental EMG data. The synthetic case has additive white noise and is compared with known components. The experimental EMG data was obtained using a custom linear electrode array designed for this study.

Wheeler, Kevin R.↗

Statistical mechanics of light elements at high pressure. VII - A perturbative free energy for arbitrary mixtures of H and He

A model free energy is presented which accurately represents results from 45 high-precision Monte Carlo calculations of the thermodynamics of hydrogen-helium mixtures at pressures of astrophysical and planetophysical interest. The free energy is calculated using free-electron perturbation theory (dielectric function theory), and is an extension of the expression given in an earlier paper in this series. However, it fits the Monte Carlo results more accurately, and is valid for the full range of compositions from pure hydrogen to pure helium. Using the new free energy, the phase diagram of mixtures of liquid metallic hydrogen and helium is calculated and compared with earlier results. Sample results for mixing volumes are also presented, and the new free energy expression is used to compute a theoretical Jovian adiabat and compare the adiabat with results from three-dimensional Thomas-Fermi-Dirac theory. The present theory gives slightly higher densities at pressures of about 10 megabars.

Hubbard, W. B.↗

A phenomenological rate of injection model for predicting fuel injection with application to mixture formation in light-duty diesel engines

Fuel injection rate laws are one of the most important pieces of information needed when modeling engine combustion with computational fluid dynamics. In this study, a simple phenomenological model of a common-rail injector was developed and calibrated for the Bosch CRI2.2 platform. The model requires three tunable parameter fits, making it relatively easy to calibrate and suitable for injector modeling when high-fidelity information about the internal injector’s geometry and electrical circuit details are not available. Each injection pulse is modeled as a sequence of up to four stages: an injection needle mechanical opening transient, a full-lift viscous flow inertial transient, a Bernoulli steady-state stage, and a needle descent transient. Parameters for each stage are obtained as polynomial fits from measured injection rate properties. The model enforces total injected mass, and the intermediate stages are only introduced if the injection pulse duration is long enough. Experimental rates of injection from two separate campaigns on the same injector were used to calibrate the model. The model was first validated against measured injection rate laws featuring pilot injections, short partially premixed combustion pulses, and conventional diesel combustion injection strategies. Then, it was employed as an input to engine computational fluid dynamics simulations, which were run to simulate experiments of mixture formation in an optically accessible light-duty diesel engine. It was found that, though simple, this model is capable of predicting both pilot and main injection pulse mass flow rates well: the simulations yielded accurate predictions of in-cylinder equivalence ratio distributions from injection strategies for both partially premixed combustion and pilot injections. Also, once calibrated, the model produced appropriate results for a wide range of injected mass and rail pressure values. Finally, it was observed that usage of such a relatively simple model can be a good choice when high-fidelity injection rate input and highly detailed information of the injector’s geometry and operation are not available, particularly as noticeable discrepancies can be present also among different experimental campaigns on similar hardware.

42 ENGINEERING↗

On Kinetics Modeling of Vibrational Energy Transfer

Two models of vibrational energy exchange are compared at equilibrium to the elementary vibrational exchange reaction for a binary mixture. The first model, non-linear in the species vibrational energies, was derived by Schwartz, Slawsky, and Herzfeld (SSH) by considering the detailed kinetics of vibrational energy levels. This model recovers the result demanded at equilibrium by the elementary reaction. The second model is more recent, and is gaining use in certain areas of computational fluid dynamics. This model, linear in the species vibrational energies, is shown not to recover the required equilibrium result. Further, this more recent model is inconsistent with its suggested rate constants in that those rate constants were inferred from measurements by using the SSH model to reduce the data. The non-linear versus linear nature of these two models can lead to significant differences in vibrational energy coupling. Use of the contemporary model may lead to significant misconceptions, especially when integrated in computer codes considering multiple energy coupling mechanisms.

Gilmore, John O.↗

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↗

Supersonic combustion of a silane/methane mixture

A test program was conducted to evaluate the ignition and combustion characteristics of a 20 volume percent silane/methane mixture using a model scramjet combustor with Mach 2 combustor entrance conditions. The test gas total temperature was varied from 1200 to 3900 R. The mixture autoignited at a total temperature of 1650 R. This autoignition temperature can be contrasted with 2330 R for hydrogen and 1350 R for a 20 percent silane/hydrogen mixture in similar hardware. Methane without the silane additive did not autoignite in this configuration at temperatures as high as 3900 R.

Northam, G. B.↗

Bayesian calibration of bubble size dynamics applied to CO 2 gas fermenters

To accelerate the scale-up of gaseous CO 2 fermentation reactors, computational models need to predict gas-to-liquid mass transfer which requires capturing the bubble size dynamics, i.e. bubble breakup and coalescence. However, the applicability of existing models beyond air–water mixtures remains to be established. Here, an inverse modeling approach, accelerated with a neural network surrogate, calibrates the breakup and coalescence closure models, that are used in class methods for population balance modeling (PBM). The calibration is performed based on experimental results obtained in a CO 2 -air–water-coflowing bubble column reactor. Bayesian inference is used to account for noise in the experimental dataset and bias in the simulation results. To accurately capture gas holdup and interphase mass transfer, the results show that the breakage rate needs to be increased by one order of magnitude. In conclusion, the inferred model parameters are then used on a separate configuration and shown to also improve bubble size distribution predictions.

09 BIOMASS FUELS↗

Thermodynamics of Thomas-Fermi screened Coulomb systems

We obtain in closed analytic form, estimates for the thermodynamic properties of classical fluids with pair potentials of Yukawa type, with special reference to dense fully ionized plasmas with Thomas-Fermi or Debye-Hueckel screening. We further generalize the hard-sphere perturbative approach used for similarly screened two-component mixtures, and demonstrate phase separation in this simple model of a liquid mixture of metallic helium and hydrogen.

Firey, B.↗

Modeling Henry's law and phase separations of water–NaCl–organic mixtures with solvation and ion-pairing

Empirical measurements of solution vapor pressure of ternary acetonitrile (MeCN) H 2 O–NaCl–MeCN mixtures were recorded, with NaCl concentrations ranging from zero to the saturation limit, and MeCN concentrations ranging from zero to an absolute mole fraction of 0.64. After accounting for speciation, the variability of the Henry's law coefficient at vapor–liquid equilibrium (VLE) of MeCN ternary mixtures decreased from 107% to 5.1%. Solute speciation was modeled using a mass action solution model that incorporates solute solvation and ion-pairing phenomena. Two empirically determined equilibrium constants corresponding to solute dissociation and ion pairing were utilized for each solute. When speciation effects were considered, the solid–liquid equilibrium of H 2 O–NaCl–MeCN mixtures appear to be governed by a simple saturation equilibrium constant that is consistent with the binary H 2 O–NaCl saturation coefficient. Further, our results indicate that the precipitation of NaCl in the MeCN ternary mixtures was not governed by changes in the dielectric constant. Our model indicates that the compositions of the salt-induced liquid–liquid equilibrium (LLE) boundary of the H 2 O–NaCl–MeCN mixture correspond to the binary plateau activity of MeCN, a range of concentrations over which the activity remains largely invariant in the binary water–MeCN system. Broader comparisons with other ternary miscible organic solvent (MOS) mixtures suggest that salt-induced liquid–liquid equilibrium exists if: (1) the solution displays a positive deviation from the ideal limits governed by Raoult's law; and (2) the minimum of the mixing free energy profile for the binary water-MOS system is organic-rich. Finally, this work is one of the first applications of speciation-based solution models to a ternary system, and the first that includes an organic solute.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Kinetic and Transport Modeling for Entry Flows in Hydrogen-Helium Atmospheres

Several different thermochemical models for H2/He mixtures are used to simulate a 1-D shock in a viscous CFD code at conditions relevant for hypersonic entry into the ice and gas giants. These models range from a 6-species, 1-temperature (1-T) model, to more complex 17 and 25-species state-to-state (StS) models. Corresponding radiance profiles for each of these models are then computed using a radiation code, and results are compared to experimental data from the NASA Ames Electric Arc Shock Tube (EAST) facility. Overall, the 6-species model with a quasi-steady state (QSS) solver for atomic H and modified ionization rates, along with the 17 and 25-species StS models, are found to reproduce the magnitude of the radiance seen in the EAST experiments reasonably accurately. Two additional changes to the modeling of species diffusion fluxes are also investigated: the exclusion of a ∇W term, as well as the inclusion of Soret/ Dufour effects. It is determined that these changes to the transport modeling have a negligible impact on the flowfield and therefore the radiance predictions for the H2/He shock conditions considered in this study.

Outer planets↗

Evaluation of possible reactivation of undetected faults during CO 2 injection

Geologic storage of carbon dioxide can efficiently contribute to reduce greenhouse gas emissions to the atmosphere. Two major hazards of this technology are leakage towards the ground surface and fluid-induced seismicity. While major faults may be detected and avoided during site characterization, undetected subseismic faults could be encountered once injection has started. This paper investigates leakage and reactivation of undetected faults through coupled thermal–hydraulic–mechanical modeling. The simulations are performed using a recently developed sequential simulator, TOUGH-Pylith, that allows to accurately account for the thermodynamics of brine-CO 2 injection mixtures, and to model faults as surfaces of discontinuity using state-of-the-art fault friction laws. The simulator is benchmarked against well-known analytical solutions and subsequently applied to investigate two cases of CO 2 injection injection close to undetected faults under normal and strike-slip faulting regimes. Furthermore, although the scenarios are generic and represent unfavorable conditions, they suggest that leakage could occur and that undetected faults could trigger minor seismic events. Therefore, careful site characterization and continuous monitoring during operations should be always performed.

58 GEOSCIENCES↗