Results for “Validated Models” — Page 58

DOE OSTI2020

Global biogeography of fungal and bacterial biomass carbon in topsoil

Bacteria and fungi, representing two major soil microorganism groups, play an important role in global nutrient biogeochemistry. Biogeographic patterns of bacterial and fungal biomass are of fundamental importance for mechanistically understanding nutrient cycling. We synthesized 1323 data points of phospholipid fatty acid-derived fungal biomass C (FBC), bacterial biomass C (BBC), and fungi:bacteria (F:B) ratio in topsoil, spanning 11 major biomes. The FBC, BBC, and F:B ratio display clear biogeographic patterns along latitude and environmental gradients including mean annual temperature, mean annual precipitation, net primary productivity, root C density, soil temperature, soil moisture, and edaphic factors. At the biome level, tundra has the highest FBC and BBC densities at 3684 (95% confidence interval: 1678–8084) mg kg -1 and 428 (237–774) mg kg -1 , respectively; desert has the lowest FBC and BBC densities at 16.92 (14.4–19.89) mg kg -1 and 6.83 (6.1–7.65) mg kg -1 , respectively. The F:B ratio varies dramatically, ranging from 1.8 (1.6–2.1) in savanna to 8.6 (6.7–11.0) in tundra. An empirical model was developed for the F:B ratio and it is combined with a global dataset of soil microbial biomass C to produce global maps for FBC and BBC in 0–30 cm topsoil. Across the globe, the highest FBC is found in boreal forest and tundra while the highest BBC is in boreal forest and tropical/subtropical forest, the lowest FBC and BBC are in shrub and desert. Global stocks of living microbial biomass C were estimated to be 12.6 (6.6–16.4) Pg C for FBC and 4.3 (0.5–10.3) Pg C for BBC in topsoil. These findings advance our understanding of the global distribution of fungal and bacterial biomass, which facilitates the incorporation of fungi and bacteria into Earth system models. The global maps of bacterial and fungal biomass serve as a benchmark for validating microbial models in simulating the global C cycle under a changing climate.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Modeling the impact of the fuel injection strategy on the combustion and performance characteristics of a heavy-duty GCI engine

Gasoline compression ignition (GCI) is a promising strategy to achieve high thermal efficiency and low emissions with limited modifications to the conventional diesel engine hardware. It is a partially premixed concept, which derives its superiority from higher volatility and longer ignition delay of gasoline-like fuels combined with higher compression ratio typical of diesel engines. Here, the present study investigates the combustion process in a GCI engine operating with different injection strategies using computational fluid dynamics (CFD). Simulations are carried out on a single cylinder of a multi cylinder heavy-duty compression ignition engine, which operates at a compression ratio of 17:1 and an engine speed of 1038 rev/min. Two different injection strategies viz., late injection (LI), and early pilot injection (EP) are investigated to understand their impact on combustion and performance of the engine. Renormalized group (RNG) k-ε model is used to describe in-cylinder turbulence and KH RT model is used to simulate the fuel spray breakup. The developed CFD methodology is validated against relevant experimental data under a wide range of operating conditions for each injection strategy. The developed CFD methodology was found to capture the engine combustion behavior quite well. Based on the validated CFD model, the differences in the progress of combustion event for the two injection strategies is highlighted. It was found that a larger pilot fuel mass fraction results in a steeper rise in the initial heat release rate which in turn influences the transition to mixing controlled combustion. In line with the experimental data, the study showed that the late pilot injection strategy with three injection pulses, results in higher performance compared to the other conditions.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Conflict Detection in Open RAN with Recurrent Neural Networks Using Geometric Manifolds

