Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conditional diffusion models”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Calculation of Propulsive Nozzle Flowfields in Multidiffusing Chemically Reacting Environments

An advanced engineering model has been developed to aid in the analysis and design of hydrogen/oxygen chemical rocket engines. The complete multispecies, chemically reacting and multidiffusing Navier-Stokes equations are modelled, including the Soret thermal diffusion and the Dufour energy transfer terms. In addition to the spectrum of multispecies aspects developed, the model developed in this study is also conservative in axisymmetric flow for both inviscid and viscous flow environments and the boundary conditions employ a viscous, chemically reacting, reference plane characteristics method. Demonstration cases are presented for a 1030:1 area ratio nozzle, a 25 lbf film cooled nozzle, and a transpiration cooled plug and spool rocket engine. The results indicate that the thrust coefficient predictions of the 1030:1 and the 25 lbf film cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements when all of the chemical reaction and diffusion terms are considered. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. The Soret thermal diffusion term is demonstrated to have a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle in both the laminar flow 1030:1 nozzle and the turbulent flow plug and spool nozzle analysis cases performed. Further, the Soret term was shown to represent an important fraction of the diffusion fluxes occurring in a transpiration cooled rocket engine.

Kacynski, Kenneth John↗
NASA NTRS1972

NASTRAN thermal analyzer: A general purpose finite element heat transfer computer program

The program not only can render temperature distributions in solids subjected to various thermal boundary conditions, including effects of diffuse-gray thermal radiation, but is fully compatible in capacity and in the finite-element model representation with that of its structural counterpart in the NASTRAN system. The development history of the finite-element approach for determining temperatures is summarized. The scope of analysis capability, program structure, features, and limitations are given with the objective of providing NASTRAN users with an overall veiw of the NASTRAN thermal analyzer.

Lee, H.↗
NASA NTRS1987

Spacecraft particulate contaminant redistribution

This paper describes analyses that were performed in support of the Hubble Space Telescope (HST)) particulate contamination control effort. The specific problems addressed include extension of available particle removal data to launch acoustic and random vibration conditions, development of an engineering model for transport of suspended particulates by airflow and in the presence of vehicle acceleration, turbulent diffusion, and migration of particulates over vibrating surfaces, and integration of the various models into a code that could be used to generate contamination level estimates for the HST primary mirror and other critical surfaces for the HST mission phases. The overall redistribution calculations were made assuming a specified initial contaminant distribution in terms of MIL STD 1246A levels, and using predicted vibration data for the various HST surfaces and mission phases. As expected, the effects of airflow were found to be significant, particularly for the larger particles. Particles smaller than about 20 microns did not participate appreciably in the redistribution.

Klavins, A.↗
NASA NTRS1988

Modeling of hydrogen-air diffusion flame

The present research objective is to determine the effects of contaminants on extinction limits of simple, well defined, counterflow Hydrogen 2-air diffusion flames, with combustion at 1 atmosphere. Results of extinction studies and other flame characterizations, with appropriate mechanistic modeling (presently underway), will be used to rationalize the observed effects of contamination over a reasonably wide range of diffusion flame conditions. The knowledge gained should help efforts to anticipate the effects of contaminants on combustion processes in Hydrogen 2-fueled scramjets.

Isaac, Kakkattukuzhy↗
NASA NTRS1992

A numerical study of Li-SF6 wick combustion - Forced and mixed convective burning

A numerical study is conducted to study Li-SF6 wick diffusion flames under mixed convective burning conditions at a pressure of 0.01 MPa. Both planar and cylindrical wicks are considered. The model is based on a conserved scalar approach. The objective of this study is to assess the effects of particular parameters on the burning rate and heat transfer. The flat-plate solution yields a fuel mass burning rate per unit surface area following the x exp -1/2 dependence of the classical similarity solution, where x is the streamwise distance. Cylindrical wick geometries yield enhanced burning rates over planar wicks. For the case of mixed convective burning, the burning rate results approach either the forced or natural convective burning limits as ambient streamwise velocity is changed. Critical Richardson numbers specifying these burning limits are determined for a given condition. Reducing gravity results in a lower burning rate because the influence of natural convection is diminished. Under reduced gravity of 1/1000 of the sea-level value, mixed convective burning nearly resembles forced convection.

