Engineering Papers⌕ Search

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.

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

At least 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

The Effect of Diffusional Transport and Surface Catalysis on the Aerothermodynamic Heating for Martian Atmospheric Entry

For the reentry heating of 70-degree blunt sphere-cones entering Mars at velocities in excess of 6 km/sec, a study was performed to determine the magnitude of both convective and catalytic heating on the cone forebody. Case studies of the peak heating conditions for the current NASA Mars Pathfinder entry probe, as well as the peak heating conditions of the proposed NASA Mars 2001 aerobraking orbiter mission were performed. Since the actual behavior of the chemical recombination of Mars atmosphere hypersonic shock layer species on heat shield materials is not known, and no experimental studies of any consequence have been done to determine surface reaction rates for the CO2 system, a parametric analysis of various reaction schemes and surface catalytic mechanistic models was performed. In many cases the actual Mars entry probe will have a heat shield composed of an ablative material which can either partially or completely preclude the existence of catalytic surface recombination phenomena. The extent of this blockage effect has not been examined at this time and is not considered in this effort. Instead only non-ablative computations are performed to separately evaluate the full extent of the surface catalysis and related diffusional phenomena. Parametric studies include peak heating point comparisons of non-catalytic, partially catalytic and fully catalytic total surface heat transfer for three difference surface recombination catalytic mechanisms (with and without simultaneous CO + O and O + O recombination) as well as with different species diffusion models. Diffusion model studies include constant Schmidt number (equal diffusivities) and the bifurcation, approximate multi-component diffusion model. Since the gas phase reaction kinetics for the Mars CO2 based atmosphere have also not been validated, the effect of two different (fast and slow) sets of gas kinetics on heat transfer is presented.

Henline, William D.↗
NASA NTRS2018

Satellite Radiation Products for Ocean Biology and Biogeochemistry: Needs, State-Of-The-Art, Gaps, Development Priorities, and Opportunities

Knowing the spatial and temporal distribution of the underwater light field, i.e., the spectral and angular structure of the radiant intensity at any point in the water column, is essential to understanding the biogeochemical processes that control the composition and evolution of aquatic ecosystems and their impact on climate and reaction to climate change. At present, only a few properties are reliably retrieved from space, either directly or via water-leaving radiance. Existing satellite products are limited to planar photosynthetically available radiation (PAR) and ultraviolet (UV) irradiance above the surface and diffuse attenuation coefficient. Examples of operational products are provided, and their advantages and drawbacks are examined. The usefulness and convenience of these products notwithstanding, there is a need, as expressed by the user community, for other products, i.e., sub-surface planar and scalar fluxes, average cosine, spectral fluxes (UV to visible), diurnal fluxes, absorbed fraction of PAR by live algae (APAR), surface albedo, vertical attenuation, and heating rate, and for associating uncertainties to any product on a pixel-by-pixel basis. Methodologies to obtain the new products are qualitatively discussed in view of most recent scientific knowledge and current and future satellite missions, and specific algorithms are presented for some new products, namely sub-surface fluxes and average cosine. A strategy and roadmap (short, medium, and long term) for usage and development priorities is provided, taking into account needs and readiness level. Combining observations from satellites overpassing at different times and geostationary satellites should be pursued to improve the quality of daily-integrated radiation fields, and products should be generated without gaps to provide boundary conditions for general circulation and biogeochemical models. Examples of new products, i.e., daily scalar PAR below the surface, daily average cosine for PAR, and sub-surface spectral scalar fluxes are presented. A procedure to estimate algorithm uncertainties in the total uncertainty budget for above-surface daily PAR, based on radiative simulations for expected situations, is described. In the future, space-borne lidars with ocean profiling capability offer the best hope for improving our knowledge of sub-surface fields. To maximize temporal coverage, space agencies should consider placing ocean-color instruments in L1 orbit, where the sunlit part of the Earth can be frequently observed.

Frouin, Robert↗
NASA NTRS2014

Xenon Fractionation, Hydrogen Escape, and the Oxidation of the Earth

