Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow conservation”

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 OSTI2025

High y + Shear-Stress Turbulence Implementation for High Flux Isotope Reactor Narrow Channel Flows

The research objective of this work was to improve the engineering predictions of the turbulence characteristics of flows in curved narrow channels. Such channel flows are commonly encountered in nuclear research and test reactors, with one of them being the high-flux isotope reactor (HFIR). Research reactors bear high heat fluxes, and the proper computing of turbulence is paramount for safe and reliable reactor operation. The study builds on the results of a previous direct numerical simulation of turbulence to inform a well-known Reynolds-averaged Navier–Stokes shear-stress turbulence model and improves its accuracy in simulating parallel channel flows. A new formulation of the loss term in the dissipation conservation equation is suggested. Combined with high wall distance computational grids, the new implementation provides a fast-running flow solution, suitable for engineering purposes. Model generalization for parallel channel flows, in a broader range of frictional Reynolds numbers, is suggested by introducing a new form of the model constants.

CFD↗
DOE OSTI2023

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗
DOE OSTI2020

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Spatial Replication Is Important for Developing Landscape Genetic Inferences for a Wetland Salamander

Habitat fragmentation is a pressing threat to wildlife populations, and maintenance of gene flow between populations is an essential goal of conservation. Resistance surfaces have emerged as an important tool for modelling connectivity and developing management strategies to mitigate effects of habitat fragmentation. However, recent studies have noted inconsistencies in the factors most strongly associated with connectivity across different landscapes. Thus, replication of genetic-based resistance surface optimisation across landscapes may be necessary for making robust conclusions about the influence of environmental variables. Accordingly, replication represents a substantive challenge and opportunity in the field of landscape genetics. In this study, we conducted replicated landscape genetic analyses across five landscapes in Tennessee and Kentucky for a threatened wetland amphibian, the four-toed salamander (Hemidactylium scutatum). We tested multiple hypotheses of how different landscape features that could directly affect small, desiccation-intolerant amphibians (e.g., canopy cover) influenced gene flow and assessed the appropriate scale at which to model different features. We found some concordance in the landscape features that influenced gene flow (e.g., a common importance of forest cover and topography), but also some differences—potentially owing to the difference in variability of predictors across landscapes. We also found discordance in the scale of effect of different features across landscapes. In conclusion, our work emphasises that flat areas of moist forest not bisected by roads may be important for H. scutatum conservation, and our replicated design allows us to identify relationships that would have been missed if only using one study site.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Abundance and Migration Success of Overshoot Steelhead in the Upper Columbia River

Abstract Summer steelhead Oncorhynchus mykiss may enter freshwater almost a year before spawning and potentially make long migrations (>1,000 km) to interior headwater habitats. However, in response to suboptimal freshwater habitat conditions (e.g., warmer water temperatures), adult summer steelhead may exhibit complex behaviors during upstream migration in the Columbia River basin. Steelhead may migrate upstream of their natal tributary (hereafter, referred to as “overshoot”) and spend days to several months before subsequently migrating downstream (hereafter, referred to as “fallback”) to their natal tributary to spawn. An expansion of an existing Bayesian patch occupancy model, derived from observations of adult steelhead that were PIT-tagged to estimate population-specific abundance upstream of the tagging location, incorporated downstream detection locations to estimate the abundance of overshoot fallbacks. Overshoot steelhead abundance at the tagging location was estimated based on the relationship between the number of known overshoot fallbacks (i.e., the number of steelhead that overshot and successfully migrated downstream to their natal tributary) and their model-estimated abundance. During the study period (2010–2017), the annual mean proportion of overshoot steelhead that successfully migrated downstream of the tagging location (Priest Rapids Dam) was 0.59 (SD = 0.14). The number of dams encountered by overshoot steelhead during their downstream migration was negatively correlated with their downstream migration success probability. Improved downstream passage survival for adult steelhead will increase the abundance of affected populations while reducing potential genetic introgression of upstream populations (i.e., strays). This is the first study to estimate the abundance of overshoot and fallback steelhead, providing the data necessary for scientists to estimate potential conservation benefits of improved downstream survival. For example, surface flow passage routes (e.g., sluiceways and temporary spillway weirs) are very effective in guiding and passing adult steelhead downstream of Columbia River hydroelectric projects and data from this assessment show that changes in dam operations throughout the downstream migration period may maximize conservation benefits.

