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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,495 records · Page 83

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Atomistic simulations of the thermal conductivity of liquids

We present a method based on sinusoidal approach to equilibrium molecular dynamics (SAEMD) to compute the thermal conductivity of liquids Similar to nonequilibrium molecular dynamics, and unlike equilibrium simulations based on the Green-Kubo formalism, the method only requires the calculation of forces and total energies. The evaluation of heat fluxes and energy densities is not necessary, thus offering the promise of efficiently implementing first principles simulations based on density functional theory or deep molecular dynamics. Our approach is a generalization of SAEMD for solids, where the thermal conductivity is computed in the steady state, instead of a transient regime, thus properly taking into account diffusive terms in the heat equation. We present results for liquid water at ambient conditions and under pressure and discuss simulation requirements to obtain converged values of the thermal conductivity as a function of size and simulation time.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Small-angle scattering of particle assemblies

Small-angle scattering formulae for crystalline assemblies of arbitrary particles are derived from powder diffraction theory using the decoupling approximation. To do so, the pseudo-lattice factor is defined, and methods to overcome the limitations of the decoupling approximation are investigated. Further, approximated equations are suggested for the diffuse scattering from various defects of the first kind due to non-ideal particles, including size polydispersity, orientational disorder and positional fluctuation about their ideal positions. Calculated curves using the formalism developed herein are compared with numerical simulations computed without any approximation. For a finite-sized assembly, the scattering from the whole domain of the assembly must also be included, and this is derived using the correlation function approach.

assembly↗

DOE OSTI2012

RELAP-7 and PRONGHORN Initial Integration Plan

Modern nuclear reactor safety codes require the ability to solve detailed coupled neutronicthermal fluids problems. For larger cores, this implies fully coupled 3-D spatial dynamics with appropriate feedback models that can provide enough resolution to accurately compute core heat generation and removal during steady and unsteady conditions. The reactor analyis code PRONGHORN is being coupled to RELAP-7 as a first step to extend RELAP's current capabilities. This report details the mathematical models, the type of coupling, and the testing that will be used to produce an integrated system. RELAP-7 is a MOOSE-based application that solves the continuity, momentum, and energy equations in 1-D for a compressible fluid. The pipe and joint capabilities enable it to model parts of the PCU system. The PRONGHORN application, also developed on the MOOSE infrastructure, solves the coupled equations that define the neutron diffusion, fluid flow, and heat transfer in a 3-D core model. Initially, the two systems will be loosely coupled to simplify the transition towards a more complex infrastructure. The integration will be tested with the OECD/NEA MHTGR-350 Coupled Neutronics-Thermal Fluids benchmark model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2012

Initial Coupling of the RELAP-7 and PRONGHORN Applications

Modern nuclear reactor safety codes require the ability to solve detailed coupled neutronic- thermal fluids problems. For larger cores, this implies fully coupled higher dimensionality spatial dynamics with appropriate feedback models that can provide enough resolution to accurately compute core heat generation and removal during steady and unsteady conditions. The reactor analysis code PRONGHORN is being coupled to RELAP-7 as a first step to extend RELAP’s current capabilities. This report details the mathematical models, the type of coupling, and the testing results from the integrated system. RELAP-7 is a MOOSE-based application that solves the continuity, momentum, and energy equations in 1-D for a compressible fluid. The pipe and joint capabilities enable it to model parts of the power conversion unit. The PRONGHORN application, also developed on the MOOSE infrastructure, solves the coupled equations that define the neutron diffusion, fluid flow, and heat transfer in a full core model. The two systems are loosely coupled to simplify the transition towards a more complex infrastructure. The integration is tested on a simplified version of the OECD/NEA MHTGR-350 Coupled Neutronics-Thermal Fluids benchmark model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2013

BISON Theory Manual: The Equations behind Nuclear Fuel Analysis

