Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A flexible nonlinear diffusion acceleration method for the S N transport equations discretized with discontinuous finite elements

This paper presents a flexible nonlinear diffusion acceleration (NDA) method that discretizes both the S N transport equation and the diffusion equation using the discontinuous finite element method (DFEM). The method is flexible in that the diffusion equation can be discretized on a coarser mesh with the only restriction that it is nested within the transport mesh and the FEM shape function orders of the two equations can be different. The consistency of the transport and diffusion solutions at convergence is defined by using a projection operator mapping the transport into the diffusion FEM space. The diffusion weak form is based on the modified incomplete interior penalty (MIP) diffusion DFEM discretization that is extended by volumetric drift, interior face, and boundary closure terms. In contrast to commonly used coarse mesh finite difference (CMFD) methods, the presented NDA method uses a full FEM discretized diffusion equation for acceleration. Suitable projection and prolongation operators arise naturally from the FEM framework. Via Fourier analysis and numerical experiments for a one-group, fixed source problem the following properties of the NDA method are established for structured quadrilateral meshes: (1) the presented method is unconditionally stable and effective in the presence of mild material heterogeneities if the same mesh and identical shape functions either of the bilinear or biquadratic type are used, (2) the NDA method remains unconditionally stable in the presence of strong heterogeneities, (3) the NDA method with bilinear elements extends the range of effectiveness and stability by a factor of two when compared to CMFD if a coarser diffusion mesh is selected. In addition, the method is tested for solving the C5G7 multigroup, eigenvalue problem using coarse and fine mesh acceleration. Finally, while NDA does not offer an advantage over CMFD for fine mesh acceleration, it reduces the iteration count required for convergence by almost a factor of two in the case of coarse mesh acceleration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Multigrid treatment of implicit continuum diffusion

Implicit treatment of diffusive terms of various differential orders common in continuum mechanics modeling, such as computational fluid dynamics, is investigated with spectral and multigrid algorithms in non-periodic 2D domains. In doubly periodic time dependent problems these terms are handled efficiently by spectral methods, but in non-periodic systems solved with distributed memory parallel computing and 2D domain decomposition, this efficiency is lost for a large number of processors. Here, we built and present here a multigrid algorithm for these types of problems that outperforms a spectral solution employing the highly optimized FFTW library. This solver is suitable for high performance computing and may be able to efficiently treat implicit diffusion of arbitrary order by introducing auxiliary equations of lower order. We test these solvers for fourth and sixth order diffusion with harmonic test functions as well as turbulent 2D magnetohydrodynamic simulations. It is also shown that an anisotropic operator without mixed-derivative terms improves model accuracy and speed, and we examine the impact that the various diffusion operators have on the energy, the enstrophy, and the qualitative aspect of a simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Perona Malik anisotropic diffusion model using Peaceman Rachford scheme on digital radiographic image

In image processing, it is important to remove noise without affecting the image structure as well as preserving all the edges. Perona Malik Anisotropic Diffusion (PMAD) is a PDE-based model which is suitable for image denoising and edge detection problems. In this paper, the Peaceman Rachford scheme is applied on PMAD to remove unwanted noise as the scheme is efficient and unconditionally stable. The capability of the scheme to remove noise is evaluated on several digital radiography weld defect images computed using MATLAB R2009a. Experimental results obtained show that the Peaceman Rachford scheme improves the image quality substantially well based on the Peak Signal to Noise Ratio (PSNR). The Peaceman Rachford scheme used in solving the PMAD model successfully removes unwanted noise in digital radiographic image.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1967

A Study to Analyze the Permeation of High Density Gases and Propellant Vapors Through Single Layer Teflon or Teflon Structure Materials and Laminations

This report contains the results of a fifteen month analytical and experimental study of the leakage rate of the pressurant gases (N2, He) and the propellant vapors (N2O4,N2H4) through bladder structures consisting of two layers of Teflon separated by a metallic foil diffusion barrier containing microscopic or larger holes. Results were obtained for the steady state leakage rate through circular holes and long rectangular openings in the barrier for arbitrary thicknesses of the two Teflon layers. The effect of hole shape and relative hole position on the leakage rate were studied. The transient problem was analyzed and it was shown that steady state calculations are adequate for estimating the leakage rate. A computer program entitled "Diffusion Analyzer Program" was developed to calculate the leakage rate, both transient and steady state. Finally, the analytical results were compared to experimentally determined values of the leakage rate through a model laminated bladder structure. The results of the analysis are in good agreement with experiment. The experimental effort (Part II of the Bladder Permeation Program) measured the solubility, diffusion coefficient and permeability of helium, nitrogen and nitrogen tetroxide vapor through Teflon TFE and FEP membranes. Data were obtained in the temperature range of 25 to 100 C at pressures ranging from near vacuum to about 20 atmospheres. Results of the experimental effort were compared with the limited data previously reported. As a verification to the applicability of results to actual bladder systems, counter diffusion tests were performed with a laminated sample containing aluminum foil with a selected group of holes.

