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 559 records · Page 31

Composition, Structure and Formation of Low-Density Planets Within 0.5 AU of Their Star

Kepler has discovered that tens of percent of sunlike stars possess one or more planets larger than Earth within 0.5 AU. Masses of dozens of small and mid-sized exoplanets have been measured using radial velocity or transit timing variations. Few transiting planets larger than 1.6 Earth radii are rocky. Rather, the planetary mass-radius diagram for planetary radius between 1.6 and 8 times that of Earth shows a tendency for mass to increase slowly with size. The scatter is larger than can be explained by observational errors and correlations with temperature, implying that the population is heterogeneous. Many Neptune-size planets have very low densities and must have a substantial fraction of their volumes occupied by hydrogen and/or helium. Models of the structure, growth and evolution of planets made of mixtures of rock and light gases will be presented and compared with observations of both individual bodies and of the characteristics of the population as a whole.

radial velocity↗

Integrated High Payoff Rocket Propulsion Technology (IHPRPT) SiC Recession Model

SiC stability and recession rates were modeled in hydrogen/oxygen combustion environments for the Integrated High Payoff Rocket Propulsion Technology (IHPRPT) program. The IHPRPT program is a government and industry program to improve U.S. rocket propulsion systems. Within this program SiC-based ceramic matrix composites are being considered for transpiration cooled injector faceplates or rocket engine thrust chamber liners. Material testing under conditions representative of these environments was conducted at the NASA Glenn Research Center, Cell 22. For the study described herein, SiC degradation was modeled under these Cell 22 test conditions for comparison to actual test results: molar mixture ratio, MR (O2:H2) = 6, material temperatures to 1700 C, combustion gas pressures between 0.34 and 2.10 atm, and gas velocities between 8,000 and 12,000 fps. Recession was calculated assuming rates were controlled by volatility of thermally grown silica limited by gas boundary layer transport. Assumptions for use of this model were explored, including the presence of silica on the SiC surface, laminar gas boundary layer limited volatility, and accuracy of thermochemical data for volatile Si-O-H species. Recession rates were calculated as a function of temperature. It was found that at 1700 C, the highest temperature considered, the calculated recession rates were negligible, about 200 m/h, relative to the expected lifetime of the material. Results compared favorably to testing observations. Other mechanisms contributing to SiC recession are briefly described including consumption of underlying carbon and pitting. A simple expression for liquid flow on the material surface was developed from a one-dimensional treatment of the Navier-Stokes Equation. This relationship is useful to determine under which conditions glassy coatings or thermally grown silica would flow on the material surface, removing protective layers by shear forces. The velocity of liquid flow was found to depend on the gas velocity, the viscosity of gas and liquid, as well as the thickness of the gas boundary layer and the liquid layer. Calculated flow rates of a borosilicate glass coating compared well to flow rates observed for this coating tested on a SiC panel in Cell 22.

Opila, E. J.↗

Numerical simulation of comet nuclei. I - Water-ice comets

A one-dimensional numerical model of pure water-ice cometary nuclei is presented, and the influence of the nuclear interior as a heat reservoir on the behavior of the nuclear surface is examined. It is shown that a number of effects, including the thermal inertia due to heat stored in the core and the release of latent heat, which goes entirely into heating the adjacent layers or into sublimation on passing through a phase transition from amorphous to crystalline ice, can help to explain such characteristics as the asymmetrical lightcurve of Comet Halley. Results are given for the cases of Comet Schwassmann-Wachmann 1 and Comet Encke. Consideration is also given to the insulating effect of an evolving dust mantle. The role of this mantle in determining the surface temperature of the ice core is studied as a function of the mass fraction of the dust in the ice-dust mixture and the thermal conductivity of the nucleus. The loose-lattice model of Mendis and Brin (1977) indicates that both high dust to ice ratios and high-core conductivities inhibit mantle blowoff.

Herman, G.↗

A subgrid mixing model for LES of nonpremixed turbulent combustion

