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 1,279 records · Page 71

jaxhps: An elliptic PDE solver built with machine learning in mind

Elliptic partial differential equations (PDEs) can model many physical phenomena, such as electrostatics, acoustics, wave propagation, and diffusion. In scientific machine learning settings, a high-throughput PDE solver may be required to generate a training dataset, run in the inner loop of an iterative algorithm, or interface directly with a deep neural network. To provide value to machine learning users, such a PDE solver must be compatible with standard automatic differentiation frameworks, scale efficiently when run on graphics processing units (GPUs), and maintain high accuracy for a large range of input parameters. We have designed the jaxhps package with these use-cases in mind by implementing a highly efficient and accurate solver for elliptic problems with native hardware acceleration and automatic differentiation support.

97 MATHEMATICS AND COMPUTING↗

Analysis of Tempered Fractional Operators

Tempered fractional operators provide an improved predictive capability for modeling anomalous effects that cannot be captured by standard partial differential equations. These effects include subdiffusion and superdiffusion (i.e. the mean square displacement in a diffusion process is proportional to a fractional power of the time), that often occur in, e.g., geoscience and hydrology. We analyze tempered fractional operators within the nonlocal vector calculus framework in order to assimilate them to the rigorous mathematical structure developed for nonlocal models. First, we show they are special instances of generalized nonlocal operators by means of a proper choice of the nonlocal kernel. Then, we present a plan for showing tempered fractional operators are equivalent to truncated fractional operators. These truncated operators are useful because they are less computationally intensive than the tempered operators.

97 MATHEMATICS AND COMPUTING↗

Test-particle method in kinetic theory of a plasma.

The introduction of a test particle into a system is considered. The system may be described by the Born-Bogoliubov-Green-Kirkwood-Yvon hierarchy. The field particles form a cloud which surrounds the test particle. The cloud is described by a conditional probability function which satisfies a certain equation. A generalization of the superposition principle reported by Rostoker (1964) to higher order correlation functions is discussed. Kinetic equations with the generalized Lenard-Balescu term are obtained, taking into account also diffusion by waves. The characteristics regarding the absorption or emission of waves by particles can be calculated.

Matsuda, K.↗

Design and calibration of a rocket-borne electron spectrometer for investigation of particle ionization in the nighttime midlatitude E region

An explanation was developed for the formation, near midnight at midlatitudes, of a broad electron density layer extending approximately from 120 to 180 km and usually referred to as the intermediate E layer. The responsible mechanism is believed to be the converging vertical ion drifts resulting from winds of the solar semidiurnal tide. Numerical solutions of the continuity equation appropriate to the intermediate layer is described for particular models of ion drift, diffusion coefficents, and ionization production. Analysis of rocket observations of the layer show that the ionization rate is highly correlated with the planetary geomagnetic index, K sub p. Particle flux measurements support the idea that energetic electrons are the principal source of this ionization. A semiconductor spectrometer experiment for investigation of the particle flux, spectrum, and angular properties was designed and successfully flown on a Nike Apache rocket. A detailed description of the theory, design, and calibration of the experiment and some preliminary results presented.

Voss, H. D.↗

Two-stream theory of spectral reflectance of snow

Spectral reflectance of snow under diffuse illumination is studied using the two-stream approximation of the radiative transfer equation. The scattering and absorption within the snowcover due to the randomly distributed ice grains are characterized by the single scattering albedo and anisotropic phase function. Geometric optics calculations are used to relate the scattering and absorption parameters to grain size and density of snow. Analytical expressions for the intensity within the snowpack and the asymptotic flux extinction coefficient are also obtained. Good agreement is shown between the theory and available experimental data on visible and near-infrared reflectance and asymptotic flux extinction coefficient. The theory also may be used to explain the observed effect of aging on the snow reflectance.

Choudhury, B. J.↗

Two-stream theory of reflectance of snow

Spectral reflectance of snow under diffuse illumination is studied using the two-stream approximation of the radiative transfer equation. The scattering and absorption parameters of the radiative transfer equation - the single scattering albedo, the optical depth, and the integrated phase function are obtained from the grain size and density of snow. Analytical expressions for the intensity within the snowpack, the reflectance, and the asymptotic flux extinction coefficient, are given. Good agreement is shown between the theory and available experimental data on visible and near-infrared reflectance, and the asymptotic flux extinction coefficient. The theory may also be used to explain the observed effect of aging on the snow reflectance.

Choudhury, B. J.↗

Viscous, resistive MHD stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one dimensional magnetohydrodynamic (MHD) quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds like numbers. Marginal stability curves, growth rates versus Reynolds like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result which appears general is that instability was found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗

Simulations of magnetic-flux transport in solar active regions

The evolution of several observed solar active regions is simulated by solving a transport equation for magnetic flux at the photosphere. The rates of rotation, meridional flow, and diffusion of the flux are determined self-consistently in the calculations. The findings are in good quantitative agreement with previous measures of the rotation rate and diffusion constant associated with photospheric magnetic fields. Although the meridional velocities are consistent in direction and magnitude with recently reported poleward flows, relatively large uncertainties in the velocity determinations make this result inconclusive.

Devore, C. R.↗

A method for calculating turbulent boundary layers and losses in the flow channels of turbomachines

An interactive inviscid core flow-boundary layer method is presented for the calculation of turbomachine channel flows. For this method, a one-dimensional inviscid core flow is assumed. The end-wall and blade surface boundary layers are calculated using an integral entrainment method. The boundary layers are assumed to be collateral and thus are two-dimensional. The boundary layer equations are written in a streamline coordinate system. The streamwise velocity profiles are approximated by power law profiles. Compressibility is accounted for in the streamwise direction but not in the normal direction. Equations are derived for the special cases of conical and two-dimensional rectangular diffusers. For these cases, the assumptions of a one-dimensional core flow and collateral boundary layers are valid. Results using the method are compared with experiment and good quantitative agreement is obtained.

Schumann, Lawrence F.↗

A method for calculating turbulent boundary layers and losses in the flow channels of turbomachines

An interactive inviscid core flow-boundary layer method is presented for the calculation of turbomachine channel flows. For this method, a one-dimensional inviscid core flow is assumed. The end-wall and blade surface boundary layers are calculated using an integral entrainment method. The boundary layers are assumed to be collateral and thus are two-dimensional. The boundary layer equations are written in a streamline coordinate system. The streamwise velocity profiles are approximated by power law profiles. Compressibility is accounted for in the streamwise direction but not in the normal direction. Equations are derived for the special cases of conical and two-dimensional rectangular diffusers. For these cases, the assumptions of a one-dimensional core flow and collateral boundary layers are valid. Results using the method are compared with experiment and good quantitative agreement is obtained.

Schumann, Lawrence F.↗

The computation of non-equilibrium chemically-reacting flows

Numerical techniques for the analysis of flows with nonequilibrium thermodynamics and chemistry are developed and demonstrated. Particular attention is given to the formulation of the thermodynamic and chemical gas models, the derivation of the governing equations, and the solution procedures. Results for typical problems involving shock tubes and supersonic diffusers are presented in extensive graphs and briefly characterized.

Grossman, B.↗

Cosmic-ray streaming in the Born approximation

The present work invokes the Born approximation to derive a more accurate expression for the streaming of cosmic rays parallel to the mean magnetic field. While all prior results pertaining to the helicity dependence of the diffusion coefficient and convection speed can be recovered as special cases from this streaming equation, it is concluded that a new set of transport parameters presented here is more appropriate for the solar modulation of galactic cosmic rays. In addition, a new parameter related to time variability, which may be a dominant cause of charge sign-dependent transport of solar particles, is introduced.

Bieber, J. W.↗

A time-accurate implicit method for chemically reacting flows at all Mach numbers

The objective of this work is to develop a unified solution algorithm capable of treating time-accurate chemically reacting flows at all Mach numbers, ranging from molecular diffusion velocities to supersonic speeds. A rescaled pressure term is used in the momentum equation to circumvent the singular behavior of pressure at low Mach numbers. A dual time-stepping integration procedure is established. The system eigenvalues become well behaved and have the same order of magnitude, even in the very low Mach number regime. The computational efficiency for moderate and high speed flow is competitive with the conventional density-based scheme. The capabilities of the algorithm are demonstrated by applying it to selected model problems including nozzle flows and flame dynamics.

Withington, J. P.↗

Simulation of the present-day atmospheric ozone, odd nitrogen, chlorine and other species using a coupled 2-D model in isentropic coordinates

A zonally averaged two-dimensional model for chemical species in the atmosphere (Olaguer et al., 1990), with coupled dynamics, radiation transfer, and chemistry, is used to study the present-day atmospheric ozone, odd nitrogen, chlorine, other species. The model utilizes all zonally averaged physical equations of momentum, energy, and mass and determines self-consistently the advective and diffusive tranasport parameters from the temperature specific to the period of observation. A comparison of the model results with observations showed good agreement.

Yang, H.↗

The implementation and validation of improved landsurface hydrology in an atmospheric general circulation model

