SEARCH · Engineering Papers

Results for “Diffusion coefficient”

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.

SourceResearch connectionFrom yearThrough year

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

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Molecular Modeling and Molecular Dynamics Simulation of a Packed and Intact Bacterial Microcompartment

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria which encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semipermeable shell proteins and are packed full of enzyme cargoes and metabolites as they fulfill their function. Coupling together recent SAXS and proteomics work, it is possible to develop molecular models for these microcompartments and interrogate enzyme and metabolite dynamics within. Our primary goal of this study is to quantify the permeability of metabolite glyceraldehyde-3-phosphate (G3P) and dihydroxyacetone phosphate (DHAP) across the BMC shell through classical molecular dynamics simulation. The Haliangium ochraceum model of BMC shell (PDB: 6MZX) was used to model an intact BMC of approximately 10 million atoms. Working at this scale presented its own challenges in managing large data sets, with multiple challenges and hardware advances discussed that facilitated this work. Over approximately 750 ns of aggregate simulation, we see multiple permeation events for these metabolites that were added at high concentration through the pores present within BMC shell tiles. When compared to independent permeability estimates for the same metabolites determined through replica exchange umbrella sampling simulations, the permeabilities varied by approximately 3 orders of magnitude. Regardless, the permeability coefficients for both G3P and DHAP are highly similar and very high, such that only very small concentration gradients can be maintained across the BMC shell between the cytosol and BMC interior. The large simulation systems also facilitated comparisons for molecular diffusivity in the crowded environment within the BMC shell. By our estimates, the viscosity within a packed BMC shell is at least 10-fold higher than it would be in neat solution and is the real driver for varying permeability estimates we obtained through simulation. These findings will be used as design inputs for future bioengineering efforts to make products from BMCs, highlighting how permeable BMC shells can be.

Diffusion↗

DOE OSTI2025

Epoxy-based vitrimeric semi-interpenetrating network/MXene nanocomposites for hydrogen gas barrier applications

Herein, we report MXene-filled epoxy-based vitrimeric nanocomposites featuring a semi-interpenetrating network (S-IPN) to develop a hydrogen gas (H 2 ) barrier coating with self-healing characteristics for compressed H 2 storage applications. The reversible epoxy network was formed by synthesizing linear epoxy chains with pendent bis-hydroxyl groups using amino diol, which were then crosslinked with 1,4-benzenediboronic acid to generate dynamic boronic ester linkages. To achieve the S-IPN-type molecular arrangement, the epoxy chains were in situ crosslinked in the presence of poly(ethylene-co-vinyl alcohol) (EVOH), giving rise to a self-healing network (EEP) with a healing efficiency of 87%. Into the S-IPN vitrimer (EEP), a 2D platelet-type nanofiller MXene was incorporated to introduce a tortuous path for H 2 gas diffusion along with improved mechanical properties. The nanocomposite coating was applied to nylon 6 liner material, which is conventionally used in all-composite H 2 storage vessels. Further, the application of a 2 wt% MXene/EEP nanocomposite coating showed a permeability coefficient of 0.062 cm 3 mm m -2 d -1 atm -1 exhibiting ~96% reduction in gas permeability compared to uncoated nylon 6. The same nanocomposite exhibited a healing efficiency of 79%. Increasing the MXene loading to 10 wt% further reduced the permeability coefficient to 0.002 cm 3 mm m -2 d -1 atm -1 ; however, the healing efficiency decreased due to restricted chain mobility. In essence, the current work highlights the potential of vitrimeric S-IPN nanocomposite coatings for H 2 gas-barrier applications, enhancing safety and performance.

42 ENGINEERING↗

DOE OSTI2020

Comparative analysis of the methods for quantitative determination of water content in skin from diffuse reflectance spectroscopy data

The predictive properties of methods aimed for estimating the water content in skin from the spectral diffuse reflection characteristics near the water absorption line in the near-IR spectral range are analysed. Numerical simulation data, experimental data on diffuse reflection from human skin phantoms, and data from the reference data set of human skin reflectance spectra are used to consider the possibility of gaining additional information about the water distribution in skin. The influence of variations in the scattering coefficient and oxyhaemoglobin concentration on the water content estimates is investigated. (paper)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Charm and bottom hadrons in hot hadronic matter

