Results for “Numerical diffusion” — Page 39

NASA NTRS1993

Alfven wave transport effects in the time evolution of parallel cosmic-ray-modified shocks

This paper presents a numerical study of the time evolution of plane, cosmic-ray modified shocks with magnetic field parallel to the shock normal, based on the diffusive shock acceleration formalism and including the effects from the finite propagation speed and energy of Alfven waves responsible for controlling the transport of the cosmic rays. The simulations discussed are based on a three-fluid model for the dynamics, but a more complete formalism is laid out for future work. The results of the simulations confirm earlier steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and the gas sound speed are comparable, i.e., when the plasma and magnetic pressures are similar. It is also clear, however, that the impact of Alfven transport effects, which tend to slow shock evolution and reduce the time asymptotic cosmic-ray pressure in the shock, is strongly dependent upon uncertain details in the transport models. Both cosmic-ray advection tied to streaming Alfven waves and dissipation of wave energy are important to include in the models. Further, Alfven transport properties on both sides of the shock are also influential.

Jones, T. W.↗

NASA NTRS1995

Numerical solution of fluid flow and heat tranfer problems with surface radiation

This paper presents a numerical scheme, based on the finite element method, to solve strongly coupled fluid flow and heat transfer problems. The surface radiation effect for gray, diffuse and isothermal surfaces is considered. A procedure for obtaining the view factors between the radiating surfaces is discussed. The overall solution strategy is verified by comparing the available results with those obtained using this approach. An analysis of a thermosyphon is undertaken and the effect of considering the surface radiation is clearly explained.

Ahuja, S.↗

NASA NTRS1997

The First National Student Conference: NASA University Research Centers at Minority Institutions

The conference includes contributions from 13 minority universities with NASA University Research Centers. Topics discussed include: leadership, survival strategies, life support systems, food systems, simulated hypergravity, chromium diffusion doping, radiation effects on dc-dc converters, metal oxide glasses, crystal growth of Bil3, science and communication on wheels, semiconductor thin films, numerical solution of random algebraic equations, fuzzy logic control, spatial resolution of satellite images, programming language development, nitric oxide in the thermosphere and mesosphere, high performance polyimides, crossover control in genetic algorithms, hyperthermal ion scattering, etc.

Daso, Endwell O.↗

NASA NTRS2004

An Investigation of Wave Propagations in Discontinuous Galerkin Method

Analysis of the discontinuous Galerkin method has been carried out for one- and two-dimensional system of hyperbolic equations. Analytical, as well as numerical, properties of wave propagation in a DGM scheme are derived and verified with direct numerical simulations. In addition to a systematic examination of the dissipation and dispersion errors, behaviours of a DG scheme at an interface of two different grid topologies are also studied. Under the same framework, a quantitative discrete analysis of various artificial boundary conditions is also conducted. Progress has been made in numerical boundary condition treatment that is closely related to the application of DGM in aeroacoustics problems. Finally, Fourier analysis of DGM for the Convective diffusion equation has also be studied in connection with the application of DG schemes for the Navier-Stokes equations. This research has resulted in five(5) publications, plus one additional manuscript in preparation, four(4) conference presentations, and three(3) departmental seminars, as summarized in part II. Abstracts of papers are given in part 111 of this report.

Hu, Fang Q.↗

NASA NTRS2018

The Two-Fluid Dynamics and Energetics of the Asymmetric Magnetic Reconnection in Laboratory and Space Plasmas

Magnetic reconnection is a fundamental process in magnetized plasma where magnetic energy is converted to plasma energy. Despite huge differences in the physical size of the reconnection layer, remarkably similar characteristics are observed in both laboratory and magnetosphere plasmas. Here we present the comparative study of the dynamics and physical mechanisms governing the energy conversion in the laboratory and space plasma in the context of two-fluid physics, aided by numerical simulations. In strongly asymmetric reconnection layers with negligible guide field, the energy deposition to electrons is found to primarily occur in the electron diffusion region where electrons are demagnetized and diffuse. A large potential well is observed within the reconnection plane and ions are accelerated by the electric field toward the exhaust region. The present comparative study identifies the robust two-fluid mechanism operating in systems over six orders of magnitude in spatial scales and over a wide range of collisionality.

