Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical diffusion”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Modeling of non-equilibrium partitioning in direct metal deposited copper–iron concentrated alloys

The direct metal deposition (DMD) additive manufacturing process produces high cooling rates within a small melt pool and can lead to high amounts of solute trapping. These high cooling rates limit diffusion and lead to the formation of non-equilibrium phases. In this work, we utilize a numerical model to calculate the degree of solute trapping, defined as non-equilibrium partitioning. A theoretical case with overall composition fixed to 50Cu-50Fe at.%. was performed to observe the influence of increasing solidification rates. We then simulate DMD of equimolar Cu-Fe powder printed on mild steel substrate and the calculated non-equilibrium phase compositions were consistent with experimental observations reported earlier on this alloy composition. For a single deposited track, cooling rates are high enough to yield significant solute trapping. Here, the degree of solute trapping is highest near the free surface and has a gradient that correlates with the cooling rate gradient.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Diffused morphotropic phase boundary in relaxor-PbTiO 3 crystals: High piezoelectricity with improved thermal stability

Solid solution ferroelectrics are the most widely used piezoelectric material for numerous electromechanical applications, including sensors, actuators, and transducers. A milestone in ferroelectric research is the discovery of the morphotropic phase boundary that was first reported in Pb(Zr x Ti 1–x )O 3 , which has been extensively solicited to improve the performance of various solid solution ferroelectrics, including those having the highest piezoelectricity known today. However, due to the inherent correlation between phase transition and thermodynamic imbalance, the efforts of building the phase boundary encounter the challenge that high performance materials are generally accompanied by property instability. In this paper, we report a comprehensive study on the crystalline symmetry and polar nano-regions in relaxor-PbTiO 3 crystals by synchrotron measurements and property characterizations. In contrast to the common belief that the morphotropic phase boundary is a narrow composition region, the morphotropic phase boundary in ternary Pb(In 1/2 Nb 1/2 )O 3 -Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 is demonstrated as an extensive region covering at least 9% of the PbTiO 3 content. Such a diffused morphotropic phase boundary is associated with the intensive interaction of polar nano-regions, leading to high piezoelectricity (>1500 pC/N) with greatly improved thermal stability, where the piezoelectric variation is ~90% over the temperature range of 273–373 K, which is about a factor of 3 lower compared to its binary counterpart Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 . This work sheds light on the fundamental understanding of nanoscale inhomogeneity and macroscopic symmetry in relaxor-PbTiO 3 , and successfully interlinks the structure and properties in complex solid solutions. The existence of a diffused morphotropic phase boundary is also expected to benefit other ferroic systems beyond ferroelectrics, such as ferromagnetic and multiferroic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Mechanistic modeling study on process optimization and precursor utilization with atmospheric spatial atomic layer deposition

Spatial atomic layer deposition (SALD) is a promising technology with the aim of combining the advantages of excellent uniformity and conformity of temporal atomic layer deposition (ALD), and an industrial scalable and continuous process. In this manuscript, an experimental and numerical combined model of atmospheric SALD system is presented. To establish the connection between the process parameters and the growth efficiency, a quantitative model on reactant isolation, throughput, and precursor utilization is performed based on the separation gas flow rate, carrier gas flow rate, and precursor mass fraction. The simulation results based on this model show an inverse relation between the precursor usage and the carrier gas flow rate. With the constant carrier gas flow, the relationship of precursor usage and precursor mass fraction follows monotonic function. The precursor concentration, regardless of gas velocity, is the determinant factor of the minimal residual time. The narrow gap between precursor injecting heads and the substrate surface in general SALD system leads to a low Péclet number. In this situation, the gas diffusion act as a leading role in the precursor transport in the small gap rather than the convection. Fluid kinetics from the numerical model is independent of the specific structure, which is instructive for the SALD geometry design as well as its process optimization.

42 ENGINEERING↗
DOE OSTI2018

Permeability Variations Associated With Fault Reactivation in a Claystone Formation Investigated by Field Experiments and Numerical Simulations

