Search NASA⌕ Search

SEARCH · Search NASA

Results for “jump 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 19 records

Estimation of the viscosities of silicate liquids at high pressure from measurements of oxygen diffusivity

The dynamics and evolution of a magma ocean depend significantly on the structure and physical properties of silicate liquids at high pressure. Viscosity is particularly important because of its effect on the dynamics of convection and crystal setting. Traditionally, the viscosities of silicate liquids have been measured at high pressure by falling sphere viscometry but, due to some limitations of this technique, few measurements have been made at pressures exceeding 2.5 GPa. An alternative approach, which has previously been proven up to 2 GPa, is to use high pressure measurements of oxygen self-diffusivity (D) to estimate viscosity (eta) from the Eyring relationship eta = kT/D lambda (where k is the Boltzmann constant, T is absolute temperature, and lambda is the diffusive jump distance). We have performed preliminary experiments on Na2Si4O9 liquid in order to explore the feasibility of this method up to pressures in excess of 10 GPa. Diffusion couples were prepared from glass starting materials with one half of the sample enriched in O-18 and the other half containing the natural abundance (0.2 wt percent). Diffusion experiments were performed at 1600 to 1800 C at pressures in the range of 2.5 to 10 GPa for times up to 6 minutes using a 1200 ton multianvil apparatus. Oxygen diffusivities were derived from the resulting O-18 concentration profiles, which were analyzed using an ion microprobe. At 1800 C, oxygen diffusivities increase continuously from 1 x 10(exp -10) m(exp -2) s(exp -1) at 2.5 Gpa to 5 x 10(exp -10) m(exp -2) s(exp -1) at 10 GPa. From the Eyring relationship, the viscosity of Na2SiO9 liquid is predicted to decrease by approx. 0.7 log units over this pressure range at 1800 C. These trends, which can be related to changes in the structure of the liquid at high pressure, are in agreement with results of molecular dynamics calculations. These results suggest that it should be possible to estimate the viscosities of silicate liquids up to at least 15 GPa and 2200 C by this method.

Rubie, D. C.↗

[Results of Simulation Studies]

Lattice Monte Carlo and off-lattice molecular dynamics simulations of h(sub 1)t(sub 4) and h(sub 4)t(sub l) (head/tail) amphiphile solutions have been performed as a function of surfactant concentration and temperature. The lattice and off-lattice systems exhibit quite different self-assembly behavior at equivalent thermodynamic conditions. We found that in the weakly aggregating regime (no preferred-size micelles), all models yield similar micelle size distributions at the same average aggregation number, albeit at different thermodynamic conditions (temperatures). In the strongly aggregating regime, this mapping between models (through temperature adjustment) fails, and the models exhibit qualitatively different micellization behavior. Incipient micellization in a model self-associating telechelic polymer solution results in a network with a transient elastic response that decays by a two-step relaxation: the first is due to a heterogeneous jump-diffusion process involving entrapment of end-groups within well-defined clusters and this is followed by rapid diffusion to neighboring clusters and a decay (terminal relaxation) due to cluster disintegration. The viscoelastic response of the solution manifests characteristics of a glass transition and entangled polymer network.

Source record↗

Simulating the Diffusion of Hydrogen in Amorphous Silicates: A ‘Jumping’ Migration Process and its Implications for Solar Wind Implanted Lunar Volatiles

We use molecular dynamics (MD) simulations to better explain the movement of atomic hydrogen in amorphous silica and quantify the planetary science implications of these findings. Previous MD simulations had a large range of predicted values and did not agree well with experiment. Our simulations sample atomic motion for a longer duration and consider a wider range of temperatures than previous simulations. In contrast to constant atomic motion, the hydrogen atoms were shown to undergo random intermittent jumps from one oxygen atom to another, the number of which increase with temperature. Predicted diffusion coefficients had a better agreement to experimental values than previous MD simulations, suggesting the importance of longer simulation durations for better statistics. The low activation energy and jumps observed at lunar temperatures do not support the theory of diurnal variations in OH content for an undamaged amorphous silica surface. Instead, we conclude that energetic solar wind impacts can induce two competing atomic hydrogen motion processes in the exposed surface: A prompt effect that induces jumps in the temperature spike volume, but also a long term effect of damage in the structure that traps atomic hydrogen. We then use SDTrimSP to quantify the damage created during exposure and MD to demonstrate the H retention and trapping near these defects. Damage was shown to be dependent on impact energy, with defects easily retaining implanted hydrogen. MD results like those presented herein on unweathered surfaces are therefore most relevant to magnetic anomalies. As a result, we demonstrate the importance of lunar volatile models to account for the damage state of the substrate when modelling hydrogen diffusion, retention, and subsequent OH/water production.