A new sub-grid mixing model for use in large eddy simulations of turbulent combustion is presented and applied to a hydrogen-air diffusion flame. The sub-grid model is based on Kerstein's Linear Eddy Model (1988). The model is first used to predict the mixing of a conserved scalar in a turbulent shear flow. The model correctly predicts the behavior of the pdf of the scalar field and displays a nonmarching peak at the preferred mixture fraction as the shear layer is traversed. It is then illustrated how a reduced chemical mechanism can be implemented within the linear eddy subgrid model formulation. The model is used to predict NO formation in hydrogen-air diffusion flame using a reduced chemical mechanism involving nine reactive scalars.

Mcmurtry, P. A.↗

Effects of Radiative Emission and Absorption on the Propagation and Extinction of Premixed Gas Flames

Premixed gas flames in mixtures of CH4, O2, N2, and CO2 were studied numerically using detailed chemical and radiative emission-absorption models to establish the conditions for which radiatively induced extinction limits may exist independent of the system dimensions. It was found that reabsorption of emitted radiation led to substantially higher burning velocities and wider extinction limits than calculations using optically thin radiation models, particularly when CO2, a strong absorber, is present in the unburned gas, Two heat loss mechanisms that lead to flammability limits even with reabsorption were identified. One is that for dry hydrocarbon-air mixtures, because of the differences in the absorption spectra of H2O and CO2, most of the radiation from product H2O that is emitted in the upstream direction cannot be absorbed by the reactants. The second is that the emission spectrum Of CO2 is broader at flame temperatures than ambient temperature: thus, some radiation emitted near the flame front cannot be absorbed by the reactants even when they are seeded with CO2 Via both mechanisms, some net upstream heat loss due to radiation will always occur, leading to extinction of sufficiently weak mixtures. Downstream loss has practically no influence. Comparison with experiment demonstrates the importance of reabsorption in CO2 diluted mixtures. It is concluded that fundamental flammability limits can exist due to radiative heat loss, but these limits are strongly dependent on the emission-absorption spectra of the reactant and product -gases and their temperature dependence and cannot be predicted using gray-gas or optically thin model parameters. Applications to practical flames at high pressure, in large combustion chambers, and with exhaust-gas or flue-gas recirculation are discussed.

Ju, Yiguang↗

Modeling the reflectance spectrum of Callisto 0.25 to 4.1 microns

Models employing Hapke's (1981) radiative transfer theory are presented for the 0.2-4.1 micron reflectance spectrum. A simultaneous intimate, as well as aereal, mixture solution of ice and dark material is envisioned, in order to satisfy both absorption band depths and reflectance levels. The models indicate that the surface's ice component is rather large-grained, and that the major spectral features above about 2.5 microns cannot be accounted for by the ice. Spectra obtained for the nonice material were similar to each other; their absorption features resemble hydrated silicates bearing both oxidation states of iron.

Calvin, Wendy M.↗

Energy transfer by radiation in non-grey atomic gases in isothermal and non-isothermal slabs

A multiband model for the absorption coefficient of atomic hydrogen-helium plasmas is constructed which includes continuum and line contributions. Emission from 28 stronger lines of 106 that have been screened is considered, of which 21 are from hydrogen and 7 belong to helium, with reabsorption due to line-line, line-continuum overlap accurately accounted for. The model is utilized in the computation of intensities and fluxes from shock-heated slabs of 85% H2-15% He mixtures for slab thicknesses from 1 to 30 cm, temperature from 10,000 to 20,000 K, and for different densities. In conjunction with the multiband model, simple numerical schemes have been devised which provide a quick and comprehensive way of computing radiative energy transfer in nonisothermal and nongrey gases.

Poon, P. T. Y.↗

The Dripping Handrail Model: Transient Chaos in Accretion Systems

We define and study a simple dynamical model for accretion systems, the "dripping handrail" (DHR). The time evolution of this spatially extended system is a mixture of periodic and apparently random (but actually deterministic) behavior. The nature of this mixture depends on the values of its physical parameters - the accretion rate, diffusion coefficient, and density threshold. The aperiodic component is a special kind of deterministic chaos called transient chaos. The model can simultaneously exhibit both the quasiperiodic oscillations (QPO) and very low frequency noise (VLFN) that characterize the power spectra of fluctuations of several classes of accretion systems in astronomy. For this reason, our model may be relevant to many such astrophysical systems, including binary stars with accretion onto a compact object - white dwarf, neutron star, or black hole - as well as active galactic nuclei. We describe the systematics of the DHR's temporal behavior, by exploring its physical parameter space using several diagnostics: power spectra, wavelet "scalegrams," and Lyapunov exponents. In addition, we note that for large accretion rates the DHR has periodic modes; the effective pulse shapes for these modes - evaluated by folding the time series at the known period - bear a resemblance to the similarly- determined shapes for some x-ray pulsars. The pulsing observed in some of these systems may be such periodic-mode accretion, and not due to pure rotation as in the standard pulsar model.

