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 901 records · Page 50

Variable stator radial turbine

A radial turbine stage with a variable area nozzle was investigated. A high work capacity turbine design with a known high performance base was modified to accept a fixed vane stagger angle moveable sidewall nozzle. The nozzle area was varied by moving the forward and rearward sidewalls. Diffusing and accelerating rotor inlet ramps were evaluated in combinations with hub and shroud rotor exit rings. Performance of contoured sidewalls and the location of the sidewall split line with respect to the rotor inlet was compared to the baseline. Performance and rotor exit survey data are presented for 31 different geometries. Detail survey data at the nozzle exit are given in contour plot format for five configurations. A data base is provided for a variable geometry concept that is a viable alternative to the more common pivoted vane variable geometry radial turbine.

Rogo, C.↗

Modulation of Galactic Cosmic Rays in the Inner Heliosphere, Comparing with PAMELA Measurements

We develop a numerical model to study the time-dependent modulation of galactic cosmic rays in the inner heliosphere. In the model, a time-delayed modified Parker heliospheric magnetic field (HMF) and a new diffusion coefficient model, NLGCE-F, from Qin and Zhang, are adopted. In addition, the latitudinal dependence of magnetic turbulence magnitude is assumed to be ∼(1+sin{sup 2}θ)/2 from the observations of Ulysses, and the radial dependence is assumed to be ∼r{sup S}, where we choose an expression of S as a function of the heliospheric current sheet tilt angle. We show that the analytical expression used to describe the spatial variation of HMF turbulence magnitude agrees well with the Ulysses, Voyager 1, and Voyager 2 observations. By numerically calculating the modulation code, we get the proton energy spectra as a function of time during the recent solar minimum, it is shown that the modulation results are consistent with the Payload for Antimatter-Matter Exploration and Light-nuclei Astrophysics measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

Computational Study of Molecular Structure and Self-Association of Tri-n-butyl Phosphates in n-Dodecane

Tri-n-butyl phosphate is an important extractant used in solvent extraction process for the recovery of uranium and plutonium from spent nuclear fuel. To understand the fundamental molecular level behavior of extracting agents in solution, an atomistic parameterization study was carried out using the AMBER force field to model TBP molecule and n-dodecane molecule, a commonly used organic solvent, for molecular dynamics simulations. For validation of the optimized force field, various thermophysical properties of pure TBP and pure n-dodecane in the bulk liquid phase such as mass density, dipole moment, self-diffusion coefficient and heat of vaporization were calculated and compared favorably with experimental values. The molecular structure of TBPs in n-dodecane at various TBP concentrations was examined based on radial distribution functions and 2D potential mean force, which was used as criteria for identifying TBP aggregates. The dimerization constant of TBP in n-dodecane was also obtained and matches the experimental value. The U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences funded the work performed by LXD.

molecular dynamics simulations↗

Counter-ion binding and mobility in the presence of hydrophobic polyions – combining molecular dynamics simulations and NMR

Counter-ion binding and mobility in aqueous solutions of partially hydrophobic ionene oligoions is studied here by a combination of all-atomic molecular dynamics (MD) simulations and NMR ({sup 19}F and {sup 81}Br nuclei) measurements. We present results for 12, 12–ionenes in the presence of different halide ions (F{sup −}, Cl{sup −}, Br{sup −} and I{sup −}), as well as their mixtures; the latter allowing us to probe counter-ion selectivity of these oligoions. We consolidate both structural and dynamic information, in particular simulated radial distribution functions and average residence times of counter-ions in the vicinity of ionenes and NMR data in the form of counter-ion chemical shift and self-diffusion coefficients. On one hand, previously reported enthalpy of dilution and mixing measurements show a reverse counter-ion sequence for 12, 12–ionenes with respect to their less hydrophobic 3, 3– and 6, 6– analogues. On the other hand, the current MD and NMR data, reflecting the counter-ion binding tendencies to the ionene chain, give evidence for the same ordering as that observed by MD for 3, 3–ionenes. This is not seen as a contradiction and can be rationalized on the basis of increasing chain hydrophobicity, which has different consequences for enthalpy and ion-binding. The latter is reflecting free energy changes and as such includes both enthalpic and entropic contributions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

MULTI-ZONE MODELING OF THE PULSAR WIND NEBULA HESS J1825-137