We studied the relation between rupture and changes in permeability within a fault zone intersecting the Opalinus Clay formation at 300 m depth in the Mont Terri Underground Research Laboratory (Switzerland). A series of water injection experiments were performed in a borehole straddle interval set within the damage zone of the main fault. A three-component displacement sensor allowed an estimation of the displacement of a minor fault plane reactivated during a succession of step rate pressure tests. The experiment reveals that the fault hydromechanical (HM) behavior is different from one test to the other with varying pressure levels needed to trigger rupture and different slip behavior under similar pressure conditions. Numerical simulations were performed to better understand the reason for such different behavior and to investigate the relation between rupture nucleation, permeability change, pressure diffusion, and rupture propagation. Here, our main findings are as follows: (i) a rate frictional law and a rate-and-state permeability law can reproduce the first test, but it appears that the rate constitutive parameters must be pressure dependent to reproduce the complex HM behavior observed during the successive injection tests; (ii) almost similar ruptures can create or destroy the fluid diffusion pathways; (iii) a too high or too low diffusivity created by the main rupture prevents secondary rupture events from occurring whereas “intermediate” diffusivity favors the nucleation of a secondary rupture associated with the fluid diffusion. However, because rupture may in certain cases destroy permeability, this succession of ruptures may not necessarily create a continuous hydraulic pathway.

58 GEOSCIENCES↗
DOE OSTI2026

Carbon fiber design improvements based on economic models

Circular fiber geometries are predominant in commercial carbon fiber material systems, but the use of this fiber shape has numerous limitations. Circular geometries have continuous symmetry, which is helpful for various processing considerations, but also have the largest possible maximum diffusion thickness for a given fiber area. This characteristic means that circular carbon fibers always have the highest material processing cost and lowest production throughput compared to any other fiber shape with the same area and tow count. To quantify material cost and other benefits for non-circular carbon fiber geometries, process models for polyacrylonitrile based carbon fiber production are developed in relationship to the carbon fiber shape, size, and tow count. For a given fiber shape, precursor production costs are shown to favor maximizing fiber size while conversion costs are minimized by the smallest fiber size. These competing cost trends result in a numerically optimal fiber size for a given shape and tow count, while both cost components are decreased by increasing tow count. Cost–performance tradeoffs for three lobe fiber geometries are studied by supplementing the cost trends with a numerical failure model to predict compressive strength for discrete shape variants. The shape selection is shown to be more sensitive to variations in cost than compressive strength while suboptimal shape designs can improve manufacturing robustness and achievable fiber volume fractions. Finally, an optimal three lobe carbon fiber is identified that balances the set of considerations while reducing costs and embodied energy and increasing production throughput compared to a commercial carbon fiber.

Carbon fiber↗
DOE OSTI2018

Quantitative soil zymography: Mechanisms, processes of substrate and enzyme diffusion in porous media

Soil membrane zymography enables 2D mapping of enzyme activities on the surface of soil samples. The method is according to diffusion of components of enzymatically-mediated reactions to/from membrane, and, hence, reflects the distribution of enzyme activities at the intact soil surface. Zymography has been already successfully implemented in numerous soil ecology applications. Here we identify two methodological aspects for further improvement and expansion of the method at micro and macro scales: first, accounting for the area of contact between the soil surface and the zymography membranes and, second, accounting for diffusion effects during the zymography procedure. We tested three methods, namely, laser-scanning, staining with a fluorescent product (e.g. MUF: 4-methylumbelliferone), and X-ray computed micro-tomography, for assessing the area of the soil surface in contact with the membranes. We quantified diffusion of MUF, enzymes and substrate between the substrate-saturated membrane and soil as well as diffusion processes during membrane zymography via HP2 software. Diffusion of the substrate from the membrane and of the MUF-product to the membrane was detected, while there were no clear evidence of enzyme diffusion to/in the membrane. According to the model simulations, the enzyme activities detected via 2D zymography probably represent only a small portion, about 20%, of the actual reactions within the soil volume that is in both direct contact and in hydrological contact with zymography membranes. This is a result of omnidirectional diffusion of reaction products. The membrane contact with the soil surface estimated by three methods ranged from 3.4 to 36.5% further signifying that only a fraction of enzymes activity is detectable in a course of 2D soil zymography.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Numerical investigation of the spatiotemporal distribution of chemical species in an atmospheric surface barrier-discharge

