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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

CFD Extraction of Heat Transfer Coefficient in Cryogenic Propellant Tanks

Current reduced-order thermal model for cryogenic propellant tanks is based on correlations built for flat plates collected in the 1950's. The use of these correlations suffers from inaccurate geometry representation; inaccurate gravity orientation; ambiguous length scale; and lack of detailed validation. This study uses first-principles based CFD methodology to compute heat transfer from the tank wall to the cryogenic fluids and extracts and correlates the equivalent heat transfer coefficient to support reduced-order thermal model. The CFD tool was first validated against available experimental data and commonly used correlations for natural convection along a vertically heated wall. Good agreements between the present prediction and experimental data have been found for flows in laminar as well turbulent regimes. The convective heat transfer between the tank wall and cryogenic propellant, and that between the tank wall and ullage gas were then simulated. The results showed that the commonly used heat transfer correlations for either vertical or horizontal plate over-predict heat transfer rate for the cryogenic tank, in some cases by as much as one order of magnitude. A characteristic length scale has been defined that can correlate all heat transfer coefficients for different fill levels into a single curve. This curve can be used for the reduced-order heat transfer model analysis.

Yang, H. Q.↗
NASA NTRS2003

A Eulerian-Lagrangian Model to Simulate Two-Phase/Particulate Flows

Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.

Apte, S. V.↗
NASA NTRS1985

Theoretical and experimental design studies for the Atmospheric General Circulation Experiment

The major criterion for the Atmospheric General Circulation Experiment (AGCE) design is that it be possible to realize strong baroclinic instability in the spherical configuration chosen. A configuration was selected in which a hemispherical shell of fluid is subjected to latitudinal temperature gradients on its spherical boundaries and the latitudinal boundaries are insulators. Work in the laboratory with a cylindrical version of this configuration revealed more instabilities than baroclinic instability. Since researchers fully expect these additional instabilities to appear in the spherical configuration also, they decided to continue the laboratory cylindrical annulus studies. Four flow regimes were identified: an axisymmetric Hadley circulation, boundary layer convection, baroclinic waves and deep thermal convection. Regime diagrams were prepared.

Fowlis, W. W.↗
NASA NTRS1987

Agreement between experimental and theoretical effects of nitrogen gas flowrate on liquid jet atomization

Two-phase flows were investigated by using high velocity nitrogen gas streams to atomize small-diameter liquid jets. Tests were conducted primarily in the acceleration-wave regime for liquid jet atomization, where it was found that the loss of droplets due to vaporization had a marked effect on drop-size measurements. In addition, four identically designed two-fluid atomizers were fabricated and tested for similarity of spray profiles. A scattered-light scanner was used to measure a characteristic drop diameter, which was correlated with nitrogen gas flowrate. The exponent of 1.33 for nitrogen gas flowrate is identical to that predicted by atomization theory for liquid jet breakup in the acceleration-wave regime. This is higher than the value of 1.2 which was previously obtained at a smapling distance of 4.4 cm downstream of the atomizer. The difference is attributed to the fact that drop-size measurements obtained at a 2.2 cm sampling distance are less affected by vaporization and dispersion of small droplets and therefore should give better agreement with atomization theory. Profiles of characteristic drop diameters were also obtained by making at least five line-of-sight measurements across the spray at several horizontal positions above and below the center line of the spray.

Ingebo, Robert D.↗
NASA NTRS1988

Space cryogenics components based on the thermomechanical (TM) effect

He II vapor-liquid phase separation (VLPS) is discussed, with emphasis on fluid-related transport phenomena. The VLPS system has been studied for both linear and nonlinear regimes, demonstrating that well-defined convection patterns exist in porous plug phase separators. In the second part, other components based on the thermomechanical effect are discussed in the limit of ideal conditions. Examples considered include the heat pipe transfer of zero net mass flow, liquid transfer pumps based on the fountain effect, mechanocaloric devices for cooling purposes, and He II vortex refrigerators.

Yuan, S. W. K.↗
NASA NTRS2013

Time-Accurate Local Time Stepping and High-Order Time CESE Methods for Multi-Dimensional Flows Using Unstructured Meshes

With the wide availability of affordable multiple-core parallel supercomputers, next generation numerical simulations of flow physics are being focused on unsteady computations for problems involving multiple time scales and multiple physics. These simulations require higher solution accuracy than most algorithms and computational fluid dynamics codes currently available. This paper focuses on the developmental effort for high-fidelity multi-dimensional, unstructured-mesh flow solvers using the space-time conservation element, solution element (CESE) framework. Two approaches have been investigated in this research in order to provide high-accuracy, cross-cutting numerical simulations for a variety of flow regimes: 1) time-accurate local time stepping and 2) highorder CESE method. The first approach utilizes consistent numerical formulations in the space-time flux integration to preserve temporal conservation across the cells with different marching time steps. Such approach relieves the stringent time step constraint associated with the smallest time step in the computational domain while preserving temporal accuracy for all the cells. For flows involving multiple scales, both numerical accuracy and efficiency can be significantly enhanced. The second approach extends the current CESE solver to higher-order accuracy. Unlike other existing explicit high-order methods for unstructured meshes, the CESE framework maintains a CFL condition of one for arbitrarily high-order formulations while retaining the same compact stencil as its second-order counterpart. For large-scale unsteady computations, this feature substantially enhances numerical efficiency. Numerical formulations and validations using benchmark problems are discussed in this paper along with realistic examples.