Murdoch, Andrew R. (ORCID:0000000324827689)↗
DOE OSTI2023

Operation Optimization using Reinforcement Learning with Integrated Artificial Reasoning Framework

In large and complex systems, operational decision-making requires a systematic analysis with a vast amount of data from both process parameters and component status monitoring. In this paper, we present an integrated artificial reasoning approach for system state transition models that can help operational decision-making with explainable and traceable reasoning. The integrated artificial reasoning framework is a physics-based approach of defining the system structure in a Bayesian network, so we leveraged it in a Markov decision process (MDP) for finding optimal operational solutions. In our proposed framework, the MDP is implemented on a dynamic Bayesian network (DBN), which represents causalities in a system. The multilevel flow modeling was utilized in order to extract these causalities in a more efficient and objective manner. Since multilevel flow modeling is based on the fundamental energy and mass conservation laws, the target system is decomposed into several mass, energy, and information structures, which serve as the basis for a DBN. The MDP consists of the processes of finding a solution for the Bellman equation, which can be derived from the conditional probability equations of the constructed DBN. System operators can capture stochastic system dynamics as multiple subsystem state transitions based on their physical relations and uncertainties coming from the component degradation process or random failures. We analyzed a simplified example system to illustrate finding an optimal operational policy with this approach.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

A mathematical model of asynchronous data flow in parallel computers *

Abstract We present a simplified model of data flow on processors in a high-performance computing framework involving computations necessitating inter-processor communications. From this ordinary differential model, we take its asymptotic limit, resulting in a model which treats the computer as a continuum of processors and data flow as an Eulerian fluid governed by a conservation law. We derive a Hamilton–Jacobi equation associated with this conservation law for which the existence and uniqueness of solutions can be proven. We then present the results of numerical experiments for both discrete and continuum models; these show a qualitative agreement between the two and the effect of variations in the computing environment’s processing capabilities on the progress of the modelled computation.

Barnard, Richard C.↗
DOE OSTI2024

Development of a Hybrid Single/Two-Phase Capillary-Based Micro-Cooler using Copper Inverse Opals Wick with Silicon 3D Manifold for High-Heat-Flux Cooling Application

Previously, we reported two-phase capillary-based cooling using narrow (200 to 1000 ..mu..m) heater bridge copper inverse opal (CIO) wicks with heat flux level of 1400 Wcm -2 and superheat ~ 10 degrees C. Here, we demonstrate the area scaling of the proposed technology to a large-area micro-cooler for high-heat-flux cooling of microprocessors and power electronics. We developed a hybrid single/two-phase micro-cooler that relies on capillary-wicking in a 25-..mu..m-thick CIO with an open silicon microchannel 3D-manifold for liquid delivery and vapor extraction, achieving a high heat flux ~ 400 Wcm -2 over a heated area of 10 x 10 mm 2 . For a range of inlet water flowrates from 5 to 60 g(min) -1 , we achieved total thermal resistances and vapor qualities of 0.68 cm 2 degrees CW -1 to 0.2 cm 2 degrees CW -1 and 0.55 to 0.12, respectively. The flowrates are 10x smaller than those of conventional single- or two-phase microchannel cooling technology. The corresponding two-phase thermal resistances ranges from 0.05 to 0.02 cm 2 degrees CW -1 with temperature superheat of 8 to 6 degrees C, respectively. While the overall performance of the large-area (10 x 10 mm 2 ) capillary-based micro-cooler degraded compared to the previous demonstration of the technology for a heated area of 5 x 5 mm 2 , however, preliminary CFD modeling indicates that an improved manifold design will be able to achieve comparable performance.

capillary flow↗
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 OSTI2020

Using Differential Forms to Find Symmetries in the Noh Problem for an Ideal Gas in a Spherical System

The Euler equations describe inviscid, compressible gas flow and are derived using the laws of conservation of mass, momentum, and energy. From these equations a common test problem, the Noh problem, can be solved for invariant solutions in planar, cylindrical, and spherical geometries. Using Lie group analysis methods, we will solve the Noh Problem for an ideal gas in a spherical system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Development of a Hybrid Single/Two-Phase Capillary-Based Micro-Cooler Using Copper Inverse Opals Wick with Silicon 3D Manifold for High-Heat-Flux Cooling Application: Preprint

