Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “criticality simulations”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Computational Fluid Dynamics (CFD) Simulation of Hypersonic Turbine-Based Combined-Cycle (TBCC) Inlet Mode Transition

Methods of computational fluid dynamics were applied to simulate the aerodynamics within the turbine flowpath of a turbine-based combined-cycle propulsion system during inlet mode transition at Mach 4. Inlet mode transition involved the rotation of a splitter cowl to close the turbine flowpath to allow the full operation of a parallel dual-mode ramjet/scramjet flowpath. Steady-state simulations were performed at splitter cowl positions of 0deg, -2deg, -4deg, and -5.7deg, at which the turbine flowpath was closed half way. The simulations satisfied one objective of providing a greater understanding of the flow during inlet mode transition. Comparisons of the simulation results with wind-tunnel test data addressed another objective of assessing the applicability of the simulation methods for simulating inlet mode transition. The simulations showed that inlet mode transition could occur in a stable manner and that accurate modeling of the interactions among the shock waves, boundary layers, and porous bleed regions was critical for evaluating the inlet static and total pressures, bleed flow rates, and bleed plenum pressures. The simulations compared well with some of the wind-tunnel data, but uncertainties in both the windtunnel data and simulations prevented a formal evaluation of the accuracy of the simulation methods.

Slater, John W.↗
NASA NTRSDate not supplied

Initial Testing the Exploration Extravehicular Mobility Unit (xEMU) in Lunar Environment Simulation at the Neutral Buoyancy Lab (NBL) in 2022

Understanding how to effectively train for Extravehicular Activities (EVAs) for Artemis missions is critical. Tasked by the new Extravehicular Activity and Human Surface Mobility Program(EHP), a team of stakeholders from engineering and flight operations worked together to test the Exploration Extravehicular Mobility (xEMU) in the Neutral Buoyancy Lab(NBL) simulating the lunar environment. The NBL is a pool in which astronauts are typically trained to complete simulated International Space Station (ISS) EVA tasks in a neutral weigh-out. The overall focus for this test series was evaluating how well the NBL could be used as an effective simulated environment for Artemis lunar EVAs. New NBL support hardware was manufactured to support this test series and update the interfaces to be more representative of xEMU hardware such as a new NBL Portable Life Support System(PLSS)mock-up and umbilical, lights and cameras, and display and control unit (DCU).To simulate the 1/6th-gravity (1/6-g) environment, extra weights were added around the suit to create an accurate center of gravity representation. The Partial Gravity Weigh-out Stand (PGWS) predicted where and how much weight would be required for the 1/6th-G simulation and predicted the center of gravity placement. After the weigh-out attempt was completed, the suited participant would complete a set of directed tasks such as kneeling, object pick-up, and walking to evaluate the weigh-out. To create the lunar environment, sand, boulders, and ramps were added to the pool floor. If time allowed after weigh-out was determined acceptable, the flight operations team would direct the suited participants to complete simulated Artemis EVA tasks using flight-like tools and tool cart. The team also evaluated turning out the overhead lights in the NBL for a portion of some of the test runs and using a solar simulator to create amore realistic lunar south pole lighting environment. The test series completed one unmanned test and nine manned tests. With four of the tests, dual xEMU suited operations were completed with subjects in both the small and large xEMU sizes. Overall, the NBL was able to create an acceptable lunar landscape and ability to have two subjects work side by side in pressurized xEMUs which currently cannot be achieved in other test environments. Achieving a stable 1/6th-G lunar weigh-out and clear, consistent dual suit communications was a challenge for this series. Future work is planned to continue to improve this simulation environment.

spacesuit↗
NASA NTRSDate not supplied

Initial Testing of the Exploration Extravehicular Mobility Unit (xEMU) in Lunar Environment Simulation at the Neutral Buoyancy Lab (NBL) in 2022

