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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Plasma diffusion at the magnetopause? The case of lower hybrid drift waves

The diffusion expected from the quasilinear theory of the lower hybrid drift instability at the Earth's magnetopause is recalculated. The resulting diffusion coefficient is in principle just marginally large enough to explain the thickness of the boundary layer under quiet conditions, based on observational upper limits for the wave intensities. Thus, one possible model for the boundary layer could involve equilibrium between the diffusion arising from lower hybrid waves and various low processes. However, some recent data and simulations seems to indicate that the magnetopause is not consistent with such a soft diffusive equilibrium model. Furthermore, investigation of the nonlinear equations for the lower hybrid waves for magnetopause parameters indicates that the quasilinear state may never arise because coalescence to large wavelengths, followed by collapse once a critical wavelengths is reached, occur on a time scale faster than the quasilinear diffusion. In this case, an inhomogeneous boundary layer is to be expected. More simulations are required over longer time periods to explore whether this nonlinear evolution really takes place at the magnetopause.

Treumann, R. A.↗
NASA NTRS1983

Ion cyclotron waves in the Io plasma torus - Polarization reversal of whistler mode noise

Because of the presence of multiple ion species in the Io plasma torus, whistler mode noise can be converted to ion cyclotron waves via a polarization reversal process at the local crossover frequency. Using whistler mode intensity measurements in the Jovian magnetosphere from Voyager 1, the pitch-angle diffusion rates that would occur if the noise is converted to ion cyclotron waves is estimated. Typical pitch-angle diffusion coefficients range from 0.000001/sec for protons resonating near the equator to 0.0001/sec for 10-keV O(+) ions resonating at high latitudes. Although complete bounce averaged diffusion coefficients have not yet been computed, preliminary estimates indicate that the energetic ion precipitation caused by these waves may be able to account for the EUV auroral emissions at the foot of the torus field lines.

Gurnett, D. A.↗
NASA NTRS1984

Turbo fluid machinery and diffusers

The general theory behind turbo devices and diffusers is explained. Problems and the state of research on basic equations of flow and experimental and measuring methods are discussed. Conventional centrifugation-type compressor and fan diffusers are considered in detail.

Sakurai, T.↗
NASA NTRS1978

Partially averaged field approach to cosmic ray diffusion

A new nonlinear technique is used to derive the kinetic equation for particles interacting with turbulent fluctuations. Difficulties associated with quasi-linear theory are avoided. The new method evaluates the effects of the fluctuations along particle orbits which themselves include the effects of a statistically averaged subset of the possible configuration of the turbulence. As an illustration, the pitch-angle diffusion coefficient is calculated for particles interacting with 'slab model' magnetic turbulence, i.e., magnetic fluctuations linearly polarized transverse to a mean magnetic field. The pitch-angle diffusion coefficient is determined in the vicinity of 90-deg pitch angles where quasi-linear theory breaks down. The spatial diffusion coefficient parallel to a mean magnetic field is evaluated by use of the calculated pitch-angle diffusion coefficient. It is suggested that the partially averaged field method is not limited to small amplitude fluctuating fields and hence is not a perturbation theory.

Jones, F. C.↗
NASA NTRS1993

Alveolar ventilation to perfusion heterogeneity and diffusion impairment in a mathematical model of gas exchange

This study describes a two-compartment model of pulmonary gas exchange in which alveolar ventilation to perfusion (VA/Q) heterogeneity and impairment of pulmonary diffusing capacity (D) are simultaneously taken into account. The mathematical model uses as input data measurements usually obtained in the lung function laboratory. It consists of two compartments and an anatomical shunt. Each compartment receives fractions of alveolar ventilation and blood flow. Mass balance equations and integration of Fick's law of diffusion are used to compute alveolar and blood O2 and CO2 values compatible with input O2 uptake and CO2 elimination. Two applications are presented. The first is a method to partition O2 and CO2 alveolar-arterial gradients into VA/Q and D components. The technique is evaluated in data of patients with chronic obstructive pulmonary disease (COPD). The second is a theoretical analysis of the effects of blood flow variation in alveolar and blood O2 partial pressures. The results show the importance of simultaneous consideration of D to estimate VA/Q heterogeneity in patients with diffusion impairment. This factor plays an increasing role in gas alveolar-arterial gradients as severity of COPD increases. Association of VA/Q heterogeneity and D may produce an increase of O2 arterial pressure with decreasing QT which would not be observed if only D were considered. We conclude that the presented computer model is a useful tool for description and interpretation of data from COPD patients and for performing theoretical analysis of variables involved in the gas exchange process.

