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 505 records · Page 28

Transport in Stochastic Media with Random Chord Length Distributions

Thermal radiation transport computations in binary Markovian random mixtures rely almost exclusively on the Levermore-Pomraning (LP) model which is obtained by applying a heuristic closure to the ensemble averaged random medium transport equation. The validity of this model has been extensively tested by comparing numerical results over a broad parameter range (material types and mixing parameters) against benchmark solutions in planar geometry. The conditions under which the LP-model provides useful results and when it breaks down are now well established, but work to date has been largely restricted to homogeneous mixing statistics, i.e., the mean chord lengths of both materials are taken to be spatially constant. In recent work, this limitation was relaxed by allowing the mean chord lengths and, in a consistent fashion, the volume fractions in the LP-model to vary spatially and in a follow-up investigation benchmark solutions were obtained by ensemble averaging results over material realizations sampled from a nonhomogeneous Poisson process (NHPP). Numerical experiments in rod geometry with specifically linear and quadratic spatial dependence of chord lengths showed that the material averaged radiation intensities vary nonmonotonically with depth into the medium, in stark contrast to solutions obtained assuming uniform chord lengths. Moreover, depending on the local optical thickness and strength of scattering, the LP-model results showed locally more nuanced deviations from the benchmark solutions than was the case with constant chord lengths. These limited numerical investigations highlight the nontrivial qualitative and quantitative consequences of nonhomogeneous mixing statistics, in particular that closure approximations may not be uniformly valid or invalid over the problem domain.

42 ENGINEERING↗

System Resilience Benefits of Dual-Fuel Capable Generators

The growing dependency on natural gas (NG) for power generation raises challenges for ensuring the resilience of power systems during extreme cold weather. Dual-fuel capable generators that can switch from burning NG to distillate fuel oil during an NG shortage over one way to alleviate these challenges. In this study, the impacts of gas unavailability on the IEEE 73-bus reliability test system (RTS) with 2019 updated generation mixture are explored. An extension to the standard production cost model formulation of unit commitment and economic dispatch is proposed to allow the representation of dual-fuel capable generators that can switch fuels between NG and distillate oil with a specified oil tank capacity and tank refueling capability. The operation of the system under gas unavailability with 0%, 25%, 50%, 75%, and 100% of combined cycle and combustion turbine generators as dual-fuel capable with a one-day supply of fuel oil are simulated using PLEXOS, a production cost model. The dual-fuel generator performance, unserved energy, and system costs are fully assessed.

dual-fuel capable generator↗

Direct Conversion of Delignified Biomass into Non-Volatile Products Using Thermophilic Bacteria

Butyric acid (BA) can be used as a platform intermediate for diesel and sustainable aviation fuels (SAFs) as well as a precursor for diverse commodity chemicals for the production of polymers, fibers, solvents, and preservatives. This work presents a co-culture-based bioprocess for the production of BA from corn stover by two thermophilic bacteria, Clostridium thermocellum, a well-known efficient degrader of insoluble and oligomeric cellulosic substrates, and Clostridium thermobutyricum, a highly efficient BA producer from monomeric sugars. First, we tested the capability of C. thermocellum to deconstruct three different biomass streams from the NREL pilot plant pretreatment process -raw, deacetylated, and deacetylated and mechanically refined corn stover (DMR)-, and up to 82% solids deconstruction and 91% carbohydrate utilization were obtained when using DMR solids. Then, the sugar utilization and by-product formation capabilities on various substrates by C. thermobutyricum were compared to evaluate its potential as a co-culture partner for C. thermocellum. Model carbon sources -glucose and xylose- were compared to a mixture of substrates present in soluble DMR hydrolysate streams, and the highest cell density as well as the maximum BA titer (8.3 g/L) were observed on DMR hydrolysate. Finally, clostridial co-cultivations on DMR solids were performed to investigate the inoculation time for both organisms that enhances BA production; while C. thermocellum was always inoculated at the beginning of the cultivation, C. thermobutyricum was inoculated at 0, 36, or 72h. It was found that BA production begins immediately upon inoculation of C. thermobutyricum regardless of the inoculation time. However, higher solids deconstruction and carbohydrate utilization results (83.7 and 95.6%, respectively) were obtained when both organisms were inoculated at the beginning of the fermentation. Under those conditions, up to 2.6 g/L of BA were obtained from DMR solids after 120h of cultivation. This work provides relevant information for the sustainable production of BA using thermophilic bacteria.

bacterial co-culture↗

Measurements of He-Ar, He-N 2 , and He-air thermal creep slip and accommodation coefficients at high temperatures