BISON is a finite element-based nuclear fuel performance code applicable to a variety of fuel forms including light water reactor fuel rods, TRISO particle fuel, and metallic rod and plate fuel. It solves the fully-coupled equations of thermomechanics and species diffusion, for either 1D spherical, 2D axisymmetric or 3D geometries. Fuel models are included to describe temperature and burnup dependent thermal properties, fission product swelling, densification, thermal and irradiation creep, fracture, and fission gas production and release. Plasticity, irradiation growth, and thermal and irradiation creep models are implemented for clad materials. Models are also available to simulate gap heat transfer, mechanical contact, and the evolution of the gap/plenum pressure with plenum volume, gas temperature, and fission gas addition. BISON is based on the MOOSE framework and can therefore efficiently solve problems using standard workstations or very large high-performance computers. This document describes the theoretical and numerical foundations of BISON.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2016

BISON Theory Manual The Equations behind Nuclear Fuel Analysis

BISON is a finite element-based nuclear fuel performance code applicable to a variety of fuel forms including light water reactor fuel rods, TRISO particle fuel, and metallic rod and plate fuel. It solves the fully-coupled equations of thermomechanics and species diffusion, for either 2D axisymmetric or 3D geometries. Fuel models are included to describe temperature and burnup dependent thermal properties, fission product swelling, densification, thermal and irradiation creep, fracture, and fission gas production and release. Plasticity, irradiation growth, and thermal and irradiation creep models are implemented for clad materials. Models are also available to simulate gap heat transfer, mechanical contact, and the evolution of the gap/plenum pressure with plenum volume, gas temperature, and fission gas addition. BISON is based on the MOOSE framework and can therefore efficiently solve problems using standard workstations or very large high-performance computers. This document describes the theoretical and numerical foundations of BISON.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2017

DCMIP2016: a review of non-hydrostatic dynamical core design and intercomparison of participating models

Abstract. Atmospheric dynamical cores are a fundamental component of global atmospheric modeling systems and are responsible for capturing the dynamical behavior of the Earth's atmosphere via numerical integration of the Navier–Stokes equations. These systems have existed in one form or another for over half of a century, with the earliest discretizations having now evolved into a complex ecosystem of algorithms and computational strategies. In essence, no two dynamical cores are alike, and their individual successes suggest that no perfect model exists. To better understand modern dynamical cores, this paper aims to provide a comprehensive review of 11 non-hydrostatic dynamical cores, drawn from modeling centers and groups that participated in the 2016 Dynamical Core Model Intercomparison Project (DCMIP) workshop and summer school. This review includes a choice of model grid, variable placement, vertical coordinate, prognostic equations, temporal discretization, and the diffusion, stabilization, filters, and fixers employed by each system.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1974

User's instructions for the 41-node thermoregulatory model (steady state version)

A user's guide for the steady-state thermoregulatory model is presented. The model was modified to provide conversational interaction on a remote terminal, greater flexibility for parameter estimation, increased efficiency of convergence, greater choice of output variable and more realistic equations for respiratory and skin diffusion water losses.

Leonard, J. I.↗

NASA NTRS1978

Closure of a kinetic model of plasma in strong turbulence by relaxation

A Fokker-Planck kinetic equation for a turbulent plasma is derived by a repeated cascade decomposition. Calculation of the propagator and the kinetic equation determine the transport coefficients (diffusivity and turbulent viscosity) by means of a closure based on a relaxation procedure governing the approach to equilibrium. The k to the minus third power spectral law is obtained, which governs the coupling between the velocity and the electrostatic field fluctuations.

Tchen, C. M.↗

NASA NTRS1979

A kinetic approach to compressible turbulence

A kinetic theory of strong turbulence is developed by decomposing fluctuation of distribution function into scaled components, or ranks, to represent three transport processes: macroscopic evolution, transport property and relaxation. Closure is obtained by formulating a relaxation process for the approach of eddy diffusivity to equilibrium. The kinetic equation is solved and three transport functions of turbulence in adiabatic gas are derived: two transfer functions describe turbulent mode transfer, and a coupling function describes the build-up of kinetic energy at the expense of internal energy. The energy balance in the inertial subrange of homogeneous and isotropic turbulence in an adiabatic gas yields a k to the -3rd law, in agreement with numerical computations from one dimensional model of compressible turbulence.

