Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow Conditions”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Assessment of Model Hub Height Wind Speed Performance Using DOE Lidar Buoy Data

As the offshore wind industry continues to gain momentum both in the U.S. and globally, reanalysis products remain essential tools for the lifecycle of a wind farm. From project siting and production estimates to determining construction and maintenance windows, reanalysis products play an important role. These regional and global long-term models also serve as inputs to higher resolution models. With such wide-ranging impacts on the offshore wind industry, the need for validation of reanalysis products is strong, particularly at typical turbine hub heights. Pacific Northwest National Laboratory (PNNL) operates two AXYS WindSentinelTM buoys for the U.S. Department of Energy (DOE) in order to collect meteorological and oceanographic data, including hub height wind speed, in areas of interest for offshore wind development. This work uses DOE buoy observations off the coasts of New Jersey and Virginia to evaluate the performance of Modern Era Retrospective Analysis for Research and Applications-2 (MERRA-2), the North American Regional Reanalysis (NARR), the European Center for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5th generation (ERA5), and the National Oceanic and Atmospheric Administration (NOAA) assimilation system Rapid Refresh (RAP). Biases and degrees of correspondence are determined for each reanalysis product in order to provide insights on the performance and uncertainty for long-term wind resource characterization. With the baseline metrics determined, this work proceeds with an investigation into the sources of large deviations between the modeled and observed hub height wind speeds during the East Coast lidar buoy deployments. Consistent sources of reanalysis model error include stably stratified flow conditions, high wind shear, flow parallel to the coastline, and tropical and winter storms. For the near shore New Jersey location, offshore winds (i.e. blowing from land to water) tend to be associated with model overestimation of observed hub height wind speed, while winds parallel to the coastline are correlated with model underestimation of wind speed. Similarly, near the Virginia buoy, during southerly winds (parallel to the coastline) models underestimate the observed wind speed due to coastal upwelling. Reanalysis model wind speeds underestimate the hub height winds with increasing stable atmospheric stratification. Reanalysis model biases are seen when the observed wind speeds approach or exceed typical turbine cut-out speeds (around 20 m s-1). Large model errors during winter and tropical storm events, high wind speed events, and ramp events are often the result of wind speed timing and magnitude offsets. Overall, ERA5 provides the most successful representation of observed offshore hub height wind speeds at the U.S. East Coast buoy locations and is therefore best suited as an input boundary condition for model case studies. ERA5 performs well for a variety of atmospheric phenomena, including storms and sea breezes, however RAP, the next most successful reanalysis model, performs best in capturing the frequency of ramp events.

17 WIND ENERGY↗
DOE OSTI2024

Active learning for the design of polycrystalline textures using conditional normalizing flows

Generative modeling has opened new avenues for solving previously intractable materials design problems. However, these new opportunities are accompanied by a drastic increase in the required amount of training data. This is in stark juxtaposition to the high expense and difficulty in curating such large materials datasets. In this work, we propose a novel framework for integrating generative models within an active learning loop. Further, this enables the training of generative models with datasets significantly smaller than what has previously been demonstrated, providing a direct route for their application in data constrained environments. The functionality of this framework is then demonstrated by addressing the challenge of designing polycrystalline textures associated with target anisotropic mechanical properties. The developed protocol exhibited a cost reduction between 14 to 18 times over a randomly sampled experimental design.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electrochemical Control for Corrosion in Molten Chlorides During CSP Plant Operation

The Liquid Pathway of the Concentrating Solar Power Generation 3 (CSP Gen3) program proposed low-cost molten chloride salt for energy storage. However, online corrosion control was identified as a remining major risk of the Liquid Pathway approach. This project addressed that risk. Electrochemical solutions for corrosion mitigation during CSP plant operation were investigated and their feasibility and scalability were evaluated. The leading cause of corrosion in molten chloride salt systems was identified as corrosive impurities that form within the salt upon exposure to trace amounts of air and moisture. Leveraging electrochemistry, reduction/oxidation reactions can be employed to remove these corrosive impurities. In Phase 1 of this project, a bench-scale batch electrochemical reactor was designed, fabricated, and used to assess the kinetics and thermodynamics of electrochemical salt purification. In Phase 2, a laboratory-scale flow reactor was designed, fabricated, and used to assess the efficacy of the electrochemical method under flowing conditions. Results show that under proposed operating conditions for the Liquid Pathway Gen3 Pilot Plant, the electrochemical method is significantly more effective at removing impurities than alternative chemical and thermal methods, and that the electrochemical method produces less harmful byproducts. A key advance made in the course of this project was the development of a 2-electrode method for electrochemical purification that is more scalable than previously developed 3 electrode methods. This novel method is based on Magnesium (Mg) electrowinning. A provisional patent based on this invention has been submitted (USPTO Application No. 63/480,355). Additional key advances made during this project include assessment of the effect of dissimilar alloys on corrosion, kinetic and thermodynamic evaluation of thermolysis reactions of impurities within the molten salt, characterization of byproducts of purification reactions, and generation of IP focused on isolating value-added products using molten salt-based electrochemistry that could be deployed to valorize the process (USPTO Application No. 63/478,806). Ultimately, this project represented a step toward feasibility of Liquid Pathway Gen3 CSP. The method developed under this project could significantly reduce capital expenses and operating costs and increase plant profitability by enabling use of less expensive alloys, decreasing maintenance, and increasing plant longevity. Key focus areas for follow-on work have been identified as 1) evaluation of the efficacy of the electrochemical method under turbulent conditions in a larger flow system, such as the FASTR loop, 2) development of methods for removal of purification byproducts, 3) modeling pilot and industrial scale performance of electrochemical salt purification during plant operation and 4) further assessment of the effect of impurities on salt vapor phase.