Understanding how to effectively train for Extravehicular Activities (EVAs) for Artemis missions is critical. Tasked by the new Extravehicular Activity and Human Surface Mobility Program(EHP), a team of stakeholders from engineering and flight operations worked together to test the Exploration Extravehicular Mobility (xEMU) in the Neutral Buoyancy Lab(NBL) simulating the lunar environment. The NBL is a pool in which astronauts are typically trained to complete simulated International Space Station (ISS) EVA tasks in a neutral weigh-out. The overall focus for this test series was evaluating how well the NBL could be used as an effective simulated environment for Artemis lunar EVAs. New NBL support hardware was manufactured to support this test series and update the interfaces to be more representative of xEMU hardware such as a new NBL Portable Life Support System(PLSS)mock-up and umbilical, lights and cameras, and display and control unit (DCU).To simulate the 1/6th-gravity (1/6-g) environment, extra weights were added around the suit to create an accurate center of gravity representation. The Partial Gravity Weigh-out Stand (PGWS) predicted where and how much weight would be required for the 1/6th-G simulation and predicted the center of gravity placement. After the weigh-out attempt was completed, the suited participant would complete a set of directed tasks such as kneeling, object pick-up, and walking to evaluate the weigh-out. To create the lunar environment, sand, boulders, and ramps were added to the pool floor. If time allowed after weigh-out was determined acceptable, the flight operations team would direct the suited participants to complete simulated Artemis EVA tasks using flight-like tools and tool cart. The team also evaluated turning out the overhead lights in the NBL for a portion of some of the test runs and using a solar simulator to create amore realistic lunar south pole lighting environment. The test series completed one unmanned test and nine manned tests. With four of the tests, dual xEMU suited operations were completed with subjects in both the small and large xEMU sizes. Overall, the NBL was able to create an acceptable lunar landscape and ability to have two subjects work side by side in pressurized xEMUs which currently cannot be achieved in other test environments. Achieving a stable 1/6th-G lunar weigh-out and clear, consistent dual suit communications was a challenge for this series. Future work is planned to continue to improve this simulation environment.

spacesuit↗
DOE OSTI2022

A Fast and Accurate Transient Stability Assessment Method Based on Deep Learning: WECC Case Study

Transient stability is one of the critical aspects of power system stability assessment. The increasing integration of inverter-based resources and the retirement of conventional synchronous generators result in the decreasing system inertia and growing complexity of system operating conditions. Using a few selected typical operating conditions cannot guarantee system transient stability in all operating conditions, and the time-domain simulation of all operating conditions requires tremendous time and is often infeasible. This paper proposes a more efficient transient stability assessment method based on deep learning. The binary search method is used to determine the critical clearing time (CCT) in creating training databased by time-domain simulation. This method is fast and accurate with 1 ms resolution. The buses whose CCTs are lower than 200 ms are considered critical buses. Buses close to each other are grouped based on their mutual admittance matrix to reduce the search space of the critical buses. This paper also proposes the generator feature normalization based on the physical model. Case study on the reduced 240-bus WECC system model demonstrates that the proposed method can predict CCT accurately and efficiently.

critical clearing time↗
DOE OSTI2023

The Local Electronic Structure of Supercritical CO 2 from X-ray Raman Spectroscopy and Atomistic-Scale Modeling

Supercritical CO 2 is encountered in several technical and natural systems related to biology, geophysics, and engineering. While the structure of gaseous CO 2 has been studied extensively, the properties of supercritical CO 2 , particularly close to the critical point, are not well-known. Here, in this work, we combine X-ray Raman spectroscopy, molecular dynamics simulations, and first-principles density functional theory (DFT) calculations to characterize the local electronic structure of supercritical CO 2 at conditions around the critical point. The X-ray Raman oxygen K-edge spectra manifest systematic trends associated with the phase change of CO 2 and the intermolecular distance. Extensive first-principles DFT calculations rationalize these observations on the basis of the 4sσ Rydberg state hybridization. X-ray Raman spectroscopy is found to be a sensitive tool for characterizing electronic properties of CO 2 under challenging experimental conditions and is demonstrated to be a unique probe for studying the electronic structure of supercritical fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2017