Understanding gas flows in capillaries has many applications in modeling the transport of gases in nano-structured, porous, or fractured media. A network of capillaries can often approximate these media, and also information on gas-surface interactions obtained from capillary experiments can be used for modeling flows in these media. Experimental data on flows of different mixtures of He-Ar, He-N 2 , and He-air in the slip regime and in stainless steel capillaries at high temperatures were obtained by using a two-bulb apparatus. An accurate expression for the thermal creep slip coefficients with the Lennard-Jones potential parameters and diffuse-specular reflection gas-surface interaction conditions were then used to obtain the accommodation coefficients for the different gases. The experimental data are best described with values of accommodation coefficients in the range of 0.1–0.3 for He and 0.5–1.0 for Ar, N 2 , and air. Further, the use of values in this range is suggested for modeling gaseous flows in capillaries and nano-structured, porous, or fractured media if other direct measured values for a particular medium are unavailable.

42 ENGINEERING↗

Ignition Delay Times and Chemical Kinetic Model Validation for Hydrogen and Ammonia Blending With Natural Gas at Gas Turbine Relevant Conditions

Ignition delay times from undiluted mixtures of natural gas (NG)/H 2 /Air and NG/NH 3 /Air were measured using a high-pressure shock tube at the University of Central Florida. The combustion temperatures were experimentally tested between 1000 and 1500 K near a constant pressure of 25 bar. As mentioned, mixtures were kept undiluted to replicate the same chemistry pathways seen in gas turbine combustion chambers. Recorded combustion pressures exceeded 200 bar due to the large energy release, hence why these were performed at the high-pressure shock tube facility. The data are compared to the predictions of the NUIGMech 1.1 mechanism for chemical kinetic model validation and refinement. An exceptional agreement was shown for stoichiometric conditions in all cases but strayed at lean and rich equivalence ratios, especially in the lower temperature regime of H 2 addition and all temperature ranges of the baseline NG mixture. Hydrogen addition also decreased ignition delay times by nearly 90%, while NH 3 fuel addition made no noticeable difference in ignition time. NG/NH 3 exhibited similar chemistry to pure NG under the same conditions, which is shown in a sensitivity analysis. Here, the reaction CH 3 + O 2 = CH 3 O + O is identified and suggested as a possible modification target to improve model performance. Increasing the robustness of chemical kinetic models via experimental validation will directly aid in designing next-generation combustion chambers for use in gas turbines, which in turn will greatly lower global emissions and reduce greenhouse effects.

33 ADVANCED PROPULSION SYSTEMS↗

Model-based compositional predictions for a differential scanning calorimetry/thermogravimetric analysis-mass spectrometry system used for heat of vaporization measurements

A Differential Scanning Calorimetry/Thermogravimetric Analysis (DSC/TGA) procedure for determination of Heat of Vaporization (HOV) of fuel and surrogate fuel mixtures has been previously utilized for evaluation of oxygenate impacts on spark ignition fuels. This analysis has been further leveraged with mass spectrometry (MS) to elucidate vapor phase compositional changes during evaporation. In this study we focus on a simple system of three compounds (ethanol, n-hexane, and n-octane) to provide a model-based understanding of the evaporation process under these experimental conditions. An Aspen Plus model was set up to provide liquid and vapor phase concentration predictions during mixture evaporation using the initial mixture composition, and measured temperature and mass loss profiles from DSC/TGA experiments. Agreements between the trends of the MS ion count profiles and model composition predictions provide a degree of validation of this modeling approach that can allow the interpretation and extrapolation of DSC/TGA measurements for more complex fuel mixtures.

33 ADVANCED PROPULSION SYSTEMS↗

Catalytic Triad-Inspired Nanozyme Catalysts for Ester Hydrolysis in Organic Solvent Mixtures

We synthesized a series of molecular nanozyme catalysts containing functional groups (imidazole, carboxylic acid, and hydroxyl groups) inspired by the catalytic triad found in natural serine hydrolases. Different structural features were incorporated using two distinct synthesis routes to investigate the influence of interactions beyond the active site on catalytic activity in designing molecular nanozymes. Molecular dynamics simulations suggest that nanozyme activity is affected by structural features that influence nanozyme hydrophilicity and the organization of the local solvent environment. The most active nanozyme showed activity comparable to the enzyme α-chymotrypsin for the hydrolysis of a model ester (4-nitrophenyl 4- hydroxybenzoate) in a 95:5 (v/v%) water/acetonitrile mixture at ambient temperature. Increasing the temperature and organic solvent content decreases the activity of α-chymotrypsin while enhancing the activity of the nanozymes. The nanozymes can be immobilized on supported metal nanoparticles using a dithiol self-assembled monolayer, which facilitates their removal from the postreaction solution. These results demonstrate the potential in creating solvent-tolerant bioinspired catalysts, thereby combining the advantages of biocatalysts and chemical catalysts as next-generation industrial catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CFD modeling of near-wall combustion and unburned methane prediction in natural gas spark ignition engines

