Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical transport”

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 397 records · Page 22

Tree-based solvers for adaptive mesh refinement code FLASH – III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources

Radiation is an important contributor to the energetics of the interstellar medium, yet its transport is difficult to solve numerically. We present a novel approach towards solving radiative transfer of diffuse sources via backwards ray tracing. Here, we focus on the radiative transfer of infrared radiation and the radiation pressure on dust. The new module, TREERAY/RADPRESSURE, is an extension to the novel radiative transfer method TREERAY implemented in the grid-based Magneto-Hydrodynamics code FLASH. In TREERAY/RADPRESSURE, every cell and every star particle is a source of infrared radiation. We also describe how gas, dust, and radiation are coupled via a chemical network. This allows us to compute the local dust temperature in thermal equilibrium, leading to a significantly improvement over the classical grey approximation. In several tests, we demonstrate that the scheme produces the correct radiative intensities as well as the correct momentum input by radiation pressure. Subsequently, we apply our new scheme to model massive star formation from a collapsing, turbulent core of 150 M ⊙ . We include the effects of both, ionizing and infrared radiation on the dynamics of the core. We find that the newborn massive star prevents fragmentation in its proximity due to radiative heating. Over time, dust and radiation temperature equalize, while the gas temperature can be either warmer due to shock heating or colder due to insufficient dust–gas coupling. Compared to gravity, the effects of radiation pressure are insignificant for the stellar mass on the simulated time-scale in this work.

(ISM:) HII regions↗

A study of efficiency in low resistivity silicon solar cells

A general device-analysis program has been utilized to study the efficiency of silicon solar cells. The analysis is applied to specific geometries of both n(+)-p and n(+)-p-p(+) solar cells, and involves a numerical solution of the basic transport and continuity equations. This approach allows solutions that are free of typical limiting assumptions involved in solving the device equations as well as solutions relating to lifetime, mobility variations, and diffused-region profiles. The analysis includes available empirical information on diffusion length, mobility, and lifetime as a function of doping, as well as a Gaussian profile for the diffused region. Results are presented which illustrate the limitations of efficiency as a function of doping. It is found that the maximum efficiencies for both types of cell converge at lower resistivities to around 16% with air-mass-zero radiation and a single-layer absorbing-SiO antireflecting film. It is also found that the minority-carrier lifetime, both in the n(+) surface and p-type bulk regions, presents serious limitations to conversion efficiency, particularly in the low-resistivity cells.

Dunbar, P. M.↗

Diurnal behaviour of water on Mars

A numerical model of the diurnal transport of water across the Martian surface is developed. The atmospheric boundary layer is modeled in terms of local radiative-convective processes, and radiative effects of ice fogs near the surface are included. Diffusion of water in the ground is treated for the cases of adsorption and condensation. The model is applied to the diurnal variation of water vapor in the atmosphere as observed by Barker (1974). The morning rise in the amount of water vapor can be explained in terms of the evaporation of ground fogs. The evening decrease is compatible with the model if adsorption dominates in the soil. The average level of vapor concentration requires that the atmosphere above the boundary layer be relatively dry. The ground fogs persist until midmorning and should be observable. Some consequences of these conclusions are discussed.

Flasar, F. M.↗

Survey of computer programs for prediction of crash response and of its experimental validation

The author seeks to critically assess the potentialities of the mathematical and hybrid simulators which predict post-impact response of transportation vehicles. A strict rigorous numerical analysis of a complex phenomenon like crash may leave a lot to be desired with regard to the fidelity of mathematical simulation. Hybrid simulations on the other hand which exploit experimentally observed features of deformations appear to hold a lot of promise. MARC, ANSYS, NONSAP, DYCAST, ACTION, WHAM II and KRASH are among some of the simulators examined for their capabilities with regard to prediction of post impact response of vehicles. A review of these simulators reveals that much more by way of an analysis capability may be desirable than what is currently available. NASA's crashworthiness testing program in conjunction with similar programs of various other agencies, besides generating a large data base, will be equally useful in the validation of new mathematical concepts of nonlinear analysis and in the successful extension of other techniques in crashworthiness.

