Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid flow regimes”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Benchmark of interfacial area concentration approaches for the two-fluid model in gas-dispersed condensing flow

Accurate predictions of interfacial area concentration (IAC) are required for closure of the condensation phase change term. The interfacial area transport equation (IATE) can dynamically predict interfacial area concentration and is increasingly being used in the benchmarking of phase change flows over static based correlations which require a flow regime map. For higher flow regimes, the two-group two-fluid IATE model can separately model the transport of spherical bubbles and cap/slug bubbles. Recently, a modified formulation of the one-group and two-group two-fluid model with interfacial area transport equation was proposed for gas-dispersed condensing flows. Three datasets have been used for validation consisting of condensing two-phase flow data in the unheated section of a vertical annulus with conditions spanning a range of pressures, inlet subcoolings and mass fluxes for benchmarking the various IAC approaches. For bubbly flows, the coupled void transport - IATE simulation results in good agreement in the predictions of the void fraction and the interfacial area concentration under nearly all sixty conditions. Considering the proposed two-group two-fluid IATE model formulation, the group-2 Nusselt number correlation is investigated to account for the appropriate heat transfer length scale for group-2 bubbles. Validation of the proposed two-group two-fluid model against the new dataset shows good agreement with experimental data with significant improvement over the one-group IATE model for conditions with substantial group-2 void fraction throughout the simulation domain.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Monitoring of Temperature Measurements for Different Flow Regimes in Water and Galinstan with Long Short-Term Memory Networks and Transfer Learning of Sensors

Temperature sensing is one of the most common measurements of a nuclear reactor monitoring system. The coolant fluid flow in a reactor core depends on the reactor power state. We investigated the monitoring and estimation of the thermocouple time series using machine learning for a range of flow regimes. Measurement data were obtained, in two separate experiments, in a flow loop filled with water and with liquid metal Galinstan. We developed long short-term memory (LSTM) recurrent neural networks (RNNs) for sensor predictions by training on the sensor’s own prior history, and transfer learning LSTM (TL-LSTM) by training on a correlated sensor’s prior history. Sensor cross-correlations were identified by calculating the Pearson correlation coefficient of the time series. The accuracy of LSTM and TL-LSTM predictions of temperature was studied as a function of Reynolds number (Re). The root-mean-square error (RMSE) for the test segment of time series of each sensor was shown to linearly increase with Re for both water and Galinstan fluids. Using linear correlations, we estimated the range of values of Re for which RMSE is smaller than the thermocouple measurement uncertainty. For both water and Galinstan fluids, we showed that both LSTM and TL-LSTM provide reliable estimations of temperature for typical flow regimes in a nuclear reactor. The LSTM runtime was shown to be substantially smaller than the data acquisition rate, which allows for performing estimation and validation of sensor measurements in real time.

Pantopoulou, Stella↗
DOE OSTI2023

A unifying criterion of the boiling crisis

Abstract We reveal and justify, both theoretically and experimentally, the existence of a unifying criterion of the boiling crisis. This criterion emerges from an instability in the near-wall interactions of bubbles, which can be described as a percolation process driven by three fundamental boiling parameters: nucleation site density, average bubble footprint radius and product of average bubble growth time and detachment frequency. Our analysis suggests that the boiling crisis occurs on a well-defined critical surface in the multidimensional space of these parameters. Our experiments confirm the existence of this unifying criterion for a wide variety of boiling surface geometries and textures, two boiling regimes (pool and flow boiling) and two fluids (water and liquid nitrogen). This criterion constitutes a simple mechanistic rule to predict the boiling crisis, also providing a guiding principle for designing boiling surfaces that would maximize the nucleate boiling performance.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Comparison of Archean and Phanerozoic granulites: Southern India and North American Appalachians

Archean granulites at the southern end of the Dharwar craton of India and Phanerozoic granulites in the southern Appalachians of North America share an important characteristic: both show continuous transitions from amphibolite facies rocks to higher grade. This property is highly unusual for granulite terranes, which commonly are bounded by major shears or thrusts. These two terranes thus offer an ideal opportunity to compare petrogenetic models for deep crustal rocks formed in different time periods, which conventional wisdom suggests may have had different thermal profiles. The salient features of the Archean amphibolite-to-granulite transition in southern India have been recently summarized. The observed metamorphic progression reflects increasing temperature and pressure. Conditions for the Phanerozoic amphibolite-to-granulite transition in the southern Appalachians were documented. The following sequence of prograde reactions was observed: kyanite = sillimanite, muscovite = sillimanite + K-feldspar, partial melting of pelites, and hornblende = orthopyroxene + clinopyroxene + garnet. The mineral compositions of low-variance assemblages in mafic and intermediate rocks are almost identical for the two granulite facies assemblages. In light of their different fluid regimes and possible mechanisms for heat flow augmentation, it seems surprising that these Archean and Phanerozoic granulite terranes were apparently metamorphosed under such similar conditions of pressure and temperature. Comparison with other terrains containing continuous amphibolite-to-granulite facies transitions will be necessary before this problem can be addressed.

