scholarly journals Numerical Simulation of Turbulent Jets With Rectangular Cross-Section

1998 ◽  
Vol 120 (2) ◽  
pp. 285-290 ◽  
Author(s):  
R. V. Wilson ◽  
A. O. Demuren

Three-dimensional turbulent jets with rectangular cross-section are simulated with a finite-difference numerical method. The full Navier-Stokes equations are solved at a low Reynolds number, whereas at a higher Reynolds number filtered forms of the equations are solved along with a sub-grid scale model to approximate effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporal discretization and a fourth-order compact scheme is used for spatial discretization. Divergence-free velocity field is obtained by solving a Poisson equation for pressure with the same spatial discretization scheme for consistent accuracy. Computations are performed for different inlet conditions which represent different types of jet forcing within the shear layer. The phenomenon of axis-switching is observed in some cases. At low Reynolds numbers, it is based on self-induction of the vorticity field, whereas at higher Reynolds numbers, the turbulent structure becomes the dominant mechanism in natural jets. Budgets of the mean streamwise velocity show that convection is balanced by gradients of the Reynolds stresses and the pressure.

2017 ◽  
Vol 836 ◽  
pp. 1009-1049 ◽  
Author(s):  
B. Viggiano ◽  
T. Dib ◽  
N. Ali ◽  
L. G. Mastin ◽  
R. B. Cal ◽  
...  

Geophysical flows occur over a large range of scales, with Reynolds numbers and Richardson numbers varying over several orders of magnitude. For this study, jets of different densities were ejected vertically into a large ambient region, considering conditions relevant to some geophysical phenomena. Using particle image velocimetry, the velocity fields were measured for three different gases exhausting into air – specifically helium, air and argon. Measurements focused on both the jet core and the entrained ambient. Experiments considered relatively low Reynolds numbers from approximately 1500 to 10 000 with Richardson numbers near 0.001 in magnitude. These included a variety of flow responses, notably a nearly laminar jet, turbulent jets and a transitioning jet in between. Several features were studied, including the jet development, the local entrainment ratio, the turbulent Reynolds stresses and the eddy strength. Compared to a fully turbulent jet, the transitioning jet showed up to 50 % higher local entrainment and more significant turbulent fluctuations. For this condition, the eddies were non-axisymmetric and larger than the exit radius. For turbulent jets, the eddies were initially smaller and axisymmetric while growing with the shear layer. At lower turbulent Reynolds number, the turbulent stresses were more than 50 % higher than at higher turbulent Reynolds number. In either case, the low-density jet developed faster than a comparable non-buoyant jet. Quadrant analysis and proper orthogonal decomposition were also utilized for insight into the entrainment of the jet, as well as to assess the energy distribution with respect to the number of eigenmodes. Reynolds shear stresses were dominant in Q1 and Q3 and exhibited negligible contributions from the remaining two quadrants. Both analysis techniques showed that the development of stresses downstream was dependent on the Reynolds number while the spanwise location of the stresses depended on the Richardson number.


1966 ◽  
Vol 24 (1) ◽  
pp. 113-151 ◽  
Author(s):  
Odus R. Burggraf

The viscous structure of a separated eddy is investigated for two cases of simplified geometry. In § 1, an analytical solution, based on a linearized model, is obtained for an eddy bounded by a circular streamline. This solution reveals the flow development from a completely viscous eddy at low Reynolds number to an inviscid rotational core at high Reynolds number, in the manner envisaged by Batchelor. Quantitatively, the solution shows that a significant inviscid core exists for a Reynolds number greater than 100. At low Reynolds number the vortex centre shifts in the direction of the boundary velocity until the inviscid core develops; at large Reynolds number, the inviscid vortex core is symmetric about the centre of the circle, except for the effect of the boundary-layer displacement-thickness. Special results are obtained for velocity profiles, skin-friction distribution, and total power dissipation in the eddy. In addition, results of the method of inner and outer expansions are compared with the complete solution, indicating that expansions of this type give valid results for separated eddies at Reynolds numbers greater than about 25 to 50. The validity of the linear analysis as a description of separated eddies is confirmed to a surprising degree by numerical solutions of the full Navier–Stokes equations for an eddy in a square cavity driven by a moving boundary at the top. These solutions were carried out by a relaxation procedure on a high-speed digital computer, and are described in § 2. Results are presented for Reynolds numbers from 0 to 400 in the form of contour plots of stream function, vorticity, and total pressure. At the higher values of Reynolds number, an inviscid core develops, but secondary eddies are present in the bottom corners of the square at all Reynolds numbers. Solutions of the energy equation were obtained also, and isotherms and wall heat-flux distributions are presented graphically.


