Unsteady Isothermal Flow Through a Staggered Tube Bundle Array

Author(s):  
Artit Ridluan ◽  
Surasing Arayangkun ◽  
Coochart Phayom

Two-dimensional Unsteady simulations of isothermal flow through a staggered tube bundle array at three different Reynolds numbers 54, 72, and 90 were investigated. The Navier-Stokes equations are numerically solved. Based on the CFD simulation results, the unsteady flow patterns were developed behind the rear row of the array, while for the other rows, the steady separated and reattached flow behaviors were observed, small, short, and closed separation bubble behind the rods. At Reynolds Number of 54, the transient flow was perfectly periodic. The complicated patterns of unsteady flow could be observed at Reynolds numbers of 72 and 90. The shedding patterns of vortices from the last rods were different and no longer periodic as found at Reynolds number of 54. The degree of chaos is increased as Reynolds number progressed.

1990 ◽  
Vol 220 ◽  
pp. 459-484 ◽  
Author(s):  
H. M. Badr ◽  
M. Coutanceau ◽  
S. C. R. Dennis ◽  
C. Ménard

The unsteady flow past a circular cylinder which starts translating and rotating impulsively from rest in a viscous fluid is investigated both theoretically and experimentally in the Reynolds number range 103 [les ] R [les ] 104 and for rotational to translational surface speed ratios between 0.5 and 3. The theoretical study is based on numerical solutions of the two-dimensional unsteady Navier–Stokes equations while the experimental investigation is based on visualization of the flow using very fine suspended particles. The object of the study is to examine the effect of increase of rotation on the flow structure. There is excellent agreement between the numerical and experimental results for all speed ratios considered, except in the case of the highest rotation rate. Here three-dimensional effects become more pronounced in the experiments and the laminar flow breaks down, while the calculated flow starts to approach a steady state. For lower rotation rates a periodic structure of vortex evolution and shedding develops in the calculations which is repeated exactly as time advances. Another feature of the calculations is the discrepancy in the lift and drag forces at high Reynolds numbers resulting from solving the boundary-layer limit of the equations of motion rather than the full Navier–Stokes equations. Typical results are given for selected values of the Reynolds number and rotation rate.


2013 ◽  
Author(s):  
Henri J. L. van der Heiden ◽  
Peter van der Plas ◽  
Arthur E. P. Veldman ◽  
Roel W. C. P. Verstappen ◽  
Roel Luppes

In offshore applications, details of viscous flow effects can become relevant when predicting e.g. drag forces on the columns of oil drilling rigs, or the flow around a semisubmersible in figure 1. This motivates a novel approach for efficiently simulating viscous flow effects at high Reynolds numbers with the CFD simulation tool ComFLOW. In ComFLOW, the Navier–Stokes equations can be solved for one-phase and for two-phase flow. The equations are discretized second-order in space, and second-order in time. An Improved Volume-of-Fluid (IVOF) algorithm is used for free-surface advection and reconstruction [1, 2]. Modeling viscous flow effects in high Reynolds number flows requires a turbulence model that provides accurate results on coarse grids. We pursue to achieve a high local grid resolution in a computationally efficient manner. Both approaches are tested for flows around a square cylinder: grid refinement at Reynolds numbers 10 and 100, and the turbulence model at Reynolds number 22,000.


Author(s):  
Justin R. Finn ◽  
Sourabh V. Apte ◽  
Brian D. Wood

The characteristics of pore scale vortical structures observed in moderate Reynolds number flow through mono-disperse packed beds of spheres are examined. Our results come from direct numerical simulations of flow through (i) a periodic, simple cubic arrangement of 54 spheres, (ii) a wall bounded, close packed arrangement of 216 spheres, and (iii) a realistic randomly packed tube containing 326 spheres with a tube diameter to sphere diameter ratio of 5.96. Pore Reynolds numbers in the steady inertial (10 ≲ Re ≲ 200) and unsteady inertial (Re ≈ 600) regimes are considered. Even at similar Reynolds numbers, the vortical structures observed in flows through these three packings are remarkably different. The interior of the arranged packings are dominated by multi-lobed vortex ring structures which align with the principal axes of the packing. The random packing and the near wall region of the close packed arrangement are dominated by helical vortices, elongated in the mean flow direction. In the simple cubic packing, unsteady flow is marked by periodic vortex shedding which occurs at a single frequency. Conversely, at a similar Reynolds number, the vortical structures in unsteady flow through the random packing oscillate with many characteristic frequencies.


