Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

An analytical solution of the solution-diffusion-electromigration equations reproduces trends in ion rejections during nanofiltration of mixed electrolytes

Here, this paper develops an analytical solution to the equations governing the solution-diffusion-electromigration transport of the ions from two salts that contain a common ion. The analytical expressions, which rely on constant ion permeances, readily enable simple spreadsheet computations of the permeate ion concentrations in nanofiltration (NF). Despite the model simplicity, calculations of ion rejections as a function of transmembrane volume fluxes, feed-solution compositions, and permeance values reproduce experimental trends, with ion permeances as the only adjustable parameters. For example, with solutions containing dissolved MA and M 2 B, when the membrane permeance to B 2– is much lower than the permeances to A – and M + , plots of A – rejection versus transmembrane volume flux show the negative minimum that is often characteristic of NF. Moreover, B 2– rejection increases with increasing MA in the feed solution, reflecting a decreasing electric field in the membrane. The analytical solution also enables spreadsheet computations of electric fields in membrane separation layers and effectively models experimental data for intrinsic rejections as a function of transmembrane volume flow. Interestingly, the model suggests that with sufficiently high permeances of A – (relative to M + ), at certain feed compositions the A – /B 2– selectivity may increase without bound as a function of volume flow. Otherwise, this selectivity does not vary greatly with feed composition and quickly reaches a limit as the volume flow increases. Despite the simplifying approximation of constant ion permeances, using only a few adjustable parameters this model readily reveals trends in NF ion rejections and may help to identify the desired ion permeances or solution compositions for specific separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

A new nonlinear diffusion formalism in a magnetized plasma - Application to space physics and astrophysics

A novel diffusion formalism that takes into account the finite width of resonances is presented. The resonance diagram technique is shown to reproduce the details of the particle orbits very accurately, and can be used to determine the acceleration/scattering in the presence of a given wave spectrum. Ways in which the nonlinear orbits can be incorporated into the diffusion equation are shown. The resulting diffusion equation is an extension of the Q-L theory to cases where the waves have large amplitudes and/or are coherent. This new equation does not have a gap at 90 deg in cases where the individual orbits can cross the gap. The conditions under which the resonance gap at 90-deg pitch angle exits are also examined.

Karimbadi, H.↗
NASA NTRS1993

The telegraph equation in charged particle transport

We present a new derivation of the telegraph equation which modifies its coefficients. First, an infinite order partial differential equation is obtained for the velocity space solid angle-averaged phase-space distribution of particles which underwent at least a few collisions. It is shown that, in the lowest order asymptotic expansion, this equation simplifies to the well-known diffusion equation. The second-order asymptotic expansion for isotropic small-angle scattering results in a modified telegraph equation with a signal propagation speed of v(5/11) exp 1/2 instead of the usual v/3 exp 1/2. Our derivation of a modified telegraph equation follows from an expansion of the Boltzmann equation in the relevant smallness parameters and not from a truncation of an eigenfunction expansion. This equation is consistent with causality. It is shown that, under steady state conditions in a convecting plasma, the telegraph equation may be regarded as a diffusion equation with a modified transport coefficient, which describes a combination of diffusion and cosmic-ray inertia.

Gombosi, T. I.↗
DOE OSTI2016

Generalized theory of diffusion based on kinetic theory

Here, we propose to use spin hydrodynamics, a two-fluid model of spin propagation, as a generalization of the diffusion equation. We show that in the dense limit spin hydrodynamics reduces to Fick's law and the diffusion equation. In the opposite limit spin hydrodynamics is equivalent to a collisionless Boltzmann treatment of spin propagation. Spin hydrodynamics avoids unphysical effects that arise when the diffusion equation is used to describe to a strongly interacting gas with a dilute corona. We apply spin hydrodynamics to the problem of spin diffusion in a trapped atomic gas. We find that the observed spin relaxation rate in the high-temperature limit [Sommer et al., Nature (London) 472, 201 (2011)] is consistent with the diffusion constant predicted by kinetic theory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Finite Volume Method Based Neutronic Solvers for Steady and Transient Analysis of Molten Salt Reactors