14 SOLAR ENERGY↗
DOE OSTI2024

Recovery of Oxpure 612C-50 Coconut Carbon Particles Using Bump Arrays

Cyanide is used to leach gold from crushed ore (solid matrix) into a gold cyanide solution. The gold is extracted from the cyanide solution by adsorption onto activated carbon. The gold extraction process occurs when the activated coconut carbon is placed into tanks that contain the gold cyanide solution either in a batch process or into a continuous flow circuit. The coconut carbon is then removed from the solution for gold recovery. The use of the mesofluidic separation system (a deterministic lateral displacement system) can significantly reduce costs and waste generated from the process of removing the coconut carbon from the solution. In the coconut carbon recovery process used at many gold mines in Nevada, the gold is attached to carbon particles that are still in the cyanide liquor. The mesofluidic system can rapidly remove the gold bearing carbon particles from the cyanide liquor quickly and with no operating costs. The benefits include: • Reduce operational costs related to filtration and hydro-cyclones • Reduced acid usage for the elution gold stripping phase as only 25% of the fluids will follow the carbon particles to the final express lane when two mesofluidic systems are used in series This alternative particle removal techniques would be advantageous. A promising technique for removing larger particles from slurries is mesofluidic separation, similar to deterministic lateral displacement arrays or “bump” arrays [1] but operates at much larger flow rates. As described by Pease, et al. [2], “Bump arrays in deterministic lateral displacement devices separate large particles from small particles using arrays of staggered posts. Large particles, defined as those with radii larger than the distance between the edge of a post and the stagnation streamline from the next downstream post, must bump toward one side of the device, whereas particles smaller than this distance slalom from entrance to exit without net lateral displacement.” Unlike filters or sieves, the posts in the arrays that cause separation do not block or occlude particle flow but work because particles go around the posts. Unlike hydrocyclones separation, separation is not driven by particle density, because gravitational forces are unimportant to mesofluidic separation. Burns, et al., and Pease, et al., have shown that particles may be separated from complex suspensions, under turbulent flow conditions, and at industrially important flow rates [3-5]. They have also recently shown that mesofluidic devices may be arranged in series to increase separation performance [6]. In this paper, we evaluate the mesofluidic system for the separation of commercial OxPure GR 612 Charred Coconut shell particles as a proof of principle test for the rapid separation using the 1500 micron cut mesofluidic separator. We first describe the experimental system and conditions. We then present the experimental results. [1] Huang, L.R., E.C. Cox, R.H. Austin, J.C. Sturm. 2004. Continuous particle separation through deterministic lateral displacement. Science, 304, 987–990. https://doi.org/10.1126/science.1094567. [2] Pease L.F., J.A. Bamberger, C.A. Burns, and M.J. Minette. 2021. Large Particle Separation from Non-Newtonian Slurries using Bump Arrays. In Proceedings of the ASME 2021 Fluids Engineering Division Summer Meeting FEDSM2021-65904, V003T08A023. New York, New York: ASME. doi:10.1115/FEDSM2021-65904 [3] Burns C.A., T.G. Veldman, J. Serkowski, R.C. Daniel, X.-Y. Yu, M.J. Minette, L.F. Pease. 2021. Mesofluidic separation versus dead-end filtration. Separation and Purification Technology, 254, 117256. https://doi.org/10.1016/j.seppur.2020.117256. [4] Pease L.F., J.E. Serkowski, T.G. Veldman, J. Williams, X.-Y. Yu, M.J. Minette, J.A. Bamberger, C.A. Burns. 2021. Can Bump Arrays Separate Particles from Turbulent Flows?. In Proceedings of the ASME 2021 Fluids Engineering Division Summer Meeting FEDSM2021-67696, V003T08A024. New York, New York: ASME. doi:10.11

