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 973 records · Page 54

Fast secant methods for the iterative solution of large nonsymmetric linear systems

A family of secant methods based on general rank-1 updates was revisited in view of the construction of iterative solvers for large non-Hermitian linear systems. As it turns out, both Broyden's good and bad update techniques play a special role, but should be associated with two different line search principles. For Broyden's bad update technique, a minimum residual principle is natural, thus making it theoretically comparable with a series of well known algorithms like GMRES. Broyden's good update technique, however, is shown to be naturally linked with a minimum next correction principle, which asymptotically mimics a minimum error principle. The two minimization principles differ significantly for sufficiently large system dimension. Numerical experiments on discretized partial differential equations of convection diffusion type in 2-D with integral layers give a first impression of the possible power of the derived good Broyden variant.

Deuflhard, Peter↗

Focused transport of energetic particles along magnetic field lines draped around a coronal mass ejection

Evidence is presented for focused transport of energetic particles along magnetic field lines draped around a coronal mass ejection. This evidence was obtained with the University of Maryland/Max-Planck-Institute experiment on the ISEE-3 spacecraft during the decay phase of the June 6, 1979, solar particle event. During the early portion of the decay phase of this event, interplanetary magnetic field lines were apparently draped around a coronal mass ejection, leading to a small focusing length on the western flank where ISEE 3 was located. A period of very slow decrease of particle intensity was observed, along with large sunward anisotropy in the solar wind frame, which is inconsistent with predictions of the standard Fokker-Planck equation models for diffusive transport. It was found possible to fit the observations, assuming that focused transport dominates and that the particle pitch angle scattering is isotropic.

Tan, L. C.↗

Relationship Between Solution Shear Viscosity and Density at the Saturation Point

Properties of supersaturated solutions such as the density, viscosity, and solute diffusivity are dependent on the solute concentration. The diffusion-boundary-layer equations are derived and solved for the natural convection case with the viscosity and density dependent on the solute concentration. The solution obtained demonstrates that, at the vicinity of the saturation concentration c(sub s), there is a non-trivial dependence of the solution viscosity eta on its density rho: eta(c(sub s)) = eta(rho(sub s)) varies as rho(sub s)(exp 1/2), where rho(sub s) = rho(c(sub s)). This result has been verified in experiments with aqueous solutions of inorganic and organic salts.

Izmailov, Alexander F.↗

A Three-Wave Model of the Stratosphere with Coupled Dynamics, Radiation and Photochemistry

A zonal mean chemistry transport model (2-D CTM) coupled with a semi-spectral dynamical model is used to simulate the distributions of trace gases in the present day atmosphere. The zonal-mean and eddy equations for the velocity and the geopotential height are solved in the semi-spectral dynamical model. The residual mean circulation is derived from these dynamical variables and used to advect the chemical species in the 2- D CTM. Based on a linearized wave transport equation, the eddy diffusion coefficients for chemical tracers are expressed in terms of the amplitude, frequency and growth rate of dynamical waves; local chemical loss rates; and a time constant parameterizing small scale mixing. The contributions to eddy flux are from the time varying wave amplitude (transient eddy), chemical reactions (chemical eddy) and small scale mixing. In spite of the high truncation in the dynamical module (only three longest waves are resolved), the model has simulated many observed characteristics of stratospheric dynamics and distribution of chemical species including ozone. Compared with the values commonly used in 2-D CTMs, the eddy diffusion coefficients for chemical species calculated in this model are smaller, especially in the subtropics. It is also found that the chemical eddy diffusion has only a small effects in determining the distribution of most slow species, including ozone in the stratosphere.

Shia, Run-Lie↗

Particle Distribution in Suspension Shear Flow

An earlier proposed constitutive relation for normal stresses originated by random particle fluctuations is used to describe a joint effect of thermal and shear-induced fluctuations on concentrational distributions in suspension flow. Averaged products of components of the fluctuation velocity are evaluated on a basis of the rational mechanics approach combined with a simple kinematic consideration. The equation of momentum conservation of the dispersed phase of a suspension closed with this constitutive relation is applied to unidirectional shear flow in the gravity field and to rotational Couette flow. Coupling of the thermal and shear-induced fluctuations results in that the ability of shear flow to suspend particles has a minimum at a certain particle size, all other things being equal. The developed model provides also for a reasonable explanation of particle distributions observed in Couette flow. The approach based on the consideration of momentum balance for the dispersed phase is proved to lead to an effective equation of convective diffusion of the suspended particles. Coefficients of mutual diffusion due to both thermal and shear-induced fluctuations are drastically different from corresponding self-diffusivities as regards both their scaling and their concentrational dependence.

Buyevich, Yu A.↗

A Theoretical Solid Oxide Fuel Cell Model for System Controls and Stability Design

