Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Comparison of greenhouse gas emission estimates from six hydropower reservoirs using modeling versus field surveys

As with most aquatic ecosystems, reservoirs play an important role in the global carbon (C) cycle and emit greenhouse gases (GHG) as carbon dioxide (CO 2 ) and methane (CH 4 ). However, GHG emissions from reservoirs are poorly quantified, especially in temperate systems, resulting in high uncertainty. We compared reservoir C emission estimates and uncertainty of diffusive, ebullitive, and degassing pathways in six hydropower reservoirs in the southeastern United States among four data sources: two field-based surveys and two models (including the GHG Reservoir “G-res” Tool). We found that CH 4 diffusion was most similar across data sources (modeled minus observed, bias = - 21 g CO 2-eq m -2 y -1 ) and had low relative uncertainty (coefficient of variation, CV = 0.98). On the other hand, CO 2 diffusion was least consistent across data sources (bias = - 518 g CO 2-eq m -2 y -1 ). Both field surveys indicated strong negative CO 2 diffusion (i.e., CO 2 uptake) at all reservoirs, while G-res estimated positive CO 2 diffusion. By extension, total C emissions showed similar discrepancies, leading to high uncertainty in upscaling and interpreting reservoir source-sink dynamics. Finally, CH 4 ebullition had the highest relative uncertainty (CV = 2.77) due to high variability across sites. We discuss limitations of field surveys and these models, including temperature-based annualization methods, varying definitions of ebullition zones, low sampling resolution, and lack of dynamism. Future field efforts focused on capturing variability in CO 2 diffusion and CH 4 ebullition will be especially valuable in reducing uncertainty and improving models to advance our understanding reservoir GHG emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗
DOE OSTI2025

Ion Size Effects on the Thermodynamic, Kinetic, and Mechanical Properties during Ion Exchange in Solid-State Electrolytes

Ion exchange offers a pathway to impose residual compressive stresses to mitigate the electro-chemo-mechanical cracking of solid-state electrolytes such as lithium lanthanum zirconium oxide. This study uses a coupled multiscale framework (integrating density functional theory (DFT), molecular dynamics (MD), and continuum modeling) to examine how exchange ion size influences stress, diffusion, fracture toughness, and electronic properties. Larger isovalent ions (Na + , Ag + , K + ) were exchanged with Li + , with DFT confirming their preference for octahedral 96h sites and a linear relationship between ion size and chemical free expansion coefficient. MD simulations reveal stress and concentration effects on exchange ion diffusivity at elevated temperatures, with Na + and Ag + maintaining favorable mobility while K + showing concentration-dependent clustering. Continuum modeling predicts the range of fracture strength improvements and the required ion exchange concentration profile. It was shown that a 5% surface exchange concentration can induce ∼0.6 GPa of surface compressive stress using Na + and ∼1.0 GPa of surface compressive stress using Ag + . On the other hand, larger ion exchange species may penalize Li + diffusivity by increasing the activation volume and activation energy. Interestingly, Na + has a negligible penalty on Li-ion diffusivity. The room temperature Li + ion diffusivity is reduced by ∼40% with Ag + ion exchange. Electronic band structure analysis shows no size-dependent change in the bandgap, though Ag + introduces localized defect states near the valence band maximum. This study highlights ion size as a key factor in optimizing LLZO properties, offering a framework to improve the solid-state battery performance.

