Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data Reasoning”

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 865 records · Page 48

Towards Pareto optimal high entropy hydrides via data-driven materials discovery

The ability to rapidly screen material performance in the vast space of high entropy alloys is of critical importance to efficiently identify optimal hydride candidates for various use cases. Given the prohibitive complexity of first principles simulations and large-scale sampling required to rigorously predict hydrogen equilibrium in these systems, we turn to compositional machine learning models as the most feasible approach to screen on the order of tens of thousands of candidate equimolar high entropy alloys (HEAs). Critically, we show that machine learning models can predict hydride thermodynamics and capacities with reasonable accuracy (e.g. a mean absolute error in desorption enthalpy prediction of ~5 kJ mol H 2 –1 ) and that explainability analyses capture the competing trade-offs that arise from feature interdependence. We can therefore elucidate the multi-dimensional Pareto optimal set of materials, i.e., where two or more competing objective properties can't be simultaneously improved by another material. This provides rapid and efficient down-selection of the highest priority candidates for more time-consuming density functional theory investigations and experimental validation. Various targets were selected from the predicted Pareto front (with saturation capacities approaching two hydrogen per metal and desorption enthalpy less than 60 kJ mol H 2 –1 ) and were experimentally synthesized, characterized, and tested amongst an international collaboration group to validate the proposed novel hydrides. Finally, additional top-predicted candidates are suggested to the community for future synthesis efforts, and we conclude with an outlook on improving the current approach for the next generation of computational HEA hydride discovery efforts.

36 MATERIALS SCIENCE↗

Numerical simulation of involute-plate research reactor flow behavior using RANS, LES and DNS

This paper investigates the flow behavior of involute-plate research reactors by performing Reynolds-Averaged Navier Stokes simulation (RANS), Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) of the channel flow between fuel plates. By modeling turbulence with different numerical approaches, this study provides data with three levels of fidelity. For the RANS simulation, three widely used turbulence models, i.e., k-ε, k-ω, Reynolds Stress Turbulence model (RST) are applied by using the commercial CFD code STAR-CCM +. For LES and DNS, the open-source CFD code, Nek5000, is used given its outstanding scalability on High Performance Computer (HPC) and high-order technique. The results from RANS simulations are compared with that from LES and DNS for benchmarking. Both macroscale parameters and turbulence statistics, such as velocity magnitude, lateral velocity and turbulence kinetic energy, are presented and analyzed. The results from RANS simulation achieve good agreement with LES and DNS on velocity and turbulence kinetic energy prediction. The RST turbulence model predicts the most similar flow pattern of lateral velocity as compared to LES and DNS. The Lambda-2 (λ2) criterion with a reasonable threshold is used to demonstrate the instantaneous vortices distribution in the involute channel from both LES and DNS calculation. The DNS simulation captures more detailed turbulence especially near the corner, which explains the discrepancy between LES and DNS results near the corner. The normalized RMS error are defined and calculated to assess the performance of those turbulence models. The RST model captures the anisotropic feature of turbulence, which enable it to outperform other turbulence models for predicting the flow behavior in an involute channel. Although some discrepancies are found between LES and DNS results in the corner, the overall deviations between LES and DNS are found to be small. In conclusion, given that the computational cost of DNS calculation is an order of magnitude higher, using LES data for benchmarking RANS model is a cost-effective approach.

DNS↗

Refined views of ancient ocean chemistry: Tracking trace element incorporation in pyrite framboids using atom probe tomography

