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 433 records · Page 24

A numerical simulation of solar energetic particle dropouts during impulsive events

This paper investigates the conditions for producing rapid variations of solar energetic particle (SEP) intensity commonly known as 'dropouts'. In particular, we use numerical model simulations based on solving the focused transport equation in the three-dimensional Parker interplanetary magnetic field to put constraints on the properties of particle transport coefficients in both directions perpendicular and parallel to the magnetic field. Our calculations of the temporal intensity profile of 0.5 and 5 MeV protons at the Earth show that the perpendicular diffusion must be small while the parallel mean free path is long in order to reproduce the phenomenon of SEP dropouts. When the parallel mean free path is a fraction of 1 AU and the observer is located at 1 AU, the perpendicular to parallel diffusion ratio must be below 10{sup –5} if we want to see the particle flux dropping by at least several times within 3 hr. When the observer is located at a larger solar radial distance, the perpendicular to parallel diffusion ratio for reproducing the dropouts should be even lower than that in the case of 1 AU distance. A shorter parallel mean free path or a larger radial distance from the source to observer will cause the particles to arrive later, making the effects of perpendicular diffusion more prominent and SEP dropouts disappear. All of these effects require the magnetic turbulence that resonates with the particles to be low everywhere in the inner heliosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Energetic charged particles in the magnetosphere of Neptune

The Voyager 2 cosmic ray system (CRS) measured significant fluxes of energetic greater than or equal to about 1 megaelectron volt (MeV) trapped electrons and protons in the magnetosphere of Neptune. The intensities are maximum near a magnetic L shell of 7, decreasing closer to the planet because of absorption by satellites and rings. In the region of the inner satellites of Neptune, the radiation belts have a complicated structure, which provides some constraints on the magnetic field geometry of the inner magnetosphere. Electron phase-space densities have a positive radial gradient, indicating that they diffuse inward from a source in the outer magnetosphere. Electron spectra from 1 to 5 MeV are generally well represented by power laws with indices near 6, which harden in the region of peak flux to power law indices of 4 to 5. Protons have significantly lower fluxes than electrons throughout the magnetosphere, with large anisotropies due to radial intensity gradients. The radiation belts resemble those of Uranus to the extent allowed by the different locations of the satellites, which limit the flux at each planet.

Stone, E. C.↗

Challenges and approaches to interpretive modeling of boundary plasma and neutral transport in a closed, pumped divertor

An experimental discharge from the DIII-D tokamak is modeled using the SOLPS-ITER code suite and compared against measurements in the pumped and relatively closed upper divertor. Uncertainties of boundary plasma simulations are identified by attempting to match code inputs to experimental conditions, including iteratively solving transport coefficients to match upstream experimental profiles using varying quantities of core particle flux, different pumping models, and various assumptions of ion thermal transport. Simulated boundary conditions for particle injection at the core interface are shown to be relevant to the plasma solution at the divertor targets, even if upstream transport is modified so that plasma profiles are comparatively similar, although seperatrix density is not held constant. When upstream plasma profiles are matched to experimental measurements by varying diffusive transport coefficients, using either poloidally symmetric or ballooning structure, the model finds a majority of injected energy being transported radially off the computational domain, in conflict with experimental radiated power measurements and heat flux measurements at the divertor target. Imposing a maximum thermal diffusivity or radially shifting the experimental separatrix location of the fitted profiles to increase power conducted to the targets by increasing the upstream electron temperature does not significantly modify this result. Including a thermalizing plenum volume in the simulation domain is shown to maintain the experimental volumetric pumping rate without knowing the neutral energy distribution incident on the pump duct a priori. By modifying transport parameters to match different assumptions for ion temperature, downstream neutral pressure changes by more than a factor of two, suggesting that attention to ion thermal transport may be a critical parameter for simulations to accurately resolve recycling and neutral transport, particularly in a closed divertor geometry. In addition to quantifying various modeling uncertainties, this work motivates both further experimental study and modeling improvements to improve predictive capabilities.

divertor↗

Diffuse ionized gas in nearby Galaxies

We discuss the distribution and spectral characteristics of diffuse ionized gas in nearby galaxies. The diffuse gas has been mapped in both edge-on and more face-on galaxies providing information on the radial and vertical distributions. Emission from diffuse ionized gas is strongest around star forming regions. The vertical distribution appears related to the radio continuum thick-disk emission. We also briefly discuss ionization mechanisms and the connection between star formation characteristics and morphology of the interstellar medium.

Walterbos, R. A. M.↗

Experimental determination of the velocity and strain rate field in a laminar H2/Air counter-flow diffusion flame via LDA

