SEARCH · Search NASA

Results for “simulation framework”

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 847 records · Page 47

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Multi-scale simulations of wind farm performance with complex terrain and weather events

Predicting the response of wind farms to local weather is necessary for optimal design and operation. In this work, simulation and analysis of a frontal passage through a utility scale wind farm is achieved for the first time using a seamless multi-scale modeling approach. A generalized actuator disk (GAD) wind turbine model is used to represent turbine wakes, and results are compared to novel radar observations during the frontal passage. The Weather Research and Forecasting (WRF) model is employed with a nested grid setup that allows for coupling between multi-scale atmospheric conditions and turbine response. Starting with mesoscale forcing, the atmosphere is dynamically downscaled to the region of interest, where the interaction between turbulent flows and individual wind turbines is simulated with the GAD at10 m resolution. Several improvements are made to the GAD model to mimic realistic turbine operation, including a yawing capability and a power output calculation. Ultimately, the model is able to capture both the dynamics of the frontal passage and the turbine response; predictions show good agreement with observed turbine wake structure, velocity deficit, and power output after accounting for a phase shift in the mesoscale forcing. This study demonstrates the utility of the WRF-GAD framework for simulating wind farm performance under complex atmospheric conditions.

17 WIND ENERGY↗

DOE OSTI2025

EMT-TS Hybrid Simulation for Large Power Grids Considering IBR-Driven Dynamics

The escalating integration of inverter-based resources (IBRs) poses new challenges to power systems by introducing fast dynamics with higher frequencies, which may need to be simulated by an electromagnetic transient (EMT) program. As an alternative to conducting EMT simulations for the entire system, which is typically time consuming, hybrid simulation between EMT and phasor-domain transient stability (TS) can greatly reduce the computational burden while preserving the detailed fast dynamics in the EMT zone. This paper establishes an EMT-TS hybrid simulation platform using open-source tools, specifically ParaEMT, GridPACK, and HELICS, which are the EMT simulator, TS simulator, and interface framework, respectively. Case studies on the 240-bus Western Electricity Coordinating Council (WECC) system demonstrate that the developed ParaEMT-HELICS-GridPACK hybrid simulator can accurately capture both slow electromechanical and fast IBR-driven dynamics with a 2.4x speedup.

electromagnetic transient↗

DOE OSTI2025

Analyzing the Impact of Future Weather Data on Energy Consumption in Weatherization Assistant

This study supports the mission of the U.S. Department of Energy’s Weatherization Assistance Program (WAP), which aims to increase the energy efficiency of dwellings and reduce their total residential expenditures. Specifically, we examine how projected future climate conditions may affect residential building energy performance by integrating future weather data into the National Energy Audit Tool (NEAT). Since WAP evaluates the cost-effectiveness of retrofit measures over lifespans of up to 30 years, accounting for evolving climate conditions is increasingly important. To reflect future household energy demands, this study replaces historically based Typical Meteorological Year (TMY3) weather inputs with Future Typical Meteorological Year (fTMY) datasets derived from global climate model (GCM) projections. A simulation-based framework was established to enable NEAT analysis under future weather conditions. This workflow involves converting EPW-format weather files into JSON inputs compatible with NEAT and generating degree-hour metrics needed for load calculations. The fTMY dataset used in this study was developed by Oak Ridge National Laboratory through downscaling of six GCMs under different emission scenarios and covers the period from 2020 to 2100. In contrast, the TMY3 dataset is based on historical weather data from 1961 to 1990. Simulations were conducted for benchmark single-family prototype buildings across ASHRAE climate zones 1–7, which cover all regions of the U.S. except the subarctic Zone 8 in northern Alaska, evaluating both heating and cooling loads under TMY3 and fTMY conditions. Four foundation types were tested, while heating systems were standardized, as NEAT does not differentiate thermal energy load by HVAC system type in its load calculations. Results show that fTMY weather input consistently yield lower heating loads and higher cooling loads across most locations, aligning with expected climate warming trends. Notably, colder regions such as zones 6A, 6B, and 7 experience marked reductions in heating load, while warmer and transitional zones, such as 2A (Lufkin, TX) and 3C (San Francisco, CA), have substantial increases in cooling loads. Although this study does not directly assess the performance of retrofit measures under future climate conditions, it provides a critical foundation for doing so. By quantifying shifts in baseline (i.e., pre-retrofit case) energy loads between historical and future weather files, the study highlights the importance of integrating climate-responsive data into audit tools. These findings will inform future efforts to evaluate the long-term effectiveness and cost-effectiveness of weatherization measures under changing climate conditions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2026