The pulsar wind nebula associated with PSR J1826-1334, HESS J1825-137, is a bright very high energy (VHE) source with an angular extent of {approx}1 Degree-Sign and spatially resolved spectroscopic TeV measurements. The gamma-ray spectral index is observed to soften with increasing distance from the pulsar, likely the result of cooling losses as electrons traverse the nebula. We describe analysis of X-ray data of the extended nebula, as well as three-dimensional time-dependent spectral energy distribution modeling, with emphasis on the spatial variations within HESS J1825-137. The multi-wavelength data place significant constraints on electron injection, transport, and cooling within the nebula. The large size and high nebular energy budget imply a relatively rapid initial pulsar spin period of 13 {+-} 7 ms and an age of 40 {+-} 9 kyr. The relative fluxes of each VHE zone can be explained by advective particle transport with a radially decreasing velocity profile with v(r){proportional_to}r{sup -0.5}. The evolution of the cooling break requires an evolving magnetic field which also decreases radially from the pulsar, B(r,t){proportional_to}r{sup -0.7} E-dot (t){sup 1/2}. Detection of 10 TeV flux {approx}80 pc from the pulsar requires rapid diffusion of high-energy particles, contrary to the common assumption of toroidal magnetic fields with strong magnetic confinement. The model predicts a rather uniform Fermi Large Area Telescope (LAT) surface brightness out to {approx}1 Degree-Sign from the pulsar, in good agreement with the recently discovered LAT source centered 0.{sup 0}5 southwest of PSR J1826-1334 with extension 0.{sup 0}6 {+-} 0.{sup 0}1.

79 ASTRONOMY AND ASTROPHYSICS↗

Performance of a small annular turbojet combustor designed for low cost

Performance investigations were conducted on a combustor utilizing several cost-reducing innovations and designed for use in a low-cost 4448-N thrust turbojet engine for commercial light aircraft. Low-cost features included simple, air-atomizing fuel injectors; combustor liners of perforated sheet; and the use of inexpensive type 304 stainless-steel material. Combustion efficiencies at the cruise and sea-level-takeoff design points were approximately 97 and 98 percent, respectively. The combustor isothermal pressure loss was 6.3 percent at the cruise-condition diffuser inlet Mach number of 0.34. The combustor exit temperature pattern factor was less than 0.24 at both the cruise and sea-level-takeoff design points. The combustor exit average radial temperature profiles at all conditions were in very good agreement with the design profile.

Fear, J. S.↗

Scale-model studies for the improvement of flow patterns of a low-speed tunnel

Parameters affecting the performance of a scale model V/STOL wind tunnel were measured. These included: (1) static pressure variation around the closed tunnel circuit as affected by the application of various screens; (2) pressure rise and radial distribution of flow through the fan; (3) variation of screen parameters and screen location affecting the flow; and (4) effects of multiple screens on the velocity distribution in the fourth diffuser. In these tests the fan was driven through a long shaft by an externally situated 11,190-W (15-HP) motor equipped to turn at 4 speeds.

Barna, P. S.↗

Calculation of Last Closed Drift Shells for the 2013 GEM Radiation Belt Challenge Events

Radiation belt behavior is often analyzed in terms of adiabatic invariants, of which the third roughly characterizes the radial distance of particle drift shells. The outermost, or last closed drift shell, can be an important boundary for numerical modeling, especially for drift-averaged treatments such as three-dimensional diffusion codes. In this work, we discuss calculation of the last closed drift shell, using the widely used International Radiation Belt Environment Modeling (IRBEM) Library, the LanlGeoMag code, and a guiding center code named AFRL-Shell, in conjunction with a variety of current magnetic field models. Furthermore, we present results for the four events in 2013, comprising the radiation belt challenge organized by the Quantitative Assessment of Radiation Belt Modeling focus group of the Geospace Environment Modeling (GEM) program.

79 ASTRONOMY AND ASTROPHYSICS↗

Construction of generalized quasilinear diffusion coefficient using neural networks with physical restrictions

The quasilinear diffusion coefficient (D QL ) derived from our machine learning framework shows comparable trends with the ground truth D QL obtained from GENRAY-CQL3D simulations. Additionally, for the strong absorption cases, the radial current drive profiles generated using the D QL from our model exhibit consistent behavior with those obtained from the original simulation. These findings indicate the potential of our surrogate modeling approach with physical restrictions to replicate key wave–plasma interaction characteristics while reducing computational costs. Traditionally, calculating D QL for wave–particle interactions relies on computationally intensive wave simulations coupled with Fokker–Planck solvers. To address this challenge, we developed a machine learning-based surrogate model with physical restrictions derived from cold plasma theory and bounce-averaged damping effects. First, we establish the propagation domain of Lower Hybrid Waves in the (N∥, ρ) space by identifying the accessibility limit and determining the upper and lower bounds of N∥ using the Potential Power Deposition (PPD) method. Subsequently, leveraging a database constructed using Latin hypercube sampling alongside the underlying physical restrictions (e.g. PPD), machine learning methods including U-Net and Recurrent Neural Networks are employed to design a physics-restricted machine learning framework capable of reconstructing D QL .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Preliminary Compendium of the Initial Conditions and Evolution of the Fundamental Interactions of SiC with Ag and Pd