Yamada, M.↗

NASA NTRS2008

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗

NASA NTRS2010

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗

NASA NTRS1981

Numerical simulation of reentry flow around the Shuttle Orbiter including real gas effects

The flight environment of a reentry vehicle is predicted from the numerical solution of fluid dynamics equations for the region between the blunt nose and the swept wings. The inviscid portion of the shock layer is modeled by the Euler equations, but the laminar viscous flow adjacent to the wall is modeled by the approximate parabolic Navier-Stokes equations. The approximations made to the axial gradients of pressure and diffusive fluxes enable the coupled inviscid and viscous equations to be solved efficiently along the body axis. The equilibrium air aftbody code contains significant improvements over its predecessor which only considers a perfect gas model and noncircular configurations. The inclusion of numerical damping either explicitly or implicitly has extended its capabilities for predicting flow field around a winged configuration at higher Machs and greater angles of attack. The results are obtained on the cylindrical coordinates and satisfactory for the Shuttle Orbiter at a free-stream Mach number of 22 and an angle of attack of 40 deg. Also discussed are the inviscid formulation and its application for wind-tunnel conditions.

Li, C. P.↗

NASA NTRS1985

A surface temperature and moisture parameterization for use in mesoscale numerical models

A modified multi-level soil moisture and surface temperature model is presented for use as in defining lower boundary conditions in mesoscale weather models. Account is taken of the hydraulic and thermal diffusion properties of the soil, their variations with soil type, and the mixing ratio at the surface. Techniques are defined for integrating the surface input into the multi-level scheme. Sample simulation runs were performed with the modified model and the original model defined by Pielke, et al. (1977, 1981). The models were applied to regional weather forecasting over soils composed of sand and clay loam. The new form of the model avoided iterations necessary in the earlier version of the model and achieved convergence at reasonable profiles for surface temperature and moisture in regions where the earlier version of the model failed.

Tremback, C. J.↗

NASA NTRS1999

A Theory of Oscillating Edge Flames

It has been known for some years that when a near-limit flame spreads over a liquid pool of fuel, the edge of the flame can oscillate relative to a frame moving with the mean speed. Each period of oscillation is characterized by long intervals of modest motion during which the edge gases radiate like those of a diffusion flame, punctuated by bursts of rapid advance during which the edge gases radiate like those in a deflagration. Substantial resources have been brought to bear on this issue within the microgravity program, both experimental and numerical. It is also known that when a near-asphyxiated candle-flame burns at zero gravity, the edge of the (hemispherical) flame can oscillate violently prior to extinction. Thus a web-surfer, turning to the NASA web-site at http://microgravity.msfc.nasa.gov, and following the trail combustion science/experiments/experimental results/candle flame, will find photographs and a description of candle burning experiments carried out on board both the Space-shuttle and the Russian space station Mir. A brief report can also be found in the proceedings of the Fourth Workshop. And recently, in a third microgravity program, the leading edge of the flame supported by injection of ethane through the porous surface of a plate over which air is blown has been found to oscillate when conditions are close to blow-off. A number of important points can be made with respect to these observations: It is the edge itself which oscillates, advancing and retreating, not the diffusion flame that trails behind the edge; oscillations only occur under near limit conditions; in each case the Lewis number of the fuel is significantly larger than 1; and because of the edge curvature, the heat losses from the reacting edge structure are larger than those from the trailing diffusion flame. We propose a general theory for these oscillations, invoking Occam's 'Law of Parsimony' in an expanded form, to wit: The same mechanism is responsible for the oscillations in all three experiments; and no new mechanism is invoked (Occam's original 'Razor'). Such a strategy eliminates Marangoni effects as the source, for these are absent in the second and third experiments. And it eliminates arguments that point to numerically predicted gas eddies as the source, a new mechanism, unelucidated. Indeed, we hypothesize that the essential driving mechanism for the instability is a combination of large Lewis number and heat losses from the reacting structure near the flame edge. Instabilities driven by these mechanisms are commonplace in 1D configurations. Chemical reactor theory, for example, leads to system responses which mimic the response of the candle flame - steady flame, oscillations, extinction. In a combustion context, oscillating instabilities were first reported for diffusion flames in a theoretical study by Kirkby and Schmitz, and here also the instabilities are associated with near-extinction conditions, large Lewis numbers, and heat losses. And deflagrations will oscillate if the Lewis number is large enough, oscillations that are exacerbated when heat losses are present, whether global or to a surface.