mesofluidic separation, slurry, bump array↗
DOE OSTI2026

Performance of pumped counterflow virtual impactors to study aerosol interactions with laboratory generated warm clouds

Pumped Counterflow Virtual Impactors (PCVI) are designed to separate aerosols based on aerodynamic diameter, which is particularly useful for isolating cloud droplets and ice crystals from smaller particles. However, the PCVI transmission efficiency (TE) values reported in the literature show considerable variability, and little information is available on the TEs for cloud liquid droplets. Here, we determined the optimal flow conditions for PCVI sampling for different activation ratios from the MTU Pi-cloud chamber, highlighting the conditions that maximize droplet residual sampling while minimizing interstitial transmission. Even for a lower limit cloud activation ratio = 1:10, an add flow of 1.5 LPM achieves a residual fraction greater than 0.80, with a droplet TE of ∼15% for typical Pi-Chamber droplets with diameters between 3.5 to 10 µm. This framework can be adapted for field cloud measurements based on the cloud conditions of interest. TE for cloud droplets under flow conditions at which ∼99% of the unactivated submicron particles were removed was substantially lower between 8 and 16% than for supermicron dry particles (∼30 to 40% for 3 µm polystyrene latex spheres). In addition, we conducted a comparison of the performances of three similar PCVI units to assess repeatability for dry aerosols and cloud droplets. All three PCVIs performed similarly in rejecting submicron particles smaller than the desired cutoff diameter. One of the units showed lower TEs due to a misalignment of the internal orifice. We discuss a procedure that improved alignment and performance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Behavior and mechanisms of Doppler wind lidar error in complex terrain: stable flow case study at Perdigão

A numerical experiment is carried out investigating the magnitude of biases in ground-based lidar measurements in complex flow conditions. Biases assessed include those arising from flow curvature and from the interaction of turbulence with the wind field reconstruction (WFR) algorithms used by a WindCube lidars and anemometers. RANS-CFD and WRF-LES simulations were performed for the Perdig˜ao Field Experiment site for a range of atmospheric conditions. Virtual anemometer and lidar data were generated for four locations: two near exposed ridge tops and two in low-speed regions in the valley. The LES data at these four locations show that the scalar inflation terms (the relation between scalar and vector averaged wind speed) for virtual lidar and virtual cups agree very well with predictions using perturbation theory. While the lidar errors vary greatly with location and height, the contribution from the flow curvature tends to be larger than the differences arising from scalar inflation. For one lidar/mast pair near the ridge top, comparisons between simulations and measurements are carried out for a resonant mountain wave event on June 14th, 2017, and for the whole duration of the Perdigão campaign for winds perpendicular to the ridges. The lidar error during the mountain wave, a period of strong stability and low inversion height, is significantly larger than the campaign average. The sensitivity of the lidar error to atmospheric stability is confirmed by the RANS simulations, which suggests strong sensitivity of flow curvature error to stability conditions and to the shape of the wind speed profile near the top of the boundary layer.

17 WIND ENERGY↗
DOE OSTI2024

Engineering a biopolymer for vascular embolization based on fundamental polymer principals

Herein we describe the engineering of a potential polymer substitute to metal coils for achieving vascular embolization. Our system is based on a methacrylate end-functionalized polyethylene oxide-polypropylene oxide-polyethylene oxide (PEO-PPO-PEO) triblock copolymer which displays thermal-reversible gelation properties and FCC ordering of the original copolymer while also being able to chemically crosslink. This system shows injectability, two-step gelation, chemical crosslinking with a five-fold increase in modulus and is highly resistant to washout under physiological flow conditions. Furthermore, we show that exposure to aqueous flow imparts resistance to swelling where less than 10 % swelling was observed under physiological flow as opposed to more than 60 % under static conditions after 30 days. Auxiliary experiments with PEO-PPO-PEO copolymers of various molecular weights indicate that this phenomenon may possibly be attributed to entanglements known to occur between adjacent micelles, where the reptation time, the time for polymer disentanglement, is long compared to the flow rate.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the interaction of a wind turbine wake with a conventionally neutral atmospheric boundary layer