Young, Karl↗

NASA Numerical and Experimental Evaluation of UTRC Low Emissions Injector

Computational and experimental analyses of a PICS-Pilot-In-Can-Swirler technology injector, developed by United Technologies Research Center (UTRC) are presented. NASA has defined technology targets for near term (called "N+1", circa 2015), midterm ("N+2", circa 2020) and far term ("N+3", circa 2030) that specify realistic emissions and fuel efficiency goals for commercial aircraft. This injector has potential for application in an engine to meet the Pratt & Whitney N+3 supersonic cycle goals, or the subsonic N+2 engine cycle goals. Experimental methods were employed to investigate supersonic cruise points as well as select points of the subsonic cycle engine; cruise, approach, and idle with a slightly elevated inlet pressure. Experiments at NASA employed gas analysis and a suite of laser-based measurement techniques to characterize the combustor flow downstream from the PICS dump plane. Optical diagnostics employed for this work included Planar Laser-Induced Fluorescence of fuel for injector spray pattern and Spontaneous Raman Spectroscopy for relative species concentration of fuel and CO2. The work reported here used unheated (liquid) Jet-A fuel for all fuel circuits and cycle conditions. The initial tests performed by UTRC used vaporized Jet-A to simulate the expected supersonic cruise condition, which anticipated using fuel as a heat sink. Using the National Combustion Code a PICS-based combustor was modeled with liquid fuel at the supersonic cruise condition. All CFD models used a cubic non-linear k-epsilon turbulence wall functions model, and a semi-detailed Jet-A kinetic mechanism based on a surrogate fuel mixture. Two initial spray droplet size distribution and spray cone conditions were used: (1) an initial condition (Lefebvre) with an assumed Rosin-Rammler distribution, and 7 degree Solid Spray Cone; and (2) the Boundary Layer Stripping (BLS) primary atomization model giving the spray size distribution and directional properties. Contour and line plots are shown in comparison with experimental data (where this data is available) for flow velocities, fuel, and temperature distribution. The CFD results are consistent with experimental observations for fuel distribution and vaporization. Analysis of gas sample results, using a previously-developed NASA NOx correlation, indicates that for sea-level takeoff, the PICS configuration is predicted to deliver an EINOx value of about three for the targeted supersonic aircraft. Emissions results at supersonic cruise conditions show potential for meeting the NASA goals with liquid fuel.

Hicks, Yolanda R.↗

NASA Numerical and Experimental Evaluation of UTRC Low Emissions Injector

Computational and experimental analyses of a PICS-Pilot-In-Can-Swirler technology injector, developed by United Technologies Research Center (UTRC) are presented. NASA has defined technology targets for near term (called "N+1", circa 2015), midterm ("N+2", circa 2020) and far term ("N+3", circa 2030) that specify realistic emissions and fuel efficiency goals for commercial aircraft. This injector has potential for application in an engine to meet the Pratt & Whitney N+3 supersonic cycle goals, or the subsonic N+2 engine cycle goals. Experimental methods were employed to investigate supersonic cruise points as well as select points of the subsonic cycle engine; cruise, approach, and idle with a slightly elevated inlet pressure. Experiments at NASA employed gas analysis and a suite of laser-based measurement techniques to characterize the combustor flow downstream from the PICS dump plane. Optical diagnostics employed for this work included Planar Laser-Induced Fluorescence of fuel for injector spray pattern and Spontaneous Raman Spectroscopy for relative species concentration of fuel and CO2. The work reported here used unheated (liquid) Jet-A fuel for all fuel circuits and cycle conditions. The initial tests performed by UTRC used vaporized Jet-A to simulate the expected supersonic cruise condition, which anticipated using fuel as a heat sink. Using the National Combustion Code a PICS-based combustor was modeled with liquid fuel at the supersonic cruise condition. All CFD models used a cubic non-linear k-epsilon turbulence wall functions model, and a semi-detailed Jet-A kinetic mechanism based on a surrogate fuel mixture. Two initial spray droplet size distribution and spray cone conditions were used: 1) an initial condition (Lefebvre) with an assumed Rosin-Rammler distribution, and 7 degree Solid Spray Cone; and 2) the Boundary Layer Stripping (BLS) primary atomization model giving the spray size distribution and directional properties. Contour and line plots are shown in comparison with experimental data (where this data is available) for flow velocities, fuel, and temperature distribution. The CFD results are consistent with experimental observations for fuel distribution and vaporization. Analysis of gas sample results, using a previously-developed NASA NOx correlation, indicates that for sea-level takeoff, the PICS configuration is predicted to deliver an EINOx value of about 3 for the targeted supersonic aircraft. Emissions results at supersonic cruise conditions show potential for meeting the NASA goals with liquid fuel.

