SEARCH · Engineering Papers

Results for “High-resolution numerical model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Improving Process Level Understanding of Boundary Layer Winds over the Northeast U.S. Shelf: The Third Wind Forecast Improvement Project (WFIP3)

The third Wind Forecast Improvement Project (WFIP3), a U.S. Department of Energy and National Oceanic and Atmospheric Administration sponsored investigation, sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the marine atmospheric boundary layer. WFIP3 focused on mesoscale and submesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high-use coastal zone, using a 3D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the marine atmospheric boundary layer over the ocean was done from an air-sea interaction flux tower and extended deployments of a large autonomous barge platform. Numerous critical forecasting phenomena were observed that are being evaluated within regional coupled and uncoupled modeling systems, including the National Oceanic and Atmospheric Administration's foundational High-Resolution Rapid Refresh forecast model.

Kirincich, Anthony↗

DOE OSTI2026

Final DOE-ASR Report for the Project “Using LASSO to bridge the gap between model and observations and to learn about atmospheric convection”

Atmospheric convection spans a wide range of spatial and temporal scales and involves complex interactions with the surrounding dynamic and thermodynamic environment, particularly over tropical continental regions. These processes remain a major source of uncertainty in weather and climate models, including persistent biases in the diurnal cycle of convective precipitation that directly affect estimates of climate sensitivity. Addressing these challenges requires the combined use of high-resolution observations and cloud-resolving modeling frameworks. In this context, the DOE Atmospheric Radiation Measurement (ARM) program’s Large-Eddy Simulation ARM Symbiotic Simulation and Observation (LASSO) activity provides a powerful platform that pairs comprehensive observations with numerical simulations to enable process-level understanding of atmospheric convection. Within this context, this Research and Development Partnership Pilot (RDPP) project was designed to initiate and expand DOE ARM/ASR research capacity at minority-serving institutions, while advancing scientific understanding of convective processes over the Amazon rainforest. Consistent with the RDPP mission, the project emphasized partnership development, training, and workforce capacity building alongside exploratory research activities. On the scientific side, the project produced two peer-reviewed journal articles, and one manuscript currently under review (see list in section 3.1). Together, these studies combine long-term ARM observations and cloud-resolving and convection-permitting modeling to investigate the environmental controls on the shallow-to-deep convective transition during the Amazon wet season. The results demonstrate the central role of early-day moisture preconditioning and large-scale dynamical forcing in regulating isolated deep convection, provide mechanistic insight into convective evolution, and establish physically informed modeling frameworks for future sensitivity experiments. These scientific outcomes are described in sections 2.1 to 2.3 and were disseminated in 8 conference presentations (see section 3.2) and 5 invited talks (see section 3.3), reflecting broad engagement with our community. Equally important, the project achieved its RDPP capacity-building objectives (see section 2.4). A sustained research partnership was established among the University of Maryland, Baltimore County (UMBC), Morgan State University (MSU), and Howard University (HU), and extended to include collaboration with Pacific Northwest National Laboratory (PNNL). The project organized multiple multi-day training events focused on ARM data, LASSO simulations, and quantitative analysis methods, directly engaging students, postdoctoral researchers, and faculty across institutions. These activities broadened participation in ASR research and led to independent adoption of LASSO workflows by students beyond the immediate project team. Finally, the project successfully positioned the participating institutions to pursue future DOE research. Preliminary scientific results, coupled with strengthened partnerships and technical capacity, enabled the submission of follow-on proposals to DOE ASR funding opportunities. In this way, the project fulfilled the RDPP goal of seeding durable research capacity and laying the foundation for larger-scale, sustained engagement with DOE ARM and ASR programs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

High-Resolution ESM Projections for Energy Applications Over the CONUS

Assessing energy resources under future scenarios requires high-resolution meteorological information that is physically consistent and suitable for regional-scale analysis. While Earth system model (ESM) projections provide valuable large-scale information, their coarse resolution and systematic biases limit direct applicability for energy system modeling and planning. In this study, we develop a high-resolution dynamical downscaling framework based on the Weather Research and Forecasting (WRF) model to translate global-scale ESM data into energy-relevant regional projections over the contiguous United States (CONUS). The framework identifies an optimized WRF configuration through numerical experiments and evaluates raw and bias-corrected ESM initial and boundary conditions, with soil moisture (SM) and soil temperature (ST) bias correction implemented as an integral part of the bias-corrected ESM forcing to improve land-atmosphere coupling prior to WRF dynamical downscaling. Using an optimized WRF configuration at 4-km resolution, we show that raw ESM forcing introduces systematic dry and cold soil biases that propagate into pronounced warm biases in near-surface air temperature and positive biases in solar irradiance, particularly during summer. Applying bias-corrected atmospheric forcing together with bias-corrected SM and ST substantially reduces these downstream biases and improves the surface energy balance and near-surface atmospheric fields. These results demonstrate that bias-aware treatment of initial conditions is critical for producing high-resolution downscaled projections suitable for energy system modeling and planning applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

