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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Effects of spatiotemporal plasma power distribution on the modeling of ignition kernel evolution in quiescent and turbulent methane/air mixtures

Abstract The present work improves a phenomenological plasma-assisted combustion model by integrating the spatiotemporal distribution of plasma power density, thereby considering the evolution of plasma streamers in the modeling, and subsequently, better predicting the ignition kernel evolution. The improved phenomenological model is validated against experiments representing the plasma discharge and post-discharge ignition kernel evolution. Specifically, the new model demonstrates a more accurate prediction of ultrafast gas heating and O 2 dissociation during the plasma discharge, compared to the original model. In addition, the new model is found to closely match the experimental pressure wave and heated channel profiles post-discharge without the need for tuning the energy deposition (unlike the original model), highlighting its accuracy of post-discharge ignition kernel dynamics. The improved phenomenological model is then employed to investigate ignition kernel evolution for a stoichiometric methane-air discharge across various discharge gap configurations. Simulations reveal a non-uniform temperature and streamer distribution progressing from the electrode tips toward the center, contrasting uniform cylindrical discharges previously described in the original model. Streamer propagation is observed to be faster for larger gaps when maintained at the same average electric field for different discharge gaps. The tendency of smaller gaps to produce detached toroidal ignition kernels is observed, while larger gaps promote cylindrical and attached ignition kernels. Interactions between successive ignition kernels from consecutive discharges varied significantly, with the smallest gap (1 mm) promoting the quenching of the preceding ignition kernel due to the initial kernel–kernel separation. The intermediate gap (2 mm) promotes detached kernel growth. In contrast, in the largest gap (4 mm), kernels consistently combine and expand attached to electrodes. The impact of homogeneous isotropic turbulence is also explored, showing the persistence of ignition kernels early on but eventually quenching due to enhanced radical and heat losses with pronounced turbulence intensity.

Johnson, Praise Noah↗
DOE OSTI2024

Uncertainty quantification of mass models using ensemble Bayesian model averaging

Developments in the description of the masses of atomic nuclei have led to various nuclear mass models that provide predictions for masses across the whole chart of nuclides. These mass models play an important role in understanding the synthesis of heavy elements in the rapid neutron capture ( r ) process. However, it is still a challenging task to estimate the size of uncertainty associated with the predictions of each mass model. In this work, a method called ensemble Bayesian model averaging (EBMA) is introduced to quantify the uncertainty of one-neutron separation energies (S 1 n ) which are directly relevant in the calculations of r -process observables. Here, this Bayesian method provides a natural way to perform model averaging, selection, and uncertainty quantification, by combining the mass models as a mixture of normal distributions whose parameters are optimized against the experimental data, employing the Markov chain Monte Carlo method using the no-u-turn sampler. The EBMA model optimized with all the experimental S 1 n from the AME2003 nuclides are shown to provide reliable uncertainty estimates when tested with the new data in the AME2020.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Calibration of the 2-Phase Bubble Tracking Model for Liquid Mercury Target Simulation with Machine Learning Surrogate Models

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory is one of the most powerful accelerator-driven neutron sources in the world. The intense protons strike on SNS's mercury target to provide bright neutron beams, which also leads to severe fluid-structure interactions inside the target. Prediction of resultant loading on the target is difficult particularly when helium gas is injected into mercury to reduce the loading and mitigate the pitting damage on vessel walls. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this multi-physics problem. However, several uncertain parameters in the R-P model require intensive simulations to determine their optimal values. With the help of machine learning and the measured target strain, we have studied the major uncertain parameters in this R-P model and developed a framework to identify optimal parameters that significantly reduce the discrepancy between simulations and experimental strains. The preliminary results show the possibility of using this mercury/helium mixture and surrogate models to predict a better match of target strain response when the helium gas is injected.

Lin, Lianshan↗
DOE OSTI2025

Inverse Thermodynamics: Designing Interactions for Targeted Phase Behavior

The traditional goal of inverse self-assembly is to design interactions that drive particles toward a desired target structure. However, achieving successful self-assembly also requires tuning the thermodynamic conditions under which the structure is stable. In this work, we extend the inverse design paradigm to explicitly address this challenge by developing a framework for inverse thermodynamics, i.e., the design of interaction potentials that realize specific thermodynamic behavior. As a step in this direction, using patchy particle mixtures as a model system, we demonstrate how precise control over both bonding topology and bond energetics enables the programming of targeted phase behavior. In particular, we establish design principles for azeotropic demixing and show how to create mixtures that exhibit azeotropy at any prescribed composition. Our predictions are validated through Gibbs-ensemble simulations [Panagiotopoulos, Mol. Phys. 1987, 61, 813−826]. These results highlight the necessity of coupling structural design with thermodynamic engineering, and provide a blueprint for controlling complex phase behavior in multicomponent systems.

