Results for “radial diffusion” — Page 34

DOE OSTI2020

Numerical predictions of cryogenic hydrogen vertical jets

Here, a comparison of Computational Fluid Dynamics (CFD) predictions with measurements is presented for cryo-compressed hydrogen vertical jets. The stagnation conditions of the experiments are characteristic of unintended leaks from pipe systems that connect cryogenic hydrogen storage tanks and could be encountered at a fuel cell refueling station. Jets with pressure up to 5 bar and temperatures just above the saturation liquid temperature were examined. Comparisons are made to the centerline mass fraction and temperature decay rates, the radial profiles of mass fraction and the contours of volume fraction. Two notional nozzle approaches are tested to model the under-expanded jet that was formed in the tests with pressures above 2 bar. In both approaches the mass and momentum balance from the throat to the notional nozzle are solved, while the temperature at the notional nozzle was assumed equal to the nozzle temperature in the first approach and was calculated by an energy balance in the second approach. The two approaches gave identical results. Satisfactory agreement with the measurements was found in terms of centerline mass fraction and temperature. However, for test with 3 and 4 bar release the concentration was overpredicted. Furthermore, a wider radial spread was observed in the predictions possibly revealing higher degree of diffusion using the k-ε turbulence model. An integral model for cryogenic jets was also developed and provided good results. Finally, a test simulation was performed with an ambient temperature jet and compared to the cold jet showing that warm jets decay faster than cold jets.

08 HYDROGEN↗

DOE OSTI2021

Turbulent impurity transport simulations in Wendelstein 7-X plasmas

A study of turbulent impurity transport by means of quasilinear and nonlinear gyrokinetic simulations is presented for Wendelstein 7-X (W7-X). The calculations have been carried out with the recently developed gyrokinetic code stella. Different impurity species are considered in the presence of various types of background instabilities: ion temperature gradient (ITG), trapped electron mode (TEM) and electron temperature gradient (ETG) modes for the quasilinear part of the work; ITG and TEM for the nonlinear results. While the quasilinear approach allows one to draw qualitative conclusions about the sign or relative importance of the various contributions to the flux, the nonlinear simulations quantitatively determine the size of the turbulent flux and check the extent to which the quasilinear conclusions hold. Although the bulk of the nonlinear simulations are performed at trace impurity concentration, nonlinear simulations are also carried out at realistic effective charge values, in order to know to what degree the conclusions based on the simulations performed for trace impurities can be extrapolated to realistic impurity concentrations. The presented results conclude that the turbulent radial impurity transport in W7-X is mainly dominated by ordinary diffusion, which is close to that measured during the recent W7-X experimental campaigns. Finally, it is also confirmed that thermodiffusion adds a weak inward flux contribution and that, in the absence of impurity temperature and density gradients, ITG- and TEM-driven turbulence push the impurities inwards and outwards, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Plasma turbulence in the scrape-off layer of tokamak devices

Plasma turbulence is explored in the scrape-off layer of tokamak devices using three-dimensional global two-fluid simulations. Two transport regimes are discussed: one in which the turbulent fluctuations saturate nonlinearly due to the Kelvin-Helmholtz instability, and another in which the fluctuations saturate due to a local flattening of the plasma gradients and associated removal of the linear instability drive. Focusing on the latter regime, analytical estimates of the cross-field transport and plasma profile gradients are obtained that display Bohm-scaling diffusion properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Binding energies and spatial structures of small carrier complexes in monolayer transition-metal dichalcogenides via diffusion Monte Carlo

Ground-state diffusion Monte Carlo is used to investigate the binding energies and intercarrier radial probability distributions of excitons, trions, and biexcitons in a variety of two-dimensional transition-metal dichalcogenide materials. We compare these results to approximate variational calculations, as well as to analogous Monte Carlo calculations performed with simplified carrier interaction potentials. Our results highlight the successes and failures of approximate approaches as well as the physical features that determine the stability of small carrier complexes in monolayer transition-metal dichalcogenide materials. In conclusion, we discuss points of agreement and disagreement with recent experiments.

