Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion equations”

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 1,387 records · Page 77

Dynamic Theory of Relativistic Electrons Stochastic Heating by Whistler Mode Waves with Application to the Earth Magnetosphere

In the Hamiltonian approach an electron motion in a coherent packet of the whistler mode waves propagating along the direction of an ambient magnetic field is studied. The physical processes by which these particles are accelerated to high energy are established. Equations governing a particle motion were transformed in to a closed pair of nonlinear difference equations. The solutions of these equations have shown there exists the energetic threshold below that the electron motion is regular, and when the initial energy is above the threshold an electron moves stochastically. Particle energy spectra and pitch angle electron scattering are described by the Fokker-Planck-Kolmogorov equations. Calculating the stochastic diffusion of electrons due to a spectrum of whistler modes is presented. The parametric dependence of the diffusion coefficients on the plasma particle density, magnitude of wave field, and the strength of magnetic field is studies. It is shown that significant pitch angle diffusion occurs for the Earth radiation belt electrons with energies from a few keV up to a few MeV.

Khazanov, G. V.↗

Snow Micro-Structure Model

PIKA is a MOOSE-based application for modeling micro-structure evolution of seasonal snow. The model will be useful for environmental, atmospheric, and climate scientists. Possible applications include application to energy balance models, ice sheet modeling, and avalanche forecasting. The model implements physics from published, peer-reviewed articles. The main purpose is to foster university and laboratory collaboration to build a larger multi-scale snow model using MOOSE. The main feature of the code is that it is implemented using the MOOSE framework, thus making features such as multiphysics coupling, adaptive mesh refinement, and parallel scalability native to the application. PIKA implements three equations: the phase-field equation for tracking the evolution of the ice-air interface within seasonal snow at the grain-scale; the heat equation for computing the temperature of both the ice and air within the snow; and the mass transport equation for monitoring the diffusion of water vapor in the pore space of the snow.

Micah Johnson, Andrew Slaughter↗

Evaluation of rodded BWR assembly pin powers with simulate

In the development of Studsvik's nodal code SIMULATE5, special attention has been given to the accurate description of pin powers. The code solves the multi-group diffusion or simplified P3 equations with high spatial resolution. A short description is given of the handling of material heterogeneities in the axial and radial directions - the axial re-homogenization and the radial sub-mesh model, respectively. Furthermore, two models which are of special importance in the presence of control rods, and hence on the pellet clad interaction (PCI) phenomenon, are described; the quarter-assembly thermal-hydraulic treatment and the depletion of absorber material. A numerical example for a real core shows that while the impact of the fine grained description is negligible on global parameters such as k{sub eff}, the effect on pin powers may be substantial. For a deeply inserted and highly depleted control rod, the difference in pin powers caused by neglecting the advanced models of SIMULATE5, may be as high as 15 %. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

High-fidelity simulation-driven model development for coarse-grained computational fluid dynamics

Nuclear reactor safety analysis requires analysis of a broad range of accident scenarios and determining their consequences. For a nuclear power plant behavior, it is impossible to obtain experimental data at sufficiently large scales to support calibration and validation of various system models. In single-phase flow convective problems, high-resolution methods of Computational Fluid Dynamics (CFD) such as Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) can provide high fidelity results under conditions where physical measurements are unavailable. However, such high-resolution simulations are computationally expensive and not suitable for simulation of long transient scenarios in nuclear reactor accidents. In this work, we investigate the use of a high fidelity simulation-driven approach to model sub-grid scale (SGS) effect in Coarse Grained Computational Fluid Dynamics CG-CFD. Notably, fine-mesh simulations are used to construct physics-informed statistical surrogate of coarse-mesh SGS model. For an initial analysis, we consider a case of turbulent natural convection in a volumetrically heated fluid layer with a thermally insulated lower boundary and isothermal upper boundary. This scenario of unstable stratification is relevant to turbulent natural convection in a severe nuclear reactor accident, as well as in containment mixing and passive cooling. For this type of flow, the SGS effect is modeled by an added turbulent diffusivity in the energy equation. It is shown that a global correction for the energy equation is sufficient to achieve a significant improvement to the CG-CFD prediction of thermal mixing in the fluid layer. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Model of the boundary layer of a vacuum-arc magnetic filter