As the aviation industry moves towards higher efficiency electrical power generation, all electric aircraft, or zero emissions and more quiet aircraft, fuel cells are sought as the technology that can deliver on these high expectations. The Hybrid Solid Oxide Fuel Cell system combines the fuel cell with a microturbine to obtain up to 70 percent cycle efficiency, and then distributes the electrical power to the loads via a power distribution system. The challenge is to understand the dynamics of this complex multi-discipline system, and design distributed controls that take the system through its operating conditions in a stable and safe manner while maintaining the system performance. This particular system is a power generation and distribution system and the fuel cell and microturbine model fidelity should be compatible with the dynamics of the power distribution system in order to allow proper stability and distributed controls design. A novel modeling approach is proposed for the fuel cell that will allow the fuel cell and the power system to be integrated and designed for stability, distributed controls, and other interface specifications. This investigation shows that for the fuel cell, the voltage characteristic should be modeled, but in addition, conservation equation dynamics, ion diffusion, charge transfer kinetics, and the electron flow inherent impedance should also be included.

Kopasakis, George↗

A Theoretical Solid Oxide Fuel Cell Model for Systems Controls and Stability Design

As the aviation industry moves toward higher efficiency electrical power generation, all electric aircraft, or zero emissions and more quiet aircraft, fuel cells are sought as the technology that can deliver on these high expectations. The hybrid solid oxide fuel cell system combines the fuel cell with a micro-turbine to obtain up to 70% cycle efficiency, and then distributes the electrical power to the loads via a power distribution system. The challenge is to understand the dynamics of this complex multidiscipline system and the design distributed controls that take the system through its operating conditions in a stable and safe manner while maintaining the system performance. This particular system is a power generation and a distribution system, and the fuel cell and micro-turbine model fidelity should be compatible with the dynamics of the power distribution system in order to allow proper stability and distributed controls design. The novelty in this paper is that, first, the case is made why a high fidelity fuel cell mode is needed for systems control and stability designs. Second, a novel modeling approach is proposed for the fuel cell that will allow the fuel cell and the power system to be integrated and designed for stability, distributed controls, and other interface specifications. This investigation shows that for the fuel cell, the voltage characteristic should be modeled but in addition, conservation equation dynamics, ion diffusion, charge transfer kinetics, and the electron flow inherent impedance should also be included.

Kopasakis, George↗

Introduction to Physical Intelligence

A slight deviation from Newtonian dynamics can lead to new effects associated with the concept of physical intelligence. Non-Newtonian effects such as deviation from classical thermodynamic as well as quantum-like properties have been analyzed. A self-supervised (intelligent) particle that can escape from Brownian motion autonomously is introduced. Such a capability is due to a coupling of the particle governing equation with its own Liouville equation via an appropriate feedback. As a result, the governing equation is self-stabilized, and random oscillations are suppressed, while the Liouville equation takes the form of the Fokker-Planck equation with negative diffusion. Non- Newtonian properties of such a dynamical system as well as thermodynamical implications have been evaluated.

Zak, Michail↗

Quasi-linear velocity space diffusion of heavy cometary pickup ions on bispherical diffusion characteristics

We consider an ion transport equation describing the source, convection, adiabatic acceleration, and quasi-linear velocity diffusion of cometary pickup ions. The equation is solved numerically along solar wind plasma flow lines for the environment of Comet Halley, to obtain distributions at positions on the Giotto spacecraft trajectory which may be compared with observations. We obtain full two-dimensional (pitch angle and velocity) numerical distributions which show pitch angle scattering about the wave scattering centers at +/- V(A), where V(A) is the Alfven wave speed in the solar wind frame. Peak ion phase-space densities approximately follow the bispherical shell geometry. The energy distributions F(v) agree well with Giotto observations at lower v, but in the high-energy tail there is evidence possibly for the first-order Fermi acceleration mechanism taking place in the Comet Halley foreshock region.

Huddleston, D. E.↗

Computations of ion diffusion coefficients from the Boltzmann-Fokker-Planck equation

The Boltzmann-Fokker-Planck equation is solved with the Chapman-Enskog method of analysis for the velocity distribution functions of helium, carbon, nitrogen, and oxygen. The analysis is a perturbation scheme based on the assumption of a collision-dominated gas, and the calculations are carried out to first order. The elements considered are treated as trace constituents in an electron-proton gas. From the resulting distribution functions, diffusion coefficients are computed which are found to be 20-30% less than those obtained by Chapman and Burgers. In addition, it is shown that the return current of cold electrons needed to maintain quasi-neutrality in a plasma with a temperature gradient contributes a term in the thermal diffusion coefficient omitted erroneously in previous works. This added term resolves the longstanding controversy over the discrepancy between the coefficients of Chapman and Burgers, which are seen to be completely equivalent in the light of this analysis. The viscosity coefficient for an electron-proton gas is also computed and found to be 7% less than that obtained by Braginskii.