Measurements of the axial and radial components of velocity on the air side of stagnation in an axisymmetric H2/Air laminar counter-flow diffusion flame are reported. Results include the two-dimensional velocity field and computed velocity gradients (strain rates) along the stagnation streamline at two 'characteristic' strain rates, below the extinction limit. The measurements generally verify the modeling assumptions appropriate to the model of Kee et al. (1988). The 'traditional' potential flow model is not consistent with the measured results.

Yeo, S. H.↗

Experimental investigation of rotor-stator interaction in a centrifugal pump with several vaned diffusers

Steady and unsteady diffuser vane pressure measurements have been conducted with a two-dimensional test impeller, in an experimental investigation of rotor-stator interaction within a centrifugal pump having several vaned diffusers, under conditions of different flow coefficients and different radial gaps between the impeller blade trailing edge and the diffuser vane leading edge. The largest pressure fluctuations on the diffuser vanes and the impeller blades were found to be of the same order of magnitude as the total pressure rise across the pump. Increasing the number of diffuser vanes was found to result in a significant decrease of impeller blade pressure fluctuations.

Arndt, N.↗

Avalanche statistics of fluctuation-induced fluxes from the SLPM and the W7-AS stellarator

Measurements of fluctuating floating potentials and ion saturation currents at different radial locations in the Santander Linear Plasma Machine (Castellanos et al 2005 Plasma Phys. Control. Fusion47 2067) and at the edge of the W7-AS stellarator by means of radially movable Langmuir probes allow to infer the corresponding fluctuation-induced radial flux temporal series. Avalanche-like transport events are identified in the time series and statistically characterized in terms of avalanche size/duration/quiet-time distributions and size-duration scaling relations. Transport is diffusive in the inner and intermediate radial region of the SLPM r < r tr ≈ 2.6 cm, undergoing a transition at r tr , becoming non-diffusive in the outermost region of the device, r > r tr . Here, the results obtained at the edge of the W7-AS stellarator are similar to those found in SLPM for r > r tr , i.e. consistent with what would be expected for scale-free, self-similar plasma transport dynamics near a critical state.

Avalanches↗

Properties of gravitationally equilibrated Yukawa systems—A molecular dynamics study

Using 2D Molecular Dynamics simulation, the equilibrium and dynamical properties of a gravitationally equilibrated Yukawa liquid are investigated. We observe that due to asymmetry introduced in one direction by gravity, several interesting features arise. For example, for a given value of coupling parameter Γ, screening parameter κ, and according to a chosen value of gravitational force g (say in y-direction), the system is seen to exhibit super-, sub- or normal diffusion. Interestingly, x-averaged density profiles, unlike a barotropic fluid, acquires sharp, free surface with scale free linear y-dependence. As can be expected for a system with macroscopic gradients, self-diffusion calculated from Green-Kubo’s formalism does not agree with that obtained from Einstein-Smoluchowski diffusion. A 2D angular-radial pair correlation function g(r, θ) clearly indicates asymmetric features induced by gravity. We observe that due to compression in y-direction, though in liquid state for all values of gravity considered, the transverse mode is found to predominant as compared to the longitudinal mode, leading to a novel Anisotropic Solid-like Yukawa liquid.

74 ATOMIC AND MOLECULAR PHYSICS↗

A NUMERICAL SIMULATION OF COSMIC-RAY MODULATION NEAR THE HELIOPAUSE

Based on a hybrid galactic cosmic-ray transport model, which incorporated MHD global heliospheric data into Parker’s cosmic-ray transport equation, we studied the behavior of the transport of galactic cosmic rays and the corresponding gradients in their flux near the heliopause (HP). We found that, (1) by increasing the ratio of the parallel diffusion coefficient to the perpendicular diffusion coefficient in the interstellar magnetic field of the outer heliosheath, the simulated radial flux near the HP increases as well. As the ratio multiplying factor reached 10{sup 10}, the radial flux experienced a sudden jump near the HP, similar to what Voyager 1 observed in 2012. (2) The effect of changing the diffusion coefficients’ ratio on the radial flux variation depends on the energy of the cosmic rays, the lower the energy, the more pronounced the effect is. (3) The magnitude of the diffusion coefficients also affect the radial flux near the HP, the modulation beyond the HP varies by adjusting the magnitude multiplying factor.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical models of solar modulation of galactic cosmic rays, including drifts