Chang, Chau-Lyan↗
NASA NTRSDate not supplied

Experimental Studies on the Effect of Inlet Liquid Subcooling on the Chilldown Characteristics of Stainless Steel Tubes

Understanding transfer line chilldown process under microgravity is important for the efficient transfer of cryogenic propellants in space fuel depots to facilitate future long duration space missions. The present work is part of the ongoing efforts to develop and test a two-phase flow chilldown test section to study the complete chilldown process under sustained microgravity conditions onboard the International Space Station. In this study, ground-based chilldown experiments are carried out on a 60 cm long SS-316 test section with PF-5060 as the working fluid. The complete chilldown curve was obtained including the film, transition and nucleate boiling regimes along with the temperature transition points. The effect of inlet liquid subcooling on the behavior of the chilldown curves are presented. Further, the chilldown and heat flux curves are analyzed to obtain re-wetting/Leidenfrost and onset of nucleate boiling temperature transition points as well as the critical heat flux values. The effect of inlet liquid subcooling on regime-specific heat flux and heat transfer coefficients are also examined.

CFD↗
NASA NTRS1977

Dissipative, forced turbulence in two-dimensional magnetohydrodynamics

The equations of motion of a turbulent two-dimensional MHD flow are solved in the presence of finite viscosity and resistivity for the case when external mechanical and/or magnetic forces act on the fluid, the goal of the study being to verify the existence of a MHD dynamo effect which is represented by a substantial back-transfer of mean square vector potential to the longest allowed Fourier wavelengths. The regime explored is that for which the mechanical and magnetic Reynolds numbers are in the range 100 to 1000. It is concluded that mechanical forcing terms alone cannot lead to dynamo action, but that dynamo action can result from either magnetic forcing terms or from both mechanical and magnetic forcing terms simultaneously.

Fyfe, D.↗
NASA NTRS1989

Solution of the hydrodynamic device model using high-order non-oscillatory shock capturing algorithms

A micron n+ - n - n+ silicon diode is simulated via the hydrodynamic model for carrier transport. The numerical algorithms employed are for the non-steady case, and a limiting process is used to reach steady state. The simulation employs shock capturing algorithms, and indeed shocks, or very rapid transition regimes, are observed in the transient case for the coupled system, consisting of the potential equation and the conservation equations describing charge, momentum, and energy transfer for the electron carriers. These algorithms, termed essentially non-oscillatory, were successfully applied in other contexts to model the flow in gas dynamics, magnetohydrodynamics, and other physical situations involving the conservation laws in fluid mechanics. The method here is first order in time, but the use of small time steps allows for good accuracy. Runge-Kutta methods allow one to achieve higher accuracy in time if desired. The spatial accuracy is of high order in regions of smoothness.

Fatemi, Emad↗
NASA NTRS2000

The Dynamics of Miscible Interfaces: A Space Flight Experiment

Experiments as well as accompanying simulations are described that serve in preparation of a space flight experiment to study the dynamics of miscible interfaces. The investigation specifically addresses the importance of both nonsolenoidal effects as well as nonconventional Korteweg stresses in flows that give rise to steep but finite concentration gradients. The investigation focuses on the flow in which a less viscous fluid displaces one of higher viscosity and different density within a narrow capillary tube. The fluids are miscible in all proportions. An intruding finger forms that occupies a fraction of the total tube diameter. Depending on the flow conditions, as expressed by the Peclet number, a dimensionless viscosity ratio, and a gravity parameter, this fraction can vary between approximately 0.9 and 0.2. For large Pe values, a quasi-steady finger forms, which persists for a time of O(Pe) before it starts to decay, and Poiseuille flow and Taylor dispersion are approached asymptotically. Depending on the specific flow conditions, we observe a variety of topologically different streamline patterns, among them some that leak fluid from the finger tip. For small Pe values, the flow decays from the start and asymptotically reaches Taylor dispersion after a time of O(Pe). Comparisons between experiments and numerical simulations based on the 'conventional' assumption of solenoidal velocity fields and without Korteweg stresses yield poor agreement as far as the Pe value is concerned that distinguishes these two regimes. As one possibility, we attribute this lack of agreement to the disregard of these terms. An attempt is made to use scaling arguments in order to evaluate the importance of the Korteweg stresses and of the assumption of solenoidality. While these effects should be strongest in absolute terms when steep concentration fronts exist, i.e., at large Pe, they may be relatively most important at lower values of Pe. We subsequently compare these conventional simulations to more complete simulations that account for nonvanishing divergence as well as Korteweg stresses. While the exact value of the relevant stress coefficients are not known, ballpark numbers do exist, and their use in the simulations indicates that these stresses may indeed be important. We plan to evaluate these issues in detail by means of comparing a space experiment with corresponding simulations, in order to extract more accurate Korteweg stress coefficients, and to confirm or deny the importance of such stresses.