Kamat, M. P.↗

Auroral electron interaction with the atmosphere in the presence of conjugate field-aligned electrostatic potentials

In the present investigation, both pitch angle and energy distributions of auroral electrons are obtained, for the first time, by rigorous numerical solutions of a multiangle transport equation at all energies in a realistic three component atmosphere in the presence of a conjugate field-aligned electrostatic potential. The results confirm all previous demonstrations of the validity of Evan's (1974, 1975, 1976) hypothesis concerning the origin of low energy auroral electrons from collisional interactions of the auroral primaries with the atmosphere. The existence of a collisionally populated field-perpendicular flux at energies near the beam energy is also demonstrated.

Prasad, S. S.↗

The Martian climate: Energy balance models with CO2/H2O atmospheres

Coupled equations are developed for mass and heat transport in a seasonal Mars model with condensation and sublimation of CO2 at the polar caps. Topics covered include physical considerations of planetary as mass and energy balance; effects of phase changes at the surface on mass and heat flux; atmospheric transport and governing equations; and numerical analysis.

Hoffert, M. I.↗

Energetic photoelectrons and the polar rain

In the daytime midlatitudes, the Low Altitude Plasma Instrument (LAPI) on board the Dynamics Explorer 2 satellite has observed photoelectrons with energies as high as 850 eV. These energetic photoelectrons are an extension of the 'classical' photoelectrons (less than 60 eV) and result from photoionization of neutrals by soft solar X-rays. Since these photoelectrons are produced wherever the solar flux is incident on the earth's atmosphere, they should be present in sunlit polar cap. But in the polar cap, over these same energies, there is a well-known electron population: the polar rain, a low intensity electron flux of magnetospheric origin. Thus, in the sunlit polar cap, an energetic population of electrons should consist of both an ionospheric (photoelectron) and a magnetospheric (polar rain) component. Using numerical solutions of an electron transport equation with appropriate boundary conditions and sunlit polar cap LAPI data, it is shown that the two populations (photoelectron and polar rain) are indeed present and are both needed to explain polar cap observations.

Decker, Dwight T.↗

Multidisciplinary optimization in aircraft design using analytic technology models

An approach to multidisciplinary optimization is presented which combines the Global Sensitivity Equation method, parametric optimization, and analytic technology models. The result is a powerful yet simple procedure for identifying key design issues. It can be used both to investigate technology integration issues very early in the design cycle, and to establish the information flow framework between disciplines for use in multidisciplinary optimization projects using much more computational intense representations of each technology. To illustrate the approach, an examination of the optimization of a short takeoff heavy transport aircraft is presented for numerous combinations of performance and technology constraints.

Malone, Brett↗

Magnetic radiation shielding - An idea whose time has returned?

One solution to the problem of shielding crew from particulate radiation in space is to use active electromagnetic shielding. Practical types of shield include the magnetic shield, in which a strong magnetic field diverts charged particles from the crew region, and the magnetic/electrostatic plasma shield, in which an electrostatic field shields the crew from positively charged particles, while a magnetic field confines electrons from the space plasma to provide charge neutrality. Advances in technology include high-strength composite materials, high-temperature superconductors, numerical computational solutions to particle transport in electromagnetic fields, and a technology base for construction and operation of large superconducting magnets. These advances make electromagnetic shielding a practical alternative for near-term future missions.

Landis, Geoffrey A.↗

A comparative study of advanced shock-capturing schemes applied to Burgers' equation

A systematic evaluation is conducted of all extant numerical schemes for nonlinear scalar transport problems, and several advanced shock-capturing schemes are used to solve the nonlinear Burgers' equation in order to characterize their ability to resolve the sharp discontinuity, expansion zone, and propagation and collision features of shocks. For discontinuous functions, the Warming-Beam scheme generates preshock wiggles, while the Lax-Wendroff scheme generates postshock ones. Such limiters as the MUSCL or the superbee are more compressive than minimod or monotonic limiters. The performance of such TVD schemes as the upwind, the symmetric, and the Roe-Sweby, resemble each other.

Yang, H. Q.↗