In this work, we investigate the dynamics of wind turbine tip-vortex breakdown in a conventionally neutral atmospheric boundary layer (ABL). To this end, high-resolution data are collected from large-eddy simulations of a wind turbine operating within a neutral ABL and studied by means of proper orthogonal decomposition (POD) and Fourier analysis. The high resolution of the generated data in both space and time allows us to gain insight into the tip-vortex breakdown mechanisms by (i) capturing the energy modes of the coherent structures, (ii) studying their contribution to the tip-vortex breakdown through their power spectra functions and mean kinetic energy (MKE) flux, and (iii) analysing the growth rate of each contributing perturbation frequency along tip vortices. Our analysis shows that under a fully turbulent scenario, the growth rate of perturbations along the tip vortices is largest for low wave numbers, i.e. long-wave perturbations. Additionally, the MKE flux reaches its highest value at two diameters downstream of the rotor plane, a behaviour that can be attributed to the coexistence of multiple interacting POD modes, with the streamwise vortex roller mode being the primary contributor to the total MKE flux budget, contributing approximately . Finally, comparisons with a laminar, uniform flow scenario subject to a single-frequency perturbation highlight the differences between the two ambient flow conditions. In the non-turbulent, uniform flow scenario, the growth rate attains its maximum value at a wave number corresponding to the out-of-phase mutual-inductance mechanism, whereas the MKE flux exhibits local minima and maxima along the wake and at different downstream locations depending on the perturbation frequency. Our analyses suggest that the breakdown of the wind turbine tip vortices under a fully turbulent neutral ABL inflow is due to complex interactions across a range of excitation frequencies, in which the mutual-inductance instability may not be the dominant one.

17 WIND ENERGY↗
DOE OSTI2025

CaloChallenge 2022: a community challenge for fast calorimeter simulation

Here, we present the results of the ‘Fast Calorimeter Simulation Challenge 2022’—the CaloChallenge. We study state-of-the-art generative models on four calorimeter shower datasets of increasing dimensionality, ranging from a few hundred voxels to a few tens of thousand voxels. The 31 individual submissions span a wide range of current popular generative architectures, including variational autoencoders (VAEs), generative adversarial networks (GANs), normalizing flows, diffusion models, and models based on conditional flow matching. We compare all submissions in terms of quality of generated calorimeter showers, as well as shower generation time and model size. To assess the quality we use a broad range of different metrics including differences in one-dimensional histograms of observables, KPD/FPD scores, AUCs of binary classifiers, and the log-posterior of a multiclass classifier. The results of the CaloChallenge provide the most complete and comprehensive survey of cutting-edge approaches to calorimeter fast simulation to date. In addition, our work provides a uniquely detailed perspective on the important problem of how to evaluate generative models. As such, the results presented here should be applicable for other domains that use generative AI and require fast and faithful generation of samples in a large phase space.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Glovebox CFD Analysis: Low-Flow, Laminar CFD Analysis of Glovebox Ventilation and Humidity Transport in Ansys Fluent 2025 R1

There are gloveboxes that frequently exceed the 30% relative humidity level required for operations. These exceedances result in alternative operations needing to be utilized. This report documents a computational fluid dynamics (CFD) analysis performed to characterize airflow and humidity transport within a representative glovebox under nominal, low-flow conditions. The analysis focuses on establishing the flow regime, identifying transport-limited regions, and evaluating the implications of laminar flow on humidity removal. A single-phase, laminar flow model with a species transport model was used to represent air and water vapor mixing. Results demonstrate that, at the nominal flow rates, airflow is laminar with weak mixing, which will result in long residence times, tens of minutes long or more, in low-velocity regions. These conclusions suggests that humidity reduction is transport-limited. This analysis provides qualitative insight into airflow structure and humidity transport trends. In addition to characterization of glovebox airflow and humidity transport, results are intended to inform experimental validation and future transient modeling efforts.

42 ENGINEERING↗
DOE OSTI2020

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

42 ENGINEERING↗
DOE OSTI2022

Additively manufactured multiplexed inertial coalescence filters

Multiphase flows often pose a significant challenge to the efficient and reliable design of thermofluidic systems. This paper describes multiplexed inertial coalescence filters composed of parallel helical pathways, designed to capture fine droplets (<40 µm) through inertial separation while maintaining a low pressure drop (<400 Pa). Filtration efficiencies for 7 µm and 30 µm droplets were characterized for varying flow conditions, with complete capture observed above a threshold flow rate. Models for filtration efficiency and pressure drop were developed and validated against experimental results to allow system design and optimization, enabled by the tunable additive manufacturing approach used to fabricate the filters. Filter quality factor was computed for varying droplet sizes, showcasing exceptional quality factor when compared to state-of-the-art filters documented in the literature. In conclusion, this multiplexed inertial coalescence filtration approach could find use in dehumidification systems, fog harvesting, chemical reactors, and microgravity droplet capture.

3D Printing↗
DOE OSTI2020

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗
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↗
DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