Engineering PapersSearch

SEARCH · Engineering Papers

Results for “anisotropic diffusion”

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 37 records · Page 2

Model atmospheres and radiation of magnetic neutron stars. I - The fully ionized case

Model neutron star atmospheres are calculated for typical cooling stars with a strong magnetic field and effective temperatures of 10 exp 5 to 10 exp 6 K. The effect of anisotropic photon diffusion in two normal modes are examined under the assumption that the opacity is due solely to the bremsstrahlung and Thomson scattering processes under conditions of LTE that are expected to prevail at the temperatures and densities obtained. The main aspects of anisotropic photon diffusion, and an original procedure for calculating model atmospheres and emitted spectra are discussed. Representative calculated spectra are given, and it is found that the hard spectral excess characterizing the nonmagnetic case, while still present, becomes less prominent in the presence of magnetic fields in the range of 10 exp 11 to 10 exp 13 G.

Shibanov, Iu. A.

Mesoscale convection in the clouds of Venus

A theory explaining the high aspect ratios for cloud level convection in the Venusian atmosphere is presented. The apparent flatness of large-scale convection cells in the Venusian clouds is accounted for by anisotropic eddy diffusion and radiative transfer effects. Horizontal eddy diffusivities must be at least ten-fold greater than vertical diffusivities. Anisotropy ratios greater than or approximately equal to 1,000 are sufficient to explain the flattening in cases where the vertical eddy diffusivity is greater than or approximately equal to 100,000 sq cm/s and the effects of radiative transfer are negligible. For vertical diffusivity greater than or approximately equal to 100,000 sq cm/s, radiative transfer contributes to the flattening. Radiative transfer alone can not account for the apparent aspect ratios, as this would require vertical heat diffusivities large enough to give a nonphysically high value of convective heat transport.

Covey, C. C.

A three dimensional model of the Venusian thermosphere with superrotation

An improved three dimensional spectral model of the thermosphere of Venus is described. The model solves the Navier-Stokes equations and includes nonlinear effects for an arbitrary number of atmospheric species. A two dimensional axisymmetric model of the superrotation of the thermosphere is also presented. This model addresses the Pioneer-Venus mission finding, which suggested the thermospheric rotation rate to be much higher than that of the planet as seen from the asymmetric distribution of hydrogen and helium. Both models include the effects of an anisotropic eddy diffusion that is consistent with atmospheric mixing length theory.

Stevens-Rayburn, D. R.

Modeling scalar flux and the energy and dissipation equations

Closure models derived from the Two-Scale Direct-Interaction Approximation were compared with data from direct simulations of turbulence. Attention was restricted to anisotropic scalar diffusion models, models for the energy dissipation equation, and models for energy diffusion.

Yoshizawa, A.

Interpretation and modeling of the averaged equations for a fluid-solid flow

A self-consistent derivation of the conservation laws is given for flows of a fluid-solid mixture. A unified analytical framework for obtaining constitutive relations is provided. This analysis uses a control volume/control surface approach that is widely used in fluid mechanics. All terms in the governing equations and the constitutive relations are written in terms of the mass-weighted averages except solid concentration. It is believed that the mass-weighted average is the natural bridge between micromechanics and constitutive relations. The derived momentum equations contain terms that differ from all existing models except that of Prosperetti and Jones (1984). However, their assumptions are not needed here. Special attention is given to the solid phase pressure. The physical basis of the previously assumed form for this pressure (Givler 1987) becomes clear. A number of related phenomena are also discussed. These include the anti-diffusion and anisotropic normal stresses. The energy equations are also different from existing models.

Shen, Hayley H.

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam

Diffuse X-rays from the galactic disk

An anisotropic feature of the diffuse hard X-ray background that tracks the concentration of interstellar hydrogen in the plane of the galaxy is reported. This feature supports a model of galactic X-ray emission by subrelativistic cosmic rays via a bremsstrahlung process. The measurement was carried out on August 9, 1971, using two multianode multilayer gas proportional counters onboard Aerobee 170 flight 13.08. A schematic diagram of the detectors is shown. This type of construction and the appropriate utilization of the signals from the many anodes result in a low detector background, a perequisite before undertaking a measurement of possible small variations in the brightness of the X-ray sky.