Damaso, R. C.↗
NASA NTRS1993

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A computer code called Proteus 3D was developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Towne, Charles E.↗
NASA NTRS1993

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 3: Programmer's reference

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Towne, Charles E.↗
NASA NTRS2000

FDNS CFD Code Benchmark for RBCC Ejector Mode Operation: Continuing Toward Dual Rocket Effects

Computational Fluid Dynamics (CFD) analysis results are compared with benchmark quality test data from the Propulsion Engineering Research Center's (PERC) Rocket Based Combined Cycle (RBCC) experiments to verify fluid dynamic code and application procedures. RBCC engine flowpath development will rely on CFD applications to capture the multi -dimensional fluid dynamic interactions and to quantify their effect on the RBCC system performance. Therefore, the accuracy of these CFD codes must be determined through detailed comparisons with test data. The PERC experiments build upon the well-known 1968 rocket-ejector experiments of Odegaard and Stroup by employing advanced optical and laser based diagnostics to evaluate mixing and secondary combustion. The Finite Difference Navier Stokes (FDNS) code [2] was used to model the fluid dynamics of the PERC RBCC ejector mode configuration. Analyses were performed for the Diffusion and Afterburning (DAB) test conditions at the 200-psia thruster operation point, Results with and without downstream fuel injection are presented.

West, Jeff↗
NASA NTRS2014

Study of Aerothermodynamic Modeling Issues Relevant to High-Speed Sample Return Vehicles

This paper examines the application of state-of-the-art coupled ablation and radiation simulations to highspeed sample return vehicles, such as those returning from Mars or an asteroid. A defining characteristic of these entries is that the surface recession rates and temperatures are driven by nonequilibrium convective and radiative heating through a boundary layer with significant surface blowing and ablation products. Measurements relevant to validating the simulation of these phenomena are reviewed and the Stardust entry is identified as providing the best relevant measurements. A coupled ablation and radiation flowfield analysis is presented that implements a finite-rate surface chemistry model. Comparisons between this finite-rate model and a equilibrium ablation model show that, while good agreement is seen for diffusion-limited oxidation cases, the finite-rate model predicts up to 50% lower char rates than the equilibrium model at sublimation conditions. Both the equilibrium and finite rate models predict significant negative mass flux at the surface due to sublimation of atomic carbon. A sensitivity analysis to flowfield and surface chemistry rates show that, for a sample return capsule at 10, 12, and 14 km/s, the sublimation rates for C and C3 provide the largest changes to the convective flux, radiative flux, and char rate. A parametric uncertainty analysis of the radiative heating due to radiation modeling parameters indicates uncertainties ranging from 27% at 10 km/s to 36% at 14 km/s. Applying the developed coupled analysis to the Stardust entry results in temperatures within 10% of those inferred from observations, and final recession values within 20% of measurements, which improves upon the 60% over-prediction at the stagnation point obtained through an uncoupled analysis. Emission from CN Violet is shown to be over-predicted by nearly and order-of-magnitude, which is consistent with the results of previous independent analyses. Finally, the coupled analysis is applied to a 14 km/s Earth entry representative of a Mars sample return. Although the radiative heating provides a larger fraction of the total heating, the influence of ablation and radiation on the flowfield are shown to be similar to Stardust.

Johnston, Christopher O.↗
NASA NTRS2002

Simultaneous Raman-Rayleigh-LIF Measurements and Numerical Modeling Results of a Lifted H2/N2 Turbulent Jet Flame in a Vitiated Coflow

