Search NASA⌕ Search

SEARCH · Search NASA

Results for “computer simulation model development”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Cluster Dynamics Modeling Needs for the Advanced Materials and Manufacturing Technologies Program

This milestone report aims to identify and assess the cluster dynamics (CD) modeling requirements within the Department of Energy's Office of Nuclear Energy (DOE-NE) Advanced Materials and Manufacturing Technologies (AMMT) program and to communicate these needs to the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The goal is to ensure NEAMS is well-informed about the CD modeling requirements to support AMMT's mission of accelerating the development, qualification, demonstration, and deployment of advanced structural materials and manufacturing for nuclear energy applications. CD modeling is an essential tool for predicting the degradation of structural materials under irradiation, which is a key component of AMMT's accelerated qualification process. The AMMT program focuses on both additively manufactured and wrought structural alloys, such as laser powder-bed fusion 316H austenitic stainless steel, alloy 709, Haynes 244, and alloy 617. These materials require a generalized CD modeling framework to facilitate rapid model development and computational simulation. A flexible, generalized CD software, similar to the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework, would enable modeling of various cluster types, including defect clusters, defect-solute clusters, and multicomponent clusters, incorporating thermodynamics and kinetics parameters. Radiation effects, microstructural feature evolution, and multi-dimensional modeling are critical considerations for the CD model. The usability of the CD code should allow for easy modification and coupling with MOOSE-based simulations. Additionally, the software should adhere to Nuclear Quality Assurance-1 standards, include a testing suite for verification and validation, and be version-controlled within a national laboratory-managed Git repository. Benchmark problems are needed to assess code predictions and performance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Surrogate modeling of Cellular-Potts agent-based models as a segmentation task using the U-Net neural network architecture

The Cellular-Potts model is a powerful and ubiquitous framework for developing computational models for simulating complex multicellular biological systems. Cellular-Potts models (CPMs) are often computationally expensive due to the explicit modeling of interactions among large numbers of individual model agents and diffusive fields described by partial differential equations (PDEs). In this work, we develop a convolutional neural network (CNN) surrogate model using a U-Net architecture that accounts for periodic boundary conditions. We use this model to accelerate the evaluation of a mechanistic CPM previously used to investigate in vitro vasculogenesis. The surrogate model was trained to predict 100 computational steps ahead (Monte-Carlo steps, MCS), accelerating simulation evaluations by a factor of 562 times compared to single-core CPM code execution on CPU. Over short timescales of up to 3 recursive evaluations, or 300 MCS, our model captures the emergent behaviors demonstrated by the original Cellular-Potts model such as vessel sprouting, extension and anastomosis, and contraction of vascular lacunae. This approach demonstrates the potential for deep learning to serve as a step toward efficient surrogate models for CPM simulations, enabling faster evaluation of computationally expensive CPM simulations of biological processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Thermal Management for Planar Package Power Electronics (CRADA Final Report)

The National Renewable Energy Laboratory (NREL) and John Deere Electronic Solutions (JDES) collaborated to develop and evaluate computer models and simulations related to the thermal performance of semiconductor device packaging for an inverter of an off-road vehicle ("Semiconductor Packaging"). The objective of the research was for NREL and JDES to develop a two or three-dimensional computer-aided design model of the Semiconductor Packaging ("Computer Model") for thermal performance evaluation in a simulation. The project developed computer models of a Semiconductor Package of silicon-carbide semiconductor devices with appropriate thermal dissipation via one or more of the following thermal features: a double-sided planar cooling configuration, an air-cooled configuration, a liquid-cooled configuration for off-road inverter applications in a relevant simulated operating environment. It is noted that surface area and prototype product embodiment consisting of air-cooled and liquid-cooled configurations may and may not be in direct contact with power semiconductor chips. The idea was to explore and develop packaging and thermal management technology for power semiconductors that was most effective in the performance yet had least burden in overall product cost for a given application.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Fault-Tolerant Decentralized Control for Large-Scale Inverter-Based Resources for Active Power Tracking