Exciton-Defect Interaction and Optical Properties from a First-Principles T-Matrix Approach

Understanding exciton-defect interactions is critical for optimizing optoelectronic and quantum information applications in many materials. However, ab initio simulations of material properties with defects are often limited to high defect density. Here, we study effects of exciton-defect interactions on optical absorption and photoluminescence spectra in monolayer MoS 2 using a first-principles T-matrix approach. We demonstrate that exciton-defect bound states can be captured by the disorderaveraged Green’s function with the T-matrix approximation and further analyze their optical properties. Our approach yields photoluminescence spectra in good agreement with experiments and provides a new, computationally efficient framework for simulating optical properties of disordered 2D materials from firstprinciples.

T-matrix↗

DOE OSTI2021

A Review of Lattice-Boltzmann Models Coupled with Geochemical Modeling Applied for Simulation of Advanced Waterflooding and Enhanced Oil Recovery Processes

To maintain economic profit and improve the oil production efficiency after the primary and secondary production phase, advanced waterflooding techniques such as low salinity waterflooding in carbonate reservoirs have been investigated in numerical simulations, laboratory experiments, and field pilot tests. Multiple underlying mechanisms have been proposed based on these studies, and they are still under debate. Various numerical modeling approaches are introduced, but there exists a lack of a pore-scale comprehensive modeling scheme to fully understand the processes. Lattice-Boltzmann method (LBM) is a type of numerical fluid flow modeling technique that shows capabilities and flexibilities in modeling pore-scale fluid flow to integrate physical–chemical processes within complex structures. The intrinsic feature of LBM makes it a promising framework for simulating advanced waterflooding due to its flexibility, accuracy, and parallel efficiency. LBM works either by itself for solving reactive transport problems or by coupling with a third-party reaction solver. This review mainly introduces the LBM fluid flow and reactive transport capabilities and the concept and modeling approaches to simulate advanced waterflooding techniques. Meanwhile, an evaluation of the coupled LBM models for enhanced oil recovery (EOR) simulations is discussed with future research challenges and directions concluded.

02 PETROLEUM↗

DOE OSTI2025

Development of interatomic potential and effect of ordering on defect properties in CrMnV

Developing materials that can withstand extreme environments, such as high radiation doses and elevated temperatures, is crucial for next-generation particle accelerators, including the 2.4 MW Long-Baseline Neutrino Facility. High-Entropy Alloys have emerged as promising candidates for beam window materials due to their superior mechanical strength, corrosion resistance, and radiation tolerance. In this study, we focus on the Cr–Mn–V alloy system, developing and employing machine-learning interatomic potentials (MLIPs) to investigate the formation of an ordered phase and its influence on defect properties. Using hybrid Monte Carlo-Molecular Dynamics simulations, we observe the formation of a B2-ordered phase at lower temperatures, consistent with Density Functional Theory (DFT) predictions. Ordered structures display a bimodal distribution of migration energies and reduced mean square displacement values, indicating suppressed vacancy diffusion. Our results also show that the migration energy barrier varies based on the atomic species, with Mn and V exhibiting the highest and lowest average barriers, respectively. These findings suggest that atomic ordering inhibits defect mobility, potentially enhancing the radiation resistance of CrMnV alloys. The validated MLIP provides a reliable framework for simulations that are faster than traditional DFT while maintaining the accuracy required to study defect and ordering properties.