The trace element chemistry of pyrite can be used to determine the origin, timing, and conditions of formation of ore deposits; as a vectoring tool for mineral exploration; and to determine the evolution of the Earth’s oceanic and atmospheric chemistry. However, little is known about whether trace elements are held with the pyrite structure or within nano-inclusions of other phases. This distinction is important for two primary reasons. First, trace element incorporation into the pyrite structure can affect the partitioning of other trace elements. Second, if trace elements are held within nano-inclusions, the partition coefficient of the mineral phases that make up the nano-inclusion, rather than pyrite, are the critical consideration in related interpretations. Previous studies addressing this topic have used laser ablation inductively coupled plasma mass spectrometry; however, the resulting data do not provide sufficient spatial resolution to delineate trace element distributions unless the inclusions are large. Further, they use these flat element profiles in time resolved laser ablation output graphs to argue that pyrite trace element content provides a direct relationship between trace element content of seawater and pyrite trace element content. To improve resolution, we have used atom probe tomography to characterize trace element distributions in pyrite framboids from the Cariaco Basin and Demerara Rise in three-dimensions at sub-nanometer resolution. Manganese was found to be concentrated in the pyrite nanocrystal part of the original framboid structure. In contrast, Ni was mostly found along the grain boundary, though it still appeared to be contained within the pyrite structure. Copper was concentrated in later pyrite overgrowths, and As varied in its location. These observations suggest that some important trace elements are incorporated into pyrite during early diagenesis, even in euxinic settings dominated by water-column pyrite formation. Statistical analysis was used to determine whether trace elements were incorporated in the lattice or within nanoscale inclusions (referred to here are nano-inclusions). We found that As, Ni, Cu, and Mn were commonly held within the pyrite structure, but As, Ni, and Cu can also be held as nano-inclusions or within grain interfaces. Incorporation of As is known to enhance the incorporation of other trace elements and in this case appears to correlated to elevated Ni and Cu concentrations in the Cariaco Basin samples. Furthermore, understanding these relationships strongly impacts our ability to utilize pyrite trace metal concentrations to analyze and quantify early ocean chemistry and its evolution through time.

54 ENVIRONMENTAL SCIENCES↗

Dynamics of Fungal and Bacterial Biomass Carbon in Natural Ecosystems: Site-level Applications of the CLM-Microbe Model

Explicitly representing microbial processes has been recognized as a key improvement to Earth system models for the realistic projections of soil carbon (C) and climate dynamics. The CLM-Microbe model builds upon the CLM4.5 and explicitly represents two major soil microbial groups, fungi and bacteria. Based on the compiled time-series data of fungal (FBC) and bacterial (BBC) biomass C from nine biomes, we parameterized and validated the CLM-Microbe model, and further conducted sensitivity analysis and uncertainty analysis for simulating C cycling. The model performance was evaluated with mean absolute error (MAE), root mean square error (RMSE), and coefficient of determination (R2) for relative change in FBC and BBC. The CLM-Microbe model is able to reasonably capture the seasonal dynamics of FBC and BBC across biomes, particularly for tropical/subtropical forest, temperate broadleaf forest, and grassland, with MAE < 0.49 for FBC and <0.36 for BBC and RMSE <0.52 FBC and <0.39 for BBC, while R2 values are relatively smaller in some biomes (e.g., shrub) due to small sample sizes. We found good consistencies between simulated and observed FBC (R 2 =0.70, P<0.001) and BBC (R 2 =0.26, P<0.05) on average across biomes, but the model is not able to fully capture the large variation in observed FBC and BBC. Sensitivity analysis shows the most critical parameters are turnover rate, carbon-to-nitrogen ratio of fungi and bacteria, and microbial assimilation efficiency. This study confirms that the explicit representation of soil microbial mechanisms enhances model performance in simulating C variables such as heterotrophic respiration and soil organic C density. The further application of the CLM-Microbe model would deepen our understanding of microbial contributions to the global C cycle.

54 ENVIRONMENTAL SCIENCES↗

Fluid structure interaction of a subaqueous pendulum: Analyzing the effect of wake correction via large eddy simulations

The dynamic behavior of a subaqueous cylindrical pendulum and corresponding flow dynamics are investigated. The objectives were twofold: (i) to examine whether the two-dimensional model equations sufficiently capture the three-dimensional dynamics and (ii) to investigate the emerging three-dimensional vortical flow structures. Large eddy simulations with two-way coupling fluid structure interaction were carried out using the immersed boundary method to simulate the motion of the pendulum and its interactions with the initially stagnant water. The resulting pendulum motion is compared against measured data obtained in a series of experimental tests to validate the simulation results and the model equations with and without wake corrections. An analysis of the flow vorticity revealed the development of a vortex ring during the first swing and the formation of tip vortices. The evolution of the vortex rings emerging from the motion of the subaqueous cylindrical pendulum was visualized using Q-criteria showing a reasonable agreement with vortical structures observed in the experiment using particle imaging velocimetry. The hydrodynamic moments acting on the simulated pendulum and the moments calculated from the model equations are analyzed. Using the insights from these numerical simulations, a modification of the wake correction is proposed to enhance the accuracy of the rate of decay and period. The transient effect of coherent flow on pendulum dynamics, especially the added mass effect, is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale

Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.

