SEARCH · Search NASA

Results for “three space dimension”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Flows in the space of interacting chiral boson theories

We study interacting theories of N left-moving and N ¯ right-moving Floreanini-Jackiw bosons in two dimensions. A parametrized family of such theories is shown to enjoy (nonmanifest) Lorentz invariance if and only if its Lagrangian obeys a flow equation driven by a function of the energy-momentum tensor. We discuss the canonical quantization of such theories along classical stress tensor flows, focusing on the case of the root- T T ¯ deformation, where we obtain perturbative results for the deformed spectrum in a certain large-momentum limit. In the special case N = N ¯ , we consider the quantum effective action for the root- T T ¯ -deformed theory by expanding around a general classical background, and we find that the one-loop contribution vanishes for backgrounds with constant scalar gradients. Our analysis can also be interpreted via dual U ( 1 ) Chern-Simons theories in three dimensions, which might be used to describe deformations of charged AdS 3 black holes or quantum Hall systems. Published by the American Physical Society 2024

Ebert, Stephen (ORCID:0000000154670155)↗

DOE OSTI2023

Microstrata

Abstract Microstrata are the non-extremal analogues of superstrata: they are smooth, non-extremal (non-BPS) solitonic solutions to IIB supergravity whose deep-throat limits approximate black holes. Using perturbation theory and numerical methods, we construct families of solutions using a consistent truncation to three-dimensional supergravity. The most general families presented here involve two continuous parameters, or amplitudes, and four quantized parameters that set the angular momenta and energy levels. Our solutions are asymptotic to the vacuum of the D1-D5 system: AdS 3 ×S 3 × 𝕋 4 . Using holography, we show that the they are dual to multi-particle states in the D1-D5 CFT involving a large number of mutually non-BPS supergravitons and we determine the anomalous dimensions of these states from the binding energies in supergravity. These binding energies are uniformly negative and depend non-linearly on the amplitudes of the states. In one family of solutions, smoothness restricts some of the fields to lie on a special locus of the parameter space. Using precision holography we show that this special locus can be identified with the multi-particle states constructed via the standard OPE of the single-particle constituents. Our numerical analysis shows that microstrata are robust at large amplitudes and the solutions can be obtained to very high precision.

Physics↗

DOE OSTI2020

Three-dimensional electronic scaffolds for monitoring and regulation of multifunctional hybrid tissues

Recently, the integration of electronic elements with cellular scaffolds has brought forth the ability to monitor and control tissue function actively by using flexible free-standing two-dimensional (2D) systems. Capabilities for electrically probing complex, physicochemical and biological three-dimensional (3D) microenvironments demand, however, 3D electronic scaffolds with well-controlled geometries and functional-component distributions. This work presents the development of flexible 3D electronic scaffolds with precisely defined dimensions and microelectrode configurations formed using a process that relies on geometric transformation of 2D precursors by compressive buckling. It demonstrates a capability to fabricate these constructs in diverse 3D architectures and/or electrode distributions aimed at achieving an enhanced level of control and regulation of tissue function relatively to that of other approaches. In addition, this work presents the integration of these 3D electronic scaffolds within engineered 3D cardiac tissues, for monitoring of tissue function, controlling tissue contraction through electrical stimulation, and initiating on-demand, local release of drugs, each through well-defined volumetric spaces. These ideas provide opportunities in fields ranging from in vitro drug development to in vivo tissue repair and many others.

3D electronic scaffolds↗

DOE OSTI2025

CLPNets: Coupled Lie–Poisson neural networks for multi-part Hamiltonian systems with symmetries