Using a one dimensional time dependent convection-reaction-diffusion model, the temporal and spatial distributions of species propagating downstream of an atmospheric pressure air surface barrier discharge was studied. It was found that the distribution of negatively charged species is more spatially spread compared to positive ions species, which is attributed to the diffusion of electrons that cool down and attach to background gas molecules, creating different negative ions downstream of the discharge region. Given the widespread use of such discharges in applications involving the remote microbial decontamination of surfaces and liquids, the transport of plasma generated reactive species away from the discharge region was studied by implementing mechanical convection through the discharge region. It was shown that increased convection causes the spatial distribution of species density to become uniform. It was also found that many species have a lower density close to the surface of the discharge as convection prevents their accumulation. While for some species, such as NO{sub 2}, convection causes a general increase in the density due to a reduced residence time close to the discharge region, where it is rapidly lost through reactions with OH. The impact of the applied power was also investigated, and it was found that the densities of most species, whether charged or neutral, are directly proportional to the applied power.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Generalizable Porous Aromatic Framework‐Included Polymer Membranes for Diffusion‐Enhanced Gas Separations

Industrial separation processes account for 10-15% of global energy consumption. Membrane-based processes are less energy-intensive than traditional gas separation technologies; however, enhanced material separation performance and stability for numerous gas mixtures are needed for widespread industrial adoption. This work presents a generalizable strategy for preparing mixed-matrix gas separation membranes exceeding the performance upper bounds of existing polymer membranes for a wide variety of industrial gases. By incorporating robust porous aromatic framework (PAF) particles into various dense commercial polymer matrices, gas diffusivity and solubility can be enhanced. For diverse gas mixtures (e.g., CO2/N2, O2/N2, He/CH4, H2/N2, and C2H4/C2H6), the resulting composite membranes exhibit enhanced gas permeabilities-by as much as 520%-and largely unchanged selectivities even after 6 years of aging under simulated flue gas conditions. These improvements arise from the ultrahigh porosity, excellent chemical compatibility, and unique physicochemical properties of the embedded PAF particles. Functionalizing the PAFs with polyamines also enables composite membranes that achieve among the highest reported performances against plasticization, a common obstacle in commercializing gas separation membranes. Significantly, the PAF-1 particles are readily dispersible in various common membrane casting solvents, suggesting their broader utility as a filler for designing high-performance membranes for many industrial gas separations.

Uliana, Adam A↗
DOE OSTI2016

Modeling and analysis of fission product transport in the AGR-3/4 experiment - 18693

In this work we describe the ongoing modeling and analysis efforts in support of the AGR-3/4 experiment. AGR-3/4 is intended to provide data to assess fission product retention and transport (e.g., diffusion coefficients) in fuel matrix and graphite materials. We describe a set of pre-test predictions that incorporate the results of detailed thermal and fission product release models into a coupled 1D radial diffusion model of the experiment, using diffusion coefficients reported in the literature for Ag, Cs, and Sr. We make some comparisons of the predicted Cs profiles to preliminary measured data and find these to be reasonable, in most cases within an order of magnitude. Our ultimate objective is to refine the diffusion coefficients using AGR-3/4 data, so we develop an analytical method for doing so and demonstrate its efficacy via a series of numerical experiments using the model predictions. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

The equilibrium-diffusion limit for radiation hydrodynamics

The equilibrium-diffusion approximation (EDA) is used to describe certain radiation-hydrodynamic (RH) environments. When this is done the RH equations reduce to a simplified set of equations. The EDA can be derived by asymptotically analyzing the full set of RH equations in the equilibrium-diffusion limit. Here, we derive the EDA this way and show that it and the associated set of simplified equations are both first-order accurate with transport corrections occurring at second order. Having established the EDA’s first-order accuracy we then analyze the grey nonequilibrium-diffusion approximation and the grey Eddington approximation and show that they both preserve this first-order accuracy. Further, these approximations preserve the EDA’s first-order accuracy when made in either the comoving-frame (CMF) or the lab-frame (LF). And while analyzing the Eddington approximation, we found that the CMF and LF radiation-source equations are equivalent when neglecting O(β 2 ) terms and compared in the LF. Of course, the radiation pressures are not equivalent. It is expected that simplified physical models and numerical discretizations of the RH equations that do not preserve this first-order accuracy will not retain the correct equilibrium-diffusion solutions. As a practical example, we show that nonequilibrium-diffusion radiative-shock solutions devolve to equilibrium-diffusion solutions when the asymptotic parameter is small.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Formation and evolution of protostellar accretion discs – I. Angular-momentum budget, gravitational self-regulation, and numerical convergence