Hicks, Yolanda R.↗

A Small Aircraft Transportation System (SATS) Demand Model

The Small Aircraft Transportation System (SATS) demand modeling is a tool that will be useful for decision-makers to analyze SATS demands in both airport and airspace. We constructed a series of models following the general top-down, modular principles in systems engineering. There are three principal models, SATS Airport Demand Model (SATS-ADM), SATS Flight Demand Model (SATS-FDM), and LMINET-SATS. SATS-ADM models SATS operations, by aircraft type, from the forecasts in fleet, configuration and performance, utilization, and traffic mixture. Given the SATS airport operations such as the ones generated by SATS-ADM, SATS-FDM constructs the SATS origin and destination (O&D) traffic flow based on the solution of the gravity model, from which it then generates SATS flights using the Monte Carlo simulation based on the departure time-of-day profile. LMINET-SATS, an extension of LMINET, models SATS demands at airspace and airport by all aircraft operations in US The models use parameters to provide the user with flexibility and ease of use to generate SATS demand for different scenarios. Several case studies are included to illustrate the use of the models, which are useful to identify the need for a new air traffic management system to cope with SATS.

Long, Dou↗

Cloud-Based Solutions for Monitoring Coastal Ecosystems and Prioritization of Restoration Efforts Across Belize

In recent years the availability of automated change detection algorithms in Google Earth Engine has permitted cloud-based processing of large time series of satellite imagery. Models such as the Continuous Change Detection and Classification (CCDC), CCDC-Spectral Mixture Analysis (CCDC-SMA), and Landsat-based Detection of Trends in Disturbance and Recovery (LandTrendr) allow users to exploit decades of Earth Observations (EO) , leveraging the Landsat archive and data from other sensors to detect disturbance in forest ecosystems. Despite the wide adoption of these methods, robust documentation and growing community of users, little research has explored their use in mangrove environments. Mangroves are dynamic environments subject to changes due to not only the natural migration of mudflats but also coastal erosion, urban expansion, aquaculture practices, etc. This work aims to identify best practices for the application of these models to identify and monitor changes in the Belizean mangroves, which experienced an estimated 5.4% decrease in extent between 1980 and 2017 (Cherrington et al. 2020). Partnering directly with the Belizean Forest Department, our team will develop a replicable, efficient methodology to annually update the country’s mangrove extent employing EO-based change detection. This collaboration supports Belize’s marine conservation and climate resilience efforts, such as its Blue Bond agreement with the Nature Conservancy, building on national processes to monitor these critical ecosystems.

Mangroves↗

Lagrange thermodynamic potential and intrinsic variables for He-3 He-4 dilute solutions