Azeotropes↗
DOE OSTI2021

Quantitative modeling of carcinogenesis induced by single beams or mixtures of space radiations using targeted and non-targeted effects

Abstract Ionizing radiations encountered by astronauts on deep space missions produce biological damage by two main mechanisms: (1) Targeted effects (TE) due to direct traversals of cells by ionizing tracks. (2) Non-targeted effects (NTE) caused by release of signals from directly hit cells. The combination of these mechanisms generates non-linear dose response shapes, which need to be modeled quantitatively to predict health risks from space exploration. Here we used a TE + NTE model to analyze data on APC (1638N/+) mouse tumorigenesis induced by space-relevant doses of protons, 4 He, 12 C, 16 O, 28 Si or 56 Fe ions, or γ rays. A customized weighted Negative Binomial distribution was used to describe the radiation type- and dose-dependent data variability. This approach allowed detailed quantification of dose–response shapes, NTE- and TE-related model parameters, and radiation quality metrics (relative biological effectiveness, RBE, and radiation effects ratio, RER, relative to γ rays) for each radiation type. Based on the modeled responses for each radiation type, we predicted the tumor yield for a Mars-mission-relevant mixture of these radiations, using the recently-developed incremental effect additivity (IEA) synergy theory. The proposed modeling approach can enhance current knowledge about quantification of space radiation quality effects, dose response shapes, and ultimately the health risks for astronauts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Application of high-speed, species-specific chemiluminescence imaging for laminar flame speed and Markstein length measurements in spherically expanding flames

A high-speed chemiluminescence imaging diagnostic system has been developed and characterized for measuring species-specific properties of spherically expanding flames, with the objective of obtaining synchronized measurements of multiple flame properties from a single flame experiment. The imaging system comprises a high-speed camera coupled to a high-speed image intensifier and optical filters to isolate emission from desired excited species. Here this study validates the diagnostic scheme specifically for measuring laminar flame speed and Markstein length from chemiluminescence at three wavelength bands corresponding to OH*, CH*, and total chemiluminescence, respectively. Mixtures of methane-air and methane/ethane-air were tested at initial conditions of atmospheric pressure and room temperature. An under-utilized methodology for computing unburned Markstein length from burned Markstein length was implemented. Flame speed modeling of each mixture was conducted in Chemkin using AramcoMech 2.0. No significant differences were observed between the results derived from the three emission wavelength bands for the flame speeds nor the Markstein lengths. For both of these properties, excellent agreement was observed between the experimental results and the values reported in the literature. Because laminar flame speed and Markstein length show no dependency on the imaging wavelength, the present results imply that multiple flame properties, such as laminar flame speed, Markstein length, flame thickness, and excited-species intensity profiles, can all be acquired simultaneously for the same flame using this method.

42 ENGINEERING↗
DOE OSTI2023

Droplet burning of n-butyl acetate/n-heptane mixtures in the standard atmosphere: experiments and numerical modeling

Blending petroleum fuels with oxygenates is a common approach to stem the depletion of crude oil while also mitigating the impact of their combustion on the environment. It has recently been considered that n-butyl acetate (BA, C6H12O2 , boiling point of 399K) could be a viable oxygenate additive to diesel fuel. In this application it is important to determine the influence of the fractional amount of BA on burning of the mixture. This presentation considers this problem from an experimental and computational approach using n-heptane (C7H14 , boiling point of 372K) as an essentially single component surrogate for diesel fuel. The burning configuration considered was ostensibly spherical symmetry as promoted by burning droplets under conditions where forced and natural convection effects were minimized. In this configuration the droplet and flame are concentric spheres and soot aggregates are trapped in a shell structure between the droplet and flame. Measurements were made of the droplet, flame and soot shell diameters through video imaging of the burning process. The experimental results were used to validate a detailed numerical model of the spherically symmetric droplet burning process that incorporated a consistent combustion kinetic mechanism for BA/heptane mixtures, comprised of 402 species and 16,872 reactions inclusive of soot chemistry. A model for soot formation was also included in the simulation to predict formation of the soot shell. The simulations agreed well with measured droplet and flame diameters. It was found both experimentally and computationally that adding BA to heptane had a minimal effect on the mixture burning rate, while the flame was positioned closer to the droplet surface and the simulated soot volume fraction decreased. The implications of these results are discussed for using BA as a potential additive to diesel fuel.