42 ENGINEERING↗

Ionization detail parameters and cluster dose: a mathematical model for selection of nanodosimetric quantities for use in treatment planning in charged particle radiotherapy

Abstract Objective . To propose a mathematical model for applying ionization detail (ID), the detailed spatial distribution of ionization along a particle track, to proton and ion beam radiotherapy treatment planning (RTP). Approach . Our model provides for selection of preferred ID parameters ( I p ) for RTP, that associate closest to biological effects. Cluster dose is proposed to bridge the large gap between nanoscopic I p and macroscopic RTP. Selection of I p is demonstrated using published cell survival measurements for protons through argon, comparing results for nineteen I p : N k , k = 2, 3, …, 10, the number of ionizations in clusters of k or more per particle, and F k , k = 1, 2, …, 10, the number of clusters of k or more per particle. We then describe application of the model to ID-based RTP and propose a path to clinical translation. Main results . The preferred I p were N 4 and F 5 for aerobic cells, N 5 and F 7 for hypoxic cells. Significant differences were found in cell survival for beams having the same LET or the preferred N k . Conversely, there was no significant difference for F 5 for aerobic cells and F 7 for hypoxic cells, regardless of ion beam atomic number or energy. Further, cells irradiated with the same cluster dose for these I p had the same cell survival. Based on these preliminary results and other compelling results in nanodosimetry, it is reasonable to assert that I p exist that are more closely associated with biological effects than current LET-based approaches and microdosimetric RBE-based models used in particle RTP. However, more biological variables such as cell line and cycle phase, as well as ion beam pulse structure and rate still need investigation. Significance . Our model provides a practical means to select preferred I p from radiobiological data, and to convert I p to the macroscopic cluster dose for particle RTP.

Engineering↗

Demonstration of μCHP in Light Commercial Hot Water Applications

Internal combustion engine driven combined heat and power (CHP) systems produce power while the waste heat is recovered and used for another purpose. For µCHP systems (less than 50kWe output) overall efficiencies can reach over 90% for products currently on the market throughout the world. On a fuel basis this offers an advantage over using grid-based power and another fuel for heat. Additionally, the cost of natural gas on an energy unit basis is significantly less expensive than electric power. This difference can lead to substantial savings for the end user of a CHP system. In some parts of the world, particularly Germany and Japan, µCHP has become a common part of the heat and distributed generation equipment mix. Market analysis conducted in 2012 showed that in North America there is potential of a 5,300 units per year market for µCHP, resulting in an annual energy savings of approximately 13.5 TBtu and more than $84 million in annual customer savings. However, µCHP still has not seen much market penetration in North America. This study helps to identify some of the reasons why, address roadblocks that still need to be overcome and provide experience to established best practices. Nine commercial sites were initially identified (fitness club, surgical center, senior community, resort and spa, commercial laundry, two multi-family housing, restaurant and manufacturing facility). Systems from two different manufacturers were installed at five sites total, but due to unforeseen circumstances, only three of those sites were finally commissioned: test data were collected by Oak Ridge National Laboratory. While not all systems were installed as intended, the project resulted in the following key findings: Micro combined heat and power is an unfamiliar technology for building owners and plumbing and electrical contractors that serve the buildings most suited for µCHP. Manufacturers will need to be aggressive in providing training and providing broad support to successfully bring the technology to market.; The proper integration of µCHP into other building systems is the most critical part of a successful installation. Manufacturers will need to be able to provide the expertise to train, advise and support the mechanical contractors during the design and installation of the heat integration system.; Regulations governing µCHP varies across the country. There are a few UL standards (UL 2200 and UL 1741) that apply nationally. Emissions are governed by EPA except in California where the California Air Resources Board (CARB) has established their own standard. Requirements that dictate installation and interconnection are dependent on local building codes and local utility processes and requirements. Some out of date building codes do not offer clear guidance for these requirements, leading to confusion. This makes developing training and determining broadly applicable best practices difficult.; Installation costs of µCHP are a key factor in making an installation successful. It is critical that these costs and all national, state and local regulations be understood and accounted for by all parties involved in the installation before a project is undertaken.; A robust customer service organization is important to meet the needs of the customers. These systems require annual maintenance and service and local capability is needed to keep costs inline.