Tchen, C. M.↗

NASA NTRS1981

Relationships between diffuse reflectance and vegetation canopy variables based on the radiative transfer theory

Out of the lengthy original expression of the diffuse reflectance formula, simple working equations were derived by employing characteristic parameters, which are independent of the canopy coverage and identifiable by field observations. The typical asymptotic nature of reflectance data that is usually observed in biomass studies was clearly explained. The usefulness of the simplified equations was demonstrated by the exceptionally close fit of the theoretical curves to two separately acquired data sets for alfalfa and shortgrass prairie canopies.

Park, J. K.↗

NASA NTRS1982

A simple finite difference procedure for the vortex controlled diffuser

A simple prediction procedure for sudden expansion incompressible flows is developed and applied to the vortex controlled diffuser. Transient Navier-Stokes equations of an incompressible fluid are solved by means of their associated finite difference equations in terms of the primitive pressure velocity variables. A computer code is developed using a laminar flow simulation with free slip or no slip wall boundary conditions. In addition, predicted results confirm that effectiveness increases with increases in duct length and bleed flow rate

Busnaina, A. A.↗

NASA NTRS1982

Beam-plasma instability in the presence of low-frequency turbulence

General equations are derived for a linear beam-plasma instability in the presence of low-frequency turbulence. Within a 'quasi-linear' statistical approximation, these equations contain Langmuir wave scattering, diffusion, resonant and nonresonant anomalous absorption, and a 'plasma laser' effect. It is proposed that naturally occurring density irregularities in the solar wind may stabilize the beam-unstable Langmuir waves which occur during type III solar emissions.

Goldman, M. V.↗

NASA NTRS1983

Transport properties of oxygen

Tables of viscosity, thermal conductivity, and thermal diffusivity of oxygen as a function of temperature and pressure from the triple point to 320 K and at pressures to 100 MPa are presented. Auxiliary tables in engineering units are also given. Viscosity and thermal conductivity are calculated from published correlations. Density and specific heat at constant pressure, required to calculate thermal diffusivity, are obtained from an equation of state. The Prandtl number can be obtained quite easily from the values tabulated.

Roder, H. M.↗

NASA NTRS1992

A burst from a thermonuclear runaway on an ONeMg white dwarf

Studies which examine the consequences of accretion, at rates of 10(exp -9) solar mass/yr and 10(exp -10) solar mass/yr, onto an ONeMg white dwarf with a mass of 1.35 solar masses are performed. In these studies, a Lagrangian, hydrodynamic, one-dimensional computer code was used. The code now includes a network with 89 nuclei up to Ca-40, elemental diffusion, new opacities, and new equation of state. The initial abundance distribution corresponded to a mixture that was enriched to either 25, 50, or 75 percent in products of carbon burning. The remaining material in each case is assumed to have a solar composition. The evolution of the thermonuclear runaway in the 1.35 solar mass white dwarf, with M = 10(exp -9) solar mass, produced peak temperatures in the shell source exceeding 300 million degrees. The sequence produced significant amounts of Na-22 from proton captures onto Ne-20 and significant amounts of Al-26 from proton captures on Mg-24. This sequence ejected 5.2 x 10(exp -6) solar mass moving with speeds from approximately 100 km/s to 2300 km/s. When the mass accretion rate was decreased to 10(exp -10) solar mass, the resulting thermonuclear runaway produced a shock that moved through the outer envelope of the white dwarf and raised the surface luminosity to L greater than 10(exp 7) solar luminosity and the effective temperature to values exceeding 10(exp 7) K. The interaction of the material expanding from off of the white dwarf with the accretion disk should produce a burst of gamma-rays.

Starrfield, S.↗

NASA NTRS1992

The structure of premixed particle-cloud flames

