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 469 records · Page 26

A meshless Galerkin method for non-local diffusion using localized kernel bases

Here, we introduce a meshless method for solving both continuous and discrete variational formulations of a volume constrained, nonlocal diffusion problem. We use the discrete solution to approximate the continuous solution. Our method is nonconforming and uses a localized Lagrange basis that is constructed out of radial basis functions. By verifying that certain inf-sup conditions hold, we demonstrate that both the continuous and discrete problems are well-posed, and also present numerical and theoretical results for the convergence behavior of the method. The stiffness matrix is assembled by a special quadrature routine unique to the localized basis. Combining the quadrature method with the localized basis produces a well-conditioned, symmetric matrix. This then is used to find the discretized solution.

97 MATHEMATICS AND COMPUTING↗

Simulation of the Dynamics of Isothermal Growth of Single-Layer Graphene on a Copper Catalyst in the Process of Chemical Vapor Deposition of Hydrocarbons

A new kinetic model of isothermal growth of single-layer graphene on a copper catalyst as a result of the chemical vapor deposition of hydrocarbons on it at a low pressure has been developed on the basis of in situ measurements of the growth of graphene in the process of its synthesis. This model defines the synthesis of graphene with regard for the chemisorption and catalytic decomposition of ethylene on the surface of a copper catalyst, the diffusion of carbon atoms in the radial direction to the nucleation centers within the thin melted near-surface copper layer, and the nucleation and autocatalytic growth of graphene domains. It is shown that the time dependence of the rate of growth of a graphene domain has a characteristic asymmetrical bell-like shape. The dependences of the surface area and size of a graphene domain and the rate of its growth on the time at different synthesis temperatures and ethylene concentrations have been obtained. Time characteristics of the growth of graphene domains depending on the parameters of their synthesis were calculated. The results obtained can be used for determining optimum regimes of synthesis of graphene in the process of chemical vapor deposition of hydrocarbons on different catalysts with a low solubility of carbon.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Thermophysical properties of FLiBe using moment tensor potentials

Fluoride salts are prospective materials for applications in some next-generation nuclear reactors and their thermophysical properties at various conditions are of interest. Experimental measurement of the properties of these salts is often difficult and, in some cases, unfeasible due to challenges from high temperatures, impurity control, and corrosivity. Therefore, accurate theoretical methods are needed for fluoride salt property prediction. In this work, we used moment tensor potentials (MTP) to approximate the potential energy surface of eutectic FLiBe (66.6% LiF – 33.3% BeF2) predicted by the ab initio (DFT-D3) method. Here, we then used the developed potential and molecular dynamics to obtain several thermophysical properties of FLiBe, including radial distribution functions, density, self-diffusion coefficients, thermal expansion, specific heat capacity, bulk modulus, viscosity, and thermal conductivity. Our results show that the MTP potential approximates the potential energy surface accurately and the overall approach yields very good agreement with experimental values. The converged fitting can be obtained with less than 600 configurations generated from DFT calculations, which data can be generated in just 1200 core hours on today's typical processors. The MTP potential is faster than many machine learning potentials and about one order of magnitude slower than widely used empirical molten salt potentials such as Tosi/Fumi.

36 MATERIALS SCIENCE↗

MULTI-WAVELENGTH VIEW OF KILOPARSEC-SCALE CLUMPS IN STAR-FORMING GALAXIES AT z {approx} 2