Integration of inverter-based resources (IBRs) which lack the intrinsic characteristics such as the inertial response of the traditional synchronous-generator (SG)-based sources presents a new challenge in the form of analyzing the grid stability under their presence. While the dynamic composition of IBRs differs from that of the SGs, the control objective remains similar in terms of tracking the desired active power. This letter presents a decentralized primal-dual-based fault-tolerant control framework for the power allocation in IBRs. Overall, a hierarchical control algorithm is developed with a lower level addressing the current control and the parameter estimation for the IBRs and the higher level acting as the reference power generator to the low level based on the desired active power profile. The decentralized network-based algorithm adaptively splits the desired power between the IBRs taking into consideration the health of the IBRs transmission lines. The proposed framework is tested through a simulation on the network of IBRs and the high-level controller performance is compared against the existing framework in the literature. The proposed algorithm shows significant performance improvement in the magnitude of power deviation and settling time to the nominal value under faulty conditions as compared to the algorithm in the literature.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Universal Solvent Viscosity Reduction via Hydrogen Bonding Disruptors

Liquid Ion Solutions LLC (DBA RoCo Global) in partnership with Carnegie Mellon University and Carbon Capture Scientific LLC, has performed lab-scale development and evaluation of novel additives that lower the viscosity of water-lean amine solvents for post-combustion carbon dioxide capture. This project focuses on developing additives that minimize the formation of long-range electrostatic and hydrogen bonding (HB) networks, decreasing the solvent viscosity, improving diffusion, and improving the process economics. The project objectives included: 1) performing computer simulation to understand the molecular interactions of the additive molecules in water-lean CO 2 capture solvents, 2) design and synthesis of HB disruptors additives, 3) performance testing with additive molecules on model amine solvents, and 4) demonstration of the effectiveness of the optimized additives in the presence of synthetic flue gas. To meet the abovementioned objectives, the project team utilized a holistic approach that combines molecular simulation, experimental testing, and economic analysis studies. The project team developed ab initio molecular model and then perform computer simulation to develop relationship between hydrogen bonding, viscosity, and performed quantitative analysis of additive on the viscosity of the solvent. The team completed computational comparative study on a range of organic functional groups such as ethers, esters, cyclic carbonates, alkanes, and ammonium salts for their effect on viscosity gaining key insights into molecular interactions and the impact of various functional groups and molecular shapes on viscosity. Assisted with molecular simulation insights, the project team conducted additive synthesis and testing, including a proof-of-concept study, additive screening, optimization, and synthetic flue gas testing. The experimental proof-of-concept study proved that the hydrogen bonding acceptors result in significant decrease of viscosities. Detailed additive screening (exploring various functionalities and molecular structures) has been performed. Several promising additives showed excellent reduction in viscosity (30-41%) at 5% additive loading, and over 50% viscosity reduction at 10% additive loading for the model solvents. The team also performed complex screening studies on additive loadings and mixing effect among additives using the design of experiments. Based on multiple screening experiments, one additive-solvent candidate was down-selected for synthetic flue gas testing. A 100-hour continuous absorption/desorption study was conducted under simulated flue gas using a lab-scale continuous capture and separation system. No degradation (for both solvent and additive) was observed based on the GC results of the solvent samples collected from the continuous study. The team conducted preliminary engineering analyses and cost-benefit analyses to quantify the potential economic benefits of the additive approach for solvent viscosity reduction. Based on the experimental data, CO 2 capture cost savings from the capital and operating cost savings are estimated at $\$$4.7/tonne and $\$$0.3/tonne CO 2 captured, respectively. Considering the additive cost, the net benefit is estimated to be between $\$$4.32~$\$$4.86/tonne CO 2 captured.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

A Transported Livengood–Wu Integral Model for Knock Prediction in Computational Fluid Dynamics Simulation

