Search NASASearch

SEARCH · Search NASA

Results for “generative diffusion model”

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 109 records · Page 6

Badhwar-O'Neill 2011 Galactic Cosmic Ray Model Update and Future Improvements

The Badhwar-O'Neill Galactic Cosmic Ray (GCR) Model based on actual GCR measurements is used by deep space mission planners for the certification of microelectronic systems and the analysis of radiation health risks to astronauts in space missions. The BO GCR Model provides GCR flux in deep space (outside the earth's magnetosphere) for any given time from 1645 to present. The energy spectrum from 50 MeV/n - 20 GeV/n is provided for ions from hydrogen to uranium. This work describes the most recent version of the BO GCR model (BO'11). BO'11 determined the GCR flux at a given time applying an emperical time delay function to past sunspot activity. We describe the GCR measurement data used in the BO'11 update - modern data from BESS, PAMELA, CAPRICE, and ACE emphasized more than the older balloon data used for the previous BO model (BO'10). We look at the GCR flux for the last 24 solar minima and show how much greater the flux was for the cycle 24 minimum in 2010. The BO'11 Model uses the traditional, steady-state Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration. It assumes a radially symmetrical diffusion coefficient derived from magnetic disturbances caused by sunspots carried outward by a constant solar wind. A more complex differential equation is now being tested to account for particle transport in the heliosphere in the next generation BO model. This new model is time-dependent (no longer a steady state model). In the new model, the dynamics and anti-symmetrical features of the actual heliosphere are accounted for so emperical time delay functions will no longer be required. The new model will be capable of simulating the more subtle features of modulation - such as the Sun's polarity and modulation dependence on gradient and curvature drift. This improvement is expected to significantly improve the fidelity of the BO GCR model. Preliminary results of its performance will be presented.

O'Neill, Pat M.

Badhwar-O'Neill 2011 Galactic Cosmic Ray Model Update and Future Improvements

The Badhwar-O'Neill Galactic Cosmic Ray (GCR) Model based on actual GR measurements is used by deep space mission planners for the certification of micro-electronic systems and the analysis of radiation health risks to astronauts in space missions. The BO GCR Model provides GCR flux in deep space (outside the earth's magnetosphere) for any given time from 1645 to present. The energy spectrum from 50 MeV/n-20 GeV/n is provided for ions from hydrogen to uranium. This work describes the most recent version of the BO GCR model (BO'11). BO'11 determines the GCR flux at a given time applying an empirical time delay function to past sunspot activity. We describe the GCR measurement data used in the BO'11 update - modern data from BESS, PAMELA, CAPRICE, and ACE emphasized for than the older balloon data used for the previous BO model (BO'10). We look at the GCR flux for the last 24 solar minima and show how much greater the flux was for the cycle 24 minimum in 2010. The BO'11 Model uses the traditional, steady-state Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration. It assumes a radially symmetrical diffusion coefficient derived from magnetic disturbances caused by sunspots carried onward by a constant solar wind. A more complex differential equation is now being tested to account for particle transport in the heliosphere in the next generation BO model. This new model is time-dependent (no longer a steady state model). In the new model, the dynamics and anti-symmetrical features of the actual heliosphere are accounted for so empirical time delay functions will no longer be required. The new model will be capable of simulating the more subtle features of modulation - such as the Sun's polarity and modulation dependence on the gradient and curvature drift. This improvement is expected to significantly improve the fidelity of the BO GCR model. Preliminary results of its performance will be presented.

O'Neill, Pat M.

Precipitation Model Validation in 3rd Generation Aeroturbine Disc Alloys

In support of application of the DARPA-AIM methodology to the accelerated hybrid thermal process optimization of 3rd generation aeroturbine disc alloys with quantified uncertainty, equilibrium and diffusion couple experiments have identified available fundamental thermodynamic and mobility databases of sufficient accuracy. Using coherent interfacial energies quantified by Single-Sensor DTA nucleation undercooling measurements, PrecipiCalc(TM) simulations of nonisothermal precipitation in both supersolvus and subsolvus treated samples show good agreement with measured gamma particle sizes and compositions. Observed longterm isothermal coarsening behavior defines requirements for further refinement of elastic misfit energy and treatment of the parallel evolution of incoherent precipitation at grain boundaries.

Olson, G. B.

Use of experimental separation limits in the theoretical design of V/STOL inlets

Experimental data from several model inlets have been used to generate two parameters which are related to the limit of operation for inlet flow separation. One parameter, called the diffusion ratio, is the ratio of the peak velocity on the inlet surface to the velocity at the diffuser exit and is related to the boundary-layer separation at low throat Mach numbers. The other parameter, the peak Mach number on the inlet surface, is related to the separation at high throat Mach numbers. These parameters are easily calculated from potential flow solutions and thus can be used as a design tool in screening proposed inlet geometries. Any of the geometric design variables can be analyzed by this technique; but, this paper is restricted to the consideration of the internal lip contraction ratio. An illustrative example of an application to an inlet design study for a tilt nacelle VTOL airplane is presented. The study will show what value of contraction ratio is required to meet the operating requirements yet allow the inlet to remain free of separation as indicated by the two separation parameters.