09 BIOMASS FUELS↗
DOE OSTI2020

First Principle based Thermal Conductivity Modeling of Metal Alloys with Applications on Uranium Alloys

Thermal conductivity is an important materials property related to heat transport, which is essential to many applications, ranging from thermoelectrics to nuclear reactor materials. High- quality thermal conductivity data is critical to these materials and their associated technologies, including the nuclear fuel materials like traditional oxide fuels and metallic uranium (U) fuels, e.g. U-Zr and U-Mo alloys. And thermal conductivity modeling is widely used to interpolate or extrapolate experimental data, to obtain high-quality thermal conductivity data over wide temperature and composition range, and to understand the impacts of different factors including defects and microstructures. However, a general thermal conductivity model for metal alloys which can work regrading different phase components is still missing. Also, for U alloys, a thermal conductivity model working with different phases is needed. Therefore, in this work, we developed thermal conductivity models for metal alloys, based on ab-initio calculations, semi-classical physics rules, and limited experimental data. The goals of these models are to help obtain high- quality thermal conductivity data within a little experimental input as possible, have models that are extendable to varying microstructures, including different types of irradiation effects, and provide mechanistic understanding of heat transfer in the modeled alloys. In this work, our model solves several challenges in the development. The DFT-BTE approach is applied to decrease the reliance on experimental data and make the model extendable to composition changes and defects. A practical and efficient way to combine the DFT inputs with physics rules is pointed out in this work too. A staged approach, which starts from simple cases of elemental metals, then extends to solid solutions, different compound phases, and multi-phase mixtures, is presented to work with the metal alloys in different phase components. Our models are demonstrated respectively on aU for elemental metal model, on U-Zr and U-Mo alloys in aU temperature range for multi-phase mixture model, on high-temperature U-Zr, U-Nb, and U-Mo alloys in the body-centered cubic phase for concentrated solid solution model, and on irradiated U-Mo alloys for irradiated metal alloy model. All models show great agreement with experimental data. In these demonstrations, our model shows its advantages compared to previous empirical fitting model. Our model requires fewer experimental data, due to the inputs from DFT. The quantitative insights into the different physical factors are provided in our model, as it incorporates the electron and phonon scattering mechanisms. Our model also can be extended to incorporate the contributions of point defects, grain boundaries, and noble gas bubbles, to integrate their effects on heat transfer. This model can serve as both a foundation for understanding the more complex thermal conductivity of realistic U alloy fuels and a useful tool to guide further modeling of thermal conductivity to aid materials and device design.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Thermodynamic modeling of double azeotropy systems

Azeotropes are boiling liquid mixtures which have identical vapor and liquid compositions. Known as polyazeotropy, azeotropic mixtures may exhibit more than one azeotrope at a given temperature or pressure. While the thermodynamic condition for azeotrope formation in vapor-liquid equilibrium is well established, thermodynamic modeling of mixtures exhibiting polyazeotropy remains challenging and prior modeling studies often failed to satisfactorily correlate the polyazeotropic behavior of associating systems. Based on the association NRTL model, this study presents successful correlation and analysis of double azeotrope formation for three diverse systems: strongly associating methanol + diethylamine binary, weakly associating ammonia + pentafluoroethane binary, and non-associating benzene + hexafluorobenzene binary systems. Furthermore, these systems exhibit polyazeotropic behavior as the temperature increases and approaches their Bancroft points.

42 ENGINEERING↗
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 OSTI2022

Altering the Supply of Shielding Gases to Fabricate Distinct Geometry in GMA Additive Manufacturing