This work describes the development of a transported Livengood–Wu (L–W) integral model for computational fluid dynamics (CFD) simulation to predict autoignition and engine knock tendency. The currently employed L–W integral model considers both single-stage and two-stage ignition processes, thus can be generally applied to different fuels such as paraffin, olefin, aromatics, and alcohol. The model implementation is first validated in simulations of homogeneous charge compression ignition (HCCI) combustion for three different fuels, showing good accuracy in prediction of autoignition timing for fuels with either single-stage or two-stage ignition characteristics. Then, the L–W integral model is coupled with G-equation model to indicate end-gas autoignition and knock tendency in CFD simulations of a direct-injection spark-ignition engine. This modeling approach is about 10 times more efficient than the ones that based on detailed chemistry calculation and pressure oscillation analysis. Two fuels with same Research Octane Number (RON) but different octane sensitivity are studied, namely, Co-Optima alkylate and Co-Optima E30. Feed-forward neural network model in conjunction with multivariable minimization technique is used to generate fuel surrogates with targets of matched RON, octane sensitivity, and ethanol content. The CFD model is validated against experimental data in terms of pressure traces and heat release rate for both fuels under a wide range of operating conditions. The knock tendency—indicated by the fuel energy contained in the autoignited region—of the two fuels at different load conditions correlates well with the experimental results and the fuel octane sensitivity, implying the current knock modeling approach can capture the octane sensitivity effect and can be applied to further investigation on composition of octane sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Real-time Nonlinear Model Predictive Control (NMPC) Strategies using Physics-Based Models for Advanced Lithium-ion Battery Management System (BMS)

Optimal operation of lithium-ion batteries requires robust battery models for advanced battery management systems (ABMS). A nonlinear model predictive control strategy is proposed that directly employs the pseudo-two-dimensional (P2D) model for making predictions. Using robust and efficient model simulation algorithms developed previously, the computational time of the nonlinear model predictive control algorithm is quantified, and the ability to use such models for nonlinear model predictive control for ABMS is established.

25 ENERGY STORAGE↗
DOE OSTI2023

Computational model development and validation of fuel dispersal phenomena

As the US nuclear industry is renewing efforts to extend the rod average burnup limits in existing pressurized water reactors, it becomes imperative to investigate the safety implications associated with high burnup fuel fragmentation and the potential dispersal of fragmented fuel into the reactor coolant system. Additionally, as there is a growing interest in transitioning to new nuclear fuel designs to accommodate higher fuel burnup levels, it becomes essential to conduct studies on fuel behavior during design basis accident scenarios, particularly during a postulated loss-of-coolant accident, in order to ensure the safe operation of light-water reactors. The present study develops a computational model to simulate the complex three-phase flow of the fuel dispersal phenomena that occur following a breach in the cladding, when the stream of solid fuel particles and carrier fission gas stream may interact with the surrounding liquid medium. Furthermore, the proposed framework considers the solid phase as discrete Lagrangian particles while the interpenetrating gas–liquid continuum is modeled using a Eulerian framework. The developed simulation framework is validated with experimental results performed in a separate-effect test facility by comparing predicted particle settlement locations with the experimental observations reported in the literature. Key highlights of this study include the interpretation of the high-pressure boundary conditions and the three-phase flow coupling strategy. Also discussed are the gas–liquid hydrodynamics and the behavior of particles as they are transported through the gas and liquid phases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A linked-scale coupled model of mass erosion and redistribution in plasma-exposed micro-foam surfaces

