Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion equation”

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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

An Equation of State and Compendium of Thermophysical Properties of Liquid Tin, a Prospective Plasma-Facing Material

In this work we collect and summarize a variety of thermophysical properties of liquid tin, which is a prospecive plasma-facing material in fusion reactors. These include the thermodynamic properties density, specific heat, sound spead, and vapor pressure; transport properties such as thermal con- ductivity, viscosity, and surface tension; and hydrogen isotope solubility and diffusivity. The thermodyamic property data are used to fit an equation of state that accurately describes these, and which has been implemented in MELCOR for fusion in order to model systems featuring liquid tin PFCs and coolant loops. Few data on hydrogen solubility could be found, and these differ by some orders of magnitude; diffusivity appears not to have been measured precisely and for subsequent analyses theoretical estimates will have to be relied upon.

Divertors↗
NASA NTRS2001

Nonlinearity Effects of Lateral Density Diffusion Coefficient on Gain-Guided VCSEL Performance

Electron and hole diffusions in the plane of semiconductor quantum wells play an important part in the static and dynamic operations of semiconductor lasers. In this paper, we apply a hydrodynamic model developed from the semiconductor Bloch equations to numerically study the effects of nonlinearity in the diffusion coefficient on single mode operation and direct modulation of a gain-guided InGaAs/GaAs multiple quantum well laser, operating not too far from threshold. We found that a small diffusion coefficient is advantageous for lowering the threshold current and increasing the modulation bandwidth. Most importantly, the effects of nonlinearity in the coefficient can be approximately reproduced by replacing the coefficient with an effective constant diffusion coefficient, which corresponds roughly to the half height density of the density distribution.

Li, Jian-Zhong↗
NASA NTRS1993

On the structure of gaseous confined laminar diffusion flames: Numerical investigation

The structure and characteristics of gaseous confined laminar diffusion flames are investigated by numerically solving the time-dependent two-dimensional axisymmetric conservation equations. The numerical model accounts for the important chemical and physical processes involved, including axial diffusion, viscous effects, radial convection, and finite-rate chemistry. The numerical results clearly show that the flame has a finite thickness and leakage of fuel vapor into the flame zone is possible. The effect of heat release is found to induce some radial flow. Predicted flame shape and dimensions are compared to the classical Burke-Schumann flame. The numerically calculated flame is observed to be about 15 percent taller and 5 percent narrower than that of the Burke-Schumann solution under the same conditions.

Mawid, M. A.↗
NASA NTRS2012

Mass and Momentum Transport in Microcavities for Diffusion-Dominant Cell Culture Applications

For the informed design of microfluidic devices, it is important to understand transport phenomena at the microscale. This letter outlines an analytically-driven approach to the design of rectangular microcavities extending perpendicular to a perfusion microchannel for microfluidic cell culture devices. We present equations to estimate the spatial transition from advection- to diffusion-dominant transport inside cavities as a function of the geometry and flow conditions. We also estimate the time required for molecules, such as nutrients or drugs to travel from the microchannel to a given depth into the cavity. These analytical predictions can facilitate the rational design of microfluidic devices to optimize and maintain long-term, physiologically-based culture conditions with low fluid shear stress.

Yew, Alvin G.↗
DOE OSTI2023

Galactic Bar Resonances with Diffusion: An Analytic Model with Implications for Bar–Dark Matter Halo Dynamical Friction

Abstract The secular evolution of disk galaxies is largely driven by resonances between the orbits of “particles” (stars or dark matter) and the rotation of non-axisymmetric features (spiral arms or a bar). Such resonances may also explain kinematic and photometric features observed in the Milky Way and external galaxies. In simplified cases, these resonant interactions are well understood: for instance, the dynamics of a test particle trapped near a resonance of a steadily rotating bar is easily analyzed using the angle-action tools pioneered by Binney, Monari, and others. However, such treatments do not address the stochasticity and messiness inherent to real galaxies—effects that have, with few exceptions, been previously explored only with complex N -body simulations. In this paper, we propose a simple kinetic equation describing the distribution function of particles near an orbital resonance with a rigidly rotating bar, allowing for diffusion of the particles’ slow actions. We solve this equation for various values of the dimensionless diffusion strength Δ, and then apply our theory to the calculation of bar–halo dynamical friction. For Δ = 0, we recover the classic result of Tremaine and Weinberg that friction ultimately vanishes, owing to the phase mixing of resonant orbits. However, for Δ > 0, we find that diffusion suppresses phase mixing, leading to a finite torque. Our results suggest that stochasticity—be it physical or numerical—tends to increase bar–halo friction, and that bars in cosmological simulations might experience significant artificial slowdown, even if the numerical two-body relaxation time is much longer than a Hubble time.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1984

