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 649 records · Page 36

Combustion of petroleum-based transportation fuels and their blends with biofuels: a new approach for developing surrogates and understanding the effects of blending

Accurate and versatile reaction mechanisms are necessary to simulate combustion of transportation fuels blended with biofuels for improving fuel efficiency, reducing fuel consumption and mitigating the formation and emissions of combustion products harmful to the atmosphere. This project involved a combined effort of simulation and experimentation of liquid fuel burning using a combustion configuration amenable to numerical modeling with a level of detail not previously achieved for blends of biofuels with gasoline certification and surrogate fuels. The configuration of an isolated droplet burning under conditions where gas transport arises solely from fuel evaporation was selected as the platform for experiments and numerical modeling. The spherical symmetry and one-dimensional gas transport that results enabled simulating the droplet burning process with an precedented level of detail. Processes associated with unsteady gas and liquid transport, formation of particulate and gaseous products, radiation, multicomponent phase equilibrium at the droplet surface and moving boundary effects from droplet evaporation, were incorporated in a single numerical framework. The numerical model was based on the open source code OPENSmoke++ (OS) adapted to incorporate the effects noted above. In addition to simulation, an experimental design was developed for the isolated droplet configuration to acquire the data needed for validating the numerical model. Several broad accomplishments of the project were the following: demonstrating generally excellent agreement between measured and simulated combustion parameters for the fuel blends examined that included mixtures of heptane and isobutanol as a model system, and surrogates comprised of up to seven miscible components mixed with ethanol or isobutanol as representative biofuel additives; using a new approach to validate reaction mechanisms of biofuel blends which incorporated fuel evaporation into the process along with developing an experimental design for acquiring data to compare with simulations; and showing that a wealth of information could be obtained on the combustion physics of biofuels from experiments requiring volumes on the order of only nanoliters at a time thus opening the way to evaluating biofuels synthesized by new processes early in development.

02 PETROLEUM↗

Cryogenic Extension of NASA Species Polynomials Using Hydrogen and Oxygen at Stoichiometry

The Helmholtz free energy for a mixture has been used in many chemistry and combustion studies for mixtures. An approach is described where polynomial fits intended for a thermally perfect ideal gases can be modified and implemented into a Helmholtz free energy model, such as those used by NIST in their REFPROP program. However, the limitation of the polynomial fits must be respected and should only used in a valid temperature interval where the flow field can be treated as a thermally perfect ideal gas when these species are present. This approach can allow a CFD simulation to have cryogenic regions of a flow field where high fidelity cryogenic species models are needed and allows for the use of species in other areas at higher temperatures where combustion species do not have such a model available or is not needed. Presented in the paper are CFD 1D detonation wave and species mass fraction simulation results which have cryogenic upstream conditions of Hydrogen and Oxygen at stoichiometry.

Combustion↗

Development of simplified model for injection rate prediction of diesel injectors during transient and steady operation

Determining the injection rate of the diesel injector during the transient short injection operations is a critical matter to enhance the combustion process by improving air/fuel mixture formation. In this regard, numerous injection models have been suggested in the literature to predict the injection rate accurately. However, such algorithms have not precisely predicted the rate of injection (ROI) in transient short injection events. In addition, most of those models require many input parameters to characterize hydraulic and mechanical subsystems in the injector, which can make them inconvenient to tune for different injector specifications, in other words not suitable for universal usage. Therefore, our study aims to develop a simple injection model that can predict injection rate precisely in the short and long injection durations considering relatively fewer input parameters. The model traces the pressure variation inside the sac based on the bulk modulus theory to predict nozzle exit velocity, thus the injection rate of the injector. In the algorithm, the bulk modulus of fuel (E) and discharge coefficients of the nozzle hole and needle seat (C d,H and C d,N ) have been considered to vary rather than set constant during the injection process since such flow parameters would have a critical impact on the flow development at the transient stage, thus affecting the injection rate of diesel injectors. Here with this modelling strategy, the accuracy of the new model on ROI prediction was confirmed with acceptable model error for both short and long injection durations by validating against experimental results. To understand the effect of the new model scheme on ROI prediction, ROI results with the new model have been discussed by comparing it with the ROI results predicted with constant flow parameters.