Numerical solutions are presented of the steady-state modulation of galactic cosmic rays for a realistic two-dimensional model which includes drifts. Calculations of energy spectra and the spatial density variation of protons are shown for a model that incorporates the solar-minimum magnetic field and nominal interplanetary propagation parameters. The calculated values are relatively insensitive to the assumed size of the modulating region and to the magnitude of the particle diffusion coefficient. The solutions show a broad interior 'plateau' in which radial gradients are small. A simple model that neglects diffusion but includes drifts and adiabatic cooling accounts for these results. It is concluded that in this model and the range of parameters, diffusion plays a role secondary to that of drifts and adiabatic cooling.

Jokipii, J. R.↗

An experimental investigation of two large annular diffusers with swirling and distorted inflow

Two annular diffusers downstream of a nacelle-mounted fan were tested for aerodynamic performance, measured in terms of two static pressure recovery parameters (one near the diffuser exit plane and one about three diameters downstream in the settling duct) in the presence of several inflow conditions. The two diffusers each had an inlet diameter of 1.84 m, an area ratio of 2.3, and an equivalent cone angle of 11.5, but were distinguished by centerbodies of different lengths. The dependence of diffuser performance on various combinations of swirling, radially distorted, and/or azimuthally distorted inflow was examined. Swirling flow and distortions in the axial velocity profile in the annulus upstream of the diffuser inlet were caused by the intrinsic flow patterns downstream of a fan in a duct and by artificial intensification of the distortions. Azimuthal distortions or defects were generated by the addition of four artificial devices (screens and fences). Pressure recovery data indicated beneficial effects of both radial distortion (for a limited range of distortion levels) and inflow swirl. Small amounts of azimuthal distortion created by the artificial devices produced only small effects on diffuser performance. A large artificial distortion device was required to produce enough azimuthal flow distortion to significantly degrade the diffuser static pressure recovery.

Eckert, W. T.↗

The Quest for Accurate Liquid Water Properties from First Principles

Developing accurate ab initio molecular dynamics (AIMD) models that capture both electronic reorganization and nuclear quantum effects associated with hydrogen bonding is key to quantitative understanding of bulk water and its anomalies, as well as its role as a universal solvent. For condensed phase simulations, AIMD has typically relied on the generalized gradient approximation (GGA) of Density Functional Theory (DFT) as the underlying model chemistry for the potential energy surface, with nuclear quantum effects (NQEs) sometimes modeled by performing classical molecular dynamics simulations at elevated temperatures. Here we show that the properties of liquid water obtained from the meta-GGA B97M-rV functional, when evaluated using accelerated path integral molecular dynamics simulations, displays accuracy comparable to a computationally expensive dispersion corrected hybrid functional, revPBE0-D 3 . We show that the meta-GGA DFT reproduces bulk water properties including radial distribution functions, self-diffusion coefficients, and infrared (IR) spectra with comparable accuracy of a much more expensive hybrid functional. In conclusion, this work demonstrates that the underlying quality of a good DFT functional requires evaluation with quantum nuclei, and that high temperature simulations are a poor proxy for properly treating nuclear quantum effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio molecular dynamics simulations of liquid water using high quality meta-GGA functionals

We have used ab initio molecular dynamics (AIMD) to characterize water properties using two meta-generalized gradient approximation (meta-GGA) functionals, M06-L-D3 and B97M-rV, and compared their performance against a standard GGA corrected for dispersion, revPBE-D3, at ambient conditions (298 K, and 1 g cm –3 or 1 atm). Simulations of the equilibrium density, radial distribution functions, self-diffusivity, the infrared spectrum, liquid dipole moments, and characterizations of the hydrogen bond network show that all three functionals have overcome the problem of the early AIMD simulations that erroneously found ambient water to be highly structured, but they differ substantially among themselves in agreement with experiment on this range of water properties. We show directly using water cluster data up through the pentamer that revPBE-D3 benefits from a cancellation of its intrinsic functional error by running classical trajectories, whereas the meta-GGA functionals are demonstrably more accurate and would require the simulation of nuclear quantum effects to realize better agreement with all cluster and condensed phase properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Communication: Quantum molecular dynamics simulation of liquid para-hydrogen by nuclear and electron wave packet approach

Liquid para-hydrogen (p-H{sub 2}) is a typical quantum liquid which exhibits strong nuclear quantum effects (NQEs) and thus anomalous static and dynamic properties. We propose a real-time simulation method of wave packet (WP) molecular dynamics (MD) based on non-empirical intra- and inter-molecular interactions of non-spherical hydrogen molecules, and apply it to condensed-phase p-H{sub 2}. The NQEs, such as WP delocalization and zero-point energy, are taken into account without perturbative expansion of prepared model potential functions but with explicit interactions between nuclear and electron WPs. The developed MD simulation for 100 ps with 1200 hydrogen molecules is realized at feasible computational cost, by which basic experimental properties of p-H{sub 2} liquid such as radial distribution functions, self-diffusion coefficients, and shear viscosities are all well reproduced.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum mechanical force field for hydrogen fluoride with explicit electronic polarization