A model is developed to describe the electrostatic boundary layer in a positively biased magnetic filter in filtered arcs with low collisionality. The set of equations used includes the electron momentum equation, with an anomalous collision term due to micro-instabilities leading to Bohm diffusion, electron mass conservation, and Poisson equation. Analytical solutions are obtained, valid for the regimes of interest, leading to an explicit expression to determine the electron density current to the filter wall as a function of the potential of the filter and the ratio of electron density at the plasma to that at the filter wall. Using a set of planar and cylindrical probes it is verified experimentally that the mentioned ratio of electron densities remains reasonably constant for different magnetic field values and probe bias, which allows to obtain a closed expression for the current. Comparisons are made with the experimentally determined current collected at different sections of a positively biased straight filter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Calculation of two-dimension radial electric field in boundary plasmas by using BOUT++

The steady state radial electric field (Er) is measured by coupling a plasma transport model with the quasi-neutrality constraint and the vorticity equation within the BOUT++ framework. Due to the experimentally measured plasma density and temperature profiles in Alcator C-Mod discharges, the effective radial particle and heat diffusivities are inferred from the set of plasma transport equations. The effective diffusivities are then extended into the scrape-off layer (SOL) to calculate the plasma density, temperature and flow profiles across the separatrix into the SOL with the electrostatic sheath boundary conditions (SBC) applied on the divertor plates. Given these diffusivities, the electric field can be calculated self-consistently across the separatrix from the vorticity equation with SBC coupled to the plasma transport equations. The sheath boundary conditions act to generate a large and positive Er in the SOL, which is consistent with experimental measurements. The effect of magnetic particle drifts is shown to play a significant role on local particle transport and Er by inducing a net particle flow in both the edge and SOL regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Model calculations for diffuse molecular clouds

A steady state isobaric cloud model is developed. The pressure, thermal, electrical, and chemical balance equations are solved simultaneously with a simple one dimensional approximation to the equation of radiative transfer appropriate to diffuse clouds. Cooling is mainly by CII fine structure transitions, and a variety of heating mechanisms are considered. Particular attention is given to the abundance variation of H2. Inhomogeneous density distributions are obtained because of the attenuation of the interstellar UV field and the conversion from atomic to molecular hyrodgen. The effects of changing the model parameters are described and the applicability of the model to OAO-3 observations is discussed. Good qualitative agreement with the fractional H2 abundance determinations has been obtained. The observed kinetic temperatures near 80 K can also be achieved by grain photoelectron heating. The problem of the electron density is solved taking special account of the various hydrogen ions as well as heavier ones.

Glassgold, A. E.↗

Preshock region acceleration of implanted cometary H(+) and O(+)

A self-consistent, three-fluid model of plasma transport and implanted ion acceleration in the unshocked solar wind is presented. The solar wind plasma is depleted by charge exchange with the expanding cometary exosphere, while implanted protons and heavy ions are produced by photoionization and charge transfer and lost by charge exchange. A generalized transport equation describing convection, adiabatic and diffusive velocity change, and the appropriate production terms is used to describe the evolution of the two cometary ion components, while the moments of the Boltzmann equation are used to calculate the solar wind density and pressure. The flow velocity is obtained self-consistently by combining the conservation equations of the three ion species. The results imply that second-order Fermi acceleration can explain the implanted spectra observed in the unshocked solar wind. Comparison of measured and calculated distribution indicates that spatial diffusion of implanted ions probably plays an important role in forming the energetic particle environment in the shock vicinity.

Gombosi, Tamas I.↗

Cosmic-ray propagation with DRAGON2: I. numerical solver and astrophysical ingredients