36 MATERIALS SCIENCE↗

NASA NTRS1968

The solar radial gradient of galactic cosmic rays.

Steady solar radial gradient of galactic cosmic ray intensity in interplanetary region noting convection, diffusion, energy losses and differential streaming

Axford, W. I.↗

NASA NTRS1989

Magnetic fields in the accretion disks of cataclysmic variables

The radial transport of magnetic flux in CV accretion disks is considered. Turbulent diffusion in a disk with finite outer radius Rd leads to decay of large-scale magnetic fields, as magnetic flux can leak out at the outer edge of the disk. Numerical computations of the decay rate and magnetic structure are presented. It is shown that the decay time is significantly shorter than the accretion time. Centrifugally driven winds may be possible from the outer parts of the disk, provided there exist efficient dynamo processes which regenerate the field.

Van Ballegooijen, A. A.↗

NASA NTRS1995

Rapid Numerical Simulation of Viscous Axisymmetric Flow Fields

A two-dimensional Navier-Stokes code has been developed for rapid numerical simulation of axisymmetric flow fields, including flow fields with an azimuthal velocity component. The azimuthal-invariant Navier-Stokes equations in a cylindrical coordinate system are mapped to a general body-fitted coordinate system, with the streamwise viscous terms then neglected by applying the thin-layer approximation. Turbulence effects are modeled using an algebraic model, typically the Baldwin-Lomax turbulence model, although a modified Cebeci-Smith model can also be used. The equations are discretized using central finite differences and solved using a multistage Runge-Kutta algorithm with a spatially varying time step and implicit residual smoothing. Results are presented for calculations of supersonic flow over a waisted body-of-revolution, transonic flow through a normal shock wave in a straight circular duct of constant cross sectional area, swirling supersonic (inviscid) flow through a strong shock in a straight radial duct, and swirling subsonic flow in an annular-to-circular diffuser duct. Comparisons between computed and experimental results are in fair to good agreement, demonstrating that the viscous code can be a useful tool for practical engineering design and analysis work.

Tweedt, Daniel L.↗

DOE OSTI2017

Confocal Raman Microscopy Investigation of Molecular Transport into Individual Chromatographic Silica Particles

Porous silica is used as a support in a variety of separation processes, including chromatographic separation and solid-phase extraction. The resolution and efficiency of these applications is significantly impacted by the kinetics of partitioning and molecular transport into the interior of the porous particles. Molecular transport in porous silica has been explored previously by measuring chromatographic elution profiles, but such measurements are limited to relatively low retention conditions, where within-particle molecular transport must be inferred from elution profiles of solutes emerging from a packed column. In this work, a measurement of within-particle molecular transport is carried out using confocal Raman microscopy to probe the time-dependent accumulation of pyrene from an aqueous mobile phase into the center of individual C 18 -chromatographic particles. The measured time constants for pyrene accumulation were much slower than diffusion-limited transport of solute in solution to the particle surface. Furthermore, the accumulation into the center of the particle did not show a time-lag characteristic of slow-transport into the particle interior. The exponential rise of pyrene concentration is, however, consistent with first-order Langmuir adsorption kinetics at low surface coverages. The linear dependence of the time-constant on particle radius indicates an adsorption barrier near the outer boundary of the particle, where the accumulation rate depends on flux across the boundary (proportional to the particle area) to satisfy the within-particle capacity at equilibrium (proportional to the particle volume). The pyrene accumulation kinetics into the porous particle, expressed as a heterogeneous rate constant, were nearly 50-times faster than the pyrene adsorption rate at a planar C 18 -functionalized silica surface, which demonstrates the impact of multiple surface encounters within the porous structure leading to much greater capture efficiency compared to a planar surface. Monte Carlo simulations of within-particle pyrene diffusion, with the adsorption efficiency estimated from the planar-surface adsorption rate, predict a diffusion-to-capture distance within the porous particle that is within 40% of that observed in the radial dependence of the pyrene within-particle accumulation results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Mutual diffusion of binary liquid mixtures containing methanol, ethanol, acetone, benzene, cyclohexane, toluene, and carbon tetrachloride

