Search NASA⌕ Search

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.

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

At least 1,243 records · Page 69

What Would a Dynamo Theorist like to Know About the Dynamics of the Solar Convection Zone?

Observations about the current status of solar dynamo theory are given. The induction equation for magnetic field is solved using assumed velocities and parametric representations of the inductive or diffusive effects of velocities. The equations of motion governing these flow are not solved in parallel. Results from global compressible convection models are discussed. Differential rotation and convection are also investigated.

Gilman, P. A.↗

Equivalence of Fluctuation Splitting and Finite Volume for One-Dimensional Gas Dynamics

The equivalence of the discretized equations resulting from both fluctuation splitting and finite volume schemes is demonstrated in one dimension. Scalar equations are considered for advection, diffusion, and combined advection/diffusion. Analysis of systems is performed for the Euler and Navier-Stokes equations of gas dynamics. Non-uniform mesh-point distributions are included in the analyses.

Wood, William A.↗

Heat transfer, diffusion, and evaporation

Although it has long been known that the differential equations of the heat-transfer and diffusion processes are identical, application to technical problems has only recently been made. In 1916 it was shown that the speed of oxidation of the carbon in iron ore depends upon the speed with which the oxygen of the combustion air diffuses through the core of gas surrounding the carbon surface. The identity previously referred to was then used to calculate the amount of oxygen diffusing to the carbon surface on the basis of the heat transfer between the gas stream and the carbon surface. Then in 1921, H. Thoma reversed that procedure; he used diffusion experiments to determine heat-transfer coefficients. Recently Lohrisch has extended this work by experiment. A technically very important application of the identity of heat transfer and diffusion is that of the cooling tower, since in this case both processes occur simultaneously.

GASES - DIFFUSION↗

Thermodynamic precursors, liquid-liquid transitions, dynamic and topological anomalies in densified liquid germania

The thermodynamic, dynamic, structural, and rigidity properties of densified liquid germania (GeO{sub 2}) have been investigated using classical molecular dynamics simulation. We construct from a thermodynamic framework an analytical equation of state for the liquid allowing the possible detection of thermodynamic precursors (extrema of the derivatives of the free energy), which usually indicate the possibility of a liquid-liquid transition. It is found that for the present germania system, such precursors and the possible underlying liquid-liquid transition are hidden by the slowing down of the dynamics with decreasing temperature. In this respect, germania behaves quite differently when compared to parent tetrahedral systems such as silica or water. We then detect a diffusivity anomaly (a maximum of diffusion with changing density/volume) that is strongly correlated with changes in coordinated species, and the softening of bond-bending (BB) topological constraints that decrease the liquid rigidity and enhance transport. The diffusivity anomaly is finally substantiated from a Rosenfeld-type scaling law linked to the pair correlation entropy, and to structural relaxation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The prediction of nozzle performance and heat transfer in hydrogen/oxygen rocket engines with transpiration cooling, film cooling, and high area ratios

An advanced engineering computational model has been developed to aid in the analysis and design of hydrogen/oxygen chemical rocket engines. The complete multi-species, chemically reacting and diffusing Navier-Stokes equations are modelled, finite difference approach that is tailored to be conservative in an axisymmetric coordinate system for both the inviscid and viscous terms. Demonstration cases are presented for a 1030:1 area ratio nozzle, a 25 lbf film cooled nozzle, and transpiration cooled plug-and-spool rocket engine. The results indicate that the thrust coefficient predictions of the 1030:1 nozzle and the film cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements when all of the chemical reaction and diffusion terms are considered. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. The Soret thermal diffusion term is demonstrated to have a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle in both the laminar flow 1030:1 nozzle and the turbulent plug-and-spool rocket engine analysis cases performed. Further, the Soret term was shown to represent a significant fraction of the diffusion fluxes occurring in the transpiration cooled rocket engine.

Kacynski, Kenneth J.↗

EFFECT OF ADIABATIC FOCUSING ON THE DIFFUSION OF ENERGETIC CHARGED PARTICLES

An approximate analytic description of a diffusion coefficient, including the effect of adiabatic focusing, has been developed. This description is formulated with the aid of stochastic differential equations and the steady perturbation solution of the Fokker–Plank transport equation. The analytical formula is based on three important assumptions. First, the pitch-angle diffusion coefficient is set to be separable from the spatial coordinate and the pitch-angle cosine. Second, the spatial dependence of the ratio between the mean-free path and focusing length is assumed to be weak. Third, the pitch-angle distribution relaxes quickly to a steady state. The new analytic formula could be applied to calculate the spatial diffusion coefficient in the interplanetary and interstellar space.

79 ASTRONOMY AND ASTROPHYSICS↗