We present version 2 of the DRAGON code designed for computing realistic predictions of the CR densities in the Galaxy. The code numerically solves the interstellar CR transport equation (including inhomogeneous and anisotropic diffusion, either in space and momentum, advective transport and energy losses), under realistic conditions. The new version includes an updated numerical solver and several models for the astrophysical ingredients involved in the transport equation. Improvements in the accuracy of the numerical solution are proved against analytical solutions and in reference diffusion scenarios. The novel features implemented in the code allow to simulate the diverse scenarios proposed to reproduce the most recent measurements of local and diffuse CR fluxes, going beyond the limitations of the homogeneous galactic transport paradigm. To this end, several applications using DRAGON2 are presented as well. This new version facilitates the users to include their own physical models by means of a modular C++ structure.

79 ASTRONOMY AND ASTROPHYSICS↗

Surface transport processes in charged porous media

Surface transport processes are important in chemistry, colloidal sciences, engineering, biology, and geophysics. Natural or externally produced charges on surfaces create electrical double layers (EDLs) at the solid-liquid interface. The existence of the EDLs produces several complex processes including bulk and surface transport of ions. In this work, a model is presented to simulate bulk and transport processes in homogeneous porous media comprising big pores. It is based on a theory for capacitive charging by ideally polarizable porous electrodes without Faradaic reactions or specific adsorption of ions. A volume averaging technique is used to derive the averaged transport equations in the limit of thin electrical double layers. Description of the EDL between the electrolyte solution and the charged wall is accomplished using the Gouy-Chapman-Stern (GCS) model. The surface transport terms enter into the average equations due to the use of boundary conditions for diffuse interfaces. Two extra surface transports terms appear in the closed average equations. One is a surface diffusion term equivalent to the transport process in non-charged porous media. The second surface transport term is a migration term unique to charged porous media. The effective bulk and transport parameters for isotropic porous media are calculated solving the corresponding closure problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling scalar flux and the energy and dissipation equations

Closure models derived from the Two-Scale Direct-Interaction Approximation were compared with data from direct simulations of turbulence. Attention was restricted to anisotropic scalar diffusion models, models for the energy dissipation equation, and models for energy diffusion.

Yoshizawa, A.↗

Diffusion and Clustering of Carbon Dioxide on Non-porous Amorphous Solid Water

Observations by ISO and Spitzer toward young stellar objects showed that CO{sub 2} segregates in the icy mantles covering dust grains. Thermal processing of the ice mixture was proposed as being responsible for the segregation. Although several laboratories studied thermally induced segregation, a satisfying quantification is still missing. We propose that the diffusion of CO{sub 2} along pores inside water ice is the key to quantify segregation. We combined Temperature Programmed Desorption and Reflection Absorption InfraRed Spectroscopy to study how CO{sub 2} molecules interact on a non-porous amorphous solid water (np-ASW) surface. We found that CO{sub 2} diffuses significantly on an np-ASW surface above 65 K and clusters are formed at well below one monolayer. A simple rate equation simulation finds that the diffusion energy barrier of CO{sub 2} on np-ASW is 2150 ± 50 K, assuming a diffusion pre-exponential factor of 10{sup 12} s{sup −1}. This energy should also apply to the diffusion of CO{sub 2} on the wall of pores. The binding energy of CO{sub 2} from CO{sub 2} clusters and CO{sub 2} from H{sub 2}O ice has been found to be 2415 ± 20 K and 2250 ± 20 K, respectively, assuming the same prefactor for desorption. CO{sub 2}–CO{sub 2} interaction is stronger than CO{sub 2}–H{sub 2}O interaction, in agreement with the experimental finding that CO{sub 2} does not wet the np-ASW surface. For comparison, we carried out similar experiments with CO on np-ASW, and found that the CO–CO interaction is always weaker than CO–H{sub 2}O. As a result, CO wets the np-ASW surface. This study should be of help to uncover the thermal history of CO{sub 2} on the icy mantles of dust grains.

