Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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

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

At least 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High-Fidelity Simulations of Plasma-Assisted Oxidation of Hydrocarbon Fuels Using Nanosecond Pulsed Discharges

Next, a robust and efficient framework for simulating NSPD in multiple dimensions is developed. The reactive Navier-Stokes equations are extended to include a drift-diffusion plasma-fluid model with a local field approximation (LFA) in a finite-volume solver, which uses an adaptive mesh refinement (AMR) strategy to address the wide separation of length scales in the problem. A two-way coupling strategy is used whereby the plasma-fluid model and reactive Navier-Stokes equations are integrated simultaneously. The oxidation of ethylene/air mixtures mediated by NSPD is simulated in a pin-to-pin configuration. All phases of the plasma discharge are simulated explicitly (including streamer ignition, propagation, and connection, as well as the subsequent spark phase), along with the evolution of the plasma during the inter-pulse period. Temporally and spatially-resolved results are presented, with an emphasis on the analysis of heating and energy deposition, as well as of the evolution of the concentration of active particles generated during the NSPD and their influence on ignition.

adaptive mesh refinement↗
DOE OSTI2014

Interface design of VSOP'94 computer code for safety analysis

Today, most software applications, also in the nuclear field, come with a graphical user interface. VSOP'94 (Very Superior Old Program), was designed to simplify the process of performing reactor simulation. VSOP is a integrated code system to simulate the life history of a nuclear reactor that is devoted in education and research. One advantage of VSOP program is its ability to calculate the neutron spectrum estimation, fuel cycle, 2-D diffusion, resonance integral, estimation of reactors fuel costs, and integrated thermal hydraulics. VSOP also can be used to comparative studies and simulation of reactor safety. However, existing VSOP is a conventional program, which was developed using Fortran 65 and have several problems in using it, for example, it is only operated on Dec Alpha mainframe platforms and provide text-based output, difficult to use, especially in data preparation and interpretation of results. We develop a GUI-VSOP, which is an interface program to facilitate the preparation of data, run the VSOP code and read the results in a more user friendly way and useable on the Personal 'Computer (PC). Modifications include the development of interfaces on preprocessing, processing and postprocessing. GUI-based interface for preprocessing aims to provide a convenience way in preparing data. Processing interface is intended to provide convenience in configuring input files and libraries and do compiling VSOP code. Postprocessing interface designed to visualized the VSOP output in table and graphic forms. GUI-VSOP expected to be useful to simplify and speed up the process and analysis of safety aspects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Derivation of the Point-Kinetics Parameters for the Coupled Reactors

The design of coupled reactors, which provide two distinct spectral zones (fast and thermal), has recently drawn attention because it can provide configurations that could not be met in current thermal neutron test reactors and thus allows versatile irradiation tasks. A coupled reactor core is characterized by the neutron multiplicity of each zone as well as by the coupling between the two zones. It is desired to evaluate the integral point kinetics parameters that characterize the neutron physics of the coupled system. This summary provides a derivation of the equations for evaluating these point kinetics parameters. Following the strong and the weak form of the multigroup diffusion equation, the weak forms for the partial forward and adjoint solutions, each for a particular active region of the reactor core, are derived. These partial solutions are then used to collapse the weak form into the coupled-reactor point kinetics equations with the approximation of equilibrium. It is noted that these parameters are not meant for transient analysis because of the equilibrium approximation. They are rather used for characterizing the coupled system. The derived reactivity equations are then compared with those in Ref. 1. The equations are implemented in Rattlesnake, the INL MOOSE-based radiation transport application for multi-physics simulations. Numerical results are presented for a simple two-group problem. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Uncovering grain and subgrain microstructure at the scale of additive manufacturing melt tracks with a scalable cellular automaton solidification model