Dynamical Downscaling of Earth System Model Data for Energy System Analysis

Assessing energy resources (e.g., solar, wind, and hydro) under future scenarios requires datasets with sufficient spatial and temporal detail to capture variability and extreme events. While global-scale Earth System Model (ESM) projections are widely used, their coarse resolution limits direct application to regional energy system analyses. Dynamical downscaling offers a robust approach to generate physically consistent, fine-scale datasets that better represent local atmospheric processes impacting energy resources. In this work, we present a two-stage approach for producing high-resolution historical and future projections over the contiguous United States (CONUS). First, we optimize the Weather Research and Forecasting (WRF) model configuration for energy-relevant variables - solar irradiance, wind speed, and precipitation - by conducting ERA5-driven simulations at 8-km and 28-km resolution. Multiple physics schemes and model configurations within the WRF are evaluated against observational datasets including the National Solar Radiation Database (NSRDB), the Parameter-elevation Regressions on Independent Slopes Model (PRISM), and the Stage IV multi-radar/multi-sensor precipitation product for the CONUS domain. Using the best-performing configuration, we dynamically downscale MPI-ESM1-2-HR simulations for 2000-2060 under SSP2-4.5 and SSP5-8.5 scenarios at 4-km spatial and hourly temporal resolution. This presentation will provide a comprehensive analysis of the results from multiple numerical experiments and high-resolution ESM projections. In addition, we will discuss potential applications of our high-resolution datasets within the energy sector and outline future research avenues dedicated to evaluating how extreme weather events influence system performance and resilience.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

Mesoscale Organization in Cumulus-Coupled Stratocumulus

Marine cloud systems cover a substantial portion of the world’s oceans. Most of these clouds form relatively close to the ocean surface, typically within one to two kilometers, a region referred to by meteorologists as the marine boundary layer. They are composed predominantly of liquid water, although ice particles can occur in mid- and high-latitude marine clouds during winter. In satellite imagery, these clouds appear bright against the darker ocean surface below, reflecting a large fraction of incoming sunlight back into space that would otherwise warm the ocean. Because marine boundary layer clouds cover such an extensive area of the ocean, they exert a significant influence on Earth’s overall transfer of solar energy absorbed by the surface and thermal energy emitted to space, a balance known as the planetary radiation budget. Marine boundary layer clouds are typically thin, and their formation and dissipation depend on a delicate balance between processes acting at the ocean surface below and the warm, dry air above. They are notoriously difficult to simulate accurately in weather forecast models, which often produce too few marine low clouds in midlatitudes and clouds in tropical regions that are excessively bright, meaning they reflect too much solar radiation. The marine boundary layer is frequently characterized by widespread overcast cloud cover that often transitions from a continuous, single-layer deck to more broken cloud fields toward the tropics. These transitions typically proceed through an intermediate stage in which shallow, broken clouds form beneath the overlying stratiform cloud deck. Once broken clouds develop below the overcast, they frequently self-organize into cloud clusters known as marine boundary layer convective complexes (MBLCCs), although the mechanisms governing the formation and organization of MBLCCs remain poorly understood. Accurately representing these transitions in long-range weather forecast models is essential because they influence the properties of air masses advected over the continental United States and Europe, and they become increasingly important for forecasts on seasonal and longer timescales. We employed two complementary approaches to investigate the processes controlling MBLCCs and their impact on marine cloud cover. Long-term observations from the U.S. Department of Energy’s Eastern North Atlantic (ENA) Observatory provided a unique dataset that allowed us to characterize fundamental properties of MBLCCs, including their typical size and frequency of occurrence. These observations were combined with high-resolution numerical simulations performed on supercomputers to examine the evolution of MBLCCs during cold-air outbreaks over the ENA region.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Geothermal Heat Pump System Showcase: Short-Term Validation of Borehole Heat Exchanger Performance from Field Data to Numerical Modeling: Preprint