36 MATERIALS SCIENCE↗

NASA NTRS2014

Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

Flexible inhibitors are generally used in solid rocket motors (SRMs) as a means to control the burning of propellant. Vortices generated by the flow of propellant around the flexible inhibitors have been identified as a driving source of instabilities that can lead to thrust oscillations in launch vehicles. Potential coupling between the SRM thrust oscillations and structural vibration modes is an important risk factor in launch vehicle design. As a means to predict and better understand these phenomena, a multidisciplinary simulation capability that couples the NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This capability is crucial to the development of NASA's new space launch system (SLS). This paper summarizes the efforts in applying the coupled software to demonstrate and investigate fluid-structure interaction (FSI) phenomena between pressure waves and flexible inhibitors inside reusable solid rocket motors (RSRMs). The features of the fluid and structural solvers are described in detail, and the coupling methodology and interfacial continuity requirements are then presented in a general Eulerian-Lagrangian framework. The simulations presented herein utilize production level CFD with hybrid RANS/LES turbulence modeling and grid resolution in excess of 80 million cells. The fluid domain in the SRM is discretized using a general mixed polyhedral unstructured mesh, while full 3D shell elements are utilized in the structural domain for the flexible inhibitors. Verifications against analytical solutions for a structural model under a steady uniform pressure condition and under dynamic modal analysis show excellent agreement in terms of displacement distribution and eigenmode frequencies. The preliminary coupled results indicate that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor. This insight could have profound implications for SRM and flexible inhibitor designs for current and future launch vehicles including SLS.

Yang, H. Q.↗

DOE OSTI2024

Using an Isotope Enabled Mass Balance to Evaluate Existing Land Surface Models

Abstract Land surface models (LSMs) play a crucial role in elucidating water and carbon cycles by simulating processes such as plant transpiration and evaporation from bare soil, yet calibration often relies on comparing LSM outputs of landscape total evapotranspiration ( ET ) and discharge with measured bulk fluxes. Discrepancies in partitioning into component fluxes predicted by various LSMs have been noted, prompting the need for improved evaluation methods. Stable water isotopes serve as effective tracers of component hydrologic fluxes, but data and model integration challenges have hindered their widespread application. Leveraging National Ecological Observation Network measurements of water isotope ratios at 16 US sites over 3 years combined with LSM‐modeled fluxes, we employed an isotope‐enabled mass balance framework to simulate ET isotope values ( δET ) within three operational LSMs (Mosaic, Noah, and VIC) to evaluate their partitioning. Models simulating δET values consistent with observations were deemed more reflective of water cycling in these ecosystems. Mosaic exhibited the best overall performance (Kling‐Gupta Efficiency of 0.28). For both Mosaic and Noah there were robust correlations between bare soil evaporation fraction and error (negative) as well as transpiration fraction and error (positive). We found the point at which errors are smallest ( x ‐intercept of the multi‐site regression) is at a higher transpiration fraction than is currently specified in the models. Which means that transpiration fraction is underestimated on average. Stable isotope tracers offer an additional tool for model evaluation and identifying areas for improvement, potentially enhancing LSM simulations and our understanding of land‐surface hydrologic processes.

58 GEOSCIENCES↗

NASA NTRS2011

Simulations of Cyclone Sidr in the Bay of Bengal with a High-Resolution Model: Sensitivity to Large-Scale Boundary Forcing