Mutual diffusion coefficients of all 20 binary liquid mixtures that can be formed out of methanol, ethanol, acetone, benzene, cyclohexane, toluene, and carbon tetrachloride without a miscibility gap are studied at ambient conditions of temperature and pressure in the entire composition range. The considered mixtures show a varying mixing behavior from almost ideal to strongly non-ideal. Predictive molecular dynamics simulations employing the Green-Kubo formalism are carried out. Radial distribution functions are analyzed to gain an understanding of the liquid structure influencing the diffusion processes. It is shown that cluster formation in mixtures containing one alcoholic component has a significant impact on the diffusion process. The estimation of the thermodynamic factor from experimental vapor-liquid equilibrium data is investigated, considering three excess Gibbs energy models, i.e., Wilson, NRTL, and UNIQUAC. It is found that the Wilson model yields the thermodynamic factor that best suits the simulation results for the prediction of the Fick diffusion coefficient. Four semi-empirical methods for the prediction of the self-diffusion coefficients and nine predictive equations for the Fick diffusion coefficient are assessed and it is found that methods based on local composition models are more reliable. Finally, the shear viscosity and thermal conductivity are predicted and in most cases favorably compared with experimental literature values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Dynamic of charged planar geometry in tilted and non-tilted frames

We investigate the dynamics of charged planar symmetry with an anisotropic matter field subject to a radially moving observer called a tilted observer. The Einstein-Maxwell field equations are used to obtain a relation between non-tilted and tilted frames and between kinematical and dynamical quantities. Using the Taub mass formalism and conservation laws, two evolution equations are developed to analyze the inhomogeneities in the tilted congruence. It is found that the radial velocity (due to the tilted observer) and the electric charge have a crucial effect on the inhomogeneity factor. Finally, we discuss the stability in the non-tilted frame in the pure diffusion case and examine the effects of the electromagnetic field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1973

Observations of diffuse galactic H alpha and forbidden N II emission

A study has been made of the intensities and radial velocities of faint galactic H alpha and forbidden N II 6584-emission lines not associated with any known bright H II regions. Although some (about 10%) of the observed radiation may be scattered galactic light from bright H II regions in the galactic plane, most of the radiation appears to be produced by an ionized component of the interstellar gas which is distributed throughout the interstellar medium within the three nearby galactic spiral arms. In the local Orion arm the emitting gas has a temperature between 3000 and 8000 K and mean square electron density approximately equal to 0.05 cm to the minus 6th power in the galactic plane, while in the Perseus and Sagittarius arms the mean square electron density is equal to about 0.1 and 0.9 cm to the minus 6th power, respectively.

Reynolds, R. J.↗

