Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equations”

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,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Large-eddy simulation of flow through a plane, asymmetric diffuser

A challenge for traditional turbulence modeling, based on the Reynolds averaged Navier-Stokes equations, remains the accurate prediction of 'mild', adverse pressure-gradient driven separation from a smooth surface. With this study we want to explore the capability of large-eddy simulation to predict the separation which occurs on the deflected wall of an asymmetric, plane diffuser with opening angle of 10 deg. The flow through the plane diffuser exhibits some additional interesting physical phenomena which make it a challenging test case. In addition to 'mild' separation about halfway down the deflected ramp, the flow is characterized by a small backflow zone with stalled fluid in the rear part of the expanding section. The turbulent flow entering the diffuser is subject to combined adverse and radial pressure gradients stemming from the convex curvature. Finally the flow recovers into a developed, turbulent channel flow in the outlet section. Obi et al. provide measurements of mean flow, Reynolds stresses, and pressure recovery, which were obtained by means of LDV in a wind tunnel. The objective of this study is to investigate whether LES with the standard dynamic model is able to accurately predict the flow in the one-sided diffuser and to explore the resolution requirements and associated costs.

Kaltenbach, Hans-Jakob↗
DOE OSTI2013

Penetration of lower hybrid current drive waves in tokamaks

Lower hybrid (LH) ray propagation in toroidal plasma is shown to be controlled by combination of the azimuthal spectrum launched by the antenna, the poloidal variation of the magnetic field, and the scattering of the waves by the drift wave fluctuations. The width of the poloidal and radial radio frequency wave spectrum increases rapidly as the rays penetrate into higher density and scatter from the drift waves. The electron temperature gradient (ETG) spectrum is particularly effective in scattering the LH waves due to its comparable wavelengths and phase velocities. ETG turbulence is also driven by the radial gradient of the electron current profile giving rise to an anomalous viscosity spreading the LH driven plasma currents. The LH wave scattering is derived from a Fokker-Planck equation for the distribution of the ray trajectories with diffusivities derived from the drift wave fluctuations. The condition for chaotic diffusion for the rays is derived. The evolution of the poloidal and radial mode number spectrum of the lower hybrid waves are both on the antenna spectrum and the spectrum of the drift waves. Antennas launching higher poloidal mode number spectra drive off-axis current density profiles producing negative central shear [RS] plasmas with improved thermal confinement from ETG transport. Core plasma current drive requires antennas with low azimuthal mode spectra peaked at m = 0 azimuthal mode numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Transport coefficients for magnetic-field evolution in inviscid magnetohydrodynamics

The magnetized resistivity and electrothermal tensors when substituted into the induction equation lead to electrothermal magnetic field generation, resistive magnetic diffusion, and magnetic field advection due to resistivity gradients, temperature gradients and currents. Here, the advection terms driven by temperature gradient and current have cross-field components (perpendicular to both the magnetic field and the driving term) that depend on significantly modified versions of Braginskii’s transport coefficients. The improved fits to Braginskii’s coefficients given by Epperlein and Haines and Ji and Held give physically incorrect results for cross-field advection at small Hall parameters (product of cyclotron frequency and collision time). The errors in Epperlein and Haines’ fits are particularly severe, giving increasing advection velocities below a Hall parameter of one when they should decrease linearly to zero. Epperlein and Haines’ fits can also give erroneous advection terms due to variations in effective atomic number. The only serious error in Braginskii’s fits is an overestimate in advection due to perpendicular resistivity. New fits for the cross-field advection terms are obtained from a direct numerical solution of the Fokker-Planck equation and Ji and Held’s higher order expansion approach that are continuous functions of the effective atomic number.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

On an origin of numerical diffusion: Violation of invariance under space-time inversion

The invariant properties of the convection equation du/dt + adu/dx = 0 (where d is the partial differential operator) with respect to spatial reflection, time reversal, and space-time inversion are studied. Generally, a finite-difference analog of this equation may possess some or none of these properties. It is shown that, under certain conditions, the von Neumann amplification factor of an analog satisfies a special relation for each invariant property this analog possesses. Particularly, an analog is neutrally stable and thus free of numerical diffusion if it possesses the invariant property related to space-time inversion. It is also explained why generally (1) an upwind scheme possesses neither the invariant property related to spatial reflection nor that related to space-time inversion, and (2) an explicit scheme possesses neither the invariant property related to time reversal nor that related to space-time inversion. Extension to the viscous case and a remarkable connection between the current work and a new numerical framework for solving conservation laws are also discussed.

Chang, Sin-Chung↗
NASA NTRS1973

Transient permeation of organic vapors through elastomeric membranes