Roussel-Dupre, R.↗

Diffusion Synthetic Acceleration for High Order S{sub N} Transport on Meshes with Curved Surfaces

Thermal radiation transport is integral to multiphysics applications such as inertial confinement fusion and astrophysics. Materials in these problems can be exceptionally optically thick and diffusive. A transport solver must be capable of accurate and efficient solutions to such problems. Previously, Woods et al. established a proof-of-concept for high order (HO) discontinuous Galerkin finite element (DGFEM) transport on meshes with curved surfaces. They numerically demonstrated the spatial convergence is order p + 1. Remarkably, these convergence rates are relatively invariant to the amount of curvature of the spatial mesh. We present here an extension of this work that demonstrates a rapidly-convergent iterative solver for the FEM transport equations on meshes with curved sides. We use the modified interior penalty (MIP) diffusion synthetic acceleration (DSA) established by Wang and Ragusa. The partially consistent MIP equations stem from the symmetric interior penalty (IP) derivation of the diffusion equation. The IP method is not stable for optically thick cells so it was combined with the diffusion conforming form (DCF), a discretized diffusion equation derived from the discretized S N equations that is not stable for intermediate optical thicknesses. The MIP equations adapt a 'switch' to the IP method in optically thinner regions and the DCF method in optically thicker regions. The literature does not reveal any Fourier analyses performed for iterative transport on meshes with curved surfaces. Wang and Ragusa performed a Fourier analysis for orthogonal meshes and our results appear to be as expected. Previous literature has reported on some of the sensitivities that affect the spectral radius of DSA schemes: polynomial order, and scattering ratios. Presently, we focus on the sensitivity of the spectral radius to polynomial order, scattering ratio, the constant C, and mesh curvature for varying cell thicknesses. We also numerically study the behavior of the spectral radius to verify our DSA implementation. A sequence of problems is designed to exercise the acceleration. As the problems become more optically thick and diffusive, the accelerated solution continues to converge toward the analytic solution. We employ the Modular Finite Element Methods (MFEM) package, developed at Lawrence Livermore National Laboratory, to generate the linear algebraic system of equations for each angle in the discrete ordinates quadrature set, and use UMFPack to solve this system of equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The combined effects of stress and moisture on the structural integrity of composites

A variational approach is used to model the moisture content and distribution profiles which cause swelling in polymers and composites which are used as aircraft structural components. The moisture intrusion stems from humidity, which makes the problem one of diffusion of mass from the environment into a semiinfinite plate with a given diffusion coefficient and thickness. Lagrangian, ordinary differential equations are defined for the mass diffusion, including driving force, mass potential and volume dissipation expressions. Account is taken of the effects of stress on diffusivity, including the dilatational stress caused by inhomogeneous swelling of the material. It is noted that material near the top and bottom of a part undergoes the most swelling, and the evolution of moisture diffusion in the materials is traced. Finally, experimental data are used to define differential equations for predicting the final moisture distribution within the material.

Hiel, C.↗

Modeling the Ionic Strength Effect on Diffusion in Clay. The DR-A Experiment at Mont Terri

Solute diffusion in compacted clays depends on ionic strength through its control on the thickness of the electrical double layer (EDL) on the charged clay surfaces. In the DR-A field experiment (Mont Terri, Switzerland), synthetic porewater was circulated through a borehole for 189 days, leading to the out-diffusion of a variety of tracers into the Opalinus Clay. The borehole solution was then replaced with a higher-salinity solution for an additional 540 days, leading to the diffusion of Cs + , Ca 2+ , Mg 2+ , and Sr 2+ back into the borehole and to an increase in the out-diffusion of anions (I - , Br - ) and 3 H. The experimental results were interpreted using the CrunchClay code, which includes a mean electrostatic potential model for the EDL. The EDL corresponds to a second continuum in addition to bulk electrically neutral porewater. Species-specific diffusion (Nernst-Planck equation) occurs through both domains. A 1D radial model considered a single pore diffusion coefficient (Dp = 10-9 m2/s) for cations and 3H in the bulk porosity, and a smaller D p (3 × 10 -10 m 2 /s) for anions. Dp values in the EDL were smaller (10 -11 m 2 /s), except for Cs + and K + (5 × 10 -10 and 2 × 10 -10 m 2 /s, respectively). The model reproduced well the experimental results and showed the capability to consider temporal changes in geochemical conditions affecting the transport and retention of potentially important radionuclide contaminants (e.g., 137 Cs + , 90 Sr 2+ , 129 I - ) in underground geological nuclear waste repositories. Finally, coupled multicomponent diffusion together with the electrostatic properties of the charged surfaces are essential in the development of predictive models for ion transport in clays.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of concentration-dependent elastic modulus on Li-ions diffusion and diffusion-induced stresses in spherical composition-gradient electrodes