In recent years, the coupled neutronics and thermal hydraulics analysis becomes more and more important, because it helps to identify the most relevant safety issues without conservative assumptions. In addition, such coupled simulations are essential for the generation IV advanced nuclear power systems such as molten salt reactor (MSR), because the presence of strong coupling between different phenomena, like the transport of delayed neutron precursors by the fluid flow?. Due to the different discretization methods used in reactor physics analysis and thermal-hydraulics analysis, it is hard to couple different kinds of codes. This work presents an attempt to solve the multigroup neutron diffusion equation with the finite volume method which is quite widely used for the numerical simulation of a variety of applications involving fluid flow and heat and mass transfer. The unstructured finite volume method (FVM) based on open source C++ library OpenFOAM was chosen as the development platform?. The objective of this paper is to concisely present the main features of the developed codes in terms of the theoretical background and numerical approach and as can be seen in the following sections, the calculation results obtained by the NDSFoam which solves the steady-state neutron diffusion equation and NDTFoam which solves the transient-state neutron diffusion equation are discussed and compared to the benchmark problems. In this paper, two new neutronic solvers using finite volume method have been developed to perform the neutron physics analysis. The solvers are based on OpenFOAM, and have been tested with several steady state and transient-state benchmark problems. By verifications of the codes, it is proved that the development of NDSFoam and NDTFoam are correct. In our future work, both of these two solvers will be coupled with the precompiled CFD solvers in OpenFOAM to analyze the molten salt fast reactors. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

The Effects of Interplanetary Transport in the Event-intergrated Solar Energetic Particle Spectra

Previous investigations on the energy spectra of solar energetic particle (SEP) events revealed that the energy spectra observed at 1 au often show double power laws with break energies from one to tens of MeV/nuc. In order to determine whether the double power-law features result from the SEP source or the interplanetary transport process from the Sun to 1 au, we separately analyze the SEP spectra in the decay phase, during which the transport effect is minimum. In this paper, we reported three events observed by the Interplanetary Monitory Platform 8 spacecraft, which occurred on 1977 September 19, November 22, and 1979 March 1. For the first two events, the event-integrated spectra of protons possess double power-law profiles with break energies in a range of several MeV to tens of MeV, while the spectra integrated in the decay (reservoir) phase yield single power laws. Moreover, a general trend from a double power law at the rising phase to a single power law at the decay phase is observed. For the third event, both the event-integrated and the reservoir spectra show double power-law features. However, the difference between the low- and high-energy power-law indices is smaller for the reservoir spectrum than the event-integrated spectrum. These features were reproduced by solving the 1D diffusion equation analytically and we suggest that the transport process, especially the diffusion process, plays an important role in breaking the energy spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A multidimensional unified gas-kinetic scheme for radiative transfer equations on unstructured mesh

In order to extend the unified gas kinetic scheme (UGKS) to solve radiative transfer equations in a complex geometry, a multidimensional asymptotic preserving implicit method on unstructured mesh is constructed in this paper. With an implicit formulation, the CFL condition for the determination of the time step in UGKS can be much relaxed, and a large time step is used in simulations. Differently from previous direction-by-direction UGKS on orthogonal structured mesh, on unstructured mesh the interface flux transport takes into account multi-dimensional effect, where gradients of radiation intensity and material temperature in both normal and tangential directions of a cell interface are included in the flux evaluation. The multiple scale nature makes the UGKS be able to capture the solutions in both optically thin and thick regions seamlessly. In the optically thick region the condition of cell size being less than photon's mean free path is fully removed, and the UGKS recovers a solver for diffusion equation in such a limit on unstructured mesh. For a distorted quadrilateral mesh, the UGKS goes to a nine-point scheme for the diffusion equation, and it naturally reduces to the standard five-point scheme for a orthogonal quadrilateral mesh. Numerical computations covering a wide range of transport regimes on unstructured and distorted quadrilateral meshes will be presented to validate the current approach. - Highlights: • A unified gas-kinetic scheme (UGKS) on unstructured mesh has been developed in all flow regimes. • The scheme solves diffusion equation in optical thick region with large cell size. • The implicit UGKS removes the CFD condition. The efficiency of the scheme has been much improved.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A Robust Numerical Treatment of Solid-Phase Diffusion in Pseudo Two-Dimensional Lithium-Ion Battery Models

Solid-phase diffusion in active materials of lithium-ion batteries significantly affects charging and safety-related behavior of lithium-ion batteries. Therefore, it is essential to develop an efficient and robust numerical algorithm for solving solid-phase diffusion equations in physics-based battery models. In this work, we discuss the origins of numerical instabilities that can occur when solving the solid-phase diffusion equations using iterative methods. Then, in order to resolve such issues, we propose a simple numerical treatment to the surface flux term of discretized solid-phase diffusion equations. To demonstrate its numerical robustness, the proposed method is implemented into a pseudo two-dimensional (P2D) physics-based battery model and simulations are conducted at wide ranges of operating conditions. Even with extremely poor initial guesses for the Li+ concentrations of the active materials, computations using the proposed method do not diverge and the their computational speeds are comparable to those with conventional initial guesses. Comprehensive tests of the proposed method are also performed with a dynamic current profile based on US06 driving profile and a multi-stage charging profile with very high initial C-rate (12C).