42 ENGINEERING↗

Low-Cost Sensor Performance Intercomparison, Correction Factor Development, and 2+ Years of Ambient PM2.5 Monitoring in Accra, Ghana

Particulate matter air pollution is a leading cause of global mortality, particularly in Asia and Africa. Addressing the high and wide-ranging air pollution levels requires ambient monitoring, but many low- and middle-income countries (LMICs) remain scarcely monitored. To address these data gaps, recent studies have utilized low-cost sensors. These sensors have varied performance, and little literature exists about sensor intercomparison in Africa. By colocating 2 QuantAQ Modulair-PM, 2 PurpleAir PA-II SD, and 16 Clarity Node-S Generation II monitors with a reference-grade Teledyne monitor in Accra, Ghana, we present the first intercomparisons of different brands of low-cost sensors in Africa, demonstrating that each type of low-cost sensor PM2.5 is strongly correlated with reference PM2.5, but biased high for ambient mixture of sources found in Accra. When compared to a reference monitor, the QuantAQ Modulair-PM has the lowest mean absolute error at 3.04 μg/m3, followed by PurpleAir PA-II (4.54 μg/m3) and Clarity Node-S (13.68 μg/m3). We also compare the usage of 4 statistical or machine learning models (Multiple Linear Regression, Random Forest, Gaussian Mixture Regression, and XGBoost) to correct low-cost sensors data, and find that XGBoost performs the best in testing (R2: 0.97, 0.94, 0.96; mean absolute error: 0.56, 0.80, and 0.68 μg/m3 for PurpleAir PA-II, Clarity Node-S, and Modulair-PM, respectively), but tree-based models do not perform well when correcting data outside the range of the colocation training. Therefore, we used Gaussian Mixture Regression to correct data from the network of 17 Clarity Node-S monitors deployed around Accra, Ghana, from 2018 to 2021. We find that the network daily average PM2.5 concentration in Accra is 23.4 μg/m3, which is 1.6 times the World Health Organization Daily PM2.5 guideline of 15 μg/m3. While this level is lower than those seen in some larger African cities (such as Kinshasa, Democratic Republic of the Congo), mitigation strategies should be developed soon to prevent further impairment to air quality as Accra, and Ghana as a whole, rapidly grow.

Humidity↗

Foundations for statistical-physical precipitation retrieval from passive microwave satellite measurements. I - Brightness-temperature properties of a time-dependent cloud-radiation model

The relationship between emerging microwave brightness temperatures (T(B)s) and vertically distributed mixtures of liquid and frozen hydrometeors was investigated, using a cloud-radiation model, in order to establish the framework for a hybrid statistical-physical rainfall retrieval algorithm. Although strong relationships were found between the T(B) values and various rain parameters, these correlations are misleading in that the T(B)s are largely controlled by fluctuations in the ice-particle mixing ratios, which in turn are highly correlated to fluctuations in liquid-particle mixing ratios. However, the empirically based T(B)-rain-rate (T(B)-RR) algorithms can still be used as tools for estimating precipitation if the hydrometeor profiles used for T(B)-RR algorithms are not specified in an ad hoc fashion.

Smith, Eric A.↗

High-precision monitoring of outgassing species in model extreme ultraviolet photoresists with a cavity ring-down spectrometer

Background Extreme ultraviolet (EUV) photoresists play a pivotal role in advancing nanopatterning technologies by balancing image quality and sensitivity. The outgassing behavior of photoresist thin films under EUV and deep ultraviolet (DUV) exposure reveals chemical details relevant to their performance. Aim Here, we focus on utilizing an analytical technique not previously used in photolithography, the tabletop cavity ring-down spectrometer, to investigate outgassing dynamics in EUV photoresists, enabling precise chemical identification and deeper insights into resist processing. Approach The spectrometer’s enhanced laser path length (∼ 20 km) and broadband absorption capabilities in the C–H overtone region allow for sensitive and temporally resolved detection of mixtures of outgassed species. Using a model resist comprising a polymer matrix with a photoacid generator and quencher, we analyzed the influence of time delays between exposure and post-exposure bake (PEB) as well as storage under varying environmental conditions. Results Suppression of isobutylene outgassing and thickness loss was observed with extended delays between exposure and PEB, potentially linked to water absorption and acid deactivation. The technique proved highly effective in distinguishing subtle chemical differences between processing stages. Conclusions Delay times and their environmental conditions, particularly humidity, reduce outgassing and thickness loss of photoresists during PEB, suggesting decreased acid-driven deprotection. This can potentially impact sensitivity, defectivity, and roughness of resist patterns, necessitating precise monitoring and control.