Natural gas-powered engines play a critical role in gas drilling, compression, and transmission sectors, but methane (CH 4 ) from engine combustion slip can be significant over their lifespan, contributing to atmospheric pollution and signaling reduced engine efficiency. Here, to address this challenge, computational fluid dynamics (CFD) simulations offer valuable insights into the in-cylinder combustion process, enabling the optimization of combustion strategies and engine designs to minimize unburned CH 4 slip. This study aims to evaluate and improve combustion models for simulating the combustion process and predicting unburned CH 4 concentrations in natural gas spark-ignition (SI) engines, including engines that are part of combined reformer-engine systems. Specifically, the performance of two flamelet-based combustion models—the Extended Coherent Flame Model (ECFM) and the G-equation model—was assessed using experimental engine data collected under varying excess-air ratio (λ) conditions and fuel compositions, including natural gas and syngas blends. In addition, to enhance the predictive capabilities of the G-equation model, a flame-wall interaction (FWI) sub-model was integrated into its framework. The effects of its model parameters, such as quenching and influence distance, on combustion behavior and unburned methane predictions were analyzed in detail. The ECFM tended to predict delayed combustion phasing under diluted mixture conditions, resulting in overprediction of unburned CH 4 concentrations. In contrast, the G-equation model provided reasonable predictions of combustion pressure, while representing higher the CH 4 reduction rate across the operating condition compared to experimental data. Incorporating the FWI sub-model—with the quenching distance calculated based on a pressure-dependent relation (P -0.48 ) and a fixed influence distance of 1.5 mm—further improved the G-equation model’s accuracy in predicting CH 4 reduction rates without compromising its ability to simulate the combustion process.

Combustion model↗

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↗

Large-Eddy Simulation of Laser-Ignited Direct Injection Gasoline Spray for Emission Control

Large-Eddy Simulations (LES) of a gasoline spray, where the mixture was ignited rapidly during or after injection, were performed in comparison to a previous experimental study with quantitative flame motion and soot formation data [SAE 2020-01-0291] and an accompanying Reynolds-Averaged Navier–Stokes (RANS) simulation at the same conditions. The present study reveals major shortcomings in common RANS combustion modeling practices that are significantly improved using LES at the conditions of the study, specifically for the phenomenon of rapid ignition in the highly turbulent, stratified mixture. At different ignition timings, benchmarks for the study include spray mixing and evaporation, flame propagation after ignition, and soot formation in rich mixtures. A comparison of the simulations and the experiments showed that the LES with Dynamic Structure turbulence were able to capture correctly the liquid penetration length, and to some extent, spray collapse demonstrated in the experiments. For early and intermediate ignition timings, the LES showed excellent agreement to the measurements in terms of flame structure, extent of flame penetration, and heat-release rate. However, RANS simulations (employing the common G-equation or well-stirred reactor) showed much too rapid flame spread and heat release, with connections to the predicted turbulent kinetic energy. With confidence in the LES for predicted mixture and flame motion, the predicted soot formation/oxidation was also compared to the experiments. The soot location was well captured in the LES, but the soot mass was largely underestimated using the empirical Hiroyasu model. An analysis of the predicted fuel–air mixture was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing.

42 ENGINEERING↗

Characterizing the contribution of bioaerosols diversity from complex aerosol particle samples.

Current global change scenarios expect to significantly increase the contribution of pollen to the total organic aerosol budget. In addition, pollen can burst into smaller fragments that are highly efficient as cloud condensation nuclei. However, current climate models still do not incorporate the effects of these particles due to detection and quantification challenges in complex aerosol mixtures. Chemical biomarkers (i.e. fructose) have commonly been used to trace pollen but this approximation can generate false positive results given the large number of shared compounds between different bioaerosols. A previous pilot EMSL experiment performed on pure pollen samples suggested that using a set of metabolites as a fingerprint of bioaerosols could serve to substantially increase detection accuracy over single biomarker approaches. In this proposal we developed novel bioinformatic strategies to characterize and quantify complex pollen mixtures present in the atmosphere. For that, we used diverse machine learning algorithms to analyze metabolic fingerprints from complex mixtures of different pollen acquired with high-resolution mass spectrometry (Orbitrap Q-Exactive). Our preliminary results demonstrated the potential to accurately discern and identify at specific relative abundances different pollen species in complex mixtures. Our research will lead to a significant advancement in the atmospheric chemistry and provide much needed data to improve the accuracy of climate models.