Xenon in Earth's atmosphere is severely mass fractionated and depleted compared to any plausible solar system source material, yet Kr is unfractionated. These observations seem to imply that Xe has escaped from Earth. But to date no process has been identified that can cause Xe, which is heavier than Kr, to escape while Kr does not. Vigorous hydrodynamic hydrogen escape can produce mass fractionation in heavy gases. The required hydrogen flux is very high but within the possible range permitted by solar Extreme Ultraviolet radiation (EUV, which here means radiation at wavelengths short enough to be absorbed efficiently by hydrogen) heating when Earth was on the order of 100 Myrs old or younger. However this model cannot explain why Xe escapes but Kr does not. Recently, what appears to be ancient atmospheric xenon has been recovered from several very ancient (3-3.5 Ga) terrestrial hydrothermal barites and cherts. What is eye-catching about this ancient Xe is that it is less fractionated that Xe in modern air. In other words, it appears that a process was active on Earth some 3 to 3.5 billion years ago that caused xenon to fractionate. By this time the Sun was no longer the EUV source that it used to. If xenon was being fractionated by escape — currently the only viable hypothesis — it had to be in the less unfamiliar context of Earth’s Archean atmosphere and under rather modest levels of EUV forcing. This requires a new model. Here we address the circumstances in which Xe, but not Kr, could escape from Earth as an ion. In a hydrodynamically escaping hydrogen wind the hydrogen is partially photo-ionized. The key concepts are that ions are much more strongly coupled to the escaping flow than are neutrals (so that a relatively modest flow of H and H+ to space could carry Xe+ along with it), and that xenon alone among the noble gases is more easily ionized than hydrogen. This sort of escape is possible if not prevented by a planetary magnetic field. The best prospects for Earth are therefore escape along the polar field lines, although a very weak or absent magnetic field would likely work as well. As applied to the Archean Earth the discussion will be constrained by diffusion-limited hydrogen escape. The extended history of hydrogen escape implicit in Xe escape in the Archean is consistent with suggestions that hydrogen escape from the anoxic Archean atmosphere was considerable, because biogenic methane is expected to have been rather abundant. Hydrogen escape plausibly played the key role in creating oxidizing condition at the surface of the Earth and setting the stage for the creation of an O2 atmosphere.

photo-ionized gases↗
NASA NTRS2015

Swarm Flyby Gravimetry

This study describes a new technology for discerning the gravity fields and mass distribution of a solar system small body, without requiring dedicated orbiters or landers. Instead of a lander, a spacecraft releases a collection of small, simple probes during a flyby past an asteroid or comet. By tracking those probes from the host spacecraft, one can estimate the asteroid's gravity field and infer its underlying composition and structure. This approach offers a diverse measurement set,equivalent to planning and executing many independent and unique flyby encounters of a single spacecraft. This report assesses a feasible hardware implementation, derives the underlying models,and analyzes the performance of this concept via simulation.In terms of hardware, a small, low mass, low cost implementation is presented, which consists of a dispenser and probes. The dispenser contains roughly 12 probes in a tube and has a total size commensurate with a 6U P-Pod. The probes are housed in disc shaped sabots. When commanded,the dispenser ejects the top-most probe using a linear motor. The ejected probe separates from its sabots and unfolds using internal springs. There are two types of probes, each designed for a particular tracking modality. The reflective probe type, tracked by a telescope, unfolds to forma diffusely reflective sphere. The retroreflector probe type, tracked by a lidar, unfolds to form a corner-cube retroreflector assembly. Both types are designed to spherical so that their attitude doesn't affect the spacecraft's tracking performance.This analysis indicates that the point-mass term of small bodies larger than roughly 500 m in diameter can be observed from a host spacecraft that tracks locally deployed probes throughout a flyby to an uncertainty of better than 5%. The conditions by which this measurement is possible depends on the characteristics of the asteroid (size, type), the flyby velocity, and the type of tracking available (angles-only or angles+ranging). For most encounters, a few (1-3) well placed probes can be very effective, with marginal improvement for additional probes. Given realistic deployment errors, an encounter may require roughly 10-12 probes to ensure that 1-3 achieve their target. Long duration tracking of probes flying by large asteroids (>5 km diameter) can sometimes provide observability of the gravity field's first spherical harmonic, J( sub 2). In summary, this method offers a feasible, affordable approach to enabling or augmenting flyby science.

Atchison, Justin↗
NASA NTRS1996