Heavy quarks, and the hadrons containing them, are excellent probes of the QCD medium formed in high-energy heavy-ion collisions, as they provide essential information on the transport properties of the medium and how quarks color-neutralize into hadrons. Large theoretical and phenomenological efforts have been dedicated thus far to assess the diffusion of charm and bottom quarks in the quark–gluon plasma and their subsequent hadronization into heavy-flavor (HF) hadrons. However, the fireball formed in heavy-ion collisions also features an extended hadronic phase, and therefore any quantitative analysis of experimental observables needs to account for the rescattering of charm and bottom hadrons. This is further reinforced by the presence of a QCD cross-over transition and the notion that the interaction strength is maximal in the vicinity of the pseudo-critical temperature. We review existing approaches for evaluating the interactions of open HF hadrons in a hadronic heat bath and the pertinent results for scattering amplitudes, spectral functions and transport coefficients. While most of the work to date has focused on D -mesons, we also discuss excited states as well as HF baryons and the bottom sector. Both the HF hadro-chemistry and bottom observables will play a key role in future experimental measurements. We also conduct a survey of transport calculations in heavy-ion collisions that have included effects of hadronic HF diffusion and assess its impact on various observables.

effective field theories↗

DOE OSTI2020

Gas Transport in the Insensitive High Explosive PBX 9502

Small-scale laboratory experiments were performed to analyze gas transport mechanics in PBX 9502, an insensitive high explosive (IHE) composition. Two independent investigations are reported here. First, gas permeametry and dynamic pycnometry techniques were used to measure the molecular flow (Knudsen) coefficient and the internally interconnected void fraction of pristine and thermally damaged samples in two uniaxial pressing orientations. The permeability of PBX 9502 was found to be unmeasurably low (of the same order of magnitude as PTFE) with gas transport being diffusion-dominated. Secondly, a pressure vessel experiment was developed to measure quasi-static and dynamic gas generation as the explosive was heated to self-ignition (cookoff). The gas generation results and the permeability/diffusivity findings provide evidence that PBX 9502 remains impermeable until seconds prior to self-ignition. At ignition, internal void-pressure drives macro-scale cracking and the sample becomes uniformly incorporated in the following deconsolidating deflagration. These results are discussed within the context of previous observations of pressure-dependent cookoff behavior of PBX 9502 and provide a more complete description of thermal damage evolution in this explosive composition.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2020

Adaptive Multirate Mass Transfer (aMMT) Model: A New Approach to Upscale Regional‐Scale Transport Under Transient Flow Conditions

Abstract The long‐term evaluation of regional‐scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady‐state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective‐dispersive contaminant transport simulated in a 3‐D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady‐state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection‐dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion‐dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Guo, Zhilin↗

DOE OSTI2022

High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces

Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Thermodynamic and Kinetic Modulation of Methylammonium Lead Bromide Crystallization Revealed by In Situ Monitoring

Hybrid organic–inorganic perovskite (HOIP) crystals are promising optoelectronic materials, but little is known about either the thermodynamic and kinetic controls on crystal growth or the underlying growth mechanism(s). Herein, we use fluid cell atomic force microscopy (AFM) and solution nuclear magnetic resonance (NMR) spectroscopy to investigate growth of the model HOIP crystal CH 3 NH 3 PbBr 3 (MAPbBr 3 ) and to determine how formic acid (HCOOH) modulates the thermodynamics and kinetics of growth. The results show that growth of MAPbBr 3 in dimethylformamide (DMF) proceeds through the classical pathway by the spreading of atomic crystal steps generated at screw dislocations on the {100} surface. Temperature dependent step velocity measurements demonstrate that with increasing concentration, HCOOH decreases both the solubility of MAPbBr 3 and the kinetic coefficient (b) of step movement. 1 H-NMR measurements indicate that HCOOH increases the lifetime of the methylammonium (MA + ) ions and promotes the association of MAPbBr 3 , thus tuning the solubility of the perovskite. HCOOH also alters the molecular tumbling motion and bulk diffusion of the MA + ions, possibly via H-bonding. Further, our findings establish a direct correlation between the mesoscale crystal growth kinetics and the molecular-scale interactions between organic additives and constituent ions, providing unprecedented insights for developing predictive syntheses of HOIP crystals with defined size, crystal habit and shape, and defect distribution.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Sectoral and growth rate control on elemental uptake by individual calcite crystals