DOE OSTI2025

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel Jana Howard1,3, Guang Yang3, Haiyan Zhao1*, Patrick Warren4, Tiankai Yao3, Steven Cavazos4 Elizabeth Sooby Wood4*, Ching-heng Shiau2 1 University of Idaho, Environmental Science, Idaho Falls, Idaho, USA 2 Boise State University, Microscopy and Characterization Suite, Idaho Falls, Idaho, USA 3 Idaho National Laboratory, Idaho Falls, Idaho, USA 4 University of Texas San Antonio, Department of Physics and Astronomy, San Antonio, Texas, USA *Corresponding author: haiyanz@uidaho.edu; elizabeth.soobywood@utsa.edu Tri-structural Isotropic (TRISO) particle fuel is the preferred choice for the newest and next-generation high-temperature nuclear reactors due to its robust construction [1]. The fuel design effectively contains fission products inside the particle at extremely high temperatures [2]. However, the effect of transition metal fission products on fuel performance and structural integrity remains largely unknown, creating uncertainties in predicting fuel performance and potential failure mechanisms [2]. For example, palladium (Pd) has a strong tendency to form intermetallic phases with uranium (U). These new phases tend to be hard and brittle which could lead to fracture of fuel kernel during irradiation [3,4]. Pd is not normally produced in high concentrations in normal fission reactions however, driving the production of Pd into the locality of a Uranium Oxycarbide (UCO) nuclear fuel provides an opportunity to closely observe the diffusion and chemical interactions. This study focuses on detailed transmission electron microscopy characterization of the intermetallic phases formed by the interaction between a UCO kernel and a Pd bar after annealing at 1100 °C for 100 hours. Figure 1 provides an overview of the UCO kernel in contact with the Pd bar. The UCO-Pd sample surface was examined using a Focused Ion Beam (FIB) Quanta 3D FEG FEI in VCD detector mode at 10kV, 83pA to reveal surface features. Several key surface features in the interaction region were observed including: (1) a small area with dendritic microstructure, (2) color contrast near the UCO kernel and Pd boundary, and (3) darker and lighter regions throughout the sample surface. To analyze diffusion behavior, FIB was used to prepare lamellae from five different areas of the UCO-Pd sample, as well as one for the as-received sample. These lamellae were examined using a Scanning Transmission Electron Microscope ThermoFisher Spectra 300 (STEM-Spectra). Energy Dispersive Spectroscopy results confirmed the diffusion between the UCO fuel and Pd bar. Figure 2 shows the identified phases of UC, UO2, and Pd dispersed throughout the interaction region. It was also found that the Pd concentration decreases as the radial distance from the UCO-Pd interface increases. The Pd concentration was 36.93% atomic fraction 11µm from the interaction region then decreased to 7.51% atomic fraction at 257 µm away from the interaction region. Figure 3 shows how concentrations of U and Pd vary across the interaction zone, highlighting the extent of diffusion. This study confirms that Pd interacts with surrounding material to form U-Pd phase and diffusion zones. These diffusion zones vary in composition depending on its radial distance from the point of contact with the Pd bar. These findings contribute to a deeper understanding of the fission product behavior in high temperature nuclear fuels, aiding in the prediction and mitigation of potential fuel degradation mechanisms. A B C Fig. 1. SEM micrographs of UCO fuel kernel in contact with solid Pd bar and the lift out locations. Figure 1A shows UCO fuel kernel in contact with solid Pd bar. Figure 1B shows a higher magnification image of the UCO fuel-Pd interaction zone. Figure 1C shows lift out sites for the lamellae. A B C Fig 2. EDS maps reveal the distribution of U, O, and Pd in location 3. Figure 2A shows the distribution of uranium. Figure 2B shows the distribution of oxygen. Figure 2C shows the distribution of palladium. A B C Fig 3. Palladium and uranium concentration across interaction zone in location 2. Figure 3A shows EDS map of lift out number 2. Figure 3B shows a zoomed in EDS map taken from the interaction zone. Figure 3C shows a line graph of uranium and palladium concentrations across the interaction zone. References: 1. B.E. Wells, N.R. Phillips, K.J. Geelhood. Pacific West Laboratory. (2021). TRISO Fuel: Properties and Failure Modes. https://www.nrc.gov/docs/ML2117/ML21175A152.pdf (Accessed January 16, 2025). 2. American Nuclear Society. TRISO Fuel Development Progresses in INL, ORNL. https://www.gen-4.org/gif/upload/docs/application/pdf/2014-03/nov13nn_fuel_reprint.pdf (Accessed January 16, 2025) 3. Clark R.A., M.A. Con

36 - MATERIALS SCIENCE↗