This paper studies the properties of kiloparsec-scale clumps in star-forming galaxies at z {approx} 2 through multi-wavelength broadband photometry. A sample of 40 clumps is identified from Hubble Space Telescope (HST)/Advanced Camera for Surveys (ACS) z-band images through auto-detection and visual inspection from 10 galaxies with 1.5 < z < 2.5 in the Hubble Ultra Deep Field, where deep and high-resolution HST/WFC3 and ACS images enable us to resolve structures of z {approx} 2 galaxies down to the kiloparsec scale in the rest-frame UV and optical bands and to detect clumps toward the faint end. The physical properties of clumps are measured through fitting spatially resolved seven-band (BVizYJH) spectral energy distribution to models. On average, the clumps are blue and have similar median rest-frame UV-optical color as the diffuse components of their host galaxies, but the clumps have large scatter in their colors. Although the star formation rate (SFR)-stellar mass relation of galaxies is dominated by the diffuse components, clumps emerge as regions with enhanced specific star formation rates, contributing individually {approx}10% and together {approx}50% of the SFR of the host galaxies. However, the contributions of clumps to the rest-frame UV/optical luminosity and stellar mass are smaller, typically a few percent individually and {approx}20% together. On average, clumps are younger by 0.2 dex and denser by a factor of eight than diffuse components. Clump properties have obvious radial variations in the sense that central clumps are redder, older, more extincted, denser, and less active on forming stars than outskirt clumps. Our results are broadly consistent with a widely held view that clumps are formed through gravitational instability in gas-rich turbulent disks and would eventually migrate toward galactic centers and coalesce into bulges. Roughly 40% of the galaxies in our sample contain a massive clump that could be identified as a proto-bulge, which seems qualitatively consistent with such a bulge-formation scenario.

79 ASTRONOMY AND ASTROPHYSICS↗

Electron kinetics in a positive column of AC discharges in a dynamic regime

Abstract We have performed hybrid kinetic-fluid simulations of a positive column in alternating current (AC) argon discharges over a range of driving frequencies f and gas pressure p for the conditions when the spatial nonlocality of the electron energy distribution function (EEDF) is substantial. Our simulations confirmed that the most efficient conditions of plasma maintenance are observed in the dynamic regime when time modulations of mean electron energy (temperature) are substantial. The minimal values of the root mean square electric field and the electron temperature have been observed at f/p values of about 3 kHz Torr −1 in a tube of radius R = 1 cm. The ionization rate and plasma density reached maximal values under these conditions. The numerical solution of a kinetic equation allowed accounting for the kinetic effects associated with spatial and temporal nonlocality of the EEDF. Using the kinetic energy of electrons as an independent variable, we solved an anisotropic tensor diffusion equation in phase space. We clarified the role of different flux components during electron diffusion in phase space over surfaces of constant total energy. We have shown that the kinetic theory uncovers a more exciting and rich physics than the classical ambipolar diffusion (Schottky) model. Non-monotonic radial distributions of excitation rates, metastable densities, and plasma density have been observed in our simulations at pR > 6 Torr cm. The predicted off-axis plasma density peak in the dynamic regime has never been observed in experiments so far. We hope our results stimulate further experimental studies of the AC positive column. The kinetic analysis could help uncover new physics even for such a well-known plasma object as a positive column in noble gases.

Physics↗

Spatially resolved electron density and electron energy distribution function in Ar magnetron plasmas used for sputter-deposition of ZnO-based thin films

Langmuir probe and trace rare gases optical emission spectroscopy were used to analyze the spatial structure of the electron density and electron energy distribution function (EEDF) in a cylindrical Ar magnetron plasma reactor used for sputter-deposition of ZnO-based thin films. While a typical Bessel (zero order) diffusion profile was observed along the radial direction for the number density of charged particles at 21 cm from the ZnO target, a significant rise of these populations with respect to the Bessel function was seen in the center of the reactor at 4 cm from the magnetron surface. As for the EEDF, it was found to transform from a more or less Maxwellian far from the target to a two-temperature Maxwellian with a depletion of high-energy electrons where magnetic field confinement effects become important. No significant change in the behavior of the electron density and EEDF across a wide range of pressures (5–100 mTorr) and self-bias voltages (115–300 V) was observed during magnetron sputtering of Zn, ZnO, and In{sub 2}O{sub 3} targets. This indicates that sputtering of Zn, In, and O atoms do not play a very significant role on the electron particle balance and electron heating dynamics, at least over the range of experimental conditions investigated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Combined Raman imaging and {sup 18}O tracer analysis for the study of zircaloy-4 high-temperature oxidation in spent fuel pool accident - 2016-0037