Allowing third-party applications on Radio Access Network (RAN) Intelligent Controllers (RICs) within the OpenRAN (O-RAN) framework introduces conflicting interactions that are often difficult to detect in advance. These conflicts, occurring between third-party applications in the Near RealTime RIC (Near-RT RIC), known as xApps, can lead to performance degradation and instability in O-RAN if not identified early. Existing conflict detection and mitigation solutions in the literature assume that the conflicts are known beforehand, which is not always accurate due to the complex and often hidden relationships between control parameters and Key Performance Indicators (KPIs). In this paper, we propose a novel Recurrent Neural Network (RNN) to detect both known and unknown conflicts in O-RAN xApps as specified in the O-RAN standards. We model the xApps, control parameters, and KPIs with nodes and edges to create graph structures and use the hidden nonEuclidean geometric properties of the Riemannian manifold to train the RNN model. The performance of this proposed model is validated using evaluation metrics and compared with benchmarks. Results demonstrate that the proposed RNN model, leveraging Riemannian geometric properties, can achieve 100% of the F1-score provided by an optimal solution in just 20 iterations.

5G↗

DOE OSTI2025

Conflict Detection in Open RAN with Recurrent Neural Networks Using Geometric Manifolds

Allowing third-party applications on Radio Access Network (RAN) Intelligent Controllers (RICs) within the OpenRAN (O-RAN) framework introduces conflicting interactions that are often difficult to detect in advance. These conflicts, occurring between third-party applications in the Near RealTime RIC (Near-RT RIC), known as xApps, can lead to performance degradation and instability in O-RAN if not identified early. Existing conflict detection and mitigation solutions in the literature assume that the conflicts are known beforehand, which is not always accurate due to the complex and often hidden relationships between control parameters and Key Performance Indicators (KPIs). In this paper, we propose a novel Recurrent Neural Network (RNN) to detect both known and unknown conflicts in O-RAN xApps as specified in the O-RAN standards. We model the xApps, control parameters, and KPIs with nodes and edges to create graph structures and use the hidden nonEuclidean geometric properties of the Riemannian manifold to train the RNN model. The performance of this proposed model is validated using evaluation metrics and compared with benchmarks. Results demonstrate that the proposed RNN model, leveraging Riemannian geometric properties, can achieve 100% of the F1-score provided by an optimal solution in just 20 iterations.

5G↗

DOE OSTI2020

SYSTEMATIC MULTI-INDEX VALIDATIONS OF SIDI ENGINE FLOW FIELD LES COMPUTATIONS USING CRANK ANGLE-RESOLVED PIV MEASUREMENTS

In-cylinder flow fields make significant impacts on the fuel atomization, fuel mixture formation, and combustion process in spark ignition direct injection (SIDI) engines. In recent years, model-based simulation approaches are preferred in regard to investigating the transient in-cylinder flow field characteristics. Most commonly, the simulation models are validated using single representative flow field at a typical crank angle measured by particle image velocimetry (PIV). However, it provides only limited knowledge about the flow field which is highly three-dimensional and of transient nature. In this study, crank angle-resolved Ply measurements are conducted on three distinct planes inside the cylinder to capture the transient process of flow field characteristics which vary with the crank angle. These three planes consist of one tumble plane through the spark plug tip, one tumble plane along two intake ports, and one swirl plane at 30 mm below the cylinder head. Large eddy simulation (LES) is employed for the numerical computations using the CONVERGE codes. On the basis of large datasets for both temporal and spatial domains, a multi-index systematic validation approach is conducted. Crank angle-resolved calculations of global indices and local indices are implemented using the flow fields velocity data obtained from both PIV and LES on select planes. Global indices reveal the trends in similarities of different crank angle degrees and locations, while local indices give the detail comparison results. In summary, with the systematic multi-index validation approach, the crucial crank angle degrees and locations for model verification will be detected. Furthermore, the corresponding critical flow features are analyzed. Practical guideline of flow field validation is proposed.

Crank angle-resolved PIV↗

DOE OSTI2025

Summary of UMI Chemical Work

Our initial work for the Urban Modeling Integration (UMI) project was to validate the QUIC model dense gas dispersion of chlorine advection up a ramp. The validation task was to use QUIC to model an experiment done in a wind tunnel. Naturally the QUIC model uses a slumping cylindrical model where inside the cylinder the gas is treated as dense but outside ideal. In trying to validate using the wind tunnel experiment, due to the scale of the wind tunnel experiment and ramp, The model formed a very clear, discrete, concentration barrier at the edge of the cylinder that likely doesn’t line up well with the real dense gas advection/diffusion. This effect is much less noticeable at a larger scale such as an instantaneous release of a 25 ton tank.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Dynamics of Gas Bubbles From a Submarine Hydrocarbon Seep Within the Hydrate Stability Zone