Metal additive manufacturing, characterized by rapid solidification, yields refined grains with a distinctive cellular subgrain microstructure that plays a pivotal role in determining material properties. Due to the significant computational expense demanded to simulate the required physics with submicron spatial resolution, their numerical simulations have been limited to proof-of-concept studies to either 2D or small subregions of a melt pool. In this study, an open-source, scalable, solidification code, muMatScale, based on the cellular automaton method, has been developed to predict the grain and the underlying subgrain microstructure over an entire melt pool. The model incorporates flexible parallelization schemes, utilizing MPI and OpenMP GPU Offloading, in addition to appropriate multi-physics specific to non-equilibrium rapid solidification in AM. The impact of nucleation parameters on grain microstructures was investigated with a focus on grain size variations and morphology transitions. With selected nucleation parameters, the simulation predicted the grain size, subgrain morphology, crystallographic orientation, and microsegregation aligned with experimental measurements. The model demonstrates that epitaxial grain growth is a dominant factor at the melt pool boundary, influencing grain size variation under different grain sizes in the build plate while maintaining consistent primary dendrite arm spacing under identical thermal conditions. Here, the highly efficient numerical model enables large-scale simulations with a spatial resolution of 100 nm or less, unveiling unprecedented insights into thermal and solutal diffusion driven grain growth, and the subgrains with microsegregation within grains in 3D across scales. muMatScale will enable the linking of submicron length-scale microstructure to part-level material behavior by investigating fundamental solidification problems at the intercellular scale in many-track and many-layer builds.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A finite element method for angular discretization of the radiation transport equation on spherical geodesic grids

Discrete ordinate (S N ) and filtered spherical harmonics (FP N ) based schemes have been proven to be robust and accurate in solving the Boltzmann transport equation but they have their own strengths and weaknesses in different physical scenarios. We present a new method based on a finite element approach in angle that combines the strengths of both methods and mitigates their disadvantages. The angular variables are specified on a spherical geodesic grid with functions on the sphere being represented using a finite element basis. A positivity-preserving limiting strategy is employed to prevent non-physical values from appearing in the solutions. Here, the resulting method is then compared with both S N and FP N schemes using four test problems and is found to perform well when one of the other methods fail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Data-Driven First-Principles Methods for the Study and Design of Alkali Superionic Conductors

We present a detailed exposition of how first-principles methods can be used to guide alkali superionic conductor (ASIC) study and design. Using the argyrodite Li 6 PS 5 Cl as a case study, we demonstrate how modern information technology (IT) infrastructure and software tools can facilitate the assessment of alkali superionic conductors in terms of various critical properties of interest such as phase and electrochemical stability and ionic conductivity. The emphasis is on well-documented, reproducible analysis code that can be readily generalized to other material systems and design problems. For our chosen Li 6 PS 5 Cl case study material, we show that Li excess is crucial to enhancing its conductivity by increasing the occupancy of interstitial sites that promote long-range Li + diffusion between cage-like frameworks. The predicted room-temperature conductivities and activation barriers are in reasonably good agreement with experimental values.

42 ENGINEERING↗
DOE OSTI2014

Relativistic electron precipitation events driven by electromagnetic ion-cyclotron waves

We adopt a canonical approach to describe the stochastic motion of relativistic belt electrons and their scattering into the loss cone by nonlinear EMIC waves. The estimated rate of scattering is sufficient to account for the rate and intensity of bursty electron precipitation. This interaction is shown to result in particle scattering into the loss cone, forming ∼10 s microbursts of precipitating electrons. These dynamics can account for the statistical correlations between processes of energization, pitch angle scattering, and relativistic electron precipitation events, that are manifested on large temporal scales of the order of the diffusion time ∼tens of minutes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Four decades of implicit Monte Carlo

In 1971, Fleck and Cummings derived a system of equations to enable robust Monte Carlo simulations of time-dependent, thermal radiative transfer problems. Denoted the “Implicit Monte Carlo” (IMC) equations, their solution remains the de facto standard of high-fidelity radiative transfer simulations. Over the course of 44 years, their numerical properties have become better understood, and accuracy enhancements, novel acceleration methods, and variance reduction techniques have been suggested. In this review, we rederive the IMC equations—explicitly highlighting assumptions as they are made—and outfit the equations with a Monte Carlo interpretation. We put the IMC equations in context with other approximate forms of the radiative transfer equations and present a new demonstration of their equivalence to another well-used linearization solved with deterministic transport methods for frequency-independent problems. We discuss physical and numerical limitations of the IMC equations for asymptotically small time steps, stability characteristics and the potential of maximum principle violations for large time steps, and solution behaviors in an asymptotically thick diffusive limit. We provide a new stability analysis for opacities with general monomial dependence on temperature. Here, we consider spatial accuracy limitations of the IMC equations and discussion acceleration and variance reduction techniques.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Calculation of the RF Conductivity and Hall Constant of a Thin Metal Film

