Modeling fast diffusion processes in time integration of stiff stochastic differential equations.
Abstract not provided.
Han, Xiaoying↗
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.
Role of charge separation and pressure diffusion in the gascap of entry objects
Diffusive transport term in hydrodynamic equations describing flow over strongly ablating entry objects, discussing charge separation and pressure diffusion
Howe, J. T.↗
Diffusion and viscosity coefficients for helium
The first order Boltzmann-Fokker-Planck equation is solved numerically to obtain diffusion and viscosity coefficients for a ternary gas mixture composed of electron, protons, and helium. The coefficients are tabulated for five He/H abundances ranging from 0.01 to 10 and for both He II and He III. Comparison with Burgers's thermal diffusion coefficients reveals a maximum difference of 9-10% for both He II and He III throughout the range of helium abundances considered. The viscosity coefficients are compared to those of Chapman and Cowling and show a maximum difference of only 5-6% for He II but 15-16% for He III. For the astrophysically important gas mixtures, it is concluded that the results of existing studies which employed Burgers's or Chapman and Cowling's coefficients will remain substantially unaltered.
Roussel-Dupre, R.↗
A Finite-Element Approach for the Partitioning of Carbon in Q&P Steel
This paper reports a Galerkin finite-element analysis of carbon partitioning from martensite into austenite during the quenching and partitioning (Q&P) processes in steels. In contrast to classical or sophisticated diffusion field models, an alternative nonlinear governing equation based on chemical potential and composition is considered. The model is applied to simulate the carbon partitioning of modified 22MnB5 alloys assuming an immobile austenite–martensite phase boundary and, in turn, is compared with experimental measurements of the volume fraction and carbon content of retained austenite. The simulations show outstanding results on the influence of Si and C content on the partitioning parameters.
36 MATERIALS SCIENCE↗
Enhancement of the droplet nucleation in a dense supersaturated Lennard-Jones vapor
The vapor–liquid nucleation in a dense Lennard-Jones system is studied analytically and numerically. A solution of the nucleation kinetic equations, which includes the elementary processes of condensation/evaporation involving the lightest clusters, is obtained, and the nucleation rate is calculated. Based on the equation of state for the cluster vapor, the pre-exponential factor is obtained. The latter diverges as a spinodal is reached, which results in the nucleation enhancement. The work of critical cluster formation is calculated using the previously developed two-parameter model (TPM) of small clusters. A simple expression for the nucleation rate is deduced and it is shown that the work of cluster formation is reduced for a dense vapor. This results in the nucleation enhancement as well. To verify the TPM, a simulation is performed that mimics a steady-state nucleation experiments in the thermal diffusion cloud chamber. The nucleating vapor with and without a carrier gas is simulated using two different thermostats for the monomers and clusters. The TPM proves to match the simulation results of this work and of other studies.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Observation of Weyl exceptional rings in thermal diffusion
A non-Hermitian Weyl equation indispensably requires a three-dimensional (3D) real/synthetic space, and it is thereby perceived that a Weyl exceptional ring (WER) will not be present in thermal diffusion given its purely dissipative nature. Here, we report a recipe for establishing a 3D parameter space to imitate thermal spinor field. Two orthogonal pairs of spatiotemporally modulated advections are employed to serve as two synthetic parameter dimensions, in addition to the inherent dimension corresponding to heat exchanges. We first predict the existence of WER in our hybrid conduction–advection system and experimentally observe the WER thermal signatures verifying our theoretical prediction. When coupling two WERs of opposite topological charges, the system further exhibits surface-like and bulk topological states, manifested as stationary and continuously changing thermal processes, respectively, with good robustness. Our findings reveal the long-ignored topological nature in thermal diffusion and may empower distinct paradigms for general diffusion and dissipation controls.
42 ENGINEERING↗
A hybridized discontinuous Galerkin method on mapped deforming domains
Here in this paper we present a hybridized discontinuous Galerkin (HDG) discretization for unsteady simulations of convection-dominated flows on mapped deforming domains. Mesh deformation is achieved through an arbitrary Lagrangian-Eulerian transformation with an analytical mapping. We present details of this transformation applied to the HDG system of equations, with focus on the auxiliary gradient equation, viscous stabilization, and output calculation. We discuss conditions under which optimal unsteady output convergence rates can be attained, and we show that both HDG and discontinuous Galerkin (DG) achieve these rates in advection-dominated flows. Results for scalar advection-diffusion and the Euler equations verify the implementation of the mesh motion mapping for both discretizations, show that HDG and DG yield similar results on a given mesh, and demonstrate differences in output convergence rates depending on the choice of HDG viscous stabilization. We note that such similar results bode well for HDG, which has fewer globally-coupled degrees of freedom compared to DG. A simulation of the unsteady compressible Navier-Stokes equations demonstrates again very similar results for HDG and DG and illustrates a pitfall of using steady-state adapted meshes for accurate unsteady simulations.
42 ENGINEERING↗
Lattice thermal conductivity of crystalline and amorphous silicon with and without isotopic effects from the ballistic to diffusive thermal transport regime
Thermal conductivity of a material is an important physical parameter in electronic and thermal devices, and as the device size shrinks down, its length-dependence becomes unable to be neglected. Even in micrometer scale devices, materials having a long mean free path of phonons, such as crystalline silicon (Si), exhibit a strong length dependence of the thermal conductivities that spans from the ballistic to diffusive thermal transport regime. In this work, through non-equilibrium molecular-dynamics (NEMD) simulations up to 17 μm in length, the lattice thermal conductivities are explicitly calculated for crystalline Si and up to 2 μm for amorphous Si. The Boltzmann transport equation (BTE) is solved within a frequency-dependent relaxation time approximation, and the calculated lattice thermal conductivities in the BTE are found to be in good agreement with the values obtained in the NEMD. The isotopic effects on the length-dependent lattice thermal conductivities are also investigated both in the crystalline and amorphous Si.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Pitch-angle scattering of energetic charged particles in nearly constant magnitude magnetic turbulence
We use a method developed by Roberts that optimizes the phase angles of an ensemble of plane waves with amplitudes determined from a Kolmogorov-like power spectrum, to construct magnetic field vector fluctuations having nearly constant magnitude and large variances in its components. This is a representation of the turbulent magnetic field consistent with that observed in the solar wind. Charged-particle pitch-angle diffusion coefficients are determined by integrating the equations of motion for a large number of charged particles moving under the influence of forces from our predefined magnetic field. We tested different cases by varying the kinetic energy of the particles (E {sub p}) and the turbulent magnetic field variance (σ{sub B}{sup 2}). For each combination of E {sub p} and σ{sub B}{sup 2}, we tested three different models: (1) the so-called “slab” model, where the turbulent magnetic field depends on only one spatial coordinate and has significant fluctuations in its magnitude (b=√(δB{sub x}{sup 2}(z)+δB{sub y}{sup 2}(z)+B{sub 0}{sup 2})); (2) the slab model optimized with nearly constant magnitude b; and (3) the slab model turbulent magnetic field with nearly constant magnitude plus a “variance-conserving” adjustment. In the last case, this model attempts to conserve the variance of the turbulent components (σ{sub Bx}{sup 2}+σ{sub By}{sup 2}), which is found to decrease during the optimization with nearly constant magnitude. We found that there is little or no effect on the pitch-angle diffusion coefficient D{sub μμ} between models 1 and 2. However, the result from model 3 is significantly different. We also introduce a new method to accurately determine the pitch-angle diffusion coefficients as a function of μ.
79 ASTRONOMY AND ASTROPHYSICS↗
Cosmic-ray transport in accelerating flows
The quasi-linear transport equation of energetic charged particles is derived, including scattering by Alfven waves propagating parallel and antiparallel to a uniform magnetic field, losses by synchrotron radiation, and acceleration of the cold background fluid supporting the waves. As in the comoving frame equations of radiative transfer, scattering, and loss terms are evaluated in the frame in which the fluid is locally at rest. The diffusion approximation is applied to the resulting equation, yielding a transport equation for the isotropic part of the distribution function and a first-order approximation to the anisotropy. The conditions under which this approximation may be applied are discussed. In the nonrelativistic regime, the standard equation of the diffusion approximation is recovered. As contributions previously investigated, new effects appear from the spatial and temporal variations of the fluid velocity as well as the synchrotron radiation terms. These are discussed in detail.
Kirk, John G.↗
PARTICLE SCATTERING OFF OF RIGHT-HANDED DISPERSIVE WAVES
Resonant scattering of fast particles off low frequency plasma waves is a major process determining transport characteristics of energetic particles in the heliosphere and contributing to their acceleration. Usually, only Alfvén waves are considered for this process, although dispersive waves are also present throughout the heliosphere. We investigate resonant interaction of energetic electrons with dispersive, right-handed waves. For the interaction of particles and a single wave a variable transformation into the rest frame of the wave can be performed. Here, well-established analytic models derived in the framework of magnetostatic quasi-linear theory can be used as a reference to validate simulation results. However, this approach fails as soon as several dispersive waves are involved. Based on analytic solutions modeling the scattering amplitude in the magnetostatic limit, we present an approach to modify these equations for use in the plasma frame. Thereby we aim at a description of particle scattering in the presence of several waves. A particle-in-cell code is employed to study wave–particle scattering on a micro-physically correct level and to test the modified model equations. We investigate the interactions of electrons at different energies (from 1 keV to 1 MeV) and right-handed waves with various amplitudes. Differences between model and simulation arise in the case of high amplitudes or several waves. Analyzing the trajectories of single particles we find no microscopic diffusion in the case of a single plasma wave, although a broadening of the particle distribution can be observed.
79 ASTRONOMY AND ASTROPHYSICS↗
Proton diffusion and hydrogen/deuterium exchange in amorphous solid water at temperatures from 114 to 134 K
The reaction coefficient for hydrogen/deuterium (H/D) exchange and the diffusion of hydrated excess protons within amorphous solid water (ASW) are characterized as a function of temperature. For these experiments, water films are deposited on a Pt(111) substrate at 108 K, and reactions with pre-adsorbed hydrogen atoms produce hydrated protons. Upon heating, protons diffuse within the water, and H/D exchange occurs when they encounter D2O probe molecules deposited in the films. The time-dependent concentration of D2O is monitored with infrared spectroscopy, and it indicates the protons diffusion from the substrate and establish an equilibrium distribution prior to significant H/D exchange for temperatures 114 K ≤T≤ 134 K. By controlling the distance between the D2O molecules and the substrate, we probe the distribution of protons within the film. It decays as x−2 for the examined range of x (12–52 nm) due to the electric field that develops between the diffusing protons and their image charges in the metal substrate. This agrees with the theoretical distance scaling for the equilibrated proton concentration in a dielectric near a metal boundary. From the proton concentration and the measured D2O decay rate, a lower bound for the proton diffusion coefficient ranging from 10−20 m2/s at 114 K to 10−18 m2/s at 134 K is estimated. The diffusion coefficient has an activation energy of 0.40 eV, which is comparable to energies reported for molecular translations and rotations of H2O, suggesting they may play a critical role in the proton diffusion mechanism within ASW.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The lunar thermal ice pump
It has long been suggested that water ice can exist in extremely cold regions near the lunar poles, where sublimation loss is negligible. The geographic distribution of H-bearing regolith shows only a partial or ambiguous correlation with permanently shadowed areas, thus suggesting that another mechanism may contribute to locally enhancing water concentrations. We show that under suitable conditions, water molecules can be pumped down into the regolith by day-night temperature cycles, leading to an enrichment of H{sub 2}O in excess of the surface concentration. Ideal conditions for pumping are estimated and found to occur where the mean surface temperature is below 105 K and the peak surface temperature is above 120 K. These conditions complement those of the classical cold traps that are roughly defined by peak temperatures lower than 120 K. On the present-day Moon, an estimated 0.8% of the global surface area experiences such temperature variations. Typically, pumping occurs on pole-facing slopes in small areas, but within a few degrees of each pole the equator-facing slopes are preferred. Although pumping of water molecules is expected over cumulatively large areas, the absolute yield of this pump is low; at best, a few percent of the H{sub 2}O delivered to the surface could have accumulated in the near-surface layer in this way. The amount of ice increases with vapor diffusivity and is thus higher in the regolith with large pore spaces.
79 ASTRONOMY AND ASTROPHYSICS↗
Reaction-diffusive dynamics of number-conserving dissipative quantum state preparation
The use of dissipation for the controlled creation of nontrivial quantum many-body correlated states is of much fundamental and practical interest. What is the result of imposing number conservation, which, in closed system, gives rise to diffusive spreading? For this work, we investigate this question for a paradigmatic model of a two-band system, with dissipative dynamics aiming to empty one band and to populate the other, which had been introduced before for the dissipative stabilization of topological states. Going beyond the mean-field treatment of the dissipative dynamics, we demonstrate the emergence of a diffusive regime for the particle and hole density modes at intermediate length- and timescales, which, interestingly, can only be excited in nonlinear response to external fields. We also identify processes that limit the diffusive behavior of this mode at the longest length- and timescales. Strikingly, we find that these processes lead to a reaction-diffusion dynamics governed by the Fisher-Kolmogorov-Petrovsky-Piskunov equation, making the designed dark state unstable towards a state with a finite particle and hole density.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Nonstationary radiation field
Radiation diffusion without redistribution by frequencies, equations for radiative transfer, and stellar spectrum
DIFFUSION WAVE↗
An analysis of the flow field near the fuel injection location in a gas core reactor.
An analytical study is presented which shows the effects of large energy release and the concurrent high acceleration of inner stream fluid on the coaxial flow field in a gas core reactor. The governing equations include the assumptions of only radial radiative transport of energy represented as an energy diffusion term in the Euler equations. The method of integral relations is used to obtain the numerical solution. Results show that the rapidly accelerating, heat generating inner stream actually shrinks in radius as it expands axially.
Weinstein, H.↗
Development of monitoring system of helium leakage from canister
This paper presents a computational method for the helium leakage from a canister. The governing equations for compressible fluids consist of mass conservation equation in Eulerian description, momentum equations and energy equation. The numerical procedures are divided into three phases, advection, diffusion and acoustic phases, and the equations of compressible fluids are discretized with a finite volume method. Thus, the mass conservation law is sufficiently satisfied in the calculation region. In particular, our computational method enables us to predict the change of the temperature distributions around the canister boundaries by calculating the governing equations for the compressible gas flows, which are leaked out from a slight crack on the canister boundary. In order to confirm the validity of our method, it was applied to the basic problem, 2-dimensional natural convection flows in a rectangular cavity. As a result, it was shown that the naturally convected flows can be reasonably simulated by our method. Furthermore, numerical experiments were conducted for the helium leakage from canister and we derived a close relationship between the inner pressure and the boundary temperature distributions.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Numerical solution of the radiation transport equation in disk geometry
An efficient numerical method for solving the problem of radiation transport in a dusty medium with two dimensional (2-D) disk geometry is described. It is a generalization of the one-dimensional quasi-diffusion method in which the transport equation is cast in diffusion form and then solved as a boundary value problem. The method should be applicable to a variety of astronomical sources, the dynamics of which are angular-momentum dominated and hence not accurately treated by spherical geometry, e.g., protoplanetary nebulae, circumstellar disks, interstellar molecular clouds, accretion disks, and disk galaxies. The computational procedure and practical considerations for implementing the method are described in detail. To illustrate the effects of 2-D radiation transport, some model results (dust temperature distributions and IR flux spectra) for externally heated, interstellar dust clouds with spherically symmetric and disk geometry are compared.
Spagna, George F., Jr.↗