On laminar flow in curved semicircular ducts

1980 ◽  
Vol 99 (3) ◽  
pp. 469-479 ◽  
Author(s):  
Jacob H. Masliyah

Calculations of the flow field under laminar conditions in a helical semicircular duct have been made by numerically solving the Navier–Stokes equations. With the flat wall of the duct being the outer wall, the solution of the momentum equations for Dean numbers below 105 gave, for the secondary flow, twin counter-rotating vortices of Taylor–Goertler type. However, above a Dean number of Dn = 105, two solutions were possible. One solution was similar to that obtained for Dn < 105. The other solution revealed four vortices for the secondary flow. For Dn > 105, convergence to either flow pattern depended on the initial guess used in the numerical solution. Flow visualization confirmed the possibility of the presence of both types of secondary flow patterns.

Proceedings ◽  
2018 ◽  
Vol 2 (11) ◽  
pp. 577 ◽  
Author(s):  
Evangelos Karvelas ◽  
Christos Liosis ◽  
Theodoros Karakasidis ◽  
Ioannis Sarris

A possible solution for water purification from heavy metals is to capture them by using nanoparticles in microfluidic ducts. In this technique, heavy metal capture is achieved by effectively mixing two streams, a nanoparticle solution and the contaminated water. In the present work, particles and water mixing is numerically studied for various inlet velocity ratios and inflow angles of the two streams. The Navier-Stokes equations are solved for the water flow while the discrete motion of particles is evaluated by a Lagrangian method. Results showed that as the velocity ratio between the inlet streams increases, by increasing the particles solution flow, the mixing of particles with the contaminated water is increased. Thus, nanoparticles are more uniformly distributed in the duct. On the other hand, angle increase between the inflow streams ducts is found to be less significant.


1996 ◽  
Vol 314 ◽  
pp. 227-246 ◽  
Author(s):  
Philip A. J. Mees ◽  
K. Nandakumar ◽  
J. H. Masliyah

Steady developing flow of an incompressible Newtonian fluid in a curved duct of square cross-section (the Dean problem) is investigated both experimentally and numerically. This study is a continuation of the work by Bara, Nandakumar & Masliyah (1992) and is focused on flow rates between Dn = 200 and Dn = 600 (Dn = Re/(R/a)1/2, where Re is the Reynolds number, R is the radius of curvature of the duct and a is the duct dimension; the curvature ratio, R/a, is 15.1).Numerical simulations based on the steady three-dimensional Navier – Stokes equations predict the development of a 6-cell secondary flow pattern above a Dean number of 350. The 6-cell state consists of two large Ekman vortices and two pairs of small Dean vortices near the outer wall that result from the primary instability that is of centrifugal nature. The 6-cell flow state develops near θ = 80° and breaks down symmetrically into a 2-cell flow pattern.The apparatus used to verify the simulations had a duct dimension of 1.27 cm and a streamwise length of 270°. At a Dean number of 453, different velocity profiles of the 6-cell flow state at θ = 90° and spanwise profiles of the streamwise velocity at every 20° were measured using a laser-Doppler anemometer. All measured velocity profiles, as well as flow visualization of secondary flow patterns, are in very good agreement with the simulations, indicating that the parabolized Navier – Stokes equations give an accurate description of the flow.Based on the similarity with boundary layer flow over a concave wall (the Görtler problem), it is suggested that the transition to the 6-cell flow state is the result of a decreasing spanwise wavelength of the Dean vortices with increasing flow rate. A numerical stability analysis shows that the 6-cell flow state is unconditionally unstable. This is the first time that detailed experiments and simulations of the development of a 6-cell flow state are reported.


1982 ◽  
Vol 119 ◽  
pp. 475-490 ◽  
Author(s):  
K. Nandakumar ◽  
Jacob H. Masliyah

The occurrence of dual solutions in curved ducts is investigated through a numerical solution of the Navier-Stokes equations in a bipolar-toroidal co-ordinate system. With the shape of duct being the region formed by the natural co-ordinate surfaces, it was possible to alter the duct geometry gradually and preserve the prevailing form of the velocity field, in a manner suggested by Benjamin (1978).In addition to the Dean number Dn = Re/Rc½, a geometrical parameter that defines the shape of the duct was also varied systematically to study the bifurcation of a two-vortex solution into a two- and four-vortex solution. Dual solutions have been found for all geometrical shapes investigated here. Of particular interest are the shapes of a full circle and a semicircle with a curved outer wall.


