Search NASA⌕ Search

SEARCH · Search NASA

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Effect of Flow on CO2 Corrosion of Self-Healing Metallic Coatings

Internal corrosion is an issue that affects natural gas pipelines, a significant part of the United States energy infrastructure. Over time, this corrosion has worn away longstanding pipelines due to the original construction materials used for the lines and impurities in the gas and liquid streams flowing through them. The main impurities in natural gas are H2O, CO2, H2S, and O2. One of the leading causes of corrosion is the CO2 dissolved in water, giving rise to carbonic acid formation, which can further dissolve the steel pipe. The National Energy Technology Laboratory (NETL) has been studying different solutions to this problem. One potential answer is a self-healing sacrificial metallic coating applied by a novel cold spray technique. Previous studies at NETL have shown that the zinc-rich self-healing coating can withstand corrosion through passivation and create a barrier to the diffusion of corrosive species. This study explored the impact of flow on corrosion with this coating and carbon steel when exposed to a saturated CO2 environment by simulating the pipeline flow profile in a small-scale lab setting. The samples were evaluated using multiple electrochemical techniques that found corrosion rates and were backed up with surface analysis to conclude the behavior of the corrosion mechanisms.

CO2 corrosion↗
DOE OSTI2017

SkyFACT: high-dimensional modeling of gamma-ray emission with adaptive templates and penalized likelihoods

We present SkyFACT (Sky Factorization with Adaptive Constrained Templates), a new approach for studying, modeling and decomposing diffuse gamma-ray emission. Like most previous analyses, the approach relies on predictions from cosmic-ray propagation codes like GALPROP and DRAGON. However, in contrast to previous approaches, we account for the fact that models are not perfect and allow for a very large number (∼> 10{sup 5}) of nuisance parameters to parameterize these imperfections. We combine methods of image reconstruction and adaptive spatio-spectral template regression in one coherent hybrid approach. To this end, we use penalized Poisson likelihood regression, with regularization functions that are motivated by the maximum entropy method. We introduce methods to efficiently handle the high dimensionality of the convex optimization problem as well as the associated semi-sparse covariance matrix, using the L-BFGS-B algorithm and Cholesky factorization. We test the method both on synthetic data as well as on gamma-ray emission from the inner Galaxy, |ℓ|<90{sup o} and | b |<20{sup o}, as observed by the Fermi Large Area Telescope. We finally define a simple reference model that removes most of the residual emission from the inner Galaxy, based on conventional diffuse emission components as well as components for the Fermi bubbles, the Fermi Galactic center excess, and extended sources along the Galactic disk. Variants of this reference model can serve as basis for future studies of diffuse emission in and outside the Galactic disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Anderson acceleration and application to the three-temperature energy equations

The Anderson acceleration method is an algorithm for accelerating the convergence of fixed-point iterations, including the Picard method. Anderson acceleration was first proposed in 1965 and, for some years, has been used successfully to accelerate the convergence of self-consistent field iterations in electronic-structure computations. Recently, the method has attracted growing attention in other application areas and among numerical analysts. Compared with a Newton-like method, an advantage of Anderson acceleration is that there is no need to form the Jacobian matrix. Thus the method is easy to implement. In this paper, an Anderson-accelerated Picard method is employed to solve the three-temperature energy equations, which are a type of strong nonlinear radiation-diffusion equations. Two strategies are used to improve the robustness of the Anderson acceleration method. One strategy is to adjust the iterates when necessary to satisfy the physical constraint. Another strategy is to monitor and, if necessary, reduce the matrix condition number of the least-squares problem in the Anderson-acceleration implementation so that numerical stability can be guaranteed. Numerical results show that the Anderson-accelerated Picard method can solve the three-temperature energy equations efficiently. Compared with the Picard method without acceleration, Anderson acceleration can reduce the number of iterations by at least half. A comparison between a Jacobian-free Newton–Krylov method, the Picard method, and the Anderson-accelerated Picard method is conducted in this paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Understanding NBI heating and fueling in LTX-β tokamak