For a two-fluid model of dilute solutions of He-3 in liquid He-4, a thermodynamic potential is constructed that provides a Lagrangian for deriving equations of motion by a variational procedure. This Lagrangian is defined for uniform velocity fields as a (negative) Legendre transform of total internal energy, and its primary independent variables, together with their thermodynamic conjugates, are identified. Here, similarities between relations in classical physics and quantum statistical mechanics serve as a guide for developing an alternate expression for this function that reveals its character as the difference between apparent kinetic energy and intrinsic internal energy. When the He-3 concentration in the mixtures tends to zero, this expression reduces to Zilsel's formula for the Lagrangian for pure liquid He-4. An investigation of properties of the intrinsic internal energy leads to the introduction of intrinsic chemical potentials along with other intrinsic variables for the mixtures. Explicit formulas for these variables are derived for a noninteracting elementary excitation model of the fluid. Using these formulas and others also derived from quantum statistical mechanics, another equivalent expression for the Lagrangian is generated.

Jackson, H. W.↗

Profile Evolution Simulation in Etching Systems Using Level Set Methods

Semiconductor device profiles are determined by the characteristics of both etching and deposition processes. In particular, a highly anisotropic etch is required to achieve vertical sidewalls. However, etching is comprised of both anisotropic and isotropic components, due to ion and neutral fluxes, respectively. In Ar/Cl2 plasmas, for example, neutral chlorine reacts with the Si surfaces to form silicon chlorides. These compounds are then removed by the impinging ion fluxes. Hence the directionality of the ions (and thus the ion angular distribution function, or IAD), as well as the relative fluxes of neutrals and ions determines the amount of undercutting. One method of modeling device profile evolution is to simulate the moving solid-gas interface between the semiconductor and the plasma as a string of nodes. The velocity of each node is calculated and then the nodes are advanced accordingly. Although this technique appears to be relatively straightforward, extensive looping schemes are required at the profile corners. An alternate method is to use level set theory, which involves embedding the location of the interface in a field variable. The normal speed is calculated at each mesh point, and the field variable is updated. The profile comers are more accurately modeled as the need for looping algorithms is eliminated. The model we have developed is a 2-D Level Set Profile Evolution Simulation (LSPES). The LSPES calculates etch rates of a substrate in low pressure plasmas due to the incident ion and neutral fluxes. For a Si substrate in an Ar/C12 gas mixture, for example, the predictions of the LSPES are identical to those from a string evolution model for high neutral fluxes and two different ion angular distributions.(2) In the figure shown, the relative neutral to ion flux in the bulk plasma is 100 to 1. For a moderately isotropic ion angular distribution function as shown in the cases in the left hand column, both the LSPES (top row) and rude's string method (bottom row) predict tapered profiles. The LSPES uses an AD with a FWHM = 13.5 degrees, and rude's model uses a ratio of sheath voltage to ion temperature of 50. The more anisotropic IADs produce profiles with more vertical sidewalls and a wider bottom surface, as shown in the right hand column. Here, the LSPES has an AD with a FWHM = 2 degrees, and rude's model uses a ratio of 500 of sheath voltage to ion temperature. The agreement between the LSPES and rude's model is excellent in both cases. We will present etching profiles generated by the LSPES, including calculations of the re-emitted fluxes of both neutrals and ions off of the profile walls. In addition, we will show the effect of geometric structures (overhangs, etc.) on the etching profiles. Other physical aspects, such as surface diffusion, will also be included in the model.

Hwang, Helen H.↗

Turbulent flame propagation in partially premixed flames

