Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equation”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Calculation of the flow field in supersonic mixed-compression inlets at angle of attack using the three-dimensional method of characteristics with discrete shock wave fitting

The influence of molecular transport is included in the computation by treating viscous and thermal diffusion terms in the governing partial differential equations as correction terms in the method of characteristics scheme. The development of a production type computer program is reported which is capable of calculating the flow field in a variety of axisymmetric mixed-compression aircraft inlets. The results agreed well with those produced by the two-dimensional method characteristics when axisymmetric flow fields are computed. For three-dimensional flow fields, the results agree well with experimental data except in regions of high viscous interaction and boundary layer removal.

Vadyak, J.↗
NASA NTRS1980

A simulation of the atmospheric cloud physics laboratory to aid in its design and the design of the experiments within the laboratory

Using the finite difference method with overrelaxation, numerical solutions of the steady-state vorticity transport equation were obtained for a continuous flow diffusion chamber of the Hudson-Squires type. The calculation neglected the effects due to temperature, gravity, and saturation. The size and shape of the manifold used to inject the aerosol laden flow were varied to obtain a design which would improve the performance of the chamber from strictly low Reynolds number (less than 20) fluid dynamical considerations.

Winchester, L. W., Jr.↗
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 NTRS1987

Simulation of the electron acoustic instability for a finite-size electron beam system

Simulations of the electron acoustic instability for a finite-size electron beam system are performed with a particle-in-cell code to investigate the heating phenomena associated with the instability and the width of the heating region. The results show that the beam radiates electrostatic electron acoustic waves whose decay time outside the beam agrees with the spatial decay length derived from the linear dispersion equation. The ambient cold electrons in a diffusion region surrounding the beam are heated to a higher temperature by absorbing the radiated electron acoustic waves, with the heating occurring mainly in the parallel direction. In the heat diffusion region, the temperature of the cold electrons decreases with distance from the beam with a temperature gradient length smaller than the decay length of the wave energy. The results are discussed with respect to the DE 1 plasma and wave observations in the polar cusp region.

Lin, C. S.↗
NASA NTRS1988

Explicit solution of the spectral radiance in integrating spheres with application to the Earth Radiation Budget Experiment ground calibration

An explicit solution of the spectral radiation leaving an arbitrary point on the wall of a spherical cavity with diffuse reflectivity is obtained. A general measurement equation is obtained that describes the output of a sensor measuring the sphere output within a given field of view and with specified angular and spectral responses. The results are applied to the Earth Radiation Budget Experiment (ERBE) integrating sphere. It is shown that by using appropriate interpretation and processing, a high-accuracy shortwave calibration of the ERBE sensors can be achieved.

Halyo, Nesim↗
NASA NTRS1989

Numerical simulation of axisymmetric turbulent flow in combustors and diffusers

Numerical studies of turbulent flow in an axisymmetric 45-deg-expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k-epsilon model were solved in a nonorthogonal curvilinear coordinate system. A zonal-grid method, where the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. Boundary overlap and interpolating techniques were used, and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite-differencing methods (hybrid, quadratic upwind, and skew upwind) were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.

Yung, Chain-Nan↗
NASA NTRS1990

Modeling of near-wall turbulence

An improved k-epsilon model and a second order closure model is presented for low Reynolds number turbulence near a wall. For the k-epsilon model, a modified form of the eddy viscosity having correct asymptotic near wall behavior is suggested, and a model for the pressure diffusion term in the turbulent kinetic energy equation is proposed. For the second order closure model, the existing models are modified for the Reynolds stress equations to have proper near wall behavior. A dissipation rate equation for the turbulent kinetic energy is also reformulated. The proposed models satisfy realizability and will not produce unphysical behavior. Fully developed channel flows are used for model testing. The calculations are compared with direct numerical simulations. It is shown that the present models, both the k-epsilon model and the second order closure model, perform well in predicting the behavior of the near wall turbulence. Significant improvements over previous models are obtained.

Shih, T. H.↗
NASA NTRS1990

Estimation of the processes controlling variability in phytoplankton pigment distributions on the southeastern U.S. continental shelf

Five coastal-zone-color-scanner images from the southeastern U.S. continental shelf are combined with concurrent moored current meter measurements to assess the processes controlling the variability in chlorophyll concentration and distribution in this region. An equation governing the space and time distribution of a nonconservative quantity such as chlorophyll is used in the calculations. The terms of the equation, estimated from observations, show that advective, diffusive, and local processes contribute to the plankton distributions and vary with time and location. The results from this calculation are compared with similar results obtained using a numerical physical-biological model with circulation fields derived from an optimal interpolation of the current meter observations and it is concluded that the two approaches produce different estimates of the processes controlling phytoplankton variability.