Previously, two-phase capillary-based cooling from narrow (200-1,000 micrometer) heater-bridge copper inverse opal (CIO) wicks with heat flux levels exceeding 1,400 watts per square cm with a low superheat of approximately 10 degrees C was demonstrated. Here, we demonstrate the area scaling of the proposed technology to large-area micro-cooler for the high-heat-flux cooling of microprocessors and power electronics. We developed a hybrid single/two-phase micro-cooler that relies on capillary wicking in 25-micrometer-thick CIOs with an open channel silicon 3D-manifold for liquid delivery and vapor extraction, to achieve a high heat flux of approximately 400 watts per square cm over a heated area of 1 square cm. For the range of inlet water (21 degrees C water temperature) flowrates from 5 to 60 mL per min, we achieved total thermal resistances and vapor qualities of 0.68-0.2 square cm Kelvin per watt and 0.55-0.12, respectively. The high heat flux levels are achieved with flowrates that are 10 times smaller than conventional single- or two-phase microchannel cooling technology. The corresponding two-phase thermal resistances are in the range of 0.05 to 0.02 square cm Kelvin per watt with temperature superheat of 8 to 6 degrees C. While the overall performance of the large-area (10 by 10 square millimeter) capillary-based micro-cooler degraded compared to previous demonstration of the technology for a heated area of 5 by 5 square millimeter, preliminary computational fluid dynamics (CFD) modeling indicates that an improved manifold design will be able to achieve comparable performance.

capillary flow↗
DOE OSTI2025

Forced flow transient safety analysis of irradiation device with adjustable orifice for research reactor fuel assemblies

The Belgium Reactor 2 (BR2) of the Belgian Nuclear Research Centre (SCK CEN) has several irradiation devices or rigs that are dedicated to the fuel performance and qualification demonstration testing of research reactor fuels. In support of the U.S. High Performance Research Reactor (USHPRR) LEU conversion project, a new flexible irradiation apparatus, MUSTANG-R, has been constructed. SCK CEN has completed the design and safety study, in cooperation with Idaho National Laboratory (INL) and Argonne National Laboratory (ANL), to allow for the irradiation testing of a full-size fuel assembly in a 200 mm diameter channel in the BR2 reactor. The moveable valve is a key design feature of the device and acts like an adjustable orifice enhancing or restricting the flow through a coolant channel inlet located in the BR2 upper plenum. This moveable valve allows the flow through the device to be adjusted prior to each BR2 cycle to obtain the necessary conditions for the fuel qualification test. This ensures accurate and representative thermal-hydraulic conditions of the fuel design are achieved. The device was designed and qualified as passively safe, implying verification by a combination of mechanical and thermal-hydraulic analysis and testing. This includes characterization of the safety margin required for a scenario where the moveable valve is assumed to be erroneously closed during irradiation. A simplified and conservative method is proposed for analyzing the corresponding forced flow transient using a critical heat flux criterion. In conclusion, this allows the required minimum valve opening to be determined for the experiments' design and safety studies.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

An immersed interface method for the 2D vorticity-velocity Navier-Stokes equations with multiple bodies

We present an immersed interface method for the vorticity-velocity form of the 2D Navier Stokes equations that directly addresses challenges posed by nonconvex immersed bodies, multiply connected domains, and the calculation of force distributions on immersed surfaces. The immersed interface method is re-interpreted as a polynomial extrapolation of flow quantities and boundary conditions into the immersed solid bodies, reducing computational cost and enabling simulations with nonconvex bodies that could not be discretized with previous immersed interface methods. In the flow, the vorticity transport equation is discretized using a conservative finite difference scheme and explicit Runge-Kutta time integration. The velocity reconstruction problem is transformed to a scalar Poisson equation that is discretized with conservative finite differences, and solved using an FFT-accelerated iterative algorithm. The use of conservative differencing throughout leads to exact enforcement of a discrete Kelvin's theorem, allowing for simulations with multiply connected domains and outflow boundaries that have challenged other immersed interface vortex methods. We also explore novel methods for recovering time-dependent pressure distributions on immersed bodies within a vorticity-based method and present a novel control volume formulation for recovering aerodynamic moments from only the vorticity and velocity fields. The method achieves second order spatial accuracy and third order temporal accuracy, and is validated on a variety of 2D flows in internal and free-space domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Conservative Finite-Volume Based Interface-Tracking Algorithm Using the Signed Distance Function