Extreme ultraviolet photoresists↗

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗

Determination of the temperature of gas mixtures by using laser Raman scattering

A theoretical model for the Raman rotational spectrum of air was formulated, yielding the spectral temperature dependence; and an error analysis of temperature measurement accuracies was performed. The theoretical spectra and the predicted temperature measurement error of + or - 20 K (7 percent) compared well with experimental data. These experimental data consisted of Raman spectra of air obtained for temperatures ranging from 243 to 313 K at a pressure of 1 atmosphere. The use of selected lines resulted in an experimental accuracy of + or - 13 K (-5 percent).

Coney, T. A.↗

Relationships of soil, grass, and bedrock over the Kaweah serpentine melange through spectral mixture analysis of AVIRIS data

A linear mixing model is used to model the spectral variability of an AVIRIS scene from the western foothills of the Sierra Nevada and calibrate these radiance data to reflectance. Five spectral endmembers from the AVIRIS data, plus an ideal 'shade' endmember were required to model the continuum reflectance of each pixel in the image. Three of the endmembers were interpreted to model the surface constituents green vegetation, dry grass, and illumination. These are the main transient surface constituents that are expected to change with shifts in land use or climatic influences and viewing conditions ('shade' only). The spectral distinction between the other three endmembers is very small, yet the spatial distributions are coherent and interpretable. These distributions cross anthropogenic and vegetation boundaries and are best interpreted as different soil types. Comparison of the fraction images to the bedrock geology maps indicates that substrate composition must be a factor contributing to the spectral properties of these endmembers. Detailed examination of the reflectance spectra of the three soil endmembers reveals that differences in the amount of ferric and ferrous iron and/or organic constituents in the soils is largely responsible for the differences in spectral properties of these endmembers.

Mustard, John F.↗

Relationships of soil, grass, and bedrock over the Kaweah Serpentinite Melange through spectral mixture analysis of AVIRIS data

A linear mixing model is used to model the spectral variability of an AVIRIS scene from the western foothills of the Sierra Nevada and calibrate these radiance data to reflectance. Five spectral endmembers from the AVIRIS data, plus an ideal 'shade' endmember were required to model the continuum reflectance of each pixel in the image. Three of the endmembers were interpreted to model the surface constituents green vegetation, dry grass, and illumination. Comparison of the fraction images to the bedrock geology maps indicates that substrate composition must be a factor contributing to the spectral properties of these endmembers. Detailed examination of the reflectance spectra of the three soil endmembers reveals that differences in the amount of ferric and ferrous iron and/or organic constituents in the soils is largely responsible for the differences in spectral properties of these endmembers.

Mustard, John F.↗

Effects of Chlorine Capture and a Proposed Density Law on the Reactivity of Plutonium Solution Systems

