Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Simulation”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

ASCR Workshop Position Paper: Challenges and Opportunities in High Energy Physics

High energy particle physics and cosmology concern themselves with estimating fundamental parameters of nature, such as the masses and interactions of fundamental particles like the Higgs boson and the rate of expansion of the universe. In doing so, they analyze exabyte-scale datasets, some of the largest in all of science, and face many challenges in subsequent data analysis. These challenges are shared between the two disciplines, but we focus on particle physics to highlight one specific domain. In particle physics, the standard method for estimating parameters involves performing Monte Carlo (MC) integration as a function of both parameters of interest and nuisance parameters using an expensive simulator, counting the number of observed collision events (i.i.d. samples) from an experiment in the corresponding integration domains, and forming a Poisson likelihood function. This likelihood function is then used in a Frequentist manner to construct a maximum likelihood point estimate (MLE) and confidence set for the parameters. To sufficiently populate the high-dimensional integration domains, simulators consume billions of CPU-hours annually and produce hundreds of petabytes of intermediate output data. Several techniques have been developed to: optimize definitions of the integration domains so as to be maximally sensitive to a particular subset of parameters, efficiently estimate the integrals, and build robust surrogate models by interpolating between integral evaluations at different parameter points. One can view this whole endeavor as classical Simulation-Based Inference (SBI).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Economic and Environmental Assessment of Biological Conversions of Agile BioFoundry (ABF) Bio-Derived Chemicals

Bio-derived chemicals are an essential part of the growing bioeconomy. They possess the potential to boost a new domestic bioproduct industry, improve the sustainability of integrated biomanufacturing, and reduce U.S. dependence on fossil energy. The Agile BioFoundry (ABF) consortium, a Department of Energy (DOE)-sponsored collaboration, is investigating biobased pathways to produce advantaged products including advanced biofuels, fuel intermediates, and bioproducts. ABF is integrating advanced computational tools for biological engineering, process design and data analysis, and economic/sustainability modeling tools into a comprehensive and dynamic platform for biomanufacturing of microbes and using them to produce key metabolic intermediates or beachhead molecules which may be derivatized into several distinct bioproducts of industrial interest. In this presentation, we first discuss a methodology to select a single exemplar product molecule to represent each beachhead pathway based on similarities with other end-molecule options (yields/titers/rates, fermentation operation mode, oxygen requirements, general separations challenges, etc.), in order to maintain a reasonable number of cases for rigorous process simulation. We then use techno-economic analysis (TEA) and life-cycle analysis (LCA) to highlight sensitivities and trends around economic performance and environmental footprint, reflecting examples for two selected ABF technology pathways to bio-derived chemicals: 1) adipic acid production via muconic acid fermentation from mixed sugars with Pseudomonas putida and 2) cineole via geranyl diphosphate with Rhodosporidium toruloides. We present multi-variable scan plots highlighting key drivers on costs and greenhouse gas emissions in order to identify the parameter space in which each pathway could ultimately achieve economic and environmental sustainability meeting or exceeding commodity product benchmarks, thus prioritizing future R&D focus areas.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

BISON Capability, Validation and Demonstration for Reactivity-Initiated Accidents