Mcclain, Charles R.↗
NASA NTRS1992

Direct Coupling Method for Time-Accurate Solution of Incompressible Navier-Stokes Equations

A noniterative finite difference numerical method is presented for the solution of the incompressible Navier-Stokes equations with second order accuracy in time and space. Explicit treatment of convection and diffusion terms and implicit treatment of the pressure gradient give a single pressure Poisson equation when the discretized momentum and continuity equations are combined. A pressure boundary condition is not needed on solid boundaries in the staggered mesh system. The solution of the pressure Poisson equation is obtained directly by Gaussian elimination. This method is tested on flow problems in a driven cavity and a curved duct.

Soh, Woo Y.↗
NASA NTRS1992

Variable gamma computations for a forward-firing Space Shuttle Main Engine

A procedure for solving an axisymmetric SSME nozzle exhausting with an equivalent mass flow of two nozzles is presented for the return-to-launch-site (RTLS) maneuver for an engine-out case. The algorithm incorporates a two-species variable gamma formulation along with diffusion terms associated with a species conservation equation. Coupling of this algorithm is accomplished explicitly to the Navier-Stokes PARC2DR and PARC3DR codes. Axisymmetric results indicate a very complicated flowfield with a bow shock extending to approximately 570 radii in front of the nozzle, and low temperatures on what would be an external tank (ET) aft of the nozzle.

Saladino, A. J.↗
NASA NTRS1994

Calculation of Propulsive Nozzle Flowfields in Multidiffusing Chemically Reacting Environments

An advanced engineering model has been developed to aid in the analysis and design of hydrogen/oxygen chemical rocket engines. The complete multispecies, chemically reacting and multidiffusing Navier-Stokes equations are modelled, including the Soret thermal diffusion and the Dufour energy transfer terms. In addition to the spectrum of multispecies aspects developed, the model developed in this study is also conservative in axisymmetric flow for both inviscid and viscous flow environments and the boundary conditions employ a viscous, chemically reacting, reference plane characteristics method. Demonstration cases are presented for a 1030:1 area ratio nozzle, a 25 lbf film cooled nozzle, and a transpiration cooled plug and spool rocket engine. The results indicate that the thrust coefficient predictions of the 1030:1 and the 25 lbf film cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements when all of the chemical reaction and diffusion terms are considered. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. The Soret thermal diffusion term is demonstrated to have a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle in both the laminar flow 1030:1 nozzle and the turbulent flow plug and spool nozzle analysis cases performed. Further, the Soret term was shown to represent an important fraction of the diffusion fluxes occurring in a transpiration cooled rocket engine.

Kacynski, Kenneth John↗
NASA NTRS1995

The applicability of turbulence models to aerodynamic and propulsion flowfields at McDonnell-Douglas Aerospace

The objective of this viewgraph presentation is to evaluate turbulence models for integrated aircraft components such as the forebody, wing, inlet, diffuser, nozzle, and afterbody. The one-equation models have replaced the algebraic models as the baseline turbulence models. The Spalart-Allmaras one-equation model consistently performs better than the Baldwin-Barth model, particularly in the log-layer and free shear layers. Also, the Sparlart-Allmaras model is not grid dependent like the Baldwin-Barth model. No general turbulence model exists for all engineering applications. The Spalart-Allmaras one-equation model and the Chien k-epsilon models are the preferred turbulence models. Although the two-equation models often better predict the flow field, they may take from two to five times the CPU time. Future directions are in further benchmarking the Menter blended k-w/k-epsilon and algorithmic improvements to reduce CPU time of the two-equation model.

Kral, Linda D.↗
NASA NTRS2001

Computational Aerothermodynamics in Aeroassist Applications

Aeroassisted planetary entry uses atmospheric drag to decelerate spacecraft from super-orbital to orbital or suborbital velocities. Numerical simulation of flow fields surrounding these spacecraft during hypersonic atmospheric entry is required to define aerothermal loads. The severe compression in the shock layer in front of the vehicle and subsequent, rapid expansion into the wake are characterized by high temperature, thermo-chemical nonequilibrium processes. Implicit algorithms required for efficient, stable computation of the governing equations involving disparate time scales of convection, diffusion, chemical reactions, and thermal relaxation are discussed. Robust point-implicit strategies are utilized in the initialization phase; less robust but more efficient line-implicit strategies are applied in the endgame. Applications to ballutes (balloon-like decelerators) in the atmospheres of Venus, Mars, Titan, Saturn, and Neptune and a Mars Sample Return Orbiter (MSRO) are featured. Examples are discussed where time-accurate simulation is required to achieve a steady-state solution.