Maxworthy, Tony↗
NASA NTRS2011

Time Accurate CFD Simulations of the Orion Launch Abort Vehicle in the Transonic Regime

Significant asymmetries in the fluid dynamics were calculated for some cases in the CFD simulations of the Orion Launch Abort Vehicle through its abort trajectories. The CFD simulations were performed steady state and in three dimensions with symmetric geometries, no freestream sideslip angle, and motors firing. The trajectory points at issue were in the transonic regime, at 0 and +/- 5 angles of attack with the Abort Motors with and without the Attitude Control Motors (ACM) firing. In some of the cases the asymmetric fluid dynamics resulted in aerodynamic side forces that were large enough that would overcome the control authority of the ACMs. MSFC's Fluid Dynamics Group supported the investigation into the cause of the flow asymmetries with time accurate CFD simulations, utilizing a hybrid RANS-LES turbulence model. The results show that the flow over the vehicle and the subsequent interaction with the AB and ACM motor plumes were unsteady. The resulting instantaneous aerodynamic forces were oscillatory with fairly large magnitudes. Time averaged aerodynamic forces were essentially symmetric.

Rojahn, Josh↗
NASA NTRS2011

Time Accurate CFD Simulations of the Orion Launch Abort Vehicle in the Transonic Regime

Significant asymmetries in the fluid dynamics were calculated for some cases in the CFD simulations of the Orion Launch Abort Vehicle through its abort trajectories. The CFD simulations were performed steady state with symmetric boundary conditions and geometries. The trajectory points at issue were in the transonic regime, at 0 and 5 angles of attack with the Abort Motors with and without the Attitude Control Motors (ACM) firing. In some of the cases the asymmetric fluid dynamics resulted in aerodynamic side forces that were large enough that would overcome the control authority of the ACMs. MSFC s Fluid Dynamics Group supported the investigation into the cause of the flow asymmetries with time accurate CFD simulations, utilizing a hybrid RANS-LES turbulence model. The results show that the flow over the vehicle and the subsequent interaction with the AB and ACM motor plumes were unsteady. The resulting instantaneous aerodynamic forces were oscillatory with fairly large magnitudes. Time averaged aerodynamic forces were essentially symmetric.

Ruf, Joseph↗
DOE OSTI2020

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗
NASA NTRSDate not supplied

Hydrodynamics of Two-Phase Flows through Porous Media in Microgravity: Packed Bed Reactor Experiment onboard of the International Space Station

The objective of the Packed Bed Reactor Experiment was to develop a fundamental understanding of the hydrodynamics of two-phase flow with no phase transition through porous media in microgravity. This work describes two experimental campaigns flown on the International Space Station (ISS) that were designed to achieve the objective. This work presents the results on flow patterns observed, two-phase flow pressure drop in porous media, and the impact of capillary effects on pressure drop at various gas and liquid flow rates. This work is the first to present predictive correlations of the two-phase friction factor for the different regimes identified based on the superficial liquid and gas velocity in microgravity, the first to hypothesize the different regimes based on the change in slope of the pressure gradient versus flow rate plots, the first to address the effects of the capillarity on the pressure gradient and the first to measure column holdup and assess the effective porosity based on the packing and the trapped gas bubbles. Experimental data on pressure drop in gas (N2) – liquid (water) flow show the presence of four different flow regimes in microgravity. Dispersed bubble flow (DB) and pulse (P) flow regimes are detected at high liquid flow rates, whereas at low liquid flow rates “large bubble” or elongated bubble (LB) and “gas channeling” (GC) regimes are observed. For these different flow regimes, different two-phase flow friction factor and pressure gradient correlations are presented for the first time as a function of the gas and liquid modified Reynolds numbers and Suratman number. Within the viscous-capillary (V-C) regime, it is found that the capillary contribution is the dominant force that contributes to the pressure drop for the wetting case (glass). However, for the non-wetting packing (Teflon), the viscous contribution dominates. It was found that the gas hold-up and pressure drop are functions of bed history at low liquid and gas flow rates with the magnitude of the hysteresis decreasing with increasing flow rates. PBRE-2 results show that the capillary force is a strong function of the superficial liquid velocity but is a much weaker function of the superficial gas velocity and varies inversely with the particle diameter. Within the Viscous-Capillary (V-C) regime, over 90% of the pressure gradient is attributed to the capillary contribution in the gas continuous regime. However, in the large bubble regime, the viscous and capillary contributions were comparable. After the completion of PBRE-2 with glass spherical beads, another bed packed with alumina was installed. Pressure gradient data for the alumina packed bed (PBRE-Water Recovery) were transferred to the PI and were not analyzed by the authors of this work. The Packed Bed Reactor Experiment concluded its on-orbit operation after two successful campaigns of testing in 2017 and in 2021. The flight hardware was brought back from orbit and is being reconditioned for a series of future experiments referred to as PBRE-WRS in support of water recovery (WRS). The WRS series consists of testing different two-phase fluid system components and packed beds to assess their pressure gradient characteristics, which will be used for designing packed bed reactors for various applications relevant to life support. Although there is so much relevance of the experimental results obtained from the two PBRE campaign in microgravity to life support, these findings also apply to other applications that involve two phase flow in porous medium such as fuel cells, transport of nutrient to plants in space and other chemical and materials processing that involve two phase flows. These systems operate differently in microgravity because, due to the lack of buoyancy, the density difference between the phases becomes irrelevant and no longer leads to phase separation.