Reactivity Initiated Accidents (RIA) are design basis accidents that could have adverse impacts on the core coolability. In the unlikely event that sufficient reactivity is inserted into the reactor core by the ejected/dropped control rod, prompt energy deposition into the fuel can occur, which when sufficiently high can lead to fuel rod failure or, at large energy deposition levels, expulsion of UO2 fragments or molten UO2 material from the fuel rod. This results in a release of fission product and fuel into the coolant and has the potential to compromise core coolability and threatening the pressure boundary of the primary coolant system. The design basis RIA is one of the industry challenging problems that the CASL aimed to address; the CASL RIA Challenge Problem Charter [1] states, ?The Pressurized Water Reactor (PWR) Rod Ejection Accident (REA) and Boiling Water Reactor (BWR) Control Rod Drop Accident (CRDA) are postulated accidents with consequences that are important to nuclear safety (fuel rod integrity and core coolability). Currently each reload core design must be analyzed to meet regulatory acceptance criteria. The goal of CASLs ModSim capability for RIA is to model the event at a higher fidelity, with validation to existing tests, to better model the transient neutronics and the progression of the fuel and cladding thermal-mechanical behavior. These improved analytical capabilities can be used to better inform reload core design, limits on fuel assembly discharge burnup, restrictions on placement of fuel in the reactor, control rod insertion limits, operating margin, and performance sensitivities.? In support of that charter, BISON, the fuel performance code has been used to demonstrate the simulation of thermal-mechanical behavior of LWR fuels during a reactivity-initiated accident. The combination of mechanical, thermal, and thermal-hydraulic phenomena all present during an RIA makes a multi-physics code such as BISON a valuable tool for modeling these scenarios. This paper will highlight many of the activities associated with BISON relevant to reactivity-initiated accident capability development and validation efforts. These efforts have been performed under both CASL and the Nuclear Energy Advanced Modeling and Simulation (NEAMS) programs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Chlorine Worth Study Nuclear Data

The Chlorine Worth Study (CWS) was a series of experiments that took place at the National Criticality Experiments Research Center (NCERC), operated by Los Alamos National Laboratory (LANL). The focus of the experiments was to develop new integral benchmarks for the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook with high sensitivity to chlorine in the thermal neutron energy region and which match sensitivities of aqueous chloride operations at LANL. This work discusses the experiment and how the experimental and simulated results using different nuclear data libraries compare to each other.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Results of South-East Flux Trap Dosimetry Measurements for the Advanced Test Reactor Critical Facility in support of Advanced Sensors and Instrumentation Development

Reactor dosimetry measurements are commonly used to validate simulation and modeling in nuclear reactor tests. Numerous standard dosimeter materials exist which are commonly utilized with sensitivities to different energy ranges of neutrons. At the Advanced Test Reactor (ATR), cobalt alloy and pure nickel wires are installed every cycle to monitor thermal- and fast-neutron fluence rates. However, there is growing interest in exploring less commonly used materials which are either more sensitive to different parts of the neutron energy spectrum or which can incorporate multiple activation paths in a single material. Epithermal and fast-neutron energies beyond the typical 1-MeV threshold are of particular interest. Two ATR-C Flux Runs took place during 2024; each flux run included supplemental dosimetry packages in the South-East Flux Trap (SEFT) Filler. The focus of the dosimetry package for flux run 23-4 was to test two novel dosimetry methods that can provide simultaneous thermal and threshold (fast) sensitivity in a single dosimeter wire. A selection of 3% Au in Cu alloyed wires were available that provided sensitivity to fast and thermal neutrons through 5 different reactions. Likewise, Fe offers multiple interaction pathways with sensitivity to both thermal and fast neutrons. The main question to be answered by these irradiations was if sufficient radioactivation would take place in the ATR-C SEFT during a nominal 20-minute irradiation at typical power levels (near 600Wth) to observe the threshold reactions that have smaller activation cross-sections than the thermal reactions without being saturated by interfering interactions and Compton continuum during the High-Purity Germanium (HPGe) measurements. The results from comparing the measurement results to anticipated activity levels provide confidence in our ability to activate both traditional and novel dosimetry materials in ATR-C, however not all the measured values matched with the predicted activities. This leaves further room for investigation both on the experimental and computational approaches for future irradiation experiments.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2026

The Role of Nuclear Data Sensitivities in Prompt α-Eigenvalue Predictions of Delayed Critical Benchmarks