Boles, M. A.

Ion-beam-induced atomic mixing

Calculations based on the diffusion model are presented of atomic mixing by ion bombardment. This mixing is assumed to have its basis, as does sputtering, in the collision cascades generated by the primary beam. Sharp interfaces within a target are seen to be smoothed by ion bombardment. Mixing may place fundamental limits on the resolution of ion microprobes.

Haff, P. K.

Low speed test of the aft inlet designed for a tandem fan V/STOL nacelle

An approximately .25 scale model of a Tandem Fan nacelle designed for a Type A V/STOL aircraft configuration was tested in a 10-by-10 foot wind tunnel. A 12 inch, tip driven, turbofan simulator was used to provide the suction source for the aft fan inlet. The front fan inlet was faired over for this test entry. Model variables consisted of a long aft inlet cowl, a short aft inlet cowl, a shaft simulator, blow-in door passages and diffuser vortex generators. Inlet pressure recovery, distortion, inlet angle of attack separation limits were evaluated at tunnel velocities from 0 to 240 knots, angles of attack from -10 to 40 degrees and inlet flow rates representative of throat Mach numbers of 0.1 to 0.6. High inlet performance and stable operation was verified at all design forward speed and angle of attack conditions. The short aft inlet configuration provided exceptionally high pressure recovery except at the highest combination of angle of attack and forward speed. The flow quality at the fan face was somewhat degraded by the addition of blow-in door passages to the long aft inlet configuration due to the pressure disturbances generated by the flow entering the diffuser through the auxiliary air passages.

Rhoades, W. W.

Additional testing of the inlets designed for a tandem fan V/STOL nacelle

The wind tunnel testing of a scale model of a tandem fan nacelle designed for a type (subsonic cruise) V/STOL aircraft configuration is discussed. The performance for the isolated front inlet and for the combined front and aft inlets is reported. Model variables include front and aft inlets with aft inlet variations of short and long aft inlet cowls, with a shaft simulator and diffuser vortex generators, cowl lip fillets, and nacelle strakes. Inlet pressure recovery, distortion, and inlet angle-to-attack separation limits were evaluated at tunnel velocity from 0 to 240 knots, angles-of-attack from -10 to +40 degrees and inlet flow rates corresponding to throat Mach number from 0.0 to 0.6. Combined nacelle pitch and yaw runs up to 30 deg. were also made.

Ybarra, A. H.

Fractal dynamics of bioconvective patterns

Biologically generated cellular patterns, sometimes called bioconvective patterns, are found to cluster into aggregates which follow fractal growth dynamics akin to diffusion-limited aggregation (DLA) models. The pattern formed is self-similar with fractal dimension of 1.66 +/-0.038. Bioconvective DLA branching results from thermal roughening which shifts the balance between ordering viscous forces and disordering cell motility and random diffusion. The phase diagram for pattern morphology includes DLA, boundary spokes, random clusters, and reverse clusters.

Noever, David A.

Inducer analysis/pump model development

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design information in a productive manner. The main goal of this study was to develop a robust and effective computational fluid dynamics (CFD) pump model for general turbopump design and analysis applications. A finite difference Navier-Stokes flow solver, FDNS, which includes an extended k-epsilon turbulence model and appropriate moving zonal interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. In the present study, three key components of the turbopump, the inducer, impeller, and diffuser, were investigated by the proposed pump model, and the numerical results were benchmarked by the experimental data provided by Rocketdyne. For the numerical calculation of inducer flows with tip clearance, the turbulence model and grid spacing are very important. Meanwhile, the development of the cross-stream secondary flow, generated by curved blade passage and the flow through tip leakage, has a strong effect on the inducer flow. Hence, the prediction of the inducer performance critically depends on whether the numerical scheme of the pump model can simulate the secondary flow pattern accurately or not. The impeller and diffuser, however, are dominated by pressure-driven flows such that the effects of turbulence model and grid spacing (except near leading and trailing edges of blades) are less sensitive. The present CFD pump model has been proved to be an efficient and robust analytical tool for pump design due to its very compact numerical structure (requiring small memory), fast turnaround computing time, and versatility for different geometries.

Cheng, Gary C.

Stream interfaces and energetic ions closer than expected: Analyses of Pioneers 10 and 11 observations