Mcsween, Harry Y., Jr.↗
NASA NTRS1992

Flux-split algorithms for hypersonic flows

This study reviews numerical techniques for the simulation of fluid flows spanning the range of reactive regimes, from local chemical equilibrium to full thermochemical non-equilibrium. In particular, characteristic-based algorithms are considered, of the flux-vector and the flux-difference type, which are becoming ever more popular due to their accurate rendition of physically complicated, shock-wave dominated flows. Consideration is given to the problems associated with modeling the thermo-chemical behavior of reactive mixtures of thermally perfect gases, and a few topics of current research are outlined.

Cinnella, P.↗
NASA NTRS1996

A New Class of Finite Difference Schemes

Fluid flows in the transitional and turbulent regimes possess a wide range of length and time scales. The numerical computation of these flows therefore requires numerical methods that can accurately represent the entire, or at least a significant portion, of this range of scales. The inaccurate representation of small scales is inherent to non-spectral schemes. This can be detrimental to computations where the energy in the small scales is comparable to that in the larger scales, e.g. large-eddy simulations of high Reynolds number turbulence. The inaccurate numerical representation of the small scales in these large-eddy simulations can result in the numerical error overwhelming the contribution of the subgrid-scale model.

Mahesh, K.↗
NASA NTRS2006

Notes on Earth Atmospheric Entry for Mars Sample Return Missions

The entry of sample return vehicles (SRVs) into the Earth's atmosphere is the subject of this document. The Earth entry environment for vehicles, or capsules, returning from the planet Mars is discussed along with the subjects of dynamics, aerodynamics, and heat transfer. The material presented is intended for engineers and scientists who do not have strong backgrounds in aerodynamics, aerothermodynamics and flight mechanics. The document is not intended to be comprehensive and some important topics are omitted. The topics considered in this document include basic principles of physics (fluid mechanics, dynamics and heat transfer), chemistry and engineering mechanics. These subjects include: a) fluid mechanics (aerodynamics, aerothermodynamics, compressible fluids, shock waves, boundary layers, and flow regimes from subsonic to hypervelocity; b) the Earth s atmosphere and gravity; c) thermal protection system design considerations; d) heat and mass transfer (convection, radiation, and ablation); e) flight mechanics (basic rigid body dynamics and stability); and f) flight- and ground-test requirements; and g) trajectory and flow simulation methods.

Rivell, Thomas↗
DOE OSTI2022

Electrostatic Relativistic Fluid Models of Electron Emission in a Warm Diode

A semi-analytic fluid model has been developed for characterizing relativistic electron emission across a warm diode gap. Here we demonstrate the use of this model in (i) verifying multi-fluid codes in modeling compressible relativistic electron flows (the EMPIRE-Fluid code is used as an example; see also Ref. 1), (ii) elucidating key physics mechanisms characterizing the influence of compressibility and relativistic injection speed of the electron flow, and (iii) characterizing the regimes over which a fluid model recovers physically reasonable solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Global weather research

Modeling, prediction, and analysis of global meteorological phenomena influencing the large scale behavior of the atmosphere are summarized. Prediction of global weather phenomena based on satellite data is discussed and models of global phenomena developed. The atmospheric general circulation model (AGCE) is reviewed, axisymmetric flow calculated, and axisymmetric states in cylindrical, spherical, three dimensional, and spin up numerical models for AGCE described. The role of latent heat release in baroclinic waves, latent heat and cyclonic systems, and a theoretical study of baroclinic flow related to the AGCE and the flow regime were studied with a simplified general circulation model. AGCE and the geophysical fluid flow cell (GFFC) instrumentation are discussed. Investigation of solar and planetary convection for GFFC is described. The utilization of satellite cloud observations to diagnose the energy state and transformations in extratropical cyclones is reviewed.

Source record↗
DOE OSTI2023

The influence of permeability anisotropy in the upper ocean crust on advective heat transport by a ridge-flank hydrothermal system