The predictability of Cyclone Sidr in the Bay of Bengal was explored in terms of track and intensity using the Advanced Research Hurricane Weather Research Forecast (AHW) model. This constitutes the first application of the AHW over an area that lies outside the region of the North Atlantic for which this model was developed and tested. Several experiments were conducted to understand the possible contributing factors that affected Sidr s intensity and track simulation by varying the initial start time and domain size. Results show that Sidr s track was strongly controlled by the synoptic flow at the 500-hPa level, seen especially due to the strong mid-latitude westerly over north-central India. A 96-h forecast produced westerly winds over north-central India at the 500-hPa level that were notably weaker; this likely caused the modeled cyclone track to drift from the observed actual track. Reducing the model domain size reduced model error in the synoptic-scale winds at 500 hPa and produced an improved cyclone track. Specifically, the cyclone track appeared to be sensitive to the upstream synoptic flow, and was, therefore, sensitive to the location of the western boundary of the domain. However, cyclone intensity remained largely unaffected by this synoptic wind error at the 500-hPa level. Comparison of the high resolution, moving nested domain with a single coarser resolution domain showed little difference in tracks, but resulted in significantly different intensities. Experiments on the domain size with regard to the total precipitation simulated by the model showed that precipitation patterns and 10-m surface winds were also different. This was mainly due to the mid-latitude westerly flow across the west side of the model domain. The analysis also suggested that the total precipitation pattern and track was unchanged when the domain was extended toward the east, north, and south. Furthermore, this highlights our conclusion that Sidr was influenced from the west side of the domain. The displacement error was significantly reduced after the domain size from the western model boundary was decreased. Study results demonstrate the capability and need of a high-resolution mesoscale modeling framework for simulating the complex interactions that contribute to the formation of tropical cyclones over the Bay of Bengal region

Kumar, Anil↗

DOE OSTI2025

Computational Modeling of Graphite Degradation due to Molten Salt Infiltration and Wear

Molten-salt reactors (MSRs) represent a promising next-generation reactor design, with graphite serving as a moderator and/or reflector in several designs. However, due to limited experimental data and operational experience, a technical understanding of the structural integrity of graphite in molten salt environments remains incomplete. This report presents a modeling-based evaluation of graphite degradation in MSR environments, focusing on the effects of salt infiltration in fuel salt-based designs and surface wear in pebble bed reactor designs. The objective of this study is to enhance understanding of the structural integrity challenges posed by these degradation mechanisms and to provide a framework for assessing graphite behavior in MSRs. The first part of the report investigates the phenomenon of molten salt infiltration into graphite. This infiltration occurs when molten salt permeates the interconnected pore structure of the graphite moderator, driven by factors such as pressure differentials and the physical properties of both the salt and graphite. The infiltration process is influenced by characteristics of the pore structure, viscosity of the molten salt, and the interfacial energies between the graphite, salt, and the atmosphere within the graphite pore. Utilizing a coupled multiphysics modeling approach with Grizzly software, the study evaluates the stress induced by internal heat sources due to infiltration, which can lead to structural concerns. This evaluation is crucial for understanding how infiltration affects the mechanical integrity of graphite components in MSRs. The study considers the Molten-Salt Reactor Experiment (MSRE) graphite stringer geometry due to the availability of relevant data. Through detailed finite element analysis, the study examines stress distributions at varying infiltration percentages, revealing that stress levels increase with higher amounts of infiltration. Rare-event simulations, using the parallel subset simulation (PSS) framework, further quantify the failure probabilities under input uncertainties, with a user-specified failure metric. The PSS framework also identifies critical input parameters that significantly affect the stress values, including infiltration amount, thermal conductivity, and power density. Additionally, considering realistic reactor scenarios, the analysis was performed to account for the combined effects of radiation and infiltration, and modeling strategies on how to analyze new reactor designs or new graphite grades are discussed. The second part of the report focuses on wear mechanisms in pebble bed-based MSRs. As graphite fuel pebbles interact with the graphite reflector block, wear can result in material loss and the formation of surface defects, which may act as stress concentrators. A similar multiphysics modeling framework is employed to assess the impact of wear on the structural integrity of graphite components. This study considers a generic fluoride-cooled high-temperature reactor (gFHR) design due to the availability of comprehensive data. Worst-case scenario dimensions of the reflector blocks were analyzed under thermal and radiation conditions. Subsequently, wear in the form of idealized pits and grooves is modeled on the inner surface of the graphite block, with the maximum stress from previous simulations. The simulations show that groove-type defects are more detrimental than pits, leading to higher stress concentrations. Considering worst-case simulation scenarios and experimental wear rates, it was determined that the formation of a surface defect critical enough to affect the stress may not be possible in a gFHR design. Overall, the findings of this research contribute to the development of robust modeling tools for predicting graphite behavior under various operational conditions in MSRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Towards Z-Next: The Integration of Theory, Experiments, and Computational Simulation in a Bayesian Data Assimilation Framework