Oxidation at high temperature (HT) of cladding materials is expected to be the primary cause of the fuel assemblies degradation in spent fuel storage pool loss of cooling accidents. Unlike a loss-of-coolant accident (LOCA) in a reactor vessel, the presence of air in the atmosphere is expected in a spent fuel pool accident and is known to be an aggravating factor because of the 'catalytic' role of nitrogen on the oxidation. In steam LOCA conditions, a low-temperature (LT) oxidation scale simulating the corrosion oxide formed during in-service use of the assemblies has been observed to have a protective effect regarding subsequent HT oxidation. In the presence of nitrogen, HT oxidation kinetic experiments show that the protective effect of an LT corrosion scale also exists but is lost much earlier than in pure steam. Oxygen transport through LT oxide layers has been studied using the {sup 18}O tracer technique. We performed experiments in {sup 18}O{sub 2} alone as well as in mixed {sup 18}O{sub 2}+N{sub 2} atmospheres, at 850 deg. C. Micro-Raman imaging, both at the specimen surfaces and on cross sections gave clear evidence for different characteristic distributions of {sup 18}O in the scales. Some of these distributions have been correlated with the presence of cracks and porosity in the LT oxide, which allows oxygen to locally penetrate in the scales. In cases in which no radial crack is present, apparent oxygen diffusion coefficients in the oxide were derived from fitting {sup 18}O diffusion profiles. Nitrogen appears to have no or only a limited influence on the oxygen diffusivity, but it is observed to reach the metal-oxide interface faster than oxygen, at which point it reacts with the metal and the sub-stoichiometric oxide to form zirconium nitride.The later conversion of the ZrN into ZrO{sub 2} might lead to the destabilization of the LT scale and to the loss of its protective effect. (authors)

36 MATERIALS SCIENCE↗

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗

Effect of rotor meridional velocity ratio on response to inlet radial and circumferential distortion

Three single transonic fan stages, each having a different meridional velocity ratio across its rotor, were tested with two magnitudes of tip radial distortion and with a 90 deg circumferential distortion imposed on the inlet flow. The rotor with the lowest meridional velocity ratio (less than 0.9 at the tip) demonstrated the least degradation of performance due to these distortions. Loss and deviation angle data (as needed for performance prediction with radial distortion) calculated along actual streamlines for radially distorted flow and correlated against diffusion factor, showed consistent agreement with data calculated along design streamlines for undistorted flow.

Sanger, N. L.↗

On solar cosmic rays

The report summarizes the major results obtained over the last four years in the following areas: (1) transport in interplanetary space which includes solar particle observations combined with theoretical efforts to develop models which simulate the dynamics of the event; (2) transport in the solar corona and release of particles in the interplanetary medium; and (3) characteristics of source spectra and acceleration models. Modelling the propagation of solar cosmic rays through interplanetary space and establishing the magnitude and the radial and energy dependence of the diffusion coefficient K are considered. An important result due to statistical observation is that the coronal propagation must be independent of or weakly dependent on both rigidity and energy. Increasing evidence is found for prolonged injection at the sun.

Van Hollebeke, M. A. I.↗

Stellar convection. III - Convection at large Rayleigh numbers

The results of a numerical, three-dimensional study of convection in a self-gravitating sphere of Boussinesq fluid with a Rayleigh number of 10-billion and a Prandtl number of 1 are presented. The velocity and temperature are computed by using spectral methods in the horizontal and finite-differencing in the radial directions. An eddy viscosity and diffusivity are needed to model the subresolution flow. For Rayleigh numbers much less than 10-billion the flows do not have well-defined inertial ranges, and an eddy viscosity and diffusivity cannot be assigned in this self-consistent manner. By computing the energy spectra as well as the detailed energy budgets as a function of wavenumber, it is shown that for Rs = 10-billion there is an inertial range for the modes corresponding to spherical harmonics with l greater than 6.

Marcus, P. S.↗

Study of spectral exponent on quiet and disturbed days

