Symmetric gaseous diffusion coefficients.
Kinetic theory development of equations for multicomponent gaseous mixtures symmetric diffusion coefficients
Curtiss, C. F.↗
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Migration of the separation point on a deforming cylinder
An iterative scheme of solving the Navier-Stokes equations for unsteady two-dimensional flow about a deforming and translating cylinder is given. In the kth iteration, the convected vorticity appears as the source term in an equation of transient diffusion of vorticity. A novel integral transform is used to reduce this transient vorticity equation to a k-dimensional heat equation. The bounded solution of this equation is obtained with a general method of superposition for problems involving a moving boundary. An equation describing the migration of the separation point on a deforming cylinder in unsteady cross flows is derived from the analytically obtained velocity field. The radial expansion of the cylinder surface is shown to hasten the separation time and to increase the separation angle. The results imply that in the steady flow past a body that is moving forward and sinking at constant rates, the locus of points at which the azimuthal component of skin friction changes sign originates on the leeward ray at a point downstream of the front tip.
Lin, S. P.↗
Reduced-order modeling for nonlocal diffusion problems
In several settings, diffusive behavior is observed to not follow the rate of spread predicted by parabolic partial differential equations (PDEs) such as the heat equation. Such behaviors, often referred to as anomalous diffusion, can be modeled using nonlocal equations for which points at a finite distance apart can interact. An example of such models is provided by fractional derivative equations. Because of the nonlocal interactions, discretized nonlocal systems have less sparsity, often significantly less, compared with corresponding discretized PDE systems. As such, the need for reduced-order surrogates that can be used to cheaply determine approximate solutions is much more acute for nonlocal models compared with that for PDEs. In this paper, we consider the construction, application, and testing of proper orthogonal decomposition (POD) reduced models for an integral equation model for nonlocal diffusion. For certain modeling parameters, the model we consider allows for discontinuous solutions and includes fractional Laplacian kernels as a special case. Preliminary computational results illustrate the potential of using POD to obtain accurate approximations of solutions of nonlocal diffusion equations at much lower costs compared with, for example, standard finite element methods. Copyright © 2016 John Wiley & Sons, Ltd.
97 MATHEMATICS AND COMPUTING↗
Improvement of the second- and third-moment modeling of turbulence: A study of Reynolds-stress closure model
Four parts of the Reynolds-stress closure modeling are reported: (1) improvement of the k and epsilon equaitons; (2) development of the third-moment transport equation; (3) formulation of the diffusion coefficient of the momentum equation by using the algebraic-stress model of turbulence; and (4) the application of the Reynolds-stress model to a heat exchanger problem. It was demonstrated that the third-moment transport model improved the prediction of the triple-velocity products in the recirculating and reattaching flow regions in comparison with the existing algebraic models for the triple-velocity products. Optimum values for empirical coefficients are obtained for the prediction of the backward-facing step flows. A functional expression is derived for the coefficient of the momentum diffusion by employing the algebraic-stress model. The second-moment closure is applied to a heat transfer problem. The computations for the flow in a corrugated-wall channel show that the second-moment closure improves the prediction of the heat transfer rates by 30% over the k - epsilon model.
Amano, R. S.↗
Resonant diffusion in strongly turbulent plasmas.
The effect of turbulent fluctuations on plasma particles is considered, and equations are derived which describe the evolution of macroscopic properties such as temperature and flow speed of the turbulent plasma. Initially, a diffusion equation for a single-particle distribution function averaged over an ensemble of plasmas is derived for an unmagnetized plasma. For the resonant diffusion in strongly turbulent plasmas, an ensemble of three-dimensional plasmas is considered with an approximately homogeneous and stationary distribution of random electromagnetic fluctuations. For each realization, the single-particle distribution function satisfies the Vlasov equation.
Birmingham, T. J.↗
Vapor Transport Within the Thermal Diffusion Cloud Chamber
A review of the equations used to determine the 1-D vapor transport in the thermal diffusion cloud chamber (TDCC) is presented. These equations closely follow those of the classical Stefan tube problem in which there is transport of a volatile species through a noncondensible, carrier gas. In both cases, the very plausible assumption is made that the background gas is stagnant. Unfortunately, this assumption results in a convective flux which is inconsistent with the momentum and continuity equations for both systems. The approximation permits derivation of an analytical solution for the concentration profile in the Stefan tube, but there is no computational advantage in the case of the TDCC. Furthermore, the degree of supersaturation is a sensitive function of the concentration profile in the TD CC and the stagnant background gas approximation can make a dramatic difference in the calculated supersaturation. In this work, the equations typically used with a TDCC are compared with very general transport equations describing the 1-D diffusion of the volatile species. Whereas no pressure dependence is predicted with the typical equations, a strong pressure dependence is present with the more general equations given in this work. The predicted behavior is consistent with observations in diffusion cloud experiments. It appears that the new equations may account for much of the pressure dependence noted in TDCC experiments, but a comparison between the new equations and previously obtained experimental data are needed for verification.
Ferguson, Frank T.↗
A closure model of diffusion transport of the Reynolds-stress equations and its application to a turbulent step flow
A study is made for the development and computations of the separating and reattaching shear flows. The highlight of the study is an attempt to predict the third-moments of turbulent velocity which is responsible for the diffusion transport of the Reynolds stresses. The present computations show that the third-moments obtained by employing the low-Reynolds number model of transport equations improve the prediction of the third-moments.
Amano, R. S.↗
Mixtures of octanol and an ionic liquid: Structure and transport
Ionic liquids (ILs) with long alkyl substituents are amphiphilic, which leads to a bicontinuous liquid structure. The strongly interacting anionic and cationic head groups form a long range charge network, with the hydrocarbon tails forming a nonpolar domain. Such nonpolar domains have been shown to dissolve a variety of neutral organic solvents. In mixtures of ILs with solvents the neutral organic molecules residing in the nonpolar domains experience different environments and friction from the charged cations and anions. Thus, the neutral molecules diffuse much faster than predicted by hydrodynamic scaling using the average viscosity of the mixture. Here, we report studies on the structure and transport properties of mixtures of 1-octanol with the IL trihexyltetradecylphosphonium bis(trifluoromethylsulfonyl)imide (P 6,6,6,14 + /NTf 2 – ). The majority of the atom fraction in the P 6,6,6,14 + cation comprises four hydrocarbon substituents. The unique amphiphilic nature of ILs with the P 6,6,6,14 + cation makes 1-octanol fully miscible with the IL at ambient temperatures. X-ray scattering experiments show that the IL structure persists in the mixtures for 1-octanol mole fractions as large as x oct = 0.90. The self-diffusion coefficients of the three molecular species in the mixtures were measured by NMR experiments. The self-diffusion of the P 6,6,6,14 + cation is well described by the Stokes–Einstein equation, while the diffusivity of the NTf 2 – anion is slightly lower than the hydrodynamic prediction. The measured diffusivities of octanol in these mixtures are 1.3–4 times higher than the hydrodynamic predictions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Modeling of Fluctuating Mass Flux in Variable Density Flows
The approach solves for both Reynolds and Favre averaged quantities and calculates the scalar pdf. Turbulent models used to close the governing equations are formulated to account for complex mixing and variable density effects. In addition, turbulent mass diffusivities are not assumed to be in constant proportion to turbulent momentum diffusivities. The governing equations are solved by a combination of finite-difference technique and Monte-Carlo simulation. Some preliminary results on simple variable density shear flows are presented. The differences between these results and those obtained using conventional models are discussed.
So, R. M. C.↗
Transportation of volatile elements in thermally evolving planetesimals: An important role of metallic iron
Ordinary chondrites are considered to have experienced thermal metamorphism in small bodies. We are interested in behaviors of volatile elements in such a kind of thermally evolving planetesimals. Volatile elements generally have high vapor pressures at high temperature. In porous bodies, with a high gas permeability, volatile elements are transported efficiently over a long range. Behavior of volatile elements transported by permeable gas flow can be handled by an equation whose form is similar to that of the equation of thermal diffusion. We can follow transportation of heats and volatile elements in planetesimals, when parameters in these equations, initial conditions and chemical behavior of volatile elements are given. Recently, we discovered that nitrogen in equilibrated H-chondrites is mainly trapped in taenite (f.c.c. Fe-Ni), probably dissolved in interstitial sites. Fegley suggests that metallic iron cannot trap nitrogen in the solar nebula gas due to its very low nitrogen partial pressure. Approximately 1 bar of nitrogen pressure is required to explain the nitrogen content in taenite. We may expect high nitrogen gas partial pressure (possibly produced by vaporization of nitrogen-bearing solids such as organic materials) at the interior of thermally evolving planetesimals. Kinetic behavior of nitrogen in taenite suggests that it can easily be equilibrated with the ambient nitrogen gas at temperatures of approximately 500 C or higher. We consider that nitrogen is trapped in taenite through a nitrogen redistribution process occurred during the thermal metamorphic event.
Hashizume, K.↗
New explicit equations for the accurate calculation of the growth and evaporation of hydrometeors by the diffusion of water vapor
The traditional explicit growth equation has been widely used to calculate the growth and evaporation of hydrometeors by the diffusion of water vapor. This paper reexamines the assumptions underlying the traditional equation and shows that large errors (10-30 percent in some cases) result if it is used carelessly. More accurate explicit equations are derived by approximating the saturation vapor-density difference as a quadratic rather than a linear function of the temperature difference between the particle and ambient air. These new equations, which reduce the error to less than a few percent, merit inclusion in a broad range of atmospheric models.
Srivastava, R. C.↗
An effective numerical algorithm for intra-granular fission gas release during non-equilibrium trapping and resolution
Both fission gas release and gaseous swelling in nuclear fuel are driven by the transport of fission gas from within the fuel grains to grain boundaries (intra-granular fission gas release). The process involves gas atom diffusion in conjunction with trapping in and resolution from intra-granular bubbles, and is described mathematically by a system of two partial differential equations (PDE). Following Speight, under the assumption of equilibrium between trapping and resolution (quasi-stationary approximation) the system can be reduced to a single diffusion equation with an effective diffusion coefficient. Numerical solutions used in engineering fuel performance calculations invariably rely on this simplification. First, we investigate the validity of the quasi-stationary approximation compared to the solution of the general system of PDEs. Results demonstrate that the approximation is valid under most conditions of practical interest, but is inadequate to model intra-granular fission gas release during rapid transients to relatively high temperatures such as postulated reactivity initiated accidents (RIA). Then, we develop a novel numerical algorithm for the solution of the general problem in time-varying conditions. We verify the PolyPole-2 algorithm against a reference finite difference solution for a large number of randomly generated operation histories including prototypical RIA transients. Results demonstrate that PolyPole-2 captures the solution of the general PDE system with a high accuracy and a low computational cost. As a result, the PolyPole-2 algorithm overcomes the quasi-stationary approximation and the concept of an effective diffusion coefficient for the solution of the intra-granular fission gas release problem in nuclear fuel analysis.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Modified Closures in Monte Carlo Algorithms for Diffusive Binary Stochastic Media Transport Problems
In a stochastic medium, the material properties at a given spatial location are known only statistically. The most common approach to solving particle transport problems involving binary stochastic media (BSM) is to use the atomic mix approximation in which the transport problem is solved using homogenized (volume-averaged) material properties. A common deterministic model developed for solving BSM particle transport problems is the Levermore-Pomraning (LP) model. Zimmerman and Adams proposed a Monte Carlo algorithm (Algorithm A) that solves the LP equations and another Monte Carlo algorithm (Algorithm B) that locally preserves the sampled material realization; we refer to these Monte Carlo algorithms as the LP and LPLRP (LP local realization preserving) algorithms, respectively. One-dimensional (1D) planar geometry benchmark studies have shown that the LPLRP algorithm is often significantly more accurate than the LP algorithm for problems with an incident angular flux as well as for problems with an interior source. The LPLRP algorithm implementation in these benchmark comparisons made explicit use of the one-dimensional nature of the problem. The LP model is derived by assuming an 'upwind' closure in the coupling term relating the two materials in a binary stochastic medium. Su and Pomraning developed a modified form of this closure by considering the small correlation length limit and requiring the modified closure to produce the correct exponential decay for a source-free half-line albedo problem in rod geometry. They concluded that the modified closure is generally not inferior to the LP closure and in some cases is significantly better. Brantley investigated the use of the Su-Pomraning (SP) closure in the Monte Carlo LP algorithm (LP-SP) for the suite of benchmark problems and concluded that 1) the LP-SP algorithm was somewhat more accurate overall than the LP algorithm and somewhat less accurate overall than the LPLRP algorithm for an incident angular flux benchmark suite, and 2) the LP-SP algorithm was generally the least accurate of the algorithms for an interior source benchmark suite. Larsen et al. performed an asymptotic analysis of the transport equation in 1D planar geometry for the situation in which the physical system is 1) a random binary stochastic medium with material macroscopic total cross sections and mean slab width values of O(1) or smaller and 2) optically thick with weak absorption and sources at each spatial point and therefore globally diffusive. The asymptotic analysis demonstrates that, under these assumptions, the transport equation limits to the conventional diffusion equation with atomically-mixed (volume-averaged) material properties. Larsen et al. also demonstrate that, under these same assumptions, the LP equations asymptotically limit to an atomically-mixed diffusion equation with a diffusion coefficient that is too large, leading to an unphysical flattening of the ensemble-averaged scalar flux distribution. Vasques and Yadav recently performed an asymptotic analysis of an adjusted Levermore-Pomraning closure in which the Markovian transition functions are rescaled such that the equations asymptotically limit to the correct atomically-mixed diffusion equation; we refer to this model as LP-VY. In this paper, we demonstrate that the adjusted LP closure proposed by Vasques and Yadav is a special case of the SP closure obtained by assuming that the absorption cross sections in the two materials are equal. We further demonstrate that the Su-Pomraning closure has the correct asymptotic limit for the diffusive physical system and is therefore an appropriate closure for these problems. We also describe how to incorporate these various closures into the Monte Carlo LP and LPLRP particle transport algorithms, and we present numerical results comparing the accuracy of the Monte Carlo LP, LP-VY, LP-SP, LPLRP, and LPLRP-SP algorithms for the set of diffusive benchmark problems. We demonstrated that the SP closure has the correct asymptotic limit for the diffusive physical system under investigation in this paper and that the VY closure is a special case of the SP closure obtained by assuming equal absorption in the two materials. Both of these closures are readily implemented in the Monte Carlo LP and LPLRP algorithms for particle transport in BSM. Through numerical comparisons to a diffusive benchmark suite, we confirmed that 1) the LP and LPLRP algorithms do not limit to the appropriate diffusion solution and 2) the LP-VY, LP-SP, and LPLRP-SP algorithms limit to the appropriate diffusion solution and are of comparable accuracy. The LPLRP algorithm has been shown to be generally more accurate than the LP algorithm for non-diffusive problems. However, the numerical results in this paper indicate that the LPLRP algorithm does not limit to the correct diffusion equation for the asymptotic limit under consideration. The work in this paper demonstrates that the SP algorithm does limit to the correct diffusion equation for the asymptotic limit considered. Future work will investigate the possibility of whether the LPLRP-SP algorithm is accurate for both diffusive and non-diffusive problems.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Nonlocal Convection-Diffusion Problems on Bounded Domains and Finite-Range Jump Processes
In this paper, a nonlocal convection-diffusion model is introduced for the master equation of Markov jump processes in bounded domains. With minimal assumptions on the model parameters, the nonlocal steady and unsteady state master equations are shown to be well-posed in a weak sense. Finally, then the nonlocal operator is shown to be the generator of finite-range nonsymmetric jump processes and, when certain conditions on the model parameters hold, the generators of finite and infinite activity Lévy and Lévy-type jump processes are shown to be special instances of the nonlocal operator.
97 MATHEMATICS AND COMPUTING↗
Accuracy of approximations to the Navier-Stokes equations
The systems of truncated differential equations that have been proposed to reduce the complexity and large computational costs of solutions to the full Navier-Stokes equations are considered. These systems are computationally efficient and capture all the physically relevant behavior. The systems follow a certain hierarchy: (1) the classical boundary-layer equations with specified edge properties (usually the streamwise pressure distribution); (2) the coupled boundary-layer/inviscid equations; (3) the so-called thin-layer equations that discard streamwise diffusion; and (4) the Navier-Stokes equations. Consideration is given to each of these approximations applied to an incompressible, laminar-separating flow at low and moderate Reynolds numbers. It is pointed out that for any flow or region of flow for which viscous-inviscid interaction effects are small, classical boundary-layer equations will provide a satisfactory description of the viscous flow at a fraction of the computational cost of any higher approximations.
Murphy, J. D.↗
The Prediction of Nozzle Performance and Heat Transfer in Hydrogen/Oxygen Rocket Engines with Transpiration Cooling, Film Cooling, and High Area Ratios
An advanced engineering computational model has been developed to aid in the analysis of chemical rocket engines. The complete multispecies, chemically reacting and diffusing Navier-Stokes equations are modelled, including the Soret thermal diffusion and Dufour energy transfer terms. Demonstration cases are presented for a 1030:1 area ratio nozzle, a 25 lbf film-cooled nozzle, and a transpiration-cooled plug-and-spool rocket engine. The results indicate that the thrust coefficient predictions of the 1030:1 nozzle and the film-cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. It is demonstrated that thermal diffusion has a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle and was shown to represent a significant fraction of the diffusion fluxes occurring in the transpiration-cooled rocket engine.
Kacynski, Kenneth J.↗
On Bi-Grid Local Mode Analysis of Solution Techniques for 3-D Euler and Navier-Stokes Equations
A procedure is presented for utilizing a bi-grid stability analysis as a practical tool for predicting multigrid performance in a range of numerical methods for solving Euler and Navier-Stokes equations. Model problems based on the convection, diffusion and Burger's equation are used to illustrate the superiority of the bi-grid analysis as a predictive tool for multigrid performance in comparison to the smoothing factor derived from conventional von Neumann analysis. For the Euler equations, bi-grid analysis is presented for three upwind difference based factorizations, namely Spatial, Eigenvalue and Combination splits, and two central difference based factorizations, namely LU and ADI methods. In the former, both the Steger-Warming and van Leer flux-vector splitting methods are considered. For the Navier-Stokes equations, only the Beam-Warming (ADI) central difference scheme is considered. In each case, estimates of multigrid convergence rates from the bi-grid analysis are compared to smoothing factors obtained from single-grid stability analysis. Effects of grid aspect ratio and flow skewness are examined. Both predictions are compared with practical multigrid convergence rates for 2-D Euler and Navier-Stokes solutions based on the Beam-Warming central scheme.
Ibraheem, S. O.↗
Concentration Dependence of Solution Shear Viscosity and Solute Mass Diffusivity in Crystal Growth from Solutions
The physical properties of a supersaturated binary solution such as its density rho, shear viscosity eta, and solute mass diffusivity D are dependent on the solute concentration c: rho = rho(c), eta = eta(c), and D = D(c). The diffusion boundary layer equations related to crystal growth from solution are derived for the case of natural convection with a solution density, a shear viscosity, and a solute diffusivity that are all depen- dent on solute concentration. The solution of these equations has demonstrated the following. (1) At the vicinity of the saturation concentration c(sub s) the solution shear viscosity eta depends on rho as eta(sub s) = eta(rho(sub s))varies as square root of rho(c(sub s)). This theoretically derived result has been verified in experiments with several aqueous solutions of inorganic and organic salts. (2) The maximum solute mass transfer towards the growing crystal surface can be achieved for values of c where the ratio of d ln(D(c)/dc) to d ln(eta(c)/dc) is a maximum.
Izmailov, Alexander F.↗