Heterogeneous distribution of trace elements (impurities) within individual calcite crystals is a phenomenon commonly observed in natural and laboratory systems. Changes in thermodynamic intensive parameters (mostly chemical potential and temperature) cannot always explain the inhomogeneous impurity patterns in calcite crystals and it has been suggested that growth rate and crystallographic orientation may exert strong effects on elemental incorporation into calcite. In addition, there are a number of experimental studies on micro-scale element (E) distribution between non-equivalent pairs of calcite vicinal faces (known as sectoral zoning); however, the variability of partition coefficients (e.g., K E = (E/Ca) calcite /(E/Ca) fluid ) within individual crystals remains undetermined. In this study, we have extended the work on elemental distribution between crystal sectors to evaluation of partition coefficients of trace and minor elements (Li, B, Mg, and Sr) in calcite crystal faces (10–14) and (01-12), whose growth rates were assessed. Growth entrapment model (GEM) and lattice strain theory were applied to explain K E heterogeneity by varying near-surface diffusivity of Mg and Sr and by varying surface enrichment factor for Li, B, Mg, and Sr. Decoupling of sectoral and growth rate effects reveals that sectoral zoning plays a key role in elemental distribution. More specifically, K Li and K B vary by more than one order of magnitude and K Mg varies by a factor of two within individual crystal faces. Strontium sectoral distribution is different from those of Li, B, and Mg and K Sr varies by up to a factor of two. Finally, these behaviors likely reflect different mechanisms of incorporation of Li, B, Mg, and Sr into calcite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Compatibility of LaFe13−x−yMnxSiyH1.6 and Eutectic Liquid GaInSn Alloy

The heat transfer rate of magnetocaloric regenerators is a topic of extensive research and the cyclability of these regenerators is critical to the operation of systems with a high coefficient of performance (e.g., potentially >22, significantly higher than typical vapor compression cooling technologies). To enable a high operating frequency that will result in a high specific cooling power, the heat transfer fluid should have high thermal conductivity and lower specific heat, i.e., higher thermal diffusivity. Eutectic metal alloys possess these qualities, such as gallium–indium–tin (Galinstan), whose thermal diffusivity has been found to be approximately an order of magnitude higher than water. For this study, the effects of eutectic liquid Galinstan exposure on the phase stability of LaFe13−x−yMnxSiyH1.6 magnetocaloric powders in an active magnetic regenerator device were investigated. The powders were characterized before and after exposure to Galinstan using X-ray diffraction, in which the phases were determined using the Rietveld refinement technique and X-ray fluorescence. It was found that after Galinstan exposure, hydrogen containing phases were present in the powder, suggesting that the hydrogen was lost from the magnetocaloric phase. The magnetocaloric phase degradation indicates that the powder was incompatible with the Galinstan metal in an environment with moisture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Viscous flow properties and hydrodynamic diameter of phenothiazine-based redox-active molecules in different supporting salt environments

We report viscous flow properties of a redox-active organic molecule, N-(2-(2-methoxyethoxy)ethyl)phenothiazine (MEEPT), a candidate for non-aqueous redox flow batteries, and two of its radical cation salts. A microfluidic viscometer enabled the use of small sample volumes in determining viscosity as a function of shear rate and concentration in the non-aqueous solvent, acetonitrile, both with and without supporting salts. All solutions tested show Newtonian behavior over shear rates of up to 30 000 s−1, which was rationalized by scaling arguments for the diffusion-based relaxation time of a single MEEPT molecule without aggregation. Neat MEEPT is flowable but with a large viscosity (412 mPa⋅s at room temperature), which is ∼1000 times larger than that of acetonitrile. MEEPT solutions in acetonitrile have low viscosities; at concentrations up to 0.5 M, the viscosity increases by less than a factor of two. From concentration-dependent viscosity measurements, molecular information was inferred from intrinsic viscosity (hydrodynamic diameter) and the Huggins coefficient (interactions). Model fit credibility was assessed using the Bayesian Information Criterion. It is found that the MEEPT and its charged cations are “flowable” and do not flocculate at concentrations up to 0.5 M. MEEPT has a hydrodynamic diameter of around 8.5 Å, which is almost insensitive to supporting salt and state of charge. This size is comparable to molecular dimensions of single molecules obtained from optimized structures using density functional theory calculations. The results suggest that MEEPT is a promising candidate for redox flow batteries in terms of its viscous flow properties.

Wang, Yilin (ORCID:0000000316102995)↗

DOE OSTI2021

Simulation of the SPARC plasma boundary with the UEDGE code