The structure of premixed flames propagating in combustible systems, containing uniformly distributed volatile fuel particles, in an oxidizing gas mixture, is analyzed. It is presumed that the fuel particles vaporize first to yield a gaseous fuel of known chemical structure, which is subsequently oxidized in the gas phase. The analysis is performed in the asymptotic limit, where the value of the characteristic Zeldovich number, based on the gas-phase oxidation of the gaseous fuel is large, and for values of phi(u) greater than or equal to 1.0, where phi(u) is the equivalence ratio based on the fuel available in the fuel particles. The structure of the flame is presumed to consist of a preheat vaporization zone where the rate of the gas-phase chemical reaction is small, a reaction zone where convection and the rate of vaporization of the fuel particles are small and a convection zone where diffusive terms in the conservation equations are small. For given values phi(u) the analysis yields results for the burning velocity and phi(g) where phi(g) is the effective equivalence ratio in the reaction zone. The analysis shows that even though phi(u) greater than or equal to 1.0, for certain cases the calculated value of phi(g) is less than unity. This prediction is in agreement with experimental observations.

Seshadri, K.↗

NASA NTRS1993

Heat Transfer from Radiatively Heated Material in a Low Reynolds Number Microgravity Environment

A mathematical model of the transient three-dimensional heat transfer between a slowly moving ambient gas stream and a thermally thick or thin flat surface heated by external radiation in a microgravity environment is presented. The problem is motivated in part by fire safety issues in spacecraft. The gas phase is represented by variable property convection-diffusion energy and mass conservation equations valid at low Reynolds numbers. The absence of gravity and low Reynolds number together permit the flow to be represented by a self-consistent velocity potential determined by the ambient velocity and the thermal expansion in the gas. The solid exchanges energy with the gas by conduction/convection and with the surroundings by surface absorption and re-emission of radiation. Heat conduction in the solid is assumed to be one dimensional at each point on the surface as a consequence of the limited times (of order of 10 seconds) of interest in these simulations. Despite the apparent simplicity of the model, the results show a complex thermally induced flow near the heated surface. The thermal exchange between the gas and solid produces an outward sourcelike flow upstream of the center of the irradiated area and a sinklike flow downstream. The responses of the temperature fields and the associated flows to changes in the intensity of the external radiation and the ambient velocity are discussed.

Yamashita, H.↗

NASA NTRS1994

Green's formula and variational principles for cosmic-ray transport with application to rotating and shearing flows

Green's theorem and Green's formula for the diffusive cosmic-ray transport equation in relativistic flows are derived. Green's formula gives the solution of the transport equation in terms of the Green's function of the adjoint transport equation, and in terms of distributed sources throughout the region R of interest, plus terms involving the particle intensity and streaming on the boundary. The adjoint transport equation describes the time-reversed particle transport. An Euler-Lagrange variational principle is then obtained for both the mean scattering frame distribution function f, and its adjoint f(dagger). Variations of the variational functional with respect to f(dagger) yield the transport equation, whereas variations of f yield the adjoint transport equation. The variational principle, when combined with Noether's theorem, yields the conservation law associated with Green's theorem. An investigation of the transport equation for steady, azimuthal, rotating flows suggests the introduction of a new independent variable H to replace the comoving frame momentum variable p'. For the case of rigid rotating flows, H is conserved and is shown to be analogous to the Hamiltonian for a bead on a rigidly rotating wire. The variable H corresponds to a balance between the centrifugal force and the particle inertia in the rotating frame. The physical interpretation of H includes a discussion of nonrelativistic and special relativistic rotating flows as well as the cases of aziuthal, differentially rotating flows about Schwarzs-child and Kerr black holes. Green's formula is then applied to the problem of the acceleration of ultra-high-energy cosmic rays by galactic rotation. The model for galactic rotation assumes an angular velocity law Omega = Omega(sub 0)(omega(sub 0)/omega), where omega denotes radial distance from the axis of rotation. Green's functions for the galactic rotation problem are used to investigate the spectrum of accelerated particles arising from monoenergetic and truncated power-law sources. We conclude that it is possible to accelerate particles beyond the knee by galactic rotation, but not in sufficient number to adequately explain the observed spectrum.

Webb, G. M.↗