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,261 records · Page 70

Satellite observations of 557.7 NM airglow as a diagnostic technique to study the dynamics of the lower thermosphere

Photometric observations of the airglow limb at 557.7 nm have been made from the ISIS 2 satellite since 1971. From an analysis of 32,000 airglow limb observations obtained during the first two years of operation, the major temporal and spatial characteristics of the night airglow have been identified. Maxima occur in mid-October and mid-April at all latitudes. There is a marked variation with latitude during the equinox periods, with peak intensities near 35 degrees in autumn and near 25 degrees in spring. It is possible to gain considerable insight into the behavior of the lower thermosphere by combining the airglow observations with a simple model of the photochemical and dynamical processes. In this way a pattern of the latitude and temporal variation of atomic oxygen concentration at the peak of the layer as well as the eddy diffusion in this altitude region have been inferred. The concentration of atomic oxygen shows midlatitude peaks near equinox, and minima about one month after solstice. Average values at midlatitudes are about 1.3 times those at the equator. The amplitude of the semiannual component is about 20% of the mean value at 35 degrees as compared to 10% at the equator, and the amplitude of 12-month component is about 12% at 35 deg and 5% at the equator. The semiannual component has a minimum that occurs approximately 3 weeks (plus or minus) after solstice at all latitudes, while the 12-month component peaks in the summer hemisphere at solstice (plus or minus 2 months).

Cogger, L. L.↗

Prediction of heat release effects on a mixing layer

A fully second-order closure model for turbulent reacting flows is suggested based on Favre statistics. For diffusion flames the local thermodynamic state is related to a single conserved scalar. The properties of pressure fluctuations are analyzed for turbulent flows with fluctuating density. Closure models for pressure correlations are discussed and modelled transport equations for Reynolds stresses, turbulent kinetic energy dissipation, density-velocity correlations, scalar moments and dissipation are presented and solved, together with the mean equations for momentum and mixture fraction. Solutions of these equations are compared with the experimental data for high heat release free mixing layers of fluorine and hydrogen in a nitrogen diluent.

Farshchi, M.↗

Implications of a strongly peaked polar magnetic field

Using the flux-transport equation in the absence of sources, the relation between a highly peaked polar magnetic field and the poleward meridional flow that concentrates it is studied. If the maximum flow speed v(m) greatly exceeds the effective diffusion speed, then the field has a quasi-equilibrium configuration in which the poleward convection of flux via meridional flow approximately balances the equatorward spreading via supergranular diffusion. Recent measurements would then give v(m) of about 7 m/s.

Sheeley, N. R., Jr.↗

Linear tearing modes of a forced current-sheet equilibrium

Linear tearing modes are studied in a nearly singular forced current-sheet equilibrium, such as could result from global magnetic forces in the solar corona. Growth rates for the tearing modes, determined by solving the linearized reduced MHD (Strauss) equations numerically, were found to scale as (gamma)tau(d) = about S exp 4/5 k(y) exp 4/5, where S is the Lundquist number, k(y) is the wavenumber, and tau(d) is the classical resistive diffusion time. This scaling is in agreement with predictions from analytical theory. Because of the faster S scaling of these modes compared to the tearing modes of a diffuse current-sheet equilibrium, the modes have much higher growth rates (by a factor of about 10,000) for coronal values of S (about 10 to the 12th). For coronal parameters, the growth times of these new modes are estimated to be on the order of several hours to days, as compared to growth times of months to years for the tearing modes in a diffuse current-shear equilibrium. The growth times are comparable to reconnection times scales required in models of coronal heating by magnetic field dissipation.

Liewer, Paulett C.↗

Numerical modeling of physical vapor transport under microgravity conditions: Effect of thermal creep and stress