The well-posedness of the Kuramoto-Sivashinsky equation

The Kuramoto-Sivashinsky equation arises in a variety of applications, among which are modeling reaction diffusion systems, flame propagation and viscous flow problems. It is considered here, as a prototype to the larger class of generalized Burgers equations: those consist of a quadratic nonlinearity and an arbitrary linear parabolic part. It is shown that such equations are well posed, thus admitting a unique smooth solution, continuously dependent on its initial data. As an attractive alternative to standard energy methods, existence and stability are derived in this case, by patching in the large short time solutions without loss of derivatives.

Tadmor, E.↗
NASA NTRS1986

The well-posedness of the Kuramoto-Sivashinsky equation

The Kuramoto-Sivashinsky equation arises in a variety of applications, among which are modeling reaction diffusion systems, flame propagation and viscous flow problems. It is considered here, as a prototype to the larger class of generalized Burgers equations: those consist of a quadratic nonlinearity and an arbitrary linear parabolic part. It is shown that such equations are well posed, thus admitting a unique smooth solution, continuously dependent on its initial data. As an attractive alternative to standard energy methods, existence and stability are derived in this case, by patching in the large short time solutions without 'loss of derivatives'.

Tadmor, E.↗
DOE OSTI2016

Multi-scale diffuse interface modeling of multi-component two-phase flow with partial miscibility

In this paper, we introduce a diffuse interface model to simulate multi-component two-phase flow with partial miscibility based on a realistic equation of state (e.g. Peng–Robinson equation of state). Because of partial miscibility, thermodynamic relations are used to model not only interfacial properties but also bulk properties, including density, composition, pressure, and realistic viscosity. As far as we know, this effort is the first time to use diffuse interface modeling based on equation of state for modeling of multi-component two-phase flow with partial miscibility. In numerical simulation, the key issue is to resolve the high contrast of scales from the microscopic interface composition to macroscale bulk fluid motion since the interface has a nanoscale thickness only. To efficiently solve this challenging problem, we develop a multi-scale simulation method. At the microscopic scale, we deduce a reduced interfacial equation under reasonable assumptions, and then we propose a formulation of capillary pressure, which is consistent with macroscale flow equations. Moreover, we show that Young–Laplace equation is an approximation of this capillarity formulation, and this formulation is also consistent with the concept of Tolman length, which is a correction of Young–Laplace equation. At the macroscopical scale, the interfaces are treated as discontinuous surfaces separating two phases of fluids. Our approach differs from conventional sharp-interface two-phase flow model in that we use the capillary pressure directly instead of a combination of surface tension and Young–Laplace equation because capillarity can be calculated from our proposed capillarity formulation. A compatible condition is also derived for the pressure in flow equations. Furthermore, based on the proposed capillarity formulation, we design an efficient numerical method for directly computing the capillary pressure between two fluids composed of multiple components. Finally, numerical tests are carried out to verify the effectiveness of the proposed multi-scale method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

First and second order operator splitting methods for the phase field crystal equation

In this paper, we present operator splitting methods for solving the phase field crystal equation which is a model for the microstructural evolution of two-phase systems on atomic length and diffusive time scales. A core idea of the methods is to decompose the original equation into linear and nonlinear subequations, in which the linear subequation has a closed-form solution in the Fourier space. We apply a nonlinear Newton-type iterative method to solve the nonlinear subequation at the implicit time level and thus a considerably large time step can be used. By combining these subequations, we achieve the first- and second-order accuracy in time. We present numerical experiments to show the accuracy and efficiency of the proposed methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Neutron matter from local chiral effective field theory interactions at large cutoffs

Neutron matter is an important many-body system that provides valuable constraints for the equation of state (EOS) of neutron stars. Neutron-matter calculations employing chiral effective field theory (EFT) interactions have been extensively used for this purpose. Among the various many-body methods, quantum Monte Carlo (QMC) methods stand out due to their nonperturbative nature and the achievable precision. However, QMC methods require local interactions as input, which leads to the appearance of stronger regulator artifacts compared to nonlocal interactions. To circumvent this, we employ large-cutoff interactions derived within chiral EFT (400 MeV ≤ Λ 𝑐 ≤ 700MeV) for studies of pure neutron matter. These interactions have been adjusted to nucleon-nucleon scattering phase shifts, the triton binding energy, as well as the triton 𝛽-decay half-life. We find that regulator artifacts significantly decrease with increasing cutoff, leading to a significant reduction of uncertainties in the neutron-matter EOS. We discuss implications for the symmetry energy and demonstrate how our new calculations lead to a reduction in the theoretical uncertainty of predicted neutron-star radii by up to 30% for low-mass stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Removing Numerical Pathologies in a Turbulence Parameterization Through Convergence Testing

