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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling fission product diffusion in TRISO fuel particles with BISON

Diffusion of fission products in intact TRISO particles depends on particle geometry, fission product source rates, time, temperature, and temperature-dependent diffusion coefficients. Simulating this diffusion process requires models for source rates and diffusion coefficients, plus computation of the temperature field if not prescribed. In addition, simulation quality depends on discretization of the geometry, appropriate time stepping, and the accuracy of the solution method. In this paper, we explore the simulation of fission product diffusion in TRISO fuel particles using the finite element method via the fuel performance code Bison. Recent material model development has occurred in Bison for each material present in tri-structural isotropic (TRISO) fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. Also, new mesh generation and fission product release fraction capabilities have been added. Diffusion capabilities are shown to converge to the correct solution via formal verification tests. A large number of code benchmarking problems are also given, with good results, showing that Bison’s computed release fractions closely match those of other software tools. Finally, a significant validation effort is detailed in which fission product release, measured as part of the AGR-1 capsule experiments, is compared to Bison outputs. Bison outputs compare very well to the experimental data and to PARFUME results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Analysis of operator splitting errors for near-limit flame simulations

High-fidelity simulations of ignition, extinction and oscillatory combustion processes are of practical interest in a broad range of combustion applications. Splitting schemes, widely employed in reactive flow simulations, could fail for stiff reaction–diffusion systems exhibiting near-limit flame phenomena. The present work first employs a model perfectly stirred reactor (PSR) problem with an Arrhenius reaction term and a linear mixing term to study the effects of splitting errors on the near-limit combustion phenomena. Analysis shows that the errors induced by decoupling of the fractional steps may result in unphysical extinction or ignition. The analysis is then extended to the prediction of ignition, extinction and oscillatory combustion in unsteady PSRs of various fuel/air mixtures with a 9-species detailed mechanism for hydrogen oxidation and an 88-species skeletal mechanism for n-heptane oxidation, together with a Jacobian-based analysis for the time scales. The tested schemes include the Strang splitting, the balanced splitting, and a newly developed semi-implicit midpoint method. Results show that the semi-implicit midpoint method can accurately reproduce the dynamics of the near-limit flame phenomena and it is second-order accurate over a wide range of time step size. For the extinction and ignition processes, both the balanced splitting and midpoint method can yield accurate predictions, whereas the Strang splitting can lead to significant shifts on the ignition/extinction processes or even unphysical results. With an enriched H radical source in the inflow stream, a delay of the ignition process and the deviation on the equilibrium temperature are observed for the Strang splitting. On the contrary, the midpoint method that solves reaction and diffusion together matches the fully implicit accurate solution. The balanced splitting predicts the temperature rise correctly but with an over-predicted peak. For the sustainable and decaying oscillatory combustion from cool flames, both the Strang splitting and the midpoint method can successfully capture the dynamic behavior, whereas the balanced splitting scheme results in significant errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Error Estimates for the Optimal Control of a Parabolic Fractional PDE

In this work, we consider the integral definition of the fractional Laplacian and analyze a linear-quadratic optimal control problem for the so-called fractional heat equation; control constraints are also considered. We derive existence and uniqueness results, first order optimality conditions, and regularity estimates for the optimal variables. To discretize the state equation we propose a fully discrete scheme that relies on an implicit finite difference discretization in time combined with a piecewise linear finite element discretization in space. We derive stability results and a novel $L^2(0,T;L^2(\Omega))$ a priori error estimate. On the basis of the aforementioned solution technique, we propose a fully discrete scheme for our optimal control problem that discretizes the control variable with piecewise constant functions, and we derive a priori error estimates for it. We illustrate the theory with one- and two-dimensional numerical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2021

Galvanostatic Intermittent Titration Technique Reinvented: Part I. A Critical Review