My work on the contract DE-SC=0023274 (July, 2022 - July, 2025) was focused on understanding the plasma fueling by NBI. On Dec. 1998 with Sergei Krasheninnikov (UCSD, San Diego, CA) we initiated the Li Wall Fusion (LiWFusion) as a new concept of magnetic fusion [1]. The new 1 2 plasma 3 4 5 6 2.1 Overview of rate coefficients, mean-free paths λH0 , and diffusions coefficients 6 2. GSV code for ⟨σv⟩ analysis of NBI fueling concept was a reaction of the lack of luck of tokamak fusion with QDT = 1 on TFTR and JET. We recognized the edge plasma cooling by recycling to be the route reason of the tokamak problems on the way to burning plasma. LiWFusion relies on plasma pumping by a lithium layer on the inner walls of the plasma chamber, combined with the plasma heating and fueling by the Neutral Beam Injection (NBI). These two innovations eliminate the route problem of the tokmak fusion. The concept became theoretically mature in 2006. In 2012, the technology of continuously Flowing Liquid Lithium (24/7-FLiLi) was invented by me for the future implementations of the LWFusion concept. At the same time, my every presentation on LiWFusion to International Symposium on Lithium Applications since 2010 was objected by some person with words “Everybody since the 1970s knows that plasma fueling by NBI is impossible”. I dropped the name of the author of this objection and ignored his views. My work on the current grant on NBI fueling of LTX-β not only clarified the issue but resolved it in an astonishing way. The NBI fueling was fully understood (thus humiliating the current dogmatic fusion community). The new future of LTX-β dedicated for decades to Li in tokamaks, as well as of the entire magnetic fusion program was envisioned, in sharp contrast with the fallure of OFES in the post-TFTR era of 21st century.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Physics-informed Karhunen-Loeve and Neural Network Approximations for Solving Inverse Differential Equation Problems

Here we present the PI-CKL-NN method for parameter estimation in differential equation (DE) models given sparse measurements of the parameters and states. In the proposed approach, the space- or time-dependent parameters are approximated by Karhunen-Loeve (KL) expansions that are conditioned on the parameters’ measurements, and the states are approximated by deep neural networks (DNNs). The unknown weights in the KL expansions and DNNs are found my minimizing the cost function that enforces the measurements of the states the DE constraint. Regularization is achieved by adding the l2 norm of the conditional KL coefficients into the loss function. Our approach assumes that the parameter fields are correlated in space or time and enforces the statistical knowledge (the mean and the covariance function) in addition to the DE constraints and measurements as opposed to the physics-informed neural network (PINN) and other similar physics-informed machine learning methods where only DE constraints and data are used for parameter estimation. We use the PI-CKL-NN method for parameter estimation in an ordinary differential equation with an unknown time-dependent parameter and the one- and two-dimensional partial differential diffusion equations with unknown space-dependent diffusion coefficients. We also demonstrate that PI-CKL-NN is more accurate than the PINN method, especially when the observations of the parameters are very sparse

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Efficient data-driven regression for reduced-order modeling of spatial pattern formation

We present an efficient data-driven regression approach for constructing reduced-order models (ROMs) of reaction-diffusion systems exhibiting pattern formation. The ROMs are learned non-intrusively from available training data of physically accurate numerical simulations. The method can be applied to general nonlinear systems through the use of polynomial model form, while not requiring knowledge of the underlying physical model, governing equations, or numerical solvers. The process of learning ROMs is posed as a low-cost least-squares problem in a reduced-order subspace identified via Proper Orthogonal Decomposition (POD). Numerical experiments on classical pattern-forming systems–including the Schnakenberg and Mimura–Tsujikawa models–demonstrate that higher-order surrogate models significantly improve prediction accuracy while maintaining low computational cost. The proposed method provides a flexible, non-intrusive model reduction framework, well suited for the analysis of complex spatio-temporal pattern formation phenomena.

Data-driven modeling↗
DOE OSTI2022

A multi-resolution approach to hydraulic fracture simulation

Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.

58 GEOSCIENCES↗
DOE OSTI2025

Joint state-parameter estimation for the reduced fracture model via the united filter

Here, in this paper, we introduce an effective United Filter method for jointly estimating the solution state and physical parameters in flow and transport problems within fractured porous media. Fluid flow and transport in fractured porous media are critical in subsurface hydrology, geophysics, and reservoir geomechanics. Reduced fracture models, which represent fractures as lower-dimensional interfaces, enable efficient multi-scale simulations. However, reduced fracture models also face accuracy challenges due to modeling errors and uncertainties in physical parameters such as permeability and fracture geometry. To address these challenges, we propose a United Filter method, which integrates the Ensemble Score Filter (EnSF) for state estimation with the Direct Filter for parameter estimation. EnSF, based on a score-based diffusion model framework, produces ensemble representations of the state distribution without deep learning. Meanwhile, the Direct Filter, a recursive Bayesian inference method, estimates parameters directly from state observations. The United Filter combines these methods iteratively: EnSF estimates are used to refine parameter values, which are then fed back to improve state estimation. Numerical experiments demonstrate that the United Filter method surpasses the state-of-the-art Augmented Ensemble Kalman Filter, delivering more accurate state and parameter estimation for reduced fracture models. This framework also provides a robust and efficient solution for PDE-constrained inverse problems with uncertainties and sparse observations.

Bayesian inference↗
DOE OSTI2017

A comparative study of discrete fracture network and equivalent continuum models for simulating flow and transport in the far field of a hypothetical nuclear waste repository in crystalline host rock