The kinetic problem of finding the rf conductivity and Hall constant of a thin metal film placed in a transverse steady magnetic field and a longitudinal variable electric field has been considered. It has been assumed that electrons diffusely reflect from the upper and lower surfaces of the film. No limitations have been imposed on the relationship between the thickness of the film and the electron free path. The dependences of the conductivity and Hall constant on dimensionless parameters, namely, electric field frequency, magnetic field induction, and thickness of the film, have been studied. Calculation results have been compared with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Calibration of the Diffusivity Predictions of Centipede Using Approximate Bayesian Computation and Applications in Nyx (Engineering Scale) and Xolotl-MARMOT (Meso-Scale) Simulations

Fission gas evolution and release in UO 2 nuclear fuel are important fuel performance metrics and occur in several distinct stages: 1) nucleation, growth and resolution of intra-granular bubbles, 2) diffusion to grain boundaries and 3) nucleation and growth of bubbles at grain boundaries, which eventually form a connected network (percolation) enabling release of gas from grain boundaries through connections to triple junctions, grain edges or free surfaces. The NE-SciDAC project is developing several computational tools to model this problem, which are connected in a hierarchical multi-scale framework. The information transfer in the multi-scale framework is a critical step that, in addition to best-estimates, should include uncertainty quantification. Despite taking a first-principles multi-scale approach, there is a need to perform parameter calibration to ensure consistency with available experimental data. In the present study, uncertainty quantification (UQ) and parameter calibration is demonstrated for one of the lower length scale codes in the multi-scale framework (Centipede) and then the results, including instances of the propagated uncertainties, are used in other codes within the framework, specifically Nyx and Xolotl-MARMOT. We calibrated the model parameters in Centipede, a computer code used to predict diffusivities of uranium (U) and xenon (Xe) in the context of the simulation of fission gas in uranium oxide (UO 2 ) nuclear fuel. The Centipede code depends on 183 parameters, all of which are subject to uncertainty. The three data sets used in our calibration effort are taken from the literature. This data is available as a set of measurements, including measurement errors. Our goal is to calibrate a statistical model that predicts both the value of the measurement and the uncertainty associated with the measurement. We perform a Bayesian calibration of the model parameters using a dedicated approximate Bayesian computation (ABC) likelihood function. To avoid excessive computational costs, we replace the expensive Centipede simulation code by a higher-order surrogate model, constructed using only the 9 most important parameters. These important parameters are identified by a preliminary global sensitivity analysis (GSA) study. Among the important parameters are T0 (the temperature at which UO 2 is perfectly stoichiometric) and Hf_pO2 (the temperature dependence of the oxygen (O) partial pressure) that should be considered as operating conditions to be estimated along with the other parameters. We consider two different cases: one where we define one set of these operating conditions for all data sets, and one where we define distinct operating condition parameters for each data set. The Xe diffusivities predicted by the latter case show distinct features that could not be observed in the former. Next, we use the diffusivity predictions by Centipede as input to Nyx, a reduced order fuel performance code focused on gas behavior alone, in order to estimate quantities associated with inter-granular bubble formation at conditions specified by the experiments. Finally, the diffusivities obtained from the calibrated Centipede runs were used in coupled Xolotl-MARMOT simulations of intra- and inter-granular gas evolution. The results are compared to simulations using the baseline diffusivities from Turnbull et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Grassmannian diffusion maps based surrogate modeling via geometric harmonics