The galvanostatic intermittent titration technique (GITT), introduced in 1977 by Weppner and Huggins, provided a readily accessible means to measuring the chemical diffusion coefficient of electrochemical electrode materials. The method continues to be widely used today, but the reported diffusivity values are highly inconsistent, ranging as much as four orders of magnitude for some Li layered oxide compositions. Even qualitative trends of diffusivity are inconsistent, suggesting significant flaws in the implementation of the method. Other variants of the GITT method also suffer from similar inconsistency problems. Here we identify numerous sources of significant error including composition-dependent reaction overpotentials, mathematical flaws in the relaxation analysis methods, finite-size and non-planar geometry effects, inter-particle inhomogeneity issues, early transient effects, and surface area uncertainties. We propose a simple relaxation analysis scheme using the time variable t relax + τ - t relax , where t relax is relaxation time and τ is the galvanostatic pulse duration. We also propose to use dense diffusion-limited samples to isolate the bulk-diffusion process in the time domain. Chemical diffusivity can be extracted much more reliably with this improved implementation of the GITT method.

25 ENERGY STORAGE↗
DOE OSTI2019

A Finite-Element Approach for the Partitioning of Carbon in Q&P Steel

This paper reports a Galerkin finite-element analysis of carbon partitioning from martensite into austenite during the quenching and partitioning (Q&P) processes in steels. In contrast to classical or sophisticated diffusion field models, an alternative nonlinear governing equation based on chemical potential and composition is considered. The model is applied to simulate the carbon partitioning of modified 22MnB5 alloys assuming an immobile austenite–martensite phase boundary and, in turn, is compared with experimental measurements of the volume fraction and carbon content of retained austenite. The simulations show outstanding results on the influence of Si and C content on the partitioning parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modulating effect of evanescent waves on thin film growth

Atomic-scale smooth thin films are keys to successful integration and proper function of many multilayer-structured devices. However, the intrinsic islandlike growth mode prevents human beings from realizing such ultrasmooth films of many important functional materials. To solve this problem, we propose a negative entropy-infusing method that employs evanescent waves to enhance the downward interlayer diffusion of adatoms and thus transform the islandlike growth mode to the layer-by-layer growth mode. In this work, the formulas of the optical force and the lowered diffusion barrier were derived, and the application of this theory on a simplified example demonstrates significantly improved surface morphologies through numerical simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Random Walks With Tweedie: A Unified View of Score-Based Diffusion Models [In the Spotlight]

We present a concise derivation for several influential score-based diffusion models that relies on only a few textbook results. Diffusion models have recently emerged as powerful tools for generating realistic, synthetic signals—particularly natural images—and often play a role in state-of-the-art algorithms for inverse problems in image processing. While these algorithms are often surprisingly simple, the theory behind them is not, and multiple complex theoretical justifications exist in the literature. Here, in this study, we provide a simple and largely self-contained theoretical justification for score-based diffusion models that is targeted towards the signal processing community. This approach leads to generic algorithmic templates for training and generating samples with diffusion models. We show that several influential diffusion models correspond to particular choices within these templates and demonstrate that alternative, more straightforward algorithmic choices can provide comparable results. This approach has the added benefit of enabling conditional sampling without any likelihood approximation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Modeling of Nonlocal Electron Conduction for Inertial Confinement Fusion

In this report Chenhall presents a new Monte Carlo (MC) transport method developed as part of his PhD thesis work funded via fellowship by Sandia National Laboratories. The model is used to simulate the electron thermal transport within inertial confinement fusion (ICF) type problems. The new model aims to improve upon the currently used implicit Schurtz, Nicolai, and Busquet (iSNB) method, in particular by using finite particle ranges in comparison to the exponential solution of a diffusion method and by improved higher order angular modeling. The new method is being developed using the 1D LILAC and 2D DRACO codes developed by the Rochester Laboratory for Laser Energetics. The model is compared to iSNB for several ICF type problems: Omega shot 60303 a shock timing experiment, Omega shot 68951 a cryo target implosion and a NIF phase plate polar direct drive design. Overall, the MC method performs at least as well as the iSNB method and appears to indicate a lower predicted preheat ahead of the shock fronts. However, the difference in overall results between iSNB and MC transport is small for the choice of test problems and further simulations will need to be run to fully understand the difference between the models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Recent Progress on Numerical Modeling for Microgravity Electric Field Flames: Preprint