Johnson, R. L.↗
NASA NTRS1989

Similarity solutions for phase-change problems

A modification of Ivantsov's (1947) similarity solutions is proposed which can describe phase-change processes which are limited by diffusion. The method has application to systems that have n-components and possess cross-diffusion and Soret and Dufour effects, along with convection driven by density discontinuities at the two-phase interface. Local thermal equilibrium is assumed at the interface. It is shown that analytic solutions are possible when the material properties are constant.

Canright, D.↗
DOE OSTI2017

High-Dimensional Intrinsic Interpolation Using Gaussian Process Regression and Diffusion Maps

This article considers the challenging task of estimating geologic properties of interest using a suite of proxy measurements. The current work recast this task as a manifold learning problem. In this process, this article introduces a novel regression procedure for intrinsic variables constrained onto a manifold embedded in an ambient space. The procedure is meant to sharpen high-dimensional interpolation by inferring non-linear correlations from the data being interpolated. The proposed approach augments manifold learning procedures with a Gaussian process regression. It first identifies, using diffusion maps, a low-dimensional manifold embedded in an ambient high-dimensional space associated with the data. It relies on the diffusion distance associated with this construction to define a distance function with which the data model is equipped. This distance metric function is then used to compute the correlation structure of a Gaussian process that describes the statistical dependence of quantities of interest in the high-dimensional ambient space. The proposed method is applicable to arbitrarily high-dimensional data sets. Here, it is applied to subsurface characterization using a suite of well log measurements. The predictions obtained in original, principal component, and diffusion space are compared using both qualitative and quantitative metrics. Considerable improvement in the prediction of the geological structural properties is observed with the proposed method.

58 GEOSCIENCES↗
NASA NTRS1988

The crack problem in bonded nonhomogeneous materials

The plane elasticity problem for two bonded half planes containing a crack perpendicular to the interface was considered. The effect of very steep variations in the material properties near the diffusion plane on the singular behavior of the stresses and stress intensity factors were studied. The two materials were thus, assumed to have the shear moduli mu(o) and mu(o) exp (Beta x), x=0 being the diffusion plane. Of particular interest was the examination of the nature of stress singularity near a crack tip terminating at the interface where the shear modulus has a discontinuous derivative. The results show that, unlike the crack problem in piecewise homogeneous materials for which the singularity is of the form r/alpha, 0 less than alpha less than 1, in this problem the stresses have a standard square-root singularity regardless of the location of the crack tip. The nonhomogeneity constant Beta has, however, considerable influence on the stress intensity factors.

Erdogan, Fazil↗
NASA NTRS1991

The crack problem in bonded nonhomogeneous materials

The plane elasticity problem for two bonded half planes containing a crack perpendicular to the interface was considered. The effect of very steep variations in the material properties near the diffusion plane on the singular behavior of the stresses and stress intensity factors were studied. The two materials were thus, assumed to have the shear moduli mu(o) and mu(o) exp (Beta x), x=0 being the diffusion plane. Of particular interest was the examination of the nature of stress singularity near a crack tip termination at the interface where the shear modulus has a discontinuous derivative. The results show that, unlike the crack problem in piecewise homogeneous materials for which the singularity is of the form r/alpha, 0 less than alpha less than 1, in this problem the stresses have a standard square-root singularity regardless of the location of the crack tip. The nonhomogeneity constant Beta has, however, considerable influence on the stress intensity factors.

Erdogan, F.↗
DOE OSTI2025

Geospatial Diffusion for Land Cover Imperviousness Change Forecasting