Liam S. Morrissey↗

Tensile Creep of Polycrystalline Near-Stoichiometric NiAl

Long term tensile creep studies were conducted on binary NiAl in the temperature range 700-1200 K with the objectives of characterizing and understanding the creep mechanisms. Inverse and normal primary creep curves were observed depending on stress and temperature. It was concluded that the creep of NiAl is limited by dislocation mobility. The stress exponent for creep, n, increased from 5.5 at 1200 K to 13.9 at 700 K. The true activation energy for creep, Qc, was constant and equal to about 400 kJ per mole between 20 and 50 MPa but decreased to a constant value of 250 kJ per mole between 50 and 110 MPa. The activation energy was observed to be stress dependent above 110 MPa. The tensile creep results reported in this investigation were compared with compression creep data reported in the literature. A detailed discussion of the probable dislocation creep mechanisms governing compressive and tensile creep of NiAl is presented. It is concluded that the non-conservative motion of jogs on screw dislocations influenced the nature of the primary creep curves, where the climb of these jogs involves either the next nearest neighbor or the six-jump cycle vacancy diffusion mechanism. The probable nature of the atom vacancy exchange that occur within the core of an edge dislocation undergoing climb in NiAl are schematically examined.

Raj, Sai V.↗

Compression and Tensile Creep of Binary NiAl

Compression creep and long term tensile creep studies were conducted on cast and extruded binary NiAl in the temperature range 700-1200 K with the objectives of characterizing and understanding the creep mechanisms. Inverse and normal primary creep curves were observed in both compression and tension creep depending on stress and temperature although an asymmetrical response was observed under these two stress states. It was concluded that the primary creep of NiAl is limited by dislocation mobility. The stress exponents, n, for compression and tensile creep were similar varying between about 5 and 14. However, there were significant differences in the stress dependence of the activation energies for compression and tensile creep. The true activation energy for tensile creep, Q(sub c), was constant and equal to about 400 kJ/mol between 20 and 50 MPa but decreased to a constant value of 250 kJ/mol between 50 and 110 MPa. The activation energy was observed to be inversely stress dependent above 110 MPa. In contrast, Q(sub c) = 300 kJ/mol for compression creep was constant between 25 and 70 MPa and inversely dependent on the true stress above 70 MPa. A detailed discussion of the probable dislocation creep mechanisms governing compressive and tensile creep of NiAl is presented. It is concluded that the non-conservative motion of jogs on screw dislocations influenced the nature of the primary creep curves, where the climb of these jogs involves either the next nearest neighbor or the six-jump cycle vacancy diffusion mechanism. The probable natures of the atom-vacancy exchange that occur within the core of an edge dislocation undergoing climb in NiAl are schematically examined.

Raj, Sai V.↗

The Hot Horizontal-Branch Stars in NGC288 - Effects of Diffusion and Stratification on Their Atmospheric Parameters*

Context. NGC288 is a globular cluster with a well developed blue horizontal branch covering the so-called u-jump which indicates the onset of diffusion. It is therefore well suited to study the effects of diffusion in blue horizontal branch (HB) stars. Aims. We compare observed abundances to predictions from stellar evolution models calculated with diffusion and from stratified atmospheric models. We verify the effect of using stratified model spectra to derive atmospheric parameters. In addition we investigate the nature of the overluminous blue HB stars around the u-jump. Methods. We define a new photometric index sz from uvby measurements that is gravity sensitive between 8 000K and 12 000 K. Using medium-resolution spectra and Stroemgren photometry we determine atmospheric parameters (Teff, logg) and abundances for the blue HB stars. We use both homogeneous and stratified model spectra for our spectroscopic analyses. Results. The atmospheric parameters and masses of the hot HB stars in NGC288 show a behaviour seen also in other clusters for temperatures between 9 000K and 14 000 K. Outside this temperature range, however, they follow rather the results found for such stars in (omega)Cen. The abundances derived from our observations are for most elements (except He and P) within the abundance range expected from evolutionary models that include the effects of atomic diffusion and assume a surface mixed mass of 10(exp −7) M. The abundances predicted by stratified model atmospheres are generally significantly more extreme than observed, except for Mg. The use of stratified model spectra to determine effective temperatures, surface gravities and masses moves the hotter stars to a closer agreement with canonical evolutionary predictions. Conclusions. Our results show definite promise towards solving the long-standing issue of surface gravity and mass discrepancies for hot HB stars, but there is still much work needed to arrive at a self-consistent solution.

horizontal branch (HB)↗

Radiative shock structure.

Shock wave structure in ionized hydrogen, using Eddington approximation to solve Rankine- Hugoniot jump equation