This paper presents the recent progress on the NASA Physical Science Informatics (PSI) project PeleLM CFD of Ion Driven Winds from Diffusion Flames in simulating the E-FIELD Flames microgravity results. The focus of this project is to comprehensively simulate the behavior of a small diffusion flame under the influence of an externally applied electric field in zero-gravity. To date, the capability of accurately simulating electric field flames has eluded researchers because the system exhibits dramatic ranges of coupled temporal and spatial scales. Moreover, in earth gravity the hot combustion products are subject to buoyancy effects that are difficult to isolate from those generated by the electric field. This work is an implementation in an existing powerful simulation framework (PeleLM) for this problem using PeleLMeX in order to validate the model and investigate the complex coupled system. The current effort includes establishing the domain of a diffusion coflow burner, examining the boundary conditions, flame geometry and ignition with gravitational forces and also with an electric field applied. A detailed model is used that includes the chemistry of charged ions and chemiluminescent flame intermediates to capture any feedback between ion-driven convection and combustion behavior, and to allow quantitative comparisons with experimental measurements.

coflow flames↗
DOE OSTI2018

Diffusion of Innovations: Interplay of Social, Economic, Technological, and Policy Drivers in the Solar Industry. Summary of UT Austin Student Capstone Research Projects

The University of Texas at Austin’s Policy Research Project (PRP), a nine-month (two semesters) capstone, is a keystone of the core curriculum at the LBJ School of Public Affairs. In PRPs, small groups of students, under the mentorship of a faculty director, take on real-world problems that require special knowledge and skill sets. PRPs expose students to challenges in formulating and executing research, and in communicating academic research and related complex data to broader stakeholder communities and decision makers. The PRP structure is an innovative and effective approach for integrating research within the teaching and training of graduate students who are preparing themselves to address important real-world problems at the intersection of society, economics, technology, and policy. The project summaries below describe seven papers developed during September 2017 – May 2018 as part of a PRP on “Diffusion of Innovations: Interplay of Social, Economic, Technological, and Policy Drivers in the Solar Industry.” Twenty graduate students, drawn from the LBJ School’s Masters in Public Affairs and Masters in Global Policy Studies programs and the Jackson School Geoscience’s Energy and Earth Resources program, participated in this PRP. Dr. Varun Rai, Associate Professor and Associate Dean for Research at the LBJ School, directed the PRP, with support from his research team including: Dr. Ariane Beck, Dr. Ashok Sekar, D. Cale Reeves, and Erik Funkhouser. Clients for the project included the U.S. Department of Energy (Casey Canfield), Lawrence Berkeley National Laboratory (Ben Hoen, Galen Barbose Joachim Seel, Naïm Darghouth, Ryan Wiser), and National Renewable Energy Laboratory (Benjamin Sigrin, Eric O’Shaughnessy). The seven projects separately addressed one of the following topics: (1) low- and middle-income PV adoption, (2) modeling economic and information intervention design, (3) evaluation of DOE’s Solar in Your Community Challenge, (4) property value impacts near large-scale solar facilities, (5) solar market maturity and evolution of business models, (6) social media data for predicting PV adoption, and (7) individual-level variation in adoption of innovations. Many of the papers relied on data collected and curated by Lawrence Berkeley National Laboratory, including data embedded within the annual Tracking the Sun and Utility-Scale Solar reports. Each of the seven teams in the PRP prepared a research paper. The PRP culminated with a full-day conference at UT Austin in May 2018 to present findings from the seven projects in this PRP to a broad audience of about 75 experts from academia, national labs, industry, and government from across the country.

14 SOLAR ENERGY↗
DOE OSTI2017

Convergence Study of LR-NDA Using Fourier Analysis