Author(s):  
Patrick H. Oosthuizen

Flow through a minichannel with a square cross-sectional shape whose center-line follows a wavy path has been numerically studied. The wavy center-line path is defined by a series of circular arcs. Now when a fluid flows through such a wavy channel unsteady flow can develop at Reynolds numbers that are far below the value at which transition to turbulence occurs. The development of this unsteadiness can lead to increases in the pressure losses in the channel and to significant increases in the heat transfer rate when the channel wall is heated or cooled. To study the conditions under which unsteady flow develops with a wavy, circular arc center-line shape attention has been given to flow through a channel system in which there is a short straight entrance section followed by a wavy channel section with four full waves followed by a short straight outlet section. It has been assumed that the flow is incompressible and that there is no slip at the channel walls. The unsteady form of the governing equations have been written in dimensionless form and solved using a commercial finite-element based software package, FIDAP. The solution has the following parameters: the Reynolds number, the Prandtl number, and the dimensionless radius of the wavy portion of the channel system. Results have only been obtained for a Prandtl number 0.7 for Reynolds numbers of between 10 and 500 for various dimensionless wave radii. The solution shows that steady exists at low Reynolds numbers but that at a critical Reynolds unsteady flow exists, this critical Reynolds number depending on the dimensionless wave radius. The major concern in this work is with defining the conditions under which this unsteady flow develops and with the effect of the development of this unsteadiness on the pressure drop in the channel system.


Materials ◽  
2003 ◽  
Author(s):  
Sayavur I. Bakhtiyarov ◽  
Ruel A. Overfelt

Numerical simulation of decomposed gases through foam pattern was conducted using finite element analysis. A new kinetic model is proposed for gaseos phase flow between molten metal and foam material. The computations were performed for a wide range of Reynolds numbers. The results of the simulations are compared with the experiemental data obtained in this study.


1985 ◽  
Vol 160 ◽  
pp. 281-295 ◽  
Author(s):  
F. A. Milinazzo ◽  
P. G. Saffman

Computations of two-dimensional solutions of the Navier–Stokes equations are carried out for finite-amplitude waves on steady unidirectional flow. Several cases are considered. The numerical method employs pseudospectral techniques in the streamwise direction and finite differences on a stretched grid in the transverse direction, with matching to asymptotic solutions when unbounded. Earlier results for Poiseuille flow in a channel are re-obtained, except that attention is drawn to the dependence of the minimum Reynolds number on the physical constraint of constant flux or constant pressure gradient. Attempts to calculate waves in Couette flow by continuation in the velocity of a channel wall fail. The asymptotic suction boundary layer is shown to possess finite-amplitude waves at Reynolds numbers orders of magnitude less than the critical Reynolds number for linear instability. Waves in the Blasius boundary layer and unsteady Rayleigh profile are calculated by employing the artifice of adding a body force to cancel the spatial or temporal growth. The results are verified by comparison with perturbation analysis in the vicinity of the linear-instability critical Reynolds numbers.


Author(s):  
Carlos Marchi ◽  
Cosmo D. Santiago ◽  
Carlos Alberto Rezende de Carvalho Junior

Abstract The incompressible steady-state fluid flow inside a lid-driven square cavity was simulated using the mass conservation and Navier-Stokes equations. This system of equations is solved for Reynolds numbers of up to 10,000 to the accuracy of the computational machine round-off error. The computational model used was the second-order accurate finite volume method. A stable solution is obtained using the iterative multigrid methodology with 8192 × 8192 volumes, while degree-10 interpolation and Richardson extrapolation were used to reduce the discretization error. The solution vector comprised five entries of velocities, pressure, and location. For comparison purposes, 65 different variables of interest were chosen, such as velocity profile, its extremum values and location, extremum values and location of the stream function. The discretization error for each variable of interest was estimated using two types of estimators and their apparent order of accuracy. The variations of the 11 selected variables are shown across 38 Reynolds number values between 0.0001 and 10,000. In this study, we provide a more accurate determination of the Reynolds number value at which the upper secondary vortex appears. The results of this study were compared with those of several other studies in the literature. The current solution methodology was observed to produce the most accurate solution till date for a wide range of Reynolds numbers.