Abstract Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e., clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, nonlinear artificial diffusion employed for improving numerical stability can introduce unintended small‐scale features into the solution of the model equations. Smoothing the limiters and using linear artificial diffusion reduces the noise and restores the expected first‐order convergence in CLUBB's solutions. These model reformulations enhance our confidence in the trustworthiness of solutions from CLUBB by eliminating the unphysical oscillations in high‐resolution simulations. The improvements in the results at coarser, near‐operational grid spacing and timestep are also seen in cumulus cloud and dry turbulence tests. In addition, convergence testing is proven to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

58 GEOSCIENCES↗
DOE OSTI2016

Report on simulation of fission gas and fission product diffusion in UO 2

In UO 2 nuclear fuel, the retention and release of fission gas atoms such as xenon (Xe) are important for nuclear fuel performance by, for example, reducing the fuel thermal conductivity, causing fuel swelling that leads to mechanical interaction with the clad, increasing the plenum pressure and reducing the fuel–clad gap thermal conductivity. We use multi-­scale simulations to determine fission gas diffusion mechanisms as well as the corresponding rates in UO 2 under both intrinsic and irradiation conditions. In addition to Xe and Kr, the fission products Zr, Ru, Ce, Y, La, Sr and Ba have been investigated. Density functional theory (DFT) calculations are used to study formation, binding and migration energies of small clusters of Xe atoms and vacancies. Empirical potential calculations enable us to determine the corresponding entropies and attempt frequencies for migration as well as investigate the properties of large clusters or small fission gas bubbles. A continuum reaction-­diffusion model is developed for Xe and point defects based on the mechanisms and rates obtained from atomistic simulations. Effective fission gas diffusivities are then obtained by solving this set of equations for different chemical and irradiation conditions using the MARMOT phase field code. The predictions are compared to available experimental data. The importance of the large Xe U3O cluster (a Xe atom in a uranium + oxygen vacancy trap site with two bound uranium vacancies) is emphasized, which is a consequence of its high mobility and high binding energy. We find that the Xe U3O cluster gives Xe diffusion coefficients that are higher for intrinsic conditions than under irradiation over a wide range of temperatures. Under irradiation the fast-­moving Xe U3O cluster recombines quickly with irradiation-induced interstitial U ions, while this mechanism is less important for intrinsic conditions. The net result is higher concentration of the Xe U3O cluster for intrinsic conditions than under irradiation. We speculate that differences in the irradiation conditions and their impact on the Xe U3O cluster can explain the wide range of diffusivities reported in experimental studies. However, all vacancy-­mediated mechanisms underestimate the Xe diffusivity compared to the empirical radiation-­enhanced rate used in most fission gas release models. We investigate the possibility that diffusion of small fission gas bubbles or extended Xe-­vacancy clusters may give rise to the observed radiation-­enhanced diffusion coefficient. These studies highlight the importance of U divacancies and an octahedron coordination of uranium vacancies encompassing a Xe fission gas atom. The latter cluster can migrate via a multistep mechanism with a rather low effective barrier, which together with irradiation-induced clusters of uranium vacancies, gives rise to the irradiation-enhanced diffusion coefficient observed in experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

ICAT: A numerical scheme to minimize numerical diffusion in advection-dispersion modeling and its application in identifying flow channeling

The advection-dispersion equation for scalar transport is essential for the numerical modeling of many fluid dynamics problems. However, solutions from numerical schemes always suffer from numerical diffusion and/or oscillation. In this study, we develop an Intra-Cell Advection Tracking (ICAT) scheme to minimize numerical diffusion and preserve monotonicity for advection-dispersion modeling. The key idea is to introduce “queues” in each discretized cell, and using a sequential transport rule and a flow distribution mechanism to track the scalar transport in these queues temporally and spatially. The capability and limitations of ICAT are first investigated through three test cases. Compared with the results obtained from other numerical schemes, the results from ICAT show substantially reduced numerical diffusion and agree better with analytical solutions. We also employ ICAT to simulate the transport process of a conservative tracer in a fracture with a highly heterogeneous aperture distribution. We report discrete flow channels in the fracture are better discerned by ICAT than by other numerical schemes, indicating the suitability of ICAT for modeling tracer transport in channelized flow fields.

15 GEOTHERMAL ENERGY↗
NASA NTRS1978

Numerical simulation of warm fog and its application to warm fog prediction and modification