During fissionable material processing, all normal and credible abnormal conditions must remain safely subcritical. Nuclear Criticality Safety (NCS) uses a number of methods to determine subcriticality, one of which is the use of neutron transport codes such as MCNP6. In order to create models for use with MCNP6, both the geometry and materials in fissionable material processes must be known, or assumptions must be made and quantified for the impact to bias. One of the systems with a significant amount of bias due to material modeling assumptions is in the area of aqueous plutonium processing. These solutions are typically plutonium nitrate solutions or plutonium chloride solutions, which are modeled as fictitious plutonium metal-water mixtures because little is known about the actual density of the solution and there is no current predictive capability approved for use at Los Alamos National Laboratory (LANL) for modeling them. This research is currently underway to fill the gap and develop an algorithm for use with MCNP6 to model the density of plutonium chloride solutions. The method is to be validated with experimental data for density, and also validated with critical experiments using MCNP6. Note that the Chlorine Worth Study (CWS) was performed in December 2021 to help bridge the gap in chlorine data for critical experiments, and is currently awaiting International Criticality Safety Benchmark Evaluation Project (ICSBEP) review. This study was performed by LANL at the National Criticality Experiments Research Center (NCERC) at the Nevada National Security Site (NNSS). Additional information regarding this experiment may currently be found in LA-UR- 22-29180. Additionally, the Chemistry-Actinide Analytical Chemistry (C-AAC) at LANL has performed a number of solution density measurements for PuCl 3 -HC 1 -H 2 O, allowing for such data be used to create a semi-empirical density via the Pitzer method. The published dataset for the measurements is documented in LA-UR-22-25454. This method has already been tested successfully for aqueous plutonium nitrate solutions in SCALE. Current solution density measurements exist of plutonium concentrations of 0-~142g/L, all at 2M HC1, at temperatures 20-40°C. Additional data was taken for HC1-corrected density values, which essentially mimics the data for a pure PuCl x -water solution. The calculations in this report aim to support the current research by demonstrating the difference in system reactivity for the current modeling method when compared to the new proposed modeling with a density law implementation, which is being written as a Python tool to be used with MCNP6.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗

Phase coexistence implications of violating Newton’s third law

Newton’s third law, action = reaction, is a foundational statement of classical mechanics. However, in natural and living systems, this law appears to be routinely violated for constituents interacting in a nonequilibrium environment. Here, we use computer simulations to explore the macroscopic phase behavior implications of breaking microscopic interaction reciprocity for a simple model system. We consider a binary mixture of attractive particles and introduce a parameter that is a continuous measure of the degree to which interaction reciprocity is broken. In the reciprocal limit, the species are indistinguishable, and the system phase separates into domains with distinct densities and identical compositions. Furthermore, increasing nonreciprocity is found to drive the system to explore a rich assortment of phases, including phases with strong composition asymmetries and three-phase coexistence. Many of the states induced by these forces, including traveling crystals and liquids, have no equilibrium analogs. By mapping the complete phase diagram for this model system and characterizing these unique phases, our findings offer a concrete path forward toward understanding how nonreciprocity shapes the structures found in living systems and how this might be leveraged in the design of synthetic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical Modeling of Glaciers in Martian Paleoclimates

Numerous geologic features suggest the presence of ice flow on the surface of mars. These features include lobate debris aprons, concentric crater fill, and lineated valley fill. The lateral extent of these features can range from 100 meters to over 20 km. Previous work has demonstrated that these features could not have formed in current Martian conditions. It has long been speculated that changes in Mars orbital properties, namely its obliquity, eccentricity, and argument of perihelion, can result in dramatic changes to climate. Recent climate model studies have shown that at periods of increased obliquity north polar water ice is mobilized southward and deposited at low ad mid latitudes. Mid latitude accumulation of ice would provide the necessary conditions for rock glaciers to form. A time-marching, finite element glacier model is used to demonstrate the ability of ice and ice-rock mixtures to flow under Martian paleoclimate conditions. Input to this model is constrained by the NASA Ames Mars General Circulation Model (MGCM).

Colaprete, A.↗

Experimental and modeling study of water time histories during H 2 S-N 2 O combustion in a shock tube

Here, with the interest in directly burning sour gas in gas turbines, and the fact that even small amounts of H 2 S or its combustion products can alter combustion characteristics, many research studies have been performed to better understand the combustion chemistry of H 2 S. In the present study, the water formation was followed by laser absorption with N 2 O as an oxidant, instead of O 2 . Nitrous Oxide being essentially consumed via N 2 O (+M) ⇌ N 2 + O (+M), the water formation via the H 2 S + O route can then be probed to further validate the models. Three H 2 S/N 2 O mixtures diluted in 98% Ar were studied to cover the following range of equivalence ratios: 0.5, 1.0, and 2.0, over a wide range of temperatures (1580–1940 K) around atmospheric pressure. A chemical kinetic model was then developed first by validating the base N 2 O kinetics mechanism, then by investigating the experimental results presented herein. The N 2 O kinetics mechanism was updated with recent work on the low- and high-pressure limits for N 2 O decomposition (N 2 O (+M) ⇌ N 2 + O (+M)) as well as a review of rate coefficients for N 2 O + H ⇌ N 2 + OH from the literature. Good agreement is shown for H 2 /N 2 O mixtures. Updates were then made to the H 2 S kinetics mechanism, specifically, an update from the literature on SO 2 + H (+M) ⇌ SO + OH (+M) and an adjustment to SO + SH ⇌ S 2 + OH. Additionally, reactions between SH and N 2 O were determined using W1BD data and transition state theory which required the addition of an NNS sub-mechanism. With these updates, the mechanism provides good agreement with the H 2 S/N 2 O experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Physical Model to Determine Snowfall over Land by Microwave Radiometry

