SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

NASA NTRS2007

Transient Numerical Modeling of Catalytic Channels

This paper presents a transient model of catalytic combustion suitable for isolated channels and monolith reactors. The model is a lumped two-phase (gas and solid) model where the gas phase is quasi-steady relative to the transient solid. Axial diffusion is neglected in the gas phase; lateral diffusion, however, is accounted for using transfer coefficients. The solid phase includes axial heat conduction and external heat loss due to convection and radiation. The combustion process utilizes detailed gas and surface reaction models. The gas-phase model becomes a system of stiff ordinary differential equations while the solid phase reduces, after discretization, into a system of stiff ordinary differential-algebraic equations. The time evolution of the system came from alternating integrations of the quasi-steady gas and transient solid. This work outlines the numerical model and presents some sensitivity studies on important parameters including internal transfer coefficients, catalytic surface site density, and external heat-loss (if applicable). The model is compared to two experiments using CO fuel: (1) steady-state conversion through an isothermal platinum (Pt) tube and (2) transient propagation of a catalytic reaction inside a small Pt tube. The model requires internal mass-transfer resistance to match the experiments at lower residence times. Under mass-transport limited conditions, the model reasonably predicted exit conversion using global mass-transfer coefficients. Near light-off, the model results did not match the experiment precisely even after adjustment of mass-transfer coefficients. Agreement improved for the first case after adjusting the surface kinetics such that the net rate of CO adsorption increased compared to O2. The CO / O2 surface mechanism came from a sub-set of reactions in a popular CH4 / O2 mechanism. For the second case, predictions improved for lean conditions with increased external heat loss or adjustment of the kinetics as in the first case. Finally, the results show that different initial surface-species distribution leads to different steady-states under certain conditions. These results demonstrate the utility of a lumped two-phase model of a transient catalytic combustor with detailed chemistry.

Struk, Peter M.↗

NASA NTRS1990

ULTRA-SHARP nonoscillatory convection schemes for high-speed steady multidimensional flow

For convection-dominated flows, classical second-order methods are notoriously oscillatory and often unstable. For this reason, many computational fluid dynamicists have adopted various forms of (inherently stable) first-order upwinding over the past few decades. Although it is now well known that first-order convection schemes suffer from serious inaccuracies attributable to artificial viscosity or numerical diffusion under high convection conditions, these methods continue to enjoy widespread popularity for numerical heat transfer calculations, apparently due to a perceived lack of viable high accuracy alternatives. But alternatives are available. For example, nonoscillatory methods used in gasdynamics, including currently popular TVD schemes, can be easily adapted to multidimensional incompressible flow and convective transport. This, in itself, would be a major advance for numerical convective heat transfer, for example. But, as is shown, second-order TVD schemes form only a small, overly restrictive, subclass of a much more universal, and extremely simple, nonoscillatory flux-limiting strategy which can be applied to convection schemes of arbitrarily high order accuracy, while requiring only a simple tridiagonal ADI line-solver, as used in the majority of general purpose iterative codes for incompressible flow and numerical heat transfer. The new universal limiter and associated solution procedures form the so-called ULTRA-SHARP alternative for high resolution nonoscillatory multidimensional steady state high speed convective modelling.

Leonard, B. P.↗

NASA NTRS2000

Design and Performance of McRas in SCMs and GEOS I/II GCMs

