Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion coefficient”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Study of cosmic-ray H and He isotopes at 23 AU

We have measured the spectra of H and He isotopes during the 1987 solar minimum with the cosmic-ray detector system (CRS) on the Voyager 2 spacecraft. By carrying out the measurement near solar minimum and at large heliospheric distances, the effects of solar modulations were reduced. In particular, the adiabatic energy losses were smaller, and these results from 23 AU over the solar minimum period of cycle 21 represent observations at energies not accessible from previous measurements near 1 AU. The modulated spectra with the diffusion coefficient constant k(sub 0) = 3.15 x 10(exp 22) sq cm/s (which corresponds to a solar modulation parameter of 360 MV at 23 AU and 500 MV at 1 AU) agree well with both our data at 23 AU and the previous solar minimum measurements at 1 AU. The measured H-1 and H-2 spectra are both consistent with the calculated spectra, using standard Galactic and heliospheric propagation models without invoking an anomalous hydrogen component. With the fixed modulation parameter of 360 MV, the mean pathlengths, source spectra, and cross sections were varied to study the effects of different input parameters on the spectra and relative abundances. At this stage of our work, we have not found any strong evidence from the low-energy H-2 and He-3 data that H-1 and He-4 should have a different propagation history, or different types of source spectra from the heavier cosmic-ray nuclei.

Seo, E. S.↗
NASA NTRS1993

A two-dimensional model with coupled dynamics, radiation, and photochemistry. 1: Simulation of the middle atmosphere

Spatial and seasonal distribution of ozone and other trace gases are simulated using a interactive two-dimensional model of the stratosphere updated to include full gas phase chemistry. The model consists of a primitive equation dynamics module, a full radiative transfer scheme, and a comprehensive gas phase chemistry module. The circulation is derived from heating rates in the stratosphere that are calculated using model-generated ozone. In the troposphere, parameterized heating rates are adopted. The eddy momentum flux divergence in the zonal mean momentum equation is given by the eddy fluxes of potential vorticity. Eddy fluxes of potential vorticity and tracers are parameterized using a set of predetermined diffusion coefficients. The adopted values for K(sub yy), show a hemispherical asymmetry in that the values in the lower stratosphere are consistently smaller in the southern hemisphere. The asymmetry in K(sub yy) and in the parameterization of the tropospheric heating rate results in an asymmetry in the circulation giving rise to unique signatures in the trace gas distributions. The model successfully simulates the observed asymmetry in the column abundance of the springtime ozone maxima between northern and southern hemisphere. Results for other trace gases are in agreement with the gross observed features although specific differences exist.

Ko, Malcolm K. W.↗
NASA NTRS1994

The Balmer-dominated northeast limb of the Cygnus loop supernova remnant