We validate a new model for mass transfer and bubble transport for natural seeps on the continental margins using an integrated observation of a seep at 883 m in the Northern Gulf of Mexico. In the model, mass transfer is assumed to transition from clean to dirty bubble mass transfer rates following a characteristic hydrate formation time that depends on the initial bubble surface area and the hydrate subcooling. We show that buoyancy-induced upwelling is negligible for the bubble stream. We initialize the model using precise data observed at the seafloor and validate the model results to optical and acoustic measurements within the bubble stream. The model accurately predicts the acoustic attenuation of the bubble flare by lateral spreading and bubble shrinkage throughout the water column. In conclusion, the model shows that up to 99.4% by mass of the released gases from this seep dissolves into the ocean within the hydrate stability zone.

03 NATURAL GAS↗

DOE OSTI2023

A parametric finite element study for determining burst strength of thin and thick-walled pressure vessels

To accurately predict the burst strength of both thin and thick-walled pressure vessels (PVs), a parametric study of PV burst strength was performed for a wide range of vessel geometries and materials using elastic-plastic finite element analysis (FEA). A valid FEA model was established through a detailed study of 2D versus 3D FEA models, the critical stress failure criterion versus the limit load criteria, and the thick-wall effect on the FEA simulations. Here, the results show that the stresses and strains at the mean diameter, rather than outside diameter, determines a more accurate burst strength for both thin and thick-walled PVs. On this basis, a parametrized FEA script using the ABAQUS Python application programming interface (API) was used to create a large database of PV burst strengths for a variety of vessel geometries and materials, demonstrating that Python scripting is a powerful technique for performing parametric studies or generating large databases. From the FEA results, using the regression method, a new burst pressure model was developed as a function of the vessel geometry (D/t ratio) and material properties (UTS and n). As validated by a large number of full-scale burst test data, the proposed burst model can very accurately predict the burst strength for both thin and thick-walled PVs.

42 ENGINEERING↗

DOE OSTI2022

Kinetic Model for the Reduction of Cu II Sites by NO + NH 3 and Reoxidation of NH 3 -Solvated Cu I Sites by O 2 and NO in Cu-SSZ-13

In this work, a kinetic model is developed for the reduction of Cu II sites by NO + NH 3 and the reoxidation of NH 3 -solvated Cu I sites by O 2 and NO in Cu-SSZ-13. Fourier transform infrared (FTIR) spectroscopy and spatially resolved capillary inlet mass spectrometry (SpaciMS) measurements during transient reactor experiments are utilized to identify the rate parameters associated with NO + NH 3 RHC (reduction half-cycle), proposed to occur via two distinct pathways involving adsorbed NH 3 and gas-phase NH 3 . The resulting NO + NH 3 RHC model is validated using spatiotemporal N 2 measurements covering a wide range of temperatures (200–450 °C) and space velocities (53,000–640,000 h –1 ). N 2 O formation is observed and modeled during NO + NH 3 RHC, with quantitative validation under standard selective catalytic reduction (SCR) conditions. Experimentally measured enthalpic and entropic changes associated with O 2 adsorption on NH 3 -solvated Cu I (ZCu(NH 3 ) 2 ) complexes [ Kamasamudram, K. Catal. Today 2010 , 151 (3–4), 212-222], along with activation energies estimated computationally for the intercage diffusion of ZCu(NH 3 ) 2 complexes [ Paolucci, C. Science 2017 , 357 (6 354), 898-903], are incorporated into a mean field kinetic model for the low-temperature oxidation half-cycle (OHC). Significant NH 3 release is observed during the isothermal oxidation of Cu I sites, attributed to desorption of NH 3 ligands from NH 3 -solvated Cu II dimers (Z 2 Cu 2 (NH 3 ) 4 O 2 ). Reduction of these dimeric complexes leads to the consumption of one NO/Cu II , contradicting the expected reduction stoichiometry. Inclusion of a global Arrhenius rate for the NO titration of Z 2 Cu 2 (NH 3 ) 4 O 2 complexes provides accurate representations of standard SCR on reduced and oxidized catalysts, predicting transient NO and NH 3 consumption between 150 and 250 °C as a function of hydrothermal aging. Deactivation of low-temperature standard SCR by NH 3 is observed at high NH 3 pressures, modeled via the formation of superoxo amino (ZCu(NH 3 ) 3 OO*) complexes during NH 3 titration of Z 2 Cu 2 (NH 3 ) 4 O 2 complexes [ Negri, C. J. Am. Chem. Soc. 2020 , 142 (37), 15884-15896]. The redox kinetic model presented here provides a foundational description of active site redox during low-temperature standard SCR, combining the recent kinetic, spectroscopic, and computational findings on the mechanism of standard SCR over Cu-SSZ-13.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

