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.

SourceResearch connectionFrom yearThrough year

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

At least 1,387 records · Page 77

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Organic synthesis in experimental impact shocks

Laboratory simulations of shocks created with a high-energy laser demonstrate that the efficacy of organic production depends on the molecular, not just the elemental composition of the shocked gas. In a methane-rich mixture that simulates a low-temperature equilibrium mixture of cometary material, hydrogen cyanide and acetylene were produced with yields of 5 x 10(17) molecules per joule. Repeated shocking of the methane-rich mixture produced amine groups, suggesting the possible synthesis of amino acids. No organic molecules were produced in a carbon dioxide-rich mixture, which is at odds with thermodynamic equilibrium approaches to shock chemistry and has implications for the modeling of shock-produced organic molecules on early Earth.

NASA Discipline Exobiology↗

NASA NTRS2007

One-Dimensional Ablation with Pyrolysis Gas Flow Using a Full Newton's Method and Finite Control Volume Procedure

The development and verification of a one-dimensional material thermal response code with ablation is presented. The implicit time integrator, control volume finite element spatial discretization, and Newton's method for nonlinear iteration on the entire system of residual equations have been implemented and verified for the thermochemical ablation of internally decomposing materials. This study is a continuation of the work presented in "One-Dimensional Ablation with Pyrolysis Gas Flow Using a Full Newton's Method and Finite Control Volume Procedure" (AIAA-2006-2910), which described the derivation, implementation, and verification of the constant density solid energy equation terms and boundary conditions. The present study extends the model to decomposing materials including decomposition kinetics, pyrolysis gas flow through the porous char layer, and a mixture (solid and gas) energy equation. Verification results are presented for the thermochemical ablation of a carbon-phenolic ablator which involves the solution of the entire system of governing equations.

Amar, Adam J.↗

NASA NTRS2018

Scattering and Radiative Properties of Morphologically Complex Carbonaceous Aerosols: A Systematic Modeling Study

This paper provides a thorough modeling-based overview of the scattering and radiative properties of a wide variety of morphologically complex carbonaceous aerosols. Using the numerically-exact superposition T-matrix method, we examine the absorption enhancement, absorption Angstroem exponent (AAE), backscattering linear depolarization ratio (LDR), and scattering matrix elements of black-carbon aerosols with 11 different model morphologies ranging from bare soot to completely embedded soot-sulfate and soot-brown carbon mixtures. Our size-averaged results show that fluffy soot particles absorb more light than compact bare-soot clusters. For the same amount of absorbing material, the absorption cross section of internally mixed soot can be more than twice that of bare soot. Absorption increases as soot accumulates more coating material and can become saturated. The absorption enhancement is affected by particle size, morphology, wavelength, and the amount of coating. We refute the conventional belief that all carbonaceous aerosols have AAEs close to 1.0. Although LDRs caused by bare soot and certain carbonaceous particles are rather weak, LDRs generated by other soot-containing aerosols can reproduce strong depolarization measured by Burton et al. for aged smoke. We demonstrate that multi-wavelength LDR measurements can be used to identify the presence of morphologically complex carbonaceous particles, although additional observations can be needed for full characterization. Our results show that optical constants of the host/coating material can significantly influence the scattering and absorption properties of soot-containing aerosols to the extent of changing the sign of linear polarization. We conclude that for an accurate estimate of black-carbon radiative forcing, one must take into account the complex morphologies of carbonaceous aerosols in remote sensing studies as well as in atmospheric radiation computations.

carbonaceous aerosols; scattering matrix; polariza↗

NASA NTRSDate not supplied

3D Phase-Field Simulations of Pattern Formation During Freeze Casting