Surface evolution due to exposure to harsh environments is of importance in many scientific and technological applications. In plasma-exposed materials, the surface receives charged particles from the plasma, leading to a series of processes that drive the system far from equilibrium and may lead to the severe deterioration of the surface properties. Although surface morphological changes are driven by atomic collisions taking place over picoseconds and nanometers, these processes result in mass loss and redistribution of matter over much larger length and time scales. This necessitates a multi-scale approach capable of capturing the range of processes linking primary atomic collision events with engineering-level surface geometry changes. In this paper, we develop a computational model to simulate the morphological evolution and effective erosion rate of micro-architected tungsten foams during low-energy plasma ion bombardment. Furthermore, the model acts on several length scales, with the energy and angular dependence of the sputtering yield of flat tungsten surfaces determined using the SRIM code based on the binary collision approximation. This information is introduced into a low-fluence, short-term Monte Carlo raytracing model, and further into a high-fluence, long-term particle transport model. In the latter, material particles representing billions of atoms are described in a digitized 3-D representation of the foam structure from X-ray tomography data, and are sputtered off and redeposited using the atomistic information. We show that the redeposition of sputtered atoms leads to partial self-healing in the bottom layers with a sharp reduction in the sputtering coefficient of low density foams, roughly 25% of the solid W value. This is in qualitative agreement with recent experiments on low porosity W structures. At high fluence, the foam structure degrades considerably as there are fewer ligaments available to recapture sputtered atoms and, consequently, the sputtering rate increases again.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multilayered Ordered Protein Arrays Self-Assembled from a Mixed Population of Virus-like Particles

Biology shows many examples of spatially controlled assembly of cells and biomacromolecules into hierarchically organized structures, to which many of the complex biological functions are attributed. While such biological structures have inspired the design of synthetic materials, it is still a great challenge to control the spatial arrangement of individual building blocks when assembling multiple types of components into bulk materials. Here, in this study, we report self-assembly of multilayered, ordered protein arrays from mixed populations of virus-like particles (VLPs). We systematically tuned the magnitude of the surface charge of the VLPs via mutagenesis to prepare four different types of VLPs for mixing. A mixture of up to four types of VLPs selectively assembled into higher-order structures in the presence of oppositely charged dendrimers during a gradual lowering of the ionic strength of the solution. The assembly resulted in the formation of three-dimensional ordered VLP arrays with up to four distinct layers including a central core, with each layer comprising a single type of VLP. A coarse-grained computational model was developed and simulated using molecular dynamics to probe the formation of the multilayered, core-shell structure. Our findings establish a simple and versatile bottom-up strategy to synthesize multilayered, ordered materials by controlling the spatial arrangement of multiple types of nanoscale building blocks in a one-pot fabrication.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

LLNL Macroscopic Anisotropic Explosives Research at INL National Security Test Range - Test Results

A select team of 23 engineers, scientists, and explosives specialists from LLNL, LANL, INL, and Marine Raiders from Marine Special Operations Command (MARSOC) and U.S. Special Operations Command (SOCOM) assembled during the second week of November at the INL National Security Test Range near Idaho Falls to investigate and demonstrate fundamental principles of explosives anisotropy. Today's explosives are isotropic in their detonation performance. That is, no matter what direction a detonation runs through bulk explosive, the performance is the same; whereas, anisotropic explosives exhibit different performance, depending on which direction the detonation wave moves through the explosive. The ANISO Team worked in subfreezing temperatures on the Snake River Plain, carrying out 55 experimental explosives shots in four days that lead to a clear understanding of the performance and behavior of an assembly of small, linerless, C4 shaped charges. These shots clearly demonstrated, for the first time, on a macroscopic scale, the principle of anisotropy in measured progression of the detonation through the explosive assembly. The outputs of nine piezo timing pins in the explosive assembly clearly showed detonation progressing through the assembly faster than nominal detonation velocity and moving slower than nominal detonation velocity in the opposite direction. Basic data from these experiments will be used to design and construct explosives assemblies that will be shot in the LLNL High Explosives Applications Facility's (HEAF). These experimental tests will provide refined basic data that will then be used by modelers to develop high explosives models. Computer simulations using these models will then be run to predict performance and design inhomogeneous, anisotropic bulk explosive charges that will be tested at LLNL.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