Aeroacoustic Simulations of a Nose Landing Gear with FUN3D: A Grid Refinement Study

A systematic grid refinement study is presented for numerical simulations of a partially-dressed, cavity-closed (PDCC) nose landing gear configuration that was tested in the University of Florida's open-jet acoustic facility known as the UFAFF. The unstructured-grid flow solver FUN3D is used to compute the unsteady flow field for this configuration. Mixed-element grids generated using the Pointwise (Registered Trademark) grid generation software are used for numerical simulations. Particular care is taken to ensure quality cells and proper resolution in critical areas of interest in an effort to minimize errors introduced by numerical artifacts. A set of grids was generated in this manner to create a family of uniformly refined grids. The finest grid was then modified to coarsen the wall-normal spacing to create a grid suitable for the wall-function implementation in FUN3D code. A hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) turbulence modeling approach is used for these simulations. Time-averaged and instantaneous solutions obtained on these grids are compared with the measured data. These CFD solutions are used as input to a FfowcsWilliams-Hawkings (FW-H) noise propagation code to compute the farfield noise levels. The agreement of the computed results with the experimental data improves as the grid is refined.

Vatsa, Veer N.↗
NASA NTRS2005

Bio-Optics of the Chesapeake Bay from Measurements and Radiative Transfer Calculations

We combined detailed bio-optical measurements and radiative transfer (RT) modeling to perform an optical closure experiment for optically complex and biologically productive Chesapeake Bay waters. We used this experiment to evaluate certain assumptions commonly used when modeling bio-optical processes, and to investigate the relative importance of several optical characteristics needed to accurately model and interpret remote sensing ocean-color observations in these Case 2 waters. Direct measurements were made of the magnitude, variability, and spectral characteristics of backscattering and absorption that are critical for accurate parameterizations in satellite bio-optical algorithms and underwater RT simulations. We found that the ratio of backscattering to total scattering in the mid-mesohaline Chesapeake Bay varied considerably depending on particulate loading, distance from land, and mixing processes, and had an average value of 0.0128 at 530 nm. Incorporating information on the magnitude, variability, and spectral characteristics of particulate backscattering into the RT model, rather than using a volume scattering function commonly assumed for turbid waters, was critical to obtaining agreement between RT calculations and measured radiometric quantities. In situ measurements of absorption coefficients need to be corrected for systematic overestimation due to scattering errors, and this correction commonly employs the assumption that absorption by particulate matter at near infrared wavelengths is zero.

Tzortziou, Maria↗
DOE OSTI2023

Comparing quantile regression forest and mixture density long short-term memory models for probabilistic post-processing of satellite precipitation-driven streamflow simulations

Abstract. Deep learning (DL) and machine learning (ML) are widely used in hydrological modelling, which plays a critical role in improving the accuracy of hydrological predictions. However, the trade-off between model performance and computational cost has always been a challenge for hydrologists when selecting a suitable model, particularly for probabilistic post-processing with large ensemble members. This study aims to systematically compare the quantile regression forest (QRF) model and countable mixtures of asymmetric Laplacians long short-term memory (CMAL-LSTM) model as hydrological probabilistic post-processors. Specifically, we evaluate their ability in dealing with biased streamflow simulations driven by three satellite precipitation products across 522 nested sub-basins of the Yalong River basin in China. Model performance is comprehensively assessed using a series of scoring metrics from both probabilistic and deterministic perspectives. Our results show that the QRF model and the CMAL-LSTM model are comparable in terms of probabilistic prediction, and their performances are closely related to the flow accumulation area (FAA) of the sub-basin. The QRF model outperforms the CMAL-LSTM model in most sub-basins with smaller FAA, while the CMAL-LSTM model has an undebatable advantage in sub-basins with FAA larger than 60 000 km2 in the Yalong River basin. In terms of deterministic predictions, the CMAL-LSTM model is preferred, especially when the raw streamflow is poorly simulated and used as input. However, setting aside the differences in model performance, the QRF model with 100-member quantiles demonstrates a noteworthy advantage by exhibiting a 50 % reduction in computation time compared to the CMAL-LSTM model with the same ensemble members in all experiments. As a result, this study provides insights into model selection in hydrological post-processing and the trade-offs between model performance and computational efficiency. The findings highlight the importance of considering the specific application scenario, such as the catchment size and the required accuracy level, when selecting a suitable model for hydrological post-processing.