An empirical study of corotating interaction regions (CIRs) observed between 3.9 AU and 5.9 AU on Pioneers 10 and 11 shows that the main corotation energetic ion population (CEIP), which is associated with the trailing reverse shock, terminates within the CIR at a definite, structural boundary, which we show here is the stream interface. This new result has significant implications for solar wind and energetic particle modeling. In particular it implies either that the reverse shock forms closer to the stream interface than models suggest or that the theories that treat the generation and transport of these energetic ions, such as preshock Fermi acceleration and cross-field diffusion must be combined or extended. We test these scenarios by comparing the CEIP intensity profiles on the two sides of the stream interface. We find that while each automatically accounts for one or two aspects of the results none of them alone can account for all of our empirical results.

Intriligator, Devrie S.

A source-sink model of the generation of plate tectonics from non-Newtonian mantle flow

A model of mantle convection which generates plate tectonics requires strain rate- or stress-dependent rheology in order to produce strong platelike flows with weak margins as well as strike-slip deformation and plate spin (i.e., toroidal motion). Here, we employ a simple model of source-sink driven surface flow to determine the form of such a rheology that is appropriate for Earth's present-day plate motions. In this model, lithospheric motion is treated as shallow layer flow driven by sources and sinks which correspond to spreading centers and subduction zones, respectively. Two plate motion models are used to derive the source sink field. As originally implied in the simpler Cartesian version of this model, the classical power law rheologies do not generate platelike flows as well as the hypothetical Whitehead-Gans stick-slip rheology (which incorporates a simple self-lubrication mechanism). None of the fluid rheologies examined, however, produce more than approximately 60% of the original maximum shear. For either plate model, the viscosity fields produced by the power law rheologies are diffuse, and the viscosity lows over strike-slip shear zones or pseudo-margins are not as small as over the prescribed convergent-divergent margins. In contrast, the stick-slip rheology generates very platelike viscosity fields, with sharp gradients at the plate boundaries, and margins with almost uniformly low viscosity. Power law rheologies with high viscosity contrasts, however, lead to almost equally favorable comparisons, though these also yield the least platelike viscosity fields. This implies that the magnitude of toroidal flow and platelike strength distributions are not necessarily related and thus may present independent constraints on the determination of a self-consistent plate-mantle rheology.

Bercovici, David

Simulation of the thermospheric tides by use of the NCAR thermospheric general circulation model

Numerical calculations of the thermospheric tidal winds and temperatures at equinox are presented. The calculations were made using the National Center for Atmospheric Research Thermospheric General Circulation Model (TGCM), which includes the effects of viscosity, conductivity, diffusion, ion drag, winds, and temperature gradients. The thermospheric diurnal and semidiurnal tides are excited in situ by solar heating and by ion-neutral momentum coupling. The semidiurnal tidal calculations also include the effects of upward propagating waves generated by heating in the lower atmosphere. This semidiurnal propagating component is modeled by use of the classical tidal perturbations as lower boundary conditions. The model is tuned by adjusting the propagating tidal forcing term until calculated semidiurnal wind and temperature fields best approximate incoherent scatter observations. The tidal TGCM results are consistent with previous theoretical work and successfully reproduce high altitude temperature and meridional velocity data, but they give significantly lower magnitudes for velocities and temperatures near 160 km than are seen by observations.

Fensen, C. G.

Gravitational effects on chemically reacting laminar boundary layer flows over a horizontal flat plate

A theoretical study of the chemically reacting laminar boundary layer flow over a horizontal flat plate with gravitationally induced buoyant force is presented. A diffusion flame sheet model was used to describe the combustion process. The effects of gravity on the purely force convection flow can be characterized by a dimensionless coordinate quantity, which is involved in the generation of the governing equations. A numerical solution of the zero and first order governing equations subject to the appropriate physical boundary conditions was obtained. It is shown that the cross stream buoyancy induced body force acts effectively to produce a streamwise pressure gradient in the fluid adjacent to the plate surface. It is concluded that buoyancy plays an important role in boundary layer diffusion flames.

Lavid, M.

Investigation of turbulent transport in an axisymmetric sudden expansion

Simultaneous two-component laser velocimeter measurements were made in the incompressible turbulent flowfield following an axisymmetric sudden expansion. Mean velocities, Reynolds stresses, and triple products were measured and are presented at axial positions ranging from x/H = 0.2-14. A balance of the turbulent kinetic energy in the flow was performed. The production, convection, and diffusion of turbulent kinetic energy were computed directly from the experimental data using central differencing. A specially designed correction lens was employed to correct for optical aberrations introduced by the circular tube. This lens system allowed the accurate simultaneous measurement of axial and radial velocities in the test section. The experimental measurements were compared to predictions generated by a code that employed the k-epsilon turbulence model. Agreement was good for mean axial velocities, turbulent kinetic energy, and turbulent shear stresses. However, the modeled turbulent normal stresses where in poor agreement with the measured values. The modeled diffusion of turbulent kinetic energy was underpredicted in the region between the shear layer and the centerline of the flow giving lower values of turbulent kinetic energy downstream of the potential core than measured.