The composition-gradient electrode material is considered as one of the most promising materials for lithium-ion batteries because of its excellent electrochemical performance and thermal stability. In this work, the effects of concentration-dependent elastic modulus on Li-ions diffusion and diffusion-induce stress in the composition-gradient electrodes were studied. The coupling equations of elasticity and diffusion under both potentiostatic charging and galvanostatic charging were developed to obtain the distributions of both the Li-ions concentration and the stress. The results indicated that the effects of the concentration-dependent elastic modulus on the Li-ions diffusion and the diffusion-induce stresses are controlled by the lithiation induced stiffening factor in the composition-gradient electrodes: a low stiffening factor at the center and a high stiffening factor at the surface lead to a significant effect, whereas a high stiffening factor at the center and a low stiffening factor at the surface result in a minimal effect. The results in this work provide guidance for the selection of electrode materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling High-Energy Gamma-Rays from the Fermi Bubbles - Oral Presentation [Slides]

In 2010, the Fermi Bubbles were discovered at the galactic center of the Milky Way. These giant gamma-ray structures, extending 50 degrees in galactic latitude and 20-30 degrees in galactic longitude, were not predicted. We wish to develop a model for the gamma-ray emission of the Fermi Bubbles. To do so, we assume that second order Fermi acceleration is responsible for the high-energy emission of the bubbles. Second order Fermi acceleration requires charged particles and irregular magnetic fields—both of which are present in the disk of the Milky Way galaxy. I use the assumption of second order Fermi acceleration in the transport equation, which describes the diffusion of particles. By solving the steady-state case of the transport equation, I compute the proton spectrum due to Fermi second order acceleration and compare this analytical solution to a numerical solution provided by Dr. P. Mertsch. Analytical solutions to the transport equation are taken from Becker, Le, & Dermer and are used to further test the numerical solution. I find that the numerical solution converges to the analytical solution in all cases. Thus, we know the numerical solution accurately calculates the proton spectrum. The gamma-ray spectrum follows the proton spectrum, and will be computed in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

Computation of turbulent boundary layer flows with an algebraic stress turbulence model

An algebraic stress turbulence model is presented, characterized by the following: (1) the eddy viscosity expression is derived from the Reynolds stress turbulence model; (2) the turbulent kinetic energy dissipation rate equation is improved by including a production range time scale; and (3) the diffusion coefficients for turbulence equations are adjusted so that the kinetic energy profile extends further into the free stream region found in most experimental data. The turbulent flow equations were solved using a finite element method. Examples include: fully developed channel flow, fully developed pipe flow, flat plate boundary layer flow, plane jet exhausting into a moving stream, circular jet exhausting into a moving stream, and wall jet flow. Computational results compare favorably with experimental data for most of the examples considered. Significantly improved results were obtained for the plane jet flow, the circular jet flow, and the wall jet flow; whereas the remainder are comparable to those obtained by finite difference methods using the standard kappa-epsilon turbulence model. The latter seems to be promising with further improvement of the expression for the eddy viscosity coefficient.

Kim, Sang-Wook↗

Transport Phenomena During Equiaxed Solidification of Alloys

Recent progress in modeling of transport phenomena during dendritic alloy solidification is reviewed. Starting from the basic theorems of volume averaging, a general multiphase modeling framework is outlined. This framework allows for the incorporation of a variety of microscale phenomena in the macroscopic transport equations. For the case of diffusion dominated solidification, a simplified set of model equations is examined in detail and validated through comparisons with numerous experimental data for both columnar and equiaxed dendritic growth. This provides a critical assessment of the various model assumptions. Models that include melt flow and solid phase transport are also discussed, although their validation is still at an early stage. Several numerical results are presented that illustrate some of the profound effects of convective transport on the final compositional and structural characteristics of a solidified part. Important issues that deserve continuing attention are identified.

Beckermann, C.↗

Energy-momentum-conserving stochastic differential equations and algorithms for the nonlinear Landau-Fokker-Planck equation

Coulomb collision is a fundamental diffusion process in plasmas that can be described by the Landau-Fokker-Planck (LFP) equation or the stochastic differential equation (SDE). While energy and momentum are conserved exactly in the LFP equation, they are conserved only on average by the conventional corresponding SDEs, suggesting that the underlying stochastic process may not be well defined by such SDEs. Here, in this study, we derive new SDEs with exact energy-momentum conservation for the Coulomb collision by factorizing the collective effect of field particles into individual particles and enforcing Newton's third law. These SDEs, when interpreted in the Stratonovich sense, have a particularly simple form that represents pure diffusion between particles without drag. To demonstrate that the new SDEs correspond to the LFP equation, we develop numerical algorithms that converge to the SDEs and preserve discrete conservation laws. Simulation results are presented in a benchmark of various relaxation processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