SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Problems in simulating the stratocumulus-topped boundary layer with a third-order closure model

The Andre et al. (1976, 1978) third-order closure model, in which the time rate of change terms, the relaxation and rapid effects for pressure-related terms, and the clipping approximation are used along with the quasi-normal closure, is invoked in the study of turbulence in a cloudy layer that is radiatively cooled from above. A spurious oscillation whose greatest amplitude lies near the inversion is shown by analysis to arise from the mean gradient and buoyancy terms of the triple-moment equations. An attempt is made to damp the oscillation through the introduction of diffusion terms into the triple-moment equations. The results obtained are noted to be sensitive to the ad hoc eddy coefficient applied in the third-moment equations.

Moeng, C.-H.↗

NASA NTRS1979

The prediction of two- and three-dimensional asymmetrical turbulent wakes - A comparision of the performance of three turbulence models for the effects of streamline curvature and rotation

Two- and three-dimensional turbulent wakes which develop under the influence of streamline curvature and rotation were calculated with three turbulence closure models. The first model comprised of transport equations for the turbulent kinetic energy and the rate of energy dissipation. The second model comprised of equations for the rate of turbulent kinetic energy dissipation and Reynolds stresses, but the effects of the convection and diffusion in the Reynolds stress transport equation were handled collectively. The third model utilizes equations of turbulent kinetic energy dissipation and Reynolds stresses in nearly exact form. All of the three models were modified for the effects of streamline curvature. The second and third models include the effects of rotation through the rotation-originated redistribution term in the transport equation of Reynolds stresses. The numerical results for the wakes of airfoils, cascades and turbomachinery rotor blades demonstrate that the second and third models provide accurate predictions, but computer time and storage can be considerably saved with the second model.

Hah, C.↗

NASA NTRS1992

A finite element method for the thermochemical decomposition of polymeric materials. I - Theory

The governing differential equations are developed to model the thermomechanical behavior of chemically decomposing, polymeric materials. These equations account for thermal and gaseous diffusion through a poroelastic, transversely isotropic solid. The Bubnov-Galerkin finite element method is applied to the governing equations to cast the coupled set into a single matrix equation. A method for solving these equations simultaneously at each time step is discussed.

Sullivan, R. M.↗

NASA NTRS1980

Velocity and concentration distribution in a Stefan diffusion tube

Numerical solutions for the diffusion- and gravity-driven flow in a cylindrical Stefan tube were obtained from the coupled diffusion and Navier-Stokes equations of Peclet numbers 0.3, 1 and 5. Distributions of binary component concentrations and velocities were calculated. The mass average velocity is parabolic in nature, except at high Peclet numbers. The solvent is not stagnant but recirculates, even in the absence of gravity. Radial concentration gradients develop which act convectively destabilizing. Consequences for the deduction of diffusion coefficients from Stefan tube experiments are discussed.

Markham, B. L.↗

DOE OSTI2013

ENERGETIC PARTICLE CROSS-FIELD PROPAGATION EARLY IN A SOLAR EVENT

Solar energetic particles (SEPs) have been observed to easily spread across heliographic longitudes, and the mechanisms responsible for this behavior remain unclear. We use full-orbit simulations of a 10 MeV proton beam in a turbulent magnetic field to study to what extent the spread across the mean field can be described as diffusion early in a particle event. We compare the full-orbit code results to solutions of a Fokker-Planck equation including spatial and pitch angle diffusion, and of one including also propagation of the particles along random-walking magnetic field lines. We find that propagation of the particles along meandering field lines is the key process determining their cross-field spread at 1 AU at the beginning of the simulated event. The mean square displacement of the particles an hour after injection is an order of magnitude larger than that given by the diffusion model, indicating that models employing spatial cross-field diffusion cannot be used to describe early evolution of an SEP event. On the other hand, the diffusion of the particles from their initial field lines is negligible during the first 5 hr, which is consistent with the observations of SEP intensity dropouts. We conclude that modeling SEP events must take into account the particle propagation along meandering field lines for the first 20 hr of the event.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1982

Solar reflection from a regular array of horizontally finite clouds

