Search NASA⌕ Search

SEARCH · Search NASA

Results for “PARTICLE 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 73 records · Page 4

Blowing in the Wind: II. Creation and Redistribution of Refractory Inclusions in a Turbulent Protoplanetary Nebula

Ca-A1 rich refractory mineral inclusions (CAIs) found at 1-6% mass fraction in primitive chondrites appear to be 1-3 million years older than the dominant (chondrule) components which were accreted into the same parent bodies. A prevalent concern is that it is difficult to retain CAIs for this long against gas-drag-induced radial drift into the sun. We reassess the situation in terms of a hot inner (turbulent) nebula context for CAI formation, using analytical models of nebula evolution and particle diffusion. We show that outward radial diffusion in a weakly turbulent nebula can overwhelm inward drift, and prevent significant numbers of CAI-size particles from being lost into the sun for times on the order of 10(exp 6) years. CAIs can form early, when the inner nebula was hot, and persist in sufficient abundance to be incorporated into primitive planetesimals at a much later time. Small (less than or approximately 0.1 mm diameter) CAIs persist for longer times than large (greater than or approximately 5mm diameter ones. To obtain a quantitative match to the observed volume fractions of CAIs in chondrites, another process must be allowed for: a substantial enhancement of the inner hot nebula in silicate-forming material, which we suggest was caused by rapid inward drift of meter-sized objects. This early in nebula history, the drifting rubble would have a carbon content probably an order of magnitude larger than even the most primitive (CI) carbonaceous chondrites. Abundant carbon in the evaporating material would help keep the nebula oxygen fugacity low, plausibly solar; as inferred for the formation environment of CAIs. The associated production of a larger than canonical amount of CO2 might also play a role in mass-independent fractionation of oxygen isotopes, leaving the gas rich in O-16 as inferred from CAIs and other high temperature condensates.

Cuzzi, Jeffrey N.↗

Trapped particle absorption by the Ring of Jupiter

The interaction of trapped radiation with the ring of Jupiter is investigated. Because it is an identical problem, the rings of Saturn and Uranus are also examined. Data from the Pioneer II encounter, deductions for some of the properties of the rings of Jupiter and Saturn. Over a dozen Jupiter magnetic field models are available in a program that integrates the adiabatic invariants to compute B and L. This program is to label our UCSD Pioneer II encounter data with the most satisfactory of these models. The expected effects of absorbing material on the trapped radiation are studied to obtain the loss rate as a function of ring properties. Analysis of the particle diffusion problem rounds out the theoretical end of the ring absorption problem. Other projects include identification of decay products for energetic particle albedo off the rings and moons of Saturn and a search for flux transfer events at the Jovian magnetopause.

Fillius, W.↗

The Role of Kinetic Instabilities and Waves in Collisionless Magnetic Reconnection

Magnetic reconnection converts magnetic field energy into particle energy by breaking and reconnecting magnetic field lines. Magnetic reconnection is a kinetic process that generates a wide variety of kinetic waves via wave-particle interactions. Kinetic waves have been proposed to play an important role in magnetic reconnection in collisionless plasmas by, for example, contributing to anomalous resistivity and diffusion, particle heating, and transfer of energy between different particle populations. These waves range from below the ion cyclotron frequency to above the electron plasma frequency and from ion kinetic scales down to electron Debye length scales. This review aims to describe the progress made in understanding the relationship between magnetic reconnection and kinetic waves. We focus on the waves in different parts of the reconnection region, namely, the diffusion region, separatrices, outflow regions, and jet fronts. Particular emphasis is placed on the recent observations from the Magnetospheric Multiscale (MMS) spacecraft and numerical simulations, which have substantially increased the understanding of the interplay between kinetic waves and reconnection. Some of the ongoing questions related to waves and reconnection are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Turbulent diffusion from a quasi-kinematical point of view

An estimate for the coefficient of self-diffusion can be derived whenever the Eulerian velocity spectrum is known. The result is more general than those previously obtained. A comparison is made with computed test-particle diffusion in an inviscid two-dimensional Navier-Stokes fluid.

Salu, Y.↗

IUTAM Symposium on Hydrodynamic Diffusion of Suspended Particles

The focus of the symposium was on multiparticle hydrodynamic interactions which lead to fluctuating motion of the particles and resulting particle migration and dispersion or diffusion. Implications of these phenomena were described for sedimentation, fluidization, suspension flows, granular flows, and fiber suspensions. Computer simulation techniques as well as experimental techniques were described. Each session had an invited leadoff talk which overviewed the session topic as well as described the speaker's own related research. Ample time for discussion was included after each talk as well as at the end of each session. The symposium started with a keynote talk on the first evening on What is so puzzling about hydrodynamic diffusion?, which set the tone for the rest of the meeting by emphasizing both recent advances and unanswered issues.

Davis, R. H.↗

Photochemistry Saturn's Atmosphere: Effects of an Influx of External Oxygen - 2

We use a one-dimensional diurnally averaged model of photochemistry and diffusion in Saturn's stratosphere to investigate the influence of extraplanetary debris on atmospheric chemistry. In particular, we consider the effects of an influx of oxygen from micrometeoroid ablation or from ring-particle diffusion; the contribution from cometary impacts, satellite debris, or ring vapor is deemed to be less important. The photochemical model results are compared directly with Infrared Space Observatory (ISO) observations to constrain the influx of extraplanetary oxygen to Saturn. From the ISO observations, we determine that the column densities of CO2 and H2O above 10 mbar in Saturn's atmosphere are (6.3 +/- 1) x 10(exp 14) and (1.4 +/- 0.4) x 10(exp 15)/ square cm, respectively; our models indicate that a globally averaged oxygen influx of (4+/-2) x 10(exp 6) O atoms /sq cm/s is required to explain these observations. Models with a locally enhanced influx of H20 operating over a small fraction of the projected area do not provide as good a fit to the ISO H2O observations. If volatile oxygen compounds comprise one-third to one-half of the exogenic source by mass, then Saturn is currently being bombarded with (3 +/- 2) x 10(exp -16) g/square cm/s of extraplanetary material. To reproduce the observed CO2/H2O ratio in Saturn's stratosphere, some of the exogenic oxygen must arrive in the form of a carbon-oxygen bonded species such as CO or CO2. An influx consistent with the composition of cometary ices fails to reproduce the high observed CO2/H2O ratio, suggesting that (i) the material has ices that are slightly more carbon-rich than is typical for comets, (ii) a contribution from an organic-rich component is required, or (iii) some of the hydrogen-oxygen bonded material is converted to carbon-oxygen bonded material without photochemistry (e.g., during the ablation process). We have also reanalyzed the 5-micron CO observations of Noll and Larson and determine that the CO lines are most sensitive to the 100- to 400-mbar column density for which we derive a range of (0.7-1.5) x 10(exp 17)/square cm; the CO observations do not allow us to distinguish between an external or internal source of CO on Saturn. If we assume that all the extraplanetary oxygen derives from a micrometeoroid source, then the unfocused dust flux at 9.5 AU must be (i) (1 +/- 0.7) x 10(exp -16) g/square cm/s if interstellar grains are the source of the external oxygen on Saturn, (ii) (4 +/- 3) x 10 (exp -17) g/ sq. cm./s if IDPs on randomly inclined, highly eccentric orbits (e.g., Halley-type comet grains) are the source of the external oxygen, or (iii) (2+/- 1.4) x 10(exp -18)/g/sq cm/s if IDPs on low inclination, low eccentricity orbits (e.g., Kuiper-belt grains) are the source of the external oxygen. These estimates can be used in combination with future Cassini dust detection data to determine the ultimate source of the dust at Saturn's distance from the Sun.

Moses, Julianne I.↗

A new continuum model for suspensions of gyrotactic micro-organisms

A new continuum model is formulated for dilute suspensions of swimming micro-organisms with asymmetric mass distributions. Account is taken of randomness in a cell's swimming direction, p, by postulating that the probability density function for p satisfies a Fokker-Planck equation analogous to that obtained for colloid suspensions in the presence of rotational Brownian motion. The deterministic torques on a cell, viscous and gravitational, are balanced by diffusion, represented by an isotropic rotary diffusivity Dr, which is unknown a priori, but presumably reflects stochastic influences on the cell's internal workings. When the Fokker-Planck equation is solved, macroscopic quantities such as the average cell velocity Vc, the particle diffusivity tensor D and the effective stress tensor sigma can be computed; Vc and D are required in the cell conservation equation, and sigma in the momentum equation. The Fokker-Planck equation contains two dimensionless parameters, lambda and epsilon; lambda is the ratio of the rotary diffusion time Dr-1 to the torque relaxation time B (balancing gravitational and viscous torques), while epsilon is a scale for the local vorticity or strain rate made dimensionless with B. In this paper we solve the Fokker-Planck equation exactly for epsilon = 0 (lambda arbitrary) and also obtain the first-order solution for small epsilon. Using experimental data on Vc and D obtained with the swimming alga, Chlamydomonas nivalis, in the absence of bulk flow, the epsilon = 0 results can be used to estimate the value of lambda for that species (lambda approximately 2.2; Dr approximately 0.13 s-1). The continuum model for small epsilon is then used to reanalyse the instability of a uniform suspension, previously investigated by Pedley, Hill & Kessler (1988). The only qualitatively different result is that there no longer seem to be circumstances in which disturbances with a non-zero vertical wavenumber are more unstable than purely horizontal disturbances. On the way, it is demonstrated that the only significant contribution to sigma, other than the basic Newtonian stress, is that derived from the stresslets associated with the cells' intrinsic swimming motions.

NASA Discipline Number 28-20↗

Reversible electron heating vs. wave-particle interactions in quasi-perpendicular shocks

The energy necessary to explain the electron heating in quasi-perpendicular collisionless shocks can be derived either from the electron acceleration in the d.c. cross shock electric potential, or by the interactions between the electrons and the waves existing in the shock. A Monte Carlo simulation has been performed to study the electron distribution function evolution through the shock structure, with and without particle diffusion on waves. This simulation has allowed us to clarify the relative importance of the two possible energy sources; in particular it has been shown that the electron parallel temperature is determined by the d.c. electromagnetic field and not by any wave-particle-induced heating. Wave particle interactions are effective in smoothing out the large gradients in phase space produced by the 'reversible' motion of the electrons, thus producing a 'cooling' of the electrons.

Veltri, P.↗

Simulation of the velocity diffusion of charged particles in turbulent magnetic fields

Numerical simulation experiments have been performed in order to study the velocity diffusion of charged particles in a static turbulent magnetic field. By following orbits of particles moving in a large ensemble of random magnetic field realizations with suitably-chosen statistical properties, a pitch-angle diffusion coefficient is derived. Results are presented for a variety of particle rigidities and r.m.s. random field strengths and a power spectrum typical of conditions found in interplanetary space.

Kaiser, T. B.↗

Trapped radiation belts of Saturn - First look

Data on the magnetosphere of Saturn obtained with the trapped radiation detector package on board the Pioneer 11 spacecraft is reported. Radiation belt profiles determined by the trapped radiation detectors on Pioneer 10 and 11 indicate that Saturn's magnetosphere is intermediate in size between those of the earth and Jupiter, with particle intensities similar to those of the earth. The outer region of the Saturn magnetosphere is found to contain particles of lower energy than the outer region, being strongly influenced by the time-varying solar wind. The moons and rings of Saturn are observed to be effective absorbers of trapped particles, confirming the discoveries of the F ring, the Pioneer ring division and the moon 1979 S 2. Particle diffusion rates are used to estimate a cross-sectional area of greater than 7 x 10 to the 13th sq cm and an opacity greater than 0.00001 for the F ring. It is suggested that cosmic-ray albedo neutron decay be studied as a possible source of energetic particles in the inner magnetosphere of Saturn.

Fillius, W.↗

Ion acceleration in impulsive solar flares

Nonrelativistic spectra of protons and ions accelerated in impulsive solar flares are derived using more realistic turbulence power spectra. The calculation is based on a particle transport equation extracted from a second step acceleration model containing stochastic acceleration. The turbulence model is generalized to waves with a small angle to the magnetic field vector and to turbulence power spectra with spectral indices s smaller than 2. Due to the occurrence of impulsive flares at low coronal heights, Coulomb losses at the dense coronal plasma and diffusive particle escape are taken into account. The ion spectra show deviations from long-duration spectra near the Coulomb barrier, where the losses become maximal. The Z-squared/A-dependence of the Coulomb losses leads to spectral variations for different ions. We present a method to estimate the turbulence parameters and injection conditions of the flare particles using ion ratios like Fe/O of impulsive flares.

Steinacker, Jurgen↗

Characteristics of the magnetospheric source of interplanetary energetic particles

The Earth's bow shock is frequently cited as an example of an astrophysical shock where particle acceleration is observed. However, because energetic particles observed upstream of the bow shock may be accelerated within the magnetosphere, it is important to understand the properties of the magnetospheric source. A first order picture of the spatial distribution of magnetospheric particles in the magnetosheath and upstream is obtained by mapping those magnetic field lines which drape over the magnetopause through the bow shock. Subsets of these field lines that connect to potential sites of magnetic merging on the magnetopause are also traced in the event that leakage occurs preferentially where normal components of the field are present across the boundary. The results can be used to determine whether the so-called diffuse particles observed upstream are accelerated locally or within the magnetosphere.

Luhmann, J. G.↗

Compression Acceleration of Protons and Heavier Ions at the Heliospheric Current Sheet

Abstract Recent observations by the Parker Solar Probe (PSP) suggest that protons and heavier ions are accelerated to high energies by magnetic reconnection at the heliospheric current sheet (HCS). By solving the energetic particle transport equation in large-scale MHD simulations, we study the compression acceleration of protons and heavier ions in the reconnecting HCS. We find that the acceleration of multispecies ions results in nonthermal power-law distributions with a spectral index consistent with the PSP observations. Our study shows that the high-energy cutoff of protons can reach E max ∼ 0.1 –1 MeV depending on the particle diffusion coefficients. We also study how the high-energy cutoff of different ion species scales with the charge-to-mass ratio E max ∝ ( Q / M ) α . When determining the diffusion coefficients from the quasi-linear theory with a Kolmogorov magnetic power spectrum, we find that α ∼ 0.4, which is somewhat smaller than α ∼ 0.7 observed by PSP.

79 ASTRONOMY AND ASTROPHYSICS↗

Transverse particle acceleration and diffusion in a planetary magnetic field

A general model of particle acceleration by plasma waves coupled with adiabatic radial diffusion in a planetary magnetic field is developed. The model assumes that a spectrum of lower hybird waves is present to resonantly accelerate ions transverse to the magnetic field. The steady state Green's function for the combined radial diffusion and wave acceleration equation is found in terms of a series expansion. The results provide a rigorous demonstration of how a quasi-Maxwellian distribution function is formed in the absence of particle collisons and elucidate the nature of turbulent heating of magnetospheric plasmas. The solution is applied to the magnetosphere of Neptune for which a number of examples are given illustrating how the spectrum of pickup N(+) ions from Triton evolves.

Barbosa, D. D.↗

Atmospheric concentrations of submicron contact-freezing nuclei

Atmospheric concentrations of contact-freezing nuclei were measured using a technique primarily sensitive to submicron aerosol particles. Diffusion and phoretic forces were relied on for the capture of nuclei by supercooled drops of distilled water exposed to the sample air. Nucleus concentrations were deduced from the rate at which the drops were observed to freeze, interpreting that rate on the basis of a theoretical prediction of aerosol capture rate for different assumed sizes of the nuclei. Measurements at Laramie, Wyoming, yielded average concentrations of contact-freezing nuclei of 1.7/L at -15 C and 3.1/L at -18 C for an assumed radius of 0.01 micron for the nucleating particles.

Deshler, Terry↗

Particle control via cryopumping and its impact on the edge plasma profiles of Alcator C-Mod

At the high n e proposed for high-field fusion reactors, it is uncertain whether ionization, as opposed to plasma transport, will be most influential in determining n e at the pedestal and separatrix. A database of Alcator C-Mod discharges is analyzed to evaluate the impact of source modification via cryopumping. The database contains similarly-shaped H-modes at fixed I P = 0.8 MA and B t = 5.4 T, spanning a large range in P net and ionization. Measurements from an edge Thomson scattering system are combined with those from a midplane-viewing Ly α camera to evaluate changes to n e and T e in response to changes in ionization rates, S ion ∙ $n^{sep}_e$ and $T^{ped}_e$ are found to be most sensitive to changes to $S^{sep}_{ion}$, as opposed to $n^{ped}_e$ and $T^{sep}_e$. Dimensionless quantities, namely α MHD and v*, are found to regulate attainable pedestal values. Select discharges at different values of P net and in different pumping configurations are analyzed further using SOLPS-ITER. It is determined that changes to plasma transport coefficients are required to self-consistently model both plasma and neutral edge dynamics. Pumping is found to modify the poloidal distribution of atomic neutral density, n 0 , along the separatrix, increasing n 0 at the active X-point. Opaqueness to neutrals from high n e in the divertor is found to play a role in mediating neutral penetration lengths and hence, the poloidal distribution of neutrals along the separatrix. Pumped discharges thus require a larger particle diffusion coefficient than that inferred purely from 1D experimental profiles at the outer midplane.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Self-consistent equilibrium and transport simulations for NSTX-U plasmas enhanced via machine learning surrogate models

The Control-Oriented Transport SIMulator (COTSIM) is an advanced equilibrium and transport code designed for simulating tokamak discharges at computational speeds suitable for control applications. COTSIM’s modular framework enables users to select models that balance accuracy with speed according to specific needs, allowing the code to operate from fast to faster-than-real-time performance levels. This work presents recent enhancements to COTSIM’s predictive accuracy for NSTX-U scenarios, achieved by integrating neural-network-based surrogate models and self-consistent equilibrium calculations. To improve source deposition predictions, a surrogate model for NUBEAM has been incorporated. Additionally, a surrogate model for the Multi-Mode Module (MMM) now supports predictions of anomalous thermal, momentum, and particle diffusivities—key factors for modeling the evolution of temperature and rotation. Each surrogate model was specifically trained for the NSTX-U operational regime to enhance COTSIM’s accuracy while maintaining computational efficiency. Moreover, COTSIM now couples fixed-boundary equilibrium solvers with its transport solvers, enabling self-consistent predictions of plasma profiles and equilibrium evolution over the discharge. Simulation results demonstrate strong agreement between COTSIM and TRANSP predictions for NSTX-U discharges. These substantial advancements expand COTSIM’s utility in model-based control applications for NSTX-U. Potential applications include simultaneous optimization of equilibrium and transport scenarios, integration into digital twins, real-time profile estimation (e.g., temperature and rotation) from limited or noisy measurements, and advanced feedback-based scenario control.

Equilibrium and transport modeling↗