We present the results of a combined experimental and phase-field modeling study of pattern formation during freeze casting. Those studies use unidirectional freezing of simple binary liquid mixtures of water and sugars (sucrose and trehalose) in a temperature gradient, which suffice to produce hierarchical templated structures similar to those observed in more complex multi-component freeze-cast systems including lamellae, undulated ridges, and more exotic “jellyfish-like” substructures. Multiscale 3D phase-field simulations reproduce remarkably well those structures quantitatively and identify key properties of the ice-water interface that control their formation. They further reveal that lamellae form as a result of a novel symmetry-breaking secondary instability of partially faceted cellular structures and pinpoint additional secondary instability mechanisms giving rise to smaller-scale substructures.

Kaihua Ji↗

NASA NTRSDate not supplied

3D Phase-Field Simulations of Pattern Formation During Freeze Casting

Kaihua Ji↗

NASA NTRS2022

A Comprehensive Machine Learning Study to Classify Precipitation Type over Land from Global Precipitation Measurement Microwave Imager (GPM-GMI) Measurements

Precipitation type is a key parameter used for better retrieval of precipitation characteristics as well as to understand the cloud–convection–precipitation coupling processes. Ice crystals and water droplets inherently exhibit different characteristics in different precipitation regimes (e.g., convection, stratiform), which reflect on satellite remote sensing measurements that help us distinguish them. The Global Precipitation Measurement (GPM) Core Observatory’s microwave imager (GMI) and dual-frequency precipitation radar (DPR) together provide ample information on global precipitation characteristics. As an active sensor, the DPR provides an accurate precipitation type assignment, while passive sensors such as the GMI are traditionally only used for empirical understanding of precipitation regimes. Using collocated precipitation type flags from the DPR as the “truth”, this paper employs machine learning (ML) models to train and test the predictability and accuracy of using passive GMI-only observations together with ancillary information from a reanalysis and GMI surface emissivity retrieval products. Out of six ML models, four simple ones (support vector machine, neural network, random forest, and gradient boosting) and the 1-D convolutional neural network (CNN) model are identified to produce 90–94% prediction accuracy globally for five types of precipitation (convective, stratiform, mixture, no precipitation, and other precipitation), which is much more robust than previous similar effort. One novelty of this work is to introduce data augmentation (subsampling and bootstrapping) to handle extremely unbalanced samples in each category. A careful evaluation of the impact matrices demonstrates that the polarization difference (PD), brightness temperature (Tc) and surface emissivity at high-frequency channels dominate the decision process, which is consistent with the physical understanding of polarized microwave radiative transfer over different surface types, as well as in snow and liquid clouds with different microphysical properties. Furthermore, the view-angle dependency artifact that the DPR’s precipitation flag bears with does not propagate into the conical-viewing GMI retrievals. This work provides a new and promising way for future physics-based ML retrieval algorithm development.

machine learning/artificial intelligence↗

DOE OSTI2021

A Simple Criterion for Feasibility of Heat Integration between Distillation Streams Based on Relative Volatilities

In a multicomponent distillation configuration, there are numerous sources and sinks of heat, and a potential way to reduce the heat duty requirement is to perform heat integration. Unfortunately, an algorithmic search of the optimal heat integration opportunities is intractable when the required temperatures of intermediate streams are computed via complex models. Instead, in this work, we introduce pressure-scaled pseudo relative volatility, a new metric to compare stream temperatures. Here, we justify the use of pseudo relative volatility by proving that this variable is a monotonically increasing function of the liquid fraction in a saturated mixture stream. Using this metric, we derive a shortcut criterion to check the feasibility of various heat integration opportunities, such as thermal coupling via heat transfer (TCH). The advantage of this approach is that it circumvents the need for explicit temperatures and instead relies on composition, component relative volatilities, and pressure—quantities that are readily available in shortcut models for optimization of distillation configurations. Leveraging this fact, we propose a new optimization framework that identifies feasible TCHs that we consider within the formulation while minimizing the total heat duty of a distillation configuration. We demonstrate, on a few examples, that our formulation can identify heat duty efficient configurations, some of which are multieffect configurations. Using this methodology, we discover configurations that are not only simpler than the fully thermally coupled (FTC) configuration but also have a much lower heat duty.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1990