1977 ◽  
Vol 83 (3) ◽  
pp. 509-527 ◽  
Author(s):  
J. A. C. Humphrey ◽  
A. M. K. Taylor ◽  
J. H. Whitelaw

Calculated values of the three velocity components and measured values of the longitudinal component are reported for the flow of water in a 90° bend of 40 x 40mm cross-section; the bend had a mean radius of 92mm and was located downstream of a 1[sdot ]8m and upstream of a 1[sdot ]2m straight section. The experiments were carried out at a Reynolds number, based on the hydraulic diameter and bulk velocity, of 790 (corresponding to a Dean number of 368). Flow visualization was used to identify qualitatively the characteristics of the flow and laser-Doppler anemometry to quantify the velocity field. The results confirm and quantify that the location of maximum velocity moves from the centre of the duct towards the outer wall and, in the 90° plane, is located around 85% of the duct width from the inner wall. Secondary velocities up to 65% of the bulk longitudinal velocity were calculated and small regions of recirculation, close to the outer corners of the duct and in the upstream region, were also observed.The calculated results were obtained by solving the Navier–Stokes equations in cylindrical co-ordinates. They are shown to exhibit the same trends as the experiments and to be in reasonable quantitative agreement even though the number of node points used to discretize the flow for the finite-difference solution of the differential equations was limited by available computer time and storage. The region of recirculation observed experimentally is confirmed by the calculations. The magnitude of the various terms in the equations is examined to determine the extent to which the details of the flow can be represented by reduced forms of the Navier–Stokes equations. The implications of the use of so-called ‘partially parabolic’ equations and of potential- and rotational-flow analysis of an ideal fluid are quantified.


2019 ◽  
Vol 16 (32) ◽  
pp. 345-361
Author(s):  
B. A. UNASPEKOV ◽  
Z. O. ZHUMADILOVA ◽  
S. S. AUELBEKOV ◽  
A. S. TAUBALDIEVA ◽  
G. B. ALDABERGENOVA

A study was made of the nature of the movement of air mass in the space between the device that warms the room and the confining outer wall of the room. The hydrodynamic motion of the air mass was simulated on the basis of the two-dimensional in space Navier-Stokes equations and the convective heat conduction equation to find out a detailed picture of the physical processes that occur. Also for some particular cases, analytical solutions were obtained for the conjugate problem, where the movement of air is caused by its thermal expansion and the action of gravity in areas with different densities (Archimedes' forces). The simulation of more complex types of air movement with the formation of vortex flows using numerical methods and the corresponding program codes written in C++. It was found that the movement of air mass in the space between the device and the fence strongly depends not only on the temperature of the device, the wall, and the outside air. It was revealed that hydrodynamics and heat transfer are significantly influenced by two geometrical parameters: the distance between the device and the wall, and the distance from the flooring to the bottom of the device. In particular, a thin layer of the downward flow of cold air along the fence and above the floor is possible. The condition for the occurrence of such a downward flow has been found, it is determined by the ratio of the temperature of the wall, the device, and the average temperature in the room.


1995 ◽  
Vol 1 (3-4) ◽  
pp. 225-235 ◽  
Author(s):  
M. J. Braun ◽  
M. Dzodzo

The flow in a hydrostatic pocket is numerically simulated using a dimensionless formulation of the 2-D Navier-Stokes equations written in primitive variables, for a body fitted coordinates system, and applied through a collocated grid. In essence, we continue the work of Braun et al. 1993a, 1993b] and extend it to the study of the effects of the pocket geometric format on the flow pattern and pressure distribution. The model includes the coupling between the pocket flow and a finite length feedline flow, on one hand, and the pocket and its adjacent lands on the other hand. In this context we shall present, on a comparative basis, the flow and the pressure patterns at the runner surface for square, ramped-Rayleigh step, and arc of circle pockets. Geometrically all pockets have the same footprint, same lands length, and same capillary feedline. The numerical simulation uses the Reynolds number based on the lid(runner) velocity and the inlet jet strengthFas the dynamic similarity parameters. The study aims at establishing criteria for the optimization of the pocket geometry in the larger context of the performance of a hydrostatic bearing.


2019 ◽  
Vol 141 (12) ◽  
Author(s):  
Filipe S. Pereira ◽  
Luís Eça ◽  
Guilherme Vaz