The design of a prognostic cloud scheme named McRAS (Microphysics of clouds with Relaxed Arakawa-Schubert Scheme) for general circulation models (GCMs) will be discussed. McRAS distinguishes three types of clouds: (1) convective, (2) stratiform, and (3) boundary-layer types. The convective clouds transform and merge into stratiform clouds on an hourly time-scale, while the boundary-layer clouds merge into the stratiform clouds instantly. The cloud condensate converts into precipitation following the auto-conversion equations of Sundqvist that contain a parametric adaptation for the Bergeron-Findeisen process of ice crystal growth and collection of cloud condensate by precipitation. All clouds convect, advect, as well as diffuse both horizontally and vertically with a fully interactive cloud-microphysics throughout the life-cycle of the cloud, while the optical properties of clouds are derived from the statistical distribution of hydrometeors and idealized cloud geometry. An evaluation of McRAS in a single column model (SCM) with the GATE Phase III and 5-ARN CART datasets has shown that together with the rest of the model physics, McRAS can simulate the observed temperature, humidity, and precipitation without many systematic errors. The time history and time mean incloud water and ice distribution, fractional cloudiness, cloud optical thickness, origin of precipitation in the convective anvil and towers, and the convective updraft and downdraft velocities and mass fluxes all show a realistic behavior. Performance of McRAS in GEOS 11 GCM shows several satisfactory features but some of the remaining deficiencies suggest need for additional research involving convective triggers and inhibitors, provision for continuously detraining updraft, a realistic scheme for cumulus gravity wave drag, and refinements to physical conditions for ascertaining cloud detrainment level.

Sud, Yogesh C.↗

NASA NTRS2017

Prediction of Turbulent Temperature Fluctuations in Hot Jets

Large-eddy simulations (LES) were used to investigate turbulent temperature fluctuations and turbulent heat flux in hot jets. A high-resolution finite-difference Navier-Stokes solver was used to compute the flow from a 2-inch round nozzle. Three different flow conditions of varying jet Mach numbers and temperature ratios were examined. The LES results showed that the temperature field behaves similar to the velocity field, but with a more rapidly spreading mixing layer. Predictions of mean, mu-bar(sub i), and fluctuating, mu'(sub i), velocities were compared to particle image velocimetry data. Predictions of mean, T-bar, and fluctuating, T', temperature were compared to data obtained using Rayleigh scattering and Raman spectroscopy. Very good agreement with experimental data was demonstrated for the mean and fluctuating velocities. The LES correctly predicts the behavior of the turbulent temperature field, but over-predicts the levels of the fluctuations. The turbulent heat flux was examined and compared to Reynolds-averaged Navier-Stokes (RANS) results. The LES and RANS simulations produced very similar results for the radial heat flux. However, the axial heat flux obtained from the LES differed significantly from the RANS result in both structure and magnitude, indicating that the gradient diffusion type model in RANS is inadequate. Finally, the LES data was used to compute the turbulent Prandtl number and verify that a constant value of 0.7 used in the RANS models is a reasonable assumption.

hot jets↗

NASA NTRS1999

Effects of Buoyancy in Hydrogen Jet Diffusion Flames

This project was carried out to understand the effects of heat release and buoyancy on the flame structure of diffusion flames. Experiments were conducted at atmospheric pressure in both normal gravity and microgravity conditions in the NASA LeRC 2.2 s drop tower. Experiments were also conducted in a variable pressure combustion facility in normal gravity to scale buoyancy and thus, to supplement the drop tower experiments. Pure H2 or H2 mixed with He was used as the jet fluid to avoid the complexities associated with soot formation. Fuel jet burning in quiescent air was visualized and quantified by the Rainbow Schlieren Deflectometry (RSD) to obtain scalar profiles (temperature, oxygen concentration) within the flame. Burner tube diameter (d) was varied from 0.3 to 1.19 mm producing jet exit Reynolds numbers ranging from 40 to 1900, and generating flames encompassing laminar and transitional (laminar to turbulent) flow structure. Some experiments were also complemented with the CFD analysis. In a previous paper, we have presented details of the RSD technique, comparison of computed and measured scalar distributions, and effects of buoyancy on laminar and transitional H2 gas-jet diffusion flames. Results obtained from the RSD technique, variable pressure combustion chamber, and theoretical models have been published. Subsequently, we have developed a new drop rig with improved optical and image acquisition. In this set up, the schlieren images are acquired in real time and stored digitally in RAM of an onboard computer. This paper deals with laminar diffusion flames of pure H2 in normal and microgravity.

Agrawal, A. K.↗

NASA NTRS2003

Jupiter System Data Analysis Program: Mechanisms, Manifestation, and Implications of Cryomagmatism on Europa