Here in this work, the UEDGE edge transport code is used to examine conditions in the SPARC divertor and edge plasma for various levels of carbon impurity and power from the core (P SOL ). A double-null magnetic configuration is simulated assuming up-down symmetry in geometry and physics. The anomalous heat and particle transport coefficients are tuned to match empirical predictions for SPARC's midplane density profiles, target plate heat flux profiles, and inner/outer divertor power sharing. Convective transport is included on the low-field side, while on the high-field side the transport is modeled as purely diffusive. Hydrogen neutrals are modeled as a fluid with inertial effects, and a carbon impurity is included using the fixed-fraction model. We find that detachment induced by impurity seeding could significantly reduce the heat flux to the divertor surfaces in the SPARC tokamak. At P SOL = 28 MW (the value predicted for SPARC's full-power H-mode scenario) cases with both divertor legs detached were obtained with a carbon impurity fraction between 0.3%–1.4%, far below Z eff limits for SPARC. When the plasma in the outer leg is detached, the peak heat flux density perpendicular to the target plate is below 1 MW m -2 , electron and ion temperatures are less than 1.5 eV, and momentum detachment is observed. However, the detachment state is found to be sensitive to the side-wall boundary conditions, the level of neutral pumping, and the target plate tilt. Finally, a broadly similar SOLPS simulation of SPARC is used to assess the appropriateness of the simpler impurity and neutral models used in UEDGE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

MELCOR Deposition Models Applied to Noble Metals

Noble metals are chemically inert elements that form particles that can circulate around the fuel loop of a liquid-fueled molten salt reactor (LF-MSR). Due to their health-hazardous abilities, it is important to track these elements and to understand their contributions to source term during a reaction, particularly for Te as it decays into iodine. This report describes three solid deposition models that were implemented into the MELCOR reactor modeling code. The models that were implemented were Dittus-Boelter for turbulent fluid deposition, Talbot thermophoretic transport equation and Richardson-Zaki gravitational deposition. These closures were used to calculate the mass transfer coefficient and was compared to those reported for the Molten Salt Reactor Experiment (MSRE). The impacts of salt, operation and particle parameters on how they affect solid deposition are described, and found that solid deposition is strongly dependent on salt flow and viscosity, as well as particle diffusion rate and size. The conclusions form the basis for activities to be carried out in FY26 and for pursuit of future activities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

The design of a photoionization front experiment using the Z-Machine as a driving source and estimated measurements

Radiation-driven heat fronts are present in the early universe during reionization, the circumstellar medium of supernovae, and in high-energy-density physics experiments. Dedicated experiments to observe and diagnose the behavior of these types of heat fronts can improve our understanding of these phenomena. A simulation study of photoionization fronts using the HELIOS-CR radiation hydrodynamics code provides an experimental design for the Z-Machine at Sandia National Laboratory using a measurement-calibrated input radiation flux to drive the photoionization front. The simulations use detailed atomic physics and non-diffusive radiation transport in 1D to determine an optimal gas pressure of 0.75 atm for an experiment in N gas as well as the effects of increasing the thickness of the window that seals the gas cell. Post-processing of these simulations demonstrates that ratios of atomic rate coefficients place the heat front in a physics regime where photoionization dominates the energy deposition. Here, to see the sensitivity of the simulations to changes in the model and spatial grid, this analysis performed resolution, atomic model detail, and radiation transport angular grid studies showing less than 10% deviation from the nominal model for increased complexity, when possible. An effort to emulate 3D geometric effects on the radiation flux using an artificial attenuation scheme has shown that, even for conservative estimates of the flux, simulations still produce a photoionization front. Estimations of a streaked, visible spectroscopy measurement using SPECT3D showed that line emission measurements are present early in time and that later in time thermal emission should become dominant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Transport Processes in a Li-ion Cell during an Internal Short-Circuit

Internal short-circuit in a lithium-ion cell causes an abrupt increase in cell temperature and triggers subsequent thermal runaway. In this work, we present a detailed electrochemical-thermal model to investigate the physical behavior during an internal short-circuit. Simulations at wide range of heat transfer coefficients and short-circuit resistances are conducted to illustrate electrochemical and thermal behavior under a wide range of conditions. The Joule heating at the shorted region promotes electrochemical reactions nearby, causing in-plane non-uniformity of electrolyte and active material transport. Furthermore, it is found that diffusion in solid active materials plays a significant role at very high shorting currents (~20 C), as electrochemical reactions rate are being controlled progressively more by availability of Li + at the interface, due to limitations in diffusion through the active material with increasing discharge rates. This diffusion limitation causes a drop in available energy, and subsequently a decrease in cell equilibrium potential and the heat generation rate at the location of the short. On the other hand, rapid depletion of lithium concentration in the electrolyte and accumulation on the electrode surface results in highly non-uniform transport properties resulting in higher heat generation rates. Hence, the heating regime shifts from "local heating" to "global heating". Based on the findings, important design parameters for battery safety are discussed.

25 ENERGY STORAGE↗