DIFFUSION WAVE↗

Effect of Static Strains on Diffusion

A theory is developed that gives the diffusion coefficient in strained systems as an exponential function of the strain. This theory starts with the statistical theory of the atomic jump frequency as developed by Vineyard. The parameter determining the effect of strain on diffusion is related to the changes in the inter-atomic forces with strain. Comparison of the theory with published experimental results for the effect of pressure on diffusion shows that the experiments agree with the form of the theoretical equation in all cases within experimental error.

Girifalco, L. A.↗

Particle drift, diffusion, and acceleration at shocks

The gradient and curvature drifts implicit in change of the ambient magnetic field at a hydromagnetic shock wave are incorporated into the diffusive theory of shock acceleration of charged particles. The conventional jump condition at the shock is modified by a term incorporating the large drift along the shock plane. This term vanished identically for one-dimensional systems, but must be included in general for shocks which are finite in transverse extent or which have transverse structure. It is found that the effect of the drift is such that the transverse drift rate is proportional to the acceleration rate, and for perpendicular shocks is exactly equal to the rate of change of energy in the V x B electric field observed in the shock frame. This establishes a connection with the 'shock drift' models which neglect diffusion.

Jokipii, J. R.↗

A thick-thin approximation for radiative transfer.

First- and second-order formulations of a thick-thin approximate technique for application to the solution of energy transport when nongray radiation is important are presented. The virtue of the technique is that it eliminates the cumbersome integral terms which usually appear. It is shown that the temperature distributions calculated with the technique agree well with the results of more rigorous calculations if a proper value of the critical optical depth, which separates the thick and thin regions of the spectrum, is used. The proper value of the critical optical depth is less than one and may be less than 0.5. It is also shown that the jump boundary conditions which are frequently employed with the Rosseland, or diffusion, approximation should incorporate the Planck mean absorption coefficient for a nongray gas. This is in contradiction to suggestions appearing in the literature.

Olstad, W. B.↗

The relation between ferroelasticity and superconductivity

The high-temperature superconductivity is explained widely by the layered crystal structure. The one- and two-dimensional subsystems and their interaction are investigated here. It is assumed that the high-T(sub c) superconductivity takes place in the two-dimensional subsystem and the increase of the phase transition temperature from 60 K up to 90 K is the consequence of turning on the influence of one-dimensional chains. The interaction between the two subsystems is transferred along the c axis by the phonons of breathing mode, which causes the hybridization of the electronic bonds between these subsystems. The experimental works indicate that the existence of both the chains Cu(1)-O and their interaction with the superconducting plane of Cu(2)-O modify the temperature of the transition to the superconducting state. It is seen from the neutron scattering data that the rates of the interatomic distance dependencies on temperature are changed around 140 K and 90 K. The 'zig-zag' order in Cu(1)-O chains has been postulated but, on the other hand, the vibrations with a large amplitude only were reported. The bi-stabilized situation of the oxygen ions can be caused by the change of distance between these ions and the Ba ions. It leads to the appearance of a two-well potential. Its parameters depend on temperature and the dynamics of the oxygen ions' movement. They can induce the antipolar order, which can be, however, more or less chaotic. The investigation of the ferroelastic properties of Y-Ba-Cu-O samples lead to the conclusion that they are related to jumps of ions inside the given chain and not to a diffusion between different sites in the ab plane. Researchers deduce thus that the fluctuating oxygen ions from these chains create dipoles in the ab plane. They can be described with the pseudo-spin formalism/ - Pauli matrices/. The system can be described with the Ising model. The pseudo-spins interact with phonons and influence the superconductivity in the second subsystem.

Molak, A.↗

The relation between ferroelasticity and superconductivity

The high-temperature superconductivity is explained widely by the layered crystal structure. The one- and two-dimensional subsystems and their interaction are investigated here. It is assumed that the high-T(sub c) superconductivity takes place in the two-dimensional subsystem and the increase of the phase transition temperature from 60 K up to 90 K is the consequence of turning on the influence of one-dimensional chains. The interaction between the two subsystems is transferred along the c axis by the phonons of breathing mode, which causes the hybridization of the electronic bonds between these subsystems. The experimental works indicate that the existence of both the chains Cu(1)-O and their interaction with the superconducting plane of Cu(2)-O modify the temperature of the transition to the superconducting state. It is seen from the neutron scattering data that the rates of the interatomic distance dependencies on temperature are changed around 240 K and 90 K. The 'zig-zag' order in Cu(1)-O chains has been postulated but, on the other hand, the vibrations with a large amplitude only were reported. The bi-stabilized situation of the oxygen ions can be caused by the change of distance between these ions and the Ba ions. It leads to the appearance of a two-well potential. Its parameters depend on temperature and the dynamics of the oxygen ions' movement. They can induce the antipolar order, which can be, however, more or less chaotic. The investigation of the ferroelastic properties of Y-Ba-Cu-O samples lead to the conclusion that they are related to jumps of ions inside the given chain and not to a diffusion between different sites in the ab plane. Researchers deduce, thus, that the fluctuating oxygen ions from these chains create dipoles in the ab plane. They can be described with the pseudo-spin formalism (- Pauli matrices). The system can be described with the Ising model. The pseudo-spins interact with phonons and influence the superconductivity in the second subsystem.