The objectives of the work completed under NASA Grant NAG5-8898 were (i) to document and characterize the low-albedo diffuse surfaces associated with triple bands and lenticulae, (ii) to determine their mechanisms of formation, and (iii) to assess the implications of these features for the resurfacing (in space and time) of the Europa and the nature of the Europan interior. The approach involved a combination of processing and analysis of Solid State Imaging data returned by the Galileo spacecraft during the primary and extended mission phases, together with numerical modeling of the physical processes interpreted to the observed features. We have modeled the formation of Europan triple explosive venting of cryoclastic material from bands and lenticulae halos by two processes: (i) a liquid layer in the Europan interior, and (ii) lag deposit formation by the thermal influence of subsurface cryomagmatic intrusions. We favor the latter hypothesis for explaining these features, and further suggest that a liquid water or brine intrusion is required to provide sufficient lateral heating of surface ice to explain the 25 km size of the largest features. (Solid ice diapirs, even under the most favorable conditions, become thermally exhausted before they heat significant lateral distances). We argue that water circulating in open fractures, or repeated cryomagmatic 'diking' events would provide sufficient thermal input to produce the observed features. Thus our work argues for the existence of a liquid beneath Europa's surface. Our results might most easily be explained by the presence of a continuous liquid layer (the putative Europan ocean); this would concur with the findings of the Galileo magnetometer team. However, we cannot rule out the possibility that discrete liquid pockets provide injections of fluid closer to the surface.

Fagents, Sarah A.↗

NASA NTRS2003

First Year Wilkinson Microwave Anisotropy Probe(WMAP)Observations: The Angular Power Spectrum

We present the angular power spectrum derived from the first-year Wilkinson Microwave Anisotropy Probe (WMAP) sky maps. We study a variety of power spectrum estimation methods and data combinations and demonstrate that the results are robust. The data are modestly contaminated by diffuse Galactic foreground emission, but we show that a simple Galactic template model is sufficient to remove the signal. Point sources produce a modest contamination in the low frequency data. After masking approximately 700 known bright sources from the maps, we estimate residual sources contribute approximately 3500 mu sq Kappa at 41 GHz, and approximately 130 mu sq Kappa at 94 GHz, to the power spectrum [iota(iota + 1)C(sub iota)/2pi] at iota = 1000. Systematic errors are negligible compared to the (modest) level of foreground emission. Our best estimate of the power spectrum is derived from 28 cross-power spectra of statistically independent channels. The final spectrum is essentially independent of the noise properties of an individual radiometer. The resulting spectrum provides a definitive measurement of the CMB power spectrum, with uncertainties limited by cosmic variance, up to iota approximately 350. The spectrum clearly exhibits a first acoustic peak at iota = 220 and a second acoustic peak at iota approximately 540, and it provides strong support for adiabatic initial conditions. Researchers have analyzed the CT(sup Epsilon) power spectrum, and present evidence for a relatively high optical depth, and an early period of cosmic reionization. Among other things, this implies that the temperature power spectrum has been suppressed by approximately 30% on degree angular scales, due to secondary scattering.

Hinshaw, G.↗

NASA NTRS2013

Chemical Enrichment and Physical Conditions in IZw18*