We present a comprehensive investigation of the Balmar-dominated northeast limb of the Cygnus Loop supernova remnant. Data presented include H alpha (O III), and X-ray images, UV and visible spectrophotometry, and high-resolution spectroscophy. The two relatively bright Balmer-dominated filaments visible on the POSS prints are seen to be part of a very smooth and regular complex of filaments. These filaments mark the current location of the blast wave and are seen to bound the sharply limb-brightened X-ray emission, including the previously reported X-ray, 'halo.' The (O III)/h beta ratio throughout the region is approximately 0.1, except for regions in which the shock is undergoing a transition from nonradiative to incomplete radiative to incomplete radiative conditions. At these locations (O III) emission from the cooling region is quite strong, while collisionally excited Balmer-line emission can be weak because of photoionization of the preshock medium by UV from the nascent cooling region. As a result (O III)/H beta is greater than 100 in some locations. The nonradiative/radiative transition is best studied along the length of the northwestern of the two brightest filaments, where the shock velocity and swept-up column go from approximately 180 km/s and 10(exp 17)/sq cm at one end to approximately 140 km/s and 8 x 10(exp 17)/cm at the other. There are also a number of locations of such incomplete radiative emission where the shock has recently encountered denser regions with characteristic sizes of approximately 10(exp 18) cm. There is a considerable amount of evidence that the shock has decelerated from approximately 400 km/s to less than 200 km/s in the last 1000 yr. We interpret this as the result of the blast wave hitting the wall of a cavity which surround the supernova precursor and succeed in matching a wide range of data with a reflected shock model in which the density ofthe cavity wall is approximately 1.2/cu cm and the density in the interior of the cavity is about 0.09/cu cm. This analysis supports previous suggestions that the Cygnus Loop was a cavity explosion. We opredict that the X-ray shell will continue to brighten for the next 1000 years or so, during which time the shock front will become radiative throughout the northeast region. The similarity of the time scales for post-shock cooling, shock deceleration, and the 1ressure evolution of the region violates assumptions of steady flow, and greatly complicated interpretation of shock specptra in the region. This may be a common and important effect in many supernova remnants. High-resolutoin spectroscopy of the H alpha line profile shows that the broad component is surprisingly narrow (approximately 130 km/s). Consistency between this line width and the shock velocity inferred from spectrophotometry requires either rapid postshock equilibration of T(sub i) and T(sub e), or that the shock front has declerated from 180 km/s to 150 km/s in the last approximately 200 yr. While we cannot unambiguously distinquish between these models, the required deceleration is consistent with extrapolation of the shock decleration rate inferred over the last 1000 years. The width of the narrow component of the h alpha line profile is approximately 3 km/s, which indicates that the preshock medium has been heated to approximately 25,000 K by a shock precursor. The most promising types of precursor are heating either by fast neutrals which overtake the shock or by cosmic rays. In the latter case the observed properties of the precursor constrain the cosmic-ray diffusion coefficient to be between about 10(exp 23) sq cm/s. The characteristics of the precursor seem to be about right to explain the presence of weak (S II) and (N II) emission previously reported by Fesen & Itoh (1985).

Hester, J. Jeff↗
NASA NTRS1994

Assimilation by lunar mare basalts: Melting of crustal material and dissolution of anorthite

We discuss techniques for calculating the amount of crustal assimilation possible in lunar magma chambers and dikes based on thermal energy balances, kinetic rates, and simple fluid mechanical constraints. Assuming parent magmas of picritic compositions, we demonstrate the limits on the capacity of such magmas to melt and dissolve wall rock of anorthitic, troctolitic, noritic, and KREEP (quartz monzodiorite) compositions. Significant melting of the plagioclase-rich crustal lithologies requires turbulent convection in the assimilating magma and an efficient method of mixing in the relatively buoyant and viscous new melt. Even when this occurs, the major element chemistry of the picritic magmas will change by less than 1-2 wt %. Diffusion coefficients measured for Al2O3 from an iron-free basalt and an orange glass composition are 10(exp -12) sq m/s at 1340 C and 10(exp -11) sq m/s at 1390 C. These rates are too slow to allow dissolution of plagioclase to significantly affect magma compositions. Picritic magmas can melt significant quantities of KREEP, which suggests that their trace element chemistry may still be affected by assimilation processes; however, mixing viscous melts of KREEP composition with the fluid picritic magmas could be prohibitively difficult. We conclude that only a small part of the total major element chemical variation in the mare basalt and volcanic glass collection is due to assimilation/fractional crystallization processes near the lunar surface. Instead, most of the chemical variation in the lunar basalts and volcanic glasses must result from assimilation at deeper levels or from having distinct source regions in a heterogeneous lunar mantle.

Finnila, A. B.↗
NASA NTRS1994

Vertical transport rates in the statosphere in 1993 from observations of CO2, N2O, and CH4