NASA Discipline Environmental Health↗
NASA NTRS1986

Diffuser investigation for advanced multi-megawatt MPD engine development

Steady-state testing of advanced, multimegawatt MPD engines requires vacuum facilities with pumping speeds which are not presently available. Gas dynamic diffusers have been proposed as a possible solution to this problem. An analytical investigation into the feasibility of using a diffuser for pumping the MPD engine exhaust is presented here. This analysis uses a detailed equation of state based on the partition function for the high temperature argon exhaust. On the basis of the electron-ion momentum exchange collision frequency in the exhaust plasma, it was concluded that the diffuser gas dynamics could be modeled to a first approximation with ordinary continuum equations. Calculations of the stagnation pressure in the diffuser, downstream of a strong normal shock, yielded pressures on the order of 10 torr suggesting that the diffuser is feasible for pumping the MPD engine exhaust.

Brophy, J. R.↗
NASA NTRS1981

Factors governing the ratio of inward to outward diffusing flux of satellite ions

The sensitivity of the flux ratio to the uncertainties in the boundary conditions and diffusion coefficient is analyzed. The principal result is that if the diffusion is driven externally, for example, by winds in Jupiter's ionosphere, the inward flux exceeds or is of the same order of magnitude as the outward flux. On the other hand, if the diffusion is driven locally by the centrifugal interchange instability, as suggested by one interpretation of Voyager plasma data, the outward flux can be one or two orders of magnitude greater than the inward flux. Fits to the Voyager in situ plasma torus data obtained through solutions to the transport equations appropriate to wind driven and centrifugally driven diffusion are compared.

Richardson, J. D.↗
NASA NTRS1994

Modification of the Two-equation Turbulence Model in NPARC to a Chien Low Reynolds Number K-epsilon Formulation

This report documents the changes that were made to the two-equation k-epsilon turbulence model in the NPARC (National-PARC) code. The previous model based on the low Reynolds number model of Speziale, was replaced with the low Reynolds number k-epsilon model of Chien. The most significant difference was in the turbulent Prandtl numbers appearing in the diffusion terms of the k and epsilon transport equations. A new inflow boundary condition and stability enhancements were also implemented into the turbulence model within NPARC. The report provides the rationale for making the change to the Chien model, code modifications required, and comparisons of the performances of the new model with the previous k-epsilon model and algebraic models used most often in PARC/NPARC. The comparisons show that the Chien k-epsilon model installed here improves the capability of NPARC to calculate turbulent flows.

Georgiadis, Nicholas J.↗
NASA NTRS1998

Comments on the Diffusive Behavior of Two Upwind Schemes

The diffusive characteristics of two upwind schemes, multi-dimensional fluctuation splitting and locally one-dimensional finite volume, are compared for scalar advection-diffusion problems. Algorithms for the two schemes are developed for node-based data representation on median-dual meshes associated with unstructured triangulations in two spatial dimensions. Four model equations are considered: linear advection, non-linear advection, diffusion, and advection-diffusion. Modular coding is employed to isolate the effects of the two approaches for upwind flux evaluation, allowing for head-to-head accuracy and efficiency comparisons. Both the stability of compressive limiters and the amount of artificial diffusion generated by the schemes is found to be grid-orientation dependent, with the fluctuation splitting scheme producing less artificial diffusion than the finite volume scheme. Convergence rates are compared for the combined advection-diffusion problem, with a speedup of 2.5 seen for fluctuation splitting versus finite volume when solved on the same mesh. However, accurate solutions to problems with small diffusion coefficients can be achieved on coarser meshes using fluctuation splitting rather than finite volume, so that when comparing convergence rates to reach a given accuracy, fluctuation splitting shows a speedup of 29 over finite volume.

Wood, William A.↗
NASA NTRS2003

Transport Equation Based Wall Distance Computations Aimed at Flows With Time-Dependent Geometry