The permeation of benzene and acetone vapors through sulfur-cured natural rubber was studied by the time-lag method. The zero concentration diffusion coefficient, was obtained by the early-time method. The Frisch time lag equation was utilized to estimate both the solubility coefficient and the additional parameter required to define the concentration-dependence of the diffusion coefficient. This form of concentration dependence was manifested by the corresponding permeability coefficient values. At low entering penetrant pressure where the transport coefficients are constant, indirect evidence was obtained that zero diffusion coefficient is the mechanically correct diffusion coefficient. The solubility coefficient values calculated for benzene vapor in natural rubber are in reasonable agreement with published equilibrium sorption data for a similar rubber compound. At higher entering penetrant pressures, average diffusion coefficients obtained at steady state tended to be larger than the corresponding average diffusion coefficients derived from the time-lags.

Curry, J. E.↗
NASA NTRS1983

A unified mathematical theory of electrophoretic processes

A mathematical theory is presented which shows that each of the four classical electrophoretic modes (zone electrophoresis, moving boundary electrophoresis, isotachophoresis, and isoelectric focusing) is based on the same general principles and can collectively be described in terms of a single set of equations. This model can predict the evolution of the four electrophoretic modes as a function of time. The model system is one-dimensional, neglecting the effects of electroosmosis, temperature gradients, and any bulk flows of liquid. The model is based on equations which express the components' dissociation equilibria, the mass transport due to electromigration and diffusion, electroneutrality, and the conservation of mass and charge. The model consists of a system of coupled partial differential and nonlinear algebraic equations which can be solved numerically by use of a computer. The versatility of this model was verified using an example of a three-component system containing cacodylate, tris hydroxylmethylaminomethane, and histidine. Results show that this model not only correctly predicts the characteristic features of each electrophoretic mode, but also gives details of the concentration, pH, and conductivity profiles not easily amenable to direct experimental measurement.

Bier, M.↗
NASA NTRS1985

Calculation of three-dimensional, viscous flow through turbomachinery blade passages by parabolic marching

The three-dimensional compressible Navier-Stokes equations are formulated in a rotating coordinate system, so as to include centrifugal and Coriolis forces. The equations are parabolized by using a previously calculated inviscid static pressure field. The thin layer Navier-Stokes approximation, which neglects streamwise diffusion, is used. A body-fitted coordinate system is used. The streamwise momentum equation is uncoupled from the cross-stream momentum equation by using contravariant momentum components, and then using the contravariant velocity components as primary unknowns. To reduce problems with small separated regions, the Reyhner and Flugge-Lotz approximation is used. The energy equation is included to allow for calculation of heat transfer. The flow may be laminar, or a simple eddy-viscosity turbulence may be used. A number of curved ducts and an axial stator were analyzed, including cases for which experimental data are available.

Katsanis, T.↗
NASA NTRS1985

Calculation of three-dimensional, viscous flow through turbomachinery blade passages by parabolic marching

The three-dimensional compressible Navier-Stokes equations are formulated in a rotating coordinate system, so as to include centrifugal and Coriolis forces. The equations are parabolized by using a previously calculated inviscid static pressure field. The thin layer Navier-Stokes approximation, which neglects streamwise diffusion, is used. A body-fitted coordinate system is used. The streamwise momentum equation is uncoupled from the cross-stream momentum equation by using contravariant momentum components, and then using the contravariant velocity components as primary unknowns. To reduce problems with small separating regions, the Reyhner and Flugge-Lotz approximation is used. The energy equation is included to allow for calculation of heat transfer. The flow may be laminar, or a simple eddy-viscosity turbulence may be used. A number of curved ducts and an axial stator were analyzed, including cases for which experimental data are available.

Katsanis, T.↗
DOE OSTI2013

Probability Density Function Method for Langevin Equations with Colored Noise

We present a novel method to derive closed-form, computable PDF equations for Langevin systems with colored noise. The derived equations govern the dynamics of joint or marginal probability density functions (PDFs) of state variables, and rely on a so-called Large-Eddy-Diffusivity (LED) closure. We demonstrate the accuracy of the proposed PDF method for linear and nonlinear Langevin equations, describing the classical Brownian displacement and dispersion in porous media.

PDF method, uncertainty quantification, Langevin e↗
DOE OSTI2015

A computational study of diffusion in a glass-forming metallic liquid

In this study, liquid phase diffusion plays a critical role in phase transformations (e.g. glass transformation and devitrification) observed in marginal glass forming systems such as Al-Sm. Controlling transformation pathways in such cases requires a comprehensive description of diffusivity, including the associated composition and temperature dependencies. In our computational study, we examine atomic diffusion in Al-Sm liquids using ab initio molecular dynamics (AIMD) and determine the diffusivities of Al and Sm for selected alloy compositions. Non-Arrhenius diffusion behavior is observed in the undercooled liquids with an enhanced local structural ordering. Through assessment of our AIMD result, we construct a general formulation for Al-Sm liquid, involving a diffusion mobility database that includes composition and temperature dependence. A Volmer-Fulcher-Tammann (VFT) equation is adopted for describing the non-Arrhenius behavior observed in the undercooled liquid. Furthermore, the composition dependence of diffusivity is found quite strong, even for the Al-rich region contrary to the sole previous report on this binary system. The model is used in combination with the available thermodynamic database to predict specific diffusivities and compares well with reported experimental data for 0.6 at.% and 5.6 at.% Sm in Al-Sm alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Identifiability analysis of rotational diffusion tensor and electronic transition moments measured in time-resolved fluorescence depolarization experiment