Low-metallicity star-forming dwarf galaxies are prime targets to understand the chemical enrichment of the interstellar medium. The H I region contains the bulk of the mass in blue compact dwarfs, and it provides important constraints on the dispersal and mixing of heavy elements released by successive star-formation episodes. The metallicity of the H I region is also a critical parameter to investigate the future star-formation history, as metals provide most of the gas cooling that will facilitate and sustain star formation. Aims. Our primary objective is to study the enrichment of the H I region and the interplay between star-formation history and metallicity evolution. Our secondary obje ctive is to constrain the spatial- and time-scales over which the HI and H II regions are enriched, and the mass range of stars responsible for the heavy element production. Finally, we aim to examine the gas heating and cooling mechanisms in the H I region. Methods. We observed the most metal-poor star-forming galaxy in the Local Universe, I Zw 18, with the Cosmic Origin Spectrograph onboard Hubble. The abundances in the neutral gas are derived from far-ultraviolet absorption-lines (H I, CIII, CIIi*, N I, OI,...) and are compared to the abundances in the H II region. Models are constructed to calculate the ionization structure and the thermal processes. We investigate the gas cooling in the HIi region through physical diagnostics drawn from the fine-structure level of C+. Results. We find that H I region abundances are lower by a factor of approx 2 as compared to the H II region. There is no differential depletion on dust between the H I and H II region. Using sulfur as a metallicity tracer, we calculate a metallicity of 1/46 Z(solar) (vs. 1/31 Z(solar) in the H II region). From the study of the C/O, [O/Fe], and N/O abundance ratios, we propose that C, N, O, and Fe are mainly produced in massive stars. We argue that the H I envelope may contain pockets of pristine gas with a metallicity essentially null. Finally, we derive the physical conditions in the H I region by investigating the C II* absorption line. The cooling rate derived from C II* is consistent with collisions with H(sup 0) atoms in the diffuse neutral gas. We calculate the star-formation rate from the CII* cooling rate assuming that photoelectric effect on dust is the dominant gas heating mechanism. Our determination is in good agreement with the values in the literature if we assume a low dust-to-gas ratio (approx 2000 times lower than the Milky Way value).

far-ultraviolet absorption-lines↗

NASA NTRS1993

Containerless measurements on liquids at high temperatures

The application of containerless techniques for measurements of the thermophysical properties of high temperature liquids is reviewed. Recent results obtained in the materials research laboratories at Intersonics are also presented. Work to measure high temperature liquid properties is motivated by both the need for reliable property data for modeling of industrial processes involving molten materials and generation of data form basic modeling of materials behavior. The motivation for this work and examples of variations in thermophysical property values from the literature are presented. The variations may be attributed to changes in the specimen properties caused by chemical changes in the specimen and/or to measurement errors. The two methods used to achieve containerless conditions were aeroacoustic levitation and electromagnetic levitation. Their qualities are presented. The accompanying slides show the layout of levitation equipment and present examples of levitated metallic and ceramic specimens. Containerless techniques provide a high degree of control over specimen chemistry, nucleation and allow precise control of liquid composition to be achieved. Effects of minor additions can thus be measured in a systematic way. Operation in reduced gravity enables enhanced control of liquid motion which can allow measurement of liquid transport properties. Examples of nucleation control, the thermodynamics of oxide contamination removal, and control of the chromium content of liquid aluminum oxide by high temperature containerless processes are presented. The feasibility of measuring temperature, emissivity, liquidus temperature, enthalpy, surface tension, density, viscosity, and thermal diffusivity are discussed in the final section of the paper.

Weber, Richard↗

NASA NTRS2003

Fundamental space radiobiology

The unique feature of the space radiation environment is the dominance of high-energy charged particles (HZE or high LET radiation) emitted by the Sun and galactic sources, or trapped in the Van Allen radiation belts. These charged particles present a significant hazard to space flight crews, and accelerator-based experiments are underway to quantify the health risks due to unavoidable radiation exposure. There are three principal properties of charged particles that distinguish them from conventional radiation, i.e. gamma rays and x-rays. First, they have a defined range in matter rather than an exponential absorption profile. Second, they undergo nuclear reactions to produce secondary particles. Third, and most important, they deposit their energy along well-defined linear paths or tracks rather than diffuse fields. The structured energy deposition pattern interacts on multiple scales with the biological structures of DNA, cells and tissues to produce correlated patterns of damage that evade repair systems. Traditional concepts of dose and its associated normalization parameter, RBE (relative biological effectiveness), break down under experimental scrutiny, and probabilistic models of risk based on the number of particle traversals per cell may be more appropriate. Unique patterns of DNA damage, gene expression, mobilization of repair proteins, activation of cytokines and remodeling of cellular microenvironment are observed following exposure to high LET radiation. At low levels of exposure the communication of bioactive substances from irradiated to unirradiated "bystander" cells can amplify the damage and cause a significant deviation from linearity in dose vs. response relations. Under some circumstances, there is even a multigenerational delay in the expression of radiation-induced genetic damage (genomic instability) which is not strictly dose dependent. These issues and the experimental evidence derived from ground based experiments at particle accelerators are presented along with speculation about how modified inertial conditions might perturb homeostatic responses to radiation to further complicate risk assessment for space flight.

