Results for “Diffusion equation” — Page 71

NASA NTRS2000

Descending Layer Variability Over Arecibo

Descending layers of ionization over Arecibo exhibit very diverse behavior from night to night that does not appear to be strongly correlated to geomagnetic activity, solar forcing, or average semidiurnal tidal winds. On some nights, three or more distinct layers are observed to form near 170 km over timescales of approx. 2 hours. Rather than descending smoothly over periods of several hours, these layers stall, abruptly disappear, or even reverse direction in the midst of their descent. The time scales for their disappearance are examined and compared to loss rates arising from diffusion and recombination. Diffusion alone is found to be too slow to account for the observations, but recombination is fast enough provided that the convergent wind shear that forms the layer is relatively weak coincident with their disappearance. The continuity equation is solved in conjunction with a time sequence of radar profiles to estimate the vertical drift and horizontal neutral wind consistent with the observed behavior. The resultant wind field is northward, has an average speed of approx. 80 m/ s, and varies significantly near the altitude where the layers are observed. These inferred winds are consistent with the presence of the observed layers, and their magnitudes as obtained from the classical continuity and momentum equations are reasonable for this altitude range.

Earle, G. D.↗

NASA NTRS1983

A finite element method for diffusion dominated unsteady viscous flows

A general conforming finite element scheme for computing viscous flows is presented which is of second-order accuracy in space and time. Viscous terms are treated implicitly and advection terms are treated explicitly in the time marching segment of the algorithm. A method for solving the algebraic equations at each time step is given. The method is demonstrated on two test problems, one of them being a plane vortex flow for which asymptotic methods are used to obtain suitable numerical boundary conditions at each time step.

Gunzburger, M. D.↗

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 NTRS1973

A study of the force-field equation for the propagation of galactic cosmic rays

A new development is given of the solution of the equation of the force-field approximation for the propagation of galactic cosmic rays in the interplanetary region. It leads to simpler methods for determining the force-field parameters. A method is given for determining the separable diffusion coefficient from observations of galactic electron spectrum and near-earth electron spectra; it is shown that this diffusion coefficient is not unique but may have a periodic-like dependence upon rigidity; and the method is used to obtain diffusion coefficients for 1965 and 1968. Approximate formulae relating small changes in intensity and diffusion coefficient are developed and some applications of these noted; in one it is shown that the form of, and changes in, diffusion coefficient deduced previously for a neutron monitor event during June-September 1969 are unnecessarily constrained and therefore probably not correct.

Gleeson, L. J.↗

NASA NTRS1974

Analytic approximations in the study of the solar modulation of electrons

Numerical solutions to the transport equation of galactic cosmic rays in the interplanetary medium have been used to investigate the applicability of commonly used approximate analytic solutions for electrons (10 MeV to 10 GeV). We find that for a given cosmic-ray diffusion coefficient the force-field approximation is in reasonable agreement with the numerical solution at energies above 200 MeV, but deviates significantly at lower energies, depending on the shape of the interstellar electron spectrum. The diffusion-convection approximation agrees generally with the numerical solution within a factor of 2 over the entire energy range.

Cummings, A. C.↗

NASA NTRS1976

Theoretical investigations of plasma processes

System analyses are presented for electrically sustained, collision dominated plasma centrifuges, in which the plasma rotates under the influence of the Lorentz forces resulting from the interaction of the current density fields with an external magnetic field. It is shown that gas discharge centrifuges are technically feasible in which the plasma rotates at speeds up to 1 million cm/sec. The associated centrifugal forces produce a significant spatial isotope separation, which is somewhat perturbed in the viscous boundary layers at the centrifuge walls. The isotope separation effect is the more pronounced. The induced magnetic fields have negligible influence on the plasma rotation if the Hall coefficient is small. In the technical realization of collision dominated plasma centrifuges, a trade-off has to be made between power density and speeds of rotation. The diffusion of sputtered atoms to system surfaces of ion propulsion systems and the deposition of the atoms are treated theoretically by means of a simple model which permits an analytical solution. The problem leads to an inhomogeneous integral equation.

Wilhelm, H. E.↗

NASA NTRS1982

Convective thinning of the lithosphere - A mechanism for the initiation of continental rifting

A model of lithospheric thinning, in which heat is convected to the base and conducted within the lithosphere, is presented. An analytical equation for determinining the amount of thinning attainable on increasing the heat flux from the asthenosphere is derived, and a formula for lithosphere thickness approximations as a function of time is given. Initial and final equilibrium thicknesses, thermal diffusivity, transition temperature profile, and plume temperature profile are all factors considered for performing rate of thinning determinations. In addition, between initial and final equilibrium states, lithospheric thinning occurs at a rate which is inversely proportional to the square root of the time. Finally, uplift resulting from thermal expansion upon lithospheric thinning is on the order of 10 to the 2nd to 10 to the 3rd m.

Spohn, T.↗

NASA NTRS1985