battery modeling↗
NASA NTRS2002

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one- dimensional convection-diffusion-reaction (CDR) equations. Accuracy, stability, conservation, and dense-output are first considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, (N = 2), additive Runge-Kutta (ARK(sub 2)) methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms of the partitioned method are of equal order to those of the elemental methods. Derived ARK(sub 2) methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z['] yields -infinity, and retain high stability efficiency in the absence of stiffness, z['] yield 0. Extrapolation-type stage- value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK(sub 2) error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK(sub 2) methods.

Kennedy, Christopher A.↗
NASA NTRS2001

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one-dimensional convection-diffusion-reaction (CDR) equations. First, accuracy, stability, conservation, and dense output are considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, N = 2, additive Runge-Kutta ARK2 methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms are of equal order to those of the elemental methods. Derived ARK2 methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z(exp [I]) goes to infinity, and retain high stability efficiency in the absence of stiffness, z(exp [I]) goes to zero. Extrapolation-type stage-value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK2 error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK2 methods while offering error control and stage-value predictors.

Kennedy, Christopher A.↗
DOE OSTI2017

Calculation of Higher Eigenpairs of the Transport Equation Using IRAM Based on Domain Decomposition

The higher eigenpairs of the neutron transport (or diffusion) equation are useful in core flux distribution reconstruction, transient analysis, stability analysis, and loading pattern optimization. There are various methods for the calculation of higher eigenpairs, such as the power iteration method with filtering technique, the Generalized Modified Power Method (GMPM), the power iteration method with innovative boundary conditions for symmetric cores, and Krylov subspace methods. In regard to Krylov subspace methods, the Implicitly Restarted Arnoldi Method (IRAM) and ORTHOMIN were implemented to calculate higher eigenpairs for both transport and diffusion equations. In Warsa's work, however, there was a significant drawback that problems with energy upscatter could not be treated efficiently. In Modak's work, the mode to which his method was likely to converge depended on the shape of initial guess flux, and only one mode could be computed each time, making it impractical for real use. Furthermore, the computational burden for higher eigenpairs is generally much larger than for only the first eigenpair, especially for neutron transport equations. Therefore, the implementation of parallel computation for higher eigenpairs is essential to improve the efficiency, and no such work has been reported in literature so far. In this summary, we propose an innovative method to calculate higher eigenpairs of the neutron transport equation based on high performance computation. In this method, 2D Matrix Method of Characteristic (MMOC) is used to discretize the neutron transport equation; IRAM is used to evaluate higher eigenpairs; spatial domain decomposition is used to realize parallel computation. The highlight of this method is using the MGMD PGMRES (Multi-Group and Multi-Domain coupled Parallel Generalized Minimal Residual Method) algorithm to simultaneously accelerate the upscatter iterations and the inner interface angular fluxes iterations, avoiding the drawback in Warsa's method and the efficiency problem of conventional Jacobi iteration. To the best of our knowledge, our work represents the first application of IRAM based on domain decomposition for higher eigenpairs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1976

A strictly Markovian expansion for plasma turbulence theory

The collision operator that appears in the equation of motion for a particle distribution function that was averaged over an ensemble of random Hamiltonians is non-Markovian. It is non-Markovian in that it involves a propagated integral over the past history of the ensemble averaged distribution function. All formal expansions of this nonlinear collision operator to date preserve this non-Markovian character term by term yielding an integro-differential equation that must be converted to a diffusion equation by an additional approximation. An expansion is derived for the collision operator that is strictly Markovian to any finite order and yields a diffusion equation as the lowest nontrivial order. The validity of this expansion is seen to be the same as that of the standard quasilinear expansion.

Jones, F. C.↗
NASA NTRS1978

A strictly Markovian expansion for plasma turbulence theory

The collision operator that appears in the equation of motion for a particle distribution function that has been averaged over an ensemble of random Hamiltonians is non-Markovian. It is non-Markovian in that it involves a propagated integral over the past history of the ensemble averaged distribution function. All formal expansions of this nonlinear collision operator to date preserve this non-Markovian character term by term yielding an integro-differential equation that must be converted to a diffusion equation by an additional approximation. In this note we derive an expansion of the collision operator that is strictly Markovian to any finite order and yields a diffusion equation as the lowest non-trivial order. The validity of this expansion is seen to be the same as that of the standard quasi-linear expansion.

Jones, F. C.↗
NASA NTRS1993

Simulations of diffusion-reaction equations with implications to turbulent combustion modeling

An enhanced diffusion-reaction reaction system (DRS) is proposed as a statistical model for the evolution of multiple scalars undergoing mixing and reaction in an isotropic turbulence field. The DRS model is close enough to the scalar equations in a reacting flow that other statistical models of turbulent mixing that decouple the velocity field from scalar mixing and reaction (e.g. mapping closure model, assumed-pdf models) cannot distinguish the model equations from the original equations. Numerical simulations of DRS are performed for three scalars evolving from non-premixed initial conditions. A simple one-step reversible reaction is considered. The data from the simulations are used (1) to study the effect of chemical conversion on the evolution of scalar statistics, and (2) to evaluate other models (mapping-closure model, assumed multivariate beta-pdf model).