54 ENVIRONMENTAL SCIENCES↗

Factorization of Binary Matrices: Rank Relations, Uniqueness and Model Selection of Boolean Decomposition

The application of binary matrices are numerous. Representing a matrix as a mixture of a small collection of latent vectors via low-rank decomposition is often seen as an advantageous method to interpret and analyze data. In this work, we examine the factorizations of binary matrices using standard arithmetic (real and nonnegative) and logical operations (Boolean and $\mathbb{Z}$ 2 ). We examine the relationships between the different ranks, and discuss when factorization is unique. In particular, we characterize when a Boolean factorization X = W$\land$H has a unique W, a unique H (for a fixed W), and when both W and H are unique, given a rank constraint. We introduce a method for robust Boolean model selection, called BMFk, and show on numerical examples that BMFk not only accurately determines the correct number of Boolean latent features but reconstruct the pre-determined factors accurately.

97 MATHEMATICS AND COMPUTING↗

Bridging confined phase behavior of CH 4 -CO 2 binary systems across scales

Phase behavior of confined fluids may deviate significantly from that of the bulk fluid due to the fluid-wall interactions being a significant portion of all intermolecular interactions under confinement. Despite recent advancements in understanding confined phase behavior of pure fluids, confined phase behavior of mixtures remains an understudied topic. In this work, we examine the confined phase behavior of a CH 4 -CO 2 binary system by combining Monte Carlo (MC) simulations, a cubic equation of state (EoS), and the lattice Boltzmann method (LBM). First, the effects of confinement on density and phase distribution in nano-size pores are established using Gibbs Ensemble MC calculations, which produce precise results of liquid and vapor confined pressures and account for the modification of the phase change location. By comparing the phase envelopes of bulk and confined mixtures at a fixed temperature, here it is observed that the phase envelopes shrink with reductions in pore size. Based on this observation, we extend a modified Peng-Robinson EoS, which was originally developed for pure fluids under confinement, to mixtures via van-der-Waals-type mixing rules and by accounting for shifts in the critical properties of confined CH 4 -CO 2 . The resulting phase envelopes are in good agreement with the MC data. In addition, a local density model is used in combination with the confined EoS to calculate adsorption isotherms of CH 4 -CO 2 mixtures and to characterize the behavior of confined matter in nanopores. Finally, we incorporate this EoS in a multicomponent multiphase LBM that uses a pseudopotential model to represent intermolecular forces. This workflow utilizes multiscale simulation techniques to bridge the behavior of multicomponent systems across scales and to shed light on the confined phase behavior of CH 4 -CO 2 binary systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of radiation heat loss on planar and spherical hydrofluorocarbon/air flames

Hydrofluorocarbons (HFC), which are mildly flammable and pose potential fire risks, have received greater attention as a viable low global warming potential alternative to traditional refrigerant and fire-suppressant compounds. Therefore, there is a demand to accurately quantify their flammability and re-activity to establish proper safety metrics. Here, this study investigates the effects of radiation heat loss on slowly-propagating HFC/air laminar flames. Planar 1-D simulations of R-32/air and R-1234yf/air flames show significant reductions in laminar flame speed due to radiative heat losses from the flame zone. Simulations of spherically expanding flames (SEF) revealed that the radiation-induced flow needs to be considered when interpreting data from experiments. To this end, a Spherical-flame RADiation-Induced Flow (SRADIF) model was developed to estimate the burned gas inward flow velocities in constant-pressure SEFs, utilizing the optically thin limit assumption to model radiation heat loss. The model was validated against results from detailed numerical simulations of SEFs, from which radiation-induced inward flow was derived using a new formulation considering both the radiation heat loss and convective flow effects. Results show that SRADIF accurately predicts the inward flow velocity for R-32/air mixtures over a range of conditions and performs significantly better compared to existing analytical models. However, the model was unable to accurately predict flow velocities for R-1234yf/air flames and the reason for this is discussed. (224 words)

42 ENGINEERING↗

Bioreactor Optimization through Multi-Phase Flow Models (CRADA Final Report)