Geology↗
NASA NTRS2017

Representing Water Scarcity in Future Agricultural Assessments

Globally, irrigated agriculture is both essential for food production and the largest user of water. A major challenge for hydrologic and agricultural research communities is assessing the sustainability of irrigated croplands under climate variability and change. Simulations of irrigated croplands generally lack key interactions between water supply, water distribution, and agricultural water demand. In this article, we explore the critical interface between water resources and agriculture by motivating, developing, and illustrating the application of an integrated modeling framework to advance simulations of irrigated croplands. We motivate the framework by examining historical dynamics of irrigation water withdrawals in the United States and quantitatively reviewing previous modeling studies of irrigated croplands with a focus on representations of water supply, agricultural water demand, and impacts on crop yields when water demand exceeds water supply. We then describe the integrated modeling framework for simulating irrigated croplands, which links trends and scenarios with water supply, water allocation, and agricultural water demand. Finally, we provide examples of efforts that leverage the framework to improve simulations of irrigated croplands as well as identify opportunities for interventions that increase agricultural productivity, resiliency, and sustainability.

Irrigated Agriculture; Crop Model; Hydrologic Mode↗
DOE OSTI2022

Reactive Transport Modeling of Mineral Precipitation and Carbon Trapping in Discrete Fracture Networks

In this study we use numerical experiments to analyze reactive flow and transport behavior in discrete intersecting fracture networks, focusing on (a) how reaction-induced changes in physical and chemical properties affect flow connectivity and (b) how fracture networks developed in the Earth's critical zone contribute to carbon sequestration via mineral weathering reactions. In the first part of the study, we used two-dimensional reactive flow and transport simulations to analyze the impacts of mixing in a natural discrete fracture network. We concluded that reaction-induced changes can substantially alter the flow connectivity, especially at fracture intersections. The second set of simulations considered the problem of natural weathering of fractured mafic and ultramafic rocks in the partially saturated Earth's critical zone as a function of infiltration rates, fracture permeability, and partially saturated flow parameters. As a model system, we considered an incongruent reaction network with dissolution of forsterite and precipitation of magnesite. The behavior is complex in terms of the rate-controlling processes because of the multicomponent nature of the system as shown by the grid Peclet number: the CO 2 behavior is gas diffusion-controlled in the partially saturated zone, while the rate of water flow via the Damkӧhler number controls Mg 2+ transport through the fracture network. The amounts of carbon that can be trapped are modest, but the naturally fractured domain considered here provides a useful “base case” against which various engineered solutions can be compared.

58 GEOSCIENCES↗
DOE OSTI2021

Computational Fluid Dynamics Simulation of Compressible Non-Newtonian Biomass in a Compression-Screw Feeder