A coupled thermo-hydro-mechanical model for simulating leakoff-dominated hydraulic fracturing with application to geologic carbon storage

A potential risk of injecting CO2 into storage reservoirs with marginal permeability (≲ 10 mD (1 mD = 10 –15 m 2 )) is that commercial injection rates could induce fracturing of the reservoir and/or the caprock. Such fracturing is essentially fluid-driven fracturing in the leakoff-dominated regime. Recent studies suggested that fracturing, if contained within the lower portion of the caprock complex, could substantially improve the injectivity without compromising the overall seal integrity. Modeling this phenomenon entails complex coupled interactions among the fluids, the fracture, the reservoir, and the caprock. Here, we develop a simple method to capture all these interplays in high fidelity by sequentially coupling a hydraulic fracturing module with a coupled thermal-hydrological-mechanical (THM) model for nonisothermal multiphase flow. The model was made numerically tractable by taking advantage of self-stabilizing features of leakoff-dominated fracturing. The model is validated against the PKN solution in the leakoff-dominated regime. Moreover, we employ the model to study thermo-poromechanical responses of a fluid-driven fracture in a field-scale carbon storage reservoir that is loosely based on the In Salah project's Krechba reservoir. The model reveals complex yet intriguing behaviors of the reservoir-caprock-fluid system with fracturing induced by cold CO 2 injection. We also study the effects of the in situ stress contrast between the reservoir and caprock and thermal contraction on the vertical containment of the fracture. The proposed model proves effective in simulating practical problems on length and time scales relevant to geological carbon storage.

58 GEOSCIENCES↗

DOE OSTI2024

Thermal and dynamic partition of dumbbell interstitials in complex concentrated alloys

Complex concentrated alloys (CCAs) are promising candidates for applications in extreme conditions, such as irradiation where interstitial mediated diffusion is important. In CCAs with N principal elements, $\frac{N(N+1)}{2}$ types of dumbbell interstitials exist. Currently, there is no way to predict the thermal partition (fractional concentration at equilibrium) and the dynamic partition (fractional time an interstitial spends during diffusion) of each type of dumbbell interstitial. To mitigate this issue, this work proposes a theoretical model for computing the equilibrium concentrations and thermal partition of dumbbell interstitials in CCAs and validates the model using grand canonical Monte Carlo simulations. Lattice kinetic Monte Carlo simulations show that the thermal partition is equivalent to the dynamic partition, and both are governed by composition and formation energies of dumbbells. As a result, the model proposed provides a foundation for understanding radiation enhanced diffusion and induced segregation in CCAs under irradiation.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Induced Polarization of Clayey Rocks and Soils: Non‐Linear Complex Conductivity Models

Abstract The past decades have witnessed the increased applications of induced polarization (IP) method in the critical zone studies with ubiquitous clay minerals. Although IP outperforms traditional electrical and electromagnetic methods through its unique ability to measure quadrature conductivity, the nonlinearity that quadrature conductivity behaves with salinities and frequencies greatly tortures IP practitioners, as (a) salinity‐dependency makes the quadrature conductivity a varyingly unstable parameter to quantitatively estimate hydraulic properties and clay content; (b) frequency‐dependent Cole‐Cole and Debye/Warburg decomposition models, although mathematically sound, physically mingle the properties of pore water and clay minerals and are empirical in nature. From basic principles, we demonstrate that quadrature conductivity remains a hybrid property involving both clay and water, and develop relevant models to distinguish them. Our models are validated by theories, experiments, simulations, and comparisons, all of which proclaim considerable advantages over previous models and offer the prospect of quantitative applications.