Here, we investigate the formation and early evolution of a protostellar disc from a magnetized prestellar core using non-ideal magnetohydrodynamic (MHD) simulations including ambipolar diffusion and Ohmic dissipation. The dynamical contraction of the prestellar core ultimately leads to the formation of a first hydrostatic core, after ambipolar diffusion decouples the magnetic field from the predominantly neutral gas. The hydrostatic core accumulates angular momentum from the infalling material, evolving into a rotationally supported torus; this ‘first hydrostatic torus’ then forms an accreting protostar and a rotationally supported disc. The disc spreads out by gravitational instability, reaching ~30 au in diameter at ~3 kyr after protostar formation. The total mass and angular momentum of the protostar–disc system are determined mainly by accretion of gas from an infalling pseudo-disc, which has low specific angular momentum because of magnetic braking; their removal from the protostar–disc system by outflow and disc magnetic braking are negligible, in part because the magnetic field is poorly coupled there. The redistribution of angular momentum within the protostar–disc system is facilitated mainly by gravitational instability; this allows formation of relatively large discs even when the specific angular momentum of infalling material is low. We argue that such discs should remain marginally unstable as they grow (with Toomre Q ~ 1–2), an idea that is broadly consistent with recent observational estimates for Class 0/I discs. We discuss the numerical convergence of our results, and show that properly treating the inner boundary condition is crucial for achieving convergence at an acceptable computational cost.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A MAGNETIC RIBBON MODEL FOR STAR-FORMING FILAMENTS

We develop a magnetic ribbon model for molecular cloud filaments. These result from turbulent compression in a molecular cloud in which the background magnetic field sets a preferred direction. We argue that this is a natural model for filaments and is based on the interplay between turbulence, strong magnetic fields, and gravitationally driven ambipolar diffusion, rather than pure gravity and thermal pressure. An analytic model for the formation of magnetic ribbons that is based on numerical simulations is used to derive a lateral width of a magnetic ribbon. This differs from the thickness along the magnetic field direction, which is essentially the Jeans scale. We use our model to calculate a synthetic observed relation between apparent width in projection versus observed column density. The relationship is relatively flat, similar to observations, and unlike the simple expectation based on a Jeans length argument.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

On the eigenvector bias of Fourier feature networks: From regression to solving multi-scale PDEs with physics-informed neural networks

Physics-informed neural networks (PINNs) are demonstrating remarkable promise in integrating physical models with gappy and noisy observational data, but they still struggle in cases where the target functions to be approximated exhibit high-frequency or multi-scale features. Here in this work we investigate this limitation through the lens of Neural Tangent Kernel (NTK) theory and elucidate how PINNs are biased towards learning functions along the dominant eigen-directions of their limiting NTK. Using this observation, we construct novel architectures that employ spatio-temporal and multi-scale random Fourier features, and justify how such coordinate embedding layers can lead to robust and accurate PINN models. Numerical examples are presented for several challenging cases where conventional PINN models fail, including wave propagation and reaction–diffusion dynamics, illustrating how the proposed methods can be used to effectively tackle both forward and inverse problems involving partial differential equations with multi-scale behavior.

42 ENGINEERING↗
DOE OSTI2025

Nonintrusive projection-based reduced order modeling using stable learned differential operators

Nonintrusive projection-based reduced order models (ROMs) are essential for dynamics prediction in multi-query applications where underlying governing equations are known but the access to the source of the underlying full order model (FOM) is unavailable; that is, FOM is a glass-box. This article proposes a learn-then-project approach for nonintrusive model reduction. In the first step of this approach, high-dimensional stable sparse learned differential operators (S-LDOs) are determined using the generated data. In the second step, the ordinary differential equations, comprising these S-LDOs, are used with suitable dimensionality reduction and low-dimensional subspace projection methods to provide equations for the evolution of reduced states. This approach allows easy integration into the existing intrusive ROM framework to enable nonintrusive model reduction while allowing the use of Petrov–Galerkin projections. The applicability of the proposed approach is demonstrated for Galerkin and LSPG projection-based ROMs through four numerical experiments: 1-D scalar advection, 1-D Burgers, 2-D scalar advection and 1-D scalar advection–diffusion–reaction equations. In conclusion, the results indicate that the proposed nonintrusive ROM strategy provides accurate and stable dynamics prediction.

42 ENGINEERING↗
DOE OSTI2020

SubTSBR to tackle high noise and outliers for data-driven discovery of differential equations