Compression-screw feeders play a critical role in biorefineries to transport lignocellulosic biomass feedstocks from the feeding hoppers to biomass-conversion reactors in order to pretreat and convert biomass to hydrocarbon liquid biofuels and other power and energy resources. One of the main challenges in the operation of screw feeder is plugging and jamming of compressed biomass with high concentration of insoluble solids. The focus of this paper is to numerically investigate the screw feeder at these challenging operating conditions and help with the optimization of the screw feeder design to avoid operation failure. In this work a customized CFD model based on open-source OpenFOAM package [1] was developed to simulate the concentrated biomass as a highly viscous non-Newtonian fluid in the screw feeder. The biomass is modeled as a single-phase compressible Bingham fluid with a plastic viscosity as well as a density-dependent yield stress. The compressibility formulation (pressure-dependent density) and the density-dependent yield stress formulation in the governing equations follow the suggestions from a recent study by Duncan et al. [2]. A pilot-scale hopper/screw feeding system at NREL [3] is used to compare the experimental observations with our simulation results. The auger is 280 mm long and tapered with outer diameter changing from 80 mm to 35 mm. The auger rotates from 10 to 60 rpm in a conical throat which contains anti-rotational bars. The simulations predicted the required torque for the screw feeder and the pressure increase at the exit for biomass feedstocks with various fluid viscosity properties and auger rotating speeds. The analysis of the stress forces helped to identify the critical conditions were the screw feeder excessive wear or jamming could occur.

biofuels↗
DOE OSTI2022

Lignin-based jet fuel and its blending effect with conventional jet fuel

Sustainable aviation fuels (SAFs) must demonstrate specific physical and chemical properties as well as material compatibility (i.e., seal swell) to be used as an aviation turbine fuel. Several alternative jet fuels incorporated in ASTM D7566 are comprised mainly of n/iso-alkanes and can only be blended up to 50 vol% due to material compatibility and density issues. Prior work illustrated the ability of cycloalkanes to replace aromatics’s role in material compatibility. Here, we report the first archival documentation of a feedstock and chemical process to yield a product composition able to complement existing SAF ASTM D7566 annexes. A lignin-based jet fuel (LJF) blend components is generated and composed of mostly C7-C18 mono, di, and tri-cycloalkanes. The neat LJF was blended with conventional jet fuel at 10 vol% (LJF blend) to simulate an ASTM “Fast Track” evaluation process. Fuel properties that are critical to engine operability (ATSM D4054 Tier 3 & 4) were either predicted or experimental tested based on the volume availability. All LJF blend’s critical properties fall within the experience range of conventional jet fuel, with o-ring swelling exceeding the typical range of conventional fuels. Here, these results in total support the potential use of this LJF pathway to complement other SAF pathways and achieve 100% drop-in SAF.

09 BIOMASS FUELS↗
DOE OSTI2020

Applying a Template of Expected Uncertainties to Updating 239 Pu(n,f) Cross-section Covariances in the Neutron Data Standards Database

Templates of uncertainties expected in specific measurement types were recently developed. One aim of these templates is to help evaluators in identifying (1) missing or suspiciously low uncertainties and (2) missing correlations between uncertainties of the same and different experiments, when estimating covariances for experimental data employed in their evaluations. These templates also provide realistic estimates of standard deviations and correlations for a particular uncertainty source and measurement type that can be used by evaluators in situations where they are not supplied by the experimenters. This information allows for a more comprehensive uncertainty analysis across all measurements considered in an evaluation and, thus, more realistic evaluated covariances. Here, in this work, we extend a template that is applicable to uncertainties expected in neutron-induced fission, (n,f), cross-section measurements. It is applied to improving covariances of 239 Pu(n,f) cross-section measurements in the database underlying the Neutron Data Standards evaluations. This particular example was chosen since this evaluation is primarily based on experimental information. Also, some uncertainties of individual 239 Pu(n,f) cross-section experiments in this database were suspected to be underestimated. The evaluated uncertainties obtained after updating the covariances in the database by means of the template indeed do increase compared to their original values. Even more importantly, the evaluated mean values change noticeably. These modified cross sections impact application calculations significantly, as is demonstrated by employing them in simulations of the effective neutron multiplication factor for a few selected critical assemblies. However, this updated evaluated 239 Pu(n,f) cross section should not be interpreted as the final one that should replace values of the current Neutron Data Standards project. Evaluations for the Neutron Data Standards of the 239 Pu(n,f) cross section must be linked to many other observables included in the associated database, most notably to cross sections for 235 U(n,f), but also to those for 10 B(n, α ), 6 Li(n,t), 238 U(n,f), and 238 U(n, γ ), because of included measurements of the 239 Pu(n,f) cross section that appear as ratios to these reactions. Some of these other reactions are correlated to further observables in the database. Hence, updating uncertainties of data sets of any of these observables can potentially impact the 239 Pu(n,f) cross section. Uncertainties for all measurements of these linked physical observables have to be updated before a comprehensive evaluation of the 239 Pu(n,f) cross section and its corresponding uncertainties can be provided.