Buckmaster, J.↗

NASA NTRS1993

Stormtime transport of ring current and radiation belt ions

A dynamical guiding-center simulation model is used to study the stormtime ion transport which leads to the formation of the ring current and diffusion in the radiation belts. Representative ions guiding-center motion in response to model storm-associated impulses in the convection electric field is traced for a range of ion mu values. The present numerical results are compared with previously formulated limiting idealization of particle transport in order to assess the limits of validity of these approximations. For ions having drift periods that exceed the duration of the main phase of the storm, their inward transport to form the stormtime ring current is appropriately described as direct convective access. For ions having drift periods comparable to the duration of the main phase of the storm, there is a transition between direct convective access and transport that resembles radial diffusion. Lower-energy ring-current ions at L of about 3 are freshly injected there from open adiabatic trajectories, whereas the higher-energy ring-current population consists of a mixture of freshly injected and previously trapped ions.

Chen, Margaret W.↗

NASA NTRS1971

Recent radial turbine research at the NASA Lewis Research Center

The high efficiencies of small radial turbines led to their application in space power systems and numerous APU and shaft power engines. Experimental and analytical work associated with these systems included examination of blade-shroud clearance, blade loading, and exit diffuser design. Results indicate high efficiency over a wide range of specific speed and also insensitivity to clearance and blade loading in the radial part of the rotor. The exit diffuser investigation indicated that a conventional conical outer wall may not provide the velocity variation consistent with minimum overall diffuser loss.

Rohlik, H. E.↗

NASA NTRS1981

A general method for computing the total solar radiation force on complex spacecraft structures

The method circumvents many of the existing difficulties in computational logic presently encountered in the direct analytical or numerical evaluation of the appropriate surface integral. It may be applied to complex spacecraft structures for computing the total force arising from either specular or diffuse reflection or even from non-Lambertian reflection and re-radiation.

Chan, F. K.↗

NASA NTRS1981

Numerical experiments on the role of radiative processes in the development and maintenance of upper level clouds

A time-dependent, two-dimensional Eulerian model is presented whose purpose is to obtain more realistic parameterizations of extended high level cloudiness, and the results of a numerical experiment using the model are reported. The model is anelastic and the Bousinesque assumption is invoked. Unresolved subgrid scale processes are parameterized as eddy diffusion processes. Two phases of water are incorporated and equilibrium between them is assumed. The effects of infrared radiative processes are parametrically represented. Two simulations were conducted with identical initial conditions; in one of them, the radiation term was never turned on. The mean values of perturbation potential temperature at each level in the domain are plotted versus height after 15, 30, and 60 minutes of simulated time. The influence of the radiative term is seen to impose a cooling trend, leading to an increased generation of ice water and an increased generation of turbulent kinetic energy in the cloud layer.

Starr, D. O'C.↗

NASA NTRS1982

Single-mode theory of diffusive layers in thermohaline convection