Propagation and Breaking at High Altitudes of Gravity Waves Excited by Tropospheric Forcing

An anelastic approximation is used with a time-variable coordinate transformation to formulate a two-dimensional numerical model that describes the evolution of gravity waves. The model is solved using a semi-Lagrangian method with monotone (nonoscillatory) interpolation of all advected fields. The time-variable transformation is used to generate disturbances at the lower boundary that approximate the effect of a traveling line of thunderstorms (a squall line) or of flow over a broad topographic obstacle. The vertical propagation and breaking of the gravity wave field (under conditions typical of summer solstice) is illustrated for each of these cases. It is shown that the wave field at high altitudes is dominated by a single horizontal wavelength; which is not always related simply to the horizontal dimension of the source. The morphology of wave breaking depends on the horizontal wavelength; for sufficiently short waves, breaking involves roughly one half of the wavelength. In common with other studies, it is found that the breaking waves undergo "self-acceleration," such that the zonal-mean intrinsic frequency remains approximately constant in spite of large changes in the background wind. It is also shown that many of the features obtained in the calculations can be understood in terms of linear wave theory. In particular, linear theory provides insights into the wavelength of the waves that break at high altitudes, the onset and evolution of breaking. the horizontal extent of the breaking region and its position relative to the forcing, and the minimum and maximum altitudes where breaking occurs. Wave breaking ceases at the altitude where the background dissipation rate (which in our model is a proxy for molecular diffusion) becomes greater than the rate of dissipation due to wave breaking, This altitude, in effect, the model turbopause, is shown to depend on a relatively small number of parameters that characterize the waves and the background state.

Prusa, Joseph M.↗
NASA NTRS1999

Characteristics of Non-Premixed Turbulent Flames in Microgravity

The momentum of the fuel (and/or air) jet is important in classifying gas-jet diffusion flame behavior. Normal-gravity data on gas-jet flames show that the flame height (non-dimensionalized with respect to an effective diameter) can be correlated to a density weighted Froude number in the buoyancy-dominated limit. In the momentum-dominated limit this non-dimensional flame height asymptotes to a constant value. The momentum-dominated limit under normal gravity conditions is usually obtained for very high injection velocities which in turn results in high values of the injection Reynolds number. This results in a complicated flame structure because of the large number of turbulence scales involved. In order to gain better insight into the structure of these flames it would be useful to reduce the injection Reynolds number while still maintaining turbulent conditions. This can be done in microgravity where momentum-dominated turbulent flames are obtained at much smaller velocities than in normal gravity. In this paper, experimental results on the effects of nozzle diameter and fuel dilution on flame height are discussed. The experimental values are compared with predictions from a numerical procedure utilizing the standard k-epsilon turbulence model. Flame height scaling with nozzle size and dilution is established. Differences between model predictions and measurements are presented. In order to explain these differences, evolutions of turbulent spectra and Taylor microscale along the flame axis are considered.

Hegde, U.↗
NASA NTRS1989

Chemoviscosity modeling for thermosetting resin systems, 4

An experimental study on the changes of chemorheological properties has been conducted and analyzed on commercial Hercules 3501-6 resin system cured under several isothermal conditions between 375 and 435 K. For the cure temperatures equal to or greater than 385 K, the storage modulus curing curves, G prime (t), exhibited abrupt changes in slope which occurred at various times depending on the curing temperatures and were attributed to the onset of gelation reactions. The crossover points between G prime (t) and G double prime (t) curves were observed for curing temperatures equal to or greater than 400 K. The gelation and the crossover points obtained from the chemorheological measurements, therefore, defined two characteristic resin states during cure. Approximately the same value for the degree of cure was reached by the advancement of the reaction at each of these states. The temperature dependency of the viscosities for the characteristic resin states and the rate constants of increase in moduli at different stages of curing were analyzed. Various G prime (t) and G double prime (t) isothermal curing curves were also shown to be capable of being superimposed on one another by the principle of time-temperature superposition. The resultant shift factors a sub t(t) and a Eta(T) were shown to follow the Arrhenius type relationship. Values of the activation energy suggested that the reaction kinetics, instead of the diffusion mechanism, was the limiting step in the overall resin advancement for the cure at temperatures equal to or greater than 385 K.

Hou, T. H.↗
NASA NTRS2012

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗
NASA NTRS2015

