Search NASA⌕ Search

SEARCH · Search NASA

Results for “radial 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 991 records · Page 55

Machine Learning Potentials with the Iterative Boltzmann Inversion: Training to Experiment

Methodologies for training machine learning potentials (MLPs) with quantum-mechanical simulation data have recently seen tremendous progress. Experimental data have a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. Here, we investigate a training procedure based on iterative Boltzmann inversion that produces a pair potential correction to an existing MLP using equilibrium radial distribution function data. By applying these corrections to an MLP for pure aluminum based on density functional theory, we observe that the resulting model largely addresses previous overstructuring in the melt phase. Interestingly, the corrected MLP also exhibits improved performance in predicting experimental diffusion constants, which are not included in the training procedure. The presented method does not require autodifferentiating through a molecular dynamics solver and does not make assumptions about the MLP architecture. Our results suggest a practical framework for incorporating experimental data into machine learning models to improve the accuracy of molecular dynamics simulations.

36 MATERIALS SCIENCE↗

Modelling enhanced confinement in drift-wave turbulence

Here, the results of modeling studies of an enhanced confinement in the drift wave turbulent plasma of the CSDX linear device are presented. The mechanism of enhanced confinement is investigated here using a reduced 1D, time-dependent model, which illustrates the exchange of enstrophy between two disparate scale structures: the mesoscale flow and profile, and the turbulence intensity fields. Mean density, mean vorticity, and turbulent potential enstrophy are the variables for this model. Total potential enstrophy is conserved in this model. Vorticity mixing occurs on a scale length related to an effective Rhines' scale of turbulence, and shrinks as both density and vorticity gradients steepen. Numerical results obtained from solution of the model agree well with the experimental data from CSDX showing: (i) a steepening of the mean density profile, indicating a radial transport barrier formation, (ii) the development of a radially sheared azimuthal flow velocity that coincides with the density steepening and initiates a turbulence quench, and (iii) negative Reynolds work values, indicating that fluctuations drive the shear flow. These observations as the magnitude of the magnetic field B increases are recovered using purely diffusive expressions for the vorticity and density fluxes. A new dimensionless turbulence parameter R DT -defined as the ratio of the integrated potential enstrophy transfer from turbulence to the flow, to the integrated potential enstrophy production due to relaxation of the density gradient is introduced as a turbulence collapse indicator that detects when the enhanced confinement state is triggered

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

VERY LARGE TELESCOPE/X-SHOOTER SPECTROSCOPY OF THE CANDIDATE BLACK HOLE X-RAY BINARY MAXI J1659-152 IN OUTBURST

We present the optical to near-infrared spectrum of MAXI J1659-152 during the onset of its 2010 X-ray outburst. The spectrum was obtained with X-shooter on the ESO Very Large Telescope early in the outburst simultaneous with high-quality observations at both shorter and longer wavelengths. At the time of the observations, the source was in the low-hard state. The X-shooter spectrum includes many broad ({approx}2000 km s{sup -1}), double-peaked emission profiles of H, He I, and He II, characteristic signatures of a low-mass X-ray binary during outburst. We detect no spectral signatures of the low-mass companion star. The strength of the diffuse interstellar bands results in a lower limit to the total interstellar extinction of A{sub V} {approx_equal} 0.4 mag. Using the neutral hydrogen column density obtained from the X-ray spectrum we estimate A{sub V} {approx_equal} 1 mag. The radial velocity structure of the interstellar Na I D and Ca II H and K lines results in a lower limit to the distance of {approx}4 {+-} 1 kpc, consistent with previous estimates. With this distance and A{sub V} , the dereddened spectral energy distribution represents a flat disk spectrum. The two 10 minute X-shooter spectra show significant variability in the red wing of the emission-line profiles, indicating a global change in the density structure of the disk, though on a timescale much shorter than the typical viscous timescale of the disk.

79 ASTRONOMY AND ASTROPHYSICS↗

Structure and dynamics of the ionosphere