03 NATURAL GAS↗

Preliminary Study of Potential Utilization of Triassic Rift Basins for Long-Term Carbon Sequestration

Triassic rift basins of the eastern United States present a potential option for long–term carbon dioxide (CO 2 ) sequestration. Exposed and buried basins are located near large point CO 2 sources. Because of their similar origins, comparable fill successions indicate a level of reproducibility that can be conveyed between basins. The thick rock successions of the exposed Culpeper and Gettysburg basins served as proxies for understanding of the basin sequences. Their study allowed recognition of five mappable assemblages of rock types, herein termed lithofacies associations. These lithofacies associations were formed by alluvial fan, braided and meandering streams and marginal and distal lake depositional processes.Stratigraphic architecture of the lithofacies associations within exposed basins suggests a vertical succession that consists of an initial coarse-grained fluvial succession that is progressively replaced upward and basinward by finer grained lacustrine deposits. Along the faulted margins of the basins alluvial fans transitioning to fluvial delta deposits reflect lateral progradation and contemporaneous lateral components of the basin filling.To test the model devised in exposed basins, facies associations were developed for the buried Taylorsville Basin. This basin preserves over 8,000 feet of Triassic rocks that are concealed beneath more than 2,000 feet of Cretaceous and Tertiary Coastal Plain sediments. Composition and stratigraphic architecture of these rocks are similar to those observed in exposed basins. Infilling of the basin was the result of vertical aggradation and lateral progradation of coarse-grained to fine-grained facies. Based upon study of outcrops in the Culpeper and Gettysburg basins, groups of recurring lithologic facies were identified as the fundamental constructs of the basin’s sediment infilling. These groups of facies, termed lithofacies associations, represent an amalgamation of lithologic components from broadly similar depositional systems. These lithologic associations were then extrapolated into the exposed and buried portions of the Taylorsville basins. This effort determined that Triassic rift basins provide several avenues for potential study of geologic sequestration of carbon. The characteristic rift basin succession in the basins presents possible reservoir targets within the marginal fluvial deposits. Furthermore, intrusive and extrusive mafic bodies present a potential source of fracture porosity that could serve as CO 2 reservoirs. Porosity values from thin section analysis of exposed basin samples are higher in alluvial fan and braided fluvial lithofacies. Porosity is highest in samples that show less compaction, due to a lack of ductile lithic fragments and/or higher stratigraphic position. Log data from the Taylorsville Basin indicates that porosity and permeability values in basin marginal fluvial strata are elevated. Near the center of the basin, porosity and permeability values are greatly reduced, primarily owing to the fine-grained character of the lacustrine deposits. Some Triassic basins also contain thick intervals of concordant extrusive and intrusive mafic igneous rocks. These igneous bodies are potential CO 2 reservoirs for several reasons. Firstly, extrusive lava flows provide potential storage in the layers of primary porosity that occur in the vesicules formed at the top of lava flows. Secondly, both lava flows and subsurface igneous sills exhibit extensive fracture porosity produced by the rapid cooling of the flows and intrusions. Thirdly, these igneous rocks are mafic in composition, and studies have shown that iron- and magnesium-rich mafic rocks provide sequestration opportunities through carbonate remineralization. Lastly, extrusive and intrusive igneous rocks are invariably preserved within fine-grained lake deposits. These lacustrine sediments can serve as a fine-grained confining layer that encases the igneous rocks both above and below.

01 COAL, LIGNITE, AND PEAT↗