Investigation of parabolic computational techniques for internal high-speed viscous flows

A feasibility study was conducted to assess the applicability of an existing parabolic analysis (ADD-Axisymmetric Diffuser Duct), developed previously for subsonic viscous internal flows, to mixed supersonic/subsonic flows with heat addition simulating a SCRAMJET combustor. A study was conducted with the ADD code modified to include additional convection effects in the normal momentum equation when supersonic expansion and compression waves were present. It is concluded from the present study that for the class of problems where strong viscous/inviscid interactions are present a global iteration procedure is required.

Anderson, O. L.↗

NASA NTRS1985

Numerical investigation of internal high-speed viscous flows using a parabolic technique

A feasibility study has been conducted to assess the applicability of an existing parabolic analysis (ADD-Axisymmetric Diffuser Duct), developed previously for subsonic viscous internal flows, to mixed supersonic/subsonic flows with heat addition simulating a SCRAMJET combustor. A study was conducted with the ADD code modified to include additional convection effects in the normal momentum equation when supersonic expansion and compression waves are present. A set of test problems with weak shock and expansion waves have been analyzed with this modified ADD method and stable and accurate solutions were demonstrated provided the streamwise step size was maintained at levels larger than the boundary layer displacement thickness. Calculations made with further reductions in step size encountered departure solutions consistent with strong interaction theory. Calculations were also performed for a flow field with a flame front in which a specific heat release was imposed to simulate a SCRAMJET combustor. In this case the flame front generated relatively thick shear layers which aggravated the departure solution problem. Qualitatively correct results were obtained for these cases using a marching technique with the convective terms in the normal momentum equation suppressed. It is concluded from the present study that for the class of problems where strong viscous/inviscid interactions are present a global iteration procedure is required.

Anderson, O. L.↗

NASA NTRS1988

Evaluation of surface energy and radiation balance systems for FIFE

The energy balance and radiation balance components were determined at six sites during the First International Satellite Land Surface Climatology Project Field Experiment (FIFE) conducted south of Manhattan, Kansas during the summer of 1987. The objectives were: to determine the effect of slope and aspect, throughout a growing season, on the magnitude of the surface energy balance fluxes as determined by the Energy Balance Method (EBM); to investigate the calculation of the soil heat flux density at the surface as calculated from the heat capacity and the thermal conductivity equations; and to evaluate the performance of the Surface Energy and Radiation Balance System (SERBS). A total of 17 variables were monitored at each site. They included net, solar (up and down), total hemispherical (up and down), and diffuse radiation, soil temperature and heat flux density, air and wet bulb temperature gradients, wind speed and direction, and precipitation. A preliminary analysis of the data, for the season, indicate that variables including net radiation, air temperature, vapor pressure, and wind speed were quite similar at the sites even though the sites were as much as 16 km apart and represented four cardinal slopes and the top of a ridge.

Fritschen, Leo J.↗

NASA NTRS1987

Fermi acceleration at shocks with arbitrary velocity profiles

Diffusive particle acceleration in modified shock fronts is considered. A method by means of which to obtain the resulting particle spectrum for an arbitrary flow profile is outlined; this method is based on a reformulation of the transport equation as a nonlinear first-order equation. Besides the shape of the flow profile, the spectral index of particles depends on the Peclet number of the shock. Results are given for some sample profiles.

Schneider, Peter↗

NASA NTRS1994

Probability density function approach for compressible turbulent reacting flows

The objective of the present work is to extend the probability density function (PDF) tubulence model to compressible reacting flows. The proability density function of the species mass fractions and enthalpy are obtained by solving a PDF evolution equation using a Monte Carlo scheme. The PDF solution procedure is coupled with a compression finite-volume flow solver which provides the velocity and pressure fields. A modeled PDF equation for compressible flows, capable of treating flows with shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed. Two super sonic diffusion flames are studied using the proposed PDF model and the results are compared with experimental data; marked improvements over solutions without PDF are observed.

Hsu, A. T.↗

NASA NTRS1993

Particle-Gas Dynamics in the Midplane of a Protoplanetary Nebula

In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimeter-to-meter sized primordial aggregates. During this stage, particles are able to settle toward the midplane into a layer of mass density comparable to or much greater than that of the gas. Nonlinear, coupled interactions between the particles and the nebula gas become significant and ultimately determine the vertical profiles of the particle density and the mean velocities of the particles and the nearby gas. This is the environment in which the earliest planetesimals probably form. Our numerical models rely on the Reynolds averaged Navier-Stokes equations for the gas and particles and are fully viscous and turbulent (and, for the particles. compressible). Our turbulence modeling uses a Prandtl local shear parametrization, validated by laboratory experiments. We have developed a new model for particle diffusion involving the profile of the nebula gas turbulence and the particle Schmidt number, which is a function of particle size and density. We model a cool, quiescent nebula at 1 AU (280 K) and 10 AU (90 K) and a higher temperature stage (1000 K) at 1 AU. Our main results include: (a) rapid accretion of planetesimals by gravitationally unstable fragmentation on an orbital time scale (the 'Goldreich-Ward instability') is unlikely to occur until objects have already accreted by some other process to the mass of the largest known meteorite samples, if at all; (b) from 'seeds' as small as 10 m, growth of 10- to 100-km planetesimals can proceed rapidly by drift-augmented accretion in the particle-rich midplane with orbital decay of about 1% for the growing planetesimals; (c) outward transport of vapor and small entrained chips can account for significant radial compositional and mineralogical mixing in primitive meteorite parent bodies. Other implications may be drawn concerning particle random collisional velocities and the time scales of midplane turbulence.