Tri-structural isotropic (TRISO) particles are the fuel of choice for very-high temperature reactors. At the core, these spherical particles consist of an inner fuel kernel of UO{sub 2}, UC{sub 2} or a combination of both, which is coated radially outward by successive layers of low density pyrolytic graphite, an inner high density pyrolytic graphite layer, silicon carbide (SiC) and an outer high density pyrolytic graphite layer. SiC is the main diffusion barrier against the release of fission products due to its hardness and high melting point. Irradiation testing of TRISO particles indicates that trace amounts of metallic fission products can damage or potentially escape through the barrier layer. Surrogate experimental studies have previously been performed by various methods seeking to shed light on and explain the diffusion behavior of Ag and Pd. Although some of these experiments have yielded realistic diffusion coefficients the discrepancy of results between different groups is such that the essential details of the diffusive behavior are still to be fully understood. This work aims to understand the fundamental physical and chemical interactions that take place at the interface of SiC and Ag/Pd, and which give rise to the observed diffusion behavior of these metallic species in TRISO particles. Theoretical simulations of the diffusion of Ag through SiC have been exhaustive but typically yield bulk diffusion coefficients that due to the nature of bulk diffusion phenomena (i.e. slow diffusion) the diffusion coefficients are calculated to be orders of magnitude lower than what experimental observations suggest. This in part is due to the difficulty in simulating realistic models that include various polymorphs of SiC, and the dynamical structural transformations that SiC undergoes in the TRISO environment (which are also to be explained). The result is a particularly challenging problem to solve given that the initial conditions of the interactions of Ag and Pd with SiC are difficult to address, both in theoretical models and surrogate experiments. In order to improve these models the approach taken here is to consider fundamental interactions between Ag and SiC that take place in short time scales by varying the initial environmental conditions, i.e. annealing temperature and time, initial amounts of the metallic species (M) Ag, Pd in SiC/M/SiC multilayers. This work aims to understand the fundamental physical and chemical interactions that take place at the interface of SiC and Ag/Pd, and which give rise to the observed diffusion behavior of these metallic species in TRISO particles. In order to do so we have synthesized SiC/Ag/SiC and SiC/Pd/SiC multilayers of various thicknesses by PLD and characterized them using SEM/EDS. Initial results show that upon annealing Ag agglomerates and ruptures the SiC layer to be released. In the case of Pd a strong chemical interaction is deduced from the slow interlayer diffusion and limited release. Further studies using photoelectron spectroscopy and RHEED will be carried out in order obtain detailed information on the chemical and physical transitions that undergo upon annealing of these systems. (authors)

36 MATERIALS SCIENCE↗

Prediction of divertor heat flux width for ITER using BOUT++ transport and turbulence module

Investigation of turbulent transport dynamics in scrape-off-layer (SOL) and divertor heat flux width prediction is performed for ITER. Both BOUT++ transport and BOUT++ turbulence codes are applied to capture the physics on different temporal scales. Simulations start with an ITER 15MA baseline scenario profile generated by CORSICA (Kim et al 2015 Paper ITER_D_R9T8J9 v1.1). In BOUT++ transport code, the plasma parameters ( n i , T i , T e ) and radial electric () profiles are evolved to steady state. Here, the initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients are calculated by inverting the plasma profiles inside the separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reducing DLOFC fuel temperatures by mixing thorium with LEU in a single-zone six-pass fuel cycle in a PBMR-DPP-400 core - 18548

Many studies have been done in neutronics and thermal-hydraulic simulation of the standard 6- pass fuel recirculation scheme for the standard 9.6 wt % enriched, 9 g per fuel sphere low-enriched uranium (LEU) fuel in the PBMR-DPP-400, using different versions of the VSOP diffusion codes. Maximum DLOFC temperatures were all below the upper limit of 1600 deg. C. The DLOFC temperature is highly dependent on the peaking factor of the power; meaning a lower maximum DLOFC temperature can be obtained by suppressing the axial power peak and by moving the radial power peak towards the external reflector. In this study the standard 6-pass fuel recirculation was retained. The improvement strategy was thus attempted by means of flattening the axial power profile by mixing substantial amounts of thorium into the LEU fuel. The addition of thorium led to breading of substantial amounts of {sup 233}U. This led to slower depletion of enrichment with burn-up, which increased fuel reactivity and power densities near the bottom of the core and thus flattened the axial power profile. The effect was a reduction in maximum DLOFC temperature by 44 deg. C. The simulations were made using the VSOP- 99/05 diffusion code. It was further shown that the results obtained are also applicable to the Chinese HTR-PM and the proposed strategies for further improvement can be expected to produce even much better results in Block Reactors than in Pebble Bed Reactors. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Role of the coupling between neck and radial degrees of freedom in evolution from dinucleus to mononucleus