Particle acceleration in cosmic plasmas; Proceedings of the Workshop, Bartol Research Inst., Newark, DE, Dec. 4-6, 1991

The present volume on particle acceleration in cosmic plasmas discusses transport theory, particle acceleration in the sun, heliosphere, and galaxy, particle acceleration at nonrelativistic shocks, and stochastic particle acceleration. Attention is given to particle acceleration in relativistic flows and shocks, simulations, composition and sources of high-energy cosmic rays, electron acceleration, and observation priorities and perspectives. Topics addressed include particle acceleration in the magnetosphere, particle transport from a turbulence perspective, numerical simulations of time-dependent cosmic ray mediated shocks, and solar flare plasma. Also discussed are relativistic shock waves and the excitation of plerions, ion acceleration at collisionless shock interactions, ultrahigh-energy cosmic rays from Fanaroff Riley, and constraints on electron acceleration in the Crab Nebula.

Zank, G. P.↗

Rates for neutron-capture reactions on tungsten isotopes in iron meteorites

High-precision W isotopic analyses by Harper and Jacobsen indicate the W-182/W-183 ratio in the Toluca iron meteorite is shifted by -(3.0 +/- 0.9) x 10(exp -4) relative to a terrestrial standard. Possible causes of this shift are neutron-capture reactions on W during Toluca's approximately 600-Ma exposure to cosmic ray particles or radiogenic growth of W-182 from 9-Ma Hf-182 in the silicate portion of the Earth after removal of W to the Earth's core. Calculations for the rates of neutron-capture reactions on W isotopes were done to study the first possibility. The LAHET Code System (LCS) which consists of the Los Alamos High Energy Transport (LAHET) code and the Monte Carlo N-Particle(MCNP) transport code was used to numerically simulate the irradiation of the Toluca iron meteorite by galactic-cosmic-ray (GCR) particles and to calculate the rates of W(n, gamma) reactions. Toluca was modeled as a 3.9-m-radius sphere with the composition of a typical IA iron meteorite. The incident GCR protons and their interactions were modeled with LAHET, which also handled the interactions of neutrons with energies above 20 MeV. The rates for the capture of neutrons by W-182, W-183, and W-186 were calculated using the detailed library of (n, gamma) cross sections in MCNP. For this study of the possible effect of W(n, gamma) reactions on W isotope systematics, we consider the peak rates. The calculated maximum change in the normalized W-182/W-183 ratio due to neutron-capture reactions cannot account for more than 25% of the mass 182 deficit observed in Toluca W.

Masarik, J.↗

Numerical Modeling of High-Temperature Corrosion Processes

Numerical modeling of the diffusional transport associated with high-temperature corrosion processes is reviewed. These corrosion processes include external scale formation and internal subscale formation during oxidation, coating degradation by oxidation and substrate interdiffusion, carburization, sulfidation and nitridation. The studies that are reviewed cover such complexities as concentration-dependent diffusivities, cross-term effects in ternary alloys, and internal precipitation where several compounds of the same element form (e.g., carbides of Cr) or several compounds exist simultaneously (e.g., carbides containing varying amounts of Ni, Cr, Fe or Mo). In addition, the studies involve a variety of boundary conditions that vary with time and temperature. Finite-difference (F-D) techniques have been applied almost exclusively to model either the solute or corrodant transport in each of these studies. Hence, the paper first reviews the use of F-D techniques to develop solutions to the diffusion equations with various boundary conditions appropriate to high-temperature corrosion processes. The bulk of the paper then reviews various F-D modeling studies of diffusional transport associated with high-temperature corrosion.

Nesbitt, James A.↗

Computation of Turbulent Recirculating Flow in Channels, and for Impingement Cooling