To accurately compute data-based prediction of Hamiltonian systems, it is essential to utilize methods that preserve the structure of the equations over time. We consider a particularly challenging case of systems with interacting parts that do not reduce to pure momentum evolution. Such systems are essential in scientific computations, such as discretization of a continuum elastic rod, which can be viewed as the group of rotations and translations $SE(3)$. The evolution involves not only the momenta but also the relative positions and orientations of the particles. The presence of Lie group-valued elements, such as relative positions and orientations, poses a problem for applying previously derived methods for data-based computing. We develop a novel method of data-based computation and complete phase space learning of such systems. We follow the original framework of SympNets (Jin et al., 2020) and LPNets (Eldred et al., 2024), building the neural network from phase space mappings that preserve the Lie–Poisson structure. We derive a novel system of mappings that are built into neural networks describing the evolution of such systems. We call such networks Coupled Lie–Poisson Neural Networks, or CLPNets. We consider increasingly complex examples for the applications of CLPNets, starting with the rotation of two rigid bodies about a common axis, progressing to the free rotation of two rigid bodies, and finally to the evolution of two connected and interacting $SE(3)$ components, describing the discretization of an elastic rod into two elements. Our method preserves all Casimir invariants to machine precision, preserves energy to high accuracy, and shows good resistance to the curse of dimensionality, requiring only a few thousand data points for all cases studied (three to eighteen dimensions). Additionally, the method is highly economical in memory requirements, requiring only about 200 parameters for the most complex case considered.

Data-based modeling↗

DOE OSTI2022

Canted antiferromagnetic phases in the candidate layered Weyl material EuMnSb 2