Modeling supersonic combustion using a fully-implicit numerical method

A fully-implicit finite-volume algorithm for two-dimensional axisymmetric flows has been coupled to a detailed hydrogen-air reaction mechanism (13 species and 33 reactions) so that supersonic combustion phenomena may be investigated. Numerical computations are compared with ballistic-range shadowgraphs of Lehr (1972) that exhibit two discontinuities caused by a blunt body as it passes through a premixed stoichiometric hydrogen-air mixture. The suitability of the numerical procedure for simulating these double-front flows is shown. The requirements for the physical formulation and the numerical modeling of these flowfields are discussed. Finally, the sensitivity of these external flowfields to changes in certain key reaction rate constants is examined.

Maccormack, Robert W.↗

DOE OSTI2026

AGR-5/6/7 Thermal Model with Non-uniform Gas Gaps

The daily as-run thermal analysis for the Advanced Gas Reactor (AGR)-5/6/7 experiment was documented in ECAR-5633 (Hawkes 2021). A finite element model was created for the entire test train, which consisted of five capsules. The fuel compacts, graphite holders, stainless-steel capsule walls, and all other major components were individually modeled, including all thermocouples (TCs). Daily heat rates for each compact and component in the model were input from daily as-run physics analyses. Daily gas compositions for each capsule were also input. The thermal conductivity of the compacts and graphite holders varied with fast neutron fluence. Gas mixture thermal conductivity was implemented using experimentally attained values from the literature. Fluence and temperature-dependent thermal conductivity was used for the graphite components and the fuel compacts. Radiation heat transfer was implemented. The model was tuned in an attempt to match the TC readings during the first cycle by adjusting the Neolube (graphite lubricant) thickness.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1991

Regenerator performance in a Vuilleumier refrigerator compared with a third-order numerical model

A three-stage Vuilleumier refrigerator was used to measure the performance of various third stage regenerators. The refrigerator operates between 2.5 and 5.0 Hz and, depending on the material used in the third stage regenerator, achieves temperatures of 8 to 20 K at the cold end of the third stage. This paper presents a comparison of regenerator performance for four regenerator materials: 229 micron diameter spheres of Pb(+)5 pct Sb, 229 micron diameter spheres of brass, 216 micron irregularly-shaped GdRh powder, and a mixture of 229 microns and 762 microns diameter spheres of Pb(+5) pct Sb. The experimental results are compared with a first-order model that neglects the void volume within the regenerator and with a third-order model that considers the effect of pressure oscillations in the regenerator void volume. Experimental results indicate that regenerator losses are dominated by the pressure oscillation in the void volume rather than the mass flow through the temperature gradient in the regenerator. These results are consistent with the third-order numerical model. This model shows that the heat capacity of the gas in the void space as well as the heat capacity of the matrix influences the regenerator performance.

Bradley, P. E.↗

NASA NTRS1990

Ice and minerals on Callisto - A reassessment of the reflectance spectra

The results of a comparison of laboratory spectral reflectance measurements of particulate mixtures of both hydrated silicates and palagonite with water ice, on the one hand, with two previously unpublished reflectance spectra of Callisto, yield direct support for the hypothesis that the measured reflectance of Callisto includes a substantial nonice component. Hapke's (1981) equations are used in a theoretical model for thorough analysis of telescopic data; a comparison of the calculation results thus obtained with measured reflectance data for Callisto indicate that three-component, ice/magnetite/serpentine mixtures are a better match for telescopic data than two-component ice mixtures.

Roush, T. L.↗

DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.'s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U parameter for the 5f electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl3) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

On Improving the Dust Lidar Ratios for Version 5 of CALIOP by Using Measurements with Source-tagged Model Results