58 GEOSCIENCES↗

DOE OSTI2021

Effect of workpiece curvature on the tool influence function during hemispherical sub-aperture tool glass polishing

The influence of workpiece curvature on the tool influence function spot during polishing of fused silica glass with cerium oxide slurry, while using a rotating hemispherical pad-foam tool for a wide variety of process conditions (tool displacement, inclination angle, and rotation rate), has been investigated. (Workpiece curvature ranged from 500 mm radius concave to 43 mm radius convex.) The TIF spot decreases in diameter and increases in the peak removal rate on more convex workpieces. In contrast, the TIF spot increases both in diameter and peak removal rate on more concave workpieces. For the range of workpiece curvatures investigated, both the spot size and the peak removal rate changed significantly, as much as 2 times. An elastic sphere–sphere contact mechanics model, which utilizes both a modified displacement (that leads to a change in the applied load) as well as a mismatch factor (that influences the pressure distribution shape), has been developed. Furthermore, the model was validated using both offline load-displacement measurements and finite-element analysis simulations. The model quantitatively describes the measured change in the relative contact diameter and relative pressure distribution, as well as semiquantitively describes the change in the relative volumetric removal rate on a large variety of TIF spots. The change in the volumetric removal rate for convex workpieces is a result of the balance between a decreasing spot size (reducing removal) and an increasing peak pressure (increasing removal), which usually results in relatively small changes in volumetric removal. In the case of concave workpieces, the volumetric removal rate change is also governed by a similar balance, but the spot size increase contribution dominates, resulting in a significant increase in volumetric removal rate. Understanding these trends can enable methods to add greater determinism during the fabrication of freeform optics by adjusting polishing parameters (such as dwell time) while the tool translates along a workpiece surface with different local curvatures.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Advanced Ceramic Membranes/Modules for Ultra Efficient Hydrogen (H2) Production/Carbon Dioxide (CO2) Capture for Coal-Based Polygeneration Plants: Fabrication, Testing, and CFD Modeling

Inorganic membrane-based systems are a promising technology for precombustion CO2 capture with simultaneous H2 production. State-of-the-art packages for high temperature and pressure service consist of multiple tube membrane bundles prepared in a "candle filter" configuration, in which the membrane tubes are open at one end and sealed at the other. This configuration is used for practical reasons, specifically the need to minimize problems due to thermal expansion mismatch between the ceramic tube bundle and the steel housing. However, the primary technical problem with the candle filter format for commercial-scale installations is the inability to purge the tube side (typically the permeate side), a feature that is crucial for high H2 recovery. In this study, the focus is to design and fabricate the first dual-end open full ceramic multiple tube membrane bundle that enables tube side (permeate) purge for gas separation applications. An additional key feature of this design is the simplified module layout, as the membrane bundles can be installed end-to-end with tube side (permeate) flow directly from one bundle to the next. This layout simplifies the membrane to housing seals and yields significant improvement in membrane packing density. Detailed focus areas in our studies include: (i) Materials development and preparation of the tube-to-tube sheet potting for the dual-ended bundle; (ii) the sealing and optimal module configuration design to minimize membrane stress upon module mounting; (iii) the demonstration, via the fabrication of CMS and Pd-alloy membranes supported on full-size, dual-ended ceramic support bundles, of the first example of a purgeable ceramic membrane and module; and (iv) development of a CFD model of the membrane module for calculation of feed flow distribution, and for use in scale-up, and capital cost estimating. The CFD model was validated using experimental data with the multi-tubular membrane system, employing He/N2 as a model gas mixture (surrogate for H2/CO2), and has been shown to be quite accurate. Employing the model, we are able to study the effects of operating pressure and temperature, feed and sweep gas flow rates, and the choice of membrane tube configuration on system performance.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2024

