Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion equation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Modelling of edge plasma dynamics with active wall boundary conditions

A self-consistent 2D model is presented for transport in boundary plasma and plasma-facing material walls. Plasma dynamics in the domain is represented by a 2D collisional plasma fluid model in the edge-plasma code UEDGE, and transport of hydrogen and heat in the wall is represented by a system of reaction–diffusion equations in the 1D wall code FACE. To account for variation of parameters along the wall, in the coupled model multiple instances of the FACE code run in parallel. Here, the coupled model provides a tool for investigating a range of dynamic plasma–material interactions phenomena in 2D. For demonstration of its capability, one application of particular interest is the role of active wall in tokamak strike point sweeping proposed for mitigation of divertor heat loads. In the present study, the coupled calculations are applied to investigation of the impact of heat and hydrogen transport in the material wall on the divertor plasma and target heat load during sweeping of the target strike point for parameters of a high-power tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Kinetic Simulations of Cu Doping in Chlorinated CdSeTe PV Absorbers

This paper reports on 1D kinetic modeling of p-type doping formation in Cu-doped chlorinated CdSeTe photovoltaic (PV) absorbers with graded Se profiles. Following the recent progress in kinetic simulations of the defect chemistry in stoichiometric CdTe films, this work extends simulation capabilities to a domain of semiconductor alloy films with graded stoichiometries. The defect formation energies and ionization levels, as well as the diffusion and reaction barriers used in the reaction-diffusion equations are calculated from the first principles. A new formalism is employed to account for the uncertainty of the defect formation energies in semiconductor alloys caused by unknown local surrounding. Finally, the developed methodology is used to model a practical fabrication process of p-type doping formation in graded CdSeTe absorbers and to study the effect of alloy stoichiometry on the doping activation.

14 SOLAR ENERGY↗
DOE OSTI2020

Learning compact physics‐aware delayed photocurrent models using dynamic mode decomposition

Abstract Radiation‐induced photocurrent in semiconductor devices can be simulated using complex physics‐based models, which are accurate, but computationally expensive. This presents a challenge for implementing device characteristics in high‐level circuit simulations where it is computationally infeasible to evaluate detailed models for multiple individual circuit elements. In this work we demonstrate a procedure for learning compact delayed photocurrent models that are efficient enough to implement in large‐scale circuit simulations, but remain faithful to the underlying physics. Our approach utilizes dynamic mode decomposition (DMD), a system identification technique for learning reduced‐order discrete‐time dynamical systems from time series data based on singular value decomposition. To obtain physics‐aware device models, we simulate the excess carrier density induced by radiation pulses by solving numerically the ambipolar diffusion equation, then use the simulated internal state as training data for the DMD algorithm. Our results show that the significantly reduced‐order delayed photocurrent models obtained via this method accurately approximate the dynamics of the internal excess carrier density—which can be used to calculate the induced current at the device boundaries—while remaining compact enough to incorporate into larger circuit simulations.

Hanson, Joshua↗
DOE OSTI2020

Nonoverlapping Localized Exponential Time Differencing Methods for Diffusion Problems

In this paper, we propose nonoverlapping localized exponential time differencing (ETD) methods for diffusion problems. The model time-dependent diffusion equation is first reformulated on subdomains based on the nonoverlapping domain decomposition, in which Neumann boundary conditions are imposed on the interfaces for the subdomain problems and Dirichlet type conditions are enforced to form a space-time interface problem. After spatial discretization by standard central finite differences and temporal integration with the first or second order ETD methods, the fully discrete interface problem is obtained. Such an interface problem is then solved iteratively either at each time step or over the whole time interval: the former involves the solution of stationary problems in each subdomain at each iteration while the latter involves the solution of time-dependent subdomain problems at each iteration. For both approaches, we prove that localized ETD solutions conserve mass exactly and converge in time to the exact space semidiscrete solution. Numerical experiments in two dimensions are also presented to illustrate the performance of the proposed methods.

