SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,693 records · Page 94

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Probability and Cumulative Density Function Methods for the Stochastic Advection-Reaction Equation

We present a cumulative density function (CDF) method for the probabilistic analysis of $d$-dimensional advection-dominated reactive transport in heterogeneous media. We employ a probabilistic approach in which epistemic uncertainty on the spatial heterogeneity of Darcy-scale transport coefficients is modeled in terms of random fields with given correlation structures. Our proposed CDF method employs a modified Large-Eddy-Diffusivity (LED) approach to close and localize the nonlocal equations governing the one-point PDF and CDF of the concentration field, resulting in a $(d + 1)$ dimensional PDE. Compared to the classsical LED localization, the proposed modified LED localization explicitly accounts for the mean-field advective dynamics over the phase space of the PDF and CDF. To illustrate the accuracy of the proposed closure, we apply our CDF method to one-dimensional single-species reactive transport with uncertain, heterogeneous advection velocities and reaction rates modeled as random fields.

advection-reaction equation↗

DOE OSTI2013

Multi-moment advection scheme in three dimension for Vlasov simulations of magnetized plasma

We present an extension of the multi-moment advection scheme [T. Minoshima, Y. Matsumoto, T. Amano, Multi-moment advection scheme for Vlasov simulations, Journal of Computational Physics 230 (2011) 6800–6823] to the three-dimensional case, for full electromagnetic Vlasov simulations of magnetized plasma. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. Similar to the two-dimensional scheme, the three-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile with little numerical dispersion or diffusion. This is a very important property for Vlasov simulations of magnetized plasma. We apply the scheme to electromagnetic Vlasov simulations. Propagation of linear waves and nonlinear evolution of the electron temperature anisotropy instability are successfully simulated with a good accuracy of the energy conservation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

On the surface photovoltaic effect in a multivalley semiconductor in an external magnetic field

A theory of the photovoltaic effect in a semi-infinite multvalley semiconductor upon the absorption of polarized light at free carriers caused by the specular reflection and diffuse scattering of electrons at a film surface is developed by us. The kinetic Boltzmann equation in the approximation of time relaxation and boundary conditions, which determine the correlation between the distribution function of electrons reflected from a semi-infinite crystal surface and the distribution function of electrons incident onto a surface for the case of both specular reflection and diffuse scattering is used. The aforementioned takes into account the fact that the distribution function of electrons diffusely scattered from a surface depends only on their energy and is determined from the condition of total electron flux vanishing at the surface. Expressions for analyzing the spectral dependence of the current which is linearly dependent on the magnetic-field strength are obtained.

36 MATERIALS SCIENCE↗

DOE OSTI2013

Ion drag force on a dust grain in a weakly ionized collisional plasma

The problem of calculating the ion drag force acting on a dust grain immersed in a weakly ionized collisional plasma is studied using an approach based on the direct numerical solution of the Vlasov-Bhatnagar-Gross-Krook kinetic equations. A uniform subthermal flow of argon plasma past a spherical dust grain is considered. The numerical computations are performed for a wide range of plasma pressures. On the basis of the obtained results, the effect of ion-neutral collisions on the ion drag force is analyzed in a wide range of ion collisionality. In the collisionless limit, our results are shown to be in good agreement with the results obtained by the binary collision approach. As the ion collisionality increases, the ion drag force is found to decrease sharply and even become negative, i.e., directed oppositely to the plasma flow. A qualitative explanation of this effect is presented and a comparison of our results with those obtained using the drift diffusion approach is discussed. The velocity dependence of the ion drag force in the highly collisional regime is examined. The relationship between the ion and the neutral drag forces in the highly collisional limit is analyzed and the possibility of a superfluid-like behavior of dust grains is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Langevin and Fokker-Planck analyses for diffusion-mediated passing of circular and discorectangular species in two-dimensional channels

The propensity for pairs of diffusing species to pass each other within narrow channels or pores is of basic interest as a first-passage-type problem. It is also of relevance for solution-phase transport in nanoporous materials, and in particular for catalytic conversion reactions where high yield requires that product species can efficiently pass reactant species to exit the pores. Here, we analyze a two-dimensional model with nonoverlapping circular and discorectangular species confined to a rectangular channel, and where passing is mediated by Brownian dynamics in an implicit solvent. For narrower channels where passing is still possible, the discorectangle must align with the channel to pass the circular species. Behavior of the passing propensity, 𝑃, can be assessed by strongly damped Langevin simulations, or within an equivalent Fokker-Planck equation (FPE) formalism. The latter corresponds to a diffusion problem in a “higher-dimensional channel” with a constriction. We assess the variation of the passing propensity, 𝑃, for a broad range of channel width including its scaling just above the threshold where passing is sterically blocked. Analysis of 𝑃 versus the rotational diffusion coefficient 𝐷 𝑟 of the discorectangle reveals a significant decrease in 𝑃 for lower 𝐷 𝑟 for moderate channel width. This prompts a direct analysis of the regime where 𝐷 𝑟 → 0 , for which the FPE can be reduced to a three-dimensional diffusion problem, precise analysis of which is facilitated by adaptive-mesh finite element methods. The dependence of 𝑃 on the aspect ratio of the discorectangle is also assessed.