Free-Energy Calculations. A Mathematical Perspective

Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.

free energy↗
NASA NTRS2021

The Response of the Amazon Ecosystem to the Photosynthetically Active Radiation Fields: Integrating Impacts of Biomass Burning Aerosol and Clouds in the NASA GEOS Earth System Model

The Amazon experiences fires every year, and the resulting biomass burning aerosols, together with cloud particles, influence the penetration of sunlight through the atmosphere, increasing the ratio of diffuse to direct photosynthetically active radiation (PAR) reaching the vegetation canopy and thereby potentially increasing ecosystem productivity. In this study, we use the NASA Goddard Earth Observing System (GEOS) model with coupled aerosol, cloud, radiation, and ecosystem modules to investigate the impact of Amazon biomass burning aerosols on ecosystem productivity, as well as the role of the Amazon’s clouds in tempering this impact. The study focuses on a seven-year period (2010-2016) during which the Amazon experienced a variety of dynamic environments (e.g., La Niña, normal years, and El Niño). The direct radiative impact of biomass burning aerosols on ecosystem productivity—called here the aerosol diffuse radiation fertilization effect —is found to increase Amazonian Gross Primary Production (GPP) by 2.6% via a 3.8% increase in diffuse PAR (DFPAR) despite a 5.4% decrease in direct PAR (DRPAR) on multiyear average during burning seasons. On a monthly basis, this increase in GPP can be as large as 9.9% (occurring in August 2010). Consequently, the net primary production (NPP) in Amazon is increased by 1.5%, or ~92 Tg C a-1– equivalent to ~37% of the average carbon lost due to Amazon fires over the seven years considered. Clouds, however, strongly regulate the effectiveness of the aerosol diffuse radiation fertilization effect. The efficiency of this fertilization effect is the highest in cloud-free conditions and linearly decreases with increasing cloud amount until the cloud fraction reaches ~0.8, at which point the aerosol-influenced light changes from being a stimulator to an inhibitor of plant growth. Nevertheless, interannual changes in the overall strength of the aerosol diffuse radiation fertilization effect are primarily controlled by the large interannual changes in biomass burning aerosols rather than by changes in cloudiness during the studied period.

biomass burning aerosols on ecosystem productivity↗
NASA NTRS2014

Behavior of Particle Depots in Molten Silicon During Float-Zone Growth in Strong Static Magnetic Fields

In the photovoltaics industry, the largest market share is represented by solar cells made from multicrystalline silicon, which is grown by directional solidification. During the growth process, the silicon melt is in contact with the silicon nitride coated crucible walls and the furnace atmosphere which contains carbon monoxide. The dissolution of the crucible coating, the carbon bearing gas, and the carbon already present in the feedstock, lead to the precipitation of silicon carbide, and silicon nitride, at later stages of the growth process. The precipitation of Si3N4 and SiC particles of up to several hundred micrometers in diameter leads to severe problems during the wire sawing process for wafering the ingots. Furthermore the growth of the silicon grains can be negatively influenced by the presence of particles, which act as nucleation sources and lead to a grit structure of small grains and are sources for dislocations. If doped with Nitrogen from the dissolved crucible coating, SiC is a semi conductive material, and can act as a shunt, short circuiting parts of the solar cell. For these reasons, the incorporation of such particles needs to be avoided. In this contribution we performed model experiments in which the transport of intentionally added SiC particles and their interaction with the solid-liquid interface during float zone growth of silicon in strong steady magnetic fields was investigated. SiC particles of 7μm and 60μm size are placed in single crystal silicon [100] and [111] rods of 8mm diameter. This is achieved by drilling a hole of 2mm diameter, filling in the particles and closing the hole by melting the surface of the rod until a film of silicon covers the hole. The samples are processed under a vacuum of 1x10(exp -5) mbar or better, to prevent gas inclusions. An oxide layer to suppress Marangoni convection is applied by wet oxidation. Experiments without and with static magnetic field are carried out to investigate the influence of melt convection on the distribution of particles and their incorporation into the crystal. The field strengths applied by a superconducting magnet are 1T, 3T, 4.5T, and 5T. The increase in field strength dampens the melt flow, and so this study provides comparative data to the crystal growth experiment to be carried out onboard the sounding rocket mission TEXUS 51, where purely diffusive growth condition will be achieved under microgravity conditions.