Autonomous Flow Electrochemistry for Accelerated Catalyst Discovery

Our objective is to develop an Autonomous Chemical Experimentation (ACE) platform that accelerates discovery of new catalytic transformations and other energy-relevant chemical reactions and processes. We intentionally designed ACE to be highly modular, both with respect to its rapid deployment to different chemistries and experimental workflows as well as incorporation of a wide range of different AI algorithms. In addition to the development of the core software architecture, initial efforts were made to incorporate Large Language Models to provide human-interpretable reasoning of the optimizer’s actions, and to develop a user-friendly graphical interface for experimental researchers. ACE was demonstrated using a flow electrocatalysis platform containing an inline FTIR spectrometer for real-time analysis and quantification of the reaction outcome. Human-in-the-loop experiments were performed in which a human researcher conducted an experiment using electrode potentials suggested by ACE, then fed the spectral data back to ACE for decision making. After confirming the successful function of the optimizer, efforts were next directed to automation of the hardware and performed full autonomy tests using three reactions: catalytic oxidation of formate, catalytic oxidation of cyclohexanol, and oxidation of hydroquinone. These studies confirm that ACE can close the loop between reaction execution, analysis, and optimization. They also reveal that more improved product detection methods will be essential for ACE to make well-informed decisions for reactions with low conversions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insights into Methodologies and Stochastic Optimization of Thermal Energy Storage-Coupled Advanced Reactor Systems: A Comparison of Methods for Accessing Long-Term Sub-System Sizing Adequacy

This paper investigates the potential of coupling Thermal Energy Storage (TES) with Advanced Reactors (ARs) to address uncertainties posed by climate change in deep decarbonized power systems. The TES Use-case Team at Idaho National Laboratory (INL) has examined the potential of storing thermal energy from ARs during low demand periods and optimizing discharge during peak-priced hours, in both steady-state and transient conditions. Building on this groundwork, this study bridges the gaps in optimal sizing of the sub-system of TES-coupled AR systems using Risk Analysis Virtual Environment (RAVEN) and Holistic Energy Resource Optimization Network (HERON), INL?s framework for grid optimization. By applying this framework, we present statistically-robust optimal charge, discharge including balance of plant (BOP), and storage sizing for the High-Temperature Gas-Cooled Reactor (HTGR) with 203 MWth output. To this end, we generated synthetic price samples for 30 years using 2018 ? 2021 real-time market data from ERCOT, PJM and MISO. Our results reveals that the TES-coupled HTGR system is highly effective in maximizing revenue from electricity sales. We observed a substantial increase of 40 % in ERCOT and a noteworthy 15 % increase in PJM and MISO when compared to the conventional BOP without TES. This improvement is achieved through regionally-tailored sub-system sizing, which ranges from 398 to 416 MWth for the discharge system and 610 to 1029 MWth for the TES. We find that the average electricity price directly impacts the overall economics, while price volatility influences storage size. Additional sensitivity analyses were performed to access the impact of key assumptions on system economics and sizing, differentiating the optimization window (i.e., 24 ? 219 hours of chronological observations) and by imposing storage continuity condition in tracking TES cycles. We observed that at the 120-hour of the optimization window, a reasonable balance between computation time and accuracy was achieved. Our analysis also highlights the significance of conducting multi-day cycle analysis (> 120-hour) for TES to capture interaction between electricity prices and storage dynamics, providing a comprehensive understanding of TES behavior that AR developers should integrate into their plant designs.

25 ENERGY STORAGE↗

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗

February 2020 Bayou Choctaw Subsidence Analysis

Subsidence monitoring is a crucial component to understanding cavern integrity of salt storage caverns. This report looks at the historical and current subsidence monitoring program and includes interpretation of the data from the Bayou Choctaw Strategic Petroleum Reserve site. The current monitoring program consists of an annual elevation survey as well as GPS and tiltmeter instruments above both Cavern 4 and Cavern 20. This year's level and rod survey indicates little subsidence across the site. In addition, the GPS and tiltmeter instruments do not indicate any substantial movement above caverns 4 and 20. As such, there is no reason to indicate any of the caverns at Bayou Choctaw have lost integrity.