An experimental and numerical investigation is presented of a H2/N2 turbulent jet flame burner that has a novel vitiated coflow. The vitiated coflow emulates the recirculation region of most combustors, such as gas turbines or furnaces. Additionally, since the vitiated gases are coflowing, the burner allows for exploration of recirculation chemistry without the corresponding fluid mechanics of recirculation. Thus the vitiated coflow burner design facilitates the development of chemical kinetic combustion models without the added complexity of recirculation fluid mechanics. Scalar measurements are reported for a turbulent jet flame of H2/N2 in a coflow of combustion products from a lean ((empty set) = 0.25) H2/Air flame. The combination of laser-induced fluorescence, Rayleigh scattering, and Raman scattering is used to obtain simultaneous measurements of the temperature, major species, as well as OH and NO. Laminar flame calculation with equal diffusivity do agree when the premixing and preheating that occurs prior to flame stabilization is accounted for in the boundary conditions. Also presented is an exploratory pdf model that predicts the flame's axial profiles fairly well, but does not accurately predict the lift-off height.

Cabra, R.↗
NASA NTRS2000

Measurements and Modeling of Nitric Oxide Formation in Counterflow, Premixed CH4/O2/N2 Flames

Laser-induced fluorescence (LIF) measurements of NO concentration in a variety of CH4/O2/N2 flames are used to evaluate the chemical kinetics of NO formation. The analysis begins with previous measurements in flat, laminar, premixed CH4/O2/N2 flames stabilized on a water-cooled McKenna burner at pressures ranging from 1 to 14.6 atm, equivalence ratios from 0.5 to 1.6, and volumetric nitrogen/oxygen dilution ratios of 2.2, 3.1 and 3.76. These measured results are compared to predictions to determine the capabilities and limitations of the comprehensive kinetic mechanism developed by the Gas Research Institute (GRI), version 2.11. The model is shown to predict well the qualitative trends of NO formation in lean-premixed flames, while quantitatively underpredicting NO concentration by 30-50%. For rich flames, the model is unable to even qualitatively match the experimental results. These flames were found to be limited by low temperatures and an inability to separate the flame from the burner surface. In response to these limitations, a counterflow burner was designed for use in opposed premixed flame studies. A new LIF calibration technique was developed and applied to obtain quantitative measurements of NO concentration in laminar, counterflow premixed, CH4/O2/N2 flames at pressures ranging from 1 to 5.1 atm, equivalence ratios of 0.6 to 1.5, and an N2/O2 dilution ratio of 3.76. The counterflow premixed flame measurements are combined with measurements in burner-stabilized premixed flames and counterflow diffusion flames to build a comprehensive database for analysis of the GRI kinetic mechanism. Pathways, quantitative reaction path and sensitivity analyses are applied to the GRI mechanism for these flame conditions. The prompt NO mechanism is found to severely underpredict the amount of NO formed in rich premixed and nitrogen-diluted diffusion flames. This underprediction is traced to uncertainties in the CH kinetics as well as in the nitrogen oxidation chemistry. Suggestions are made which significantly improve the predictive capability of the GRI mechanism in near-stoichiometric, rich, premixed flames and in atmospheric-pressure, diffusion flames. However, the modified reaction mechanism is unable to model the formation of NO in ultra-rich, premixed or in high-pressure, nonpremixed flames, thus indicating the need for additional study under these conditions.

Thomsen, D. Douglas↗
NASA NTRS2001

Numerical Study of Buoyancy and Different Diffusion Effects on the Structure and Dynamics of Triple Flames