Jauss, T.↗
NASA NTRS2000

State Relationships of Laminar Permanently-Blue Opposed-Jet Hydrocarbon-Fueled Diffusion Flames

The structure and state relationships of laminar soot-free (permanently-blue) diffusion flames at various strain rates were studied experimentally using an opposed-jet configuration, motivated by the importance of soot-free hydrocarbon-fueled diffusion flames for many practical applications. Measurements of gas velocities, temperatures and compositions were carried out along the stagnation stream line. Flame conditions studied included propylene- and 1,3-butadiene-fueled opposed-jet diffusion flames having a stoichiometric mixture fractions of 0.7 and strain rates of 60-240/s at normal temperature and pressure. It was found that oxygen leakage to fuel-rich conditions and carbon monoxide leakage to fuel-lean conditions both increased as strain rates increased. Furthermore, increased strain rates caused increased fuel concentrations near the flame sheet decreased peak gas temperatures, and decreased concentrations of carbon dioxide and water vapor throughout the flames. State relationships for major gas species and gas temperatures for these flames were found to exist over broad ranges of strain rates. In addition, current measurements, as well as previous measurements and predictions of ethylene-fueled permanently-blue diffusion flames, all having a stoichiometric mixture fraction of 0.7, were combined to establish generalized state relationships for permanently-blue diffusion flames for this stoichiometric mixture fraction. The combined measurements and predictions support relatively universal generalized state relationships for N2, CO2, H2O and fuel over a broad range of strain rates and fuel types. State relationships for O2 in the fuel-rich region, and for CO in the fuel-lean region, however, are functions of strain rate and fuel type. State relationships for H2 and temperature exhibit less universality, mainly due to the increased experimental uncertainties for these variables. The existence of state relationships for soot-free hydrocarbon-fueled diffusion flames provides potential for significant computational simplifications for modeling purposes in many instances, allowing for effects of finite-rate chemistry while avoiding time-consuming computations of Arrhenius expressions.

Lin, K.-C.↗
NASA NTRS2000

State Relationships of Laminar Permanently-Blue Opposed-Jet Hydrocarbon-Fueled Diffusion Flames

The structure and state relationships of laminar soot-free (permanently-blue) diffusion flames at various strain rates were studied experimentally using an opposed-jet configuration, motivated by the importance of soot-free hydrocarbon-fueled diffusion flames for many practical applications. Measurements of gas velocities, temperatures and compositions were carried out along the stagnation stream line. Flame conditions studied included propylene- and 1,3-butadiene-fueled opposed-jet diffusion flames having a stoichiometric mixture fractions of 0.7 and strain rates of 60-240 s (exp -1) at normal temperature and pressure. It was found that oxygen leakage to fuel-rich conditions and carbon monoxide leakage to fuel-lean conditions both increased as strain rates increased. Furthermore, increased strain rates caused increased fuel concentrations near the flame sheet, decreased peak gas temperatures, and decreased concentrations of carbon dioxide and water vapor throughout the flames. State relationships for major gas species and gas temperatures for these flames were found to exist over broad ranges of strain rates. In addition, current measurements, as well as previous measurements and predictions of ethylene-fueled permanently-blue diffusion flames, all having a stoichiometric mixture fraction of 0.7, were combined to establish generalized state relationships for permanently-blue diffusion flames for this stoichiometric mixture fraction. The combined measurements and predictions support relatively universal generalized state relationships for N2, CO2, H2O and fuel over a broad range of strain rates and fuel types. State relationships for O2 in the fuel-rich region, and for CO in the fuel-lean region, however, are functions of strain rate and fuel type. State relationships for H2 and temperature exhibit less universality, mainly due to the increased experimental uncertainties for these variables. The existence of state relationships for soot-free hydrocarbon-fueled diffusion flames provides potential for significant computational simplifications for modeling purposes in many instances, allowing for effects of finite-rate chemistry while avoiding time-consuming computations of Arrhenius expressions.

Lin, K.-C.↗
NASA NTRS1995

Velocity Fields of Axisymmetric Hydrogen-Air Counterflow Diffusion Flames from LDV, PIV, and Numerical Computation