DOE OSTI2025

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel Jana Howard1,3, Guang Yang3, Haiyan Zhao1*, Patrick Warren4, Tiankai Yao3, Steven Cavazos4 Elizabeth Sooby Wood4*, Ching-heng Shiau2 1 University of Idaho, Environmental Science, Idaho Falls, Idaho, USA 2 Boise State University, Microscopy and Characterization Suite, Idaho Falls, Idaho, USA 3 Idaho National Laboratory, Idaho Falls, Idaho, USA 4 University of Texas San Antonio, Department of Physics and Astronomy, San Antonio, Texas, USA *Corresponding author: haiyanz@uidaho.edu; elizabeth.soobywood@utsa.edu Tri-structural Isotropic (TRISO) particle fuel is the preferred choice for the newest and next-generation high-temperature nuclear reactors due to its robust construction [1]. The fuel design effectively contains fission products inside the particle at extremely high temperatures [2]. However, the effect of transition metal fission products on fuel performance and structural integrity remains largely unknown, creating uncertainties in predicting fuel performance and potential failure mechanisms [2]. For example, palladium (Pd) has a strong tendency to form intermetallic phases with uranium (U). These new phases tend to be hard and brittle which could lead to fracture of fuel kernel during irradiation [3,4]. Pd is not normally produced in high concentrations in normal fission reactions however, driving the production of Pd into the locality of a Uranium Oxycarbide (UCO) nuclear fuel provides an opportunity to closely observe the diffusion and chemical interactions. This study focuses on detailed transmission electron microscopy characterization of the intermetallic phases formed by the interaction between a UCO kernel and a Pd bar after annealing at 1100 °C for 100 hours. Figure 1 provides an overview of the UCO kernel in contact with the Pd bar. The UCO-Pd sample surface was examined using a Focused Ion Beam (FIB) Quanta 3D FEG FEI in VCD detector mode at 10kV, 83pA to reveal surface features. Several key surface features in the interaction region were observed including: (1) a small area with dendritic microstructure, (2) color contrast near the UCO kernel and Pd boundary, and (3) darker and lighter regions throughout the sample surface. To analyze diffusion behavior, FIB was used to prepare lamellae from five different areas of the UCO-Pd sample, as well as one for the as-received sample. These lamellae were examined using a Scanning Transmission Electron Microscope ThermoFisher Spectra 300 (STEM-Spectra). Energy Dispersive Spectroscopy results confirmed the diffusion between the UCO fuel and Pd bar. Figure 2 shows the identified phases of UC, UO2, and Pd dispersed throughout the interaction region. It was also found that the Pd concentration decreases as the radial distance from the UCO-Pd interface increases. The Pd concentration was 36.93% atomic fraction 11µm from the interaction region then decreased to 7.51% atomic fraction at 257 µm away from the interaction region. Figure 3 shows how concentrations of U and Pd vary across the interaction zone, highlighting the extent of diffusion. This study confirms that Pd interacts with surrounding material to form U-Pd phase and diffusion zones. These diffusion zones vary in composition depending on its radial distance from the point of contact with the Pd bar. These findings contribute to a deeper understanding of the fission product behavior in high temperature nuclear fuels, aiding in the prediction and mitigation of potential fuel degradation mechanisms. A B C Fig. 1. SEM micrographs of UCO fuel kernel in contact with solid Pd bar and the lift out locations. Figure 1A shows UCO fuel kernel in contact with solid Pd bar. Figure 1B shows a higher magnification image of the UCO fuel-Pd interaction zone. Figure 1C shows lift out sites for the lamellae. A B C Fig 2. EDS maps reveal the distribution of U, O, and Pd in location 3. Figure 2A shows the distribution of uranium. Figure 2B shows the distribution of oxygen. Figure 2C shows the distribution of palladium. A B C Fig 3. Palladium and uranium concentration across interaction zone in location 2. Figure 3A shows EDS map of lift out number 2. Figure 3B shows a zoomed in EDS map taken from the interaction zone. Figure 3C shows a line graph of uranium and palladium concentrations across the interaction zone. References: 1. B.E. Wells, N.R. Phillips, K.J. Geelhood. Pacific West Laboratory. (2021). TRISO Fuel: Properties and Failure Modes. https://www.nrc.gov/docs/ML2117/ML21175A152.pdf (Accessed January 16, 2025). 2. American Nuclear Society. TRISO Fuel Development Progresses in INL, ORNL. https://www.gen-4.org/gif/upload/docs/application/pdf/2014-03/nov13nn_fuel_reprint.pdf (Accessed January 16, 2025) 3. Clark R.A., M.A. Con