The structure of the Venus ionosphere and the major processes occurring within it are summarized. The daytime ionosphere is created by solar EUV radiation incident on the thermosphere; it is in photochemical equilibrium near its peak at about 142 km, where O2(+) is the major ion, and near diffusive equilibrium in its upper regions, where the major ion is O(+). The day-to-night plasma pressure gradient across the terminator drives a nightward ion flow which, together with electron precipitation, contributes to the formation of the nighttime ionosphere. Large-scale radial holes or plasma depletions extending downwards to nearly the ionization peak in the antisolar region are also observed which are associated with regions of strong radial magnetic fields. The ionopause is a highly dynamic and complex surface, extending from an average altitude of 290 km at the subsolar point to about 1000 km at the terminator and from 200 to over 3000 km on the nightside. A variety of solar wind interaction products are observed in the mantle, a transition region between the ionospheric plasma and the flowing shocked solar wind.

Nagy, A. F.↗

Onion-shell model for cosmic ray electrons and radio synchrotron emission in supernova remnants

The spectrum of cosmic ray electrons, accelerated in the shock front of a supernova remnant (SNR), is calculated in the test-particle approximation using an onion-shell model. Particle diffusion within the evolving remnant is explicity taken into account. The particle spectrum becomes steeper with increasing radius as well as SNR age. Simple models of the magnetic field distribution allow a prediction of the intensity and spectrum of radio synchrotron emission and their radial variation. The agreement with existing observations is satisfactory in several SNR's but fails in other cases. Radiative cooling may be an important effect, especially in SNR's exploding in a dense interstellar medium.

Beck, R.↗

Extensional flow convecting a reactant undergoing a first order homogeneous reaction and diffusional mass transfer from a sphere at low to intermediate Peclet and Damkohler numbers

Forced convective diffusion-reaction is considered for viscous axisymmetric extensional convecting velocity in the neighborhood of a sphere. For Peclet numbers in the range 0.1 less than or equal to Pe less than or equal to 500 and for Damkohler numbers increasing with increasing Pe but in the overall range 0.02 less than or equal to Da less than or equal to 10, average and local Sherwood numbers have been computed. By introducing the eigenfunction expansion c(r, Theta) = Sum of c(n)(r)P(n)(cos Theta) into the forced convective diffusion equation for the concentration of a chemical species undergoing a first order homogeneous reaction and by using properties of the Legendre functions Pn(cos Theta), the variable coefficient PDE can be reduced to a system of N + 1 second order ODEs for the radial functions c(sub n)(r), n = 0, 1, 2,..., N. The adaptive grid algorithm of Pereyra and Lentini can be used to solve the corresponding 2(N + 1) first order differential equations as a two-point boundary value problem on 1 less than or equal to r less than or equal to r(sub infinity). Convergence of the expansion for a specific value of N can thus be established and provides 'spectral' behavior as well as the full concentration field c(r, Theta).

Shah, N. Y.↗

Near-infrared extinction due to cool supernova dust in Cassiopeia A

We present the results of extinction measurements toward the main ejecta shell of the Cassiopeia A supernova (SN) remnant using the flux ratios between the two near-infrared (NIR) [Fe ii] lines at 1.26 and 1.64 μm. We find a clear correlation between the NIR extinction (E(J−H)) and the radial velocity of ejecta knots, showing that redshifted knots are systematically more obscured than blueshifted ones. This internal “self-extinction” strongly indicates that a large amount of SN dust resides inside and around the main ejecta shell. At one location in the southern part of the shell, we measure E(J−H) by the SN dust of 0.23 ± 0.05 mag. By analyzing the spectral energy distribution of thermal dust emission at that location, we show that there are warm (∼100 K) and cool (∼40 K) SN dust components and that the latter is responsible for the observed E(J−H). We investigate the possible grain species and size of each component and find that the warm SN dust needs to be silicate grains such as MgSiO{sub 3}, Mg{sub 2}SiO{sub 4}, and SiO{sub 2}, whereas the cool dust could be either small (≲0.01 μm) Fe or large (≳0.1 μm) Si grains. We suggest that the warm and cool dust components in Cassiopeia A represent grain species produced in diffuse SN ejecta and in dense ejecta clumps, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Dielectric-barrier-discharge plasma-assisted hydrogen diffusion flame. Part 1: Temperature, oxygen, and fuel measurements by one-dimensional fs/ps rotational CARS imaging