Jagad, Harsh D. [Brown Univ., Providence, RI (Unit↗
DOE OSTI2024

Impact of Variable Perpendicular Transport Coefficients in WEST Simulations Using SolEdge-HDG

Plasma–wall interaction is one of the key research topics on the way to controlled fusion. To study the best operational designs with reduced heat and particle fluxes onto tokamak plasma facing components (PFCs) comprehensive plasma simulations are required. A recent implementation of a hybridized discontinuous Galerkin scheme into a new version of SolEdge code has the advantage of using magnetic equilibrium-free mesh. This allows us to conduct pioneering 2-D transport simulations of a full discharge in the WEST tokamak. In this work, we implemented plasma transport coefficients as functions of coordinate in the poloidal plane and neutral diffusion as a function of neutral mean free path. Moreover, the perpendicular convection flux terms were added to the code. Using the new features, a few test cases were investigated. Finally, the influence of nonconstant transport coefficients on the simulated particle and heat fluxes onto the WEST tokamak PFCs are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modularization of Ceramic Hollow Fiber Membrane Technology for Air Separation

This proposed project is aimed at studying high performance and economically competitive ceramic membrane technology for air separation and high-purity oxygen production using hollow fiber ceramic membrane stack and module technology. The design of the single permeate membrane is a hollow fiber substrate-supported thin tri-layer structure. Radially well-aligned micro-channels are embedded in the thick substrate and open at the inner surface of the substrate, enabling facile air/gas diffusion. The tri-layer structure of thin dense membrane layer (~ 10 µm) sandwiched by porous surface layer on either side is then built on the substrate using advanced fabrication process. The tri-layer structure design allows different materials to be used in different layers, where the materials of surface layers have high surface exchange coefficients while the material of dense layer has high bulk diffusivity. Such a synergetic combination leads to high permeation performance. The focus of the proposed project will be on membrane stack/module development using the developed single hollow fiber membranes, including: 1) fabrication and characterization of novel single hollow fiber membranes; 2) membrane stack design and assembly using fabricated single membranes; 3) stack modeling and analysis to guide membrane stack designs; 4) permeation performance testing and characterization of membrane stacks. The hollow fiber feature and simple sealing requirement enable very compact design of membrane stacks, which have excellent flexibility for further modularizations at different scales. The operations of such membrane stack and module may employ the exhaust heat from other components of Integrated Gasification Combined Cycle and oxy-combustion systems. Therefore, modularization of such an air separation membrane technology can be incorporated into the DOE’s REMS (radically engineered modular systems)-gasification skid and support the oxidant feed of an oxygen-blown REMS gasifier scaled to different ranges.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Hydrogen transport in yttrium hydride under asymmetric heat

Metal hydrides are a promising moderator material for high temperature fission reactors. Yttrium hydride can be loaded to a high hydrogen density with relatively high hydrogen stability at temperatures up to 800°C. This makes yttrium hydride a potential moderator material for microreactors as foreseen in the fission surface power program. However, during the operation of such advanced reactors temperature gradients are expected which can change the local hydrogen density in the moderator. Hydrogen diffusion in metals is driven by a concentration gradient (Fick’s law) and thermal diffusion (Soret diffusion). Thermal diffusion is the transport of hydrogen, or other species, due to a temperature gradient. For example, hydrogen might migrate from the hot side of a sample to the cold side of a sample. Measuring Fickian diffusion is achieved through various permeation or absorption experiments, however measuring thermal diffusion is challenging and has rarely been performed. The Hydrogen Experimental Apparatus for Thermal Diffusion (HEATD) experiment is designed to induce thermal diffusion in samples and quench those samples so that the hydrogen distribution can be analyzed using hot vacuum extraction (HVE). One side of the sample was heated to a high temperature e.g., 800°C, while the other side of the sample is at a lower temperature. The sample was held under the applied temperature gradient for a given time until the anticipated hydrogen diffusion has occurred. The actual time depends depend on the sample composition and hydrogen concentration. The results from HVE showed that thermal diffusion took place in the specimen and the Soret coefficient was calculated.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unraveling Interdiffusion Phenomena and the Role of Nanoscale Diffusion Barriers in the Copper–Gold System

Diffusion is one of the most fundamental concepts in materials science, playing a pivotal role in materials synthesis, forming, and degradation. Of particular importance is solid state interdiffusion of metals which defines the usable parameter space for material combinations in the form of alloys. This parameter space can be explored on the macroscopic scale by using diffusion couples. However, this method reaches its limit when going to low temperatures, small scales, and when testing ultrathin diffusion barriers. Therefore, this work transfers the principle of the diffusion couples to small scales by using core–shell nanowires and in situ heating. This allows us to delve into the interdiffusion dynamics of copper and gold, revealing the interplay between diffusion and the disorder–order phase transition. Our in situ TEM experiments in combination with chemical mapping reveal the interdiffusion coefficients of Cu and Au at low temperatures and highlight the impact of ordering processes on the diffusion behavior. The formation of ordered domains within the solid-solution is examined using high-resolution imaging and nanodiffraction including strain mapping. In addition, we examine the effectiveness of ultrathin Al 2 O 3 barrier layers to control interdiffusion of the diffusion couple. Our findings indicate that a 5 nm thick layer serves as an efficient diffusion barrier. Furthermore, this research provides valuable insights into the interdiffusion behavior of Cu and Au on the nanoscale, offering potential applications in the development of miniaturized integrated circuits and nanodevices.

alloys↗
DOE OSTI2026

Variations in lithium vapor cave performance predictions due to radial transport and recycling assumptions

