Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow analysis”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

ENSO-Induced Teleconnection: Process-Oriented Diagnostics to Assess Rossby Wave Sources and Ambient Flow Properties in Climate Models

Climate model fidelity in representing ENSO-induced teleconnection is assessed with process-oriented diagnostics that examine a chain of processes, from equatorial Pacific precipitation to the midlatitude circulation pattern over the Pacific–North American regions. Such processes are rarely addressed during model development. Using an upper-tropospheric divergent level, local vorticity gradient of the ambient zonal flow (∂ 2 U¯/∂y 2 ) and a restoring force for Rossby waves (β * ) are estimated, the equivalent barotropic vorticity equation is solved, and an anomalous Rossby wave source (RWS') quantified. The analysis is applied to AMIP5 and AMIP6 simulations. For a realistic circulation response representation, the hypothesis that models accurately represent the strength and location of RWS', and spatial variations in β * is tested. Compared to AMIP5, in AMIP6 there are clear improvements in representing RWS' and β * . To validate the hypothesis, the analysis identifies two metrics: spatially coherent RWS' in the subtropical North Pacific, and longitudes of negative β * over the western-central North Pacific. By projecting these metrics in two and three-dimensional views, improvements or degradations in model versions are apparent. If a model’s fidelity in representing ∂ 2 U¯/∂y 2 and RWS' are compromised, then radiated Rossby waves are reflected more equatorward, resulting in zonally elongated circulation anomalies over the central North Pacific. Thus, during climate model development, applying this analysis frequently will keep a regular check on the fidelity of the modeled response to anomalous El Niño convection in conjunction with changing model ambient flow dependencies. As a result, this analysis is intended to form a process-oriented diagnostics package, a community contribution to the NOAA Model Diagnostics Task Force.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Analysis of the linear and nonlinear stability of Alfven eigenmodes and fish-bones in JET DT discharges: mode identification and shear flows generation

The plasma in future nuclear fusion reactors will be heated by neutral beam injectors (NBIs) and high frequency electromagnetic waves as well as fusion born alpha particles. Energetic particles (EPs), with energies up to two orders of magnitude larger than the thermal plasma, can trigger EP driven modes and induce harmful EP losses, reducing the plasma heating efficiency and the economical viability of the reactor. The present study is dedicated to analyze the Alfven Eigenmode (AE) activity in JET D–T discharges, the closest experiment to reactor-like operation performed until now. There, EP driven modes are induced by the combined effect of tangential NBIs and ion cyclotron resonance heating (ICRH) driven EP. Linear and nonlinear simulations are performed with the gyro-fluid FAR3d code to analyze the AE activity observed in the discharge 99896. The linear simulations reproduce the unstable n = 3 to 5 toroidal AEs (TAE) at the inner plasma region observed in the experiment, triggered by highly energetic passing deuterium populations injected by the tangential NBIs, further accelerated by the effect of the ICRH up to 1 MeV. In addition, fish-bones triggered by energetic trapped hydrogen induced by the ICRH are also reproduced. On the other hand, the alpha particles density is too small to destabilize AEs in the experiment. Nonetheless, increasing artificially the alpha density by one order of magnitude, an n = 1 beta induced AE can be destabilized in the inner plasma region. Nonlinear simulations indicate the generation of zonal structures during the AE/fish-bone saturation phase. TAE and fish-bones causes a rather weak increase of the passing D and trapped H EP (around 2%), respectively. Shear flows and zonal currents are generated during the saturation of TAE and fish-bones. Nonlinear simulations performed for D–T and pure deuterium thermal plasma indicate AE/fish-bone activity is weaker and shear flows are less intense in the pure deuterium case, trends consistent with the experimental observations that also indicates a deterioration of the thermal plasma confinement. Therefore, both numerical studies and experimental evidence indicate the generation of shear flows by AE/fish-bones could be connected with an improvement of the thermal plasma confinement.

AE↗
DOE OSTI2025

Dimensionless Regime Mapping and Chaos Analysis of Solids Volume Fraction Data in a Circulating Fluidized Bed Riser

Presentation of experimental results of solids volume fraction data measurements in a 12-inch diameter cold flow circulating fluidized bed in relation to a dimensionless fluidization regime map. Statistical and chaotic system analysis results are shown to illustrate differences in hydrodynamic behavior at different operating conditions. Presentation was made at the ASTFE 10th Thermal and Fluids Engineering Conference, March 09-12, 2025.

chaos↗
DOE OSTI2025

Stability Analysis of Parallel Connected Bidirectional WPT System

This paper presents a stability analysis of parallel-connected bi-directional series-series resonant network wireless power transfer (WPT), optimized for Electric Vehicle (EV) charging and vehicle-to-grid (V2G) applications. The study addresses critical stability challenges in systems integrated with diverse distributed energy resources (DERs), including photovoltaics, fuel cells, wind turbines, energy storage systems, and the AC grid. The stability of such integrated DC grid systems is paramount for ensuring reliable operation, particularly under varying power flow conditions and dynamic interactions between parallel WPT systems. The analysis included system impedance characterization, state-space modeling, and open and closed-loop stability evaluations. The results demonstrated that the integration of a robust control architecture effectively mitigates instability risks and supports scalable, efficient operation. This work underscores the converter's adaptability and its potential for large-scale deployment in wireless EV charging infrastructures and integrated DC grid systems.

Asa, Erdem [ORNL] (ORCID:0000000190884812)↗
DOE OSTI2023

Exact solution of burst pressure for thick-walled pipes using the flow theory of plasticity

This paper develops an exact solution of burst pressure for defect-free, thick-walled pipes using the flow theory of plasticity in terms of the average shear stress yield criterion (or Zhu-Leis yield criterion). The pipe steel is assumed to obey the power-law strain hardening rule, and large plastic deformation is described by the finite strain theory. On this basis, internal pressure is obtained as a power series function of the effective strains on the inside and outside surfaces of the thick-walled pipe with use of the second Bernoulli numbers. At burst failure, the Zhu-Leis flow solution of burst pressure is determined as a power series solution, and the burst effective strains, the burst effective stresses, and the burst pressure are functions of the diameter ratio (D o / D i ), strain hardening exponent (n), and ultimate tensile strength (UTS). Similarly, the Tresca and von Mises flow solutions of burst pressure are also determined as a power series solution in terms of the Tresca and von Mises yield criteria. A general closed-from exact solution of burst pressure is then proposed for the three yield criteria, and the results showed that the proposed exact solution matches well with the power series solution for each yield criterion. Moreover, the von Mises flow solution is an upper bound prediction, the Tresca flow solution is a lower bound prediction, and the Zhu-Leis flow solution is an intermediate prediction that agrees well with the finite element analysis results of burst pressure for thick-walled pipes. Furthermore, two datasets of full-scale burst tests are then utilized to evaluate and validate the proposed flow solutions of burst pressure for both thin and thick-walled pipes.

42 ENGINEERING↗
DOE OSTI2024

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Benchmarking FFTF LOFWOS Test# 13 using SAM code: Baseline model development and uncertainty quantification

The development and deployment of advanced reactors, such as the sodium-cooled fast reactor (SFR), relies on sophisticated modeling tools to ensure the safety of the design under various transients. The predictive capability of these advanced modeling tools requires validation to garner trust in supporting the licensing of the advanced reactors. For this reason, the International Atomic Energy Agency (IAEA) initiated a coordinated research project (CRP) in 2018 for the analysis of the Fast Flux Test Facility (FFTF) Loss of Flow Without Scram (LOFWOS) Test #13.In this study, we present and discuss the benchmarking efforts of the modern system code SAM on the FFTF LOFWOS Test #13. Further, the SAM baseline model was developed according to the benchmark specification, which included a detailed core model with reactivity feedback. Generally, good agreement was observed between the baseline results and benchmark measurements; however, discrepancies persisted, particularly in predicted fuel assembly coolant outlet temperatures. Utilizing the baseline model, uncertainty quantification (UQ) and sensitivity analysis (SA) were conducted with the assistance of various statistical learning and machine learning methods, including kernel density estimation, Gaussian processes, and Sobol indices. Following the baseline model prediction and UQ and SA results, we discuss the reasons for the simulation discrepancies and propose further improvements to the model. This benchmarking effort adheres to the best-estimate plus uncertainty approach and can serve as a valuable example for supporting risk-informed licensing of advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Design and Analysis of Liquid Lithium Plasma Facing Components

In this study, liquid metal (LM) on plasma-facing components (PFCs) creates a renewable protective cover enhancing power exhaust and protecting the solid substrate. Additionally, energy confinement improvement through particle pumping can be achieved in the case of liquid lithium (LL). A new variant of PFCs was recently introduced at Princeton Plasma Physics Laboratory (PPPL) where a porous wall is used to stabilize the LM surface, while using magnetohydrodynamic (MHD) drive to push the LM flow inside the component. This arrangement allows efficient heat exhaust, and enhanced control of the LM surface temperature, leading to spatial control of evaporation and condensation of LL on the plasma interface. This feature is particularly attractive when vapor shielding is introduced to allow heat flux redistribution. This system has the advantage that as the heat flux increases, the evaporation rate will increase, while decreasing heat flux will decrease the evaporation rate of lithium, ideally creating a feedback effect which could self-regulate the amount of lithium evaporated. Analytical and numerical models for LL PFCs were developed in PPPL. To calculate the target temperature distribution for the case of evaporation from the divertor, we apply an iterative process allowing two-way coupling between the fluid-kinetic analysis of plasma using SOLPS-ITER code and the flow and heat transfer analysis of the PFC using an analytical model. In the last stage, the results are validated using computational fluid dynamics (CFD) analysis with customized version of the CFX code from ANSYS. CFX was modified at PPPL to allow MHD analysis at the high magnetic field typical for fusion applications. Results of the analysis for NSTX-U tokamak conditions will be presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Deep Learning for Subsurface Flow: A Comparative Study of U‐Net, Fourier Neural Operators, and Transformers in Underground Hydrogen Storage

Subsurface flow research is essential for the sustainable management of natural resources and the environment. Deep learning (DL) has significantly advanced this field by developing efficient and accurate surrogate models to replace computationally expensive physics‐based simulations. These surrogate models are commonly used to predict the spatiotemporal evolution of state variables, such as gas saturation and reservoir pressure, in heterogeneous geological formations. Despite the various DL models applied to this task, there is a lack of studies systematically comparing their performance. This absence of comparative analysis leads to somewhat arbitrary DL model selection in subsurface flow research, resulting in suboptimal performance and potentially inaccurate predictions. To bridge this gap, we conduct a systematic comparison study of three popular DL architectures—U‐Net, Fourier Neural Operators (FNO), and Segmentation Transformer (SETR)—in surrogate modeling of underground hydrogen storage (UHS). We focus on UHS due to its promise of enhancing clean energy resilience and its cyclic operational conditions that represent common scenarios in various subsurface applications. We evaluate the models based on accuracy, training cost, and inference speed. The comparison shows that U‐Net achieves the highest accuracy, followed by SETR and FNO. Despite its lower accuracy, FNO has the highest inference speed. SETR offers competitive accuracy with the least training memory usage, demonstrating the potential of transformers in learning subsurface flow. Our results provide guidance for selecting DL models for surrogate modeling in a wide range of subsurface flow problems.

42 ENGINEERING↗
DOE OSTI2023

Multi-Sourced Collaboration for the Production and Refining of Rare Elements and Critical Metals (Final Technical Report)

The project objective was to develop a feasible and cost-effective method for recovering rare earth elements (REEs) and critical materials (CMs) from coal and coal byproducts, resulting in high-purity individually separated REEs and CMs. The targeted REEs included Y, Pr, Nd, Gd, Dy, and Sm, with a purity of over 99.5%, while the CMs included Co, Mn, Ga, Sr, Li, Ni, Zn, and Ge, with a purity of over 90%. The project aimed to design a prototype facility capable of producing 1-3 tonnes/day of high-purity REO mixes. The work was divided into four designated circuits: 1) REE extraction and concentration, 2) REE separation and purification, 3) RE metal production, and 4) CM production. To achieve these goals, the project involved 11 tasks, including technology reviews, research, process flow diagram development, mass balance estimation, and preliminary technical-economic analysis. The project team included researchers from the University of Kentucky, University of Alabama and Virginia Tech as well as process specialists from Argonne National Laboratory. MP Materials provided technical support regarding rare earth markets and processing while Alliance Coal performed resource assessment. The project included a market analysis for Nd/Pr, Tb, Dy, Gd, Y, Co, Mn, Li, Sr, Ga, Ni, Zn, and Ge. These analyses provided insights into the supply and demand trends as well as historic and future projections of market price relative to purity requirements for these elements. Two coal resources were selected for the project: the West Kentucky No. 13 (Baker) Seam and an undisclosed lignite resource in the Illinois coal basin. The estimated quantities of REEs in these resources were calculated based on production samples and drilling data. It was estimated that there is adequate supply for an operation producing one metric ton daily of higher purity mixed rare earth oxides (MREO) for approximately 20 years at a site located in western Kentucky. In Circuit 1, project data was obtained from a pilot heap leach and REE concentration facility. It was concluded that the existing circuit, which generated a MREO concentrate, two types of CM mixed products, and Li- and Sr-containing waters, would be suitable feed for circuits 2-4. Data from the first-of-its-kind coal coarse refuse heap leach pilot pad played a crucial role in estimating reliable elemental concentrations of the pregnant leaching solution (PLS). The average total REE concentration in the PLS was found to be 28.6 ppm. In Circuit 2, several concepts were explored including a novel process referred to as solvent-assisted chromatography (SAC). This concept involved a novel columnar reactor that incorporated multiple mixer/settlers, thereby enabling the operation of counter-flowing aqueous and organic phases. Unfortunately, due to project time constraints, a complete fundamental modeling analysis could not be completed to fully evaluate the technology. Molten salt electrowinning was considered as an alternative for circuit 3 following circuit 2 purification circuit utilizing the novel SAC process. A mass and energy balance of Nd reduction to metal in a fluoride containing molten salt electrolyte was conducted. Comparisons were made with the current state of Asian molten salt electrorefining, and potential improvements in siphoning rare earth metals (REM) from the reactor were presented. A cost estimate was performed for the production of 1 tonne per day, which yielded a total of $2.29 million for the nine electrowinning (EW) cells required. The selected option for circuits 2 and 3 was a plasma distillation process, which initially separates rare earth elements (REEs) from other elements. This is followed by selective electrowinning in various ionic liquids. The selection was made on the basis of thermodynamic modeling and experimental data previously published by a project partner. The combination offers an innovative approach to integrated refining and RE metal production. For Circuit 4, an extensive literature review was conducted for the processing of the CMs. The ultimate decision was to utilize a combined plasma and ionic liquid process as well to produce individual high-purity concentrates of Zn, Ni, Co, Mn, and Mg. A separate flowsheet for Li and Sr was recommended, which would yield carbonates of these elements. Due to the lack of suitable experimental data at this time, a process recommendation could not be provided but several methods have been proposed for consideration. Lastly, a techno-economic analysis (TEA) was conducted to assess the effectiveness of the proposed process for further investigation. The TEA results revealed a capital expense (CapEx) of $737 million and an annual operational expense (OpEx) of $220 million. Due to the selected elements, the hypothetical heap leach pad can produce 1 metric tonne per day of REO equivalent, but a conscious decision was made to only treat targeted REEs, resulting in the production of 0.4 metric tonne of REM. An estimated annual revenue of $90.87 million was projected based on standard market pricing information provided by the funding agency. During the TEA, ten different modules were evaluated for costing purposes. The precipitation circuit was identified as the largest single operational expense, followed by the Mg/Mn process due to the amount of treated metal. In terms of capital expenditures, the heap leach process incurred the highest cost, followed by the Mg/Mn process. The scalability of the plasma process is a crucial consideration since the reactors cannot be scaled beyond the largest demonstrated size due to their reliance on surface area of the slag and vapor phase. The purity estimate for the REEs are generally 98%±2% to produce a metal. The purity level being lower than the project objective was due to the lack of specific experimental data needed to tighten the tolerance of the estimates. Based on literature and previous experience, the CMs are estimated as follows; Ga (95%+, metal), Sr (95%+, carbonate), Li (95%+, carbonate), Ni (98%±2%, metal), Zn (95%+, metal sponge), Ge (95%+, metal), Co (98%±2%, metal), and Mn (98%±2%, metal).

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