Mass conservation↗
DOE OSTI2026

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗
DOE OSTI2018

An asymmetric backward problem for the inhomogeneous parabolic equation with time-dependent diffusivity

In this paper, we deal with an asymmetric case of the non-homogeneous backward parabolic problem associated with time-dependent diffusivity in polar coordinates which arises in describing the heat transfer in cylinder. In general, this problem is severely ill-posed by the Hadamard instability. To subdue the instability of this problem, we apply the modified quasi-boundary value method. According to some a priori assumptions on the exact solution, we get an explicit error estimate of Hölder type for all t∈(0,T]. In addition, a numerical experiment is given to illustrate the efficiency and flexibility of our method.

42 ENGINEERING↗
DOE OSTI2017

Pore-network model of evaporation-induced salt precipitation in porous media: The effect of correlations and heterogeneity

Salt transport and precipitation in porous media constitute a set of complex and fascinating phenomena that are of considerable interest to several important problems, ranging from storage of CO 2 in geological formations, to soil fertility, and protection of pavements and roads, as well as historical monuments. The phenomena occur at the pore scale and are greatly influenced by the heterogeneity of the pore space morphology. In this paper, we present a pore-network (PN) model to study the phenomena. Vapor diffusion, capillary effect at the brine-vapor interface, flow of brine, and transport of salt and its precipitation in the pores that plug the pores partially or completely are all accounted for. The drying process is modeled by the invasion percolation, while transport of salt in brine is accounted for by the convective-diffusion equation. We demonstrate that the drying patterns, the clustering and connectivity of the pore throats in which salt precipitation occurs, the saturation distribution, and the drying rate are all strongly dependent upon the pore-size distribution, the correlations among the pore sizes, and the anisotropy of the pore space caused by stratification that most natural porous media contain. In particular, if the strata are more or less parallel to the direction of injection of the gas that dries out the pore space (air, for example) and/or causes salt precipitation (CO 2 , for example), the drying rate increases significantly. Moreover, salt tends to precipitate in clusters of neighboring pores that are parallel to the open surface of the porous medium.

42 ENGINEERING↗
DOE OSTI2018

Control rod treatment for FEM based radiation transport methods

This paper presents a novel control rod treatment for finite element method (FEM) based reactor physics simulations. The new method falls into the class of homogenization methods which treat material interfaces within mesh elements by defining a representative, homogenized nuclear cross section set. In contrast to existing homogenization based cusping treatments that attempt to approximate flux-volume weighting of the cross sections, we propose to exactly integrate terms appearing in the variational form of the transport or diffusion equation over mixed material elements. This goal is achieved by defining smoothly varying cross sections in elements that straddle material boundaries. It is shown that the smoothly varying cross sections can be obtained by projecting the piecewise set of cross sections onto the set of Legendre polynomials up to a finite order. Implementation of smoothly varying cross sections is supported within the Rattlesnake radiation multiphysics code that is currently under development at the Idaho National Laboratory. Homogenization methods traditionally suffer from cusping effects that originate from overestimation of absorption cross sections in strongly neutron-absorbing materials as found in control rods. We demonstrate that the presented method is (1) highly effective in removing cusping effects and (2) provides accurate results based on the C5G7-TD benchmark, a time-dependent version of the C5G7 benchmark.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Orders of magnitude reduction in the thermal conductivity of polycrystalline diamond through carbon, nitrogen, and oxygen ion implantation