The coupled Langevin equations in two-dimensional collective space are used to study the dynamics of nuclear neck growth. Special attention is paid to the effects of coupling between neck and radial degrees of freedom on the evolution from dinucleus to mononuleus. The dynamic model is applied for the study of neck evolution of the mass asymmetric system {sup 50}Ti + {sup 249}Cf. In order to estimate the effects of the coupling, we use the equations in the coupled and uncoupled cases. Our results show that the coupling between neck and radial motions reduces the neck growth velocity and delays the transition from dinuclear to mononuclear regimes. In addition, by solving these dynamic equations we get the probability distributions of radial degree of freedom at the injection point s{sub in} in the asymmetric fission valley. In this way, s{sub in} is no longer an adjustable parameter in the fusion-by-diffusion model. The distributions obtained are located at positions very close to s=0. Moreover, the coupling significantly reduces the fluctuation in the s{sub in} space. Our results show that the quasifissionlike events are unlikely to take place during the transition period for the mass asymmetric systems. Based on the s{sub in} distribution obtained, the evaporation residue (ER) cross sections for 3n and 4n evaporation channels in the {sup 50}Ti + {sup 249}Cf reaction leading to formation of {sup 296}120 and {sup 295}120 isotopes are calculated. The maximum ER cross sections in 3n and 4n channels with the s{sub in} distributions are equal to 0.1 and 0.065 pb, respectively, which are more than two times larger than those obtained by the fusion-by-diffusion model with the s{sub in}=0 assumption.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Anisotropic Azimuthal Power and Temperature distribution on FuelRod: Impact on Hydride Distribution

The degradation of the zirconium cladding may limit nuclear fuel performance. In the high temperature environment of a reactor, the zirconium in the cladding corrodes, releasing hydrogen in the process. Some of this hydrogen is absorbed by the cladding in a highly inhomogeneous manner. The distribution of the absorbed hydrogen is extremely sensitive to temperature and stress concentration gradients. The absorbed hydrogen tends to concentrate near lower temperatures. This hydrogen absorption and hydride formation can cause cladding failure. This project set out to improve the hydrogen distribution prediction capabilities of the BISON fuel performance code. The project was split into two primary sections, first was the use of a high fidelity multi-physics coupling to accurately predict temperature gradients as a function of r, θ , and z, and the second was to use experimental data to create an analytical hydrogen precipitation model. The Penn State version of thermal hydraulics code COBRA-TF (CTF) was successfully coupled to the DeCART neutronics code. This coupled system was verified by testing and validated by comparison to FRAPCON data. The hydrogen diffusion and precipitation experiments successfully calculated the heat of transport and precipitation rate constant values to be used within the hydrogen model in BISON. These values can only be determined experimentally. These values were successfully implemented in precipitation, diffusion and dissolution kernels that were implemented in the BISON code. The coupled output was fed into BISON models and the hydrogen and hydride distributions behaved as expected. Simulations were conducted in the radial, axial and azimuthal directions to showcase the full capabilities of the hydrogen model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Faint Object Camera imaging and spectroscopy of NGC 4151

We describe ultraviolet and optical imaging and spectroscopy within the central few arcseconds of the Seyfert galaxy NGC 4151, obtained with the Faint Object Camera on the Hubble Space Telescope. A narrowband image including (O III) lambda(5007) shows a bright nucleus centered on a complex biconical structure having apparent opening angle approximately 65 deg and axis at a position angle along 65 deg-245 deg; images in bands including Lyman-alpha and C IV lambda(1550) and in the optical continuum near 5500 A, show only the bright nucleus. In an off-nuclear optical long-slit spectrum we find a high and a low radial velocity component within the narrow emission lines. We identify the low-velocity component with the bright, extended, knotty structure within the cones, and the high-velocity component with more confined diffuse emission. Also present are strong continuum emission and broad Balmer emission line components, which we attribute to the extended point spread function arising from the intense nuclear emission. Adopting the geometry pointed out by Pedlar et al. (1993) to explain the observed misalignment of the radio jets and the main optical structure we model an ionizing radiation bicone, originating within a galactic disk, with apex at the active nucleus and axis centered on the extended radio jets. We confirm that through density bounding the gross spatial structure of the emission line region can be reproduced with a wide opening angle that includes the line of sight, consistent with the presence of a simple opaque torus allowing direct view of the nucleus. In particular, our modelling reproduces the observed decrease in position angle with distance from the nucleus, progressing initially from the direction of the extended radio jet, through our optical structure, and on to the extended narrow-line region. We explore the kinematics of the narrow-line low- and high-velocity components on the basis of our spectroscopy and adopted model structure.