Land-use and land-cover (LULC) has a significant effect on several Earth system processes. For example, impervious surfaces reduce infiltration and speed water flow, impacting regional hydrology and flood risk. While Earth System models have improved forecasting hydrologic and atmospheric processes at higher resolutions, the ability to forecast LULC change has lagged behind. In this paper, we propose a new paradigm exploiting Generative AI (GenAI) for land cover change forecasting by framing it as a data synthesis problem conditioned on historical and auxiliary data-sources. To demonstrate the feasibility of our methodology, we perform experiments where a diffusion model is trained for decadal forecasting of imperviousness change across the entire United States. We find that our model yields MAE lower than a no-change baseline for resolutions ≥ 0.7 X 0.7km2 on average, demonstrating its ability to capture and project accurate spatiotemporal patterns. Finally, we discuss future research to incorporate Earth's physical properties and enabling scenario simulations via driver variables.

Varshney, Debvrat [ORNL] (ORCID:0000000188981736)↗
NASA NTRS1997

Spherical Diffusion Flames: Structure and Dynamics

The spherical geometry is the most suitable one to study fundamental issues concerning the structure and the dynamics of diffusion flames. From a theoretical point of view, it is the only geometry that permits the existence of a truly one dimensional stationary diffusion flame. A stationary planar diffusion flame with the fuel supplied upstream, at x approaches -infinity say, and the oxidant downstream, at x approaches +infinity is not possible. For a steady diffusion flame to exist, one must have nonzero fluxes of fuel and oxidant towards the flame. However, in the unlimited region behind the planar flame the only bounded solutions to the reaction-free convective-diffusive operator are constants. Hence the oxidant concentration behind the flame remains constant and there is no mechanism to generate the necessary flux towards the flame. A one-dimensional problem can be formulated if the reactants are supplied at finite locations; but the boundary conditions in this case introduce unnecessary complications and do not appropriately model the physical reality. Indeed, a planar diffusion flame can be established in the stagnation-point flow of two opposed jets but the flame in this case is stretched and the flow is essentially two-dimensional. The only stationary one-dimensional diffusion flame in an unlimited environment is therefore the spherical flame.

Matalon, Moshe↗
NASA NTRS1993

Applying and validating the RANS-3D flow-solver for evaluating a subsonic serpentine diffuser geometry

Subsonic inlet ducts for advanced, high-performance aircraft are evolving towards complex three-dimensional shapes for reasons of overall integration and weight. These factors lead to diffuser geometries that may sacrifice inlet performance, unless careful attention to design details and boundary layer management techniques are employed. The ability of viscous computational fluid dynamic (CFD) analysis of such geometries to aid the aircraft configurator in this complex design problem is herein examined. The RANS-3D Reynolds-Averaged Navier-Stokes solver is applied to model the complex flowfield occurring in a representative diffuser geometry and the solutions are compared to experimental results from a static test of the inlet duct. The computational results are shown to compare very favorably with experimental results over a range of mass flow rates, including those involving large amounts of separation in the diffuser. In addition, a novel grid topology is presented, and two turbulence models are evaluated in this study as part of the RANS-3D code.

Fletcher, Michael J.↗
DOE OSTI2014

A cell-local finite difference discretization of the low-order quasidiffusion equations for neutral particle transport on unstructured quadrilateral meshes

We present a quasidiffusion (QD) method for solving neutral particle transport problems in Cartesian XY geometry on unstructured quadrilateral meshes, including local refinement capability. Neutral particle transport problems are central to many applications including nuclear reactor design, radiation safety, astrophysics, medical imaging, radiotherapy, nuclear fuel transport/storage, shielding design, and oil well-logging. The primary development is a new discretization of the low-order QD (LOQD) equations based on cell-local finite differences. The accuracy of the LOQD equations depends on proper calculation of special non-linear QD (Eddington) factors from a transport solution. In order to completely define the new QD method, a proper discretization of the transport problem is also presented. The transport equation is discretized by a conservative method of short characteristics with a novel linear approximation of the scattering source term and monotonic, parabolic representation of the angular flux on incoming faces. Analytic and numerical tests are used to test the accuracy and spatial convergence of the non-linear method. All tests exhibit O(h{sup 2}) convergence of the scalar flux on orthogonal, random, and multi-level meshes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗
DOE OSTI2018

A Phase-Field Lattice-Boltzmann Study on Dendritic Growth of Al-Cu Alloy Under Convection