Alpha (α) eigenvalues, which describe the logarithmic time derivative of the neutron population in a multiplying system, are integral to time-dependent behavior and diagnostic applications. However, uncertainties in the evaluated nuclear data can significantly impact the accuracy of transport simulations for such quantities. This work explores the use of machine learning models to predict two key outputs, α-eigenvalues and keff bias, using input features derived from α-eigenvalue sensitivities to nuclear data. The criticality safety benchmark models used in this study come from the International Handbook of Evaluated Criticality Safety Benchmark Experiments. Three models, random forest, XGBoost, and NGBoost, are trained on both energy-resolved and energy-summed α sensitivities. For the α-eigenvalue bias prediction, NGBoost achieved the highest R 2 (0.9476) using energy-resolved features, while XGBoost performed best using summed sensitivities. In contrast, when predicting the keff bias, all the models showed moderate predictive capability (best R 2 ≈ 0.72), as the mapping from the static α-sensitivities to the static keff bias was less direct. SHAP (SHapley Additive exPlanations) analysis was used to interpret the model predictions. Across both prediction tasks, the features associated with neutron capture [H-1 (n, γ)], uranium scattering reactions (such as 235 U elastic/inelastic), and actinide capture/fission reactions (such as 239 Pu and 234 U) were consistently identified as the most impactful. This highlights the key role of specific nuclear reactions and energy ranges in shaping both time-dependent and steady-state criticality behavior. These results demonstrated that α-sensitivities, despite being computed for time-dependent metrics, can provide valuable insights for predicting both α-eigenvalues and the keff bias. Moreover, machine learning models offer a promising pathway for uncovering important nuclear data dependencies and guiding future data evaluation efforts.

Nuclear data↗
DOE OSTI2023

A semi-analytic study of self-interacting dark-matter haloes with baryons

ABSTRACT We combine the isothermal Jeans model and the model of adiabatic halo contraction into a semi-analytic procedure for computing the density profile of self-interacting dark-matter (SIDM) haloes with the gravitational influence from the inhabitant galaxies. The model agrees well with cosmological SIDM simulations over the entire core-forming stage up to the onset of gravothermal core-collapse. Using this model, we show that the halo response to baryons is more diverse in SIDM than in CDM and depends sensitively on galaxy size, a desirable feature in the context of the structural diversity of bright dwarfs. The fast speed of the method facilitates analyses that would be challenging for numerical simulations – notably, we quantify the SIDM halo response as functions of the baryonic properties, on a fine mesh grid spanned by the baryon-to-total-mass ratio, Mb/Mvir, and galaxy compactness, r1/2/Rvir; we show with high statistical precision that for typical Milky-Way-like systems, the SIDM profiles are similar to their CDM counterparts; and we delineate the regime of core-collapse in the Mb/Mvir − r1/2/Rvir space, for a given cross section and concentration. Finally, we compare the isothermal Jeans model with the more sophisticated gravothermal fluid model, and show that the former yields faster core formation and agrees better with cosmological simulations. We attribute the difference to whether the target CDM halo is used as a boundary condition or as the initial condition for the gravothermal evolution, and thus comment on possible improvements of the fluid model. We have made our model publicly available at https://github.com/JiangFangzhou/SIDM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Using Radiative Convective Equilibrium to Explore Clouds and Climate in the Community Atmosphere Model

Abstract Characteristics of, and fundamental differences between, the radiative‐convective equilibrium (RCE) climate states following the Radiative‐Convective Equilibrium Model Intercomparison Project (RCEMIP) protocols in the Community Atmosphere Model version 5 (CAM5) and version 6 (CAM6) are presented. This paper explores the characteristics of clouds, moisture, precipitation and circulation in the RCE state, as well as the tropical response to surface warming, in CAM5 and CAM6 with different parameterizations. Overall, CAM5 simulates higher precipitation rates that result in larger global average precipitation, despite lower outgoing longwave radiation compared to CAM6. Differences in the structure of clouds, particularly the amount and vertical location of cloud liquid, exist between the CAM versions and can, in part, be related to distinct representations of shallow convection and boundary layer processes. Both CAM5 and CAM6 simulate similar peaks in cloud fraction, relative humidity, and cloud ice, linked to the usage of a similar deep convection parameterization. These anvil clouds rise and decrease in extent in response to surface warming. More generally, extreme precipitation, aggregation of convection, and climate sensitivity increase with warming in both CAM5 and CAM6. This analysis provides a benchmark for future studies that explore clouds, convection, and climate in CAM with the RCEMIP protocols now available in the Community Earth System Model. These results are discussed within the context of realistic climate simulations using CAM5 and CAM6, highlighting the usefulness of a hierarchical modeling approach to understanding model and parameterization sensitivities to inform model development efforts.