One-dimensional hybrid fs/ps CARS imaging provides single-laser-shot measurements of temperature, oxygen, and hydrogen in a plasma-assisted hydrogen diffusion flame. The coaxial dielectric-barrier-discharge burner collapses the Re ~50 hydrogen diffusion flame to within ~5 mm of the burner surface at an applied AC potential of 8.75 kV at 18 kHz, coinciding nicely with the full spatial extent of the 1D CARS measurements. Translating the burner through the measurement volume allowed for measurements at numerous radial locations in increments of 1 mm with a resolution of 140 µm × 30 µm × 600 µm, sufficient to resolve spatial gradients in this unsteady flame. Longer probe delays, required for improved dynamic range in regions of high temperature fluctuations, proved difficult to model as a result of a nontrivial decay in the O 2 Raman coherence arising from complexities associated with the triplet ground electronic state of the O 2 molecule. Oxygen linewidths were treated empirically using the observed O 2 coherence decay in spectra acquired from the product gases of lean, near-adiabatic H 2 /air flames stabilized on a Hencken flat-flame burner. While still leading to errors up to 10% at worst, the empirically determined Raman linewidth factors eliminated any systematic error in the O 2 /N 2 measurements with probe delay. Temperature measurements in the Hencken Burner flames proved to be insensitive to probe pulse delay, providing robust thermometry. Here, demonstration of this technique in both the canonical Hencken burner flames and a new DBD burner validates its effectiveness in producing multiple spatially resolved measurements in combustion environments. Measurements in the DBD burner revealed an unsteady, counterflow flattened flame structure near the fuel orifice which became unsteady as the reaction zone curves towards the surface for larger radial positions. Lastly, fluctuations in the fuel concentration were largest at the source, as the large, plasma-generated, unsteady external toroidal vortex that dominates the transport in this flame provides enhanced ventilation of the flame surface in close proximity to the fuel tube.

42 ENGINEERING↗

Evolution of plasma loops in a semi-toroidal pinch experiment

The FlareLab experiment is a pulsed-power discharge generating magnetized plasma loops similar to a pinch experiment in a semi-toroidal configuration. After gas breakdown along a circular magnetic guide field, the structure expands in its major radius as the plasma becomes highly conductive and the discharge current rises. Photographs, current and electron density measurements reveal a significant broadening in the lateral direction leading to an increasing departure from radial symmetry of plasma parameters in the cross section. It is shown that the luminosity is related to both high electron density and high current density. Simultaneous measurements of current density and electric field reveal a high parallel resistivity of the plasma leading to fast diffusion across the magnetic field. Indications for anomalous resistivity are found by comparison with the Spitzer formula. As the experiment differs from a z-pinch experiment only by the semi-circular shape of the current path, the observed evolution is unexpected and might be of more fundamental significance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First principles-based multiparadigm, multiscale strategy for simulating complex materials processes with applications to amorphous SiC films

Progress has recently been made in developing reactive force fields to describe chemical reactions in systems too large for quantum mechanical (QM) methods. In particular, ReaxFF, a force field with parameters that are obtained solely from fitting QM reaction data, has been used to predict structures and properties of many materials. Important applications require, however, determination of the final structures produced by such complex processes as chemical vapor deposition, atomic layer deposition, and formation of ceramic films by pyrolysis of polymers. This requires the force field to properly describe the formation of other products of the process, in addition to yielding the final structure of the material. We describe a strategy for accomplishing this and present an example of its use for forming amorphous SiC films that have a wide variety of applications. Extensive reactive molecular dynamics (MD) simulations have been carried out to simulate the pyrolysis of hydridopolycarbosilane. The reaction products all agree with the experimental data. After removing the reaction products, the system is cooled down to room temperature at which it produces amorphous SiC film, for which the computed radial distribution function, x-ray diffraction pattern, and the equation of state describing the three main SiC polytypes agree with the data and with the QM calculations. Extensive MD simulations have also been carried out to compute other structural properties, as well the effective diffusivities of light gases in the amorphous SiC film.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conditions for acceleration of energetic ions greater than 30 keV associated with the earth's bow shock