Abstract A novel surrogate model based on the Grassmannian diffusion maps (GDMaps) and utilizing geometric harmonics (GH) is developed for predicting the response of complex physical phenomena. The method utilizes GDMaps to obtain a low‐dimensional representation of the underlying behavior of physical/mathematical systems with respect to uncertain input parameters. Using this representation, GH, an out‐of‐sample extension technique, is employed to create a global map from the input parameter space to a Grassmannian diffusion manifold. GH is further employed to locally map points on the diffusion manifold onto the tangent space of a Grassmann manifold. The exponential map is then used to project the points in the tangent space onto the Grassmann manifold, where reconstruction of the full solution is performed. The performance of the proposed surrogate model is verified with three examples. The first problem is a toy example used to illustrate the technique. In the second example, errors associated with the various mappings are assessed by studying response predictions of the electric potential of a dielectric cylinder in a homogeneous electric field. The last example applies the method for uncertainty prediction in the strain field evolution in a model amorphous material using the shear transformation zone theory of plasticity.

42 ENGINEERING↗
DOE OSTI2022

First-passage time statistics on surfaces of general shape: Surface PDE solvers using Generalized Moving Least Squares (GMLS)

Here, we develop numerical methods for computing statistics of stochastic processes on surfaces of general shape with drift-diffusion dynamics d X t = a (X t ) dt + b(X t ) d W t . We formulate descriptions of Brownian motion and general drift-diffusion processes on surfaces. We consider statistics of the form u (x) = E x [$∫^{τ}_{0}$ g (X t ) dt ] + E x [ f (X τ )] for a domain Ω and the exit stopping time τ = inf t { t >0 | X i Ω}, where f , g are general smooth functions. For computing these statistics, we develop high-order Generalized Moving Least Squares (GMLS) solvers for associated surface PDE boundary-value problems based on Backward- Kolmogorov equations. We focus particularly on the mean First Passage Times (FPTs) given by the case f = 0, g = 1 where u (x) = E x [τ]. We perform studies for a variety of shapes showing our methods converge with high-order accuracy both in capturing the geometry and the surface PDE solutions. We then perform studies showing how statistics are influenced by the surface geometry, drift dynamics, and spatially dependent diffusivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Preliminary Compendium of the Initial Conditions and Evolution of the Fundamental Interactions of SiC with Ag and Pd

Tri-structural isotropic (TRISO) particles are the fuel of choice for very-high temperature reactors. At the core, these spherical particles consist of an inner fuel kernel of UO{sub 2}, UC{sub 2} or a combination of both, which is coated radially outward by successive layers of low density pyrolytic graphite, an inner high density pyrolytic graphite layer, silicon carbide (SiC) and an outer high density pyrolytic graphite layer. SiC is the main diffusion barrier against the release of fission products due to its hardness and high melting point. Irradiation testing of TRISO particles indicates that trace amounts of metallic fission products can damage or potentially escape through the barrier layer. Surrogate experimental studies have previously been performed by various methods seeking to shed light on and explain the diffusion behavior of Ag and Pd. Although some of these experiments have yielded realistic diffusion coefficients the discrepancy of results between different groups is such that the essential details of the diffusive behavior are still to be fully understood. This work aims to understand the fundamental physical and chemical interactions that take place at the interface of SiC and Ag/Pd, and which give rise to the observed diffusion behavior of these metallic species in TRISO particles. Theoretical simulations of the diffusion of Ag through SiC have been exhaustive but typically yield bulk diffusion coefficients that due to the nature of bulk diffusion phenomena (i.e. slow diffusion) the diffusion coefficients are calculated to be orders of magnitude lower than what experimental observations suggest. This in part is due to the difficulty in simulating realistic models that include various polymorphs of SiC, and the dynamical structural transformations that SiC undergoes in the TRISO environment (which are also to be explained). The result is a particularly challenging problem to solve given that the initial conditions of the interactions of Ag and Pd with SiC are difficult to address, both in theoretical models and surrogate experiments. In order to improve these models the approach taken here is to consider fundamental interactions between Ag and SiC that take place in short time scales by varying the initial environmental conditions, i.e. annealing temperature and time, initial amounts of the metallic species (M) Ag, Pd in SiC/M/SiC multilayers. This work aims to understand the fundamental physical and chemical interactions that take place at the interface of SiC and Ag/Pd, and which give rise to the observed diffusion behavior of these metallic species in TRISO particles. In order to do so we have synthesized SiC/Ag/SiC and SiC/Pd/SiC multilayers of various thicknesses by PLD and characterized them using SEM/EDS. Initial results show that upon annealing Ag agglomerates and ruptures the SiC layer to be released. In the case of Pd a strong chemical interaction is deduced from the slow interlayer diffusion and limited release. Further studies using photoelectron spectroscopy and RHEED will be carried out in order obtain detailed information on the chemical and physical transitions that undergo upon annealing of these systems. (authors)