239Pu Neutron-induced fission cross-section↗
DOE OSTI2020

Risk-Oriented PMU Placement Approach in Electric Power Systems

The objective of traditional PMU placement approaches is to use the minimal number of PMUs to cover system complete observability. Very few consider the probability and the consequence severity of observability loss. For example, following line trip, grid topology changes, which may results in observability loss and overloaded lines. As long as the overloaded lines are still observable, then we can take action in time to correct the system; and the consequence of the observability loss is not severe, although some parts of the system is not observable. Otherwise, if the overloaded lines become unobservable, then we lose the situation awareness of the system and the consequence of the observability loss is severe. To address the above points in PMU placement, this paper proposes a risk-oriented PMU placement approach. The proposed approach uses two indices, namely the probability of observability loss and its corresponding consequence of severity to design the optimal PMU placement strategy. We show that our proposed approach is able to capture the critical events while the traditional approaches fail to do so. Extensive simulation results carried out on the IEEE 39-bus and 118-bus systems demonstrate its effectiveness.

observability, risk of observability loss, PMUs pl↗
DOE OSTI2026

Direct numerical simulations of activation and deactivation in turbulent atmospheric clouds

Significant knowledge gaps remain in our understanding of turbulence–cloud–aerosol interactions in the Earth's atmosphere, and direct numerical simulation (DNS) has increasingly become an indispensable tool to fill such critical knowledge gaps. Here, this study is an extension of our previous DNS model [Gao et al., J. Geophys. Res.: Atmos., 123(4), 2194–2214 (2018)], with a focus on the activation of aerosol particles into cloud droplets and deactivation of cloud droplets into aerosol particles in a microscale cloud environment. The effects of turbulence intensity, particle curvature, and solute, as well as the initial distributions of the aerosol particles (monodisperse or polydisperse) are investigated. The governing equations for the flow of air, temperature, and water vapor mixing ratio are solved numerically in the Eulerian fashion, assuming homogeneous and isotropic turbulence. The dynamics of the aerosol and cloud particles are calculated with the Lagrangian particle tracking method. The results show that the deviations of the thermodynamic variables from their respective means are significantly reduced, the activation process appears to be delayed, and the deactivation process occurs more rapidly, as the turbulence intensity is increased. The inclusion of particle curvature and solute effects, as well as polydispersity, tends to retard the activation of aerosols into cloud droplets. It is also observed that fluctuations in supersaturation broaden the spread of particle radii, and the broadening is followed by a narrowing as turbulent homogenization reduces thermodynamic fluctuations over time.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Measurements of Density and Sound Speed in Mixtures Relevant to Supercritical CO2 Cycles

Abstract The objective of this research is to validate properties of mixtures relevant to supercritical carbon dioxide (sCO2) power cycles. Direct-fired sCO2 cycles are a promising technology for the future power generation systems. The working fluid of sCO2 cycles will be near and above critical point of CO2. One of the challenges is that the simulation of mixtures should consider real gas behavior. Expected operating conditions of Allam cycles reach up to 300 bar and 1000 °C. Characterizing the mixtures at the extreme conditions is an important issue in current researches and industrial applications. Thermophysical properties of mixtures may be beyond the valid range of the widely used database such as NIST REFPROP. Experimental data of mixture properties in the literature are limited which is necessary to develop high-fidelity design tools for sCO2 power cycles. We measured the density and sound speed of several multi-component mixtures. A temperature-controlled high-pressure test cell was used for the density measurements. Sound speed was measured by resonant frequency detection using an external speaker and a piezoelectric pressure sensor. Mixtures studied in this work include carbon dioxide, methane, oxygen, and water vapor. Properties of pure CO2 were measured to show the validity of our technique. Compositions were selected to be close to frozen mixtures at the inlet, mid-progress, and exhaust conditions of a model sCO2 combustor in the previous numerical simulation work. Corresponding reaction progress variables (RPV) were RPV = 0, 0.5, and 1. Temperature and pressure conditions of experiments are 310–450 K and 0–150 bar. In our study, density and sound speed from the NIST REFPROP database agree with experimental measurements within the range of our measurement uncertainties.