Serlemitsos, P. J.

Remote measurement of in-plane diffusivity components in plates

A method of determining thermal diffusivity in thin plates is presented. The method, using infrared images of evolving thermal patterns previously injected with a laser, is noncontacting, one-sided, and remote. It does not require independent estimates of either the emissivity of the sample or the sample thickness. With a line-segment pattern for thermal input, it yields the inplane components of the diffusivity tensor in anisotropic materials and also the rate of heat loss to the environment of the plate. Two methods of data analysis are presented, one corresponding to a heating line of general cross section and the other considering a Gaussian cross section, thereby saving considerable computer time. Both methods produce a statistical evaluation of measurement quality as well as estimates of diffusivity and loss rate. Results are shown for plates of metals and graphite-epoxy composite materials. Principal components and orientation for the diffusivity tensor are obtained in the anisotropic graphite-epoxy sample.

Welch, Christopher S.

THEMIS Observations of the Magnetopause Electron Diffusion Region: Large Amplitude Waves and Heated Electrons

We present the first observations of large amplitude waves in a well-defined electron diffusion region based on the criteria described by Scudder et al at the subsolar magnetopause using data from one Time History of Events and Macroscale Interactions during Substorms (THEMIS) satellite. These waves identified as whistler mode waves, electrostatic solitary waves, lower hybrid waves, and electrostatic electron cyclotron waves, are observed in the same 12 s waveform capture and in association with signatures of active magnetic reconnection. The large amplitude waves in the electron diffusion region are coincident with abrupt increases in electron parallel temperature suggesting strong wave heating. The whistler mode waves, which are at the electron scale and which enable us to probe electron dynamics in the diffusion region were analyzed in detail. The energetic electrons (approx. 30 keV) within the electron diffusion region have anisotropic distributions with T(sub e(right angle))/T(sub e(parallel)) > 1 that may provide the free energy for the whistler mode waves. The energetic anisotropic electrons may be produced during the reconnection process. The whistler mode waves propagate away from the center of the "X-line" along magnetic field lines, suggesting that the electron diffusion region is a possible source region of the whistler mode waves.

magnetopause and boundary layers

Self-Diffusion of Drops in a Dilute Sheared Emulsion

Self-diffusion coefficients that describe cross-flow migration of non-Brownian drops in a dilute sheared emulsion were obtained by trajectory calculations. A boundary integral formulation was used to describe pairwise interactions between deformable drops; interactions between undeformed drops were described with mobility functions for spherical drops. The results indicate that drops have large anisotropic self-diffusivities which depend strongly on the drop viscosity and modestly on the shear-rate. Pairwise interactions between drops in shear-flow do not appreciably promote drop breakup.

Loewenberg, Michael

Morphological evolution of growing crystals - A Monte Carlo simulation

The combined effects of nutrient diffusion and surface kinetics on the crystal morphology were investigated using a Monte Carlo model to simulate the evolving morphology of a crystal growing from a two-component gaseous nutrient phase. The model combines nutrient diffusion, based on a modified diffusion-limited aggregation process, with anisotropic surface-attachment kinetics and surface diffusion. A variety of conditions, ranging from kinetic-controlled to diffusion-controlled growth, were examined. Successive transitions from compact faceted (dominant surface kinetics) to open dendritic morphologies (dominant volume diffusion) were obtained.

Xiao, Rong-Fu

Short term prediction of E greater than or equal to 10 MeV proton fluxes from solar flares

Both the anisotropic and isotropic diffusion theories can be used to extrapolate proton fluxes for E greater than or equal to 10 meV for over 50% of the particle events. The isotropic diffusion theory uses a diffusion coefficient: D = Mr sup beta. It was found that M and beta tended to be functions of flare position on the solar disk. A measurement of the interplanetary flux in near earth space gives a good indication of the polar cap fluxes. It was found that the 30 MHz absorption over the poles during a PCA is proportional to the square root of the integral proton flux E greater than or equal to 11 meV in interplanetary space, J = KA squared, with K = 8 plus or minus 2 and J in protons/sq cm-sec-ster.