Effects of convection (forced and natural) on dendritic evolution of the Al-Cu alloy were investigated using a phase-field lattice-Boltzmann approach. The non-linear coupled equations were solved by applying a parallel and adaptive mesh refinement algorithm. Important physical aspects including dendritic fragmentation, splitting, and formation of solute plumes were simulated. Results showed that the dendritic growth patterns under convection exhibited remarkable difference from those without convection. The presence of flow led to variation of solute diffusion and upstream–downstream dendritic growth difference, which further influenced the development of dendritic arms and multi-dendritic competitive growth. When the convection intensity was magnified, the convection-induced anisotropy became dominated, and the growth patterns changed accordingly to accommodate the local thermodynamic variation.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Equilibration in fermionic systems

Highlights: •Local equilibration in finite fermion systems is based on a nonlinear diffusion equation. •The fermionic diffusion equation is solved analytically in closed form. •The exact solutions agree with the numerical results. •Particle-number and energy conservation are shown to be fulfilled. •The model is applied to hadrons at low energies, and relativistic quarks. -- Abstract: The time evolution of a finite fermion system towards local statistical equilibrium is investigated using analytical solutions of a nonlinear partial differential equation that had been derived earlier from the Boltzmann collision term. The solutions of this fermionic diffusion equation are rederived in closed form, evaluated exactly for simplified initial conditions, and applied to hadron systems at low energies in the MeV-range, as well as to quark systems at relativistic energies in the TeV-range where antiparticle production is abundant. Conservation laws for particle number including created antiparticles, and for the energy are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

An asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations

The solutions of radiative transport equations can cover both optical thin and optical thick regimes due to the large variation of photon's mean-free path and its interaction with the material. In the small mean free path limit, the nonlinear time-dependent radiative transfer equations can converge to an equilibrium diffusion equation due to the intensive interaction between radiation and material. In the optical thin limit, the photon free transport mechanism will emerge. In this paper, we are going to develop an accurate and robust asymptotic preserving unified gas kinetic scheme (AP-UGKS) for the gray radiative transfer equations, where the radiation transport equation is coupled with the material thermal energy equation. The current work is based on the UGKS framework for the rarefied gas dynamics [14], and is an extension of a recent work [12] from a one-dimensional linear radiation transport equation to a nonlinear two-dimensional gray radiative system. The newly developed scheme has the asymptotic preserving (AP) property in the optically thick regime in the capturing of diffusive solution without using a cell size being smaller than the photon's mean free path and time step being less than the photon collision time. Besides the diffusion limit, the scheme can capture the exact solution in the optical thin regime as well. The current scheme is a finite volume method. Due to the direct modeling for the time evolution solution of the interface radiative intensity, a smooth transition of the transport physics from optical thin to optical thick can be accurately recovered. Many numerical examples are included to validate the current approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Enabling practical nanoparticle electrodeposition from aqueous nanodroplets

The rapid rise of technology in the modern world has led to an increased demand for energy. Consequently, it is essential to increase the efficiency of current energy-producing systems due to the poor activity of their catalysts. Nanoparticles play a significant role in energy storage and conversion; however, electrodeposition of nanoparticles is difficult to achieve due to surface heterogeneities, nanoparticle diffusion layer overlap, and the inability to electrodeposit multi-metallic nanoparticles with stoichiometric control. These problems can be solved through nanodroplet-mediated electrodeposition, a technique where water nanodroplets are filled with metal salt precursors that form stable nanoparticles when they collide with a negatively-biased electrode. Further, this method has demonstrated control over nanoparticle size and morphology, displaying a wide variety of applications for the generation of materials with excellent catalytic properties. Historically, the cost of nanodroplet-mediated electrodeposition experimentation is prohibitive because practitioners use 0.1 M to 0.5 M tetrabutylammonium perchlorate (TBAP) dissolved in the oil phase (∼10 mL). Such high concentrations of electrolytes have been used to lower ohmic drop and provide ions to maintain charge balance during electrodeposition. Here, we show that supporting electrolyte is not necessary for the oil phase. In fact, one can use a suitable salt (such as lithium perchlorate) in the aqueous phase to achieve nanoparticle electrodeposition. This simple change, grounded in an understanding of ion transfer, drives down the cost per experiment by nearly three orders of magnitude, representing a necessary step forward in enabling practical nanoparticle electrodeposition from water nanodroplets. This approach is a promising procedure for future cost-effective energy conversion systems relying on electrocatalytic nanoparticles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