A multi-dimensional dataset for two-phase instability in low pressure natural circulation based on direct transient local measurement

Thermal hydraulic instability can challenge the control and the safety of nuclear reactors. This experimental research captures and investigates the instability in a low-pressure (1–10 atmospheres) natural circulation loop with an adiabatic chimney installed above the heated section. Direct transient local measurement of two-phase flow is performed at multiple axial locations along the flow channel with thermocouples and conductivity probes traversed by customized high-precision mechanisms. Based on the periodicity of oscillations, a post-processing method is proposed to align the data captured non-simultaneously to reconstruct the instability behavior in one representative period. The averaging operation over multiple oscillation cycles is performed to reduce the random errors. Cyclic phenomena are clearly demonstrated by the post-processed data, and rapid void fraction change induced by the flashing is observed to be significant in the driving mechanism for low pressure oscillations. Parametric effects of inlet subcooling, heat flux, inlet flow restriction and system pressure are studied within the experiment range, and qualitative results are explained and compared with existing literature. The ratio between the oscillation period and the travelling time of the fluid particle through the test section is compared with the theoretical value for the density wave oscillation. The new dataset contains directly measured transient multi-dimensional two-phase data along the flow channel for natural circulation under density wave oscillations. The generated dataset is valuable for the validation of system analysis codes and CFD codes in the simulation of natural circulation phenomena and dynamic two-phase flow.