The subject of this paper is studies of the deterministic identifiability of molecular parameters, such as rotational diffusion tensor components and orientation of electronic transition moments, resulting from the time-resolved fluorescence anisotropy experiment. In the most general case considered, a pair of perpendicularly polarized emissions enables the unique determination of all the rotational diffusion tensor's principal components. The influence of the tensor's symmetry and the associated degeneration of its eigenvalues on the identifiability of the electronic transitions moments is systematically investigated. The analysis reveals that independently of the rotational diffusion tensor's symmetry, the transition moments involved in photoselection and emission processes cannot be uniquely identified without a priori information about their mutual orientation or their orientation with respect to the principal axes of the tensor. Moreover, it is shown that increasing the symmetry of the rotational diffusion tensor deteriorates the degree of the transition moments identifiability. To obtain these results analytically, a novel approach to solve bilinear system of equations for Markov parameters is applied. The effect of the additional information, obtained from fluorescence measurements for different molecular mobilities, to improve the identifiability at various levels of analysis is shown. The effectiveness and reliability of the target analysis method for experimental determination of the molecular parameters is also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2013

DNS, Enstrophy Balance, and the Dissipation Equation in a Separated Turbulent Channel Flow

The turbulent flows through a plane channel and a channel with a constriction (2-D hill) are numerically simulated using DNS and RANS calculations. The Navier-Stokes equations in the DNS are solved using a higher order kinetic energy preserving central schemes and a fifth order accurate upwind biased WENO scheme for the space discretization. RANS calculations are performed using the NASA code CFL3D with the komega SST two-equation model and a full Reynolds stress model. Using DNS, the magnitudes of different terms that appear in the enstrophy equation are evaluated. The results show that the dissipation and the diffusion terms reach large values at the wall. All the vortex stretching terms have similar magnitudes within the buffer region. Beyond that the triple correlation among the vorticity and strain rate fluctuations becomes the important kinematic term in the enstrophy equation. This term is balanced by the viscous dissipation. In the separated flow, the triple correlation term and the viscous dissipation term peak locally and balance each other near the separated shear layer region. These findings concur with the analysis of Tennekes and Lumley, confirming that the energy transfer terms associated with the small-scale dissipation and the fluctuations of the vortex stretching essentially cancel each other, leaving an equation for the dissipation that is governed by the large-scale motion.

Balakumar, Ponnampalam↗
NASA NTRS1978

Molecular evolution of contracting clouds - Basic methods and initial results

The relationship between the dynamics of the interstellar gas and the thermal and chemical effects associated with interstellar molecules and dust is investigated. The evolution of a rather massive isolated initially diffuse cloud under self-gravity is studied, using the equations of hydrodynamics; only radial motions are considered, and the heat, chemical, and radiative-transfer equations are solved simultaneously with the hydrodynamic equations. The relevant chemistry is described along with the thermal model, the radiative-transfer process, and the numerical methods employed. Results for a contracting cloud are discussed in terms of the problem of initial conditions, the dynamical evolution of the cloud, its chemical and thermal evolution, time scales, and column densities. It is shown that the chemical evolution of a massive contracting diffuse cloud is sensitive to such physical properties as temperature and ion abundances, that warm and cool versions of a typical cloud evolve differently, and that the physical origin of this effect is the level of heating due to H2 formation on interstellar dust grains.

Gerola, H.↗
NASA NTRS1999

Relativistic Electrons Observed at UARS and the Interpretation of their Storm-Associated Intensity Variations