The lithium vapor cave is a detached divertor design that uses a single private flux region baffle to contain a dense cloud of lithium vapor to dissipate heat flux. Plasma flows are created via fuel gas puffing in order to minimize lithium contamination of the main plasma. Significant modeling using the 2D edge code SOLPS-ITER has already been performed, predicting that sufficient target heat flux reductions (q$^{max}_{Target}$ MW m −2 ) with ( $n$ Li /$n$ e ) LCFS < 0.05 is possible in a case with 90 MW m −2 unmitigated heat flux in NSTX-U. Low heat flux and low upstream concentration was found with a variety of combinations of divertor geometries, target recycling coefficients, upstream plasma parameters, lithium evaporation locations and deuterium fueling locations, with variations in performance found for each design choice. However, the most universal uncertainty of SOLPS-ITER simulations has until now remained unaddressed systematically, namely the cross-field anomalous particle and heat diffusivities. This article aims to bound the uncertainty in lithium concentration prediction as a result of the assumed deuterium cross-field transport and recycling. For the simulations presented here, a factor of 2.1 increase in the upstream lithium density prediction across a factor of four decrease to the assumed deuterium radial particle diffusivity is found. This result is compared across different assumed deuterium recycling coefficients, known to reduce with lithium injection. Upstream lithium density is found to vary by a factor of 2.4 across a feasible range of recycling coefficients at similar $n$$^{OMP,sep}_{e}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗
DOE OSTI2025

Measurement of turbulent velocity and bounds for thermal diffusivity in laser shock compressed foams by x-ray photon correlation spectroscopy

Experimental benchmarking of transport coefficients under extreme conditions is required for validation of differing theoretical models. To date, measurement of transport properties of dynamically compressed samples remains a challenge with only a limited number of studies able to quantify transport in high pressure and temperature matter. x-ray photon correlation spectroscopy utilizes coherent x-ray sources to measure time correlations of density fluctuations, thus providing measurements of length and timescale-dependent transport properties. Here, we present a first-of-a-kind experiment to conduct x-ray photon correlation spectroscopy in laser shock compression experiments. We report measurement of the turbulent velocity in the wake of a laser driven supersonic shock and place an upper bound on thermal diffusivity in a solid density plasma on nanosecond timescales.

High-energy-density plasmas↗
DOE OSTI2025

Optimizing energy yield of monolithic perovskite/silicon tandem solar cells in real-world Conditions: The impact of luminescent coupling

Efficient light management is key to maximizing power conversion efficiency (PCE) in monolithic perovskite/silicon tandem solar cells. Achieving peak efficiency requires closely matched current generation in all junctions, especially in integrated configurations. However, real-world conditions vary significantly due to factors such as sunlight spectrum, diffuse-to-direct sunlight ratio, angular distribution of light, subcell temperature coefficients, and ground reflection. This study introduces a comprehensive optical and device simulation to optimize perovskite/silicon tandem cells, considering experimental luminescent coupling (LC) efficiency and its dependence on working conditions, alongside variations in radiative recombination, effect of temperature on absorptivity spectra, and cloud cover. Our results show potential energy yield improvements of up to 1.4 % with LC, based on current perovskite radiative recombination records, and up to 4 % with direct bandgap materials. Although radiative recombination's dependence on excitation intensity reduces output power and requires thicker absorbers, LC compensates for these losses. LC also lowers the optimized bandgap for the perovskite top cell from 1.72 eV to 1.64–1.68 eV, or even lower in regions with redshifted irradiance. Additionally, optimization revealed that thinner silicon bottom cells require a higher perovskite top cell bandgap, impacting the balance between fabrication cost and cell stability.

14 SOLAR ENERGY↗
DOE OSTI2024

Multinuclear PFGSTE NMR description of 39 K, 23 Na, 7 Li, and 1 H specific activation energies governing diffusion in alkali nitrite solutions

While pulsed field gradient stimulated echo nuclear magnetic resonance (PFGSTE NMR) spectroscopy has found widespread use in the quantification of self-diffusivity for many NMR-active nuclei, extending this technique to uncommon nuclei with unfavorable NMR properties remains an active area of research. Potassium-39 ( 39 K) is an archetypical NMR nucleus exhibiting an unfavorable gyromagnetic ratio combined with a very low Larmor frequency. Despite these unfavorable properties, this work demonstrates that 39 K PFGSTE NMR experiments are possible in aqueous solutions of concentrated potassium nitrite. Further, analysis of the results indicates that 39 K NMR diffusometry is feasible when the nuclei exhibit spin–lattice and spin–spin relaxation coefficients on the order of 60–100 ms and 50–100 ms, respectively. The diffusivity of 39 K followed Arrhenius behavior, and comparative 23 Na, 7 Li, and 1 H PFGSTE NMR studies of equimolal sodium nitrite and lithium nitrite solutions led to correlations between the enthalpy of hydration with the activation energy governing self-diffusion of the cations and also of water. Realizing the feasibility of 39 K PFGSTE NMR spectroscopy has a widespread impact across energy sciences because potassium is a common alkali element in energy storage materials and other applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗
DOE OSTI2025