Here, in this study, we highlight the importance of permeability anisotropy on the hydrogeological regime of a ridge-flank hydrothermal system. Our study site, North Pond, is a marine sediment pond on ~8 Ma seafloor in the North Atlantic, and represents a low-temperature, end-member ridge-flank hydrothermal system. Previous simulations of North Pond elucidated long-standing hypotheses concerning hydrothermal fluid and heat transport in the upper volcanic crust but failed to fully explain observed patterns of seafloor heat flux in this area. Here we use variography, a geostatistical method, to quantify relations between seafloor heat-flux measurements, and coupled numerical simulations of fluid and heat flow to simulate the hydrogeologic regime. Directional variography shows that heat-flux observations are correlated along-strike of the regional crustal fabric. Three-dimensional simulations that include permeability anisotropy are able to replicate seafloor heat-flux patterns across North Pond. The simulations that result in the best match to thermal data incorporate permeability anisotropy in the horizontal plane. We find that the feedback between permeability anisotropy and the asymmetric geometry of North Pond combine to promote advective removal of heat and mass within the crustal aquifer. These findings suggest that permeability anisotropy in the oceanic crust may influence ridge-flank hydrothermal circulation more broadly.

58 GEOSCIENCES↗
NASA NTRS1986

Semi-implicit and fully implicit shock-capturing methods for hyperbolic conservation laws with stiff source terms

Some numerical aspects of finite-difference algorithms for nonlinear multidimensional hyperbolic conservation laws with stiff nonhomogenous (source) terms are discussed. If the stiffness is entirely dominated by the source term, a semi-implicit shock-capturing method is proposed provided that the Jacobian of the soruce terms possesses certain properties. The proposed semi-implicit method can be viewed as a variant of the Bussing and Murman point-implicit scheme with a more appropriate numerical dissipation for the computation of strong shock waves. However, if the stiffness is not solely dominated by the source terms, a fully implicit method would be a better choice. The situation is complicated by problems that are higher than one dimension, and the presence of stiff source terms further complicates the solution procedures for alternating direction implicit (ADI) methods. Several alternatives are discussed. The primary motivation for constructing these schemes was to address thermally and chemically nonequilibrium flows in the hypersonic regime. Due to the unique structure of the eigenvalues and eigenvectors for fluid flows of this type, the computation can be simplified, thus providing a more efficient solution procedure than one might have anticipated.

Yee, H. C.↗
NASA NTRS1987

Semi-implicit and fully implicit shock-capturing methods for hyperbolic conservation laws with stiff source terms

Some numerical aspects of finite-difference algorithms for nonlinear multidimensional hyperbolic conservation laws with stiff nonhomogeneous (source) terms are discussed. If the stiffness is entirely dominated by the source term, a semi-implicit shock-capturing method is proposed provided that the Jacobian of the source terms possesses certain properties. The proposed semi-implicit method can be viewed as a variant of the Bussing and Murman point-implicit scheme with a more appropriate numerical dissipation for the computation of strong shock waves. However, if the stiffness is not solely dominated by the source terms, a fully implicit method would be a better choice. The situation is complicated by problems that are higher than one dimension, and the presence of stiff source terms further complicates the solution procedures for alternating direction implicit (ADI) methods. Several alternatives are discussed. The primary motivation for constructing these schemes was to address thermally and chemically nonequilibrium flows in the hypersonic regime. Due to the unique structure of the eigenvalues and eigenvectors for fluid flows of this type, the computation can be simplified, thus providing a more efficient solution procedure than one might have anticipated.

Yee, H. C.↗
NASA NTRSDate not supplied

Predictions of Line Chilldown Boiling Regime Transitions by a Coupled CFD-Sub-Grid Boiling Model Validated against 1G LN2 Experiments