One of the major challenges of simulating flow and transport in the far field of a geologic repository in crystalline host rock is related to reproducing the properties of the fracture network over the large volume of rock with sparse fracture characterization data. Various approaches have been developed to simulate flow and transport through the fractured rock. The approaches can be broadly divided into Discrete Fracture Network (DFN) and Equivalent Continuum Model (ECM). The DFN explicitly represents individual fractures, while the ECM uses fracture properties to determine equivalent continuum parameters. In this paper, we compare DFN and ECM in terms of upscaled observed transport properties through generic fracture networks. The major effort was directed on making the DFN and ECM approaches similar in their conceptual representations. This allows for separating differences related to the interpretation of the test conditions and parameters from the differences between the DFN and ECM approaches. The two models are compared using a benchmark test problem that is constructed to represent the far field (1 × 1 × 1 km 3 ) of a hypothetical repository in fractured crystalline rock. The test problem setting uses generic fracture properties that can be expected in crystalline rocks. The models are compared in terms of the: 1) effective permeability of the domain, and 2) nonreactive solute breakthrough curves through the domain. The principal differences between the models are mesh size, network connectivity, matrix diffusion and anisotropy. We demonstrate how these differences affect the flow and transport. Finally, we identify the factors that should be taken in consideration when selecting an approach most suitable for the site-specific conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Positivity enhancements to the P n transport equations

One of the common methods for solving the radiation transport equation is to use a polynomial expansion for the angle variable(s). Here, while each solution technique has its disadvantages, in some instances this method can result in oscillations large enough to give a negative radiation energy density, which is nonphysical. Some authors have recast the closure of the transport equations to guarantee positivity. Unfortunately, this results in messy nonlinear equations. Other authors post process each time step’s solution to improve positivity. In this work, the P n equations are modified to be more diffusive, less oscillatory, more positive, and more accurate. No new solution techniques are required, just the scaling of some terms with constant factors and the addition of constant filter opacities; therefore, the solutions are still linear. This new solution method is called stretched and filtered P n . Test problems are presented in one and two dimensions. Comparisons to the D n and TP n closures are made in one dimension. Modifications to D n and TP n are suggested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Hybrid super homogenization and discontinuity factor method for continuous finite element diffusion

This paper presents a novel homogenization equivalence technique aiming to simultaneously leveragethe simplicity of the Super Homogenization (SPH) method to reproduce reference reaction rates andthe ability to preserve reference leakage rates at desired surfaces through Discontinuity Factors (DF).The need for this new class of methods arises from the inability of the current state-of-the-art SPH technology to properly reproduce the reactor multiplication factor for problems with significant leakage. This work shows that this defect lies in the use of normalization factors in the SPH algorithm: while they are necessary with purely reflecting problems to ensure uniqueness of the solution, they introduce homogenization inconsistencies if at least one vacuum boundary condition is present. Two solutions to this problem are presented in this work: (i) to simply remove the normalization factor or (ii) to introduce additional degrees of freedom in the form of DFs in such a way that the normalization factors can still constrain the problem. While the former clearly offers unrivaled simplicity and can lead to very satisfactory results for leakage-dominated cores or if spatial restriction of the SPH regions is applied – as demonstrated for the High Temperature Test Reactor, the latter is more robust and does not require to scale the fission terms with the reference eigenvalue.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Residence time distribution in a structured packing unit for monitoring aerosol emissions

Solvent absorption in a structured packed column has been demonstrated as a viable option for post-combustion carbon capture. Amine emissions in the form of aerosol particles pose a major challenge that severely impacts the separation efficiency of the packed column. The aerosol particles quickly achieve the motion of the gas flow. In this context, residence time distribution (RTD) of the gaseous phase will provide insights for managing aerosol particles, thereby controlling amine loss in the packed column. Computational fluid dynamics (CFD) modeling of such columns is a multiscale problem because of a wide range of disparity in length scales. Accordingly, three-dimensional CFD simulations in the representative elementary unit of a structured packing are conducted using OpenFOAM. The RTD is analyzed in terms of residence time distribution function, mean residence time, and variance. Effects of the diffusivity and flow rates on RTD are extensively studied, and the dispersion coefficient is calculated using a one-dimensional axial dispersion model. Furthermore, RTDs computed in a smaller domain are extrapolated to a larger domain using a convolution integral. The computed and corresponding CFD models derived from the convolution integral match well at different flow rates and domain sizes. A comparison of RTD among different packing designs also is presented at a fixed number of repeating units and different flow rates.

structured packings, residence time distribution↗
DOE OSTI2025

A Comparison of Electronic Structure Methods for Predicting the Hydrogenation Energies of Candidate Molecules for Hydrogen Storage