Reed, Kevin A.↗
DOE OSTI2021

Exploring Western North Pacific Tropical Cyclone Activity in the High‐Resolution Community Atmosphere Model

Abstract High‐resolution climate models (∼28 km grid spacing) can permit realistic simulations of tropical cyclones (TCs), thus enabling their investigation in relation to the climate system. On the global scale, previous works have demonstrated that the Community Atmosphere Model (CAM) version 5 presents a reasonable TC climatology under prescribed present‐day (1980–2005) forcing. However, for the Western North Pacific (WNP) region, known biases in simulated TC genesis frequency and location under‐represent the basin's dominant share in observations. This study addresses these model biases in WNP by evaluating WNP TCs in a decadal simulation, and exploring potential improvements through nudging experiments. Among the major environmental controls of TC genesis, the lack of mid‐level moisture is identified as the leading cause of the deficit in simulated WNP TC genesis over the Pacific Warm Pool. Subsequent seasonal experiments explore the effect of constraining the large‐scale environment on TC development by nudging WNP temperature field toward reanalysis at various strengths. Temperature nudging elicits a significant response in TC genesis and intensity development, as well as in moisture and convection over the Warm Pool. These responses are sensitive to the choice of nudging timescale. Overall, the nudging experiments demonstrate that improvements in the large‐scale environment can lead to improvements in simulated TCs, suggesting future model developments in relation to model physics. In this way, the potential improvements in model fidelity will contribute to the understanding of how the mean state of current or future climates may give rise to extremes such as TCs.

58 GEOSCIENCES↗
DOE OSTI2021

An Integrative Model for Soil Biogeochemistry and Methane Processes: I. Model Structure and Sensitivity Analysis

Abstract Environmental changes are anticipated to generate substantial impacts on carbon cycling in peatlands, affecting terrestrial‐climate feedbacks. Understanding how peatland methane (CH 4 ) fluxes respond to these changing environments is critical for predicting the magnitude of feedbacks from peatlands to global climate change. To improve predictions of CH 4 fluxes in response to changes such as elevated atmospheric CO 2 concentrations and warming, it is essential for Earth system models to include increased realism to simulate CH 4 processes in a more mechanistic way. To address this need, we incorporated a new microbial‐functional group‐based CH 4 module into the Energy Exascale Earth System land model (ELM) and tested it with multiple observational data sets at an ombrotrophic peatland bog in northern Minnesota. The model is able to simulate observed land surface CH 4 fluxes and fundamental mechanisms contributing to these throughout the soil profile. The model reproduced the observed vertical distributions of dissolved organic carbon and acetate concentrations. The seasonality of acetoclastic and hydrogenotrophic methanogenesis—two key processes for CH 4 production—and CH 4 concentration along the soil profile were accurately simulated. Meanwhile, the model estimated that plant‐mediated transport, diffusion, and ebullition contributed to ∼23.5%, 15.0%, and 61.5% of CH 4 transport, respectively. A parameter sensitivity analysis showed that CH 4 substrate and CH 4 production were the most critical mechanisms regulating temporal patterns of surface CH 4 fluxes both under ambient conditions and warming treatments. This knowledge will be used to improve Earth system model predictions of these high‐carbon ecosystems from plot to regional scales.

58 GEOSCIENCES↗
DOE OSTI2022

Characterization of nonlinear ultrasonic waves behavior while interacting with poor interlayer bonds in large-scale additive manufactured materials