Propellant storage and transfer during future long-duration missions will involve fuel depot operations in which a donor depot tank is used to fill a receiver spacecraft tank with a cryogenic propellant. Prior to the on-Orbit filling operation both the receiver tank and the transfer line must be cooled. The line chilldown process involves transition between boiling regimes in microgravity that will be quite different from their 1g ground-based counterparts. Since the cryogenic propellant itself will be used to perform the chilldown process, the time constants to cool the wall and the amount of fuel that will be used become important design considerations. In this light, the focus of the present work is to capture the flow boiling characteristics of a cryogenic fluid, namely, liquid nitrogen, during the chill-down of a transfer line using CFD modelling and simulations. The cryogenic chill-down process involves different flow boiling regimes: film boiling, transition boiling and nucleate boiling. The prediction of transition between these regimes in a CFD framework is new and challenging. The present work addresses this challenge by employing a volume-of-fluid (VOF) based methodology with Lee phase change model to predict the film boiling regime of the chill-down process in ANSYS Fluent®. The transition and nucleate boiling regimes are predicted by incorporating a sub-grid model that accounts for bubble nucleation, growth, shedding frequency, and departure diameter. The sub-grid model is implemented into Fluent via a user-defined function for wall-fluid heat flux calculations. The sub-grid model is similar in formulation to the well-known Rensselaer Polytechnic Institute (RPI) boiling model. The CFD model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Predicted results shows good agreement of wall temperature, rewetting temperature, and transition between film and nucleate boiling with the experimental measurements for several different LN2 flowrates in the vertical orientation. The strengths and weaknesses of this modeling approach are presented and discussed.

Cryogenic Fluid Management↗
NASA NTRS1993

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗
NASA NTRSDate not supplied

Validation of CFD Model Prediction of Flow Boiling Regime Transitions during LN2 Line Chilldown

Propellant storage and transfer during future long-duration missions will involve fuel depot operations in which a donor depot tank is used to fill a receiver spacecraft tank with a cryogenic propellant. Prior to the on-Orbit filling operation both the receiver tank and the transfer line must be cooled. The line chilldown process involves transition between boiling regimes in microgravity that will be quite different from their 1g ground-based counterparts. Since the cryogenic propellant itself will be used to perform the chilldown process, the time constants to cool the wall and the amount of fuel that will be used become important design considerations. In this light, the focus of the present work is to understand the flow boiling characteristics of a cryogenic fluid, namely, liquid nitrogen, during the chill-down of a transfer line using CFD modelling and simulations. The cryogenic chill-down process involves different flow boiling regimes: film boiling, transition boiling and nucleate boiling. The prediction of transition between these regimes in a CFD framework is new and challenging. The present work addresses this challenge by employing a volume-of-fluid (VOF) based methodology with Lee phase change model to predict the film boiling regime of the chill-down process in ANSYS Fluent®. The transition and nucleate boiling regimes are predicted by incorporating a sub-grid model developed by Craft Tech that accounts for bubble nucleation, growth, shedding frequency, and departure diameter. The sub-grid model is implemented into Fluent via a user-defined function for wall-fluid heat flux calculations. The sub-grid model is similar in formulation to the well-known Rensselaer Polytechnic Institute (RPI) boiling model. The model constants are tested for different operating conditions and designated values are reported. The CFD model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Predicted results showing good agreement of wall temperature, rewetting temperature, and transition between film and nucleate boiling with the experimental measurements are presented and discussed.

Boiling regimes↗
NASA NTRS1976

Fundamentals of fluid sealing

The fundamentals of fluid sealing, including seal operating regimes, are discussed and the general fluid-flow equations for fluid sealing are developed. Seal performance parameters such as leakage and power loss are presented. Included in the discussion are the effects of geometry, surface deformations, rotation, and both laminar and turbulent flows. The concept of pressure balancing is presented, as are differences between liquid and gas sealing. Mechanisms of seal surface separation, fundamental friction and wear concepts applicable to seals, seal materials, and pressure-velocity (PV) criteria are discussed.

Zuk, J.↗
NASA NTRS1996

Directional Solidification of Mercury Cadmium Telluride During the Second United States Microgravity Payload Mission (USMP-2)

As a solid solution semiconductor having, a large separation between liquidus and solidus, mercury cadmium telluride (MCT) presents a formidable challenge to crystal growers desiring an alloy of high compositional uniformity. To avoid constitutional supercooling during Bridgman crystal growth it is necessary to solidify slowly in a high temperature gradient region. The necessary translation rate of less than 1 mm/hr results in a situation where fluid flow induced by gravity on earth is a significant factor in material transport. The Advanced Automated Directional Solidification Furnace (AADSF) is equipped to provide the stable thermal environment with a high gradient, and the required slow translation rate needed. Ground based experiments in AADSF show clearly the dominance of flow driven transport. The first flight of AADSF in low gravity on USMP-2 provided an opportunity to test theories of fluid flow in MCT and showed several solidification regimes which are very different from those observed on earth. Residual acceleration vectors in the orbiter during the mission were measured by the Orbital Acceleration Research Experiment (OARE), and correlated well with observed compositional differences in the samples.

Gillies, D. C.↗