Making reliable predictions in the presence of uncertainty is critical to high-consequence modeling and simulation activities, such as those encountered at Sandia National Laboratories. Surrogate or reduced-order models are often used to mitigate the expense of performing quality uncertainty analyses with high-fidelity, physics-based codes. However, phenomenological surrogate models do not always adhere to important physics and system properties. This project develops surrogate models that integrate physical theory with experimental data through a maximally-informative framework that accounts for the many uncertainties present in computational modeling problems. Correlations between relevant outputs are preserved through the use of multi-output or co-predictive surrogate models; known physical properties (specifically monotoncity) are also preserved; and unknown physics and phenomena are detected using a causal analysis. By endowing surrogate models with key properties of the physical system being studied, their predictive power is arguably enhanced, allowing for reliable simulations and analyses at a reduced computational cost.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Dark Energy Survey Year 3 results: galaxy–halo connection from galaxy–galaxy lensing

ABSTRACT Galaxy–galaxy lensing is a powerful probe of the connection between galaxies and their host dark matter haloes, which is important both for galaxy evolution and cosmology. We extend the measurement and modelling of the galaxy–galaxy lensing signal in the recent Dark Energy Survey Year 3 cosmology analysis to the highly non-linear scales (∼100 kpc). This extension enables us to study the galaxy–halo connection via a Halo Occupation Distribution (HOD) framework for the two lens samples used in the cosmology analysis: a luminous red galaxy sample (redmagic) and a magnitude-limited galaxy sample (maglim). We find that redmagic (maglim) galaxies typically live in dark matter haloes of mass log10(Mh/M⊙) ≈ 13.7 which is roughly constant over redshift (13.3−13.5 depending on redshift). We constrain these masses to ${\sim}15{{\ \rm per\ cent}}$, approximately 1.5 times improvement over the previous work. We also constrain the linear galaxy bias more than five times better than what is inferred by the cosmological scales only. We find the satellite fraction for redmagic (maglim) to be ∼0.1−0.2 (0.1−0.3) with no clear trend in redshift. Our constraints on these halo properties are broadly consistent with other available estimates from previous work, large-scale constraints, and simulations. The framework built in this paper will be used for future HOD studies with other galaxy samples and extensions for cosmological analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Parallel derivative-free optimization for simulation-based design of behind-the-meter energy systems

In this work, the integrated design and dispatch of behind-the-meter or distributed resources (e.g. stationary battery storage and solar PV generation) is considered. A simulation-based framework is employed, generating high-fidelity results with closed-loop predictive control at a fine resolution, at the expense of high computational cost (several minutes to a few hours per design point). To address this challenge, parallel derivative-free design methods are considered. Four methods are compared, including state-of-the-art surrogate-based methods (Radial-Basis Functions and Gaussian processes) and sampling strategies, an evolutionary-based method, and a simple sequential grid refinement method. As a case study, two types of design problem with increasing complexity are considered, namely, the design of behind-the-meter resources (three design variables) and the inclusion of grid capacity (four design variables). The second yields a constrained design problem for which violations can only be determined after solving the computationally expensive simulation. For the three-dimensional case, all methods present a good performance, achieving a solution within 1% of the optimum after the first iteration, with the sequential grid refinement exhibiting the fastest convergence and achieving the best final objective value. This indicates that the parallel evaluation of multiple sampling points may be more important than the choice of method for small decision spaces. For the four-dimensional constrained case, the Genetic Algorithm presents the best tradeoff between performance and computational effort, while the rough objective function terrain generated by constraint violation penalties reduces the performance of surrogate-based methods. Contour plots with flat regions indicate flexibility in the optimal design and highlight the importance of characterizing the solution space.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2018