Laminar fuel-air counterflow diffusion flames (CFDFs) were studied using axisymmetric convergent-nozzle and straight-tube opposed jet burners (OJBs). The subject diagnostics were used to probe a systematic set of H2/N2-air CFDFs over wide ranges of fuel input (22 to 100% Ha), and input axial strain rate (130 to 1700 Us) just upstream of the airside edge, for both plug-flow and parabolic input velocity profiles. Laser Doppler Velocimetry (LDV) was applied along the centerline of seeded air flows from a convergent nozzle OJB (7.2 mm i.d.), and Particle Imaging Velocimetry (PIV) was applied on the entire airside of both nozzle and tube OJBs (7 and 5 mm i.d.) to characterize global velocity structure. Data are compared to numerical results from a one-dimensional (1-D) CFDF code based on a stream function solution for a potential flow input boundary condition. Axial strain rate inputs at the airside edge of nozzle-OJB flows, using LDV and PIV, were consistent with 1-D impingement theory, and supported earlier diagnostic studies. The LDV results also characterized a heat-release hump. Radial strain rates in the flame substantially exceeded 1-D numerical predictions. Whereas the 1-D model closely predicted the max I min axial velocity ratio in the hot layer, it overpredicted its thickness. The results also support previously measured effects of plug-flow and parabolic input strain rates on CFDF extinction limits. Finally, the submillimeter-scale LDV and PIV diagnostics were tested under severe conditions, which reinforced their use with subcentimeter OJB tools to assess effects of aerodynamic strain, and fueVair composition, on laminar CFDF properties, including extinction.

Pellett, Gerald L.↗
NASA NTRS2014

An Experimental and Computational Study on Soot Formation in a Coflow Jet Flame Under Microgravity and Normal Gravity

Upon the completion of the Structure and Liftoff in Combustion Experiment (SLICE) in March 2012, a comprehensive and unique set of microgravity coflow diffusion flame data was obtained. This data covers a range of conditions from weak flames near extinction to strong, highly sooting flames, and enabled the study of gravitational effects on phenomena such as liftoff, blowout and soot formation. The microgravity experiment was carried out in the Microgravity Science Glovebox (MSG) on board the International Space Station (ISS), while the normal gravity experiment was performed at Yale utilizing a copy of the flight hardware. Computational simulations of microgravity and normal gravity flames were also carried out to facilitate understanding of the experimental observations. This paper focuses on the different sooting behaviors of CH4 coflow jet flames in microgravity and normal gravity. The unique set of data serves as an excellent test case for developing more accurate computational models.Experimentally, the flame shape and size, lift-off height, and soot temperature were determined from line-of-sight flame emission images taken with a color digital camera. Soot volume fraction was determined by performing an absolute light calibration using the incandescence from a flame-heated thermocouple. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the chemically reacting flow, and the soot evolution was modeled by the sectional aerosol equations. The governing equations and boundary conditions were discretized on an axisymmetric computational domain by finite differences, and the resulting system of fully coupled, highly nonlinear equations was solved by a damped, modified Newtons method. The microgravity sooting flames were found to have lower soot temperatures and higher volume fraction than their normal gravity counterparts. The soot distribution tends to shift from the centerline of the flame to the wings from normal gravity to microgravity.

microgravity↗
NASA NTRS2000

An Overview of the NCC Spray/Monte-Carlo-PDF Computations

This paper advances the state-of-the-art in spray computations with some of our recent contributions involving scalar Monte Carlo PDF (Probability Density Function), unstructured grids and parallel computing. It provides a complete overview of the scalar Monte Carlo PDF and Lagrangian spray computer codes developed for application with unstructured grids and parallel computing. Detailed comparisons for the case of a reacting non-swirling spray clearly highlight the important role that chemistry/turbulence interactions play in the modeling of reacting sprays. The results from the PDF and non-PDF methods were found to be markedly different and the PDF solution is closer to the reported experimental data. The PDF computations predict that some of the combustion occurs in a predominantly premixed-flame environment and the rest in a predominantly diffusion-flame environment. However, the non-PDF solution predicts wrongly for the combustion to occur in a vaporization-controlled regime. Near the premixed flame, the Monte Carlo particle temperature distribution shows two distinct peaks: one centered around the flame temperature and the other around the surrounding-gas temperature. Near the diffusion flame, the Monte Carlo particle temperature distribution shows a single peak. In both cases, the computed PDF's shape and strength are found to vary substantially depending upon the proximity to the flame surface. The results bring to the fore some of the deficiencies associated with the use of assumed-shape PDF methods in spray computations. Finally, we end the paper by demonstrating the computational viability of the present solution procedure for its use in 3D combustor calculations by summarizing the results of a 3D test case with periodic boundary conditions. For the 3D case, the parallel performance of all the three solvers (CFD, PDF, and spray) has been found to be good when the computations were performed on a 24-processor SGI Origin work-station.