Numerical solution of neutron transport equations is challenging. It often requires acceleration techniques to improve computational efficiency. In the last two decades, various acceleration methods have been developed. Among those methods, the coarse-mesh nonlinear diffusion acceleration (CM-NDA) technique, or commonly known as coarse-mesh finite difference (CMFD), becomes one of the most popular methods due to its effectiveness and simplicity in recent years. However, CM-NDA will lose its effectiveness and could become unstable when the mesh optical thickness becomes large. There exists a number of stabilization techniques to improve the stability of CM-NDA, such as pCMFD and most recent odCMFD. Wang et al., proposed a new CM-NDA scheme, LR-NDA, which employs local refinement on the coarse mesh. The new scheme can significantly improve the effectiveness and stability of CM-NDA. In this paper, we perform a convergence study of LR-NDA using Fourier analysis. The spectral radius calculated from Fourier analysis is consistent with the numerical results for a 1-D fixed source model problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Excitation Energies of Localized Correlated Defects via Quantum Monte Carlo: A Case Study of Mn 4+ -Doped Phosphors

Accurate excitation energies of localized defects have been a long-standing problem for electronic structure calculation methods. Using Mn 4+ -doped solids as our proof of principle, we show that diffusion quantum Monte Carlo (DMC) is able to predict phosphorescence emission energies within statistical error. To demonstrate the generality of our DMC approach for other possible localized defects, we conduct charge density analyses using DMC and density functional theory (DFT). We also identify a new material with an emission energy of 1.97(8) eV, which is close to the optimum of 2.03 eV for a red-emitting phosphor. To our knowledge, our work is the first report on studying excitation energies of a transition metal impurity using an ab initio many-body electronic structure method. In contrast, semilocal and hybrid-DFT largely underestimates, and fails to reproduce, some of the trends in the emission energies. In conclusion, our work underscores the importance of an accurate account of exchange, correlation, and excitonic effects for localized excitations in defective solids.

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

Analytical Sensitivity Analysis of a Spent Nuclear Fuel Cask

Here, we report nuclear science and engineering is a field increasingly dominated by computational studies resulting from increasingly powerful computational tools. As a result, analytical studies, which previously pioneered nuclear engineering, are increasingly viewed as secondary or unnecessary. However, analytical solutions to reduced-fidelity models can provide important information concerning the underlying physics of a problem and aid in guiding computational studies. Similarly, there is increased interest in sensitivity analysis studies. These studies commonly use computational tools. However, providing a complementary sensitivity study of relevant analytical models can lead to a deeper analysis of a problem. This work provides the analytical sensitivity analysis of the one-dimensional (1D) cylindrical mono-energetic neutron diffusion equation using the forward sensitivity analysis procedure (FSAP) developed by Cacuci. Further, these results are applied to a reduced-fidelity model of a spent nuclear fuel cask, demonstrating how computational analysis might be improved with a complementary analytic sensitivity analysis.

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

Theoretical Analysis of the Entrainment–Mixing Process at Cloud Boundaries. Part II: Motion of Cloud Interface

The entrainment–mixing process leads to formation of a cloud–environment interface zone, consisting of the cloud dilution zone, where the liquid water content decreases as the distance from the cloud interior increases, and a humid shell with higher humidity in comparison with the dry environment. The time evolution of the location and width of these zones is analyzed by solving the diffusion–evaporation equation for the open region in the vicinity of the cloud–dry air interface. Upon normalization, the problem is reduced to a one-parametric one, the governing parameter being the potential evaporation parameter R < 0, which is proportional to the ratio of saturation deficit in the dry air to the available liquid water content in the cloud air. It is shown that the widths of the dilution zone and the humid shell increase with time. At R < −1, the interface between the dilution zone and the humid shell (i.e., the cloud edge) moves toward the cloud interior (i.e., the cloud dissipates). If R > −1, the cloud edge moves outward; that is, the mixing leads to an increase in the cloud volume. It is shown that the relative humidity remains high within most of the cloud dilution zone, and the effective radius changes only slightly, especially in growing clouds, where the effective radius changes by less than 10% within more than 80% of the cloud dilution zone. The strong decrease in the LWC is accompanied by a corresponding decrease in the droplet concentration. The study is illustrated by means of two examples of clouds: a small cumulus (Cu) typical of the trade wind zone and a deep convective cloud typical of the monsoon period over India. It is shown that while the small Cu tends to dissipate, the deep convective cloud expands. The characteristic sizes of the dilution zone and the humid shell are several tens of meters for the small Cu and several hundred meters for the deep convective cloud.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Polarization Tomography of Residual Stresses in Cylindrical Gradient-Index Lenses