Despite the exceptional thermal and mechanical functionalities of diamond, its superlative properties are highly subject to the presence of point defects, dislocations, and interfaces. Here, polycrystalline diamond is ion implanted with C 3+ , N 3+ , and O 3+ ions at an energy of 16.5 MeV, producing an amorphous layer at the projected range and a damaged crystalline region between the surface and amorphous layer. Using time-domain thermoreflectance in combination with thermal penetration depth calculations based upon the multilayer heat diffusion equation, it is determined that reductions in the thermal conductivity can span nearly two orders of magnitude while still maintaining a polycrystalline structure within the regions thermally probed. Dynamical diffraction simulations of high-resolution x-ray diffraction measurements demonstrate the formation of a strained layer localized at the end of range, with much lower levels of strain near the surface. Furthermore, within the polycrystalline region above the amorphous layer, the average number of displacements-per-atom from the ion irradiation is found to be < 1%, with mass impurity concentrations much less than 1%. These low defect concentrations within the thermally probed region demonstrate the remarkably large impact that dilute levels of defects from the ion implantation can have on the thermal conductivity of diamond.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Multi-scale topology optimization of multi-material structures with controllable geometric complexity — Applications to heat transfer problems

This study presents a topology optimization method to design assemblies of periodic cellular materials with controllable geometric complexity. The framework is based on a novel multi-scale and multi-material design model in which the structure, the layout of the material subdomains, and their micro-structures are optimized concurrently. To allow a tight control over its geometrical complexity, the layout at the macro-scale is described by level-set fields, parameterized by geometric primitives. The micro-scale geometry is represented through a density approach. A nonlinear programming algorithm drives the optimization process using design sensitivities computed by the discrete adjoint method. The proposed design framework is studied with heat transfer problems. Practical design problems, such as a heat sink and a thermal storage unit with phase change are discussed. The macro-scale analysis model relies on a generalized transient diffusion equation. At the micro-scale, homogenization is used to compute equivalent material properties. The numerical examples in this study show that the optimized multi-scale multi-material layouts outperform the corresponding mono-scale structures while maintaining geometric simplicity at the macro-scale level.

42 ENGINEERING↗
DOE OSTI2020

An open source fast fluid dynamics model for data center thermal management

Although computational fluid dynamics (CFD) has been widely adopted to improve data center thermal management, the high computational demand limits its applications, such as multivariate optimal design and operation. Fast fluid dynamics (FFD), which has been applied for fast airflow simulation, shows great potential. However, few research applied FFD for optimal design and operation of data center thermal management. This research improves the FFD model for data centers and conducts a comprehensive evaluation and demonstration. First, the FFD model is improved by solving the advection and diffusion equations together using an upwind scheme instead of a semi-Lagrangian advection solver in the conventional FFD model. Second, new features for data centers are added, such as a pressure correction method to simulate plenum airflow and dynamic boundary conditions for IT racks. The new FFD model is first validated with two indoor environment cases and the results show that the new FFD model has slightly better overall prediction accuracy and faster speed compared to the conventional FFD model. It is also observed that both FFD models achieve acceptable accuracy, except for a few localized disparities with experimental data, which might be due to simplified handling of turbulence viscosity near the boundaries. Furthermore, validation with a real data center shows that the FFD model achieves a similar level of accuracy as CFD when compared to the experimental measurements with some level of uncertainties. It is then demonstrated for data center optimal design and operation, which saves 53.4–58.8% of annual energy while still meeting the thermal requirements. In conclusion, with a much faster speed and comparable accuracy compared to CFD, the FFD model parallelized on a graphics processing unit is promising for practical model-based data center early design and operation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2019

TRANSP-based closed-loop simulations of current profile optimal regulation in NSTX-Upgrade

Active control of the toroidal current density profile is critical for the upgraded National Spherical Torus eXperiment device (NSTX-U) to maintain operation at the desired high-performance, MHD-stable, plasma regime. Initial efforts towards current density profile control have led to the development of a control-oriented, physics-based, plasma-response model, which combines the magnetic diffusion equation with empirical correlations for the kinetic profiles and the non-inductive current sources. The developed control-oriented model has been successfully tailored to the NSTX-U geometry and actuators. Moreover, a series of efforts have been made towards the design of model-based controllers, including a linear-quadratic-integral optimal control strategy that can regulate the current density profile around a prescribed target profile while rejecting disturbances. In this work, the tracking performance of the proposed current-profile optimal controller is tested in numerical simulations based on the physics-oriented code TRANSP. Furthermore, these high-fidelity closed-loop simulations, which are a critical step before experimental implementation and testing, are enabled by a flexible framework recently developed to perform feedback control design and simulation in TRANSP.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Pressure falloff testing to characterize CO 2 plume and dry-out zone during CO 2 injection in saline aquifers