36 - MATERIALS SCIENCE↗

NASA NTRS1992

Nature of fluid flows in differentially heated cylindrical container filled with a stratified solution

Semiconductor crystals such as Hg(1-x)Cd(x)Te grown by unidirectional solidification Bridgmann method have shown compositional segregations in both the axial and radial directions. Due to the wide separation between the liquidus and the solidus of its pseudobinary phase diagram, there is a diffusion layer of higher HgTe content built up in the melt near the melt-solid interface which gives a solute concentration gradient in the axial direction. Because of the higher thermal conductivity in the melt than that in the crystal there is a thermal leakage through the fused silica crucible wall near the melt-solid interface. This gives a thermal gradient in the radial direction. Hart (1971), Thorpe, Hutt and Soulsby (1969) have shown that under such condition a fluid will become convectively unstable as a result of different diffusivities of temperature and solute. It is quite important to understand the effects of this thermosolute convection on the compositional segregation in the unidirectionally solidified crystals. To reach this goal, we start with a simplified problem. We study the nature of fluid flows of a stratified solution in a cylindrical container with a radial temperature gradient. The cylindrical container wall is considered to be maintained at a higher temperature than that at the center of the solution and the solution in the lower gravitational direction has higher solute concentration which decrease linearly to a lower concentration and then remain constant to the top of the solution. The sample solution is taken to be salt water.

Wang, Jai-Ching↗

DOE OSTI2014

Magnetohydrodynamics and deep mixing in evolved stars. I. Two- and three-dimensional analytical models for the asymptotic giant branch

The advection of thermonuclear ashes by magnetized domains emerging near the H shell was suggested to explain asymptotic giant branch (AGB) star abundances. Here we verify this idea quantitatively through exact MHD models. Starting with a simple two-dimensional (2D) geometry and in an inertia frame, we study plasma equilibria avoiding the complications of numerical simulations. We show that below the convective envelope of an AGB star, variable magnetic fields induce a natural expansion, permitted by the almost ideal MHD conditions, in which the radial velocity grows as the second power of the radius. We then study the convective envelope, where the complexity of macroturbulence allows only for a schematic analytical treatment. Here the radial velocity depends on the square root of the radius. We then verify the robustness of our results with 3D calculations for the velocity, showing that for both studied regions the solution previously found can be seen as a planar section of a more complex behavior, in which the average radial velocity retains the same dependency on the radius found in 2D. As a final check, we compare our results to approximate descriptions of buoyant magnetic structures. For realistic boundary conditions, the envelope crossing times are sufficient to disperse in the huge convective zone any material transported, suggesting magnetic advection as a promising mechanism for deep mixing. The mixing velocities are smaller than for convection but larger than for diffusion and adequate for extra mixing in red giants.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Multiple emission components in the Cygnus cocoon detected from Fermi -LAT observations