Fully elliptic forms of the transport equations have been solved numerically for two flow configurations. The first is turbulent flow in a channel with transverse rectangular ribs, and the second is impingement cooling of a plane surface. Both flows are relevant to proposed designs for active cooling of hypersonic vehicles using supercritical hydrogen as the coolant. Flow downstream of an abrupt pipe expansion and of a backward-facing step were also solved with various near-wall turbulence models as benchmark problems. A simple form of periodicity boundary condition was used for the channel flow with transverse rectangular ribs. The effects of various parameters on heat transfer in channel flow with transverse ribs and in impingement cooling were investigated using the Yap modified Jones and Launder low Reynolds number k-epsilon turbulence model. For the channel flow, predictions were in adequate agreement with experiment for constant property flow, with the results for friction superior to those for heat transfer. For impingement cooling, the agreement with experiment was generally good, but the results suggest that improved modelling of the dissipation rate of turbulence kinetic energy is required in order to obtain improved heat transfer prediction, especially near the stagnation point. The k-epsilon turbulence model was used to predict the mean flow and heat transfer for constant and variable property flows. The effect of variable properties for channel flow was investigated using the same turbulence model, but comparison with experiment yielded no clear conclusions. Also, the wall function method was modified for use in the variable properties flow with a non-adiabatic surface, and an empirical model is suggested to correctly account for the behavior of the viscous sublayer with heating.

Chang, Byong Hoon↗

A New Axisymmetric MHD Model of the Interaction of the Solar Wind with Venus

A new two-dimensional axisymmetric MHD model is used to study the interaction of the solar wind with Venus under conditions where the interplanetary field is approximately aligned with the solar wind velocity. This numerical model solves the MHD transport equations for density, velocity, pressure, and magnetic field on an adaptively refined, unstructured grid system. This use of an adaptive grid allows high spatial resolution in regions of large density/velocity gradients and yet can be run on a workstation. The actual grid sizes vary from about 0.06 R(sub v) near the bowshock to 2 R(sub v) in the unperturbed solar wind. The results of the calculations are compared with observed magnetic field values obtained from the magnetometer on the Pioneer Venus Orbiter, at a time when the angle between the solar wind velocity vector and the interplanetary magnetic field (IMF) was only 7.6 deg. Good qualitative agreement between the observed and calculated field behavior is found. The overall results suggest that the induced magnetotail disappears when the IMF is radial for an extended time period and implies that it weakens when the field rotated through a near-radial orientation.

DeZeeuw, Darren L.↗

Development of the Joint NASA/NCAR General Circulation Model

The Data Assimilation Office at NASA/Goddard Space Flight Center is collaborating with NCAR/CGD in an ambitious proposal for the development of a unified climate, numerical weather prediction, and chemistry transport model which is suitable for global data assimilation of the physical and chemical state of the Earth's atmosphere. A prototype model based on the NCAR CCM3 physics and the NASA finite-volume dynamical core has been built. A unique feature of the NASA finite-volume dynamical core is its advanced tracer transport algorithm on the floating Lagrangian control-volume coordinate. The model currently has a highly idealized ozone production/loss chemistry derived from the observed 2D (latitude-height) climatology of the recent decades. Nevertheless, the simulated horizontal wave structure of the total ozone is in good qualitative agreement with the observed (TOMS). Long term climate simulations and NWP experiments have been carried out. Current up to date status and futur! e plan will be discussed in the conference.

Lin, S.-J.↗

A Sectored-One-Dimensional Model for Simulating Combustion Instabilities in Premix Combustors

A one-dimensional, CFD based combustor simulation has been developed that exhibits self-excited, thermoacoustic oscillations in premixed combustor geometries that typically have large, abrupt changes in cross sectional area. The combustor geometry is approximated by dividing it into a finite number of one-dimensional sectors. Within each sector, the equations of motion are integrated numerically, along with a species transport and a reaction equation. Across the sectors, mass and energy are conserved, and momentum loss is prescribed using appropriately compatible boundary conditions that account for the area change. The resulting simulation and associated boundary conditions essentially represent a one-dimensional, multi-block technique. Details of the simulation code are presented herein. Results are then shown comparing experimentally observed and simulated operation of a particular combustor rig that exhibited different instabilities at different operating points. It will be shown that the simulation closely matched the rig data in oscillation amplitudes, frequencies, and operating points at which the instabilities occurred. Finally, advantages and limitations of the simulation technique are discussed.

Paxson, Daniel E.↗