LDRD poster - Identifying the speciation of salt-based actinides in the presence of contaminants

Molten salt reactors (MSRs) are among the most promising Generation IV nuclear reactors. The salts used in these reactors exhibit a high affinity for moisture and oxygen, and heating does not eliminate the oxides that are formed. Additionally, due to the high operating temperatures, reactor parts experience pronounced corrosion, and the resulting corrosion products integrate into the fuel salt system. It is hypothesized that the speciation of f-elements in molten salts would change in the presence of oxides and corrosion products. These changes in local structure can alter the rheological properties, which are critical to the operational and safety parameters of the reactors. Therefore, an in-depth understanding of how corrosion products influence lanthanide/actinide speciation is crucial. To test this hypothesis, lanthanides (as surrogates for actinides) and actinide chlorides were mixed with oxides (such as nickel oxide and lithium oxide) and/or corrosion products (e.g., nickel chloride) and reacted in molten salts. The resulting phases were analyzed using techniques such as powder X-ray diffraction, Raman spectroscopy, UV-Visible spectroscopy, and nuclear magnetic resonance spectroscopy. In the presence of oxides, due to the high oxophilicity of f-elements, these elements capture oxygen from the oxides, forming lanthanide/actinide oxide or oxychloride phases. Specifically, with cerium chloride, the introduction of lithium oxide to the formation of various final products depending on the oxide concentration. At lower oxide concentrations, cerium oxychloride (CeOCl) in the +3 oxidation state was observed, while an increase in the oxide concentration led to the formation of cerium dioxide (CeO2), where cerium is in the +4 oxidation state. In contrast, with nickel oxide, CeOCl formation was consistently observed, regardless of the oxide concentration. This project successfully determined that the speciation of f-elements change in the presence of oxides and different oxides in molten salts result in different final products. When the influence of nickel chloride was examined, changes in f-element speciation were not observed. However, it reacted with lithium chloride to generate a lithium-nickel-chloride (Li6NiCl8) phase, which was also observed when nickel oxide was reacted. Computational models were developed to simulate the phases anticipated under experimental conditions.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2021

A Poisson equation method for prescribing fully developed non-Newtonian inlet conditions for computational fluid dynamics simulations in models of arbitrary cross-section

Prescribing inlet boundary conditions for computational fluid dynamics (CFD) simulations of internal flow in complex geometries such as anatomical vascular models is challenging. In the absence of patient-specific inlet velocity data, a common approach for long blood vessels is to assume that the inlet flow is fully developed. In vessels of irregular cross section, however, prescribing fully developed conditions is complicated due to the lack of a general closed-form analytical solution. In this study, we develop a simple Poisson equation method for prescribing fully developed inlet conditions for the flow of either Newtonian or non-Newtonian fluids in CFD models of arbitrary cross-section. We first derive the generalized Poisson equation for fully developed flow of a non-Newtonian fluid and we then develop and verify a methodology for numerically computing the solution on any planar boundary domain. In addition, we develop a simple extension of the method for prescribing a non-orthogonal inlet velocity that represents fully developed flow from an upstream tube that is connected to the CFD inlet at a non-orthogonal angle. This may be used to investigate a common source of uncertainty in CFD simulations of internal flow that is due to a lack of information concerning the exact streamwise flow direction at the inlets. Comparison to several Newtonian and non-Newtonian benchmark verification solutions shows the method to be extremely accurate. As a practical demonstration case, we use the method to prescribe fully developed conditions on multiple non-circular inlets for the non-Newtonian flow of blood in a patient-specific model of the inferior vena cava (IVC). Finally, we further demonstrate the utility of the method by performing a sensitivity study using the patient-specific IVC model, wherein we investigate the influence of inlet velocity flow direction on the non-Newtonian IVC hemodynamics. Given its simplicity and computational efficiency, the method is shown to be far superior to alternative approaches for prescribing fully developed inlet conditions in such complicated geometries. In conclusion, to facilitate the adoption of our Poisson equation method, we have distributed our OpenFOAM source code and the associated test cases from this study as open-source software.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Robust Combined Heat and Hybrid Power (CHHP) for High Electrical Efficiency Cogeneration