Burial-driven methane recycling in marine gas hydrate systems

Natural gas hydrate may be buried with sediments until it is no longer stable at a given pressure and temperature, resulting in conversion of hydrate into free gas. This gas may migrate upward and recycle back into the hydrate stability zone to form hydrate. As of yet, however, no quantitative description of the methane recycling process has been developed using multiphase flow simulations to model burial-driven gas hydrate recycling. In this study, we present a series of 1D multiphase transport simulations to investigate the methane recycling process in detail. By invoking the effects of capillary phenomena on hydrate and gas formation in pores of varying size, we find that a free gas phase can migrate a significant distance above the bulk base of hydrate stability. Since the top of the free gas occurrence is often identified as the base of the hydrate stability zone from seismic data, our results demonstrate that not only could this assumption mischaracterize a hydrate system, but that under recycling conditions the highest hydrate saturations can occur beneath the top of the free gas occurrence. We show that the presence of pore size distributions requires a replacement zone through which hydrate saturations progressively decrease with depth and are replaced with free gas. This replacement zone works to buffer against significant gas buildup that could lead to fracturing of overlying sediments. In conclusion, this work provides a framework for simulating flow and transport of methane within the 3-phase stability zone from a mass conservation perspective.

58 GEOSCIENCES↗

DOE OSTI2020

2D Retrieval Frameworks for Hot Jupiter Phase Curves

Spectroscopic phase curves provide unique access to the three-dimensional properties of transiting exoplanet atmospheres. However, a modeling framework must be developed to deliver accurate inferences of atmospheric properties for these complex data sets. Here, we develop an approach to retrieve temperature structures and molecular abundances from phase curve spectra at any orbital phase. In the context of a representative hot Jupiter with a large day–night temperature contrast, we examine the biases in typical one-dimensional (1D) retrievals as a function of orbital phase/geometry, compared to two-dimensional (2D) models that appropriately capture the disk-integrated phase geometry. We guide our intuition by applying our new framework on a simulated Hubble Space Telescope (HST)+Spitzer phase curve data set in which the “truth” is known, followed by an application to the spectroscopic phase curve of the canonical hot Jupiter, WASP-43b. We also demonstrate the retrieval framework on simulated James Webb Space Telescope (JWST) phase curve observations. We apply our new geometric framework to a joint fit of all spectroscopic phases, assuming longitudinal molecular abundance homogeneity, resulting in an a factor of 2 improvement in abundances precision when compared to individual phase constraints. With a 1D retrieval model on simulated HST+Spitzer data, we find strongly biased molecular abundances for CH{sub 4} and CO{sub 2} at most orbital phases. With 2D, the day and night profiles retrieved from WASP-43b remain consistent throughout the orbit. JWST retrievals show that a 2D model is strongly favored at all orbital phases. Based on our new 2D retrieval implementation, we provide recommendations on when 1D models are appropriate and when more complex phase geometries involving multiple TP profiles are required to obtain an unbiased view of tidally locked planetary atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Large-eddy simulations for hill terrains: validation with wind-tunnel and field measurements

Development and validation of large-eddy simulation (LES) framework to study atmospheric boundary layer flows over complex terrains is reported here. The LES model uses the fourth-order time-accurate Runge–Kutta scheme and a fractional step method. The inflow boundary conditions are generated by using the so-called recycling (or mapping) technique. Evaluation of potential in-land wind park locations is the main application of the validated model. In this work, LES are carried out for turbulent boundary-layer flows over both simple and complex hill geometries. The prediction of the flow separation and reattachment length for a steeper wind-tunnel hill was closer to the measurements than previous numerical studies reported by other authors for the same hill geometry. For the complex hill case, the LES model showed good agreement with full-scale measurements.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

The Monolithic Heat Pipe Microreactor Reference Plant Model