The High Energy Particle Spectrometer (HEPS) instrument on the Upper Atmosphere Research Satellite (UARS) provides a database of electron intensities well resolved in energy and pitch-angle. Because of its 57 deg. orbital inclination, UARS encounters with magnetic shells L greater than 2 occur quite far off-equator (B/B (sub 0) greater than 9), corresponding to equatorial pitch angle alpha (sub 0) greater than 20 deg. Data acquired by HEPS (October 1991 through September 1994) span the declining phase of Solar Cycle 22. To reveal the storm-associated time dependence of relativistic electron intensities over the wide range of energies (50 keV to 5 MeV) covered by HEPS, we divide the daily average of the measured spectrum at a given L value (bin width = 0.25) by the corresponding 500-day average and plot the results with a color scale that spans only 2.5 decades. The data show that our off-equatorial electron intensities typically increase with time after the end of recovery phase (not during main phase or recovery phase) of each geomagnetic storm. The delay in off-equatorial energetic electron response and the subsequent lifetime of the corresponding electron flux enhancement seem to increase with particle energy above 300 keV. The trend below 300 keV seems to be opposite, such that the delay varies inversely with electron energy. Our working hypothesis for interpretation is that stormtime radial transport tends to increase the phase-space densities of trapped relativistic electrons but typically leads to a flux increases at specified energies only as the current (as indicated by Dst) decays. Flux enhancements in early recovery phase are greatest for equatorially mirroring electrons, and to pitch-angle anisotropies are initially large. Subsequent pitch-angle diffusion broadens the flux enhancement to particles that mirror off equator, thus gradually increasing low-altitude electron intensities (as detected by HEPS/UARS) on time scales equal to about 20% of corresponding lifetimes against diffusion into the loss cone. Alternative interpretations will also be examined.

Pesnell, W. D.↗
DOE OSTI2015

Hybrid approaches for multiple-species stochastic reaction–diffusion models

Reaction–diffusion models are used to describe systems in fields as diverse as physics, chemistry, ecology and biology. The fundamental quantities in such models are individual entities such as atoms and molecules, bacteria, cells or animals, which move and/or react in a stochastic manner. If the number of entities is large, accounting for each individual is inefficient, and often partial differential equation (PDE) models are used in which the stochastic behaviour of individuals is replaced by a description of the averaged, or mean behaviour of the system. In some situations the number of individuals is large in certain regions and small in others. In such cases, a stochastic model may be inefficient in one region, and a PDE model inaccurate in another. To overcome this problem, we develop a scheme which couples a stochastic reaction–diffusion system in one part of the domain with its mean field analogue, i.e. a discretised PDE model, in the other part of the domain. The interface in between the two domains occupies exactly one lattice site and is chosen such that the mean field description is still accurate there. In this way errors due to the flux between the domains are small. Our scheme can account for multiple dynamic interfaces separating multiple stochastic and deterministic domains, and the coupling between the domains conserves the total number of particles. The method preserves stochastic features such as extinction not observable in the mean field description, and is significantly faster to simulate on a computer than the pure stochastic model. - Highlights: • A novel hybrid stochastic/deterministic reaction–diffusion simulation method is given. • Can massively speed up stochastic simulations while preserving stochastic effects. • Can handle multiple reacting species. • Can handle moving boundaries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

The diffusion of contaminant through an accretion disc

The manner in which a trace contaminant diffuses through an accretion disk is calculated. For a steady disk in which the ratio of diffusivity to viscosity, zeta, is taken to be constant through the disk and for which the surface density distribution is taken to be a power law, analytic solutions to the time-dependent equations governing contaminant evolution are presented. The amount of contaminant which can diffuse its way outwards against the bulk inward accretion flow is calculated, and it is shown that this is a sensitive function of a and zeta.

Clarke, C. J.↗
NASA NTRS1991

Quasi-linear pitch angle and energy diffusion of pickup ions near Comet Halley

The process of ion pickup in the environment of Halley's comet is studied in order to see if velocity diffusion driven by the observed level of turbulence can explain the observed development of the implanted ion distribution. The theoretical description used is based on a quasi-linear approach and considers the implantation and transport of cometary ions along solar wind flow lines. To make such a study requires some way of extrapolating the measurements on the Giotto trajectory into the upstream region; models for mass loading and turbulence are used. A simplified kinetic equation describing the source, convection, and quasi-linear velocity diffusion of the heavy cometary ions is solved numerically along flow lines parallel to the sun-comet line. Full two-dimensional (pitch angle and velocity) distributions are obtained at positions along the Giotto trajectory, which are compared with measurements. This study finds that quasi-linear theory, with the empirical model for the observed turbulence level, produces the right order of pitch angle diffusion.

Huddleston, D. E.↗
NASA NTRS1992

Turbulent transport of a passive-scalar field by using a renormalization-group method

A passive-scalar field is considered to evolve under the influence of a turbulent fluid governed by the Navier-Stokes equation. Turbulent-transport coefficients are calculated by small-scale elimination using a renormalization-group method. Turbulent processes couple both the viscosity and the diffusivity. In the absence of any correlation between the passive-scalar fluctuations and any component of the fluid velocity, the renormalized diffusivity is essentially the same as if the fluid velocity were frozen, although the renormalized equation does contain higher-order nonlinear terms involving viscosity. This arises due to the nonlinear interaction of the velocity with itself. In the presence of a finite correlation, the turbulent diffusivity becomes coupled with both the velocity field and the viscosity. There is then a dependence of the turbulent decay of the passive scalar on the turbulent Prandtl number.

Hossain, Murshed↗