The explicit polarization (X-Pol) theory is a fragment-based quantum chemical method that explicitly models the internal electronic polarization and intermolecular interactions of a chemical system. X-Pol theory provides a framework to construct a quantum mechanical force field, which we have extended to liquid hydrogen fluoride (HF) in this work. The parameterization, called XPHF, is built upon the same formalism introduced for the XP3P model of liquid water, which is based on the polarized molecular orbital (PMO) semiempirical quantum chemistry method and the dipole-preserving polarization consistent point charge model. We introduce a fluorine parameter set for PMO, and find good agreement for various gas-phase results of small HF clusters compared to experiments and ab initio calculations at the M06-2X/MG3S level of theory. In addition, the XPHF model shows reasonable agreement with experiments for a variety of structural and thermodynamic properties in the liquid state, including radial distribution functions, interaction energies, diffusion coefficients, and densities at various state points.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using stepped anvils to make even insulation layers in laser-heated diamond-anvil cell samples

Here, we describe a method to make even insulation layers for high-pressure laser-heated diamond-anvil cell samples using stepped anvils. Moreover, the method works for both single-sided and double-sided laser heating using solid or fluid insulation. The stepped anvils are used as matched pairs or paired with a flat culet anvil to make gasket insulation layers and not actually used at high pressures; thus, their longevity is ensured. We also compare the radial temperature gradients and Soret diffusion of iron between self-insulating samples and samples produced with stepped anvils and find that less pronounced Soret diffusion occurs in samples with even insulation layers produced by stepped anvils.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Striations in moderate pressure dc driven nitrogen glow discharge

Plasma stratification has been studied for more than a century. Despite the many experimental studies reported on this topic, theoretical analyses and numerical modeling of this phenomenon have been mostly limited to rare gases. In this work, a one-dimensional fluid model with detailed kinetics of electrons and vibrationally excited molecules is employed to simulate moderate-pressure (i.e. a few Torrs) dc discharge in nitrogen in a 15.5 cm long tube of radius 0.55 cm. The model also considers ambipolar diffusion to account for the radial loss of ions and electrons to the wall. The proposed model predicts self-excited standing striations in nitrogen for a range of discharge currents. The impact of electron transport parameters and reaction rates obtained from a solution of local two-term and a multi-term Boltzmann equation on the predictions are assessed. In-depth kinetic analysis indicates that the striations result from the undulations in electron temperature caused due to the interaction between ionization and vibrational reactions. Furthermore, the vibrationally excited molecules associated with the lower energy levels are found to influence nitrogen plasma stratification and the striation pattern strongly. A balance between ionization processes and electron energy transport allows the formation of the observed standing striations. Simulations were conducted for a range of discharge current densities from ~0.018 to 0.080 mA cm –2 , for an operating pressure of 0.7 Torr. Parametric studies show that the striation length decreases with increasing discharge current. Furthermore, the predictions from the model are compared against experimental measurements and are found to agree favorably.

42 ENGINEERING↗

FY24 Progress Report on Viscosity and Thermal Conductivity Measurements of Nuclear Industry Relevant Chloride Salts: An Experimental and Computational Study

As presented in this report, experimental and computational techniques were performed to assess the viscosity and thermal conductivity of key alkali and actinide chloride mixtures for molten salt reactor developers. These mixtures were pure LiCl, NaCl-KCl, LiCl-NaCl, LiCl-KCl, LiCl-NaCl-KCl, and NaCl-UCl 3 . Experimental measurements of viscosity were performed with a rolling ball viscometer, whereas experimental measurements of thermal conductivity were performed with a variable gap apparatus. Additional benchmarking work was performed using both property measurement systems to prepare for x-ray radiography in stainless-steel crucibles for viscosity and to ensure that calibration methods were accurate for thermal conductivity before assessing the NaCl-UCl 3 system. Validation data for the NaCl-UCl 3 in literature are minimal. Details on the calibration methods, salt measurement processes, and sources of error and uncertainty are discussed in detail for both property measurements. The computational methods described herein involved ab-initio molecular dynamics (AIMD) calculations using CP2K. The calculations were performed for the LiCl-KCl-NaCl and NaCl-UCl 3 systems. These calculations not only provided thermophysical property estimations for comparison to experimental data, but they also allowed for the determination of diffusion coefficients, coordination numbers, and radial distribution functions to provide insight into ion mobility and local coordination environments, which is linked to macroscopic property trends.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