Data-driven discovery of differential equations has been an emerging research topic. We propose a novel algorithm subsampling-based threshold sparse Bayesian regression (SubTSBR) to tackle high noise and outliers. The subsampling technique is used for improving the accuracy of the Bayesian learning algorithm. It has two parameters: subsampling size and the number of subsamples. When the subsampling size increases with fixed total sample size, the accuracy of our algorithm goes up and then down. When the number of subsamples increases, the accuracy of our algorithm keeps going up. We demonstrate how to use our algorithm step by step and compare our algorithm with threshold sparse Bayesian regression (TSBR) for the discovery of differential equations. We show that our algorithm produces better results. We also discuss the merits of discovering differential equations from data and demonstrate how to discover models with random initial and boundary condition as well as models with bifurcations. The numerical examples are: (1) predator-prey model with noise, (2) shallow water equations with outliers, (3) heat diffusion with random initial and boundary condition, and (4) fish-harvesting problem with bifurcations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Microstructure evolution and mechanical performance of nickel based superalloy C1023 at elevated temperatures

Highlights: • Growth and strengthening mechanisms of carbide in superalloy C1023 were investigated. • Carbide dendrite primary arm grows selectively based on the vertex of blocky carbonitride. • Monolayer or multilayer carbide dendrite arm has certain crystal relationship with γ′ and γ. • Plastic deformation of γ hinders interior crack propagating from MC to matrix at 650 °C. • Oxidation of GBs significantly weakens the tensile strength of the superalloy at 850 °C. - Abstract: Microstructure evolution of polycrystalline nickel based superalloy C1023 casting was investigated by comparing as cast superalloy and heat treated superalloy in the paper. Influences of strengtheners and carbides on the tensile behaviors of heat treated superalloy at different temperatures (20 °C, 650 °C, and 850 °C) were studied. Resulting data from deep corrosion suggested that growth mode of scripted carbide included nucleation on the vertex of blocky carbonitride, expansion of dendrite arm in C-rich area, and enlargement of dendrite cap. Decomposition of MC carbide generated numerous octahedral M{sub 23}C{sub 6} carbide as well as nanoscale Cr-Mo-rich nuclei. Diffusion of Mo from γ to precipitates such as γ', MC, and M{sub 23}C{sub 6}, occurred during the heat treatment while irregular γ' (4 nm–500 nm) transformed into cuboidal particles. At 650 °C, no significant deterioration of tensile strength was found while plastic deformation ability of γ base was improved. Oxidation of grain boundaries (GBs) dominated the crack behavior of the superalloy at 850 °C.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Thermodynamics–Structure–Dynamics Correlations and Nonuniversal Effects in the Elastically Collective Activated Hopping Theory of Glass-Forming Liquids

We employ the microscopic Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated dynamics in combination with crystal-avoiding simulations to study four inter-related questions for metastable monodisperse hard sphere fluids. The first is how significantly improved integral equation theory structural input (Modified-Verlet (MV) closure) changes the dynamical predictions of ECNLE theory. The main consequence is a modest enhancement of the importance of the collective elastic barrier relative to its local cage contribution, which increases the alpha relaxation time and fragility relative to prior results based on the Percus–Yevick closure. Second, ECNLE-MV theory predictions for the alpha time and self-diffusion constant in the metastable regime are quantitatively compared to our new simulations. The small adjustment of a numerical prefactor that enters the collective elastic barrier leads to quantitative agreement over three decades. Third, using the more accurate MV structural input, ECNLE theory is shown to predict thermodynamics–structure–dynamics “correlations” based on various long and short wavelength scalar properties all related to static two-point collective density fluctuations. The logarithm of the alpha relaxation time scales as a power law with these scalar metrics with an exponent that is significantly lower in the less dense noncooperative activated regime compared to the very dense highly cooperative regime. However, the discovered correlation of activated relaxation with a thermodynamic property (dimensionless compressibility) is not causal in ECNLE theory, but rather reflects a strong connection between the local structural quantities that quantify kinetic constraints in the theory with the amplitude of long wavelength density fluctuations. Fourth, the consequences of chemically specific nonuniversalities associated with the onset condition and relative importance of collective elasticity are studied. In conclusion, the predicted thermodynamics–structure–dynamics correlations are found to be robust, albeit with nontrivial shifts of the onset condition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