FY24 Advanced Computing HPC Annual Report for allocation "vtocei"

This project aims to elucidate the mechanisms underlying the formation of the cathodeelectrolyte interphase (CEI) and its impact on the performance of Li-ion batteries during electrochemical cycling. Our approach integrates advanced multiscale modeling with multimodal characterization to offer a comprehensive understanding of CEI dynamics. We employ a validated multiscale modeling framework to analyze microstructure-dependent transport properties, alongside joint theory-experiment protocols for detailed resolution of interfacial chemistry via spectroscopy.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Advanced Ceramic Membranes/Modules for Ultra Efficient Hydrogen (H2) Production/Carbon Dioxide (CO2) Capture for Coal-Based Polygeneration Plants: Fabrication, Testing and CFD Modeling

Inorganic membrane-based systems are a promising technology for precombustion CO2 capture with simultaneous H2 production. State-of-the-art packages for high temperature and pressure service consist of multiple tube membrane bundles prepared in a "candle filter" configuration, in which the membrane tubes are open at one end and sealed at the other. This configuration is used for practical reasons, specifically the need to minimize problems due to thermal expansion mismatch between the ceramic tube bundle and the steel housing. However, the primary technical problem with the candle filter format for commercial-scale installations is the inability to purge the tube side (typically the permeate side), a feature that is crucial for high H2 recovery. In this study, the focus is to design and fabricate the first dual-end open full ceramic multiple tube membrane bundle that enables tube side (permeate) purge for gas separation applications. An additional key feature of this design is the simplified module layout, as the membrane bundles can be installed end-to-end with tube side (permeate) flow directly from one bundle to the next. This layout simplifies the membrane to housing seals and yields significant improvement in membrane packing density. Detailed focus areas in our studies include: (i) Materials development and preparation of the tube-to-tube sheet potting for the dual-ended bundle; (ii) the sealing and optimal module configuration design to minimize membrane stress upon module mounting; (iii) the demonstration, via the fabrication of CMS and Pd-alloy membranes supported on full-size, dual-ended ceramic support bundles, of the first example of a purgeable ceramic membrane and module; and (iv) development of a CFD model of the membrane module for calculation of feed flow distribution, and for use in scale-up, and capital cost estimating. The CFD model was validated using experimental data with the multi-tubular membrane system, employing He/N2 as a model gas mixture (surrogate for H2/CO2), and has been shown to be quite accurate. Employing the model, we are able to study the effects of operating pressure and temperature, feed and sweep gas flow rates, and the choice of membrane tube configuration on system performance.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2020

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Post-DTL Beam Delivery

Ensuring that the beam delivered from the upgraded Front-End (FE) meets the Key Performance Parameters (KPPs) at each user facility is critical to the success of the LANSCE Accelerator Modernization Project (LAMP). For a high-intensity, multi-user facility like LANSCE, compliance with beam loss and radiation thresholds is as important as the charge delivered to each target. While early LAMPF/LANSCE operations relied on iterative tuning to minimize losses from beam halo and tail particles, the new FE may introduce different beam distributions and loss modes—making predictive modeling essential. To manage this, the F2E (Front-End to End) effort is developing detailed particle-tracking models that reflect realistic beamline conditions, including halo formation and expected diagnostic readings. These "snapshot" simulations aim to benchmark live machine performance at a given moment. This will help quantify how beam quality from the new FE will propagate downstream through the facility. Only by validating these models can we confidently assess and mitigate the potential impacts of the LAMP FE on beam delivery. Post-DTL, the beam splits to serve five major user facilities. Historically, low-energy beam transport has been modeled using TRACE, and higher-energy sections with TRANSPORT. These have now been unified into MAD-X format and validated with codes such as Elegant, pyOrbit, XSuite, Impact-Z, and HPSim. The primary focus now is on accurate modeling of full particle distributions (including beam halo) as they traverse the accelerator and beamlines to each experimental station. All models are at various stages of validation with empirical data.

43 PARTICLE ACCELERATORS↗