Methods for tracking an interface between two fluid phases are developed to ensure desired fluid properties, conservation, and stability are preserved in a finitevolume (FV) discretization. Here, the interface is tracked using a level set method where the signed distance function implicitly defines the interface. Marching methods are used to evaluate the value of the signed distance function, including a novel initialization method to initialize any implicit function to the signed distance function around sharp corners in the level set. Global conservation and consistency with a set of governing equations is enforced by a compression coefficient that measures the volumetric compression or expansion due to inaccuracies in the level set evaluation. A redistribution method is integrated into the volume correction to eliminate the small-cell instability while maintaining global conservation. This suite of methods is implemented and tested using static uniform velocity, and potential flow cases with multiple interface geometries. Results show these methods achieve up to second order accuracy, and are conservative. The application for these methods is intended to track the interface of a 3D printing filament in a finite-volume discretization of the all-speed Navier-Stokes equations.

42 ENGINEERING↗
DOE OSTI2023

Advanced two-phase subchannel method via non-linear iteration

A fast-running, robust two-phase flow, sub-channel model is presented based on non-linear solution of the steady-state subchannel fluid flow equations. The drift-flux model solves for conservation of liquid and vapor mass, mixture energy, and axial and transverse mixture momentum as part of an efficient planar marching scheme and nonlinear, nested outer and inner iteration. Here, models based on mechanistic subcooled boiling, two-phase turbulent void mixing, and drift are included. Solution verification and mesh convergence studies were performed for modern GE 10 × 10 fuel geometry and are shown to have excellent convergence behavior. Run time performance for a 50 axial mesh model showed 2.2 seconds on a single CPU core to tightly converge all 3D distributions (flow, void, pressure) for the GE 10 × 10 fuel geometry, supporting its efficient use within the Virtual Environment for Reactor Applications boiling water reactor framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Scale interactions and anisotropy in Rayleigh–Taylor turbulence

We study energy scale transfer in Rayleigh–Taylor (RT) flows by coarse graining in physical space without Fourier transforms, allowing scale analysis along the vertical direction. Two processes are responsible for kinetic energy flux across scales: baropycnal work Λ, due to large-scale pressure gradients acting on small scales of density and velocity; and deformation work Π, due to multiscale velocity. Our coarse-graining analysis shows how these fluxes exhibit self-similar evolution that is quadratic-in-time, similar to the RT mixing layer. We find that Λ is a conduit for potential energy, transferring energy non-locally from the largest scales to smaller scales in the inertial range where Π takes over. In three dimensions, Π continues a persistent cascade to smaller scales, whereas in two dimensions Π rechannels the energy back to larger scales despite the lack of vorticity conservation in two-dimensional (2-D) variable density flows. This gives rise to a positive feedback loop in 2-D RT (absent in three dimensions) in which mixing layer growth and the associated potential energy release are enhanced relative to 3-D RT, explaining the oft-observed larger α values in 2-D simulations. Despite higher bulk kinetic energy levels in two dimensions, small inertial scales are weaker than in three dimensions. Moreover, the net upscale cascade in two dimensions tends to isotropize the large-scale flow, in stark contrast to three dimensions. Furthermore, our findings indicate the absence of net upscale energy transfer in three-dimensional RT as is often claimed; growth of large-scale bubbles and spikes is not due to ‘mergers’ but solely due to baropycnal work Λ.

42 ENGINEERING↗
DOE OSTI2022

Hydrogen Plus Other Alternative Fuels Risk Assessment Models (HyRAM+) Version 4.1 Technical Reference Manual

The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, methane, and propane systems. HyRAM+ is designed to facilitate the use of state-of-the-art models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impacts on people. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards, and to provide developed models in a publicly-accessible toolkit usable by all stakeholders. This document provides a description of the methodology and models contained in HyRAM+ version 4.1. The two most significant changes for HyRAM+ version 4.1 from HyRAM+ version 4.0 are direct incorporation of unconfined overpressure into the QRA calculations and modification of the models for cryogenic liquid flow through an orifice. In QRA mode, the user no longer needs to input peak overpressure and impulse values that were calculated separately; rather, the unconfined overpressure is estimated for the given system inputs, leak size, and occupant location. The orifice flow model now solves for the maximum mass flux through the orifice at constant entropy while conserving energy, which does not require a direct speed of sound calculation. This does not affect the mass flow for all-gaseous releases; the method results in the same speed of sound for choked flow. However, this method does result in a higher (and more realistic) mass flow rate for a given leak size for liquid releases than was previously calculated.

08 HYDROGEN↗