Cuzzi, Jeffrey N.↗

NASA NTRS2011

3D Modeling of Ultrasonic Wave Interaction with Disbonds and Weak Bonds

Ultrasonic techniques, such as the use of guided waves, can be ideal for finding damage in the plate and pipe-like structures used in aerospace applications. However, the interaction of waves with real flaw types and geometries can lead to experimental signals that are difficult to interpret. 3-dimensional (3D) elastic wave simulations can be a powerful tool in understanding the complicated wave scattering involved in flaw detection and for optimizing experimental techniques. We have developed and implemented parallel 3D elastodynamic finite integration technique (3D EFIT) code to investigate Lamb wave scattering from realistic flaws. This paper discusses simulation results for an aluminum-aluminum diffusion disbond and an aluminum-epoxy disbond and compares results from the disbond case to the common artificial flaw type of a flat-bottom hole. The paper also discusses the potential for extending the 3D EFIT equations to incorporate physics-based weak bond models for simulating wave scattering from weak adhesive bonds.

Leckey, C.↗

NASA NTRS2020

Helicity in the Large-Scale Galactic Magnetic Field

We search for observational signatures of magnetic helicity in data from all-sky radio polarization surveys of the Milky Way Galaxy. Such a detection would help confirm the dynamo origin of the field and may provide new observational constraints for its shape. We compare our observational results to simulated observations for both a simple helical field, and for a more complex field that comes from a solution to the dynamo equation. Our simulated observations show that the large-scale helicity of a magnetic field is reflected in the large-scale structure of the fractional polarization derived from the observed synchrotron radiation and Faraday depth of the diffuse Galactic synchrotron emission. Comparing the models with the observations provides evidence for the presence of a quadrupolar magnetic field with a vertical component that is pointing away from the observer in both hemispheres of the Milky Way Galaxy. Since there is no reason to believe that the Galactic magnetic field is unusual when compared to other galaxies, this result provides further support for the dynamo origin of large-scale magnetic fields in galaxies.

J. L. West↗

NASA NTRS2015

A Numerical and Experimental Study of Coflow Laminar Diffusion Flames: Effects of Gravity and Inlet Velocity

In this work, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting behavior of laminar coflow methane-air diffusion flames was investigated both computationally and experimentally. A series of flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) was assessed numerically under microgravity and normal gravity conditions with the fuel stream CH4 mole fraction ranging from 0.4 to 1.0. Computationally, the MC-Smooth vorticity-velocity formulation of the governing equations was employed to describe the reactive gaseous mixture; the soot evolution process was considered as a classical aerosol dynamics problem and was represented by the sectional aerosol equations. Since each flame is axisymmetric, a two-dimensional computational domain was employed, where the grid on the axisymmetric domain was a nonuniform tensor product mesh. The governing equations and boundary conditions were discretized on the mesh by a nine-point finite difference stencil, with the convective terms approximated by a monotonic upwind scheme and all other derivatives approximated by centered differences. The resulting set of fully coupled, strongly nonlinear equations was solved simultaneously using a damped, modified Newton's method and a nested Bi-CGSTAB linear algebra solver. Experimentally, the flame shape, size, lift-off height, and soot temperature were determined by flame emission images recorded by a digital camera, and the soot volume fraction was quantified through an absolute light calibration using a thermocouple. For a broad spectrum of flames in microgravity and normal gravity, the computed and measured flame quantities (e.g., temperature profile, flame shape, lift-off height, and soot volume fraction) were first compared to assess the accuracy of the numerical model. After its validity was established, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting tendency of laminar coflow methane-air diffusion flames was explored further by examining quantities derived from the computational results.

microgravity↗

NASA NTRS1977

The NASA Ames Research Center one- and two-dimensional stratospheric models. Part 1: The one-dimensional model

A one-dimensional model of stratospheric trace constituents is described in detail. Specifically, the numerical solution of the species continuity equations, including a technique for treating the stiff differential equations representing the chemical kinetic terms, and an appropriate method for simulating the diurnal variations of the species concentrations, are discussed. A specialized treatment of atmospheric photodissociation rates is outlined in the text. The choice of a vertical eddy diffusivity profile and its success in predicting the vertical tracer distributions (carbon 14, methane, and nitrous oxide) are also discussed.

Turco, R. P.↗