Measurements of CO2, N2O, and CH4 are analyzed to define hemispheric average vertical exchange rates in the lower stratosphere from November 1992 to October 1993. Effective vertical diffusion coefficients were small in summer, less than or equal to 1 sq m/s at altitudes below 25 km; values were similar near the tropopause in winter, but increased markedly with altitude. The analysis suggests possible longer residence times for exhaust from stratospheric aircraft, and more efficient transport from 20 km to the middle stratosphere, than predicted by many current models. Seasonally-resolved measurements of stratospheric CO2 and N2O provide significant new constraints on rates for global-scale vertical transport.

Wofsy, Steven C.↗
NASA NTRS1995

Time-dependent simulation of oblique MHD cosmic-ray shocks using the two-fluid model

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic-ray (CR)-modified plane shock evolution. Most of the calculations are done with a two-fluid model for diffusive shock acceleration, but we provide also comparisons between a typical shock computed that way against calculations carried out using the more complete, momentum-dependent, diffusion-advection equation. We also illustrate a test showing that these simulations evolve to dynamical equilibria consistent with previously published steady state analytic calculations for such shocks. In order to improve understanding of the dynamical role of magnetic fields in shocks modified by CR pressure we have explored for time asymptotic states the parameter space of upstream fast mode Mach number, M(sub f), and plasma beta. We compile the results into maps of dynamical steady state CR acceleration efficiency, epsilon(sub c). We have run simulations using constant, and nonisotropic, obliquity (and hence spatially) dependent forms of the diffusion coefficient kappa. Comparison of the results shows that while the final steady states achieved are the same in each case, the history of CR-MHD shocks can be strongly modified by variations in kappa and, therefore, in the acceleration timescale. Also, the coupling of CR and MHD in low beta, oblique shocks substantially influences the transient density spike that forms in strongly CR-modified shocks. We find that inside the density spike a MHD slow mode wave can be generated that eventually steepens into a shock. A strong layer develops within the density spike, driven by MHD stresses. We conjecture that currents in the shear layer could, in nonplanar flows, results in enhanced particle accretion through drift acceleration.

Frank, Adam↗
NASA NTRS1993

Large-scale isentropic mixing properties of the Antarctic polar vortex from analyzed winds

Winds derived from analyzed geopotential height fields are used to study quasi-horizontal mixing by the large-scale flow in the lower stratosphere during austral spring. This is the period when the Antarctic ozone hole appears and disappears. Trajectories are computed for large ensembles of particles initially inside and outside the main polar vortex. Mixing and transport are diagnosed through estimates of finite time Lyapunov exponents and Lagrangian dispersion statistics of the tracer trajectories. At 450 K and above prior to the vortex breakdown: Lyapunov exponents are a factor of 2 smaller inside the vortex than outside; diffusion coefficients are an order of magnitude smaller inside than outside the vortex; and the trajectories reveal little exchange of air across the vortex boundary. At lower levels (425 and 400 K) mixing is greater, and there is substantial exchange of air across the vortex boundary. In some years there are large wave events that expel small amounts of vortex air into the mid-latitudes. At the end of the spring season during the vortex breakdown there is rapid mixing of air across the vortex boundary, which is evident in the mixing diagnostics and the tracer trajectories.

Bowman, Kenneth P.↗
NASA NTRS1994

Focused interplanetary transport of approximately 1 MeV solar energetic protons through self-generated Alfven waves