42 ENGINEERING↗
DOE OSTI2022

Catalytic Activation of Polyethylene Model Compounds Over Metal-Exchanged Beta Zeolites

Decomposition of polymers by heterogeneous catalysts presents a promising approach for reuse of waste plastics. Here, we demonstrated non-hydrogenative decomposition of model polyolefins over proton-form and metal (Cu, Ni) ion-exchanged beta zeolites at moderate temperatures (ca. 300 degrees C). Near complete polyolefin decomposition was observed in batch reactions monitored by thermogravimetric analysis, while decomposition at partial conversion was studied in flow reactions. Ni-exchanged zeolites produced H 2 at substantially higher rates (> 10x) than other catalysts while also uniquely resisting deactivation over time. Application of the delplot formalism offered insights into the reaction network for polyolefin decomposition over Ni/BEA, most notably that H 2 is solely a primary product. We deduce that H 2 production is catalyzed by activation of C-H bonds at ionic Ni sites, and H 2 prevents buildup of polyaromatic coke species in Ni-exchanged zeolites that deactivate Cu-exchanged and protonic zeolites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Quantitative concentration measurements in a turbulent helium jet using rainbow schlieren deflectometry

Rainbow schlieren deflectometry (RSD) is applied to acquire quantitative concentration measurements in a turbulent environment for the first time. RSD methodology is developed and validated using Rayleigh scattering data available in the literature. Experiments were performed in an initially laminar, momentum-driven helium jet injected into ambient air. Full-field measurements are acquired that span the laminar, transition, and fully turbulent regions of the jet. Furthermore, the RSD results are shown to provide accurate measurements in each region and also confirm the self-similar nature of the fully turbulent jet. Analysis shows that the transition from laminar to turbulent flow occurs across an axial span that decreases with increasing Reynolds number. Experiments with a long exposure time were conducted to evaluate the limitations of hardware averaging of the RSD signal in comparison with software averaging.