The considered theoretical model describes the evolution of potential temperature, water vapor content, liquid water content, and horizontal and vertical winds as determined by the processes of vertical turbulent transfer and horizontal advection for momentum, energy, and moisture, as well as radiation cooling, growth of water droplets based on microphysical processes, and drop sedimentation. The mathematical model is two-dimensional in the X-Z plane. The diffusivity coefficient is the same for liquid water droplets as for vapor. The fundamental equations governing the macrophysical processes of the evolution of wind components, water vapor content, liquid water content, and potential temperature under the influences of vertical turbulent diffusion transfer, turbulent momentum transfer, and turbulent energy transfer are expressed by three sets of conservation equations.

Hung, R. J.↗
DOE OSTI2023

Thermo4PFM: Facilitating Phase-field simulations of alloys with thermodynamic driving forces

Phase-field modeling is a popular front-tracking approach used to model solidification. Its time-evolution equations are often coupled to alloy composition and/or thermal diffusion in high-resolution multiphysics approaches. Materials thermodynamic properties tabulated in CALPHAD databases can be used for phase-field modeling to parameterize bulk energies of alloys. In addition, they can be naturally integrated into models such as the Kim-Kim-Suzuki (KKS) model where driving forces depend on the differences between chemical potentials of co-existing phases. In that case, a small system of coupled nonlinear equations needs to be solved at every point in space where the phase-field order parameter is to be updated and evolved in time. Here we present Thermo4PFM, a solver for the KKS equations for binary and ternary alloys, with two or three phases, and parameterized with CALPHAD models. Thermo4PFM is open source, written in C++, and can take advantage of Graphics Processing Units (GPU) accelerators. Using OpenMP offload capabilities for C++ classes, an excellent performance is demonstrated on GPU using the LLVM compiler. CALPHAD data is read from simple JSON files using an open source parser from the boost library.

36 MATERIALS SCIENCE↗
DOE OSTI2015

The hysteresis response of soil CO 2 concentration and soil respiration to soil temperature

Diurnal hysteresis between soil temperature (T s ) and both CO 2 concentration ([CO 2 ]) and soil respiration rate (R s ) were reported across different field experiments. However, the causes of these hysteresis patterns remain a subject of debate, with biotic and abiotic factors both invoked as explanations. Here, to address these issues, a CO 2 gas transport model is developed by combining a layer-wise mass conservation equation for subsurface gas phase CO 2 , Fickian diffusion for gas transfer, and a CO 2 source term that depends on soil temperature, moisture, and photosynthetic rate. Using this model, a hierarchy of numerical experiments were employed to disentangle the causes of the hysteretic [CO 2 ]-T s and CO 2 flux T s (i.e., F-T s ) relations. Model results show that gas transport alone can introduce both [CO 2 ]-T s and F-T s hystereses and also confirm prior findings that heat flow in soils lead to [CO 2 ] and F being out of phase with T s , thereby providing another reason for the occurrence of both hystereses. The area (A hys ) of the [CO 2 ]-T s hysteresis near the surface increases, while the A hys of the Rs-Ts hysteresis decreases as soils become wetter. Moreover, a time-lagged carbon input from photosynthesis deformed the [CO 2 ]-T s and R s -T s patterns, causing a change in the loop direction from counterclockwise to clockwise with decreasing time lag. An asymmetric 8-shaped pattern emerged as the transition state between the two loop directions. Lastly, tracing the pattern and direction of the hysteretic [CO 2 ]-T s and R s -T s relations can provide new ways to fingerprint the effects of photosynthesis stimulation on soil microbial activity and detect time lags between rhizospheric respiration and photosynthesis.

58 GEOSCIENCES↗
DOE OSTI2016

Modeling the Capacitive Deionization Process in Dual-Porosity Electrodes

In many areas of the world, there is a need to increase water availability. Capacitive deionization (CDI) is an electrochemical water treatment process that can be a viable alternative for treating water and for saving energy. A model is presented to simulate the CDI process in heterogeneous porous media comprising two different pore sizes. It is based on a theory for capacitive charging by ideally polarizable porous electrodes without Faradaic reactions or specific adsorption of ions. A two steps volume averaging technique is used to derive the averaged transport equations in the limit of thin electrical double layers. A one-equation model based on the principle of local equilibrium is derived. The constraints determining the range of application of the one-equation model are presented. The effective transport parameters for isotropic porous media are calculated solving the corresponding closure problems. The source terms that appear in the average equations are calculated using theoretical derivations. The global diffusivity is calculated by solving the closure problem.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Fluctuating Hydrodynamics and Debye-Hückel-Onsager Theory for Electrolytes

We apply fluctuating hydrodynamics to strong electrolyte mixtures to compute the concentration corrections for chemical potential, diffusivity, and conductivity. We show these corrections to be in agreement with the limiting laws of Debye, Hückel, and Onsager. We compute explicit corrections for a symmetric ternary mixture and find that the co-ion Maxwell-Stefan diffusion coefficients can be negative, in agreement with experimental findings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