Eikonal, Hamilton-Jacobi and Poisson equations can be used for economical nearest wall distance computation and modification. Economical computations may be especially useful for aeroelastic and adaptive grid problems for which the grid deforms, and the nearest wall distance needs to be repeatedly computed. Modifications are directed at remedying turbulence model defects. For complex grid structures, implementation of the Eikonal and Hamilton-Jacobi approaches is not straightforward. This prohibits their use in industrial CFD solvers. However, both the Eikonal and Hamilton-Jacobi equations can be written in advection and advection-diffusion forms, respectively. These, like the Poisson s Laplacian, are commonly occurring industrial CFD solver elements. Use of the NASA CFL3D code to solve the Eikonal and Hamilton-Jacobi equations in advective-based forms is explored. The advection-based distance equations are found to have robust convergence. Geometries studied include single and two element airfoils, wing body and double delta configurations along with a complex electronics system. It is shown that for Eikonal accuracy, upwind metric differences are required. The Poisson approach is found effective and, since it does not require offset metric evaluations, easiest to implement. The sensitivity of flow solutions to wall distance assumptions is explored. Generally, results are not greatly affected by wall distance traits.

Tucker, Paul G.↗
NASA NTRS2003

Transport Equation Based Wall Distance Computations Aimed at Flows With Time-Dependent Geometry

Eikonal, Hamilton-Jacobi and Poisson equations can be used for economical nearest wall distance computation and modification. Economical computations may be especially useful for aeroelastic and adaptive grid problems for which the grid deforms, and the nearest wall distance needs to be repeatedly computed. Modifications are directed at remedying turbulence model defects. For complex grid structures, implementation of the Eikonal and Hamilton-Jacobi approaches is not straightforward. This prohibits their use in industrial CFD solvers. However, both the Eikonal and Hamilton-Jacobi equations can be written in advection and advection-diffusion forms, respectively. These, like the Poisson's Laplacian, are commonly occurring industrial CFD solver elements. Use of the NASA CFL3D code to solve the Eikonal and Hamilton-Jacobi equations in advective-based forms is explored. The advection-based distance equations are found to have robust convergence. Geometries studied include single and two element airfoils, wing body and double delta configurations along with a complex electronics system. It is shown that for Eikonal accuracy, upwind metric differences are required. The Poisson approach is found effective and, since it does not require offset metric evaluations, easiest to implement. The sensitivity of flow solutions to wall distance assumptions is explored. Generally, results are not greatly affected by wall distance traits.

Tucker, Paul G.↗
NASA NTRS1998

Diffusion Characteristics of Upwind Schemes on Unstructured Triangulations

The diffusive characteristics of two upwind schemes, multi-dimensional fluctuation splitting and dimensionally-split finite volume, are compared for scalar advection-diffusion problems. Algorithms for the two schemes are developed for node-based data representation on median-dual meshes associated with unstructured triangulations in two spatial dimensions. Four model equations are considered: linear advection, non-linear advection, diffusion, and advection-diffusion. Modular coding is employed to isolate the effects of the two approaches for upwind flux evaluation, allowing for head-to-head accuracy and efficiency comparisons. Both the stability of compressive limiters and the amount of artificial diffusion generated by the schemes is found to be grid-orientation dependent, with the fluctuation splitting scheme producing less artificial diffusion than the dimensionally-split finite volume scheme. Convergence rates are compared for the combined advection-diffusion problem, with a speedup of 2-3 seen for fluctuation splitting versus finite volume when solved on the same mesh. However, accurate solutions to problems with small diffusion coefficients can be achieved on coarser meshes using fluctuation splitting rather than finite volume, so that when comparing convergence rates to reach a given accuracy, fluctuation splitting shows a 20-25 speedup over finite volume.

Wood, William A.↗
NASA NTRS2023

A Data-driven, Physics-based Transport Model of Solar Energetic Particles Accelerated by Coronal Mass Ejection Shocks Propagating through the Solar Coronal and Heliospheric Magnetic Fields

In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.

Solar energetic particles↗
NASA NTRS1989

Hopf bifurcation in the driven cavity

Incompressible two dimensional calculations are reported for the impulsively started lid driven cavity with aspect ratio two. The algorithm is based on the time dependent streamfunction equation, with a Crank-Nicolson differencing scheme for the diffusion terms, and with an Adams-Bashforth scheme for the convection terms. A multigrid method is used to solve the linear implicit equations at each time step. Periodic asymptotic solutions have been found for Re = 10000 and for Re = 5000. The Re = 5000 results are validated by grid refinement calculations. The solutions are shown to be precisely periodic, and care is taken to demonstrate that asymptotic states were reached. A discussion is included about the indicators that are used to show that an asymptotic state was reached, and to show that the asymptotic state is indeed periodic.

Goodrich, John W.↗