This work introduces a reference plant model for a generic monolithic heat-pipe-cooled microreactor. The model will serve as a springboard to develop future evaluation models in the licensing process of similar microreactor designs at the U.S. Nuclear Regulatory Commission. This model has been developed with the Comprehensive Reactor Analysis Bundle and its specifications are based on open literature publications for the eVinci TM design. BlueCRAB is the U.S. Nu- clear Regulatory Commission non-light-water reactor analysis system based on the Multiphysics Object-Oriented Simulation Environment framework, which can couple the Griffin, BISON, and Sockeye applications to resolve the various physics that are essential for the safety analysis of this type of reactor system. The core specifications includes tristructural isotropic fuel, graphite monolith, graphite reflectors, and drums composed of graphite and B 4 C. No moderator or burnable poison pins are used in the design. The fuel enrichment is reduced to control excess reactivity in the core. This core design is not optimized and only serves for testing purposes, since the primary objective of this work is to exercise the multiphysics coupling for this type of reactor system. A three dimensional (3D) core heterogeneous Griffin discrete ordinates (SN) transport model allows the precise calculation of the flux distribution and pin powers. Griffin transfers the power density distribution and obtains a temperature distribution to and from BISON. The BISON model com- putes the 3D core temperature distribution and is coupled to 876 Sockeye subapplications running a heat pipe model. This 3D conduction model is coupled to the various heat pipes via heat flux boundary conditions. The model includes a small gap between the heat pipe and the monolith. Convective heat transfer boundaries with either ambient temperature or condenser temperature as heat sinks are imposed at the model boundaries. The 2D Sockeye heat pipe model uses a vapor- only methodology, which provides the needed resolution for transient calculations and allows the determination of various heat pipe limits. This approach is superior to the superconductor model traditionally used in steady-state calculations. BlueCRAB computes steady-state power and temperature distributions that serve as the initial condition for a loss-of-heat-sink transient simulation. The steady-state results show significant peaking due to the position of the control drum, but this is a characteristic of the particular design used, which is not optimized at this stage. The transient results show the reactor power slowly stabilizing towards a 3% power level after the partial loss of secondary heat removal. Several recriticalities are observed due to cooling through the secondary system but the reactor is self-stabilizing and behaves as expected.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Toward engineering lattice structures with the material point method (MPM)

This study examines the potential of two variants of the material point method—the generalized interpolation material point (GIMP) and dual domain material point (DDMP) methods—in developing a robust computational framework for engineering lattice structures under different loading conditions. The study begins with assessing the ability of the two methods in predicting elastic buckling phenomena using column geometries with and without initial geometric imperfections. The results indicate that both methods effectively capture buckling phenomena when initial geometric imperfections are introduced. After this verification step, we create several models of tetrahedral lattice structures with varying strut diameter and orientation and subject them to quasi-static loading. We then validate the numerical results using laboratory test results. The results show that, while both methods accurately predict load–displacement curves in the pre-buckling regime, their predictive capabilities diminish in the post-buckling regime. Through visual comparison between the numerical and experimental deformed shapes, it appears that the discrepancies between model and experimental results are attributed to initial geometric imperfections in the lattices that occurred during 3D printing. We then establish a second set of lattice models where different types of initial geometric imperfections are considered. The results from these models show that imperfections have a negligible influence in the pre-buckling regime but affect the behavior considerably in the post-buckling regime. As a final step in this work, we subject the lattice models to impact loading and employ hypothetical soft and stiff materials. These results show that the lattice stiffness, which depends on material stiffness, strut diameter, and orientation, significantly influences the ability of a lattice structure to resist impact. In particular, we find that a stiffer lattice (i.e., one made with a stiff material and thicker struts) is capable of absorbing more energy than a softer one during impact. Although material nonlinearities, inelasticity, and detailed contact formulations are not considered in this study, the findings obtained herein lay the groundwork for engineering lattice structures under extreme loading conditions through a simulation-driven framework based on particle-based methods.

97 MATHEMATICS AND COMPUTING↗