NASA Discipline Radiation Health↗

NASA NTRS2008

Multiscale Modeling of Damage Processes in Aluminum Alloys: Grain-Scale Mechanisms

This paper has two goals related to the development of a physically-grounded methodology for modeling the initial stages of fatigue crack growth in an aluminum alloy. The aluminum alloy, AA 7075-T651, is susceptible to fatigue cracking that nucleates from cracked second phase iron-bearing particles. Thus, the first goal of the paper is to validate an existing framework for the prediction of the conditions under which the particles crack. The observed statistics of particle cracking (defined as incubation for this alloy) must be accurately predicted to simulate the stochastic nature of microstructurally small fatigue crack (MSFC) formation. Also, only by simulating incubation of damage in a statistically accurate manner can subsequent stages of crack growth be accurately predicted. To maintain fidelity and computational efficiency, a filtering procedure was developed to eliminate particles that were unlikely to crack. The particle filter considers the distributions of particle sizes and shapes, grain texture, and the configuration of the surrounding grains. This filter helps substantially reduce the number of particles that need to be included in the microstructural models and forms the basis of the future work on the subsequent stages of MSFC, crack nucleation and microstructurally small crack propagation. A physics-based approach to simulating fracture should ultimately begin at nanometer length scale, in which atomistic simulation is used to predict the fundamental damage mechanisms of MSFC. These mechanisms include dislocation formation and interaction, interstitial void formation, and atomic diffusion. However, atomistic simulations quickly become computationally intractable as the system size increases, especially when directly linking to the already large microstructural models. Therefore, the second goal of this paper is to propose a method that will incorporate atomistic simulation and small-scale experimental characterization into the existing multiscale framework. At the microscale, the nanoscale mechanics are represented within cohesive zones where appropriate, i.e. where the mechanics observed at the nanoscale can be represented as occurring on a plane such as at grain boundaries or slip planes at a crack front. Important advancements that are yet to be made include: 1. an increased fidelity in cohesive zone modeling; 2. a means to understand how atomistic simulation scales with time; 3. a new experimental methodology for generating empirical models for CZMs and emerging materials; and 4. a validation of simulations of the damage processes at the nano-micro scale. With ever-increasing computer power, the long-term ability to employ atomistic simulation for the prognosis of structural components will not be limited by computation power, but by our lack of knowledge in incorporating atomistic models into simulations of MSFC into a multiscale framework.

Hochhalter, J. D.↗

NASA NTRS1986

Ground testing in a simulated auroral environment

Shuttle polar orbit missions are being planned in which astronaut extravehicular activities (EVA) will eventually include external operations or repair and refurbishment of satellites. Recent satellite data and model calculations show that significant spacecraft charging and differential charging occurs in that environment and could endanger equipment and personnel involved in the EVA. These spacecraft charging events take place in response to an intense flux of high energy electrons accompanied by a large drop in the ambient plasma density. The extent of which the auroral environment can be simulated in a large vacuum chamber are determined so that conditions which lead to significant charging and discharging may be observed and controlled. A large chamber would be required in order to simulate the effects of an aurora-like, large cross-section, energetic electron precipitation on a target system the size of an astronaut with a life support system. The precipitating electrons can be simulated with an array of multipactor sources and the background ionization produced by the beam can be limited to the range of the auroral density if the chamber pressure is not significantly above 1 x .000001 Torr. Some small chamber work has shown that the plasma in the beam will diffuse radially at a rate determined by the Bohm diffusion coefficient.

Mcintyre, Bernard↗

NASA NTRS1963

Performance Capability of Single-Cavity Vortex Gaseous Nuclear Rockets