One of the most promising applications of microgravity (micro-g) environments is the manufacture of exotic and high-quality crystals in closed cylindrical ampoules using physical vapor transport (PVT) processes. The quality enhancements are believed to be due to the absence of buoyant convection in the weightless environment - resulting in diffusion-limited transport of the vapor. In a typical experiment, solid-phase sample material is initially contained at one end of the ampoule. The sample is made to sublime into the vapor phase and deposit onto the opposite end by maintaining the source at an elevated temperature with respect to the deposit. Identification of the physical factors governing both the rates and uniformity of crystal growth, and the optimization of the micro-g technology, will require an accurate modeling of the vapor transport within the ampoule. Previous micro-g modeling efforts have approached the problem from a 'classical' convective/diffusion formulation, in which convection is driven by the action of buoyancy on thermal and solutal density differences. The general conclusion of these works have been that in low gravity environments the effect of buoyancy on vapor transport is negligible, and vapor transport occurs in a diffusion-limited mode. However, it has been recently recognized than in the non-isothermal (and often low total pressure) conditions encountered in ampoules, the commonly-assumed no-slip boundary condition to the differential equations governing fluid motion can be grossly unrepresentative of the actual situation. Specifically, the temperature gradients can give rise to thermal creep flows at the ampoule side walls. In addition, temperature gradients in the vapor itself can, through the action of thermal stress, lead to bulk fluid convection.

Mackowski, Daniel W.↗

COBE DIRBE near-infrared polarimetry of the zodiacal light: Initial results

This Letter describes near-infrared polarimetry of the zodiacal light at 2.2 micrometers, measured with the Diffuse Infrared Background Experiment (DIRBE) aboard the Cosmic Background Explorer (COBE) spacecraft. The polarization is due to scattering of sunlight. The polarization vector is perpendicular to the scattering plane, and its observed amplitude on the ecliptic equator at an elongation of 90 deg and ecliptic longitude of 10 deg declines from 12.0 +/- 0.4% at 1.25 micrometers to 8.0 +/- 0.6% at 3.5 micrometers (cf. 16% in the visible); the principal source of uncertainty is photometric noise due to stars. The observed near-infrared colors at this location are redder than Solar, but at 3.5 micrometers this is due at least in part to the thermal emission contribution from the interplanetary dust. Mie theory calculations show that both polarizations and colors are important in constraining models of interplanetary dust.

Berriman, G. B.↗

Symbolic Computational Approach to the Marangoni Convection Problem With Soret Diffusion

A recently reported solution for stationary stability of a thermosolutal system with Soret diffusion is re-derived and examined using a symbolic computational package. Symbolic computational languages are well suited for such an analysis and facilitate a pragmatic approach that is adaptable to similar problems. Linearization of the equations, normal mode analysis, and extraction of the final solution are performed in a Mathematica notebook format. An exact solution is obtained for stationary stability in the limit of zero gravity. A closed form expression is also obtained for the location of asymptotes in relevant parameter, (Sm(sub c), Mac(sub c)), space. The stationary stability behavior is conveniently examined within the symbolic language environment. An abbreviated version of the Mathematica notebook is given in the Appendix.

Skarda, J. Raymond↗

Physics of Reconnection and MMS Mission

Reconnection is the most important process driving the Earth's magnetosphere. Key to the success of the MMS science plan is the coupling of theory and observation. Determining the kinetic processes occurring in the diffusion region and physical parameters that control the rate of magnetic reconnection are among primary objectives of the MMS mission. Analysis of the role played by particle inertial effects in the diffusion region where the plasma is unmagnetized will be presented. The reconnection electric field in he diffusion region is supported primarily by particle non-gyrotropic effects. At the quasi-steady stage the reconnection electric field serves to accelerate and heat the incoming plasma population to maintain the current flow in the diffusion region the pressure balance. The primary mechanism controlling the dissipation in the vicinity of the reconnection site is incorporated into the fluid description in terms of non-gyrotropic corrections to the. induction and energy equations. The results of kinetic and fluid simulations illustrating the physics of magnetic reconnection will be presented. We will dem~:tistrate that kinetic nongyrotropic effects can significantly alter the global magnetosphere evolution and location of reconnection sites.

Kuznetsova, M. M.↗

Electromagnetic Scattering by Spheroidal Volumes of Discrete Random Medium