EuMnSb 2 is a candidate topological material which can be tuned towards a Weyl semimetal, but there are differing reports for its antiferromagnetic (AFM) phases. The coupling of bands dominated by pure Sb layers hosting topological fermions to Mn and Eu magnetic states provides a potential path to tune the topological properties. In this study we present single-crystal neutron diffraction, magnetization, and heat-capacity data as well as polycrystalline 151 Eu Mössbauer data which show that three AFM phases exist as a function of temperature, and we present a detailed analysis of the magnetic structure in each phase. The Mn magnetic sublattice orders into a C-type AFM structure below T N Mn = 323 (1) K with the ordered Mn magnetic moment μ Mn lying perpendicular to the layers. AFM ordering of the Eu sublattice occurs below T N Eu1 = 23 (1) K with the ordered Eu magnetic moment μ Eu canted away from the layer normal and μ Mn retaining its higher temperature order. μ Eu is ferromagnetically aligned within each Eu layer but exhibits a complicated AFM layer stacking. Both of these higher-temperature phases are described by magnetic space group (MSG) $\mathcal{Pn'm'a'}$ with the chemical and magnetic unit cells having the same dimensions. Cooling below T N Eu2 =9 (1) K reveals a third AFM phase where μ Mn remains unchanged but μ Eu develops an additional substantial in-plane canting. This phase has MSG $\mathcal{P11}$ $\frac{2_1}{a′}$. We also find some evidence of short-range magnetic correlations associated with the Eu between 12 K≲ T≲ 30 K . Using the determined magnetic structures, we postulate the signs of nearest-neighbor intralayer and interlayer exchange constants and the magnetic anisotropy within a general Heisenberg model. We then discuss implications of the various AFM states in EuMnSb 2 and their potential for tuning topological properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Theory of criticality for quantum ferroelectric metals

A variety of compounds, for example, doped paraelectrics and polar metals, exhibit both ferroelectricity and correlated electronic phenomena such as low-density superconductivity and anomalous transport. Characterizing such properties is tied to understanding the quantum dynamics of inversion symmetry breaking in the presence of itinerant electrons. Here, we present a comprehensive analysis of the properties of a metal near a quantum critical transition to a ferroelectric state, in both two and three dimensions. In this study, starting from a minimal model of electrons coupled to a transverse polar phonon via a Rashba-type spin-orbit interaction, we compute the dynamical response of both electrons and phonons. We find that the system can evince both Fermi and non-Fermi liquid phases, as well as enhanced pairing in both singlet and triplet channels. Furthermore, we systematically compute corrections to one-loop theory and find a tendency to quantum order-by-disorder, leading to a phase diagram that can include second-order, first-order, and finite-momentum phase transitions. Finally, we show that the entire phase diagram can be controlled via application of external strain, either compressive or volume-preserving. Our results provide a map of the dynamical and thermodynamical phase space of quantum ferroelectic metals, which can serve in characterizing existing materials and in seeking applications for quantum technologies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Gauge-invariant renormalization of four-quark operators

We study the renormalization of four-quark operators in one-loop perturbation theory. We employ a coordinate-space gauge-invariant renormalization scheme (GIRS), which can be advantageous compared to other schemes, especially in nonperturbative lattice investigations. From our perturbative calculations, we extract the conversion factors between GIRS and the modified minimal subtraction scheme ( MS ¯ ) at the next-to-leading order. As a by-product, we also obtain the relevant anomalous dimensions in the GIRS scheme. A formidable issue in the study of the four-quark operators is that operators with different Dirac matrices mix among themselves upon renormalization. We focus on both parity-conserving and parity-violating four-quark operators, which change flavor numbers by two units ( Δ F = 2 ). The extraction of the conversion factors entails the calculation of two-point Green’s functions involving products of two four-quark operators, as well as three-point Green’s functions with one four-quark and two bilinear operators. The significance of our results lies in their potential to refine our understanding of QCD phenomena, offering insights into the precision of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and shedding light on the nonperturbative treatment of complex mixing patterns associated with four-quark operators. Published by the American Physical Society 2024

Constantinou, M. (ORCID:0000000269881745)↗

DOE OSTI2020

Anomalies in the space of coupling constants and their dynamical applications I

It is customary to couple a quantum system to external classical fields. One application is to couple the global symmetries of the system (including the Poincaré symmetry) to background gauge fields (and a metric for the Poincaré symmetry). Failure of gauge invariance of the partition function under gauge transformations of these fields reflects ’t Hooft anomalies. It is also common to view the ordinary (scalar) coupling constants as background fields, i.e. to study the theory when they are spacetime dependent. We will show that the notion of ’t Hooft anomalies can be extended naturally to include these scalar background fields. Just as ordinary ’t Hooft anomalies allow us to deduce dynamical consequences about the phases of the theory and its defects, the same is true for these generalized ’t Hooft anomalies. Specifically, since the coupling constants vary, we can learn that certain phase transitions must be present. We will demonstrate these anomalies and their applications in simple pedagogical examples in one dimension (quantum mechanics) and in some two, three, and four-dimensional quantum field theories. An anomaly is an example of an invertible field theory, which can be described as an object in (generalized) differential cohomology. We give an introduction to this perspective. Also, we use Quillen’s superconnections to derive the anomaly for a free spinor field with variable mass. In a companion paper we will study four-dimensional gauge theories showing how our view unifies and extends many recently obtained results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Dispersion curve regression – assisted wideband local wavenumber analysis for characterizing three-dimensional (3D) profile of hidden corrosion damage

Corrosion as common damage in civil, petrochemical, nuclear, and aerospace structures affects the integrity and safety of the structures and may lead to catastrophic failures. This paper presents dispersion curve regression-assisted local wavenumber analysis method, which can analyze the time–space wavefield containing wideband information of wave-damage interaction and further extract the structural information carried by such wavefield for characterizing hidden corrosion damage in an isotropic plate. To acquire the time–space wavefield, a noncontact scanning laser Doppler vibrometer is used. In our analysis method, the acquired time–space wavefield is firstly processed to generate a local wavenumber function. The derived wavenumber function is further analyzed with a dispersion curve regression step, which searches in a set of theoretical frequency-wavenumber dispersion curves for different plate thicknesses and identifies the optimal plate thickness whose theoretical curve best matches the frequency-wavenumber relation contained in the local wavenumber function. By this means, a 3D profile including both in plane and thickness dimensions of the structure can be constructed for corrosion visualization and quantification. The experimental study demonstrates that our method can quantify the profile of hidden uniform corrosion in metal plates in terms of its in-plane shape and size as well as its out-of-plane depth in the subwavelength scale. Moreover, compared to previous Lamb wave-based corrosion inspection methods, our technique allows for noncontact 3D characterization of hidden corrosion from a far distance to the structure. Finally, we believe this work will inspire new 3D damage quantification methods that are based on wavefield analysis as well as enable potential applications for the quantification of hidden corrosion in civil, petrochemical, nuclear, and aerospace structures.

42 ENGINEERING↗

DOE OSTI2020

Bosonic crystalline symmetry protected topological phases beyond the group cohomology proposal

It is demonstrated by explicit construction that three-dimensional bosonic crystalline symmetry protected topological (cSPT) phases are classified by $H^5_\phi(G; \mathbb{Z})⊕H^1_\phi(G; \mathbb{Z})$ for all 230 space groups $\textit{G}$, where $H^n_\phi(G; \mathbb{Z})$ denotes the $\textit{n}$th twisted group cohomology of $\textit{G}$ with $\mathbb{Z}$ coefficients, and $\phi$ indicates that $\textit{g ε G}$ acts nontrivially on coefficients by sending them to their inverses if $\textit{g}$ reverses space-time orientation and acts trivially otherwise. The previously known summand $H^5_\phi(G;\mathbb{Z})$ corresponds only to crystalline phases built without the $E_8$ state or its multiples on 2-cells of space. It is the crystalline analog of the “group cohomology proposal” for classifying bosonic symmetry protected topological (SPT) phases, which takes the form $H^{d + 2}_\phi(G; \mathbb{Z}) ≅ H^{d + 1}_\phi(G; \text{U}(1))$ for finite internal symmetry groups in d spatial dimensions. The new summand $H^1_\phi(G; \mathbb{Z})$ classifies possible configurations of $E_8$ states on 2-cells that can be used to build crystalline phases beyond the group cohomology proposal. The completeness of our classification and the physical meaning of $H^1_\phi(G; \mathbb{Z})$ are established through a combination of dimensional reduction, surface topological order, and explicit cellular construction. The value of $H^1_\phi(G; \mathbb{Z})$ can be easily read off from the international symbol for $\textit{G}$. In this work, our classification agrees with the prediction of the “generalized cohomology hypothesis,” which concerns the general structure of the classification of SPT phases, and therefore provides strong evidence for the validity of the said hypothesis in the realm of crystalline symmetries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Benchmark of the Kilowatt Reactor Using Stirling TechnologY (KRUSTY) Component Critical Configurations

Kilowatt Reactor Using Stirling TechnologY (KRUSTY) was a prototype for the U.S. National Aeronautics and Space Administration’s Kilopower Program. KRUSTY has a highly enriched uranium–molybdenum alloy (with 7.65 wt% molybdenum) annular core reflected by beryllium oxide with an outer stainless steel shield. Five configurations from the experimental campaign were chosen to be evaluated as benchmark cases. Uncertainties were evaluated in five categories: (1) criticality measurement, (2) mass and density, (3) dimensions, (4) material compositions, and (5) positioning. The largest contribution to the overall uncertainty in each case was from the radial alignment of the movable platen. A simplified model was created to increase computational efficiency, and an average bias of –16 pcm was calculated due to the simplifications. Sample calculations were completed for each case using MCNP6.2, COG, and MC21, all with ENDF/B-VIII.0 nuclear data. For MCNP6.2, the average difference (absolute value) between the calculated and experimental $k_{eff}$ for the five configurations was 14 pcm for both the detailed and the simplified models. The $k_{eff}$ results from all three codes are within 1σ of the benchmark values. KRUSTY’s value as a benchmark is due to its sensitivity to beryllium and molybdenum. For beryllium, KRUSTY adds an 18th benchmark with a total cross-section sensitivity greater than 0.05%/%/(unit lethargy). For molybdenum, KRUSTY adds a 9th benchmark with a total cross-section sensitivity greater than 0.004%/%/(unit lethargy).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

The Virasoro minimal string

We introduce a critical string theory in two dimensions and demonstrate that this theory, viewed as two-dimensional quantum gravity on the worldsheet, is equivalent to a double-scaled matrix integral. The worldsheet theory consists of Liouville CFT with central charge c ≥ ≥ 25 coupled to timelike Liouville CFT with central charge 26-c. The double-scaled matrix integral has as its leading density of states the universal Cardy density of primaries in a two-dimensional CFT, thus motivating the name Virasoro minimal string. The duality holds for any value of the continuous parameter c c and reduces to the JT gravity/matrix integral duality in the large central charge limit. It thus provides a precise stringy realization of JT gravity. The main observables of the Virasoro minimal string are quantum analogues of the Weil-Petersson volumes, which are computed as absolutely convergent integrals of worldsheet CFT correlators over the moduli space of Riemann surfaces. By exploiting a relation of the Virasoro minimal string to three-dimensional gravity and intersection theory on the moduli space of Riemann surfaces, we are able to give a direct derivation of the duality. We provide many checks, such as explicit numerical - and in special cases, analytic - integration of string diagrams, the identification of the CFT boundary conditions with asymptotic boundaries of the two-dimensional spacetime, and the matching between the leading non-perturbative corrections of the worldsheet theory and the matrix integral. As a byproduct, we discover natural conformal boundary conditions for timelike Liouville CFT.

Physics↗

DOE OSTI2021

A machine learning approach for efficient multi-dimensional integration

Many physics problems involve integration in multi-dimensional space whose analytic solution is not available. The integrals can be evaluated using numerical integration methods, but it requires a large computational cost in some cases, so an efficient algorithm plays an important role in solving the physics problems. We propose a novel numerical multi-dimensional integration algorithm using machine learning (ML). After training a ML regression model to mimic a target integrand, the regression model is used to evaluate an approximation of the integral. Then, the difference between the approximation and the true answer is calculated to correct the bias in the approximation of the integral induced by ML prediction errors. Because of the bias correction, the final estimate of the integral is unbiased and has a statistically correct error estimation. Three ML models of multi-layer perceptron, gradient boosting decision tree, and Gaussian process regression algorithms are investigated. The performance of the proposed algorithm is demonstrated on six different families of integrands that typically appear in physics problems at various dimensions and integrand difficulties. The results show that, for the same total number of integrand evaluations, the new algorithm provides integral estimates with more than an order of magnitude smaller uncertainties than those of the VEGAS algorithm in most of the test cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate

Characterizing the properties of groundwater aquifers is essential for predicting aquifer response and managing groundwater resources. In this work, we develop a high-dimensional scalable Bayesian inversion framework governed by a three-dimensional quasi-static linear poroelastic model to characterize lateral permeability variations in groundwater aquifers. In this work, we determine the maximum a posteriori (MAP) point of the posterior permeability distribution from centimeter-level surface deformation measurements obtained from Interferometric Synthetic Aperture Radar (InSAR). The scalability of our method to high parameter dimension is achieved through the use of adjoint-based derivatives, inexact Newton methods to determine the MAP point, and a Mat´ern class sparse prior precision operator. Together, these guarantee that the MAP point is found at a cost, measured in number of forward/adjoint poroelasticity solves, that is independent of the parameter dimension. We apply our methodology to a test case for a municipal well in Mesquite, Nevada, in which InSAR and GPS surface deformation data are available. We solve problems with up to 320,824 state variable degrees of freedom (DOFs) and 16,896 parameter DOFs. A consistent treatment of noise level is employed so that the aquifer characterization result does not depend on the pixel spacing of surface deformation data. Our results show that the use of InSAR data significantly improves characterization of lateral aquifer heterogeneity, and the InSAR-based aquifer characterization recovers complex lateral displacement trends observed by independent daily GPS measurements.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Agentic AI vs ML-Based Autotuning: A Comparative Study for Loop Reordering Optimization

High Performance Computing (HPC) applications rely heavily on code optimizations to achieve good performance on modern CPU and GPU architectures. Traditional Machine Learning auto-tuning approaches have demonstrated success in exploring high-dimensional spaces, but they often require expensive compile-run evaluations and lack adaptability for large HPC applications. The recent advances in Large Language Models (LLMs) and Agentic AI systems raise intriguing questions about the potential of these approaches to address specific optimization methodologies. This work aims to answer an essential question for the HPC community: “How Agentic AI Systems Compare to Traditional ML Autotuning Techniques?” To address this question, we present a comparative analysis between a traditional ML-based optimization approach and an Agentic AI system, evaluating their respective capabilities and limitations for loop-level optimization. In addition, we introduced a new Agentic AI system named LoopGen-AI using three different Large Language Models: GPT-4.1, Claude 4.0, and Gemini 2.5. A key finding is that LoopGen-AI achieves competitive per-formance with only a few program runs, the reasoning logs from the agents revealed that their decisions rely heavily on the combination of semantic understanding of the target kernel with dynamic feedback from the environment, highlighting a promising new dimension in performance tuning. In contrast, ML-based autotuners focus on statistical exploration, and require orders of magnitude more runs to reach peak performance. Additionally, our analysis shows that prompt engineering, particularly using Persona + Context Manager patterns, significantly impacts the effectiveness of Agentic AI. Our results indicate that while Agentic AI systems are not yet a complete replacement for ML-based autotuners, it can effectively complement traditional methods.

Rosas, Miguel Romero↗

DOE OSTI2020

Thermodynamics and magnetism in the two-dimensional to three-dimensional crossover of the Hubbard model

The realization of antiferromagnetic (AF) correlations in ultracold fermionic atoms on an optical lattice is a significant achievement. Experiments have been carried out in one, two, and three dimensions, and have also studied anisotropic configurations with stronger tunneling in some lattice directions. Such anisotropy is relevant to the physics of cuprate superconductors and other strongly correlated materials. Moreover, this anisotropy might be harnessed to enhance AF order. Here we numerically investigate, using the determinant quantum Monte Carlo method, a simple realization of anisotropy in the three-dimensional (3D) Hubbard model in which the tunneling between planes, $t_⊥$, is unequal to the intraplane tunneling $\textit{t}$. This model interpolates between the three-dimensional isotropic ($t_⊥ = t$) and two-dimensional (2D;$t_⊥$ = 0) systems. In this work, we show that at fixed interaction strength to tunneling ratio ($\textit{U/t}$), anisotropy can enhance the magnetic structure factor relative to both 2D and 3D results. However, this enhancement occurs at interaction strengths below those for which the Néel temperature $T_{\text{Néel}}$ is largest, in such a way that the structure factor cannot be made to exceed its value in isotropic 3D systems at the optimal $\textit{U/t}$. We characterize the 2D-3D crossover in terms of the magnetic structure factor, real space spin correlations, number of doubly occupied sites, and thermodynamic observables. An interesting implication of our results stems from the entropy's dependence on anisotropy. As the system evolves from 3D to 2D, the entropy at a fixed temperature increases. Correspondingly, at fixed entropy, the temperature will decrease going from 3D to 2D. This suggests a cooling protocol in which the dimensionality is adiabatically changed from 3D to 2D.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Breaking the mold: Overcoming the time constraints of molecular dynamics on general-purpose hardware

The evolution of molecular dynamics (MD) simulations has been intimately linked to that of computing hardware. For decades following the creation of MD, simulations have improved with computing power along the three principal dimensions of accuracy, atom count (spatial scale), and duration (temporal scale). Since the mid-2000s, computer platforms have, however, failed to provide strong scaling for MD, as scale-out central processing unit (CPU) and graphics processing unit (GPU) platforms that provide substantial increases to spatial scale do not lead to proportional increases in temporal scale. Important scientific problems therefore remained inaccessible to direct simulation, prompting the development of increasingly sophisticated algorithms that present significant complexity, accuracy, and efficiency challenges. While bespoke MD-only hardware solutions have provided a path to longer timescales for specific physical systems, their impact on the broader community has been mitigated by their limited adaptability to new methods and potentials. In this work, we show that a novel computing architecture, the Cerebras wafer scale engine, completely alters the scaling path by delivering unprecedentedly high simulation rates up to 1.144 M steps/s for 200 000 atoms whose interactions are described by an embedded atom method potential. This enables direct simulations of the evolution of materials using general-purpose programmable hardware over millisecond timescales, dramatically increasing the space of direct MD simulations that can be carried out. In this paper, we provide an overview of advances in MD over the last 60 years and present our recent result in the context of historical MD performance trends.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