Vertically resolved retrievals of aerosol backscatter and extinction have been available since 2006 from global measurements acquired by the space borne lidar CALIOP onboard the CALIPSO satellite. These retrievals depend crucially upon the ratio of extinction to backscatter (i.e., lidar ratio) which typically has to be assumed for the various aerosol species for an elastic backscatter lidar like CALIOP. For all CALIPSO data releases up until now (version 4.2), these lidar ratios have been taken as constant both spatially and temporally. However, it is well established that the aerosol lidar ratios can vary significantly regionally, seasonally, and also as a function of age and relative humidity. In particular, the dust lidar ratios show significant variation depending upon their source regions, which is related to the mineral composition of dust in that region. In this paper, we shall present a scheme to obtain the regionally and seasonally dependent lidar ratios of dust. In this scheme, which may be implemented in a forthcoming version of CALIOP aerosol extinction retrievals, we first derive lidar ratios of aerosol layers globally, classified as dust by CALIOP, by using the backscatter measurements constrained by collocated MODIS aerosol optical depths. Since this dust will frequently be a mixture of dust from various source regions, we obtain the corresponding regional contributions from the dust simulations using the GEOS model tagged to nine source regions. Finally, combining information from both measurements and models, we obtain the lidar ratios of dust corresponding to the nine source regions. The lidar ratio of the mixed dust anywhere can now be obtained from the tagged model fractions and the corresponding regional lidar ratio. We shall show the application of this scheme for various regions and seasons.

CALIPSO↗

NASA NTRS1975

Evidence for meteoritic sulfur in the lunar regolith

Abundance and isotopic data are presented for a suite of 19 soils from the Apollo 16 mission. It is concluded that the S systematics of the lunar regolith support a model in which S from a number of indigenous rock types is mixed with a time-dependent, extralunar component of probably meteoritic origin. This mixture is processed on the lunar surface by a time-dependent mechanism which removes S-32 enriched S from the moon, enriching the residue in S-34. The loss of S by this mechanism approximately equals the meteoritic input, producing the apparent coincidence between rock abundances and abundances in soils which has led to adoption of simplified mixing models. In reality, at least for the Apollo 16 site, addition of meteoritic S slightly outweighs loss from the moon, resulting in a positive trend.

Kerridge, J. F.↗

DOE OSTI2023

A High-Pressure Shock Tube Study of Hydrogen and Ammonia Addition to Natural Gas for Reduced Carbon Emissions in Power Generation Gas Turbines

Ignition delay times from undiluted mixtures of natural gas (NG)/H2/Air and NG/NH3/Air were measured using a high-pressure shock tube at the University of Central Florida. The combustion temperatures were experimentally tested between 1000-1500 K near a constant pressure of 25 bar. 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 is 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 H2 addition and all temperature ranges of the baseline NG mixture. Hydrogen addition also decreased ignition delay times by nearly 90%, while NH3 fuel addition made no noticeable difference in ignition delay time. NG/NH3 exhibited similar chemistry to pure NG under the same conditions, which is shown in a sensitivity analysis, demonstrating hydrogen chemistry to be dominant in NG/H2 mixtures and hydrocarbon chemistry to be dominant in NG/NH3 mixtures. The reaction CH3 + O2 = CH3O + 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.

Pierro, Michael↗

DOE OSTI2021

The Influence of the Instantaneous Collapse of Tailings Pond on Downstream Facilities

To evaluate the evolutionary processes guiding the formation of the tailings‐water mixtures produced by the instantaneous collapse of tailings ponds and the influence of these on downstream facilities, a 2D simulation model with reasonable boundary and working conditions derived from actual engineering practice was built in this study, and the relationship between dam‐break elevation and impact on downstream facilities was also analyzed to determine the relevant mechanism of influence. Computational results indicated that lowering the dam‐break elevation caused the maximum velocity and flooding depth, along with the flooded area at monitoring points, to gradually increase. The occurrence times of maximum velocity and flooding depth were also gradually moved forward as the breaking elevation was reduced; this effect is directly related to the increase in the total potential energy at the lower break elevations. Further simulations of sand‐prevent dams with different heights located downstream from a tailings pond were carried out to identify methods for mitigating the impact of dam failure. The results revealed that increasing the height of the sand‐prevent dam reduced the production of tailings mixtures. Based on the results, the construction of a sand‐prevent dam with a crest elevation equal to that of the starter dam was recommended.