Quasilinear theory for inhomogeneous plasma

Here, this paper presents quasilinear theory (QLT) for a classical plasma interacting with inhomogeneous turbulence. The particle Hamiltonian is kept general; for example, relativistic, electromagnetic and gravitational effects are subsumed. A Fokker–Planck equation for the dressed ‘oscillation-centre’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, interaction with the background fields and ponderomotive effects simultaneously. The local diffusion coefficient is manifestly positive-semidefinite. Waves are allowed to be off-shell (i.e. not constrained by a dispersion relation), and a collision integral of the Balescu–Lenard type emerges in a form that is not restricted to any particular Hamiltonian. This operator conserves particles, momentum and energy, and it also satisfies the H -theorem, as usual. As a spin-off, a general expression for the spectrum of microscopic fluctuations is derived. For on-shell waves, which satisfy a quasilinear wave-kinetic equation, the theory conserves the momentum and energy of the wave–plasma system. The action of non-resonant waves is also conserved, unlike in the standard version of QLT. Dewar's oscillation-centre QLT of electrostatic turbulence is proven formally as a particular case and given a concise formulation. Also discussed as examples are relativistic electromagnetic and gravitational interactions, and QLT for gravitational waves is proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Influence of the interplanetary magnetic field on cometary and primordial dust orbits - Applications of Lorentz Scattering

The equations describing the change in orbital elements of interplanetary dust due to Lorentz-force accelerations are presented in a simplified form. Such accelerations depend on the charge state of the dust; results of theoretical calculations for five possible dust materials are presented. Under present-day conditions, it is possible that semiconducting material such as graphite might carry a net voltage near zero, compared with a roughly 10-V charge expected for other grains. The scattering of dust by a randomly changing magnetic field can be viewed analogously to the dust diffusing in space; the equations presented thus can be used to interpret observations of the present distribution of dust in terms of its possible sources and sinks. The stronger magnetic fields of the early solar system would have led to more vigorous scattering of the dust; particles as large as 1 mm could have been significantly transported by Lorentz scattering during this time.

Consolmagno, G. J.↗

An asymptotic description of transient settling and ultrafiltration of colloidal dispersions

The present model for the sedimentation of colloidal systems includes a diffusion term in the governing equation which, in the regions above the sediment, acts as small perturbation to the Kynch (1952) theory. Within the sediment, diffusion is comparable to convection due to the the high solid volume fraction. An application of the method of matched asymptotic expansions to the conservation equation allows a complete description of the settling process to be formulated, with specific attention to volume-fraction evolution in the sediment.

Davis, K. E.↗

Remobilizing the Interfaces of Thermocapillary Driven Bubbles Retarded by the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.

Palaparthi, Ravi↗

Diffusive cosmic-ray acceleration at relativistic shock waves with magnetostatic turbulence

The analytical theory of diffusive cosmic-ray acceleration at parallel stationary shock waves with magnetostatic turbulence is generalized to arbitrary shock speeds V{sub s}=β{sub 1}c, including, in particular, relativistic speeds. This is achieved by applying the diffusion approximation to the relevant Fokker–Planck particle transport equation formulated in the mixed comoving coordinate system. In this coordinate system, the particle's momentum coordinates p and μ=p{sub ∥}/p are taken in the rest frame of the streaming plasma, whereas the time and space coordinates are taken in the observer's system. For magnetostatic slab turbulence, the diffusion–convection transport equation for the isotropic (in the rest frame of the streaming plasma) part of the particle's phase space density is derived. For a step-wise shock velocity profile, the steady-state diffusion–convection transport equation is solved. For a symmetric pitch-angle scattering Fokker–Planck coefficient, D{sub μμ}(−μ)=D{sub μμ}(μ), the steady-state solution is independent of the microphysical scattering details. For nonrelativistic mono-momentum particle injection at the shock, the differential number density of accelerated particles is a Lorentzian-type distribution function, which at large momenta approaches a power-law distribution function N(p⩾p{sub c})∝p{sup −ξ} with the spectral index ξ(β{sub 1})=1+[3/(Γ{sub 1}√(r{sup 2}−β{sub 1}{sup 2})−1)(1+3β{sub 1}{sup 2})]. For nonrelativistic (β{sub 1}≪1) shock speeds, this spectral index agrees with the known result ξ(β{sub 1}≪1)≃(r+2)/(r−1), whereas for ultrarelativistic (Γ{sub 1}≫1) shock speeds the spectral index value is close to unity.

79 ASTRONOMY AND ASTROPHYSICS↗

X-ray Studies of Interfacial Strontium–Extractant Complexes in a Model Solvent Extraction System