02 PETROLEUM↗
DOE OSTI2020

Characteristics of flow through randomly packed impermeable and permeable particles using pore resolved simulations

Pore resolved simulations are performed to study the mean flow characteristics in porous media formed by random distribution of impermeable and permeable particles. Here, permeability of the medium and reactive surface areas are analyzed for a range of solid fractions, packing patterns and flow rates. The results are compared with models presented in the literature and with simulations of flow through ordered packings. Permeability of a randomly packed porous medium is found to depend on the heterogeneous distribution of porosity which results from particle agglomeration and is typically ignored by simpler models. The reactive surface area is dependent on the flow field at pore scale in addition to geometry. Moreover, the wake formed downstream of the particles at moderately high flow rates also affects the reactive surface area. Similar analysis is also performed by assuming particles to be permeable to account for dual scale porosity relevant for various engineering applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Interaction between dissolution and precipitation during olivine carbonation: Implications for CO 2 mineralization

Large-scale carbonation of olivine is considered a promising approach for in situ mineral carbonation, offering a permanent and stable method for CO 2 storage. A critical aspect of this process is understanding how dissolution and precipitation interact, as this could drive fracturing and enhance further reactions. In this study, we conducted carbonation experiments on olivine using CO 2 -saturated aqueous solutions of NaHCO 3 and NaCl. Two experimental setups were used: one representing an open geochemical system and the other a closed system, corresponding to reaction-limited and flow-limited scenarios, respectively. Further, post-reaction textural analysis using scanning electron microscopy (SEM) revealed surface coatings of reaction products in the closed system, while etch pits and etch channels were prevalent in the open system. Although no direct evidence of reaction-driven fracturing was observed, etch pits and etch channels may serve as initiation points for subcritical crack formation and growth, potentially maintaining permeability and exposing new unreacted surfaces. Using linear elastic fracture mechanics (LEFM) model, we estimate that microcracks could propagate under a pressure of 0.1 GPa if reaction products accumulate within the etch pits. Our findings offer new insights into the mechanisms governing olivine carbonation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗
DOE OSTI2026

