Study of Fully Developed Incompressible Flow in Curved Ducts, Using a Multi-Grid Technique

1987 ◽  
Vol 109 (3) ◽  
pp. 226-236 ◽  
Author(s):  
K. N. Ghia ◽  
U. Ghia ◽  
C. T. Shin
Keyword(s):  
Cross Sections ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Cross Flow ◽  
Internal Flow ◽  
Navier Stokes ◽  
Dean Number ◽  
Similarity Parameter ◽  
Fine Grid ◽  
Curved Ducts

Fully developed flows inside curved ducts of rectangular as well as polar cross sections have been analyzed using the Navier-Stokes equations in terms of the axial velocity and vorticity and the cross-flow stream function. Numerical solutions of the three second-order coupled elliptic partial differential equations governing this flow have been obtained using efficient numerical schemes. For curved-duct flows, the similarity parameter of significance is the Dean number K, rather than the Reynolds number Re. Results have been obtained for curved ducts with square cross sections for K up to 900 which, in the present study, corresponds to Re = 9,000 for this internal flow configuration. The fine-grid calculations show that, for square cross-section ducts, Dean’s instability occurs at K ≈ 125 and, further, that this phenomenon does not disappear even for K = 900. In ducts of polar cross sections, which are geometrically more representative of turbomachinery cascade passages, the phenomenon of Dean’s instability is not seen to occur for K up to 600.

Laminar Incompressible Viscous Flow in Curved Ducts of Regular Cross-Sections

1977 ◽  
Vol 99 (4) ◽  
pp. 640-648 ◽  
Author(s):  
K. N. Ghia ◽  
J. S. Sokhey
Keyword(s):  
Cross Sections ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Navier Stokes ◽  
Dean Number ◽  
Navier Stokes Equations ◽  
Flow Parameters ◽  
Incompressible Viscous Flow ◽  
Square Ducts ◽  
Curved Ducts

The laminar three-dimensional flow in curved ducts has been analyzed for an incompressible viscous fluid. The mathematical model is formulated using three-dimensional parabolized Navier-Stokes equations. The equations are generalized using two indices which permit the choice of Cartesian or cylindrical coordinate systems and straight or curved ducts. The solutions are obtained numerically using an ADI method for a number of duct geometries and flow parameters. The study presents detailed results for developing laminar flow in rectangular curved ducts; also, the effect of longitudinal curvature on secondary flow is fully analyzed. An investigation is made of the occurrence of Dean’s instability and, for curved square ducts, it is found to first appear at Dean number ≃ 143.

Some numerical experiments on developing laminar flow in circular-sectioned bends

1985 ◽  
Vol 154 ◽  
pp. 357-375 ◽  
Author(s):  
J. A. C. Humphrey ◽  
H. Iacovides ◽  
B. E. Launder
Keyword(s):  
Shear Stress ◽  
Laminar Flow ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Inviscid Flow ◽  
Streamwise Velocity ◽  
Navier Stokes ◽  
Dean Number ◽  
Navier Stokes Equations ◽  
Numerical Predictions

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.

Complex solutions of the Dean equations and non-uniqueness at all Reynolds numbers

2017 ◽  
Vol 818 ◽  
pp. 241-259 ◽  
Author(s):  
F. A. T. Boshier ◽  
A. J. Mestel
Keyword(s):  
Numerical Solutions ◽  
Stokes Equations ◽  
Reynolds Numbers ◽  
Navier Stokes ◽  
Dean Number ◽  
Navier Stokes Equations ◽  
Vortex Solution ◽  
Vortex Solutions ◽  
Complex Solutions ◽  
Unknown Solution

Steady incompressible flow down a slowly curving circular pipe is considered, analytically and numerically. Both real and complex solutions are investigated. Using high-order Hermite–Padé approximants, the Dean series solution is analytically continued outside its circle of convergence, where it predicts a complex solution branch for real positive Dean number, $K$. This is confirmed by numerical solution. It is shown that other previously unknown solution branches exist for all $K>0$, which are related to an unforced complex eigensolution. This non-uniqueness is believed to be generic to the Navier–Stokes equations in most geometries. By means of path continuation, numerical solutions are followed around the complex $K$-plane. The standard Dean two-vortex solution is shown to lie on the same hypersurface as the eigensolution and the four-vortex solutions found in the literature. Elliptic pipes are considered and shown to exhibit similar behaviour to the circular case. There is an imaginary singularity limiting convergence of the Dean series, an unforced solution at $K=0$ and non-uniqueness for $K>0$, culminating in a real bifurcation.

Bifurcation in steady laminar flow through curved tubes

1982 ◽  
Vol 119 ◽  
pp. 475-490 ◽  
Author(s):  
K. Nandakumar ◽  
Jacob H. Masliyah
Keyword(s):  
Stokes Equations ◽  
Outer Wall ◽  
Dual Solutions ◽  
Navier Stokes ◽  
Dean Number ◽  
Navier Stokes Equations ◽  
Vortex Solution ◽  
Duct Geometry ◽  
Flow Through ◽  
Curved Ducts

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.

scholarly journals A RAYLEIGH–RITZ METHOD FOR NAVIER–STOKES FLOW THROUGH CURVED DUCTS

2019 ◽  
Vol 61 (1) ◽  
pp. 1-22 ◽  
Author(s):  
BRENDAN HARDING
Keyword(s):  
Cross Sections ◽  
Cell Separation ◽  
Ritz Method ◽  
Cross Flow ◽  
Axial Flow ◽  
Navier Stokes ◽  
Cross Sectional ◽  
Straightforward Method ◽  
Flow Through ◽  
Curved Ducts