A statistical analysis of particles (greater than 30 keV/charge) upstream of the earth's bow shock is conducted and shows that the rate of occurrence of upstream particle events is relative to the angle between the magnetic field and the shock normal at the shock intersection point as well as relative to the angle between the magnetic field and the radial direction (i.e., the sun-earth line). In addition, the occurrence rate of upstream particle events relative to the bow shock connection time of a field line convected with the solar wind is presented for a model bow shock. A linear dependence of the diffusion coefficient on energy per charge is apparent with the value of the mean free path of a 30-keV proton found to be about 4 earth radii, and the free escape boundary to be at about 30 earth radii in front of the bow shock.

Scholer, M.↗

Forced cocurrent smoldering combustion

An analytical model of cocurrent smoldering combustion through a very porous solid fuel is developed. Smoldering is initiated at the top of a long radially insulated uniform fuel cylinder, so that the smolder wave propagates downward, opposing an upward-forced flow of oxidizer, with the solid fuel and the gaseous oxidizer entering the reaction zone from the same direction (hence, cocurrent). Radiative heat transfer was incorporated using a diffusion approximation, and smoldering was modeled using a one-step reaction mechanism. The results indicate that, for a given fuel, the final temperature depends only on the initial oxygen mass flux, increasing logarithmically with the mass flux. The smolder velocity is linearly dependent on the initial oxygen mass flux, and, at a fixed value of the flux, increases with initial oxygen mass fraction. The mathematical relationship determining the conditions for steady smolder propagation is presented.

Dosanjh, Sudip S.↗

The Acceleration of Energetic Particles at Coronal Shocks and Emergence of a Double Power-law Feature in Particle Energy Spectra

We present numerical modeling of particle acceleration at coronal shocks propagating through a streamer-like magnetic field by solving the Parker transport equation with spatial diffusion both along and across the magnetic field. We show that the location on the shock where the high-energy particle intensity is the largest, depends on the energy of the particles and on time. The acceleration of particles to more than 100 MeV mainly occurs in the shock-streamer interaction region, due to perpendicular shock geometry and the trapping effect of closed magnetic fields. A comparison of the particle spectra to that in a radial magnetic field shows that the intensity at 100 MeV (200 MeV) is enhanced by more than one order (two orders) of magnitude. This indicates that the streamer-like magnetic field can be an important factor in producing large solar energetic particle events. We also show that the energy spectrum integrated over the simulation domain consists of two different power laws. Further analysis suggests that it may be a mixture of two distinct populations accelerated in the streamer and open field regions, where the acceleration rate differs substantially. Our calculations also show that the particle spectra are affected considerably by a number of parameters, such as the streamer tilt angle, particle spatial diffusion coefficient, and shock compression ratio. While the low-energy spectra agree well with standard diffusive shock acceleration theory, the break energy ranges from ~1 MeV to ~90 MeV and the high-energy spectra can extend to ~1 GeV with a slope of ~2–3.

97 MATHEMATICS AND COMPUTING↗

Keyhole fluctuation and pore formation mechanisms during laser powder bed fusion additive manufacturing

Keyhole porosity is a key concern in laser powder-bed fusion (LPBF), potentially impacting component fatigue life. However, some keyhole porosity formation mechanisms, e.g., keyhole fluctuation, collapse and bubble growth and shrinkage, remain unclear. Using synchrotron X-ray imaging we reveal keyhole and bubble behaviour, quantifying their formation dynamics. The findings support the hypotheses that: (i) keyhole porosity can initiate not only in unstable, but also in the transition keyhole regimes created by high laser power-velocity conditions, causing fast radial keyhole fluctuations (2.5-10 kHz); (ii) transition regime collapse tends to occur part way up the rear-wall; and (iii) immediately after keyhole collapse, bubbles undergo rapid growth due to pressure equilibration, then shrink due to metal-vapour condensation. Concurrent with condensation, hydrogen diffusion into the bubble slows the shrinkage and stabilises the bubble size. The keyhole fluctuation and bubble evolution mechanisms revealed here may guide the development of control systems for minimising porosity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spatiotemporal dynamics of ionic reorganization near biological membrane interfaces