Landsurface hydrological parameterizations are implemented in the NASA Goddard Institute for Space Studies (GISS) General Circulation Model (GCM). These parameterizations are: (1) runoff and evapotranspiration functions that include the effects of subgrid scale spatial variability and use physically based equations of hydrologic flux at the soil surface, and (2) a realistic soil moisture diffusion scheme for the movement of water in the soil column. A one dimensional climate model with a complete hydrologic cycle is used to screen the basic sensitivities of the hydrological parameterizations before implementation into the full three dimensional GCM. Results of the final simulation with the GISS GCM and the new landsurface hydrology indicate that the runoff rate, especially in the tropics is significantly improved. As a result, the remaining components of the heat and moisture balance show comparable improvements when compared to observations. The validation of model results is carried from the large global (ocean and landsurface) scale, to the zonal, continental, and finally the finer river basin scales.

Johnson, Kevin D.↗

Mixing and mixing enhancement in supersonic reacting flowfields

Propulsion systems planned for use late in this century and beyond will require appropriate physical models for describing supersonic combustion and numerical techniques for solving the model governing equations. A computer program to study these flows is reported which considers the multicomponent diffusion and convection of important chemical species, the finite-rate reaction of these species, and the resulting interaction of the field mechanics and the chemistry. The application of the program to a spatially developing and reacting mixing layer, which serves an an excellent physical model for the mixing and reaction processes that take place in a scramjet combustor, is reported. Several techniques to enhance the fuel-air mixing and growth of that layer and improve its overall combustion efficiency are considered.

Drummond, J. P.↗

The implementation and validation of improved land-surface hydrology in an atmospheric general circulation model

New land-surface hydrologic parameterizations are implemented into the NASA Goddard Institute for Space Studies (GISS) General Circulation Model (GCM). These parameterizations are: 1) runoff and evapotranspiration functions that include the effects of subgrid-scale spatial variability and use physically based equations of hydrologic flux at the soil surface and 2) a realistic soil moisture diffusion scheme for the movement of water and root sink in the soil column. A one-dimensional climate model with a complete hydrologic cycle is used to screen the basic sensitivities of the hydrological parameterizations before implementation into the full three-dimensional GCM. Results of the final simulation with the GISS GCM and the new land-surface hydrology indicate that the runoff rate, especially in the tropics, is significantly improved. As a result, the remaining components of the heat and moisture balance show similar improvements when compared to observations. The validation of model results is carried from the large global (ocean and land-surface) scale to the zonal, continental, and finally the regional river basin scales.

Johnson, Kevin D.↗

Forced Convection and Sedimentation Past a Flat Plate

The steady laminar flow of a well-mixed suspension of monodisperse solid spheres, convected steadily past a horizontal flat plate and sedimenting under the action of gravity, is examined. It is shown that, in the limit as Re approaches infinity and epsilon approaches 0, where Re is the bulk Reynolds number and epsilon is the ratio of the particle radius a to the characteristic length scale L, the analysis for determining the particle concentration profile has several aspects in common with that of obtaining the temperature profile in forced-convection heat transfer from a wall to a fluid stream moving at high Reynolds and Prandtl numbers. Specifically, it is found that the particle concentration remains uniform throughout the O(Re(exp -1/2)) thick Blasius boundary layer except for two O(epsilon(exp 2/3)) thin regions on either side of the plate, where the concentration profile becomes non-uniform owing to the presence of shear-induced particle diffusion which balances the particle flux due to convection and sedimentation. The system of equations within this concentration boundary layer admits a similarity solution near the leading edge of the plate, according to which the particle concentration along the top surface of the plate increases from its value in the free stream by an amount proportional to X(exp 5/6), with X measuring the distance along the plate, and decreases in a similar fashion along the underside. But, unlike the case of gravity settling on an inclined plate in the absence of a bulk flow at infinity considered earlier, here the concentration profile remains continuous everywhere. For values of X beyond the region near the leading edge, the particle concentration profile is obtained through the numerical solution of the relevant equations. It is found that, as predicted from the similarity solution, there exists a value of X at which the particle concentration along the top side of the plate attains its maximum value phi(sub m) and that, beyond this point, a stagnant sediment layer will form that grows steadily in time. This critical value of X is computed as a function of phi(sub s), the particle volume fraction in the free stream. In contrast, but again in conformity with the similarity solution, for values of X sufficiently far removed from the leading edge along the underside of the plate, a particle-free region is predicted to form adjacent to the plate. This model, with minor modifications, can be used to describe particle migration in other shear flows, as, for example, in the case of crossflow microfiltration.

Pelekasis, Nikolaos A.↗