Optimizing district energy systems by integrating Borehole Thermal Energy Storage Using a Mixed-Integer Linear Programming g-function framework with a Multi-Timescale Rolling Horizon method

Shallow geothermal has gained increasing attention in recent years; however, a reliable framework for its accurate incorporation into large-scale energy system optimization remains lacking. This study proposes a Mixed-Integer Linear Programming (MILP) framework combined with the g-function approach to integrate Borehole Thermal Energy Storage (BTES) technology into energy system optimization. Validation against a Modelica-based reservoir network simulation demonstrates that the proposed framework effectively captures the ground thermal response under varying energy loads and accurately estimates the borefield energy supply. To enhance scalability, a Rolling Horizon with Multi-Timescale (RH-MTS) method is further introduced, reducing computational time by 73 % for the 1-year optimization model with only minor loss of optimality. The framework is demonstrated through the case study of the UC Berkeley campus. Results indicate that BTES is a cost-effective and low-carbon solution: two borefields comprising 382 boreholes can meet 8.0 % and 6.6 % of the total campus heating and cooling demand, respectively, at an average energy rate of 0.70–0.77 USD/kWh and carbon intensity of 0.54 kg-CO2/kWh. Short-term analysis reveals a 35%–65% decline in BTES energy flow after 3–6 months of continuous heating/cooling operation, while long-term simulation shows that annual energy production of BTES can vary by up to 12.0 % after four years before stabilizing. Overall, this study develops a novel optimization framework that couples physics-based g-function method with MILP optimization framework, thereby advancing methodological development for shallow-geothermal integration and providing actionable guidance for BTES deployment in district-energy systems.

Yang, Jiahui↗
DOE OSTI2024

Development of a high-temperature Inconel 625 heat exchanger by model design and binder jetting additive manufacturing

A nickel-based Inconel 625 superalloy heat exchanger for high-temperature applications was developed via binder jetting additive manufacturing. The material properties were characterized first on printed and sintered parts. Two sintering temperatures were used to investigate the effects of the temperature on densification and microstructure. For channel geometry design, the heat transfer capabilities of the heat exchanger were optimized for the cross-section geometry of fluid flow channels in a counterflow configuration, and a stress analysis was conducted to investigate the effects of channel geometry. After headers were incorporated into the heat exchanger as a one-piece component, a prototype was printed. A complete depowdering was achieved via compressed air blowing, and sintering was performed with the developed profile. Ultimately, a one-piece heat exchanger with a nearly full density was obtained. Heat transfer tests were performed on this unit, and the results were compared with those from the simulations. In this study, systematic processes were developed for an additively manufactured Inconel-based heat exchanger for high-pressure, high-temperature heat transfer applications.

36 MATERIALS SCIENCE↗