We use the superposition T-matrix method to compare the far-field scattering matrices generated by spheroidal and spherical volumes of discrete random medium having the same volume and populated by identical spherical particles. Our results fully confirm the robustness of the previously identified coherent and diffuse scattering regimes and associated optical phenomena exhibited by spherical particulate volumes and support their explanation in terms of the interference phenomenon coupled with the order-of-scattering expansion of the far-field Foldy equations. We also show that increasing non-sphericity of particulate volumes causes discernible (albeit less pronounced) optical effects in forward and backscattering directions and explain them in terms of the same interference/multiple-scattering phenomenon.

Electromagnetic scattering; Multi-particle groups;↗

A tutorial solution to scattering of radiation in a thin atmosphere bounded below by a diffusely reflecting, absorbing surface

A simple tutorial method, based on a photon tracking procedure, is described to determine the spherical albedo for a thin atmosphere overlying a reflecting surface. This procedure is used to provide a physical structure with which to interpret the more detailed but highly mathematical analyses presented. The final equations are shown to be in good numerical agreement with more exact solutions for thin atmospheres.

Buglia, J. J.↗

Numerical Study of Buoyancy and Different Diffusion Effects on the Structure and Dynamics of Triple Flames

Numerical simulations of 2-D triple flames under gravity force have been implemented to identify the effects of gravity on triple flame structure and propagation properties and to understand the mechanisms of instabilities resulting from both heat release and buoyancy effects. A wide range of gravity conditions, heat release, and mixing widths for a scalar mixing layer are computed for downward-propagating (in the same direction with the gravity vector) and upward-propagating (in the opposite direction of the gravity vector) triple flames. Results of numerical simulations show that gravity strongly affects the triple flame speed through its contribution to the overall flow field. A simple analytical model for the triple flame speed, which accounts for both buoyancy and heat release, is developed. Comparisons of the proposed model with the numerical results for a wide range of gravity, heat release and mixing width conditions, yield very good agreement. The analysis shows that under neutral diffusion, downward propagation reduces the triple flame speed, while upward propagation enhances it. For the former condition, a critical Froude number may be evaluated, which corresponds to a vanishing triple flame speed. Downward-propagating triple flames at relatively strong gravity effects have exhibited instabilities. These instabilities are generated without any artificial forcing of the flow. Instead disturbances are initiated by minute round-off errors in the numerical simulations, and subsequently amplified by instabilities. A linear stability analysis on mean profiles of stable triple flame configurations have been performed to identify the most amplified frequency in spatially developed flows. The eigenfunction equations obtained from the linearized disturbance equations are solved using the shooting method. The linear stability analysis yields reasonably good agreements with the observed frequencies of the unstable triple flames. The frequencies and amplitudes of disturbances increase with the magnitude of the gravity vector. Moreover, disturbances appear to be most amplified just downstream of the premixed branches. The effects of mixing width and differential diffusion are investigated and their roles on the flame stability are studied.

Chen, Jyh-Yuan↗

Reformulation of quasi-linear theory.

Standard plasma quasi-linear theory is reformulated on the basis of a classical quantum derivation proceeding from the Vlasov equation and dealing only with frequency, wavenumber, and velocity. The wave amplitudes are assumed to be weakly time-dependent, and no distinction is made between growing and decaying waves. The proposed method leads to no negative diffusivity of 'fake' diffusion. By appropraite treatment of nonresonant interaction, expressions are obtained for wave energy and momentum.

Kaufman, A. N.↗

Langevin equation with stochastic damping - Possible application to critical binary fluid

We solve the familiar Langevin equation with stochastic damping to represent the motion of a Brownian particle in a fluctuating medium. A connection between the damping and the random driving forces is proposed which preserves quite generally the Einstein relation between the diffusion and mobility coefficients. We present an application to the case of a Brownian particle in a critical binary mixture.

Jasnow, D.↗

Verification of reflectance models in turbid waters

Inherent optical parameters of very turbid waters were used to evaluate existing water reflectance models. Measured upwelling radiance spectra and Monte Carlo simulations of the radiative transfer equations were compared with results from models based upon two flow, quasi-single scattering, augmented isotropic scattering, and power series approximation. Each model was evaluated for three separate components of upwelling radiance: (1) direct sunlight; (2) diffuse skylight; and (3) internally reflected light. Limitations of existing water reflectance models as applied to turbid waters and possible applications to the extraction of water constituent information are discussed.