Computational Design of Cost-Effective High-Entropy Thermal/Environmental Barrier Coatings

Developing cost-effective thermal/environmental barrier coatings (TEBC) requires balance among various properties including low thermal conductivity, matching coefficient of thermal expansion (CTE), high thermal stability, high fracture toughness, and high recession resistance while being affordable. Low oxygen diffusivity is desirable as it can slow down oxygen transport to reach the underlying bond coating and hence delay oxidation of the bond coating. This project aims to design low-cost high-performance TEBC based on high entropy rare earth disilicates to protect SiC-based ceramic matrix composites from chemical and thermal attack for better performance of components in the hot section of gas turbine engines. To accelerate the TEBC design, we utilize first-principles density functional theory to predict key properties including phase stability, CTE, lattice thermal conductivity, temperature-dependent elastic constants, and oxygen diffusivity. Alloying elements including Yb, Y, Er, Eu, Gd, Lu, La, and Ce are considered, and modeling prediction are compared with available experimental results.

coefficient of thermal expansion↗
DOE OSTI2025

Flux Hypothesis for Odd Transport Phenomena

Onsager's regression hypothesis makes a fundamental connection between macroscopic transport phenomena and the average relaxation of spontaneous microscopic fluctuations. This relaxation, however, is agnostic to odd transport phenomena, in which fluxes run orthogonal to the gradients driving them. To account for odd transport, we generalize the regression hypothesis, postulating that macroscopic linear constitutive laws are, on average, obeyed by microscopic fluctuations, whether they contribute to relaxation or not. From this "flux hypothesis," Green-Kubo and reciprocal relations follow, elucidating the separate roles of broken time-reversal and parity symmetries underlying various odd transport coefficients. As an application, we derive and verify the Green-Kubo relation for odd collective diffusion in chiral active matter, first in an analytically tractable model and subsequently through molecular dynamics simulations of concentrated active spinners.

Hargus, Cory↗
DOE OSTI2026

Statistical inference of anomalous thermal transport with uncertainty quantification for interpretive 2D SOL models

The critical task of inferring anomalous cross-field transport coefficients is addressed in simulations of boundary plasmas with fluid models. A workflow for parameter inference in the UEDGE fluid code is developed using Bayesian optimization with parallelized sampling and integrated uncertainty quantification. In this workflow, transport coefficients are inferred by maximizing their posterior probability distribution, which is generally multidimensional and non-Gaussian. Uncertainty quantification is integrated throughout the optimization within the Bayesian framework that combines diagnostic uncertainties and model limitations. As a concrete example, we infer the anomalous electron thermal diffusivity $\chi_\perp$ from an interpretive 2D model describing electron heat transport in the conduction-limited region with radiative power loss. The workflow is first benchmarked against synthetic data and then tested on H-, L-, and I-mode discharges to match their midplane temperature and divertor heat flux profiles. We demonstrate that the workflow efficiently infers diffusivity and its associated uncertainty, generating 2D profiles that match 1D measurements. Future efforts will focus on incorporating more complicated fluid models and analyzing transport coefficients inferred from a large database of experimental results.

Bayesian optimization↗
DOE OSTI2024

A high-order finite difference method for moving immersed domain boundaries and material interfaces

Here, we present a high-order sharp treatment of immersed moving domain boundaries and material interfaces, and apply it to the advection-diffusion equation in two and three dimensions. The spatial discretization combines dimension-split finite difference schemes with an immersed boundary treatment based on a weighted least-squares reconstruction of the solution, providing stable discretizations with up to sixth order accuracy for diffusion terms and third order accuracy for advection terms. The temporal discretization relies on a novel strategy for maintaining high-order temporal accuracy in problems with moving boundaries that minimizes implementation complexity and allows arbitrary explicit or diagonally-implicit Runge-Kutta schemes. The approach is broadly compatible with popular PDE-specialized Runge-Kutta time integrators, including low-storage, strong stability preserving, and diagonally implicit schemes. Through numerical experiments we demonstrate that the full discretization maintains high-order spatial and temporal accuracy in the presence of complex 3D geometries and for a range of boundary conditions, including Dirichlet, Neumann, and flux conditions with large jumps in coefficients.

97 MATHEMATICS AND COMPUTING↗