2021 ◽  
Author(s):  
Matthew Searle ◽  
Arnab Roy ◽  
James Black ◽  
Doug Straub ◽  
Sridharan Ramesh

Abstract In this paper, experimental and numerical investigations of three variants of internal cooling configurations — dimples only, ribs only and ribs with dimples have been explored at process conditions (96°C and 207bar) with sCO2 as the coolant. The designs were chosen based on a review of advanced internal cooling features typically used for air-breathing gas turbines. The experimental study described in this paper utilizes additively manufactured square channels with the cooling features over a range of Reynolds number from 80,000 to 250,000. Nusselt number is calculated in the experiments utilizing the Wilson Plot method and three heat transfer characteristics — augmentation in Nusselt number, friction factor and overall Thermal Performance Factor (TPF) are reported. To explore the effect of surface roughness introduced due to additive manufacturing, two baseline channel flow cases are considered — a conventional smooth tube and an additively manufactured square tube. A companion computational fluid dynamics (CFD) simulation is also performed for the corresponding cooling configurations reported in the experiments using the Reynolds Averaged Navier Stokes (RANS) based turbulence model. Both experimental and computational results show increasing Nusselt number augmentation as higher Reynolds numbers are approached, whereas prior work on internal cooling of air-breathing gas turbines predict a decay in the heat transfer enhancement as Reynolds number increases. Comparing cooling features, it is observed that the “ribs only” and “ribs with dimples” configurations exhibit higher Nusselt number augmentation at all Reynolds numbers compared to the “dimples only” and the “no features” configurations. However, the frictional losses are almost an order of magnitude higher in presence of ribs.


2021 ◽  
Vol 930 ◽  
Author(s):  
Kartik P. Iyer ◽  
Katepalli R. Sreenivasan ◽  
P.K. Yeung

Using direct numerical simulations performed on periodic cubes of various sizes, the largest being $8192^3$ , we examine the nonlinear advection term in the Navier–Stokes equations generating fully developed turbulence. We find significant dissipation even in flow regions where nonlinearity is locally absent. With increasing Reynolds number, the Navier–Stokes dynamics amplifies the nonlinearity in a global sense. This nonlinear amplification with increasing Reynolds number renders the vortex stretching mechanism more intermittent, with the global suppression of nonlinearity, reported previously, restricted to low Reynolds numbers. In regions where vortex stretching is absent, the angle and the ratio between the convective vorticity and solenoidal advection in three-dimensional isotropic turbulence are statistically similar to those in the two-dimensional case, despite the fundamental differences between them.


Author(s):  
Patricia Streufert ◽  
Terry X. Yan ◽  
Mahdi G. Baygloo

Local turbulent convective heat transfer from a flat plate to a circular impinging air jet is numerically investigated. The jet-to-plate distance (L/D) effect on local heat transfer is the main focus of this study. The eddy viscosity V2F turbulence model is used with a nonuniform structured mesh. Reynolds-Averaged Navier-Stokes equations (RANS) and the energy equation are solved for axisymmetric, three-dimensional flow. The numerical solutions obtained are compared with published experimental data. Four jet-to-plate distances, (L/D = 2, 4, 6 and 10) and seven Reynolds numbers (Re = 7,000, 15,000, 23,000, 50,000, 70,000, 100,000 and 120,000) were parametrically studied. Local and average heat transfer results are analyzed and correlated with Reynolds number and the jet-to-plate distance. Results show that the numerical solutions matched experimental data best at low jet-to-plate distances and lower Reynolds numbers, decreasing in ability to accurately predict the heat transfer as jet-to-plate distance and Reynolds number was increased.


Sign in / Sign up

Export Citation Format

Share Document