Gnoffo, Peter A.↗
NASA NTRS2014

Comptonization in Ultra-Strong Magnetic Fields: Numerical Solution to the Radiative Transfer Problem

We consider the radiative transfer problem in a plane-parallel slab of thermal electrons in the presence of an ultra-strong magnetic field (B approximately greater than B(sub c) approx. = 4.4 x 10(exp 13) G). Under these conditions, the magnetic field behaves like a birefringent medium for the propagating photons, and the electromagnetic radiation is split into two polarization modes, ordinary and extraordinary, that have different cross-sections. When the optical depth of the slab is large, the ordinary-mode photons are strongly Comptonized and the photon field is dominated by an isotropic component. Aims. The radiative transfer problem in strong magnetic fields presents many mathematical issues and analytical or numerical solutions can be obtained only under some given approximations. We investigate this problem both from the analytical and numerical point of view, provide a test of the previous analytical estimates, and extend these results with numerical techniques. Methods. We consider here the case of low temperature black-body photons propagating in a sub-relativistic temperature plasma, which allows us to deal with a semi-Fokker-Planck approximation of the radiative transfer equation. The problem can then be treated with the variable separation method, and we use a numerical technique to find solutions to the eigenvalue problem in the case of a singular kernel of the space operator. The singularity of the space kernel is the result of the strong angular dependence of the electron cross-section in the presence of a strong magnetic field. Results. We provide the numerical solution obtained for eigenvalues and eigenfunctions of the space operator, and the emerging Comptonization spectrum of the ordinary-mode photons for any eigenvalue of the space equation and for energies significantly lesser than the cyclotron energy, which is on the order of MeV for the intensity of the magnetic field here considered. Conclusions. We derived the specific intensity of the ordinary photons, under the approximation of large angle and large optical depth. These assumptions allow the equation to be treated using a diffusion-like approximation.

acceleration of particles↗
NASA NTRS2014

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

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

Dissolution of multicomponent bubbles

The behavior of an isolated, stationary, multicomponent gas bubble in a glassmelt containing several dissolved gases is considered. The relevant mass-transport equations are formulated and calculations are performed for the case of two diffusing gases using a quasi-stationary model and a numerical solution of the exact mass-transfer equations. The results obtained from these two approaches are compared. The factors which govern the dissolution or growth of a bubble are thermodynamic and kinetic in origin. The tendency of a bubble to grow or shrink at long times is controlled by departure from overall equilibrium, whereas the short-time bubble dynamics may be dominated by kinetic effects. As a result of the existence of these dual influences, maxima and/or minima occur in the functional dependence of the bubble radius on time.

Weinberg, M. C.↗
NASA NTRS2007

Multi-Scale Modeling of Magnetospheric Reconnection

One of the major challenges in modeling the magnetospheric magnetic reconnection is to quantify the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions near reconnection sites. There is still considerable debate as to what degree microphysical processes on kinetic scales affect the global evolution and how important it is to substitute numerical dissipation and/or ad hoc anomalous resistivity by a physically motivated model of dissipation. Comparative studies of magnetic reconnection in small scale geometries demonstrated that MHD simulations that included non-ideal processes in terms of a resistive term $\eta J$ did not produce the fast reconnection rates observed in kinetic simulations. For a broad range of physical parameters in collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. We utilize the global MHD code BATSRUS and incorporate nongyrotropic effects in diffusion regions in terms of corrections to the induction equation. We developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated selfconsistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites is updated during the simulations. To clarify the role of nongyrotropic effects in diffusion region on the global magnetospheric dynamic we perform simulations with steady southward IMF driving of the magnetosphere. Ideal MHD simulations with magnetic reconnection supported by numerical resistivity produce steady configuration with almost stationary near-earth neutral line (NENL). Simulations with non-gyrotropic corrections demonstrate dynamic quasi-periodic response to the steady driving condition. The loading/unloading cycle in non-gyrotropic MHD results has a non-stationary reconnection site in the magnetotail, with the retreating during the stretching phase and then a new NENL forming in the resulting thin plasma sheet. We expect that this model will lead to improved representations of space weather event in the magnetosphere.

Kuznetsova, M. M.↗