The reflected flux from a regular array of 2- and 3-D clouds has been computed to estimate the effect of fractional cloud cover on albedos and the solar flux available to heat the earth's surface. The broken clouds are represented by a regular array of identical cuboids for the 3-D problem and equally spaced, infinitely long, bars for the 2-D problem. A diffusion approximation to the radiative transfer equation is used to compute the fluxes leaving each face of the cloud. Interaction between clouds is simulated by assuming diffuse exitance from the cloud faces and applying angle factors to obtain modified boundary conditions on each cloud face.

Weinman, J. A.↗

NASA NTRS1986

Surface mass transport of alumina

Measurements were made of the kinetics of thermal grooving at the intersection of rhombohedral twin boundaries with the (1 0 -1 0) plane in alumina, at temperatures from 1773 to 2273 K. An analysis of the values obtained on the basis of the model of Mullins (1957) indicated that surface diffusion was the dominant mass transport mechanism. Comparisons are made with the results of other investigations of alumina, and an equation is derived for the surface diffusion value.

Gaddipati, A. R.↗

DOE OSTI2022

On the Diffusivity of Moist Static Energy and Implications for the Polar Amplification Response to Climate Warming

Energy balance models (EBMs) have been widely used in a range of climate problems, but the assumption of constant diffusivity in the parameterization of the moist static energy (MSE) flux can be hardly justified. We demonstrate in this study that the diffusive MSE flux can be derived from the basic energy balance equation with a few tolerable assumptions. The estimated diffusivity is both spatially and seasonally dependent, and its midlatitude average is then tested against several scaling theories for the midlatitude eddy diffusivity. The result supports the diffusivity theory of Held and Larichev (1996) modified for the moist atmosphere, affording a dynamics-based parameterization of MSE diffusivity. The implementation of the parameterization in an EBM leads to an interactive MSE diffusivity that accounts for the midlatitude eddy response to climate forcing perturbations. Under a uniform radiative forcing, the EBM with a diffusivity so parameterized produces a weakening of the midlatitude diffusivity and a modestly polar-amplified surface temperature response as an inevitable outcome under the dual constraints of the nonlinear Clausius-Clapeyron relation and the temperature gradient-dependent diffusivity, even in the absence of any poleward amplifying radiative feedbacks. As the consequence of more isothermal temperature and reduced diffusivity, the variance of the midlatitude surface temperature also decreases with warming.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1981

Numerical simulation of reentry flow around the Shuttle Orbiter including real gas effects

The flight environment of a reentry vehicle is predicted from the numerical solution of fluid dynamics equations for the region between the blunt nose and the swept wings. The inviscid portion of the shock layer is modeled by the Euler equations, but the laminar viscous flow adjacent to the wall is modeled by the approximate parabolic Navier-Stokes equations. The approximations made to the axial gradients of pressure and diffusive fluxes enable the coupled inviscid and viscous equations to be solved efficiently along the body axis. The equilibrium air aftbody code contains significant improvements over its predecessor which only considers a perfect gas model and noncircular configurations. The inclusion of numerical damping either explicitly or implicitly has extended its capabilities for predicting flow field around a winged configuration at higher Machs and greater angles of attack. The results are obtained on the cylindrical coordinates and satisfactory for the Shuttle Orbiter at a free-stream Mach number of 22 and an angle of attack of 40 deg. Also discussed are the inviscid formulation and its application for wind-tunnel conditions.

Li, C. P.↗

DOE OSTI2015

Anomalous law of cooling

We analyze the temperature relaxation phenomena of systems in contact with a thermal reservoir that undergoes a non-Markovian diffusion process. From a generalized Langevin equation, we show that the temperature is governed by a law of cooling of the Newton’s law type in which the relaxation time depends on the velocity autocorrelation and is then characterized by the memory function. The analysis of the temperature decay reveals the existence of an anomalous cooling in which the temperature may oscillate. Despite this anomalous behavior, we show that the variation of entropy remains always positive in accordance with the second law of thermodynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Microscopic Theory of Long-Time Center-of-Mass Self-Diffusion and Anomalous Transport in Ring Polymer Liquids