Georgia Tech (Prime Recipient), the University of Texas at El Paso (UTEP) and the National Energy Technology Laboratory (NETL) investigated a hybrid fuel cell/ gas turbine system concept as a combined heat and hybrid power (CHHP) system for both robust and high power-to-process heat ratio cogeneration. The novelty of the proposed system entailed the distinct, elevated electrical efficiencies it maintains while simultaneously supporting a broad span of heating needs (e.g., supply temperatures) demanded across variable heat loads. The scope included: 1) leveraging a pioneering national lab facility configured for dynamic system operability development of hybrid fuel cell/gas turbine cycles; 2) enabling technology development to adjust and modulate the quality and quantity of thermal supply to bottoming heat loads via novel extreme temperature gas bypass valves. Hybrid fuel cell/gas turbine systems have primarily been reduced-to-practice in a constrained (e.g., initial proof-of-concept) manner and have still demonstrated considerable electrical efficiencies. However, these pre-pilot systems have focused upon electrical efficiencies and electrical power generation as the exclusive energy demand. Such hybrid systems had not been extensively researched or developed for flexible and variable operation consisting of both power and heat demands; however, these variable combined power and heat demands are characteristic of many types of manufacturers such as animal/poultry processing, bakeries and milk/flour/pastry manufacturing, textile mills, and electrochemical processing. Commercially, developing the system into a working combined heat and power system benefits these types of manufacturers by allowing them to meet their power and heat demands at a lower cost, higher efficiency, and/or through onsite generation. Therefore, the technical scope of this project was largely to study and facilitate these hybrid systems as combined heat and hybrid power (CHHP) systems that include dynamic operability for variable heat and power loads and/or grid dynamics for various types of manufacturers. Simulation results were used to predict the performance of the CHHP system and conceptually develop it to achieve desired dynamic operability. Experimentally, the primary goal was to design, manufacture, and experiment upon a high-temperature bypass valve. Experimental data included air mass flow rates through the valve orifice when the valve was changed to variable extent between fully closed and fully open. The experimental data was then used to create a semi-empirical computational model of the bypass valve. Concluded simulation goals for the research included developing computational heat exchanger models for the hybrid system inclusive of the bottoming heat exchanger and the recuperative heat exchanger, and then combining the computational recuperator model with the computational valve model. Afterwards, the computational models were then integrated to predict the dynamic operation of hybrid fuel cell/gas turbine cycles throughout a design space and reporting such. The scope stated in the preceding paragraph was packaged into five specific goals: 1) enabling the simulation of dynamic combined heat and power through the creation of computational, modular heat exchanger models; 2) simulation and exploration of the CHHP system’s performance by integrating the heat exchanger models with the national lab’s pre-existing hybrid system (computational) simulation, but without the recuperator bypass valve concept in order to initially determine how the (baseline) system behaves and can be controlled in order to meet variable heat and power demands; 3) development and initial deployment of the high-temperature recuperator bypass valve technology in order to confirm and characterize the approach; 4) usage of the experimental data for the valve to create a semi-empirical computational model for the bypass valve which could then be combined with the heat exchanger computational models; 5) repeat of the second task of simulating and exploring the system’s performance, but this time including the bypass valve to resolve its efficacy. Tasks were successfully completed, and the general notion of flexibly operating, high electrical efficiency CHHP was further corroborated. Supportive details are provided in the report.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Data-model files associated with the manuscript "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS)