The fluid mobilities and the lateral extent of the CO 2 plume and the dry-out zone are among important unknowns required to effectively manage CO 2 storage in saline aquifers projects. These unknowns can be derived from pressure falloff testing which has been widely used to characterize the subsurface. This paper provides a three-region analytical solution for falloff pressure during the infinite-acting period of CO 2 injection in saline aquifers. The solution is derived using Laplace transformation technique of the governing diffusivity equations and the accompanying initial and boundary conditions. The exact solution is validated against numerical results obtained from a black-oil reservoir simulator for two cases. The solution is approximated into simple forms in real-time domain considering different time periods. Then, a graphical semi-log analysis technique based on the approximate solutions is provided as an interpretation tool of falloff tests where the fluid mobilities and the extents of dry-out zone and CO 2 plume can be inferred. The interpretation technique is applied to synthetic falloff test data to investigate its potential and reliability.

02 PETROLEUM↗
DOE OSTI2021

Pulsating Poiseuille flow of a cement slurry

In this paper, we investigate the pulsatile flow of a cement slurry in a pipe. The constitutive relation for the viscous stress tensor is based on the power-law model, where the shear viscosity not only depends on the shear rate but also the volume fraction of the cement particles. To solve for the volume fraction field, a convection-diffusion equation is used. The dimensionless form of the governing equations and the boundary conditions are solved numerically at different pulsatile cycles. Further, a parametric study is conducted to study the effect of different dimensionless numbers at different cycles on the velocity and the volume fraction profiles.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Impact of parametric uncertainty on estimation of the energy deposition into an irradiated brain tumor

Predictions of the total energy deposited into a brain tumor through X-ray irradiation are notoriously error-prone. We investigate how this predictive uncertainty is affected by uncertainty in both the location of the region occupied by a dose-enhancing iodinated contrast agent and the agent's concentration. This is done within the probabilistic framework in which these uncertain parameters are modeled as random variables. We employ the stochastic collocation (SC) method to estimate statistical moments of the deposited energy in terms of statistical moments of the random inputs, and the global sensitivity analysis (GSA) to quantify the relative importance of uncertainty in these parameters on the overall predictive uncertainty. A nonlinear radiation–diffusion equation dramatically magnifies the coefficient of variation of the uncertain parameters, yielding a large coefficient of variation for the predicted energy deposition. This demonstrates that accurate prediction of the energy deposition requires a proper treatment of even small parametric uncertainty. Our analysis also reveals that SC outperforms standard Monte Carlo, but its relative efficiency decreases as the number of uncertain parameters increases from one to three. A robust GSA ameliorates this problem by reducing this number.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Self-induced spin-orbit torques in metallic ferromagnets

Herein we present a phenomenological theory of spin-orbit torques in a metallic ferromagnet with spin-relaxing boundaries. The model is rooted in the coupled diffusion of charge and spin in the bulk of the ferromagnet, where we account for the anomalous Hall effects as well as the anisotropic magnetoresistance in the corresponding constitutive relations for both charge and spin sectors. The diffusion equations are supplemented with suitable boundary conditions reflecting the spin-sink capacity of the environment. In inversion-asymmetric heterostructures, the uncompensated spin accumulation exerts a dissipative torque on the order parameter, giving rise to a current-dependent linewidth in the ferromagnetic resonance with a characteristic angular dependence. We compare our model to recent spin-torque ferromagnetic resonance measurements, illustrating how rich self-induced spin-torque phenomenology can arise even in simple magnetic structures.

36 MATERIALS SCIENCE↗