Star-forming regions may play an important role in the life cycle of Galactic cosmic rays (CRs), notably as home to specific acceleration mechanisms and transport conditions. Gamma-ray observations of Cygnus X have revealed the presence of an excess of hard-spectrum gamma-ray emission, possibly related to a cocoon of freshly accelerated particles. We seek an improved description of the gamma-ray emission from the cocoon using ~13 yr of observations with the Fermi-Large Area Telescope (LAT) and use it to further constrain the processes and objects responsible for the young CR population. We developed an emission model for a large region of interest, including a description of interstellar emission from the background population of CRs and recent models for other gamma-ray sources in the field. Thus, we performed an improved spectro-morphological characterisation of the residual emission including the cocoon. The best-fit model for the cocoon includes two main emission components: an extended component FCES G78.74+1.56, described by a 2D Gaussian of extension r 68 = 4.4° ± 0.1° -0.1° +0.1° and a smooth broken power law spectrum with spectral indices 1.67 ± 0.05 -0.01 +0.02 and 2.12 ± 0.02 -0.01 +0.00 below and above 3.0 ± 0.6 -0.2 +0.0 GeV, respectively; and a central component FCES G80.00+0.50, traced by the distribution of ionised gas within the borders of the photo-dissociation regions and with a power law spectrum of index 2.19 ± 0.03 -0.01 +0.00 that is significantly different from the spectrum of FCES G78.74+1.56. An additional extended emission component FCES G78.83+3.57, located on the edge of the central cavities in Cygnus X and with a spectrum compatible with that of FCES G80.00+0.50, is likely related to the cocoon. For the two brightest components FCES G80.00+0.50 and FCES G78.74+1.56, spectra and radial-azimuthal profiles of the emission can be accounted for in a diffusion-loss framework involving one single population of non-thermal particles with a flat injection spectrum. Particles span the full extent of FCES G78.74+1.56 as a result of diffusion from a central source, and give rise to source FCES G80.00+0.50 by interacting with ionised gas in the innermost region. For this simple diffusion-loss model, viable setups can be very different in terms of energetics, transport conditions, and timescales involved, and both hadronic and leptonic scenarios are possible. The solutions range from long-lasting particle acceleration, possibly in prominent star clusters such as Cyg OB2 and NGC 6910, to a more recent and short-lived release of particles within the last 10–100 kyr, likely from a supernova remnant. The observables extracted from our analysis can be used to perform detailed comparisons with advanced models of particle acceleration and transport in star-forming regions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Effect on plasma rotation of lower hybrid (LH) waves in Alcator C-Mod

The injection of LH waves for current drive into a tokamak changes the ion toroidal rotation. In Alcator C-Mod, the direction of the steady state rotation change due to LH waves depends on the plasma current and the density. The change in rotation can be estimated by balancing the external torque of lower hybrid waves with the turbulent radial transport of the momentum. For high plasma current, the turbulent pinch and diffusion of the injected counter-current momentum are sufficient to explain the rotation change. However, for low plasma current, the change in the the intrinsic momentum transport (residual stress) for a non-rotating state is required to explain the co-current rotation change. Accordingly, we investigate the intrinsic momentum transport for the non-rotating state when diamagnetic flow and ExB flow cancel each other. The change in the intrinsic momentum transport due to lower hybrid waves is significant when the plasma current is low, which may explain the rotation reversal for low plasma current. The effect of changed q (safety factor) profile by lower hybrid on the intrinsic momentum transport is estimated by gyrokinetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Structural, electronic, and dynamical properties of liquid water by ab initio molecular dynamics based on SCAN functional within the canonical ensemble

We perform ab initio molecular dynamics (AIMD) simulation of liquid water in the canonical ensemble at ambient conditions using the strongly constrained and appropriately normed (SCAN) meta-generalized-gradient approximation (GGA) functional approximation and carry out systematic comparisons with the results obtained from the GGA-level Perdew-Burke-Ernzerhof (PBE) functional and Tkatchenko-Scheffler van der Waals (vdW) dispersion correction inclusive PBE functional. We analyze various properties of liquid water including radial distribution functions, oxygen-oxygen-oxygen triplet angular distribution, tetrahedrality, hydrogen bonds, diffusion coefficients, ring statistics, density of states, band gaps, and dipole moments. We find that the SCAN functional is generally more accurate than the other two functionals for liquid water by not only capturing the intermediate-range vdW interactions but also mitigating the overly strong hydrogen bonds prescribed in PBE simulations. We also compare the results of SCAN-based AIMD simulations in the canonical and isothermal-isobaric ensembles. Our results suggest that SCAN provides a reliable description for most structural, electronic, and dynamical properties in liquid water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