Over the past decades, researchers have developed several nonlinear ultrasonic techniques for quality control of materials commonly used in different applications. Owing to the superior sensitivity of nonlinear ultrasound waves to small defects such as micro-cracks, their applicability in different nondestructive testing (NDT) problems has been investigated in numerous studies. These studies utilize frequency domain analysis to detect the generation of higher harmonics because of the formation of defects in the inspected medium. Frequency domain analysis based on the Fourier transform is a significant approach used in linear systems; however, it may not perform adequately on nonlinear systems. Hence, studies on nonlinear dynamics and physics consider analyzing systems' behavior in the phase-space domain. In contrast to the frequency domain analysis, which can result in information loss, analysis in the phase-space domain retains all the information regarding a system's states. Here, we investigate the nonlinearities induced by poor interlayer bonding in polymer-based additive manufactured parts in the phase-space domain. It is convenient to characterize the nonlinearity in the phase-space domain because it provides a geometrical representation of a system's states. Two types of low quality interlayer bond are considered. The first type is simulated artificially while the second type is manufactured by reducing the bond quality during the printing process. The analysis verified that the received ultrasonic signals exhibit classical nonlinear behavior in the phase-space domain while interacting with simulated poor interlayer bonds. In addition, the results showed that the behavior of ultrasonic waves is amplitude-dependent and evolves into models that have not been previously reported. Furthermore, Largest Lyapunov Exponent (LLE) is used to quantify the behavior of nonlinear ultrasonic waves while interacting with poor interlayer bonds. Using LLE, it was observed that the divergence rate of the phase-space trajectories depends on the amplitude of the excitation. This observation quantitatively proves that nonlinear behavior of ultrasound while interacting with poor interlayer bonds can be amplitude-dependent. The results of both simulated and inherent poor interlayer bond cases showed that LLE can be used as a reliable quantitative damage-sensitive feature to detect and potentially characterize weak bonds, which are difficult to detect using conventional approaches. Additionally, the reported results in the phase-space domain provide a basis for proposing a new mathematical model for ultrasonic waves interacting with poor interlayer bonds.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Improved Modeling of Waste Isolation Pilot Plant Disposal Room Porosity

Disposal rooms at the Waste Isolation Pilot Plant (WIPP) contain waste and gas, and their porosity evolves over time. This report presents several improvements to the disposal room porosity model and presents new porosity predictions for use in future WIPP Performance Assessment activities. The improvements pertain to three sub-models: the geomechanical model, the waste compaction model, and the gas pressurization model. The impacts of each major improvement were quantified and the new porosity predictions were shown to be both mesh and domain size converged. Also, sensitivity studies on the disposal room horizon, clay seam friction coefficients, and homogenized waste representation were performed to support assumptions in the disposal room porosity model. To compare the legacy and new porosity predictions, the simulation results were plotted as a response surface, where gas pressure and time are inputs and porosity is the output. The new porosity response surface is insensitive to pressures beneath lithostatic pressure and highly sensitive to pressures above lithostatic pressure. The legacy porosity response surface, on the other hand, has moderate porosity gradients over all pressures. The new porosity response surface has a stronger scientific foundation than the legacy surface and may now be used for Compliance Decision Analyses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Uncertainty Quantification in High-Low Dynamic System Coupling using RAVEN and TRANSFORM

This work demonstrates new functionality and applications stemming from the development of high-fidelity to low-fidelity (high-low) coupling for system simulations and to further explore the capabilities of the Risk Analysis Virtual Environment (RAVEN) in the performance of uncertainty quantification in this kind of high-low coupled system models. The work builds from previous work on high-low coupling that utilized COBRA-TF (CTF), the high-fidelity subchannel analysis code, with a low fidelity model built in ORNL’s TRANSFORM, the system analysis code, utilizing the Functional Mock-Up Interface (FMI). Steady-state and transient analysis examples using the high/low coupled models generated from CTF and TRANSFORM/FMI are investigated. The workflows for both steady-state and transient coupled simulations are described. A steady-state parameter sweep and uncertainty analysis of the primary flow rates and reactor power are demonstrated. Likewise, a transient pump trip and power ramp sensitivity studies are also demonstrated. This work elucidates some of the potential benefits and future needs of using RAVEN for high/low system coupling analysis of energy systems. It also shows some of the difficulties that can be encountered in coupling system simulations.

Williams, Wesley↗
DOE OSTI2022