Since 2011, a geothermal heat pump (GHP) system has been operating to provide space heating and cooling for the Solar Radiation and Research Laboratory building at the National Laboratory of the Rockies (NLR) in Golden, Colorado. The system consists of 23 vertical boreholes, each extending to a depth of 300 ft (91 m), connected to 11 water-to-air heat pump units and four circulation pumps. Between fiscal years 2023 and 2025, additional power meters and temperature sensors were retrofitted to support detailed system performance assessment and model development. This study presents preliminary monitoring results and the development of an initial numerical model of the borehole heat exchanger field. The model incorporated site-specific geometry, ground thermal properties derived from thermal response tests, and ambient temperatures, and simulated system behavior over a representative operating day in September. Model predictions of outlet temperatures were compared against corresponding field measurements. Results showed that modeling initialized with a simplified linear subsurface temperature gradient presents systematic discrepancies in outlet temperature, whereas incorporating depth-resolved borehole temperature measurements for initialization yields substantially improved agreement with observations. The findings highlight the sensitivity of short-term predictive modeling to the representation of initial subsurface thermal conditions and underscore the value of high-resolution field measurements for model calibration and validation. These preliminary results inform ongoing efforts to extend the modeling framework to longer time horizons and to refine monitoring and modeling strategies that support the design guidance and operational optimization of GHP systems in research and commercial buildings.

15 GEOTHERMAL ENERGY↗

DOE OSTI2026

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗

DOE OSTI2026

Hierarchical-embedding autoencoder with a predictor as efficient architecture for learning time-evolution in multi-scale turbulent flows

We introduce a scale-aware, data-driven deep learning modeling framework for accurately predicting the time evolution of multi-scale turbulent plasma and liquid flows. The approach is motivated by the idea of scale separation. Structures of vastly different length scales emerge in these systems, and interactions between these structures occur only locally. To exploit this structure, the flow state is transformed by a hierarchical, fully convolutional autoencoder, not into a single embedding layer as in conventional convolutional surrogate models, but into a series of embedding layers. A stepwise training strategy ensures that fine-scale features are encoded on a high-resolution grid, while larger structures are represented on progressively coarser layers. The time evolution predictor advances all embedding layers in sync, capturing local interactions between features at the same scale as well as between all scales. This approach enables efficient modeling of multi-scale systems since negligible interactions between distant, small-scale structures do not need to be directly modeled. Our hierarchical-embedding autoencoder with a predictor framework is evaluated on canonical examples of multi-scale turbulence: two-dimensional Kolmogorov flow and Hasegawa–Wakatani plasma turbulence. In both cases, the proposed framework significantly improves predictive accuracy relative to conventional convolutional network architectures. A significant improvement in prediction accuracy was observed for crucial statistical characteristics of the Hasegawa–Wakatani plasma as well as for individual trajectories of the Kolmogorov flow turbulence. Importantly, the model's rollout for the Hasegawa–Wakatani problem demonstrates a four-order-of-magnitude speedup compared to traditional numerical solvers.