Raju, M. S.↗
NASA NTRS1992

Second order closure modeling of turbulent buoyant wall plumes

Non-intrusive measurements of scalar and momentum transport in turbulent wall plumes, using a combined technique of laser Doppler anemometry and laser-induced fluorescence, has shown some interesting features not present in the free jet or plumes. First, buoyancy-generation of turbulence is shown to be important throughout the flow field. Combined with low-Reynolds-number turbulence and near-wall effect, this may raise the anisotropic turbulence structure beyond the prediction of eddy-viscosity models. Second, the transverse scalar fluxes do not correspond only to the mean scalar gradients, as would be expected from gradient-diffusion modeling. Third, higher-order velocity-scalar correlations which describe turbulent transport phenomena could not be predicted using simple turbulence models. A second-order closure simulation of turbulent adiabatic wall plumes, taking into account the recent progress in scalar transport, near-wall effect and buoyancy, is reported in the current study to compare with the non-intrusive measurements. In spite of the small velocity scale of the wall plumes, the results showed that low-Reynolds-number correction is not critically important to predict the adiabatic cases tested and cannot be applied beyond the maximum velocity location. The mean and turbulent velocity profiles are very closely predicted by the second-order closure models. but the scalar field is less satisfactory, with the scalar fluctuation level underpredicted. Strong intermittency of the low-Reynolds-number flow field is suspected of these discrepancies. The trends in second- and third-order velocity-scalar correlations, which describe turbulent transport phenomena, are also predicted in general, with the cross-streamwise correlations better than the streamwise one. Buoyancy terms modeling the pressure-correlation are shown to improve the prediction slightly. The effects of equilibrium time-scale ratio and boundary condition are also discussed.

Zhu, Gang↗
NASA NTRS1993

A model of concurrent flow flame spread over a thin solid fuel

A numerical model is developed to examine laminar flame spread and extinction over a thin solid fuel in lowspeed concurrent flows. The model provides a more precise fluid-mechanical description of the flame by incorporating an elliptic treatment of the upstream flame stabilization zone near the fuel burnout point. Parabolic equations are used to treat the downstream flame, which has a higher flow Reynolds number. The parabolic and elliptic regions are coupled smoothly by an appropriate matching of boundary conditions. The solid phase consists of an energy equation with surface radiative loss and a surface pyrolysis relation. Steady spread with constant flame and pyrolysis lengths is found possible for thin fuels and this facilitates the adoption of a moving coordinate system attached to the flame with the flame spread rate being an eigen value. Calculations are performed in purely forced flow in a range of velocities which are lower than those induced in a normal gravity buoyant environment. Both quenching and blowoff extinction are observed. The results show that as flow velocity or oxygen percentage is reduced, the flame spread rate, the pyrolysis length, and the flame length all decrease, as expected. The flame standoff distance from the solid and the reaction zone thickness, however, first increase with decreasing flow velocity, but eventually decrease very near the quenching extinction limit. The short, diffuse flames observed at low flow velocities and oxygen levels are consistent with available experimental data. The maximum flame temperature decreases slowly at first as flow velocity is reduced, then falls more steeply close to the quenching extinction limit. Low velocity quenching occurs as a result of heat loss. At low velocities, surface radiative loss becomes a significant fraction of the total combustion heat release. In addition, the shorter flame length causes an increase in the fraction of conduction downstream compared to conduction to the fuel. These heat losses lead to lower flame temperatures, and ultimately, extinction. This extinction mechanism differs from that of blowoff, where the flame is unable to be stabilized due to the high flow velocity.

Ferkul, Paul V.↗