We construct a microscopic theory at the level of segment-scale correlated space–time intermolecular forces for the long-time center-of-mass (CM) diffusion constant and intermediate-time non-Fickian transport in dense solutions and melts of ring polymers. The approach combines ideas of polymer, colloid, and liquid-state statistical mechanics to quantify how the multifractal intra-ring conformational structure and inter-ring packing correlations determine dynamic caging constraints and time-dependent friction. Breakdown of Rouse theory is predicted to occur due to length scale-dependent temporal correlation of forces exerted on pairs of tagged ring segments from surrounding polymers. At large enough degrees of polymerization (N), a stronger scaling of the CM diffusion constant (D ∝ N –2 ) is predicted, with a crossover N D proportional to the product of the system-specific macromolecular volume fraction and dimensionless compressibility. In analogy with caging effects in glass-forming fluids, the theory does appear to begin to fail at sufficiently high N/N D for the center-of-mass diffusivity, likely due to another crossover to an even slower activated transport regime. However, use of N/N D with the theoretically predicted N D , in conjunction with dynamic blob scaling ideas for local Rouse friction, collapses semidilute and concentrated solution simulation data onto a master curve. Based on the same physical ideas employed to predict the diffusion constant, a generalized Langevin equation description is formulated for intermediate-time CM transport. It predicts two subdiffusive regimes that emerge due to the self-similar nature of internal ring structure. Analytic and numerical predictions for the apparent non-Fickian exponent as a function of time, N/N D and dimensionless compressibility, the evolution of the maximum degree of subdiffusive motion with N/N D , and the time scale for recovering Fickian diffusion are made, all of which are in good accordance with melt simulations. Furthermore, the present work sets the stage to address activated dynamics and glass formation on the macromolecular scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Equation of state and transport properties of warm dense helium via quantum molecular dynamics simulations

The equation of state, self-diffusion, and viscosity coefficients of helium have been investigated by quantum molecular dynamics (QMD) simulations in the warm dense matter regime. Our simulations are validated through the comparison with the reliable experimental data. The calculated principal and reshock Hugoniots of liquid helium are in good agreement with the gas-gun data. On this basis, we revisit the issue for helium, i.e., the possibility of the instabilities predicted by chemical models at around 2000 GPa and 10 g/cm{sup 3} along the pressure isotherms of 6309, 15 849, and 31 623 K. Our calculations show no indications of instability in this pressure-temperature region, which reconfirm the predictions of previous QMD simulations. The self-diffusion and viscosity coefficients of warm dense helium have been systematically investigated by the QMD simulations. We carefully test the finite-size effects and convergences of statistics, and obtain numerically converged self-diffusion and viscosity coefficients by using the Kubo-Green formulas. The present results have been used to evaluate the existing one component plasma models. Finally, the validation of the Stokes-Einstein relationship for helium in the warm dense regime is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Nucleation and growth by diffusion under Ostwald-Freundlich boundary condition

The critical radius of a nucleus grown by diffusion in a solution is studied thermodynamically as well as kinetically. The thermodynamic growth equation called Zeldovich equation of classical nucleation theory and the kinetic diffusional growth equation combined with the Ostwald-Freundlich boundary condition lead to the same critical radius. However, it should be pointed out that the diffusional equation may lead to a kinetic critical radius that is different from the thermodynamic critical radius, thus indicating the possibility of kinetically controlling the critical radius of a nucleus.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Quasilinear theory of general electromagnetic fluctuations in unmagnetized plasmas

The general quasilinear Fokker-Planck kinetic equation for the plasma particle distribution functions in unmagnetized plasmas is derived, making no restrictions on the frequency of the electromagnetic fluctuations. The derived kinetic particle equation complements our earlier study of the general fluctuation's kinetic equation. For collective plasma eigenmodes and gyrotropic particle distribution functions, the two coupled kinetic equations describe the self-consistent dynamical evolution of the plasma. The limit of weakly damped collective modes correctly reproduces the well-known textbook kinetic particle equation with longitudinal Langmuir and ion-acoustic fluctuations, demonstrating, in particular, the resonant nature of parallel momentum diffusion of particles. In the limit of aperiodic modes, the Fokker-Planck equation contains the nonresonant diffusion of particles in momentum and the parallel and perpendicular momentum drag coefficients. As an application these drag and diffusion coefficients are calculated for extragalactic cosmic ray particles propagating in the unmagnetized intergalactic medium. Whereas for all cosmic rays, the perpendicular momentum diffusion in intergalactic aperiodic fluctuations is negligibly small; cosmic ray protons with energies below 10{sup 5 }GeV are affected by the plasma drag.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Stochastic control with rough paths