Because microwave brightness temperatures emitted by snow covered surfaces are highly variable, snowfall above such surfaces is difficult to observe using window channels that occur at low frequencies (v less than 100 GHz). Furthermore, at frequencies v less than or equal to 37 GHz, sensitivity to liquid hydrometeors is dominant. These problems are mitigated at high frequencies (v greater than 100 GHz) where water vapor screens the surface emission and sensitivity to frozen hydrometeors is significant. However the scattering effect of snowfall in the atmosphere at those higher frequencies is also impacted by water vapor in the upper atmosphere. This work describes the methodology and results of physically-based retrievals of snow falling over land surfaces. The theory of scattering by randomly oriented dry snow particles at high microwave frequencies appears to be better described by regarding snow as a concatenation of equivalent ice spheres rather than as a sphere with the effective dielectric constant of an air-ice mixture. An equivalent sphere snow scattering model was validated against high frequency attenuation measurements. Satellite-based high frequency observations from an Advanced Microwave Sounding Unit (AMSU-B) instrument during the March 5-6, 2001 New England blizzard were used to retrieve snowfall over land. Vertical distributions of snow, temperature and relative humidity profiles were derived from the Pennsylvania State University-National Center for Atmospheric Research (PSU-NCAR) fifth-generation Mesoscale Model (MM5). Those data were applied and modified in a radiative transfer model that derived brightness temperatures consistent with the AMSU-B observations. The retrieved snowfall distribution was validated with radar reflectivity measurements obtained from the National Oceanic and Atmospheric Administration (NOAA) National Weather Service (NWS) ground-based radar network.

Skofronick-Jackson, G.↗

Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels

Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.

10 SYNTHETIC FUELS↗

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields

The primary accomplishments of the project were as follows: (1) From an overall standpoint, the primary accomplishment of this research was the development of a complete gasdynamic-radiatively coupled nonequilibrium viscous shock layer solution method for axisymmetric blunt bodies. This method can be used for rapid engineering modeling of nonequilibrium re-entry flowfields over a wide range of conditions. (2) Another significant accomplishment was the development of an air radiation model that included local thermodynamic nonequilibrium (LTNE) phenomena. (3) As part of this research, three electron-electronic energy models were developed. The first was a quasi-equilibrium electron (QEE) model which determined an effective free electron temperature and assumed that the electronic states were in equilibrium with the free electrons. The second was a quasi-equilibrium electron-electronic (QEEE) model which computed an effective electron-electronic temperature. The third model was a full electron-electronic (FEE) differential equation model which included convective, collisional, viscous, conductive, vibrational coupling, and chemical effects on electron-electronic energy. (4) Since vibration-dissociation coupling phenomena as well as vibrational thermal nonequilibrium phenomena are important in the nonequilibrium zone behind a shock front, a vibrational energy and vibration-dissociation coupling model was developed and included in the flowfield model. This model was a modified coupled vibrational dissociation vibrational (MCVDV) model and also included electron-vibrational coupling. (5) Another accomplishment of the project was the usage of the developed models to investigate radiative heating. (6) A multi-component diffusion model which properly models the multi-component nature of diffusion in complex gas mixtures such as air, was developed and incorporated into the blunt body model. (7) A model was developed to predict the magnitude and characteristics of the shock wave precursor ahead of vehicles entering the Earth's atmosphere. (8) Since considerable data exists for radiating nonequilibrium flow behind normal shock waves, a normal shock wave version of the blunt body code was developed. (9) By comparing predictions from the models and codes with available normal shock data and the flight data of Fire II, it is believed that the developed flowfield and nonequilibrium radiation models have been essentially validated for engineering applications.

Carlson, Leland A.↗