An analysis was made to determine the maximum powerplant thrust-to-weight ratio possible with a single-cavity vortex gaseous reactor in which all the hydrogen propellant must diffuse through a fuel-rich region. An assumed radial temperature profile was used to represent conduction, convection, and radiation heat-transfer effects. The effect of hydrogen property changes due to dissociation and ionization was taken into account in a hydrodynamic computer program. It is shown that, even for extremely optimistic assumptions of reactor criticality and operating conditions, such a system is limited to reactor thrust-to-weight ratios of about 1.2 x 10(exp -3) for laminar flow. For turbulent flow, the maximum thrust-to-weight ratio is less than 10(exp -3). These low thrusts result from the fact that the hydrogen flow rate is limited by the diffusion process. The performance of a gas-core system with a specific impulse of 3000 seconds and a powerplant thrust-to-weight ratio of 10(exp -2) is shown to be equivalent to that of a 1000-second advanced solid-core system. It is therefore concluded that a single-cavity vortex gaseous reactor in which all the hydrogen must diffuse through the nuclear fuel is a low-thrust device and offers no improvement over a solid-core nuclear-rocket engine. To achieve higher thrust, additional hydrogen flow must be introduced in such a manner that it will by-pass the nuclear fuel. Obviously, such flow must be heated by thermal radiation. An illustrative model of a single-cavity vortex system employing supplementary flow of hydrogen through the core region is briefly examined. Such a system appears capable of thrust-to-weight ratios of approximately 1 to 10. For a high-impulse engine, this capability would be a considerable improvement over solid-core performance. Limits imposed by thermal radiation heat transfer to cavity walls are acknowledged but not evaluated. Alternate vortex concepts that employ many parallel vortices to achieve higher hydrogen flow rates offer the possibility of sufficiently high thrust-to-weight ratios, if they are not limited by short thermal-radiation path lengths.

Ragsdale, Robert G.↗

NASA NTRS2000

Numerical Simulation of a Chemically Reacting Sorbent Bed for LSS Applications

A detailed numerical model of a chemisorption bed has been developed. The model is based on the constant pressure mass transport equation for gaseous flow through a packed bed, and the equation for diffusion and reaction within a spherical particle. Because there is a wealth of data from the NASA and the Navy bodies of literature, the LiOH-H2O-CO2 system is chosen for application of the model and interpretation of results. Prior models of this system from the life support literature are limited. The current model incorporates many of the features of elaborate models developed for investigation of industrial systems or energy applications (e.g., coal, desulphurization): it distinguishes bulk convection and bed dispersion; mass transport to the particle surface, transport within the particle, and reaction. It uses the nonsteady (not pseudo-steady state) form of the equations. The chemistry is modeled as a multi-step, reversible reaction with evolving solid structure. The resulting system of equations is large. The ODEPACK family of solvers is used to integrate the system. Reaction coefficients are determined by experiment. Typical results of the model are illustrated with mission input parameters. Using the model, an explanation is offered for 1) the varied performance results found after pre-breathing (or after simulated pre-breathe conditions), 2) interrupted use and 3) low temperature use. In addition, options for a reusable canister are explored. The computational resource implications of adding energy equations are discussed briefly, as are applicability to other relevant space and undersea systems.

Luna, Bernadette↗

NASA NTRS2003

Modeling Candle Flame Behavior In Variable Gravity

The burning of a candle, as typical non-propagating diffusion flame, has been used by a number of researchers to study the effects of electric fields on flame, spontaneous flame oscillation and flickering phenomena, and flame extinction. In normal gravity, the heat released from combustion creates buoyant convection that draws oxygen into the flame. The strength of the buoyant flow depends on the gravitational level and it is expected that the flame shape, size and candle burning rate will vary with gravity. Experimentally, there exist studies of candle burning in enhanced gravity (i.e. higher than normal earth gravity, g(sub e)), and in microgravity in drop towers and space-based facilities. There are, however, no reported experimental data on candle burning in partial gravity (g < g(sub e)). In a previous numerical model of the candle flame, buoyant forces were neglected. The treatment of momentum equation was simplified using a potential flow approximation. Although the predicted flame characteristics agreed well with the experimental results, the model cannot be extended to cases with buoyant flows. In addition, because of the use of potential flow, no-slip boundary condition is not satisfied on the wick surface. So there is some uncertainty on the accuracy of the predicted flow field. In the present modeling effort, the full Navier-Stokes momentum equations with body force term is included. This enables us to study the effect of gravity on candle flames (with zero gravity as the limiting case). In addition, we consider radiation effects in more detail by solving the radiation transfer equation. In the previous study, flame radiation is treated as a simple loss term in the energy equation. Emphasis of the present model is on the gas-phase processes. Therefore, the detailed heat and mass transfer phenomena inside the porous wick are not treated. Instead, it is assumed that a thin layer of liquid fuel coated the entire wick surface during the burning process. This is the limiting case that the mass transfer process in the wick is much faster than the evaporation process at the wick surface.