We study a class of controlled differential equations driven by rough paths (or rough path realizations of Brownian motion) in the sense of Lyons. It is shown that the value function satisfies a HJB type equation; we also establish a form of the Pontryagin maximum principle. Deterministic problems of this type arise in the duality theory for controlled diffusion processes and typically involve anticipating stochastic analysis. We make the link to old work of Davis and Burstein (Stoch Stoch Rep 40:203–256, 1992) and then prove a continuous-time generalization of Roger’s duality formula [SIAM J Control Optim 46:1116–1132, 2007]. The generic case of controlled volatility is seen to give trivial duality bounds, and explains the focus in Burstein–Davis’ (and this) work on controlled drift. Our study of controlled rough differential equations also relates to work of Mazliak and Nourdin (Stoch Dyn 08:23, 2008).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Theory of Current Transients in Planar Semiconductor Devices: Insights and Applications to Organic Solar Cells

Time-domain current measurements are widely used to characterize semiconductor material properties, such as carrier mobility, doping concentration, carrier lifetime, and the static dielectric constant. It is therefore critical that these measurements be theoretically understood if they are to be successfully applied to assess the properties of materials and devices. In this paper, we derive generalized relations for describing current-density transients in planar semiconductor devices at uniform temperature. By spatially averaging the charge densities inside the semiconductor, we are able to provide a rigorous, straightforward, and experimentally relevant way to interpret these measurements. The formalism details several subtle aspects of current transients, including how the electrode charge relates to applied bias and internal space charge, how the displacement current can alter the apparent free-carrier current, and how to understand the integral of a charge-extraction transient. We also demonstrate how the formalism can be employed to derive the current transients arising from simple physical models, like those used to describe charge extraction by linearly increasing voltage (CELIV) and time-of-flight experiments. In doing so, we find that there is a nonintuitive factor-of-2 reduction in the apparent free-carrier concentration that can be easily missed, for example, in the application of charge-extraction models. Finally, to validate our theory and better understand the different current contributions, we perform a full time-domain drift-diffusion simulation of a CELIV trace and compare the results to our formalism. As expected, our analytic equations match precisely with the numerical solutions to the drift-diffusion, Poisson, and continuity equations. Furthermore, overall, our formalism provides a straightforward and general way to think about how the internal space-charge distribution, the electrode charge, and the externally applied bias translate into a measured current transient in a planar semiconductor device.

36 MATERIALS SCIENCE↗

NASA NTRS1989

An improvement of convection fidelity in Euler calculations

A new solution procedure is presented to solve the Euler equations for steady, compressible, rotational, inviscid flows. The approach is aimed at achieving real inviscid solutions in Euler calculations by eliminating numerical diffusion. The variables in the Euler equations are divided into elliptic and convective quantities, using the Clebsch velocity decomposition. The convective quantities are then transported without numerical contamination using an efficient convection operator, while the elliptic quantities are integrated with a relaxation procedure. This approach provides a generalization of the full potential formulation to rotation Euler physics by allowing variations of convective quantities. Results are demonstrated for several transonic flows in two dimensions.

Lee, K. D.↗

DOE OSTI2013

DYNAMO EFFECTS IN MAGNETOROTATIONAL TURBULENCE WITH FINITE THERMAL DIFFUSIVITY

We investigate the saturation level of hydromagnetic turbulence driven by the magnetorotational instability in the case of vanishing net flux. Motivated by a recent paper of Bodo et al., we here focus on the case of a non-isothermal equation of state with constant thermal diffusivity. The central aim of the paper is to complement the previous result with closure parameters for mean-field dynamo models, and to test the hypothesis that the dynamo is affected by the mode of heat transport. We perform computer simulations of local shearing-box models of stratified accretion disks with approximate treatment of radiative heat transport, which is modeled via thermal conduction. We study the effect of varying the (constant) thermal diffusivity, and apply different vertical boundary conditions (BCs). In the case of impenetrable vertical boundaries, we confirm the transition from mainly conductive to mainly convective vertical heat transport below a critical thermal diffusivity. This transition is however much less dramatic when more natural outflow BCs are applied. Similarly, the enhancement of magnetic activity in this case is less pronounced. Nevertheless, heating via turbulent dissipation determines the thermodynamic structure of accretion disks and clearly affects the properties of the related dynamo. This effect, however, may have been overestimated in previous work, and a careful study of the role played by boundaries will be required.

79 ASTRONOMY AND ASTROPHYSICS↗