02 PETROLEUM↗

Two-stage formation-energy correction (NbZr, TaZr, VZr)

This bundle contains the scripts, the raw and corrected per-structure data, and the manuscript plots for the NbZr / TaZr / VZr BCC binary formation energies and the associated RMSDs. Why a two-stage correction is necessary: The "raw" formation energy of every relaxed VASP configuration is computed in the usual way, FE_raw(c) = E_alloy(c) - sum_i x_i * E_pure_i , where E_pure_i are the per-atom total energies of the pure-element reference structures (Nb, Ta, V, Zr in the same BCC supercell, with identical INCAR / KPOINTS / PAW choices). With perfectly consistent reference runs the raw FE should vanish at the two pure-element endpoints (x = 0 and x = 1) by construction. In practice this does not hold for two reasons that are present in our dataset: 1. Reference-energy inconsistency (composition-dependent bias). Even with identical input parameters, the pure-element runs (stored in `corrected_DFT_pure_element_runs/`) differ slightly from the values that would be implied by the alloy runs at near-pure compositions (a few meV/atom). This bias is approximately linear in concentration, because the residual error in E_pure_Nb (or E_pure_Ta / E_pure_V) propagates into FE_raw(c) as (1 - x) * dE_pure_1, and the corresponding error in E_pure_Zr propagates as x * dE_pure_2. Left uncorrected, this produces a non-physical "tilt" of FE_raw(x) and shifts the entire FE-vs-x cloud away from zero at the endpoints. 2. Endpoint anchoring against the audited true endpoints. The strict endpoint values (FE_x0_meVatom, FE_x1_meVatom in `corrected_fe_strict_endpoints_20260518/strict_endpoint_check_20260518.csv`) were re-derived from an independent cross-check of the pure-element runs. After stage 1 removes the linear bias, the near-pure compositions in the alloy dataset still extrapolate to values that differ slightly from these audited endpoints — because stage 1 is fit from a few near-end alloy bins, not from the audited pure-element references themselves. The README.txt file discusses how these issues are addressed by the two-stage correction, and describes folder layout, pipeline summary, and how to re-run.

36 MATERIALS SCIENCE↗

Permeate fluxes from desalination of brines and produced waters: A reactive transport modeling study

The increasing interest in the use of membrane systems to desalinate inland brackish water, agricultural drainage, and industrially produced wastewater demands improved means of predicting desalination system performance under variable feedwater compositions. The interaction among water flow, solute transport, and chemical composition in these systems impacts permeate flux evolution. Here, an established multicomponent reactive transport simulator that accounts for these coupled processes is applied to compute osmotic pressure and permeate fluxes in reverse osmosis (RO) systems. The model is first validated by predicting permeate fluxes for a set of benchtop crossflow experiments subject to a range of feed flow rates and compositions, under fouling and non-fouling conditions. Results compare favorably with measured data that show that solutions with similar total dissolved solids concentrations but different compositions result in different permeate fluxes. The model is then applied to predict permeate fluxes from the desalination of produced waters using a commercial spiral wound RO module. For NaCl-dominant brines, at total dissolved salt concentrations (TDS) below about 70 g/L, permeate fluxes are inversely proportional to water mole fraction as the latter is a reasonable approximation of water activity (i.e. ideal mixing). In the case of Ca–Cl-, Na–CO3- and Na–SO4-dominant brines below about 70 g/L TDS, this relationship does not hold as well and tends to overpredict osmotic pressure and thus underpredict permeate fluxes. However, the opposite becomes true at higher TDS values for typical produced waters. The scaling potential of these waters is also computed by allowing the precipitation of minerals above their saturation limit on the RO membrane. This work demonstrates how reactive transport models developed for the analysis of waters from geological systems can be extended to improve process design, optimization, and control in desalination systems from produced waters and beyond.

Molins, Sergi↗

Multimuons in cosmic-ray events as seen in ALICE at the LHC