Kuck, G. A.

Solar particles /observations, relationship to the sun acceleration, interplanetary medium/

The major features of the propagation of flare particles in the interplanetary medium are discussed in terms of the classic well-behaved flare having unique impulsive injection and a smooth time profile. Topics include flare events, their frequency of occurrence, development of a typical event, energy spectra, proton and electron types, charge and isotopic composition, solar flares and particle accelerations, radio and X-ray observations, the Fermi mechanism, the betatron mechanism, acceleration models, plasma instabilities, two-stage acceleration, propagation mechanisms, the anisotropic stage, the diffusive stage, and the convection and energy loss stage.

Mcdonald, F. B.

Pitch Angle Distributions of Solar Energetic Particles and the Local Scattering Properties of the Interplanetary Magnetic Field

An approximate solution of the Fokker-Planck equation containing pitch angle scattering and adiabatic focusing is discussed. For modest focusing effects the omnidirectional density obeys an ordinary diffusion equation with a modified diffusion coefficient. The anisotropic part of the distribution function is properly normalized and split into an even and an odd part. The even part is determined by the ratio between the scattering mean free path and the focusing length and by the degree of polarization of the magnetic field fluctuations. The odd part is determined by the deviation of the pitch angle scattering from isotropic scattering. The method supplies a powerful tool to obtain the local characteristics of pitch angle scattering. It is insensitive to long lasting solar injections and to moderate radial variations of the mean free path. The method is applied to solar particle events observed on Helio-1 and -2.

Beeck, J.

Inertial Effects in Suspension Dynamics

The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

J. F. Brady

Thermophysical Property Measurements of Indium Iodide Crystals

InI single crystals are a promising room temperature detector material for X-rays and γ-rays. To improve crystal growth of the material by simulations, knowledge of thermophysical properties is essential, and since InI is orthorhombic, the anisotropy has to be taken into account. The temperature dependence of several thermophysical properties have been measured for InI, including the anisotropic thermal expansion, specific heat, and the thermal diffusivity in the b direction. The anisotropic thermal expansion coefficients, determined by X-ray diffraction, were α11=1.03·10-5 K-1, α22= 3.77·10-5 K-1, and α33 = 6.26·10-5 K-1. The specific heat, measured by DSC, was 0.226 J g-1 K-1 at 335K, with a temperature dependence of 9.582·10-5 J g-1 K-2. In the course of the X-ray diffraction and DSC measurements, it was shown that supposed phase changes, reported in older literature, are actually not phase changes but oxidation effects. The thermal diffusivity in the b direction, measured by the Xenon Flash method, was 0.288·10-6 m2s-1 at RT, decreasing to 0.253·10-6 m2 s-1 at 450K. In addition, the volume increase upon melting and the thermal expansion of the melt have been determined.

Martin Volz

Mass loading and velocity diffusion models for heavy pickup ions at comet Grigg-Skjellerup

We compare model predictions of cometary water group ion densities and the solar wind slow down with measurements made by the Giotto Johnstone plasma analyzer implanted ion sensor at the encounter with comet Grigg-Skjellerup (G-S) on July 10, 1992. The observed slope of the ion density profile on approach to the comet is unexpectedly steep. Possible explanations for this are discussed. We present also a preliminary investigation of the quasilinear velocity-space diffusion of the implanted heavy ion population at G-S using a transport equation including souce, convection, adiabatic compression, and velocity diffusion terms. Resulting distributions are anisotropic, in agreement with observations. We consider theoretically the waves that may be generated by the diffusion process for the observed solar wind conditions. At initial ion injections, waves are generated at omega approximately Omega(sub i) the ion gyrofrequency, and lower frequencies are predicted for diffusion toward a bispherical shell.

Huddleston, D. E.