Gould, Richard D.

Dynamics and structure of turbulent premixed flames

In earlier work (Mantel & Bilger, 1994) the structure of the turbulent premixed flame was investigated using statistics based on conditional averaging with the reaction progress variable as the conditioning variable. The DNS data base of Trouve and Poinsot (1994) was used in this investigation. Attention was focused on the conditional dissipation and conditional axial velocity in the flame with a view to modeling these quantities for use in the conditional moment closure (CMC) approach to analysis of kinetics in premixed flames (Bilger, 1993). Two remarkable findings were made: there was almost no acceleration of the axial velocity in the flame front itself; and the conditional scalar dissipation remained as high, or higher, than that found in laminar premixed flames. The first finding was surprising since in laminar flames all the fluid acceleration occurs through the flame front, and this could be expected also for turbulent premixed flames at the flamelet limit. The finding gave hope of inventing a new approach to the dynamics of turbulent premixed flames through use of rapid distortion theory or an unsteady Bernoulli equation. This could lead to a new second order closure for turbulent premixed flames. The second finding was contrary to our measurements with laser diagnostics in lean hydrocarbon flames where it is found that conditional scalar dissipation drops dramatically below that for laminar flamelets when the turbulence intensity becomes high. Such behavior was not explainable with a one-step kinetic model, even at non-unity Lewis number. It could be due to depletion of H2 from the reaction zone by preferential diffusion. The capacity of the flame to generate radicals is critically dependent on the levels of H2 present (Bilger, et al., 1991). It seemed that a DNS computation with a multistep reduced mechanism would be worthwhile if a way could be found to make this feasible. Truly innovative approaches to complex problems often come only when there is the opportunity to work close at hand with the (in this case numerical) experimental data. Not only can one spot patterns and relationships in the data which could be important, but one can also get to know the limitations of the technique being used, so that when the next experiment is being designed it will address resolvable questions. A three-year grant from the Australian Research Council has enabled us to develop a small capability at the University of Sydney to work on DNS of turbulent reacting flow, and to analyze data bases generated at CTR. Collaboration between the University of Sydney and CTR is essential to this project and finding a workable modus operandum for this collaboration, given the constraints involved, has been a major objective of the past year's effort. The overall objectives of the project are: (1) to obtain a quantitative understanding of the dynamics of turbulent premixed flames at high turbulence levels with a view to developing improved second order closure models; and (2) to carry out new DNS experiments on turbulent premixed flames using a carefully chosen multistep reduced mechanism for the chemical kinetics, with a view to elucidating the laser diagnostic findings that are contrary to the findings for DNS using one-step kinetics. In this first year the objectives have been to make the existing CTR data base more accessible to coworkers at the University of Sydney, to make progress on understanding the dynamics of the flame in this existing CTR data base, and to carefully construct a suitable multistep reduced mechanism for use in a new set of DNS experiments on turbulent premixed flames.

Bilger, R. W.

Shock Wave Dynamics in Weakly Ionized Plasmas

An investigation of the dynamics of shock waves in weakly ionized argon plasmas has been performed using a pressure ruptured shock tube. The velocity of the shock is observed to increase when the shock traverses the plasma. The observed increases cannot be accounted for by thermal effects alone. Possible mechanisms that could explain the anomalous behavior include a vibrational/translational relaxation in the nonequilibrium plasma, electron diffusion across the shock front resulting from high electron mobility, and the propagation of ion-acoustic waves generated at the shock front. Using a turbulence model based on reduced kinetic theory, analysis of the observed results suggest a role for turbulence in anomalous shock dynamics in weakly ionized media and plasma-induced hypersonic drag reduction.

Johnson, Joseph A., III

The role of gravity wave generated advection and diffusion in transport of tracers in the mesosphere

The roles of advection and turbulent diffusion in controlling tracer transport in the meridional plane of the mesosphere are compared using a version of the global primitive equation model for the middle atmosphere developed by Holton and Wehrbein (1980). The results of model sensitivity tests demonstrated that transport in the mesosphere for long-lived species with large-scale heights, such as H2O, CO, and NO(x), is primarily an advective process; on a global scale, small-scale turbulent diffusion produced by gravity wave breaking plays a secondary role. It is concluded that, for such tracers, the eddy diffusion coefficient used in one-dimensional chemical models of the mesosphere should be regarded as a parameterization of the transport by the meridional circulation, not as a parameterization of mixing by small-scale diffusion.

Holton, James R.

Profile Evolution Simulation in Etching Systems Using Level Set Methods

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

Hwang, Helen H.