Packed Bed↗
DOE OSTI2021

Preconditioners for multiphase poromechanics with strong capillarity

This paper aims to enhance the performance of Newton–Krylov solvers for coupled poromechanical problems with two-phase flow. In particular, we investigate the impact of capillary pressure on preconditioning strategies. Capillarity complicates the coupling between the solid deformation and fluid pressure degrees of freedom, as well as increases the nonlinearity of the system. Depending on the capillary pressure relation used in the constitutive formulation, the flow equations may exhibit a spectrum of advection-dominated to diffusion-dominated behavior. We propose preconditioning approaches that account for this behavior and lead to robust numerical performance within a broad range of regimes.

42 ENGINEERING↗
NASA NTRS1988

Experimental study of thermocapillary flows in a thin liquid layer with heat fluxes imposed on the free surface

To study thermocapillary flows in a two-dimensional thin liquid layer with heat fluxes imposed on the free surface experimentally, a long tray configuration was employed to simulate the infinite layer. The surface temperature distribution due to thermocapillary convective for different flow regimes was measured and compared with theoretical predictions. A short tray configuration was also employed to study the end wall effects (insulating or conducting). The results show that for a strong convection flow with an insulating wall as the boundary the surface temperature distribution became quite uniform. Consequently, the thermocapillary driving force was greatly reduced. On the other hand, a strong fluid motion always existed adjacent to the conducting wall because of the large surface temperature gradient near the wall.

Lai, Chun-Liang↗
NASA NTRS1988

Experimental study of thermocapillary flows in a thin liquid layer with heat fluxes imposed on the free surface

To study thermocapillary flows in a two-dimensional thin liquid layer with heat fluxes imposed on the free surface experimentally, a long tray configuration was employed to simulate the infinite layer. The surface temperature distribution due to thermocapillary convection for different flow regimes was measured and compared with theorectical predictions. A short tray configuration was also employed to study the end wall effects (insulating or conducting). The results show that, for a strong convection flow with an insulating wall as the boundary, the surface temperature distribution became quite uniform. Consequently, the thermocapillary driving force was greatly reduced. On the other hand, a strong fluid motion always existed adjacent to the conducting wall because of the large surface temperature gradient near the wall.

Lai, Chun-Liang↗
NASA NTRS1991

A computer code for multiphase all-speed transient flows in complex geometries. MAST version 1.0

The operation of the MAST code, which computes transient solutions to the multiphase flow equations applicable to all-speed flows, is described. Two-phase flows are formulated based on the Eulerian-Lagrange scheme in which the continuous phase is described by the Navier-Stokes equation (or Reynolds equations for turbulent flows). Dispersed phase is formulated by a Lagrangian tracking scheme. The numerical solution algorithms utilized for fluid flows is a newly developed pressure-implicit algorithm based on the operator-splitting technique in generalized nonorthogonal coordinates. This operator split allows separate operation on each of the variable fields to handle pressure-velocity coupling. The obtained pressure correction equation has the hyperbolic nature and is effective for Mach numbers ranging from the incompressible limit to supersonic flow regimes. The present code adopts a nonstaggered grid arrangement; thus, the velocity components and other dependent variables are collocated at the same grid. A sequence of benchmark-quality problems, including incompressible, subsonic, transonic, supersonic, gas-droplet two-phase flows, as well as spray-combustion problems, were performed to demonstrate the robustness and accuracy of the present code.

Chen, C. P.↗