36 MATERIALS SCIENCE↗
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 OSTI2017

Surface tension models for a multi-material ALE code with AMR

A number of surface tension models have been implemented in a 3D multi-physics multi-material code, ALE–AMR, which combines Arbitrary Lagrangian Eulerian (ALE) hydrodynamics with Adaptive Mesh Refinement (AMR). ALE–AMR is unique in its ability to model hot radiating plasmas, cold fragmenting solids, and most recently, the deformation of molten material. The surface tension models implemented include a diffuse interface approach with special numerical techniques to remove parasitic flow and a height function approach in conjunction with a volume-fraction interface reconstruction package. These surface tension models are benchmarked with a variety of test problems. In conclusion, based on the results, the height function approach using volume fractions was chosen to simulate droplet dynamics associated with extreme ultraviolet (EUV) lithography.

42 ENGINEERING↗
DOE OSTI2017

Axial magnetic field injection in magnetized liner inertial fusion

MagLIF is a fusion concept using a Z-pinch implosion to reach thermonuclear fusion. In current experiments, the implosion is driven by the Z-machine using 19 MA of electrical current with a rise time of 100 ns. MagLIF requires an initial axial magnetic field of 30 T to reduce heat losses to the liner wall during compression and to confine alpha particles during fusion burn. This field is generated well before the current ramp starts and needs to penetrate the transmission lines of the pulsed-power generator, as well as the liner itself. Consequently, the axial field rise time must exceed hundreds of microseconds. Any coil capable of being submitted to such a field for that length of time is inevitably bulky. The space required to fit the coil near the liner, increases the inductance of the load. In turn, the total current delivered to the load decreases since the voltage is limited by driver design. Yet, the large amount of current provided by the Z-machine can be used to produce the required 30 T field by tilting the return current posts surrounding the liner, eliminating the need for a separate coil. However, the problem now is the field penetration time, across the liner wall. This paper discusses why skin effect arguments do not hold in the presence of resistivity gradients. Numerical simulations show that fields larger than 30 T can diffuse across the liner wall in less than 60 ns, demonstrating that external coils can be replaced by return current posts with optimal helicity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Role of Ion Kinetic Physics in the Interaction of Magnetic Flux Ropes

To explain many natural magnetized plasma phenomena, it is crucial to understand how rates of collisionless magnetic reconnection scale in large magnetohydrodynamic (MHD) scale systems. Simulations of isolated current sheets conclude such rates are independent of system size and can be reproduced by the Hall-MHD model, but neglect sheet formation and coupling to MHD scales. Here in this study, it is shown for the problem of flux-rope merging, which includes this formation and coupling, that the Hall-MHD model fails to reproduce the kinetic results. The minimum sufficient model must retain ion kinetic effects, which set the ion diffusion region geometry and give time-averaged rates that reduce significantly with system size, leading to different global evolution in large systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Resilient Reinforcement in Secure State Estimation against Sensor Attacks with a priori Information

This paper deals with the state estimation problem in the presence of malicious sensor attacks. Unlike the existing work, in this paper, we address the problem with a priori information of the state to be estimated. Specifically, we address three prior knowledge: the sparsity information, (alpha, n_0)-sparsity information, and side information, and in each scenario, we show that the state can be reconstructed even if more sensors are compromised. This implies that the prior information reinforces the system resilience against sensor attacks. Then, an estimator under sensor attacks considering the information is developed and, under a certain condition, the estimator can be relaxed into a tractable convex optimization problem. Further, we extend this analysis to systems in presence of measurement noises, and it is shown that the prior information reduces the state-estimation error caused by the noise. Here, the numerical simulations in a diffusion process finally illustrate the reinforcement and error-reduction results with the information.

24 POWER TRANSMISSION AND DISTRIBUTION↗