Numerical simulations of 2-D triple flames under gravity force have been implemented to identify the effects of gravity on triple flame structure and propagation properties and to understand the mechanisms of instabilities resulting from both heat release and buoyancy effects. A wide range of gravity conditions, heat release, and mixing widths for a scalar mixing layer are computed for downward-propagating (in the same direction with the gravity vector) and upward-propagating (in the opposite direction of the gravity vector) triple flames. Results of numerical simulations show that gravity strongly affects the triple flame speed through its contribution to the overall flow field. A simple analytical model for the triple flame speed, which accounts for both buoyancy and heat release, is developed. Comparisons of the proposed model with the numerical results for a wide range of gravity, heat release and mixing width conditions, yield very good agreement. The analysis shows that under neutral diffusion, downward propagation reduces the triple flame speed, while upward propagation enhances it. For the former condition, a critical Froude number may be evaluated, which corresponds to a vanishing triple flame speed. Downward-propagating triple flames at relatively strong gravity effects have exhibited instabilities. These instabilities are generated without any artificial forcing of the flow. Instead disturbances are initiated by minute round-off errors in the numerical simulations, and subsequently amplified by instabilities. A linear stability analysis on mean profiles of stable triple flame configurations have been performed to identify the most amplified frequency in spatially developed flows. The eigenfunction equations obtained from the linearized disturbance equations are solved using the shooting method. The linear stability analysis yields reasonably good agreements with the observed frequencies of the unstable triple flames. The frequencies and amplitudes of disturbances increase with the magnitude of the gravity vector. Moreover, disturbances appear to be most amplified just downstream of the premixed branches. The effects of mixing width and differential diffusion are investigated and their roles on the flame stability are studied.

Chen, Jyh-Yuan↗
NASA NTRS1997

Temperature and Radiative Heat Flux Measurements in Microgravity Jet Diffusion Flames

The objective of this project is to provide detailed measurements and modeling analyses of local soot concentration, temperature and radiation heat flux distributions in laminar and turbulent jet diffusion flames under normal (1-g) and reduced gravity (0-g) conditions. Results published to date by these co-PI's and their co-workers include: 1. thermophoretic sampling and size and morphological analyses of soot aggregates in laminar flames under normal and reduced gravity conditions; 2. full-field absorption imaging for soot volume fraction maps in laminar and turbulent flames under normal and reduced gravity conditions; 3. an accurate solver module for detailed radiation heat transfer in nongray nonhomogeneous media; 4. a complete model to include flame structure, soot formation and an energy equation to couple with radiation solver.

Ku, Jerry C.↗
NASA NTRS1993

Meso-Beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

The Control-B simulation experiment was designed to improve on the Control-A experiment performed with the GMASS model. This experiment addressed several inadequacies with the first smooth terrain numerical simulation by including: (1) increased nested-grid resolution to better define the simulated gravity waves, (2) increased horizontal diffusion to remove outflow boundary condition noise, and (3) the use of reanalyzed rawinsonde data and surface observations in the initial state to increase the definition of the observed jet streak as well as other low-level features. A smoothed-terrain dry simulation with the nested-grid GMASS model has revealed many important aspects of the processes which resulted in the generation of gravity waves in the region and time when and where they were observed. However, the vertical structure, number, and characteristics of the waves are still quite different from observed waves as diagnosed thus necessitating future improved simulations. However, this control simulation has produced substantial insight into processes which occur on many spatial scales over a 30 hour time period thus allowing one to draw promising inferences as to the mechanisms for the complex process which occurred in nature during the CCOPE case study. The theoretical aspects of the project have focus on understanding the nature of the ageostrophic circulations which are produced in idealized models of the atmosphere in which the troposphere is modeled in one of two ways. The first model assumes that the lower atmosphere can be represented as a single layer of homogeneous fluid whose upper surface is free to exhibit vertical displacement. Two-dimensional internal convergence (divergence) occurring during the adjustment to an asymptotic equilibrium state from an ageostrophic initial state whose momentum structure is representative of a midlatitude localized zonal wind anomaly will cause the free upper surface of the homogeneous atmosphere to rise (fall), and therefore the response can be viewed as being physically three-dimensional. The second model assumes that the troposphere can be represented by an unbounded continuously stratified Boussinesq fluid of constant Brunt-Vaisala frequency N, where the vertical gradient of the basic state potential temperature profile allows for the existence of vertically propagating internal inertia-gravity waves.