The development of novel energy materials and fuels is required to expand current available energy sources. Aiming to reach this goal, there is growing interest in using molecular hydrogen as an energy carrier due to its abundance and high energy density. Liquid organic hydrogen carriers (LOHCs) are a promising route to the large-scale storage and transport of hydrogen for use in the energy economy. The search for thermodynamically viable LOHC molecules for real world use has led to a set of constraints on the dehydrogenation enthalpy and the minimum gravimetric hydrogen capacity. These constraints allow one to formulate the search for an ideal LOHC candidate molecule as an optimization problem well suited to the strengths of machine learning and artificial intelligence computational approaches. A critical barrier to a large-scale, high-throughput screening of LOHC candidate molecules is the lack of reliable training data. Computational electronic structure methods including density functional theory, coupled cluster approximations, and diffusion Monte Carlo can be used to provide training data where experimental data are either unreliable or do not exist. In this work, we use these methods to calculate the dehydrogenation energies and enthalpies of candidate LOHC molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Robust RF Photocathodes (Phase II Final Report)

This program addressed requirements for photocathodes with high quantum efficiency (QE) and long lifetimes. Currently, photocathodes with QEs above 1% exhibit lifetimes measured in a few hours. This is unacceptable for commercial devices using photocathodes as the electron source. This includes accelerators for industrial, medical, defense, and educational applications as well as emerging applications, such as time resolved electron microscopy. Calabazas Creek Research, Inc. (CCR) and the University of Maryland (UMD) attempted to develop reservoir photocathodes that dramatically increased lifetime and increased the ability to recover from contaminating events, such as loss of system vacuum. The Phase I program demonstrated the viability of the reservoir approach for higher QE cathodes. The Phase II program focused on design, fabrication, and delivery. This included increasing understanding of photocathode physics, investigating design options, addressing performance issues, and developing fabrication and packaging techniques leading to commercial products. The initial task identified specific photocathode materials that provide high QE and are compatible with the reservoir approach. The goal was that these photocathodes be manufacturable in a production environment, understanding that the “environment” is highly specialized. The photocathode geometry was modified to address several photogun performance issues related primarily to emittance. A goal was to reduce or eliminate thermal emittance issues arising from the reservoir heater. The program developed fabrication, handling, and transport techniques and facilities to allow commercial photocathode production for certain materials. Not only was complex chemistry involved, but practical issues of transporting and installing the cathodes into user systems were addressed. A critical goal was demonstration of reservoir photocathodes in operating RF guns, where performance could be directly measured. Success depended on clear demonstration that the cathodes could be produced, delivered, installed, and provide the required performance. The principal effort in this program was to develop reservoir cesium antimonide photocathodes. Cesium antimonide cathodes provide a work function in excess of the 1% goal, while the current generation exhibits lifetimes measured in hours. The reservoir in CCRs’ controlled porosity reservoir (CPR) cathodes provides a mechanism for replacing cesium evaporating from the surface by transporting replacement cesium from the reservoir. This is accomplished by gently heating the reservoir to induce cesium diffusion through micron-scale pores in the reservoir cap and photocathode emitter. Experiments at the UMD determined that cesium diffused from the reservoir would only react with antimony on the surface if the temperature exceeded 300C. At this temperature, however, the surface cesium evaporated faster than it could be replenished from the reservoir. Consequently, the QE decreased each time recessiation was attempted. This appears to be a fundamental problem with this approach. The sections below provide information on accomplishments for each task in the program. Once the issue with cesium antimonide cathodes was identified, CCR refocused the program to begin investigation of graphene coated photocathodes in collaboration with SLAC National Accelerator Laboratory (SLAC). The photocathode test chamber developed in the program was transferred to SLAC for use in the graphene research. This is described below in a task added to the program.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Transport properties of disordered two‐dimensional complex plasma crystal

In this study, we numerically investigate the transport properties of a two‐dimensional ( 2D ) complex plasma crystal using diffusion of coplanar dust lattice waves. In the limit where the Hamiltonian interactions can be decoupled from the non‐Hamiltonian effects, we identify two distinct types of wave transport: Anderson‐type delocalization and long‐distance excitation. We use a recently developed spectral approach to evaluate the contribution of the Anderson problem and compare it to the results of the simulation. The benefit of our approach to transport problems is twofold. First, we employ a highly tuneable macroscopic hexagonal crystal, which exhibits many‐body interactions and allows for the investigation of transport properties at the kinetic level. Second, the analysis of the transport problem in 2D is provided using an innovative spectral approach, which avoids the use of scaling and boundary conditions. The comparison between the analytically predicted and numerically observed wave dynamics allows for the study of important characteristics of this open system. In our simulations, we observe long‐distance lattice excitation, which occurs around lattice defects even when the initial perturbation does not spread from the centre to the exterior of the crystal. In the decoupled Hamiltonian regime, this many‐body effect can be attributed to the dust lattice interaction with the plasma environment.

Kostadinova, E. G.↗