We present a model of the focused transport of approximately 1 MeV solar energetic protons through interplanetary Alfven waves that the protons themselves amplify or damp. It is based on the quasi-linear theory but with a phenomenological pitch angle diffusion coefficient in the 'resonance gap.' For initial Alfven wave distributions that give mean free paths greater than approximately 0.5 AU for approximately 1 MeV protons in the inner heliosphere, the model predicts greater than roughly an order of magnitude amplification (damping) in the outward (inward) propagating resonant Alfven waves at less than or approximately equal to o.3 AU heliocentric distance. As the strength of proton source is increased, the peak differential proton intensity at approximately 1 MeV at 1 AU increases to a maximum of approximately 250 particles (/(sq cm)(s)(sr)(MeV)) and then decreases slowly. It may be attenuated by a factor of 5 or more relative to the case without wave evolution, provided that the proton source is sufficiently intense that the resulting peak differential intensity of approximately 1 MeV protons at 1 AU exceeds approximately 200 particles (/(sq cm)(s)(sr)(MeV)). Therefore, in large solar proton events, (1) one may have to take into account self-amplified waves in studying solar particle propagation, (2) the number of accelerated protons escaping from a flare or interplanetary shock may have been underestimated in past studies by a significant factor, and (3) accelerated protons escaping from a traveling interplanetary shock at r less than or approximately equal to 0.3 AU should amplify the ambient hydromagnetic waves siginificantly to make the shock an efficient accelerator, even if initially the mean free path is greater than or approximately equal to 1 AU.

Ng, C. K.↗
NASA NTRS1994

Kinetic Alfven waves as a source of plasma transport at the dayside magnetopause

As the shocked solar wind with variable plasma density and magnetic field impinges on the dayside magnetopause, it is likely to generate large-scale Alfven waves at the solar wind magnetosphere interface. However, large gradients in the density and magnetic field at the magnetopause boundary effectively couple large-scale Alfven waves with kinetic Alfven waves may play an important role in plasma transport at the dayside magnetopause and in electron acceleration along field lines. The transport can occur because, unlike the magnetohydrodynamic (MHD) shear Alfven wave, the kinetic Alfven wave has an associated parallel electric field which breaks down the 'frozen-in' condition and decouples the plasma from field lines. We calculate the average deviation of the plasma from the field line from which we estimate the diffusion coefficient associated with these 'bundles' of lead to acceleration of electrons along field lines in the magnetopause boundary and may possibly provide an explanation for observed counterstreaming electron beams characterized by energies of 50-200 eV.

Lee, L. C.↗
NASA NTRS1995

Photochemical modeling of Titan's atmosphere

We have developed a new photochemical model of Titan's atmosphere which includes all the important compounds and reactions in spherical geometry from the surface to 1240 km. Compared to the previous model of Yung et al. (1984), the most significant recent change in the reactions used is the updated methane dissociation scheme (Mordaunt et al. 1993). Moreover, the transfer of the solar radiation in the atmosphere and the photolysis rates have been calculated by using a Monte Carlo code. Finally, the eddy diffusion coefficient profile is adjusted in order to fit the mean vertical distribution of HCN retrieved from millimeter groundbased observations of Tanguy et al. (1990); using new values for the boundary flux of atomic nitrogen (Strobel et al. 1992). We have run the model in both steady-state and diurnal modes, with 62 species involved in 249 reactions. There is little difference between diurnal and steady-state results. Overall our results are in a closer agreement with the abundances inferred from the Voyager infrared measurements at the equator than the Yung et al. results. We find that the catalytic scheme for H recombination invoked by Yung et al. only slightly improves the model results and we conclude that this scheme is not essential to fit observations.

Toublanc, D.↗
NASA NTRS1995

Improved Optics For Quasi-Elastic Light Scattering

Improved optical train devised for use in light-scattering measurements of quasi-elastic light scattering (QELS) and laser spectroscopy. Measurements performed on solutions, microemulsions, micellular solutions, and colloidal dispersions. Simultaneous measurements of total intensity and fluctuations in total intensity of light scattered from sample at various angles provides data used, in conjunction with diffusion coefficients, to compute sizes of particles in sample.

Cheung, Harry Michael↗
NASA NTRS1994

Energetic electrons in the magnetosphere of Saturn