Tanis, F. J.↗

Numerical analysis of confined recirculating gas-solid turbulent flows

Particle dispersion in confined recirculating turbulent flows has been investigated numerically. The present computational scheme utilizes Eulerian forms of the governing equations and allows two-way coupling between gas and solid phases. A recently developed two-phase closure model based on the multiple-scale turbulence model was used for the estimation of turbulent viscosities and diffusivities. For the particle size and loading considered in this study, the effect of particulate phase on the conveying gas is small, the nonequilibrium between the two phases is significant. Effects of recirculation, expanded chamber size and secondary annular jet momentum on the particle mixing rate are also investigated. In general, the present numerical results are in reasonably good agreement with the available experimental data.

Chen, C. P.↗

Experimental thermal conductivity, thermal diffusivity, and specific heat values for mixtures of nitrogen, oxygen, and argon

Experimental measurements of thermal conductivity and thermal diffusivity obtained with a transient hot-wire apparatus are reported for three mixtures of nitrogen, oxygen, and argon. Values of the specific heat, Cp, are calculated from these measured values and the density calculated with an equation of state. The measurements were made at temperatures between 65 and 303 K with pressures between 0.1 and 70 MPa. The data cover the vapor, liquid, and supercritical gas phases for the three mixtures. The total reported points are 1066 for the air mixture (78.11 percent nitrogen, 20.97 percent oxygen, and 0.92 percent argon), 1058 for the 50 percent nitrogen, 50 percent oxygen mixture, and 864 for the 25 percent nitrogen, 75 oxygen mixture. Empirical thermal conductivity correlations are provided for the three mixtures.

Perkins, R. A.↗

Diffusion Of Mass In Evaporating Multicomponent Drops

Report summarizes study of diffusion of mass and related phenomena occurring in evaporation of dense and dilute clusters of drops of multicomponent liquids intended to represent fuels as oil, kerosene, and gasoline. Cluster represented by simplified mathematical model, including global conservation equations for entire cluster and conditions on boundary between cluster and ambient gas. Differential equations of model integrated numerically. One of series of reports by same authors discussing evaporation and combustion of sprayed liquid fuels.

Bellan, Josette↗

Evolution of dynamo-generated magnetic fields in accretion disks around compact and young stars

Geometrically thin, optically thick, turbulent accretion disks are believed to surround many stars. Some of them are the compact components of close binaries, while the others are throught to be T Tauri stars. These accretion disks must be magnetized objects because the accreted matter, whether it comes from the companion star (binaries) or from a collapsing molecular cloud core (single young stars), carries an embedded magnetic field. In addition, most accretion disks are hot and turbulent, thus meeting the condition for the MHD turbulent dynamo to maintain and amplify any seed field magnetic field. In fact, for a disk's magnetic field to persist long enough in comparison with the disk viscous time it must be contemporaneously regenerated because the characteristic diffusion time of a magnetic field is typically much shorter than a disk's viscous time. This is true for most thin accretion disks. Consequently, studying magentic fields in thin disks is usually synonymous with studying magnetic dynamos, a fact that is not commonly recognized in the literature. Progress in studying the structure of many accretion disks was achieved mainly because most disks can be regarded as two-dimensional flows in which vertical and radial structures are largely decoupled. By analogy, in a thin disk, one may expect that vertical and radial structures of the magnetic field are decoupled because the magnetic field diffuses more rapidly to the vertical boundary of the disk than along the radius. Thus, an asymptotic method, called an adiabatic approximation, can be applied to accretion disk dynamo. We can represent the solution to the dynamo equation in the form B = Q(r)b(r,z), where Q(r) describes the field distribution along the radius, while the field distribution across the disk is included in the vector function b, which parametrically depends on r and is normalized by the condition max (b(z)) = 1. The field distribution across the disk is established rapidly, while the radial distribution Q(r) evolves on a considerably longer timescale. It is this evolution that is the subject of this paper.

Stepinski, Tomasz F.↗