Girimaji, Sharath S.↗
DOE OSTI2022

The “Fresnel Equations” for Diffuse radiation on Inclined photovoltaic Surfaces (FEDIS)

Here the well-known Fresnel equations solve for the reflection and transmission of light for precise incident angles. The transmission of diffuse radiation incident on a planar or domed surface is often needed for real-world applications. Due to the complexity of the Fresnel equations, the analytical solution of the integration has hitherto been unobtainable over the last centuries. Therefore, this problem was numerically solved by integrating the angular transmittances in space often leading to substantial computing burden and bias in the results. To efficiently estimate the solar energy resource for a glass-covered photovoltaic (PV) module, we derive an analytical solution of diffuse transmission based on the rigorous integration of an alternate form of the Fresnel equations. The approach leads to a simple yet accurate relative transmittance model that reconciles the solar energy sensed by pyranometers and PV panels. With limited and clearly stated approximations, the complex mathematical derivation resulted in an elegant solution. An experiment using 1-year of data at the National Renewable Energy Laboratory's (NREL's) Solar Radiation Research Laboratory (SRRL) shows that the new model dramatically decreases the disparity between the solar radiation measurements by a Kipp and Zonen CM Pyranometer 22 (CMP22) and an IMT reference cell on a 1-axis tracking system. The solution in this paper can be widely used in scientific and engineering research, development, and applications wherever the Fresnel equations are used.

14 SOLAR ENERGY↗
DOE OSTI2020

Efficient smoothed particle radiation hydrodynamics I: Thermal radiative transfer

This work presents efficient solution techniques for radiative transfer in the smoothed particle hydrodynamics discretization. Two choices that impact efficiency are how the material and radiation energy are coupled, which determines the number of iterations needed to converge the emission source, and how the radiation diffusion equation is solved, which must be done in each iteration. The coupled material and radiation energy equations are solved using an inexact Newton iteration scheme based on nonlinear elimination, which reduces the number of Newton iterations needed to converge within each time step. During each Newton iteration, the radiation diffusion equation is solved using Krylov iterative methods with a multigrid preconditioner, which abstracts and optimizes much of the communication when running in parallel. The code is verified for an infinite medium problem, a one-dimensional Marshak wave, and a two and three-dimensional manufactured problem, and exhibits first-order convergence in time and second-order convergence in space. For these problems, the number of iterations needed to converge the inexact Newton scheme and the diffusion equation is independent of the number of spatial points and the number of processors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Kinetic Models for Topological Nearest-Neighbor Interactions

We consider systems of agents interacting through topological interactions. These have been shown to play an important part in animal and human behavior. Precisely, the system consists of a finite number of particles characterized by their positions and velocities. At random times a randomly chosen particle, the follower, adopts the velocity of its closest neighbor, the leader. We study the limit of a system size going to infinity and, under the assumption of propagation of chaos, show that the limit kinetic equation is a non-standard spatial diffusion equation for the particle distribution function. We also study the case wherein the particles interact with their K closest neighbors and show that the corresponding kinetic equation is the same. Finally, we prove that these models can be seen as a singular limit of the smooth rank-based model previously studied in Blanchet and Degond (J Stat Phys 163:41–60, 2016). The proofs are based on a combinatorial interpretation of the rank as well as some concentration of measure arguments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Lumping Techniques for Mixed Finite Element Diffusion Discretizations

We are interested in positivity preserving spatial discretizations applicable to arbitrary order mixed finite element spatial discretizations of the radiative diffusion equations. Though not guaranteed to yield strict positivity, matrix lumping is appealing as a method of increasing solution positivity because unlike ad hoc fixups or strictly non-negative, non-linear solution representations, matrix lumping keeps the mixed finite element diffusion equations as a linear system of equations. Self-lumping schemes have previously been shown to be a generalization of more traditional matrix lumping techniques. Self-lumping schemes restrict integration quadrature points to the interpolation points used to define interpolatory finite element basis functions. As shown previously, the positivity and order of accuracy of self-lumping schemes are directly related to the choice of finite element interpolation point. Here, we develop self-lumping schemes for arbitrary order, mixed finite element spatial discretizations of the diffusion equation in Cartesian geometry. After introducing several possible self-lumping schemes appropriate for mixed finite elements, we test the ability of each scheme to increase solution positivity relative to an unlumped discretization using a simple test problem designed to induce negative solutions, evaluating efficacy as a function of finite element order and finite element interpolation point type. The method of manufactured solutions is then used to evaluate order of convergence for promising self-lumping schemes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