Wire arc additive manufacturing (WAAM) is the process by which large, metallic structures are built, layer-by-layer, using a welding arc to melt wire feedstock. In this process, the proper selection of the shielding gas plays a vital role in the achievement of structurally acceptable part geometries and quality surface finishes. In this study, the authors used either a ternary mix (He, Ar and CO 2 ) or a binary mix (Ar and CO 2 ) of shielding gases to deposit wall geometries using an open loop-controlled WAAM system developed at Oak Ridge National Laboratory’s Manufacturing Demonstration Facility. The binary blend produced a wider and shorter geometry, while the ternary blend resulted in a narrower build that was more equivalent to the CAD geometry. The data indicated that the binary blend provided a higher oxygen concentration in the weld as compared to that of the ternary blend. The results imply that the arc characteristics and heat input had a significantly higher impact on the weld penetration than the surface tension effect of surface active elements. This was further verified by developing and applying a high-fidelity computational fluid dynamics (CFD) model of the thermophysical properties of gas mixtures. The results from the model showed that, while the influence of increased oxygen concentration on the surface tension for the binary blend led to a deeper penetration, the ternary blend gave rise to heat flux to the workpiece.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Data-Driven Many-Body Models for Molecular Fluids: CO 2 /H 2 O Mixtures as a Case Study

In this study, we extend the scope of the many-body TTM-nrg and MB-nrg potential energy functions (PEFs), originally introduced for halide ion–water and alkali-metal ion–water interactions, to the modeling of carbon dioxide (CO2) and water (H2O) mixtures as prototypical examples of molecular fluids. Both TTM-nrg and MB-nrg PEFs are derived entirely from electronic structure data obtained at the coupled cluster level of theory and are, by construction, compatible with MB-pol, a many-body PEF that has been shown to accurately reproduce the properties of water. Although both TTM-nrg and MB-nrg PEFs adopt the same functional forms for describing permanent electrostatics, polarization, and dispersion, they differ in the representation of short-range contributions, with the TTM-nrg PEFs relying on conventional Born–Mayer expressions and the MB-nrg PEFs employing multidimensional permutationally invariant polynomials. By providing a physically correct description of many-body effects at both short and long ranges, the MB-nrg PEFs are shown to quantitatively represent the global potential energy surfaces of the CO 2 –CO 2 and CO 2 –H 2 O dimers and the energetics of small clusters, as well as to correctly reproduce various properties in both gas and liquid phases. Furthermore, building upon previous studies of aqueous systems, our analysis provides further evidence for the accuracy and efficiency of the MB-nrg framework in representing molecular interactions in fluid mixtures at different temperature and pressure conditions.

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

Uncertain characterization of reservoir fluids due to brittleness of equation of state regression

Equations of state (EoS) play a central role in modeling the phase equilibrium of fluid mixtures. Their parameterization involves fitting a model to experimental data, i.e., solving a nonlinear, non-convex, multivariate optimization problem. The latter requires one to select design variables, domains of definition for each variable, and weights assigned to individual measurements. We demonstrate that subjective choices of an optimization algorithm and an initial guess also impact the regression process. Consequently, EoS predictions are fundamentally uncertain even after the EoS tuning to a limited set of experimental data points. We demonstrate this observation for two hydrocarbon reservoir fluids, in which five properties of the heaviest carbon fraction are treated as design variables. While all the optimization algorithms and initial guesses match experimental data for the gas and liquid properties, the resulting EoS parameterizations lead to dramatically different predictions of the fluid’s thermophysical behavior in the unsampled pressure and temperature regions. In conclusion, we propose the probabilistic treatment of design variables to quantify the predictive uncertainty of the resulting fluid models.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Rapid Coal-Ash Characterization using Geophysical Methods & Machine Learning

Coal combustion products (CCP) are challenging to delineate in heterogeneous field settings. Conventional methods (test pits, coring, and laboratory analyses) are labor-intensive, slow, invasive, and provide sparse spatial coverage. This study evaluates whether rapid non-invasive geophysical screening methods—induced polarization (IP), magnetic susceptibility, and nuclear magnetic resonance (NMR) —combined with surface colorimetry (RGB_24), can discriminate CCP-soil mixtures and provide reliable estimates of CCP content. Laboratory measurements were collected on five CCP-soil mixtures (series) and modeled using (i) a linear baseline, (ii) a calibrated non-linear (power-mean) model, and (iii) a machine-learning (ML) Random Forest approach, with validation via leave-one-series-out and site-specific tests. Across the five series, individual signals—particularly IP and magnetic susceptibility—were strongly predictive of ash content but were consistently outperformed by combined models. The pooled calibrated non-linear and ML models captured the observed non-linearity and achieved high accuracy and precision, improving on linear fits. Colorimetry showed the weakest direct relationship with ash content for the tested samples but improved performance when included in multi-signal models. At pre-selected 3.5% decision threshold, calibrated and ML approaches yielded near-perfect classification (Matthews correlation coefficient ˜ 1), suggesting strong practical operability for field screening. Additionally, field-analog tests highlighted the role of endmembers—accuracy declined without access to end-member measurements but was largely recovered by collecting a minimal labeled pair for local recalibration. With end members, accuracy remained high. Globally trained models performed well on three operational unknowns; however, series-specific refits provided the most accurate predictions. Overall, these results highlight the potential of combining rapid geophysics and minimal local calibration for improved coal-ash delineation.