We present a Rayleigh–Ritz method for the approximation of fluid flow in a curved duct, including the secondary cross-flow, which is well known to develop for nonzero Dean numbers. Having a straightforward method to estimate the cross-flow for ducts with a variety of cross-sectional shapes is important for many applications. One particular example is in microfluidics where curved ducts with low aspect ratio are common, and there is an increasing interest in nonrectangular duct shapes for the purpose of size-based cell separation. We describe functionals which are minimized by the axial flow velocity and cross-flow stream function which solve an expansion of the Navier–Stokes model of the flow. A Rayleigh–Ritz method is then obtained by computing the coefficients of an appropriate polynomial basis, taking into account the duct shape, such that the corresponding functionals are stationary. Whilst the method itself is quite general, we describe an implementation for a particular family of duct shapes in which the top and bottom walls are described by a polynomial with respect to the lateral coordinate. Solutions for a rectangular duct and two nonstandard duct shapes are examined in detail. A comparison with solutions obtained using a finite-element method demonstrates the rate of convergence with respect to the size of the basis. An implementation for circular cross-sections is also described, and results are found to be consistent with previous studies.

Receptivity to surface roughness near a swept leading edge

1999 ◽  
Vol 380 ◽  
pp. 141-168 ◽  
Author(s):  
S. SCOTT COLLIS ◽  
SANJIVA K. LELE
Keyword(s):  
Surface Roughness ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Convex Surface ◽  
Cross Flow ◽  
Leading Edge ◽  
Surface Curvature ◽  
Parallel Flow ◽  
Navier Stokes ◽  
Swept Wing

The formation of stationary cross flow vortices in a three-dimensional boundary layer due to surface roughness located near the leading edge of a swept wing is investigated using numerical solutions of the compressible Navier–Stokes equations. The numerical solutions are used to evaluate the accuracy of theoretical receptivity predictions which are based on the parallel-flow approximation. By reformulating the receptivity theory to include the effect of surface curvature, it is shown that convex surface curvature enhances receptivity. Comparisons of the parallel-flow predictions with Navier–Stokes solutions demonstrate that non-parallel effects strongly reduce the initial amplitude of stationary cross flow vortices. The curvature and non-parallel effects tend to counteract one another; but, for the cases considered here, the non-parallel effect dominates leading to significant over-prediction of receptivity by parallel-flow receptivity theory. We conclude from these results that receptivity theories must account for non-parallel effects in order to accurately predict the amplitude of stationary crossflow instability waves near the leading edge of a swept wing.

A numerical study of the interaction between unsteay free-stream disturbances and localized variations in surface geometry

1989 ◽  
Vol 209 ◽  
pp. 285-308 ◽  
Author(s):  
R. J. Bodonyi ◽  
W. J. C. Welch ◽  
P. W. Duck ◽  
M. Tadjfar
Keyword(s):  
Numerical Solutions ◽  
Free Stream ◽  
Stokes Equations ◽  
Numerical Study ◽  
Navier Stokes ◽  
Surface Geometry ◽  
Navier Stokes Equations ◽  
S Waves ◽  
Tollmien Schlichting

A numerical study of the generation of Tollmien-Schlichting (T–S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite–difference and spectral methods. The nonlinear steady flow is of the viscous–inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier–Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T–S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T–S waves.

NUMERICAL SOLUTIONS OF 2D AND 3D SLAMMING PROBLEMS

2021 ◽  
Vol 153 (A2) ◽  
Author(s):  
Q Yang ◽  
W Qiu
Keyword(s):  
Free Surface ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Type Equation ◽  
The Body ◽  
Navier Stokes ◽  
Time Step ◽  
Navier Stokes Equations ◽  
Highly Nonlinear ◽  
2D And 3D

Slamming forces on 2D and 3D bodies have been computed based on a CIP method. The highly nonlinear water entry problem governed by the Navier-Stokes equations was solved by a CIP based finite difference method on a fixed Cartesian grid. In the computation, a compact upwind scheme was employed for the advection calculations and a pressure-based algorithm was applied to treat the multiple phases. The free surface and the body boundaries were captured using density functions. For the pressure calculation, a Poisson-type equation was solved at each time step by the conjugate gradient iterative method. Validation studies were carried out for 2D wedges with various deadrise angles ranging from 0 to 60 degrees at constant vertical velocity. In the cases of wedges with small deadrise angles, the compressibility of air between the bottom of the wedge and the free surface was modelled. Studies were also extended to 3D bodies, such as a sphere, a cylinder and a catamaran, entering calm water. Computed pressures, free surface elevations and hydrodynamic forces were compared with experimental data and the numerical solutions by other methods.

Iterative method for solving fully-coupled fluid solid systems

2010 ◽  
pp. 653-670
Author(s):  
Elisabeth Longatte
Keyword(s):  
Stokes Equations ◽  
Cross Flow ◽  
Elastic Cylinder ◽  
Navier Stokes ◽  
Arbitrary Lagrangian Eulerian ◽  
Navier Stokes Equations ◽  
Ale Formulation ◽  
Volume Method ◽  
Fully Coupled ◽  
Solid System

This work is concerned with the modelling of the interaction of a fluid with a rigid or a flexible elastic cylinder in the presence of axial or cross-flow. A partitioned procedure is involved to perform the computation of the fully-coupled fluid solid system. The fluid flow is governed by the incompressible Navier-Stokes equations and modeled by using a fractional step scheme combined with a co-located finite volume method for space discretisation. The motion of the fluid domain is accounted for by a moving mesh strategy through an Arbitrary Lagrangian-Eulerian (ALE) formulation. Solid dyncamics is modeled by a finite element method in the linear elasticity framework and a fixed point method is used for the fluid solid system computation. In the present work two examples are presented to show the method robustness and efficiency.

Sign in / Sign up

Export Citation Format

Share Document