A two-layer configuration of thermohaline convection is studied, with the principal aim of explaining the observed independence of the buoyancy-flux ratio on the stability parameter when the latter is large. Temperature is destabilizing and salinity is stabilizing, so diffusive interfaces separate the convecting layers. The convection is treated in the single-mode approximation, with a prescribed horizontal planform and wavenumber. Surveys of numerical solutions are presented for a selection of Rayleigh numbers R, stability parameters lambda and horizontal wavenumbers. The solutions yield a buoyancy flux ratio chi that is insensitive to lambda, in accord with laboratory experiments. However chi increases with increasing R, in contradiction to laboratory observations.

Gough, D. O.↗

NASA NTRS1985

Features of a reattaching turbulent shear layer in divergent channel flow

Experimental data have been obtained in an incompressible turbulent flow over a rearward-facing step in a diverging channel flow. Mean velocities, Reynolds stresses, and triple products that were measured by a laser Doppler velocimeter are presented for two cases of tunnel wall divergence. Eddy viscosities, production, convection, turbulent diffusion, and dissipation (balance of kinetic energy equation) terms are extracted from the data. These data are compared with various eddy-viscosity turbulence models. Numerical calculations incorporating the k-epsilon and algebraic-stress turbulence models are compared with the data. When determining quantities of engineering interest, the modified algebraic-stress model (ASM) is a significant improvement over the unmodified ASM and the unmodified k-epsilon model; however, like the others, it dramatically overpredicts the experimentally determined dissipation rate.

Driver, D. M.↗

NASA NTRS1986

Middle Atmosphere Program. Handbook for MAP, volume 20

Various topics related to investigations of the middle atmosphere are discussed. Numerical weather prediction, performance characteristics of weather profiling radars, determination of gravity wave and turbulence parameters, case studies of gravity-wave propagation, turbulence and diffusion due to gravity waves, the climatology of gravity waves, mesosphere-stratosphere-troposphere radar, antenna arrays, and data management techniques are among the topics discussed.

Bowhill, S. A.↗

NASA NTRS1992

Navier-Stokes analysis and experimental data comparison of compressible flow within ducts

Many aircraft employ ducts with centerline curvature or changing cross-sectional shape to join the engine with inlet and exhaust components. S-ducts convey air to the engine compressor from the intake and often decelerate the flow to achieve an acceptable Mach number at the engine compressor by increasing the cross-sectional area downstream. Circular-to-rectangular transition ducts are used on aircraft with rectangular exhaust nozzles to connect the engine and nozzle. To achieve maximum engine performance, the ducts should minimize flow total pressure loss and total pressure distortion at the duct exit. Changes in the curvature of the duct centerline or the duct cross-sectional shape give rise to streamline curvature which causes cross stream pressure gradients. Secondary flows can be caused by deflection of the transverse vorticity component of the boundary layer. This vortex tilting results in counter-rotating vortices. Additionally, the adverse streamwise pressure gradient caused by increasing cross-sectional area can lead to flow separation. Vortex pairs have been observed in the exit planes of both duct types. These vortices are due to secondary flows induced by pressure gradients resulting from streamline curvature. Regions of low total pressure are produced when the vortices convect boundary layer fluid into the main flow. The purpose of the present study is to predict the measured flow field in a diffusing S-duct and a circular-to-rectangular transition duct with a full Navier-Stokes computer program, PARC3D, and to compare the numerical predictions with new detailed experimental measurements. The work was undertaken to extend previous studies and to provide additional CFD validation data needed to help model flows with strong secondary flow and boundary layer separation. The S-duct computation extends the study of Smith et al, and Harloff et al, which concluded that the computation might be improved by using a finer grid and more advanced turbulence models. The present study compares results for both the Baldwin-Lomas and k-epsilon turbulence models and is conducted with a refined grid. For the transition duct, two inlet conditions were considered, the first with straight flow and the second with swirling flow. The first case permits examination of the effects of the geometric transition on the flow field, while the second case includes the rotational flow effect characteristic of a gas turbine engine.

Harloff, G. J.↗