Sensitivity analyses of a homogeneous model and a RZ model of the MYRRHA reactor in its critical configuration

Two homogenized models of the MYRRHA (Multi-purpose hybrid Research Reactor for High-tech Applications) reactor in its critical configuration are presented in this study. Both models will be verified by comparing the main reactor parameters (K{sub eff}, β{sub eff} and Λ) with those of the heterogeneous model and a sensitivity analysis will be performed. A sensitivity analysis of some important parameters related to reactivity, namely the vacuum coefficient, the Doppler and the power peaking factors, is also presented. The study of the vacuum coefficient has been carried out by analyzing two possible scenarios: 50% vacuum in the coolant (lead bismuth eutectic) and 100%, obtaining a list of the 10 reactions and nuclides that most affect the k{sub eff} value. With respect to the Doppler coefficient, it will be seen that a 300 degree increase in fuel temperature (starting at 800 K) results in a reduction of the total reactivity of the system by 134 pcm.. Finally, the map of the power peaking factors will be shown, as well as a sensitivity analysis. For making the sensitivity calculations Serpent-2 reactor physics Monte Carlo code has been used and JEFF-3.3 library has been chosen in order to continue with the OECD/NEA WPEC SG46 work evaluating the nuclear data of this library. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Heat pulse testing at monitoring wells to estimate subsurface fluid velocities in geological CO 2 storage

Monitoring the injected CO 2 during geological CO 2 storage (GCS) is essential to assure containment and identify CO 2 leakage. Here in this work, a new approach is introduced to estimate the evolution of the downhole fluid velocity at a monitoring well and identify CO 2 arrival time using in-well heat pulse/tracer test. The proposed technique involves using a downhole heater to generate a series of heat pulses and measuring their corresponding temperature response. The surface temperature of the downhole heater is controlled by the supplied electrical power and the heat loss by convection to the surroundings. Convective heat transfer is well described using Newton's law of cooling in which the temperature difference between the heater and the surrounding fluids drives the heat transfer, for which the convection heat transfer coefficient (h) controls the magnitude of heat loss. Among various factors that control h, it depends on the type of the flowing fluid and its velocity. Through analyzing the measured temperature at different heat pulses, the changes in h - due to mobilization of the in-situ brine or CO 2 arrival - can be estimated. Consequently, the velocity of the flowing fluid across the heater can be obtained. Since heat transfer by convection is sensitive to the type of the surrounding fluid, intrusion of CO 2 can be detected from the relatively higher surface temperature obtained at CO 2 arrival. Churchill and Bernstein (1977)'s correlation is adopted to estimate the change of fluid velocity in terms of the change in h. To demonstrate the validity of the proposed technique, the results are applied and validated against those of COMSOL Multiphysics simulation tool for single-phase brine (before CO 2 arrival) and single-phase CO 2 (after CO 2 arrival). The observed temperature heating is sensitive to the flowing fluid velocity and fluid type. The temperature signal observed at CO 2 arrival is large and easily detectable using temperature monitoring tool which provides reliable indication for tracking CO 2 arrival at monitoring wells compared with passive temperature monitoring. The results obtained using the proposed technique agree very well with the numerical results obtained from the simulation tool with a maximum estimation error of 7 percent.

02 PETROLEUM↗
DOE OSTI2021

Vadose and Saturated Zone Flow and Transport Model Package Report for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington

This model package report (MPR) documents the development of the integrated vadose and saturated zone flow and transport model developed for the performance assessment of Trenches 31 and 34 in the low-level burial grounds (LLBGs). This modeling capability is intended for use in addressing the analysis requirements outlined in DOE O 435.1, Chg 1, Radioactive Waste Management1. The overall objective of the modeling effort is to provide a basis for making informed disposal decisions pertinent to Trenches 31 and 34. The purpose of the MPR is to document the development of the three-dimensional numerical vadose zone and saturated flow and transport model test case and evaluate its adequacy to support the LLBG performance assessment. The purpose is not to present results for DOE 435.1 decision making. The use of the model to perform base case and sensitivity analysis for the LLBG performance assessment, including inputs and results, is documented separately in subsequent environmental calculation files. This report discusses the development and translation of the conceptual model for flow and contaminant transport into the LLBG performance assessment three-dimensional numerical flow and transport model evaluated using the Subsurface Transport Over Multiple Phases (STOMP©) simulator. The development of representative geologic framework is described along with the implementation of waste release models used to represent contaminant release from waste disposed in the trenches. The report also provides the technical basis for specific model parameters and boundary conditions, along with description of modeling assumptions. This MPR includes certain calculations that are necessary to demonstrate the soundness of the model. Results provided by the model include vadose zone and saturated zone flow fields, and estimates of the possible future concentration in groundwater of technetium-99 and iodine-129 as example test cases. As an evaluation of the test cases, the model estimates of current vadose conditions are compared to available and analogous field and laboratory data, and compared to the results from the original performance assessment model (WHC-EP-0645, Performance Assessment for the Disposal of Low-Level Waste in the 200 West Area Burial Grounds 2 ). The features, events, and processes applicable to vadose zone and saturated zone flow and transport model are identified, and representative initial estimates for various parameters are documented. Note that the parameter estimates presented in this MPR are for illustration purposes, and may or may not reflect values selected for eventual performance assessment. Several key topical discussions (i.e., basis for recharge estimates, basis for vadose zone modeling, basis for saturated zone model development, and calibration) are included, which serve as the groundwork for confidence building for the groundwater pathway modeling and results. Numerical simulation results based on an example test case are included to illustrate the use of the combined saturated-unsaturated model in performance assessment calculations. Sensitivity and uncertainty results are not included in this MPR. Those results will be included in future environmental calculation files. The inclusion of the trapezoidal trench geometry and construction details in the finite difference grid introduces some gross simplifications regarding the trench and liner systems. The model test case presented in this MPR includes the assumption that the trench liner system does not affect flow through or around the trenches after the liner system is assumed to fail. This MPR also includes results of other test cases that involve alternate assumptions about how to incorporate the hydraulic effects of the trenches into the vadose zone of the model. These alternate cases do not attempt to account for the presence of the liner system after its assumed failure either. None of these cases is considered to be the base case at this time. Analysis and alternate cases that attempt to account for the presence of the liner system in greater detail, and its effect on flow through and around it, are to be documented in subsequent environmental calculation files. Depending on the results of those analyses and alternate cases, the eventual base case may involve more detailed inclusion of the effects of the liner system hydraulics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Conceptual Design Report for the MATHUSLA Long-Lived Particle Detector near CMS

We present the Conceptual Design Report (CDR) for the MATHUSLA (MAssive Timing Hodoscope for Ultra-Stable neutraL pArticles) long-lived particle detector at the HL-LHC, covering the design, fabrication and installation at CERN Point 5. MATHUSLA is a 40 m-scale detector with an air-filled decay volume that is instrumented with scintillator tracking detectors, to be located near CMS. Its large size, close proximity to the CMS interaction point and about 100 m of rock shielding from HL-LHC backgrounds allows it to detect LLP production rates and lifetimes that are one to two orders of magnitude beyond the ultimate sensitivity of the HL-LHC main detectors for many highly motivated LLP signals. Data taking is projected to commence with the start of HL-LHC operations. We present a new 40m design for the detector: its individual scintillator bars and wavelength-shifting fibers, their organization into tracking layers, tracking modules, tower modules and the veto detector; define a high-level design for the supporting electronics, DAQ and trigger system, including supplying a hardware trigger signal to CMS to record the LLP production event; outline computing systems, civil engineering and safety considerations; and present preliminary cost estimates and timelines for the project. We also conduct detailed simulation studies of the important cosmic ray and HL-LHC muon backgrounds, implementing full track/vertex reconstruction and background rejection, to ultimately demonstrate high signal efficiency and $\ll 1$ background event in realistic LLP searches for the main physics targets at MATHUSLA. This sensitivity is robust with respect to detector design or background simulation details. Appendices provide various supplemental information.

Aitken, Branden [Victoria U.]↗