The energy spectra and angular distributions of electrons observed by Pioneer 11 as a function of radial distance in the inner magnetosphere of Saturn are reanalyzed and phase space densities are then calculated. The radial dependence of phase space density requires a distributed loss process. The loss is greatest in the region of the E ring (5.5 less than L less than 8.5) and is attributed to collisions with the ring particles in agreement with earlier work by Van Allen et al. (1980). Quantitative analysis yields the following properties of the E ring: the particle radii are in the range of 4 x 10(exp -5) to 3.2 x 10(exp -4) cm and the thickness of the ring is approximately 3 R(sub s). Between the inner edge of the E ring (5.5 R(sub s)) and the outer edge of the A ring (2.3 R(sub s)) there are more energetic electrons than can be supplied by radial diffusion from an external source. Detailed calculations show that a cosmic ray albedo neutron decay (CRAND) source in the A and B rings is a plausible source for this excess. The radial diffusion coefficient required to explain the E ring absorption and CRAND source for electrons is 1 x 10(exp -12) greater than D(sub 0) greater than 3 x 10(exp -12) R(exp 2, sub s)/s, assuming that D(sub LL) = D(sub 0)L(exp 3). As part of the reanalysis program, a method for the deconvolution of pitch angle distributions observed by simple detectors on a rotating spacecraft is developed. This process removes the instrumental response and rotational smear due to finite sampling periods and yields true angular distributions.

Randall, B. A.↗
NASA NTRS1995

The observed North-South Asymmetry of IMF spiral

We appraise the finding, reported in the literature, that a small but finite north-south asymmetry (NSA) exists in the interplanetary magnetic field (IMF) spiral at Earth's orbit. We have analyzed the data available on the Omnitape for the 1963 to 1993 period. The coverage is very uneven, ranging from less than 40% to greater than 80%. The magnitude of NSA fluctuates considerably during the period of our analysis. This is true even if one considers the period 1967 to 1982 when the coverage is greater than 50%. The values of NSA derived from 27-day averages of the hourly data points range from greater than +50 deg to less than -40 deg. If one arranges the data according to the magnetic polarity epochs of the solar polar field, the epoch averages gives the magnitude of NSA less than approximately 2 deg. This is also true, if one considers the average magnitude of NSA for the 1965 to 1993 period, when the coverage is greater than 25%. A genuine, persistent, NSA of IMF spiral is likely to affect the cosmic ray modulation, on either side of the current sheet, by introducing a corresponding change in the radial diffusion coefficient of energetic particle transport in the heliosphere. The annual mean values of the observed NSA of IMF spiral are compared with the observed off-ecliptic contributions to cosmic ray modulation.

Ahluwalia, H. S.↗
NASA NTRS1996

A Note on Substructuring Preconditioning for Nonconforming Finite Element Approximations of Second Order Elliptic Problems

In this paper an algebraic substructuring preconditioner is considered for nonconforming finite element approximations of second order elliptic problems in 3D domains with a piecewise constant diffusion coefficient. Using a substructuring idea and a block Gauss elimination, part of the unknowns is eliminated and the Schur complement obtained is preconditioned by a spectrally equivalent very sparse matrix. In the case of quasiuniform tetrahedral mesh an appropriate algebraic multigrid solver can be used to solve the problem with this matrix. Explicit estimates of condition numbers and implementation algorithms are established for the constructed preconditioner. It is shown that the condition number of the preconditioned matrix does not depend on either the mesh step size or the jump of the coefficient. Finally, numerical experiments are presented to illustrate the theory being developed.

Maliassov, Serguei↗
NASA NTRS1994

Assimilation by Lunar Mare Basalts: Melting of Crustal Material and Dissolution of Anorthite