Author(s):  
S. P. Vanka ◽  
C. M. Winkler ◽  
J. Coffman ◽  
E. Linderman ◽  
S. Mahjub ◽  
...  

We present two new designs of compact mixers that can provide good mixing at low Reynolds numbers encountered in many microfluidic devices. The new designs benefit from curvature induced cross-stream vortices to enhance mixing of two co-flowing streams of fluids arranged side by side. One of the designs is a spiral of rectangular cross-section, while the other is a series of concentric circular channels arranged as a labyrinth. Both utilize the formation of sustained secondary flows to enhance mixing between two streams. Currently, the devices are fabricated in aluminum using standard machining techniques. However, they can be reduced further in size using standard microfabrication techniques. Mixing experiments were conducted in these channels at a Reynolds number of 6.8 using two sucrose solutions, one of which was laced with Rhodamine 6G dye. Compared to a experiment in an equivalent straight channel, a significant enhancement in the mixing of the two streams, as indicated by the intensity of the second fluid’s color, was observed. The present designs provide a compact and easy-to-fabricate alternative to various other concepts proposed in literature.


2019 ◽  
Vol 868 ◽  
pp. 141-175 ◽  
Author(s):  
Haithem Taha ◽  
Amir S. Rezaei

The application of the Kutta condition to unsteady flows has been controversial over the years, with increased research activities over the 1970s and 1980s. This dissatisfaction with the Kutta condition has been recently rejuvenated with the increased interest in low-Reynolds-number, high-frequency bio-inspired flight. However, there is no convincing alternative to the Kutta condition, even though it is not mathematically derived. Realizing that the lift generation and vorticity production are essentially viscous processes, we provide a viscous extension of the classical theory of unsteady aerodynamics by relaxing the Kutta condition. We introduce a trailing-edge singularity term in the pressure distribution and determine its strength by using the triple-deck viscous boundary layer theory. Based on the extended theory, we develop (for the first time) a theoretical viscous (Reynolds-number-dependent) extension of the Theodorsen lift frequency response function. It is found that viscosity induces more phase lag to the Theodorsen function particularly at high frequencies and low Reynolds numbers. The obtained theoretical results are validated against numerical laminar simulations of Navier–Stokes equations over a sinusoidally pitching NACA 0012 at low Reynolds numbers and using Reynolds-averaged Navier–Stokes equations at relatively high Reynolds numbers. The physics behind the observed viscosity-induced lag is discussed in relation to wake viscous damping, circulation development and the Kutta condition. Also, the viscous contribution to the lift is shown to significantly decrease the virtual mass, particularly at high frequencies and Reynolds numbers.


2010 ◽  
Vol 668 ◽  
pp. 33-57 ◽  
Author(s):  
LAURA GUGLIELMINI ◽  
R. RUSCONI ◽  
S. LECUYER ◽  
H. A. STONE

In recent microfluidic experiments with solutions of bacteria we observed the formation of biofilms in the form of thread-like structures, called ‘streamers’, which float in the middle plane of the channel and are connected to the side walls at the inner corners. Motivated by this observation, we discuss here the pressure-driven low-Reynolds-number flow around a corner bounded by the walls of a channel with rectangular cross-section. We numerically solve the flow field in a channel of constant cross-section, which exhibits 90° sharp corners, or turns with constant curvature, or portions with slowly changing curvature along the flow direction, for finite, but small, values of the Reynolds numbers and including the limit of vanishingly small Reynolds numbers. In addition, we develop a matched asymptotic expansion solution for the flow around two boundaries intersecting at an angle α and spanning the small gap h between two horizontal plates. We illustrate the basic features of the flow in these channel geometries by describing the three-dimensional velocity field and the distribution of streamwise vorticity and helicity, and comparing the numerical solutions with predictions based on the asymptotic approach. We demonstrate that near a corner or a change in the curvature of the side wall the flow is three-dimensional and pairs of counter-rotating vortical structures are present, as identified by Balsa (J. Fluid Mech., vol. 372, 1998, p. 25). Finally, we discuss how this secondary flow depends on the significant geometric parameters, the aspect ratio of the channel cross-section, the radius of curvature of the turn and, more generally, the variation of the curvature of the channel side boundary. We believe that these three-dimensional secondary flow structures are relevant to transport problems where accumulation of material at the boundary is possible.