Khrabry, Alexander I. [Princeton Univ., NJ (United↗

DOE OSTI2026

Assessing shellfish water exposure to fecal bacteria pollution in Salish Sea: three-dimensional modeling and implications for monitoring

Fecal bacteria (FB) contamination poses significant risks to shellfish safety and management in coastal and estuarine waters. Despite extensive pollution identification and correction efforts, FB contamination in shellfish-growing areas persists in the Salish Sea, highlighting the need to identify overlooked sources and better understand FB transport from riverine and shoreline inputs to shellfish beds. To address this, a high-resolution three-dimensional hydrodynamic model coupled with FB kinetics was developed and applied to a case study site in Salish Sea—Portage Bay—to simulate freshwater plume circulation, flushing dynamics, and bacterial transport. Daily FB loading from the major freshwater inflow—Nooksack River was generated by both linear interpolation and integrating a machine learning approach (XGBoost), trained on historical hydrological and meteorological data. The model successfully reproduced both the magnitude and seasonal variation of FB concentrations in Portage Bay for the year of 2021, demonstrating that simplified FB kinetics with first-order decay due to mortality was effective in this dynamic coastal environment with short flushing time. Model results identified the Nooksack River as the dominant far-field FB source, while scenario simulations showed that near-field coastal stormwater outfalls elevated local FB levels following rainfall, particularly under low-flow conditions. The XGBoost prediction provided comparable or superior accuracy to linear interpolation, particularly during periods of missing observational data, by capturing short-term variability and event-driven loading more effectively. Integrating data-driven riverine FB inputs with mechanistic coastal numerical modeling provides a robust framework for operational forecasting of shellfish bed exposure risk and supports adaptive monitoring and management of shellfish growing areas in the Salish Sea and similar coastal systems.

Salish Sea↗

DOE OSTI2026

High-Resolution South American Wind Resource Data Downscaled with Generative Machine Learning Conditioned on Near-Surface Observations

High-resolution historical wind data was developed for the entirety of South America using the innovative Super-Resolution for Renewable Resource Data (sup3r) machine learning framework. The publicly available Sup3rWind South America dataset represents a significant advancement in wind resource data generation, leveraging generative machine learning conditioned on near-surface observations from the Meteorological Assimilation Data Ingest System (MADIS) to efficiently and accurately downscale coarse reanalysis data from the European Centre for Medium-Range Weather Forecasts (ERA5). This approach produces fine-scale, spatially and temporally coherent wind and meteorological fields hundreds of times more computationally efficient than traditional numerical weather modeling methods, enabling access to high-fidelity wind information across both continental and offshore regions. Sup3rWind South America builds on the earlier Sup3rWind Ukraine dataset through improvements in model architecture and outputs conditioned on near-surface observation inputs. As with the Ukraine data release, this dataset includes wind speed, wind direction, temperature, relative humidity, and pressure at a horizontal resolution of ~2 km, representing a 15x spatial enhancement relative to the 31 km ERA5 grid. Wind speed and direction are provided at 5-minute resolution, a 12x temporal refinement compared to the hourly ERA5 data, while temperature, relative humidity, and pressure remain at hourly resolution. The data covers all years from 2005 to 2024. Before downscaling, ERA5 inputs were bias-corrected using long-term monthly means and a limited number of quality-controlled observations to align large-scale statistics with regional conditions. The resulting dataset is the first publicly available high-resolution timeseries wind record that provides full spatial coverage of South America. Model validation demonstrates strong agreement with observations across several statistical metrics, consistent with other state-of-the-art high-resolution wind resource datasets. The potential applications of Sup3rWind South America span renewable energy resource assessment, energy system modeling, and grid resilience analysis. The 20-year record and high spatial and temporal resolution support accurate estimation of long-term energy yield and the economic feasibility of potential wind development sites. Continuous coverage across both continental and offshore regions enables comprehensive site prospecting within exclusive economic zones. The 2 km, 5-minute resolution data provide the spatial and temporal variability required for power system simulation, operational planning, and regional risk assessments.

17 WIND ENERGY↗

NASA NTRS2004

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda↗

DOE OSTI2026

CoURAGE KAZR b1 Processing: Corrections, Calibrations, and Processing Report

The U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) User Facility supports atmospheric and earth system research through a comprehensive network of fixed and mobile observatories. These facilities provide long-term and intensive campaign-based observations of clouds, aerosols, precipitation, radiation, and meteorological state variables. ARM observations are designed to improve the physical understanding and numerical representation of atmospheric processes in earth system models, with particular emphasis on cloud-radiation interactions and precipitation processes. The Coast-Urban-Rural Atmospheric Gradient Experiment (CoURAGE) deploys one of the ARM Mobile Facilities (AMF) to the Mid-Atlantic region surrounding Baltimore, Maryland, for the period 1 December 2024 through 30 November 2025. This deployment focuses on characterizing atmospheric structure, cloud properties, and precipitation processes across strong land-use and surface heterogeneity gradients associated with urban, rural, and coastal (Chesapeake Bay) environments. The CoURAGE deployment complements the Baltimore Social-Environmental Collaborative (BSEC), a DOE Urban Integrated Field Laboratory (UIFL), by providing high-quality atmospheric observations needed to connect urban surface processes, emissions, and meteorology to cloud and precipitation responses. In addition to the central urban site, ancillary observing sites were deployed to rural Maryland northwest of Baltimore and to an island site in Chesapeake Bay. These measurements further complement a long-term atmospheric observatory operated in Beltsville, Maryland, by Howard University in collaboration with the Maryland Department of the Environment. Together, these assets form a four-node regional atmospheric observatory network representing Baltimore and its three primary surrounding environments—urban, rural, and coastal/bay. This coordinated observational strategy enables investigation of spatial gradients in boundary-layer structure, cloud occurrence, precipitation evolution, and aerosol-cloud interactions across complex surface regimes. Within this network, vertically pointing cloud radars play a critical role by providing continuous, high-resolution measurements of cloud and precipitation vertical structure.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

High-Burnup LOCA Burst Susceptibility BISON Analysis in PWRs and BWRs

Accurately assessing high-burnup fuel behavior during loss-of-coolant accidents (LOCAs) is essential for understanding fuel fragmentation, relocation, and dispersal (FFRD) risks across the US light-water reactor fleet. This work updates previous Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program multiphysics LOCA analyses for a pressurized water reactor (PWR) and a boiling water reactor (BWR) by incorporating recent model and material property advancements in the BISON fuel performance code, including a high-burnup structure (HBS) model, revised cladding burst criteria, and updated thermal–mechanical correlations. This update was needed to support ongoing industry initiatives and upcoming regulatory changes. Full-core, rod-resolved operating histories generated using Virtual Environment for Reactor Analysis (VERA) and system-level LOCA conditions obtained from TRACE were applied to statistically representative rod samples in BISON to evaluate burst behavior and FFRD susceptibility. These calculations used two cladding burst correlations and three fuel pulverization models so that the predictions of these models could be compared. The updated PWR simulations show markedly improved numerical stability as the number of crashed simulations decreased by 95% compared to the previous study, and hence higher confidence in results. The updated PWR simulations predicted cladding bursts exclusively among once-burned, high-power rods, with two different cladding burst models identifying the same burst-susceptible population. Resulting FFRD susceptibility estimates are significantly reduced compared with earlier studies, driven by cooler predicted fuel and plenum temperatures, lower hoop strains, and reduced fission gas release in the updated models. In contrast, none of the BWR rods were predicted to burst under either burst criterion, reaffirming minimal BWR FFRD susceptibility even with updated HBS and material models. Comparisons between the PWR and BWR end-of-cycle predictions are made. Comparison with prior work highlights significant shifts in PWR fuel performance metrics and confirmation of earlier BWR conclusions. Overall, the updated results underscore the importance of having high-resolution detailed modeling capability and continuously integrating evolving material models and physics into high-resolution multiphysics simulations. The unified assessment presented here strengthens confidence in predicting high-burnup LOCA behavior by improving agreement between different cladding burst correlations. These results also provide an improved foundation for future BISON model development, FFRD susceptibility calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

The Role of Three-dimensional Effects on Ion Injection and Acceleration in Perpendicular Shocks

Understanding the conditions that enable particle acceleration at nonrelativistic collisionless shocks is essential to unveil the origin of cosmic rays. We employ 2D and 3D hybrid simulations (with kinetic ions and fluid electrons) to explore particle acceleration and magnetic field amplification in nonrelativistic perpendicular shocks, focusing on the role of shock drift acceleration and its dependence on the shock Mach number. We perform an analysis of the ion injection process and demonstrate why efficient acceleration is only observed in 3D. In particular, we show that ion injection critically depends on the “porosity” of the magnetic turbulence in the downstream region near the shock, a property describing how easily the postshock region allows particles to traverse it and return upstream without being trapped. This effect can only be properly captured in 3D. Additionally, we explore the impact of numerical resolution on ion energization, highlighting how resolving small-scale turbulence—on scales below the thermal ion gyroradius—is essential for accurately modeling particle injection. Overall, our results emphasize the necessity of high-resolution 3D simulations to capture the fundamental microphysics driving particle acceleration at perpendicular shocks.

Orusa, Luca [Princeton Univ., NJ (United States); ↗

DOE OSTI2026

Distortion of a Relativistic Jet Echoing a Magnetic Flux Eruption

Magnetized accretion onto spinning black holes (BHs) can accumulate a large magnetic flux across the event horizon and launch a pair of relativistic jets via the Blandford–Znajek mechanism. In the magnetically saturated (arrested) state, excess magnetic flux is ejected from the BH in episodic magnetic flux eruptions, which result in a significant yet temporary reduction of jet power. We analyze results of a high-resolution 3D general-relativistic magnetohydrodynamic numerical simulation of geometrically thick magnetically saturated accretion onto a high-spin Kerr BH for a single cycle of magnetic flux eruption and accumulation. We show that following an eruption, a weakened jet develops a strong helical distortion with a distinct structure of magnetic fields—the poloidal field along the jet core is unaffected by the eruption; while toroidal field lines, ejected from the BH during the eruption and later readvected onto it, form poloidal “bypasses” along the inner jet sheath. Such a distortion may appear in sources fed by geometrically thick accretion flows as an asymmetric superluminal knot, strongly interacting with the jet sheath along an oblique working surface. The jet section repowered by magnetic flux reaccumulated on the BH is tilted by a few degrees, implying significant variations in radiation boost towards observers of BL Lac blazars. The intrinsic structure of the jet spine is consistent with axisymmetric semi-analytical models.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRSDate not supplied

LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Secondary Flows↗

NASA NTRSDate not supplied

A GPU-based Approach for Turbomachinery Application

Gas Turbine↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