The spectral exponents of the diurnal and semidiurnal anisotropies are studied on quiet and disturbed days for a period from 1964 to 75. It has been observed for the diurnal anisotropy that the spectral exponent for quiet days is 0.0 from 1964 to 70, however during 1971 and onwards its value is found different from zero and is -0.4; for disturbed days the spectral exponent is observed to be 0.0 from 1964 to 75 with an exception of the year 1966. Further, it has been observed for the semi-diurnal anisotropy that the best value for the spectral exponent on quiet days being + 0.6 from 1964 to 74; however on disturbed days the value of spectral exponent is sometimes negative or positive but different from zero. It was derived from the present investigation that superimposed on the corotational anisotropy, expected from a balance between outward radial convection and inward field aligned diffusion, there is an additional component operating during the period under consideration.

Kumar, S.↗

2D/1D approximations to the 3D neutron transport equation. I: Theory

A new class of '2D/1D' approximations is proposed for the 3D linear Boltzmann equation. These approximate equations preserve the exact transport physics in the radial directions x and y and diffusion physics in the axial direction z. Thus, the 2D/1D equations are more accurate approximations of the 3D Boltzmann equation than the conventional 3D diffusion equation. The 2D/1D equations can be systematically discretized, to yield accurate simulation methods for 3D reactor core problems. The resulting solutions will be more accurate than 3D diffusion solutions, and less expensive to generate than standard 3D transport solutions. In this paper, we (i) show that the simplest 2D/1D equation has certain desirable properties, (ii) systematically discretize this equation, and (iii) derive a stable iteration scheme for solving the discrete system of equations. In a companion paper [1], we give numerical results that confirm the theoretical predictions of accuracy and iterative stability. (authors)

97 MATHEMATICS AND COMPUTING↗

A time-dependent diffusion-convection model for the long-term modulation of cosmic rays

The model considered incorporates, to first order, the direct effects of the time-dependent diffusive propagation of interstellar cosmic rays in a slowly changing interplanetary medium. A one-dimensional model is discussed along with a three-dimensional radial wind. An evaluation is conducted of average diffusive propagation times. It is found that the effects of time-dependent diffusive propagation can be quite significant. A rigidity-dependent time delay or 'lag' in the modulated spectra is predicted.

Ogallagher, J. J.↗

Modulation and diffusion theory of cosmic rays

Papers given on the modulation of galactic cosmic rays by the solar wind at the 14th International Cosmic Ray Conference are reviewed. Some of the topics treated in this review are Pioneer and Helios radial gradient measurement, heliolatitude effects in modulation models, diffusion (scattering) theory - including pitch angle diffusion diagrams - solar flare particle transport theory and cosmic ray fluctuations at neutron monitor energies.

Forman, M. A.↗

The radial variation of solar flare proton anisotropies observed in deep space on Pioneers 10 and 11

The anisotropy of solar-flare protons (greater than 3 MeV) was evaluated at the time of maximum for each of 24 events measured by Pioneers 10 and 11 and was observed to decrease gradually between 1 and 4.7 AU. Since diffusive streaming is related directly to the effective diffusion coefficient parallel to the interplanetary spiral magnetic field, these observations provide a measure of the radial dependence of interplanetary diffusion parameters. Analysis of these results suggests that the spiral field increases gradually with heliocentric distance. An expression is given for the radial dependence of that field.

Mccarthy, J.↗

The stability of the contact interface of cylindrical and spherical shock tubes

The stability of the contact interface for radial shock tubes is investigated as a model for explosive dispersal. The advection upstream splitting method with velocity and pressure diffusion (AUSM+-up) is used to solve for the radial base flow. To explore the stability of the resulting contact interface, perturbed governing equations are derived assuming harmonic modes in the transverse directions. The perturbed harmonic flow is solved by assuming an initial disturbance and using a perturbed version of AUSM+-up derived in this paper. The intensity of the perturbation near the contact interface is computed and compared to theoretical results obtained by others. Despite the simplifying assumptions of the theoretical analysis, very good agreement is confirmed. Not only can the magnitude of the instability be predicted during the initial expansion, but also remarkably the agreement between the numerical and theoretical results can be maintained through the collision between the secondary shock and the contact interface. Since the theoretical results only depend upon the time evolution of the base flow, the stability of various modes could be quickly studied without explicitly solving a system of partial differential equations for the perturbed flow.

42 ENGINEERING↗