Boksenberg, A.↗

Interchange Slip-Running Reconnection and Sweeping SEP-Beams

We present a new model to explain how particles, accelerated at a reconnection site that is not magnetically connected to the Earth, could eventually propagate along the well-connected open flux tube. Our model is based on the results of a low-beta resistive magnetohydrodynamics simulation of a three-dimensional line-tied and initially current-free bipole, that is embedded in a non-uniform open potential field. The topology of this configuration is that of an asymmetric coronal null-point, with a closed fan surface and an open outer spine. When driven by slow photospheric shearing motions, field lines, initially fully anchored below the fan dome, reconnect at the null point, and jump to the open magnetic domain. This is the standard interchange mode as sketched and calculated in 2D. The key result in 3D is that, reconnected open field lines located in the vicinity of the outer spine, keep reconnecting continuously, across an open quasi-separatrix layer, as previously identified for non-open-null-point reconnection. The apparent slipping motion of these field lines leads to form an extended narrow magnetic flux tube at high altitude. Because of the slip-running reconnection, we conjecture that if energetic particles would be travelling through, or be accelerated inside, the diffusion region, they would be successively injected along continuously reconnecting field lines that are connected farther and farther from the spine. At the scale of the full Sun, owing to the super-radial expansion of field lines below 3 solar radius, such energetic particles could easily be injected in field lines slipping over significant distances, and could eventually reach the distant flux tube that is well-connected to the Earth.

Masson, S.↗

Push-pull tracer tests: Their information content and use for characterizing non-Fickian, mobile-immobile behavior: INFORMATION CONTENT OF PUSH-PULL TESTS

Path reversibility and radial symmetry are often assumed in push-pull tracer test analysis. In reality, heterogeneous flow fields mean that both assumptions are idealizations. In this paper, to understand their impact, we perform a parametric study which quantifies the scattering effects of ambient flow, local-scale dispersion, and velocity field heterogeneity on push-pull breakthrough curves and compares them to the effects of mobile-immobile mass transfer (MIMT) processes including sorption and diffusion into secondary porosity. We identify specific circumstances in which MIMT overwhelmingly determines the breakthrough curve, which may then be considered uninformative about drift and local-scale dispersion. Assuming path reversibility, we develop a continuous-time-random-walk-based interpretation framework which is flow-field-agnostic and well suited to quantifying MIMT. Adopting this perspective, we show that the radial flow assumption is often harmless: to the extent that solute paths are reversible, the breakthrough curve is uninformative about velocity field heterogeneity. Our interpretation method determines a mapping function (i.e., subordinator) from travel time in the absence of MIMT to travel time in its presence. A mathematical theory allowing this function to be directly “plugged into” an existing Laplace-domain transport model to incorporate MIMT is presented and demonstrated. Algorithms implementing the calibration are presented and applied to interpretation of data from a push-pull test performed in a heterogeneous environment. A successful four-parameter fit is obtained, of comparable fidelity to one obtained using a million-node 3-D numerical model. In conclusion, we demonstrate analytically and numerically how push-pull tests quantifying MIMT are sensitive to remobilization, but not immobilization, kinetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multilevel well modeling in aggregation-based nonlinear multigrid for multiphase flow in porous media

A full approximation scheme (FAS) nonlinear multigrid solver for two-phase flow and transport problems driven by wells with multiple perforations is developed here. It is an extension to our previous work on FAS solvers for diffusion and transport problems. The solver is applicable to discrete problems defined on unstructured grids as the coarsening algorithm is aggregation-based and algebraic. To construct coarse basis that can better capture the radial flow near wells, coarse grids in which perforated well cells are not near the coarse-element interface are desired. This is achieved by an aggregation algorithm proposed in this paper that makes use of the location of well cells in the cell-connectivity graph. Numerical examples in which the FAS solver is compared against Newton's method on benchmark problems are given. In particular, for a refined version of the SAIGUP model, the FAS solver is at least 35% faster than Newton's method for time steps with a CFL number greater than 10.

58 GEOSCIENCES↗