ALICE is a large experiment at the CERN Large Hadron Collider. Located 52 meters underground, its detectors are suitable to measure muons produced by cosmic-ray interactions in the atmosphere. In this paper, the studies of the cosmic muons registered by ALICE during Run 2 (2015–2018) are described. The analysis is limited to multimuon events defined as events with more than four detected muons (N μ > 4) and in the zenith angle range 0° < θ < 50°. The results are compared with Monte Carlo simulations using three of the main hadronic interaction models describing the air shower development in the atmosphere: QGSJET-II-04, EPOS-LHC, and SIBYLL 2.3d. The interval of the primary cosmic-ray energy involved in the measured muon multiplicity distribution is about 4 × 10 15 < E prim < 6 × 10 16 eV. In this interval none of the three models is able to describe precisely the trend of the composition of cosmic rays as the energy increases. However, QGSJET-II-04 is found to be the only model capable of reproducing reasonably well the muon multiplicity distribution, assuming a heavy composition of the primary cosmic rays over the whole energy range, while SIBYLL 2.3d and EPOS-LHC underpredict the number of muons in a large interval of multiplicity by more than 20% and 30%, respectively. The rate of high muon multiplicity events (N μ > 100) obtained with QGSJET-II-04 and SIBYLL 2.3d is compatible with the data, while EPOS-LHC produces a significantly lower rate (55% of the measured rate). For both QGSJET-II-04 and SIBYLL 2.3d, the rate is close to the data when the composition is assumed to be dominated by heavy elements, an outcome compatible with the average energy E prim ∼ 10 17 eV of these events. This result places significant constraints on more exotic production mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Pebble Tanker Model for Nuclear Criticality Safety Needs

This report documents a study performed to investigate the requirements for criticality safety benchmark experiments for high-assay, low-enriched uranium (HALEU) fuel in transportation applications. In this work, an exploratory application model, the “Pebble Tanker,” was developed to represent TRISO fuel in a transportation scenario for an analysis of the validation basis in industrial quantities. An aspect of the criticality validation process involves assessing the “similarity” between application and experimental benchmark systems through an integral index parameter evaluation. Here, this includes propagating nuclear data uncertainties and calculating a correlation coefficient (hereinafter referred to as “c k ”) to evaluate the similarity of benchmark experiments compared with the application Pebble Tanker model. Finding sufficient critical benchmark experiments allows for the evaluation of bias and bias uncertainty, thus determining the upper subcritical limit (USL) of the transportation package. A target k eff of ~0.94 was used in this work to establish appropriate modeling conditions, reflecting a reasonable estimate for a USL. Two container models were investigated: one with the Hermes-type pebble and one with the Pebble Bed Modular Reactor (PBMR)–type pebble. The models were simplified, considering only fuel, containment structure, and either water or air. This allows a focus on the underlying physics of applications involving TRISO fuel pebbles using the Pebble Tanker model. A crucial consideration is the transport package's ability to safely hold pebbles while flooded, maintaining subcritical conditions. Tools available in the SCALE 6.3.1 suite—the CSAS6-Shift, TSUNAMI-3D-Shift, and TSUNAMI-IP sequences—were employed for neutronics and sensitivity and uncertainty (S/U) analysis of the Pebble Tanker. Findings demonstrated sufficient available critical experiment benchmarks to perform a validation of the Pebble Tanker in the most reactive state, i.e., when the Tanker is flooded.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Detecting Quantum Critical Points of Correlated Systems by Quantum Convolutional Neural Network Using Data from Variational Quantum Eigensolver

Machine learning has been applied to a wide variety of models, from classical statistical mechanics to quantum strongly correlated systems, for classifying phase transitions. The recently proposed quantum convolutional neural network (QCNN) provides a new framework for using quantum circuits instead of classical neural networks as the backbone of classification methods. We present the results from training the QCNN by the wavefunctions of the variational quantum eigensolver for the one-dimensional transverse field Ising model (TFIM). We demonstrate that the QCNN identifies wavefunctions corresponding to the paramagnetic and ferromagnetic phases of the TFIM with reasonable accuracy. The QCNN can be trained to predict the corresponding ‘phase’ of wavefunctions around the putative quantum critical point even though it is trained by wavefunctions far away. The paper provides a basis for exploiting the QCNN to identify the quantum critical point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