Molak, A.↗

Steady state and dynamical structure of a cosmic-ray-modified termination shock

A hydrodynamic model is developed for the structure of a cosmic-ray-modified termination shock. The model is based on the two-fluid equations of diffuse shock acceleration (Drury and Volk, 1981). Both the steady state structure of the shock and its interaction with outer heliospheric disturbances are considered. Under the assumption that the solar wind is decelerated by diffusing interstellar cosmic rates, it is shown that the natural state of the termination shock is a gradual deceleration and compression, followed by a discontinuous jump to a downstream state which is dominated by the pressure contribution of the cosmic rays. A representative model is calculated for the steady state which incorporates both interstellar cosmic ray mediation and diffusively accelerated anomalous ions through a proposed thermal leakage mechanism. The interaction of large-scale disturbances with the equilibrium termination shock model is shown to result in some unusual downstream structure, including transmitted shocks and cosmic-ray-modified contact discontinuities. The structure observed may be connected to the 2-kHz outer heliospheric radio emission (Cairns et al., 1992a, b). The time-dependent simulations also demonstrate that interaction with solar wind compressible turbulence (e.g., traveling interplanetary shocks, etc.) could induce the termination shock to continually fluctuate between cosmic-ray-dominated and gas-dynamic states. This fluctuation may represent a partial explanation of the galactic cosmic ray modulation effect and illustrates that the Pioneer and Voyager satellites will encounter an evolving shock whose structure and dynamic properties are strongly influence by the mediation of interstellar and anomalous cosmic rays.

Donohue, D. J.↗

Computer simulations of cosmic-ray diffusion near supernova remnant shock waves

A plasma simulation model was used to study the resonant interactions between streaming cosmic-ray ions and a self-consistent spectrum of Alfven waves, such as might exist in the interstellar medium upstream of a supernova remnant shock wave. The computational model is a hybrid one, in which the background interstellar medium is an MHD fluid and the cosmic-rays are discrete kinetic particles. The particle sources for the electromagnetic fields are obtained by averaging over the fast cyclotron motions. When the perturbed magnetic field is larger than 10 percent of the background field, the macro- and microphysics are no longer correctly predicted by quasi-linear theory. The particles are trapped by the waves and show sharp jumps in their pitch-angles relative to the background magnetic field, and the effective ninety-degree scattering time for diffusion parallel to the background magnetic field is reduced to between 5 and 30 cyclotron periods. Simulation results suggest that Type 1 supernova remnants may be the principal sites of cosmic ray acceleration.

Max, C. E.↗

Wall boundary equations with slip and catalysis for multicomponent, nonequilibrium gas flows

Boundary equations obtained for a low Reynolds number, high enthalpy gas flow in regions of velocity slip and temperature jump are presented. The formulation treats, through a first-order distribution function used to include multicomponent diffusion, a multicomponent gas mixture that may be in nonequilibrium with finite-rate catalytic recombination occurring on the wall. In the boundary equations, which are obtained for use in flow-field calculations applicable to low-density flow regimes, a simplified gas/wall interaction is assumed wherein individual atoms or molecules either reflect specularly off the wall or stick and are fully accommodated. Fluxes in terms of evaluated integrals over the distribution function and integrals necessary for determining fluxes are given.

Scott, C. D.↗

Edge Diffusion Flame Propagation and Stabilization Studied

In most practical combustion systems or fires, fuel and air are initially unmixed, thus forming diffusion flames. As a result of flame-surface interactions, the diffusion flame often forms an edge, which may attach to burner walls, spread over condensed fuel surfaces, jump to another location through the fuel-air mixture formed, or extinguish by destabilization (blowoff). Flame holding in combustors is necessary to achieve design performance and safe operation of the system. Fires aboard spacecraft behave differently from those on Earth because of the absence of buoyancy in microgravity. This ongoing in-house flame-stability research at the NASA Glenn Research Center is important in spacecraft fire safety and Earth-bound combustion systems.

Takahashi, Fumiaki↗