Lin, Yuh-Lang↗
NASA NTRS1985

The perturbed neutral circulation in the vicinity of a symmetric stable auroral arc

In the present investigation, a sophisticated high-resolution time-dependent numerical model is used to study the neutral response to changes in the E region electrodynamical forcing and particle heating associated with the appearance of a stable auroral arc. The initial condition is a nighttime ionosphere in the presence of a background precipitation of electrons associated with diffuse auroras. A simulation study is conducted of the perturbation in the neutral circulation. The simulations were performed on the basis of a model which describes nonlinear, nonhydrostatic viscous flow on a rotating atmosphere.

Walterscheid, R. L.↗
NASA NTRS1972

The effects of boundary condition asymmetries on the interplanetary magnetic field-moon interaction.

Boundary condition asymmetries inherent in the solar wind flow past the moon are included in a cylindrical model of the interplanetary magnetic field-moon interaction. Numerical examinations of the sunward side response of this model are compared in the frequency domain with those of symmetrically excited spherical and cylindrical models and two characteristic differences are observed: the response of the asymmetric model is depressed at low frequencies due to magnetic diffusion around a conducting core, and is flattened at high frequencies because of the finite application time of the incident interplanetary magnetic field. The diffusion of field lines around the core is also evident in the time response of the model in the antisolar cavity. The above features of the lunar response resulting from boundary condition asymmetries are shown to be evident in observational measurements.

Reisz, A. C.↗
NASA NTRS1972

Investigation of bubbles in arterial heat pipes

The behavior of gas occlusions in arterial heat pipes has been studied experimentally and theoretically. Specifically, the gas-liquid system properties, solubility and diffusivity, have been measured from -50 to 100 C for helium and argon in ammonia, Freon-21 (CHC12F), and methanol. Properties values obtained were then used to experimentally test models for gas venting from a heat pipe artery under isothermal conditions (i.e., no-heat flow), although the models, as developed, are also applicable to heat pipes operated at power, with some minor modifications. Preliminary calculations indicated arterial bubbles in a stagnant pipe require from minutes to days to collapse and vent. It has been found experimentally that a gas bubble entrapped within an artery structure has a very long lifetime in many credible situations. This lifetime has an approximately inverse exponential dependence on temperature, and is generally considerably longer for helium than for argon. The models postulated for venting under static conditions were in general quantitative agreement with experimental data. Factors of primary importance in governing bubble stability are artery diameter, artery wall thickness, noncondensible gas partial pressure, and the property group (the Ostwald solubility coefficient multiplied by the gas/liquid diffusivity).

Saaski, E. W.↗
NASA NTRSDate not supplied

Ice Growth Model for Polar Water Capture Systems

In-situ resource utilization (ISRU) is essential to NASA’s goal of sustainable and cost-effective missions to the Moon and beyond. Water has become a particular source of interest since its discovery in the permanently shadowed regions (PSRs) of the Moon. This water exists as ice bound within the regolith, so one mode of capture is heating to sublimation temperatures and capturing at very low pressure. It then can be stored in ice form and transported out of the PSR and to a sunlit ridge, where processes such as electrolysis can be performed[2]. An understanding of frost growth dynamics under rarefied conditions is necessary for successful water capture and storage due to transient thermal properties of frost. A Diffusion-Limited Aggregation (DLA) approach was taken to analytically model bulk thermal conductivity changes of the growing frost layer under varying wall temperatures and frost layer porosities. The model utilizes well-known correlations for effective thermal conductivity and records the values as the frost layer thickness increases. This thermal conductivity can be applied to a known temperature change and be used to evaluate heat flux through the combined frost layer and tank wall for either a flat plate or cylindrical geometry. An experiment is described that will evaluate this model using one-dimensional frost growth on a horizontal, flat plate under low-pressure. The ice growth will be measured using a millimeter-scale photogrammetry camera system and a load cell measured transient mass growth.

In situ Resource Utilization↗