The importance of the turbulence closure to the modeling accuracy of the partially-averaged Navier–Stokes equations (PANS) is investigated in prediction of the flow around a circular cylinder at Reynolds number of 3900. A series of PANS calculations at various degrees of physical resolution is conducted using three Reynolds-averaged Navier–Stokes equations (RANS)-based closures: the standard, shear-stress transport (SST), and turbulent/nonturbulent (TNT) k–ω models. The latter is proposed in this work. The results illustrate the dependence of PANS on the closure. At coarse physical resolutions, a narrower range of scales is resolved so that the influence of the closure on the simulations accuracy increases significantly. Among all closures, PANS–TNT achieves the lowest comparison errors. The reduced sensitivity of this closure to freestream turbulence quantities and the absence of auxiliary functions from its governing equations are certainly contributing to this result. It is demonstrated that the use of partial turbulence quantities in such auxiliary functions calibrated for total turbulent (RANS) quantities affects their behavior. On the other hand, the successive increase of physical resolution reduces the relevance of the closure, causing the convergence of the three models toward the same solution. This outcome is achieved once the physical resolution and closure guarantee the precise replication of the spatial development of the key coherent structures of the flow.


1985 ◽  
Vol 154 ◽  
pp. 357-375 ◽  
Author(s):  
J. A. C. Humphrey ◽  
H. Iacovides ◽  
B. E. Launder

The paper reports numerical solutions to a semi-elliptic truncation of the Navier–Stokes equations for the case of developing laminar flow in circular-sectioned bends over a range of Dean numbers. The ratios of bend radius to pipe radius are 7:1 and 20:1, corresponding with the configurations examined experimentally by Talbot and his co-workers in recent years. The semi-elliptic treatment facilitates a much finer grid than has been possible in earlier studies. Numerical accuracy has been further improved by assuming radial equilibrium over a thin sublayer immediately adjacent to the wall and by re-formulating the boundary conditions at the pipe centre.Streamwise velocity profiles at Dean numbers of 183 and 565 are in excellent agreement with laser-Doppler measurements by Agrawal, Talbot & Gong (1978). Good, albeit less complete, accord is found with the secondary velocities, though the differences that exist may be mainly due to the difficulty of making these measurements. The paper provides new information on the behaviour of the streamwise shear stress around the inner line of symmetry. Upstream of the point of minimum shear stress, our numerical predictions display a progressive shift towards the result of Stewartson, Cebici & Chang (1980) as the Dean number is successively raised. Downstream of the minimum, however, in contrast with the monotonic approach to an asymptotic level reported by Stewartson, the numerical solutions display a damped oscillatory behaviour reminiscent of those from Hawthorne's (1951) inviscid-flow calculations. The amplitude of the oscillation grows as the Dean number is raised.


2015 ◽  
Vol 772 ◽  
Author(s):  
L. Kahouadji ◽  
N. Périnet ◽  
L. S. Tuckerman ◽  
S. Shin ◽  
J. Chergui ◽  
...  

We report the first simulations of the Faraday instability using the full three-dimensional Navier–Stokes equations in domains much larger than the characteristic wavelength of the pattern. We use a massively parallel code based on a hybrid front-tracking/level-set algorithm for Lagrangian tracking of arbitrarily deformable phase interfaces. Simulations performed in square and cylindrical domains yield complex patterns. In particular, a superlattice-like pattern similar to those of Douady & Fauve (Europhys. Lett., vol. 6, 1988, pp. 221–226) and Douady (J. Fluid Mech., vol. 221, 1990, pp. 383–409) is observed. The pattern consists of the superposition of two square superlattices. We conjecture that such patterns are widespread if the square container is large compared with the critical wavelength. In the cylinder, pentagonal cells near the outer wall allow a square-wave pattern to be accommodated in the centre.


Author(s):  
Joon Sang Lee ◽  
Xiaofeng Xu ◽  
R. H. Pletcher

Flow in an annular pipe with and without a wall rotating about its axis was investigated at moderate Reynolds numbers. The compressible filtered Navier-Stokes equations were solved using a second order accurate finite volume method. Low Mach number preconditioning was used to enable the compressible code to work efficiently at low Mach numbers. A dynamic subgrid-scale stress model accounted for the subgrid-scale turbulence. When the outer wall rotated, a significant reduction of turbulent kinetic energy was realized near the rotating wall and the intensity of bursting effects appeared to decrease. This modification of the turbulent structures was related to the vortical structure changes near the rotating wall. It has been observed that the wall vortices were pushed in the direction of rotation and their intensity increased near the non-rotating wall. The consequent effect was to enhance the turbulent kinetic energy and increased the intensity of the heat transfer rate there.


Sign in / Sign up

Export Citation Format

Share Document