This data package contains the model inputs and outputs used in "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS). The data.zip file contains the data used to drive the model simulations. The model.zip file contains the XML input file for ATS. The notebook.zip file contains the Jupyter notebooks for pre- and post- processing model results. The figures.zip file contains the raw figures associated with the manuscript.Meteorological forcing plays a critical role in accurately simulating the watershed hydrological cycle. With the advancement of high-performance computing and the development of integrated watershed models, simulating the watershed hydrological cycle at high temporal (hourly to daily) and spatial resolution (10s of meters) has become efficient and computationally affordable. These hyperresolution watershed models require high resolution of meteorological forcing as model input to ensure the fidelity and accuracy of simulated responses. In this study, we utilized the Advanced Terrestrial Simulator (ATS), an integrated watershed model, to simulate surface and subsurface flow and land surface processes using unstructured meshes at the Coal Creek Watershed near Crested Butte (Colorado). We compared simulated watershed hydrologic responses including streamflow, and distributed variables such as evapotranspiration, snow water equivalent (SWE), and groundwater table driven by three publicly available, gridded meteorological forcings (GMFs) -- Daily Surface Weather and Climatological Summaries (Daymet), Parameter-elevation Regressions on Independent Slopes Model (PRISM), and North American Land Data Assimilation System (NLDAS). By comparing various spatial resolutions (ranging from 400 m to 4 km) of PRISM, the simulated streamflow only becomes marginally worse when spatial resolution of meteorological forcing is coarsened to 4 km (or 30% of the watershed area). However, the 4 km resolution has much worse performance than finer resolution in spatially distributed variables such as SWE. Using temporally disaggregated PRISM, we compared models forced by different temporal resolutions (hourly to daily), sub-daily resolution preserves the dynamic watershed responses (e.g., diurnal fluctuation of streamflow) that are absent in results forced by daily resolution. Conversely, the simulated streamflow shows better performance using daily resolution compared to that using sub-daily resolution. Our findings suggest that the choice of GMF and its spatiotemporal resolution depends on the quantity of interest and its spatial and temporal scale, which may have important implications on model calibration and watershed management decisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Assessment of Fast Reactor Hot Channel Factor Calculation Capability in Griffin and NekRS

The DOE-NE Advanced Modeling and Simulation (NEAMS) program develops, demonstrates, and deploys predictive computational modeling tools for nuclear reactor design and analysis. Under the Fast Reactor Applications activity within the Application Drivers Technical Area of NEAMS, the ability to perform high fidelity “hot channel factor” (HCF) coupled multiphysics simulations has been identified as a valuable simulation capability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Cross-scale modeling of storm surge, tide, and inundation in Mid-Atlantic Bight and New York City during Hurricane Sandy, 2012

Driven by high-resolution NAM (North America Mesoscale) and ECMWF (European Centre for Medium-Range Weather Forecasts) atmospheric fields, a 3-D unstructured grid SCHISM (Semi-implicit Cross-scale Hydroscience Integrated System Model) was applied to simulate total water level in the Mid-Atlantic Bight during Hurricane Sandy in 2012. The simulated storm surge, tide, significant wave height, and wave peak period results compared favorably with NOAA observations along US East Coast, Long Island Sound, and New York Harbor. The maximum total water level at The Battery in the New York Harbor was accurately simulated with an absolute error of less than 0.08 m (out of 2.9 m peak surge) and a timing difference less than 10 min, which includes the effects on the order of 0.1–0.3 m (5–10%) from the coastal setup induced by the surface waves. The scenario comparisons of (1) “NAM” versus “ECMWF”, (2) “2-D” versus “3-D”, and (3) “with” versus “without” wind wave model were examined. The 3-D barotropic model forced by ECMWF including the effects of wind waves performs the best, attributed to wave-induced radiation stress and reduction of bottom stress when the 3-D version is used. Simultaneously, the ELCIRC-Sub model, using a 5-m topography/bathymetry sub-grid and a regular 200-m resolution finite volume computational model grid, was developed to simulate the street-level inundation in New York City. The momentum and mass fluxes calculated by the coarser base grid model were effectively coupled with the sub-grid so that running a full-blown high-resolution base model is not required. The ELCIRC-Sub model uses MPI (Message Passing Interface) parallel computing to enlarge the coverage of the land surface and an efficient non-linear solver to improve the accuracy of the wetting-and-drying processes. The temporal comparisons of modeled water level with NOAA's tidal stations and USGS' rapid-deployed gauges showed overall performance with an average error of 0.1 m. Particularly, ELCIRC-Sub captured the profile of the highest peak surge (3.9 m), which rose 2 m within 3 hour at Kings Point, NY. The spatial comparisons of the modeled peak water level at 80 surveyed locations showed an average error less than 0.13 m. The modeled maximum inundation extent also matched well with an 80% hit rate against FEMA's Hurricane Sandy maximum flooding extent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Fuel Injection Dynamics and Composition Effects on RDE Performance