The interfacial behavior of a model solvent extraction liquidliquid system, consisting of solutions of dihexadecyl phosphate (DHDP) in dodecane and SrCl2 in water, was studied to determine the structure of the interfacial ionextractant complex and its variation with pH. Previous experiments on a similar extraction system with ErCl3 demonstrated that the kinetics of the extraction process could be greatly retarded by cooling through an adsorption transition, thus providing a method to immobilize ionextractant complexes at the interface and further characterize them with X-ray interface-sensitive techniques. Here, we use this same method to study the SrCl2 system. X-ray reflectivity and fluorescence near total reflection measured the molecular-scale interfacial structure above and below the adsorption transition for a range of pH. Below the transition, DHDP molecules form a homogeneous monolayer at the interface with Sr2+ coverage increasing from zero to saturation (one Sr2+ per two DHDP) within a narrow range of pH. Experimental values of Sr2+ interfacial density determined from fluorescence measurements are larger than those from reflectivity measurements. Although both techniques probe Sr2+ bound to DHDP, only the fluorescence provides adequate sensitivity to Sr2+ in the diffuse double layer. A Stern equation determines the Sr2+ binding constant from the reflectivity measurements and the additional Sr2+ measured in the diffuse double layer is accounted for by GouyChapman theory. Above the transition temperature, a dilute concentration of DHDPSr complexes resides at the interface, even for temperatures far above the transition. A comparison is made of the structure of the interfacial ionextractant complex for this divalent metal ion to recent results on trivalent Er3+ metal ions, which provides insight into the role of metal ion charge on the structure of interfacial ionextractant complexes, as well as implications for extraction of these two differently charged ions

Bu, Wei↗

Peridynamic thermal diffusion

This study presents the derivation of ordinary state-based peridynamic heat conduction equation based on the Lagrangian formalism. The peridynamic heat conduction parameters are related to those of the classical theory. An explicit time stepping scheme is adopted for numerical solution of various benchmark problems with known solutions. It paves the way for applying the peridynamic theory to other physical fields such as neutronic diffusion and electrical potential distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Luminescent coupling in planar opto-electronic devices

Effects of luminescent coupling are observed in monolithic 5 V, five-junction GaAs phototransducers. Power conversion efficiency was measured at 61.6% ± 3% under the continuous, monochromatic illumination for which they were designed. Modeling shows that photon recycling can account for up to 350 mV of photovoltage in these devices. Drift-diffusion based simulations including a luminescent coupling term in the continuity equation show a broadening of the internal quantum efficiency curve which agrees well with experimental measurements. Luminescent coupling is shown to expand the spectral bandwidth of the phototransducer by a factor of at least 3.5 for devices with three or more junctions, even in cases where multiple absorption/emission events are required to transfer excess carriers into the limiting junction. We present a detailed description of the novel luminescent coupling modeling technique used to predict these performance enhancements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Weakly dynamic dark energy via metric-scalar couplings with torsion

We study the dynamical aspects of dark energy in the context of a non-minimally coupled scalar field with curvature and torsion. Whereas the scalar field acts as the source of the trace mode of torsion, a suitable constraint on the torsion pseudo-trace provides a mass term for the scalar field in the effective action. In the equivalent scalar-tensor framework, we find explicit cosmological solutions representing dark energy in both Einstein and Jordan frames. We demand the dynamical evolution of the dark energy to be weak enough, so that the present-day values of the cosmological parameters could be estimated keeping them within the confidence limits set for the standard LCDM model from recent observations. For such estimates, we examine the variations of the effective matter density and the dark energy equation of state parameters over different redshift ranges. In spite of being weakly dynamic, the dark energy component differs significantly from the cosmological constant, both in characteristics and features, for e.g. it interacts with the cosmological (dust) fluid in the Einstein frame, and crosses the phantom barrier in the Jordan frame. We also obtain the upper bounds on the torsion mode parameters and the lower bound on the effective Brans-Dicke parameter. The latter turns out to be fairly large, and in agreement with the local gravity constraints, which therefore come in support of our analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Coarse-Graining Conformational Dynamics with Multidimensional Generalized Langevin Equation: How, When, and Why

A data-driven ab initio generalized Langevin equation (AIGLE) approach is developed to learn and simulate high-dimensional, heterogeneous, coarse-grained (CG) conformational dynamics. Constrained by the fluctuation–dissipation theorem, the approach can build CG models in dynamical consistency (DC) with all-atom molecular dynamics. Here, we also propose practical criteria for AIGLE to enforce long-term DC. Case studies of a toy polymer, with 20 CG sites, and the alanine dipeptide, with two dihedral angles, elucidate why one should adopt AIGLE or its Markovian limit for modeling CG conformational dynamics in practice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