In most cases, cylindrical structures with a gradient refractive index (GRIN) are formed on the basis of diffusion technology. The refractive index and the residual stresses in them are axially symmetrical. The distortion determined by diffusion and the thermal internal stresses in them are simulated by the fictive temperature of residual stresses. The use of the integral photoelasticity for the reconstruction of residual stresses in GRIN lenses is usually limited to the case of a planar deformed state. In this study, this algorithm is extended to the case of the axial change of residual stresses. It is assumed that the illumination procedure is carried out in the sample cross section. The optical problem is solved within the quasi-isotropic approximation. The reconstruction is based on measuring the characteristic parameters of polarized light passed through the sample and on the combined solution of the thermoelasticity problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning

The reaction-diffusion master equation (RDME) is a model that allows for efficient on-lattice simulation of spatially resolved stochastic chemical kinetics. Compared to off-lattice hard-sphere simulations with Brownian dynamics or Green’s function reaction dynamics, the RDME can be orders of magnitude faster if the lattice spacing can be chosen coarse enough. However, strongly diffusion-controlled reactions mandate a very fine mesh resolution for acceptable accuracy. It is common that reactions in the same model differ in their degree of diffusion control and therefore require different degrees of mesh resolution. This renders mesoscopic simulation inefficient for systems with multiscale properties. Mesoscopic-microscopic hybrid methods address this problem by resolving the most challenging reactions with a microscale, off-lattice simulation. However, all methods to date require manual partitioning of a system, effectively limiting their usefulness as “black-box” simulation codes. In this paper, we propose a hybrid simulation algorithm with automatic system partitioning based on indirect a priori error estimates. We demonstrate the accuracy and efficiency of the method on models of diffusion-controlled networks in 3D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Acceleration of diffusive molecular dynamics simulations through mean field approximation and subcycling time integration

Highlights: • Two numerical methods are presented to accelerate DMD simulations. • The methods greatly reduce the computation time without much loss of accuracy. • The accelerated DMD model is capable of simulating long-term mass transport at atomistic length scales. • Demonstrated for the propagation of an atomically sharp phase boundary over 30 seconds. - Abstract: Diffusive Molecular Dynamics (DMD) is a class of recently developed computational models for the simulation of long-term diffusive mass transport at atomistic length scales. Compared to previous atomistic models, e.g., transition state theory based accelerated molecular dynamics, DMD allows the use of larger time-step sizes, but has a higher computational complexity at each time-step due to the need to solve a nonlinear optimization problem at every time-step. This paper presents two numerical methods to accelerate DMD based simulations. First, we show that when a many-body potential function, e.g., embedded atom method (EAM), is employed, the cost of DMD is dominated by the computation of the mean of the potential function and its derivatives, which are high-dimensional random variables. To reduce the cost, we explore both first- and second-order mean field approximations. Specifically, we show that the first-order approximation, which uses a point estimate to calculate the mean, can reduce the cost by two to three orders of magnitude, but may introduce relatively large error in the solution. We show that adding an approximate second-order correction term can significantly reduce the error without much increase in computational cost. Second, we show that DMD can be significantly accelerated through subcycling time integration, as the cost of integrating the empirical diffusion equation is much lower than that of the optimization solver. To assess the DMD model and the numerical approximation methods, we present two groups of numerical experiments that simulate the diffusion of hydrogen in palladium nanoparticles. In particular, we show that the computational framework is capable of capturing the propagation of an atomically sharp phase boundary over a time window of more than 30 seconds. The effects of the proposed numerical methods on solution accuracy and computation time are also assessed quantitatively.

36 MATERIALS SCIENCE↗