Alsairafi, A.↗

NASA NTRS1994

Nonaxisymmetric evolution in protostellar disks

We present a two-dimensional, multigridded hydrodynamical simulation of the collapse of an axisymmetric, rotating, 1 solar mass protostellar cloud, which forms a resolved, hydrotastic disk. The code includes the effects of physical viscosity, radiative transfer and radiative acceleration but not magnetic fields. We examine how the disk is affected by the inclusion of turbulent viscosity by comparing a viscous simulation with an inviscid model evolved from the same initial conditions, and we derive a disk evolutionary timescale on the order of 300,000 years if alpha = 0.01. Effects arising from non-axisymmetric gravitational instabilities in the protostellar disk are followed with a three-dimensional SPH code, starting from the two-dimensional structure. We find that the disk is prone to a series of spiral instabilities with primary azimulthal mode number m = 1 and m = 2. The torques induced by these nonaxisymmetric structures elicit material transport of angular momentum and mass through the disk, readjusting the surface density profile toward more stable configurations. We present a series of analyses which characterize both the development and the likely source of the instabilities. We speculate that an evolving disk which maintains a minimum Toomre Q-value approximately 1.4 will have a total evolutionary span of several times 10(exp 5) years, comparable to, but somewhat shorter than the evolutionary timescale resulting from viscous turbulence alone. We compare the evolution resulting from nonaxisymmetric instabilities with solutions of a one-dimensional viscous diffusion equation applied to the initial surface density and temperature profile. We find that an effective alpha-value of 0.03 is a good fit to the results of the simulation. However, the effective alpha will depend on the minimum Q in the disk at the time the instability is activated. We argue that the major fraction of the transport characterized by the value of alpha is due to the action of gravitational torques, and does not arise from inherent viscosity within the smoothed particle hydrodynamics method.

Laughlin, Gregory↗

NASA NTRS1999

Experiments on Diffusion Flame Structure of a Laminar Vortex Ring

The study of flame-vortex interactions provides one of the means to better understand turbulent combustion, and allows for canonical configurations that contain the fundamental elements found in turbulent flames, These include concentrated vorticity, entrainment and mixing, strain and nonequilibrium phenomena, diffusion and differential diffusion, partial premixing and diluent effects, and heat release effects. In flame- vortex configurations, these fundamental elements can be studied under more controlled conditions than is possible in direct investigations of turbulent flames. Since the paper of Marble, the problem of the flame-vortex interaction has received considerable attention theoretically, numerically and experimentally. Several configurations exist for study of the premixed flame/vortex ring interaction but more limited results have been obtained to date for the diffusion flame/vortex ring case. The setup of Chen and Dahm, which is conceptually similar to that of Karagozian and Manda and Karagozian, Suganuma and Strom where the ring is composed of fuel and air and combustion begins during the ring formation process, is used in the current study. However, it is essential to conduct the experiments in microgravity to remove the asymmetries caused by buoyancy and thus obtain highly symmetric and repeatable interactions. In previous studies it was found that the flame structure of the vortex ring was similar to that obtained analytically by Karagozian and Manda. Dilution of propane with nitrogen led mainly to a reduction in flame luminosities, flame burnout times were affected by both fuel volumes and amount of dilution, and a simple model of the burnout times was developed. In this paper, a discussion on reacting ring displacement and flame burnout time will be given, and the flame structures of vortex rings containing ethane and air will be compared to those of propane reacting in air.

Chen, Shin-Juh↗

NASA NTRS1994

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.↗