Rotating detonation engines (RDEs) provide a promising route to substantially increasing cycle efficiency in stationary gas turbines. Much of this increase relies on the ability to achieve consistent pressure gain within the combustor. In particular, the design of injectors that feed fuel and air into the detonation channel plays a crucial role. Such injectors have to ensure proper mixing of fuel and oxidizer, while minimizing backflow of detonation products into the feed plenums, and reduce susceptibility to the complex wave structures that exist within the combustor. From a practical perspective, such RDEs also need to operate with variable fuel composition. When fuel mixtures with components that possess vastly different oxidation pathways and time-scales are used, there could be additional losses through deflagrative burning instead of detonation-driven heat release. Such sensitivity to the complete flow path is akin to the physics of thermoacoustic instabilities in conventional gas turbines. In this sense, RDEs pose a unique research challenge: the performance of the device relies on the small-scale heat release process, which is highly dependent on the flow interactions within the full-scale system. As a result, canonical flow configurations, instrumented with detailed diagnostics or modeled using high-fidelity tools, but only focus on the small-scale processes will not contain the key system-level interactions. At the same time, macroscopic measurements and models that only capture system-level performance will not provide insight into the key sources of pressure losses. These couplings and sensitivities provide a formidable challenge to both experimental and simulation studies of the effects Thus, a joint experimental/computational program designed specifically to address these challenges was undertaken in this program. The focus of this program was on two key topics: a) the interaction between injector flow and the overall wave dynamics within the combustor, and b) the deflagration/detonation structure in multi-component fuels that are of practical interest. Both topics involve interaction of small-scale heat release processes with the geometry-dependent wave structure. Studies focused on the study of full-scale RDE systems, based on a 6-inch conventional annular geometry. Experimentally RDEs were studied using a combination of diagnostics. A combination of optical diagnostics and aero-thermo-acoustic analysis based on a combination of spectral and mode decomposition analysis was used to identify the dynamics of the detonation wave and other secondary waves that exist in the system. These studies have helped the identification and investigation of injector and detonation dynamics arising from coupling, and how they affect RDE mixing, detonation structure, operability and performance. Performance of RDEs was investigated through thrust stand measurements, which was used to evaluate the effective pressure gain generated by the system through the concept of equivalent available pressure. Optical diagnostics were developed and implemented to investigate the distribution of heat release, across the detonation wave. Novel optical diagnostics of NIR imaging was also developed and applied to investigate the high temperature / high pressure distribution across the detonation wave. In order to complement the experiments, the computational tools were geared to simulate the full experimental setup. GPU-based acceleration of the models and computations were developed to enable rapid simulation of the full system. In addition, the use of adaptive mesh refinement, and unstructured grid formulation, enabled the investigation of realistic geometries studied in the laboratory. The simulations produced a wealth of detail on the structure of the detonation wave under different operating conditions. Emphasis was placed on quantifying mixture pre-burning and the impact on wave propagation and structure.

03 NATURAL GAS↗