Energy & Fuels↗
DOE OSTI2020

A constitutive model for sheared dense suspensions of rough particles

In a concentrated suspension, particles come into contact due to the presence of asperities on their surfaces. As a result, the contact forces and interparticle friction become one of the important factors governing the rheology of rigid particle suspensions at high concentrations. We show that a load-dependent friction model can be used to reproduce the experimentally observed shear thickening [ST—continuous and discontinuous (DST)] behavior with an increase in the asperity size. Increasing the particle surface roughness size leads to a decrease in the critical shear rate above which shear-thickening takes place, a reduction in the critical volume fraction for DST and an increase in the viscosity jump across non-shear-thickened to shear-thickened regime. In this paper, we propose a constitutive model to quantify the effect of increasing the roughness size on the rheology of dense ST suspensions as well as on the critical shear rate for ST and the critical volume fraction for DST. We fit this model to our simulation data for stress controlled shear flow of dense rough particle suspensions. Once the fitting is complete, these equations are used to predict exact volume fractions and shear stress values for transitions between three regimes on the shear stress-shear rate flow state diagram for different roughness values. Furthermore, the results of this study can be used to tune the particle surface roughness for manipulating the dense suspension rheology according to different applications.

42 ENGINEERING↗
DOE OSTI2025

Openpronghorn

OpenPronghorn is a simulation tool specifically tailored for modeling thermal-hydraulic phenomena in advanced nuclear reactors. It is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source platform that facilitates the development of high-performance scientific computing applications. OpenPronghorn solves the Navier-Stokes equations, which describe the conservation of mass, momentum, and energy in fluid flows, using the finite volume numerical method. The code supports a wide range of fluid flow conditions that are applicable to nuclear reactors, including incompressible and weakly compressible flows, as well as single-phase and multiphase flows. It is capable of modeling diverse flow regimes, including laminar and turbulent flows, using various turbulence models such as the standard k-epsilon models, the v2f model, and the mixing length model. For multiphase flows, OpenPronghorn employs a mixture a Eulerian modeling approach with mixture, drift-flux, and full Eulerian models, and includes open-sourced interfacial transfer correlations for drag, exchange, and heat transfer coming from the scientific literature. OpenPronghorn's modular design allows it to handle multiscale simulations, ranging from detailed Reynolds-Averaged Navier Stokes (RANS) simulations to coarse-mesh and lumped parameter models. This flexibility enables users to perform high-fidelity simulations of specific reactor components as well as system-level analyses of entire reactor circuits. The code can be coupled with other MOOSE-based tools using the MultiApp system, allowing for the transfer of coupling quantities such as mass flow rates, heat fluxes, and boundary conditions between different simulation scales. One of the main features of OpenPronghorn is the it includes built-in validation cases from the open-source scientific literature and supports the implementation of user-defined models and correlations through MOOSE's FunctorMaterial system. OpenPronghorn is designed to be computationally efficient, leveraging the SIMPLE projection method for large-scale problems, and can be run on high-performance computing systems to handle the extensive computational demands of detailed reactor simulations. Overall, OpenPronghorn is a versatile and robust tool that provides critical insights into the thermal-hydraulic behavior of advanced nuclear reactors, supporting the design, safety, and optimization of next-generation nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