Peshtani, Klaudio↗
DOE OSTI2024

Assessment of simulated and observed cavitation-induced erosion damage in Spallation Neutron Source target vessels

Cavitation-induced erosion damage in different Spallation Neutron Source (SNS) target designs are simulated using explicit finite element–based techniques and compared with observations of erosion in targets after operation. The efficacy of the previously developed method, called saturation time, was evaluated using erosion-damaged samples from new target designs. A new metric called maximum bubble size was implemented under the rationale that larger cavitation bubbles will collapse more intensely. The maximum cavitation bubble size over 1 ms of simulated time was calculated based on the Rayleigh–Plesset equation for each element integration point and presented as a contour map at the vessel surface for assessing with erosion observations. SNS targets are now operated with helium gas injection to reduce cavitation damage. A simulation method using a material model for the mixture of mercury and gas bubbles was recently developed and used to account for the effect of small gas bubbles on the structural response of the target vessel. Furthermore, this work compares the new method's results with observed cavitation damage. Maps of the calculated maximum bubble size for targets operated with and without gas injection were compared with photographs of erosion damage observed in SNS targets. The patterns in maximum bubble size maps correlated well with observations of erosion patterns in target vessels after service. Advantages and challenges of the maximum bubble size simulation technique are provided, and differences between results from the previous and the newly proposed metric are discussed.

Jiang, Hao↗
DOE OSTI2020

Capturing Deuteration Effects in a Molecular Mechanics Force Field: Deuterated THF and the THF–Water Miscibility Gap

Deuteration is a common chemical modification used in conjunction with experiments such as neutron scattering, NMR, and Fourier-transform infrared for the study of molecular systems. Under the Born–Oppenheimer (BO) approximation, while the underlying potential energy surface remains unchanged by isotopic substitutions, isotopic substitution still alters intramolecular vibrations, which in turn may alter intermolecular interactions. Molecular mechanics (MM) force fields used in classical molecular dynamics simulations are assumed to represent local approximations of the BO potential energy surfaces, and hence, MD simulations using simple isotopic mass substitutions should capture BO-compatible isotope effects. However, standard MM force-field parameterizations do not directly fit to the local harmonic quantum mechanical (QM) Hessian that describes the BO surface, but rather to QM normal-modes and/or mass-dependent internal-coordinate derived distortion energies. In this study, using tetrahydrofuran (THF)–water mixtures as our model system, we show that not only does a simple mass-substitution approach fail to capture an experimentally characterized deuteration effect (the loss of the closed-loop miscibility gap associated with the complete deuteration of THF) but also it is necessary to generate new MM force-field parameters that correctly describe isotopic dependent vibrations to capture the experimental deuteration effect. We show that the origin of this failure is a result of using mass-dependent features to fit the THF MM force field, which unintentionally biases the bonded terms of the force field to represent only the isotopologue used during the original force-field parameterization. In addition, we make use of our isotopologue-corrected force field for D 8 THF to examine the molecular origins of the isotope-dependent loss of the THF–water miscibility gap.

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

Molecular-scale origins of wettability at petroleum–brine–carbonate interfaces

Wettability control of carbonates is a central concept for enhanced petroleum recovery, but a mechanistic understanding of the associated molecular-scale chemical processes remains unclear. We directly probe the interface of calcium carbonate (calcite) with natural petroleum oil, synthetic petroleum analogues, and aqueous brines to understand the molecular scale behavior at this interface. The calcite–petroleum interface structure is similar whether or not calcite was previously exposed to an aqueous brine, and is characterized by an adsorbed interfacial layer, significant structural changes within the calcite surface, and increased surface roughness. No evidence for an often-assumed thin-brine wetting layer at the calcite–petroleum interface is observed. These features differ from those observed at the calcite–brine interface, and for parallel measurements using model synthetic petroleum mixtures (consisting of representative components, including dodecane, toluene, and asphaltene). Changes to the interface after petroleum displacement by aqueous brines are also discussed.

02 PETROLEUM↗