1994 ◽  
Vol 47 (8) ◽  
pp. 307-365 ◽  
Author(s):  
Mohamed Gad-el-Hak ◽  
Promode R. Bandyopadhyay

This paper reviews the state of the art of Reynolds number effects in wall-bounded shear-flow turbulence, with particular emphasis on the canonical zero-pressure-gradient boundary layer and two-dimensional channel flow problems. The Reynolds numbers encountered in many practical situations are typically orders of magnitude higher than those studied computationally or even experimentally. High-Reynolds number research facilities are expensive to build and operate and the few existing are heavily scheduled with mostly developmental work. For wind tunnels, additional complications due to compressibility effects are introduced at high speeds. Full computational simulation of high-Reynolds number flows is beyond the reach of current capabilities. Understanding of turbulence and modeling will continue to play vital roles in the computation of high-Reynolds number practical flows using the Reynolds-averaged Navier-Stokes equations. Since the existing knowledge base, accumulated mostly through physical as well as numerical experiments, is skewed towards the low Reynolds numbers, the key question in such high-Reynolds number modeling as well as in devising novel flow control strategies is: what are the Reynolds number effects on the mean and statistical turbulence quantities and on the organized motions? Since the mean flow review of Coles (1962), the coherent structures, in low-Reynolds number wall-bounded flows, have been reviewed several times. However, the Reynolds number effects on the higher-order statistical turbulence quantities and on the coherent structures have not been reviewed thus far, and there are some unresolved aspects of the effects on even the mean flow at very high Reynolds numbers. Furthermore, a considerable volume of experimental and full-simulation data have been accumulated since 1962. The present article aims at further assimilation of those data, pointing to obvious gaps in the present state of knowledge and highlighting the misunderstood as well as the ill-understood aspects of Reynolds number effects.


Author(s):  
Biswash Shrestha ◽  
Nawraj Bhattarai

This study aims to achieve an improved airfoil performance at low Reynolds number, and to determine the optimum position and size of rectangular cross-section burst control plate (BCP) to suppress stall in airfoil. The type of airfoil used in the present study is NACA0015 (National Advisory Committee for Aeronautics) airfoil with 200 mm of chord (c) length. Here, rectangular cross-section burst control plates with different sizes and at different locations are investigated numerically at the low Reynolds number of 1.6×105. Total of three positions (0.05c, 0.1c and 0.2c from the leading edge of airfoil), and four sizes (with heights 0.3 mm, 0.7mm, 1mm and 1.5 mm, and constant width 4 mm) of rectangular BCPs are simulated in ANSYS Fluent software using Transition SST model. The results indicate that the rectangular cross-section burst control plate is an effective device in the suppression of airfoil stall. For 0.7 mm and 1 mm height BCPs, the stall angle is postponed by 2° for all positions, while for 0.3 mm and 1.5 mm height BCPs, the reduction in sudden fall of lift can be observed but at the cost of reduction in maximum lift coefficient. Among various configurations, the 1mm height BCP located at 0.2c position is found to be most effective in the suppression of stall.


2010 ◽  
Vol 1 (1-2) ◽  
pp. 15-20 ◽  
Author(s):  
B. Bolló

Abstract The two-dimensional flow around a stationary heated circular cylinder at low Reynolds numbers of 50 < Re < 210 is investigated numerically using the FLUENT commercial software package. The dimensionless vortex shedding frequency (St) reduces with increasing temperature at a given Reynolds number. The effective temperature concept was used and St-Re data were successfully transformed to the St-Reeff curve. Comparisons include root-mean-square values of the lift coefficient and Nusselt number. The results agree well with available data in the literature.


1991 ◽  
Vol 113 (4) ◽  
pp. 608-616 ◽  
Author(s):  
H. M. Jang ◽  
J. A. Ekaterinaris ◽  
M. F. Platzer ◽  
T. Cebeci

Two methods are described for calculating pressure distributions and boundary layers on blades subjected to low Reynolds numbers and ramp-type motion. The first is based on an interactive scheme in which the inviscid flow is computed by a panel method and the boundary layer flow by an inverse method that makes use of the Hilbert integral to couple the solutions of the inviscid and viscous flow equations. The second method is based on the solution of the compressible Navier–Stokes equations with an embedded grid technique that permits accurate calculation of boundary layer flows. Studies for the Eppler-387 and NACA-0012 airfoils indicate that both methods can be used to calculate the behavior of unsteady blade boundary layers at low Reynolds numbers provided that the location of transition is computed with the en method and the transitional region is modeled properly.


Sign in / Sign up

Export Citation Format

Share Document