79 ASTRONOMY AND ASTROPHYSICS↗

Diffusion, Fluxes, Friction Forces, and Joule Heating in Two-Temperature Multicomponent Magnetohydrodynamics

The relationship between Joule heating, diffusion fluxes, and friction forces has been studied for both total and electron thermal energy equations, using general expressions for multicomponent diffusion in two-temperature plasmas with the velocity dependent Lorentz force acting on charged species in a magnetic field. It is shown that the derivation of Joule heating terms requires both diffusion fluxes and friction between species which represents the resistance experienced by the species moving at different relative velocities. It is also shown that the familiar Joule heating term in the electron thermal energy equation includes artificial effects produced by switching the convective velocity from the species velocity to the mass-weighted velocity, and thus should not be ignored even when there is no net energy dissipation.

Chang, C. H.↗

A computational method for subsonic compressible flow in diffusers

Integral methods are used to model subsonic flow in planar diffusers operating with or without stall. Turbulent boundary layers are computed with a lag-entrainment method employing new correlations for shape factor and skin friction, which are based on a Coles-Van Driest velocity profile. The inviscid flow is modeled as one dimensional and is easily coupled with the boundary layer method to yield equations that can be solved by integrating in the downstream direction. Analysis of the equations describing flow in a diffuser shows that the boundary layer separation singularity may not be eliminated by viscous-inviscid interaction; a modified interaction scheme is used to guarantee a non-singular equation set. The nature of this singularity offers an explanation for the apparent universality of asymmetric stall in planar diffusers. Results for several flows demonstrate the accuracy and range of applicability of the method.

Childs, R. E.↗

Fully coupled simulation of cosmic reionization. I. numerical methods and tests

Here, we describe an extension of the Enzo code to enable fully coupled radiation hydrodynamical simulation of inhomogeneous reionization in large similar to(100 Mpc)(3) cosmological volumes with thousands to millions of point sources. We solve all dynamical, radiative transfer, thermal, and ionization processes self-consistently on the same mesh, as opposed to a postprocessing approach which coarse-grains the radiative transfer. But, we employ a simple subgrid model for star formation which we calibrate to observations. The numerical method presented is a modification of an earlier method presented in Reynolds et al. differing principally in the operator splitting algorithm we use to advance the system of equations. Radiation transport is done in the gray flux-limited diffusion (FLD) approximation, which is solved by implicit time integration split off from the gas energy and ionization equations, which are solved separately. This results in a faster and more robust scheme for cosmological applications compared to the earlier method. The FLD equation is solved using the hypre optimally scalable geometric multigrid solver from LLNL. By treating the ionizing radiation as a grid field as opposed to rays, our method is scalable with respect to the number of ionizing sources, limited only by the parallel scaling properties of the radiation solver. We test the speed and accuracy of our approach on a number of standard verification and validation tests. We show by direct comparison with Enzo's adaptive ray tracing method Moray that the well-known inability of FLD to cast a shadow behind opaque clouds has a minor effect on the evolution of ionized volume and mass fractions in a reionization simulation validation test. Finally, we illustrate an application of our method to the problem of inhomogeneous reionization in a 80 Mpc comoving box resolved with 3200(3) Eulerian grid cells and dark matter particles.

79 ASTRONOMY AND ASTROPHYSICS↗

Sensitivity of the interpretation of the experimental ion thermal diffusivity to the determination of the ion conductive heat flux

A moments equation formalism for the interpretation of the experimental ion thermal diffusivity from experimental data is used to determine the radial ion thermal conduction flux that must be used to interpret the measured data. It is shown that the total ion energy flux must be corrected for thermal and rotational energy convection, for the work done by the flowing plasma against the pressure and viscosity, and for ion orbit loss of particles and energy, and expressions are presented for these corrections. Each of these factors is shown to have a significant effect on the interpreted ion thermal diffusivity in a representative DIII-D [J. Luxon, Nucl. Fusion 42, 614 (2002)] discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