We discuss techniques for calculating the amount of crustal assimilation possible in lunar magma chambers and dikes based on thermal energy balances, kinetic rates, and simple fluid mechanical constraints. Assuming parent magmas of picritic compositions, we demonstrate the limits on the capacity of such magmas to melt and dissolve wall rock of anorthitic, troctolitic, noritic, and KREEP (quartz monzodiorite) compositions. Significant melting of the plagioclase-rich crustal lithologies requires turbulent convection in the assimilating magma and an efficient method of mixing in the relatively buoyant and viscous new melt. Even when this occurs, the major element chemistry of the picritic magmas will change by less than 1-2 wt %. Diffusion coefficients measured for Al2O3 from an iron-free basalt and an orange glass composition are 10(exp -12) m(exp 2) s(exp -1) at 1340 C and 10(exp -11) m(exp 2) s(exp -1) at 1390 C. These rates are too slow to allow dissolution of plagioclase to significantly affect magma compositions. Picritic magmas can melt significant quantities of KREEP, which suggests that their trace element chemistry may still be affected by assimilation processes; however, mixing viscous melts of KREEP composition with the fluid picritic magmas could be prohibitively difficult. We conclude that only a small part of the total major element chemical variation in the mare basalt and volcanic glass collection is due to assimilation/fractional crystallization processes near the lunar surface. Instead, most of the chemical variation in the lunar basalts and volcanic glasses must result from assimilation at deeper levels or from having distinct source regions in a heterogeneous lunar mantle.

Finnila, A. B.↗
NASA NTRS1994

Mathematical Model for the Mineralization of Bone

A mathematical model is presented for the transport and precipitation of mineral in refilling osteons. One goal of this model was to explain calcification 'halos,' in which the bone near the haversian canal is more highly mineralized than the more peripheral lamellae, which have been mineralizing longer. It was assumed that the precipitation rate of mineral is proportional to the difference between the local concentration of calcium ions and an equilibrium concentration and that the transport of ions is by either diffusion or some other concentration gradient-dependent process. Transport of ions was assumed to be slowed by the accumulation of mineral in the matrix along the transport path. The model also mimics bone apposition, slowing of apposition during refilling, and mineralization lag time. It was found that simple diffusion cannot account for the transport of calcium ions into mineralizing bone, because the diffusion coefficient is two orders of magnitude too low. If a more rapid concentration gradient-driven means of transport exists, the model demonstrates that osteonal geometry and variable rate of refilling work together to produce calcification halos, as well as the primary and secondary calcification effect reported in the literature.

Martin, Bruce↗
NASA NTRS1994

Vertical Transport Rates in the Stratosphere in 1993 from Observations of CO2, N2O and CH4

Measurements of CO2, N2O and CH4 are analyzed to define hemispheric average vertical exchange rates in the lower stratosphere from November 1992 to October 1993. Effective vertical diffusion coefficients were small in summer, less than or equal to 1 m(exp 2)/sec at altitudes below 25 km; values were similar near the tropopause in winter, but increased markedly with altitude. The analysis suggests possibly longer residence times for exhaust from stratospheric aircraft, and more efficient transport from 20 km to the middle stratosphere, than predicted by many current models. Seasonally-resolved measurements of stratospheric CO2 and N2O provide significant new constraints on rates for global-scale vertical transport.

Wofsy, Steven C.↗
NASA NTRS1994

Mathematical Model for the Mineralization of Bone

A mathematical model is presented for the transport and precipitation of mineral in refilling osteons. One goal of this model was to explain calcification 'halos,' in which the bone near the haversian canal is more highly mineralized than the more peripheral lamellae, which have been mineralizing longer. It was assumed that the precipitation rate of mineral is proportional to the difference between the local concentration of calcium ions and an equilibrium concentration and that the transport of ions is by either diffusion or some other concentration gradient-dependent process. Transport of ions was assumed to be slowed by the accumulation of mineral in the matrix along the transport path. ne model also mimics bone apposition, slowing of apposition during refilling, and mineralization lag time. It was found that simple diffusion cannot account for the transport of calcium ions into mineralizing bone, because the diffusion coefficient is two orders of magnitude too low. If a more rapid concentration gradient-driven means of transport exists, the model demonstrates that osteonal geometry and variable rate of refilling work together to produce calcification halos, as well as the primary and secondary calcification effect reported in the literature.

Martin, Bruce↗