Wang, Sha (ORCID:0000000324875445)↗

NASA NTRSDate not supplied

Community Radiative Transfer Model: Implementing A New Cloud Scattering Database and Developing the Forward Radar Module

The Mie theory is used by many fast RT models to estimate the optical properties of single particles. The Mie theory assumes spherical shapes for ice or snow particles with mixture of air and ice. However, hydrometeors scattering radiation at microwave frequencies have different shapes, sizes, and orientations. Therefore, using Mie theory to determine their optical properties leads to large uncertainties in all-sky radiative transfer calculations. The discrete dipole approximation (DDA) which approximates the optical properties of large objects in terms of discrete dipoles has shown promises in calculating the scattering properties of particles with different shapes in the microwave frequencies. The goal of the research was to enhance the CRTM scattering calculations for frozen hydrometeors in the microwave frequencies using the DDA technique. Given that such optical properties cannot be practically calculated on the fly, pre-computed look-up tables need to be implemented into the fast RT models to calculate the scattering properties of these particles using the DDA technique and the inputs provided by the users. Therefore, we implemented such pre-computed DDA databases into CRTM. In addition to using stand-alone CRTM calculations using collocated ATMS and reanalysis profiles, the data assimilation experiments conducted using the NOAA FV3GFS forecast system will be used to evaluate the scattering improvements. Additionally, we used the backscattering information from the DDA database to implement a radar simulator into CRTM. The radar operators takes advantage of CRTM different modules to calculate clouds absorption and scattering properties. In addition to the forward model both adjoint and tangent linear of the radar simulator are implemented and evaluated as well. The radar simulator is currently being tested within the JEDI/GEOS data assimilation framework to facilitate the assimilation of radar measurements such as CloudSat CPR and GPM DPR into the NASA GEOS model.

Isaac Moradi↗

DOE OSTI2025

CFD modeling of non-catalytic, partial-oxidation engine reformer for flare mitigation

Flaring associated natural gas is commonly employed in the oil and gas industry to reduce methane (CH 4 ) emissions but generates carbon dioxide (CO 2 ) and harmful pollutants, significantly contributing to air pollution and posing risks to public health. To mitigate this impact, M2X Energy Inc. has developed a small-scale, modular gas-to-methanol system. This system features an engine reformer that performs fuel-rich partial oxidation of wellhead gas to produce syngas—a mixture of carbon monoxide (CO) and hydrogen (H 2 )—followed by a downstream reactor for methanol synthesis. This study focused on computational fluid dynamics (CFD) modeling of the engine reformer to simulate partial oxidation chemistry, predict the rich-burn operating limit, and assess syngas quality, ultimately aiding in design and operational optimization. The CFD model, developed within a Reynolds-Averaged Navier-Stokes (RANS) turbulence framework, incorporated sub-models for turbulent combustion, a chemical mechanism with polycyclic aromatic hydrocarbon (PAH) pathways, and soot emissions to accurately capture the fuel-rich, turbulent jet ignition and combustion processes. Model validation against experimental data showed good agreement across pre- and main-chamber pressures, apparent heat release rates, and exhaust gas concentrations of key species (H 2 , CO, CO 2 , CH 4 ) for varying intake equivalence ratios. Here, the model identified a rich-burn operating limit near a fuel-air equivalence ratio of 2.35, consistent with experimental observations. Furthermore, syngas quality analysis revealed that extending the rich-burn limit through engine reformer optimization could enhance syngas production, contributing to higher methanol synthesis efficiency.

Computational Fluid Dynamics↗