Chemical manufacturing uses 29% of energy in the United States and produces 925 million metric tons of CO2 annually. Biomanufacturing offers the potential to leverage America’s rich agricultural resources to produce critical chemicals such as lubricants, pharmaceutical precursors, and components of energetic materials that today are sourced extensively from overseas. The bioreactors used in biomanufacturing applications, such as one developed by Capra Biosciences, involve multiphase flow of biofilm-coated solid support particles that are continuously circulated in a fluidized state within the reactor along with a constant supply of oxygen via an aeration mechanism. In this project, Capra Biosciences and LBNL developed a multiscale modeling framework to simulate the multiphase flows of solid particles in a liquid-gas bubble mixture that occurs in the bioreactor using the current MFIX-Exa software, an opensource multiphase flow solver developed and maintained at LBNL and NETL. By leveraging HPC capabilities, this high-fidelity multiscale model was used to inform design decisions for bioreactor architecture to make them operationally efficient.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

N 2 vibrational excitation in atmospheric pressure ns pulse and RF plasma jets

Time-resolved N 2 vibrational temperature and translational–rotational temperature in quasi-two-dimensional atmospheric pressure plasma jets sustained by ns pulse and RF discharges in nitrogen/noble gas mixtures are measured by the broadband vibrational Coherent Anti-Stokes Raman Scattering (CARS) . The results indicate a much stronger vibrational excitation in the RF plasma jet, due to the lower reduced electric field and higher discharge power. In a ns pulse discharge in N 2 /He, N 2 vibrational temperature is significantly lower compared to that in N 2 /Ar, due to the more rapid vibration–translation (V–T) relaxation of nitrogen by helium atoms. In the RF plasma jets in N 2 /Ne and N 2 /Ar, the vibrational excitation increases considerably as the nitrogen fraction in the mixture is reduced. The experimental data in the RF plasma jet in N 2 /Ar jet are compared with the kinetic modeling predictions. The results indicate that nitrogen vibrational excitation in N 2 /Ar plasma jets with a small N 2 fraction in the mixture (several percent) is controlled primarily by electron impact, anharmonic vibration–vibration (V–V) pumping, and V–T relaxation by N atoms. In comparison, V–V energy transfer from the vibrationally excited molecules in the first excited electronic state, N 2 (A 3 Σ u + , v), which are generated primarily by the energy transfer from the metastable Ar atoms, has a minor effect on the vibrational populations of the ground electronic state, N 2 (X 1 Σ g + , v). Although the discharge energy fraction going to electronic excitation is significant, the predicted quasi-steady-state N 2 (A 3 Σ u + ) number density, controlled by the energy pooling and quenching by N atoms, remains relatively low. Because of this, the net rate of N 2 (X 1 Σ g + ) vibrational excitation by the V–V energy transfer from N 2 (A 3 Σ u + ) is much lower compared to that by the direct electron impact. Furthermore, the results show that atmospheric pressure RF plasma jets can be used as sources of highly vibrationally excited N 2 molecules and N atoms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mixed Material Scintillator Systems Particle ID Modelling Report

AMixed-Material Scintillator System (MMSS) is a radiation detector using a scintillator made of a heterogeneous, structured mixture of two or more materials. MMSSs are designed so that the structured mixing of materials encodes properties of the radiation detected in the scintillator. This allows for new radiation detectors with advantages over detectors using traditional homogeneous scintillators. This report contains results from this project’s second task, to model and quantify the promise of the MMSSs invented earlier in the project. Those inventions have been divided into two categories, and this report focuses on the first category: MMSS inventions the use zones of differently-colored scintillators to detect and characterize neutron sources. We call this category particle-ID, or PID, detectors. Inventions using gradients designed to encode the position of radiation interactions using gradients will be discussed in a future report.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Corrosion of Rare-Earth Disilicate Environmental Barrier Coatings in Gas Turbine Environments—Part 2: Microstructure Model

A model for the volatization of ytterbium disilicate (YbDS) by water vapor is developed. The new model is based on the Deal & Grove theory of the thermal oxidation of silicon [1], generalized here to account for the multi-step, second-order chemical reactions for silica depletion, multi-component mixture effects, and the unsteady porous microstructural evolution of the corroded inner zone. With the new model, the fundamental kinetic rate and standard porous media parameters are determined using the atmospheric steam jet impingement experiment of Ridley & Opila [2] and compared to results from the original Deal & Grove model. Fixing the input parameters in each respective model, a validation study (blind prediction) is presented for the high Reynolds number (1 × 105) and high pressure (14.7 atm) burner rig experiment of Wan et al. [3].

Cha, Chong M [Rolls Royce]↗