Turbulent premixed flame propagation is essential in many practical devices. In the past, fundamental and modeling studies of propagating flames have generally focused on turbulent flame propagation in mixtures of homogeneous composition, i.e. a mixture where the fuel-oxidizer mass ratio, or equivalence ratio, is uniform. This situation corresponds to the ideal case of perfect premixing between fuel and oxidizer. In practical situations, however, deviations from this ideal case occur frequently. In stratified reciprocating engines, fuel injection and large-scale flow motions are fine-tuned to create a mean gradient of equivalence ratio in the combustion chamber which provides additional control on combustion performance. In aircraft engines, combustion occurs with fuel and secondary air injected at various locations resulting in a nonuniform equivalence ratio. In both examples, mean values of the equivalence ratio can exhibit strong spatial and temporal variations. These variations in mixture composition are particularly significant in engines that use direct fuel injection into the combustion chamber. In this case, the liquid fuel does not always completely vaporize and mix before combustion occurs, resulting in persistent rich and lean pockets into which the turbulent flame propagates. From a practical point of view, there are several basic and important issues regarding partially premixed combustion that need to be resolved. Two such issues are how reactant composition inhomogeneities affect the laminar and turbulent flame speeds, and how the burnt gas temperature varies as a function of these inhomogeneities. Knowledge of the flame speed is critical in optimizing combustion performance, and the minimization of pollutant emissions relies heavily on the temperature in the burnt gases. Another application of partially premixed combustion is found in the field of active control of turbulent combustion. One possible technique of active control consists of pulsating the fuel flow rate and thereby modulating the equivalence ratio (Bloxsidge et al. 1987). Models of partially premixed combustion would be extremely useful in addressing all these questions related to practical systems. Unfortunately, the lack of a fundamental understanding regarding partially premixed combustion has resulted in an absence of models which accurately capture the complex nature of these flames. Previous work on partially premixed combustion has focused primarily on laminar triple flames. Triple flames correspond to an extreme case where fuel and oxidizer are initially totally separated (Veynante et al. 1994 and Ruetsch et al. 1995). These flames have a nontrivial propagation speed and are believed to be a key element in the stabilization process of jet diffusion flames. Different theories have also been proposed in the literature to describe a turbulent flame propagating in a mixture with variable equivalence ratio (Muller et al. 1994), but few validations are available. The objective of the present study is to provide basic information on the effects of partial premixing in turbulent combustion. In the following, we use direct numerical simulations to study laminar and turbulent flame propagation with variable equivalence ratio.

Poinsot, T.↗

Two-dimensional model of a Space Station Freedom thermal energy storage canister

The Solar Dynamic Power Module being developed for Space Station Freedom uses a eutectic mixture of LiF-CaF2 phase change salt contained in toroidal canisters for thermal energy storage. Results are presented from heat transfer analyses of the phase change salt containment canister. A 2-D, axisymmetric finite difference computer program which models the canister walls, salt, void, and heat engine working fluid coolant was developed. Analyses included effects of conduction in canister walls and solid salt, conduction and free convection in liquid salt, conduction and radiation across salt vapor filled void regions and forced convection in the heat engine working fluid. Void shape, location, growth or shrinkage (due to density difference between the solid and liquid salt phases) were prescribed based on engineering judgement. The salt phase change process was modeled using the enthalpy method. Discussion of results focuses on the role of free-convection in the liquid salt on canister heat transfer performance. This role is shown to be important for interpreting the relationship between ground based canister performance (in l-g) and expected on-orbit performance (in micro-g). Attention is also focused on the influence of void heat transfer on canister wall temperature distributions. The large thermal resistance of void regions is shown to accentuate canister hot spots and temperature gradients.

Kerslake, Thomas W.↗

A Study of Soil and Duricrust Models for Mars

An analysis of soil and duricrust formation mechanisms on Mars is presented. Soil analog mixtures have been prepared, characterized and tested through wet/dry cycling experiments; results are compared with Mars Pathfinder soil data (spectral, chemical and magnetic). Additional information is contained in the original extended abstract.

Bishop, J. L.↗

Two-dimensional model of a Space Station Freedom thermal energy storage canister

The Solar Dynamic Power Module being developed for Space Station Freedom uses a eutectic mixture of LiF-CaF2 phase change salt contained in toroidal canisters for thermal energy storage. Results are presented from heat transfer analyses of the phase-change salt containment canister. A 2-D, axisymmetric finite-difference computer program which models the canister walls, salt, void, and heat engine working fluid coolant was developed. Analyses included effects of conduction in canister walls and solid salt, conduction and free convection in liquid salt, conduction and radiation across salt vapor filled void regions, and forced convection in the heat engine working fluid. Void shape, location, and growth or shrinkage (due to density difference between the solid and liquid salt phases) were prescribed based on engineering judgement. The salt phase change process was modeled using the enthalpy method. Discussion of results focuses on the role of free-convection in the liquid salt on canister heat transfer performance. This role is shown to be important for interpreting the relationship between groundbased canister performance (in 1-g) and expected on-orbit performance (in micro-g). Attention is also focused on the influence of void heat transfer on canister wall temperature distributions. The large thermal resistance of void regions is shown to accentuate canister hot spots and temperature gradients.

Kerslake, Thomas W.↗