Rahman, Md Khaledur [Ames Laboratory (AMES), Ames,↗

DOE OSTI2014

Temperature dependent electrical transport studies of self-aligned ZnO nanorods/Si heterostructures deposited by sputtering

Self-aligned ZnO nanorods (NRs) were grown on n-Si(100) substrate by RF sputtering techniques. The NRs are uniformly grown on 2-inch wafer along [0001] direction. Single-crystalline wurtzite structure of ZnO NRs was confirmed by X-ray diffraction. The average diameter, height, and density of NRs are found 48 nm, 750 nm, and 1.26 × 10{sup 10} cm{sup −2}, respectively. The current-voltages (I-V) characteristics of ZnO NRs/Si heterojunction (HJ) were studied in the temperature range of 120–300 K and it shows a rectifying behavior. Barrier height (ϕ{sub B}) and ideality factor (η) were estimated from thermionic emission model and found to be highly temperature dependent in nature. Richardson constant (A{sup *}) was evaluated using Richardson plot of ln(I{sub o}/T{sup 2}) versus q/kT plot by linear fitting in two temperature range 120–180 K and 210–300 K. Large deviation in Richardson constant from its theoretical value of n-Si indicates the presence of barrier inhomogeneities at HJ. Double Gaussian distribution of barrier height with thermionic equation gives mean barrier heights of 0.55 ± 0.01 eV and 0.86 ± 0.02 eV for two different temperature regions 120–180 K and 210–300 K, respectively. Modified Richardson plot provided two values of Richardson constant for two temperature regions. However, for higher temperature range (210–300 K), the calculated value of Richardson constant ∼123 A cm{sup −2} K{sup −2} was close to the ideal Richardson constant for n-Si.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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

Tadmor, E.↗

DOE OSTI2013

Photoacoustic elastic bending in thin film—Substrate system

Theoretical model for optically excited two-layer elastic plate, which includes plasmaelastic, thermoelastic, and thermodiffusion mechanisms, is given in order to study the dependence of the photoacoustic (PA) elastic bending signal on the optical, thermal, and elastic properties of thin film—substrate system. Thin film-semiconductor sample (in our case Silicon) is modeled by simultaneous analysis of the plasma, thermal, and elastic wave equations. Multireflection effects in thin film are included in theoretical model and analyzed. Relations for the amplitude and phase of electronic and thermal elastic bending in the optically excited two-layer mechanically-supported circular plate are derived. Theoretical analysis of the thermodiffusion, plasmaelastic, and thermoelastic effects in a sample-gas-microphone photoacoustic detection configuration is given. Two normalization procedures of the photoacoustic elastic bending signal in function of the modulation frequency of the optical excitation are established. Given theoretical model can be used for various photoacoustic detection configurations, for example, in the study of optical, thermal, and elastic properties of the dielectric-semiconductor or metal-semiconductor structure, etc., Theoretical analysis shows that it is possible to develop new noncontact and nondestructive experimental method—PA elastic bending method for thin film study, with possibility to obtain the optical, thermal, and elastic parameters of the film thinner than 1 μm.

36 MATERIALS SCIENCE↗

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 OSTI2014

Quantum molecular dynamics simulation of shock-wave experiments in aluminum

We present quantum molecular dynamics calculations of principal, porous, and double shock Hugoniots, release isentropes, and sound velocity behind the shock front for aluminum. A comprehensive analysis of available shock-wave data is performed; the agreement and discrepancies of simulation results with measurements are discussed. Special attention is paid to the melting region of aluminum along the principal Hugoniot; the boundaries of the melting zone are estimated using the self-diffusion coefficient. Also, we make a comparison with a high-quality multiphase equation of state for aluminum. Independent semiempirical and first-principle models are very close to each other in caloric variables (pressure, density, particle velocity, etc.) but the equation of state gives higher temperature on the principal Hugoniot and release isentropes than ab initio calculations. Thus, the quantum molecular dynamics method can be used for calibration of semiempirical equations of state in case of lack of experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