Electrical signals in excitable cells involve spatially localized ionic fluxes through ion channels and pumps on cellular lipid membranes. Common approaches to understand how these fluxes spread assume that the membrane and the surrounding electrolyte comprise an equivalent circuit of capacitors and resistors, which ignores the localized nature of transmembrane ion transport, the resulting ionic gradients and electric fields, and their spatiotemporal relaxation. Here, we consider a model of localized ion pumping across a lipid membrane, and use theory and simulation to investigate how the electrochemical signal propagates spatiotemporally in and out of plane along the membrane. The localized pumping generates long-ranged electric fields with three distinct scaling regimes along the membrane: a constant potential near-field region, an intermediate monopolar region, and a far-field dipolar region. Upon sustained pumping, the monopolar region expands radially in plane with a steady speed that is enhanced by the dielectric mismatch and the finite thickness of the lipid membrane. For unmyelinated lipid membranes in physiological settings, we find remarkably fast propagation speeds of ∼ 40 m / s , allowing faster ionic reorganization compared to bare diffusion. Together, our paper shows that transmembrane ionic fluxes induce transient long-ranged electric fields in electrolyte solutions, which may play hitherto unappreciated roles in biological signaling. Published by the American Physical Society 2025

Row, Hyeongjoo (ORCID:000000033623512X)↗

The spiral wave instability induced by a giant planet. I. Particle stirring in the inner regions of protoplanetary disks

We have recently shown that spiral density waves propagating in accretion disks can undergo a parametric instability by resonantly coupling with and transferring energy into pairs of inertial waves (or inertial-gravity waves when buoyancy is important). In this paper, we perform inviscid three-dimensional global hydrodynamic simulations to examine the growth and consequence of this instability operating on the spiral waves driven by a Jupiter-mass planet in a protoplanetary disk. We find that the spiral waves are destabilized via the spiral wave instability (SWI), generating hydrodynamic turbulence and sustained radially alternating vertical flows that appear to be associated with long wavelength inertial modes. In the interval 0.3 R{sub p}⩽R⩽0.7R{sub p}, where R {sub p} denotes the semimajor axis of the planetary orbit (assumed to be 5 au), the estimated vertical diffusion rate associated with the turbulence is characterized by α{sub diff}∼ (0.2--1.2)×10{sup −2}. For the disk model considered here, the diffusion rate is such that particles with sizes up to several centimeters are vertically mixed within the first pressure scale height. This suggests that the instability of spiral waves launched by a giant planet can significantly disperse solid particles and trace chemical species from the midplane. In planet formation models where the continuous local production of chondrules/pebbles occurs over Myr timescales to provide a feedstock for pebble accretion onto these bodies, this stirring of solid particles may add a time constraint: planetary embryos and large asteroids have to form before a gas giant forms in the outer disk, otherwise the SWI will significantly decrease the chondrule/pebble accretion efficiency.

79 ASTRONOMY AND ASTROPHYSICS↗

Diffusive redistribution of water vapor in the solar nebula revisited

Stevenson and Lunine presented a model for enhancing the abundance of solid material in the region of the solar nebula at the water condensation point. This was used to provide a means to produce a much more rapid formation of Jupiter than the standard solar nebula models. However, they underestimated the drag induced sun-ward radial drift of the planetesimals of interest. Reanalysis reveals that these particles would spread over the inner solar system and might influence the formation of the asteroids.

Sears, William D.↗

SELF-GRAVITY, RESONANCES AND ORBITAL DIFFUSION IN STELLAR DISKS

Fluctuations in a stellar system's gravitational field cause the orbits of stars to evolve. The resulting evolution of the system can be computed with the orbit-averaged Fokker–Planck equation once the diffusion tensor is known. We present the formalism that enables one to compute the diffusion tensor from a given source of noise in the gravitational field when the system's dynamical response to that noise is included. In the case of a cool stellar disk we are able to reduce the computation of the diffusion tensor to a one-dimensional integral. We implement this formula for a tapered Mestel disk that is exposed to shot noise and find that we are able to explain analytically the principal features of a numerical simulation of such a disk. In particular the formation of narrow ridges of enhanced density in action space is recovered. As the disk's value of Toomre's Q is reduced and the disk becomes more responsive, there is a transition from a regime of heating in the inner regions of the disk through the inner Lindblad resonance to one of radial migration of near-circular orbits via the corotation resonance in the intermediate regions of the disk. The formalism developed here provides the ideal framework in which to study the long-term evolution of all kinds of stellar disks.

79 ASTRONOMY AND ASTROPHYSICS↗