The effect of Reynolds number on boundary layer turbulence

A new facility for studying high Reynolds number incompressible turbulent boundary layer flows has been constructed. It consists of a moderately sized wind tunnel, completely enclosed by a pressure vessel, which can raise the ambient air pressure in and around the wind tunnel to 8 atmospheres. This results in a Reynolds number range of about 20:1, while maintaining incompressible flow. Results are presented for the zero pressure gradient flat plate boundary layer over a momentum thickness Reynolds number range 1500–15?000. Scaling issues for high Reynolds number non-equilibrium boundary layers are discussed, with data comparing the three-dimensional turbulent boundary layer flow over a swept bump at Reynolds numbers of 3800 and 8600. It is found that successful prediction of these types of flows must include length scales which do not scale on Reynolds number, but are inherent to the geometry of the flow.     

后台阶流动的数值模拟

访述了大涡模拟的基本思想，指出大涡模拟的效率主要取决于四个因素，即流动中须有大尺度涡存在、合理的计算格式、合适的滤波器和亚格子应力模型。在深入考虑粘性不可压缩流Navier—Stokes方程各个子项作用的基础上，提出二阶全展开Euler—Taylor—Galerkin有限元方法作为大涡模拟的离散格式，并采用Gauss滤波器，对典型算例——后台阶处的流动进行大涡模拟，计算结果与相关文献符合的很好。从计算结果还可以看出大涡模拟与二阶全展开ETG有限元方法的结合在捕捉涡系及反映涡动时变过程方面具有明显的优势，说明大涡模拟适合于边界几何形状复杂区域流动的模拟。同时应用二阶全展开ETG有限元方法对低雷诺数粘性不可压缩后台阶流动进行了计算，得到与相关文献符合良好的计算结果，即该方法也可独立用于对低雷诺数粘性不可压缩流动的计算。     

A mixture differential quadrature method for solving two-dimensional incompressible Navier-Stokes equations

Ｔｈｅｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄ（ＤＱＭ）ｐｒｏｐｏｓｅｄｂｙＲ．Ｂｅｌｌｍａｎ［１,２］ｈａｓｂｅｅｎｓｕｃｃｅｓｓｆｕｌｌｙｅｍｐｌｏｙｅｄｉｎｎｕｍｅｒｉｃａｌｃｏｍｐｕｔａｔｉｏｎｓｏｆｐｒｏｂｌｅｍｓｉｎｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｈｙｓｉｃａｌｓｃｉｅｎｃｅ．ＢｅｃａｕｓｅｔｈｅｉｎｆｏｒｍａｔｉｏｎｏｎａｌｌｇｒｉｄｐｏｉｎｔｓｉｓｕｓｅｄｔｏｆｉｔｔｈｅｄｅｒｉｖａｔｉｖｅｓａｔｇｒｉｄｐｏｉｎｔｓｉｎｔｈｅＤＱＭ,ｉｔｉｓｅｎｏｕｇｈｔｏｏｂｔａ…     

Fractional step method for solution of incompressible Navier-Stokes equations on unstructured triangular meshes

In this paper an implicit fractional step method for the solution of the two-dimensional, time-dependent, incompressible Navier-Stokes equations is presented. The current method was developed for use on an unstructured grid made up of triangles. The basic principles of this method are that the evaluation of the time evolution is split into intermediate steps and that for the spatial discretization of the flow equations a finite volume discretization on an unstructured triangular mesh is used. The present approach has been used to simulate viscous, laminar flows for various Reynolds numbers in test cases such as a backward-facing step, a square cavity and a channel with wavy boundaries.     

Numerical simulation of high‐Reynolds number flow around circular cylinders by a three‐step FEM–BEM model

An innovative computational model, developed to simulate high‐Reynolds number flow past circular cylinders in two‐dimensional incompressible viscous flows in external flow fields is described in this paper. The model, based on transient Navier–Stokes equations, can solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the projection method. The pressure is assumed to be zero at infinite boundary and the external flow field is simulated using a direct boundary element method (BEM) by solving a pressure Poisson equation. A three‐step finite element method (FEM) is used to solve the momentum equations of the flow. The present model is applied to simulate high‐Reynolds number flow past a single circular cylinder and flow past two cylinders in which one acts as a control cylinder. The simulation results are compared with experimental data and other numerical models and are found to be feasible and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical solution of viscous flows using integral equation methods

A formulation of the boundary element method for the solution of non-zero Reynolds number incompressible flows in which the non-linear terms are lumped together to form a forcing function is presented. Solutions can be obtained at low to moderate Reynolds numbers. The method was tested using the flow of a fluid in a two-dimensional converging channel (Hamel flow) for which an exact solution is available. An axisymmetric formulation is demonstrated by examining the drag experienced by a sphere held stationary in uniform flow. Performance of the method was satisfactory. New results for an axisymmetric free jet at zero Reynolds number obtained using the boundary element method are also included. The method is ideal for this type of free-surface problem.     

On the use of adaptive hierarchical meshes for numerical simulation of separated flows

This paper describes the use of adaptive hierarchical grids to predict incompressible separated flow at low Reynolds number. The grids consist of a quadtree system of hierarchical Cartesian meshes which are generated by recursive subdivision about seeding points. The governing equations are discretized in collocated primitive variable form using finite volumes and solved using a pressure correction scheme. The mesh is locally adapted at each time step, with panel division or removal dependent on the vorticity magnitude. The resulting grids have fine local resolution and are economical in array size. Results are presented for unidirectional, impulsively started flow past a circular and a square cylinder at various Reynolds numbers up to 5000 and 250 respectively. It is clear that hierarchical meshes may offer gains in efficiency when applied to complex flow domains or strongly sheared flows. However, as expected, the stepped approximation to curved boundaries resulting from the Cartesian quadtree representation adversely affects the accuracy of the results for flow past a circular cylinder. © 1998 John Wiley & Sons, Ltd.     

不可压缩湍流非对称脱体涡流场数值分析

采用拟压缩性方法，Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因式分解格式数值求解三维定常不可压Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程。对Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，采用Ｄｅｇａｎｉ－ｓｃｈｉｆｆ修正。计算绕尖头正切拱型旋成柱体的大迎角大雷诺数脱体涡流场，计算结果中的非对称脱体涡与实验相符．     

Implicit preconditioned high‐order compact scheme for the simulation of the three‐dimensional incompressible Navier–Stokes equations with pseudo‐compressibility method

This paper combines the pseudo‐compressibility procedure, the preconditioning technique for accelerating the time marching for stiff hyperbolic equations, and high‐order accurate central compact scheme to establish the code for efficiently and accurately solving incompressible flows numerically based on the finite difference discretization. The spatial scheme consists of the sixth‐order compact scheme and 10th‐order numerical filter operator for guaranteeing computational stability. The preconditioned pseudo‐compressible Navier–Stokes equations are marched temporally using the implicit lower–upper symmetric Gauss–Seidel time integration method, and the time accuracy is improved by the dual‐time step method for the unsteady problems. The efficiency and reliability of the present procedure are demonstrated by applications to Taylor decaying vortices phenomena, double periodic shear layer rolling‐up problem, laminar flow over a flat plate, low Reynolds number unsteady flow around a circular cylinder at Re = 200, high Reynolds number turbulence flow past the S809 airfoil, and the three‐dimensional flows through two 90°curved ducts of square and circular cross sections, respectively. It is found that the numerical results of the present algorithm are in good agreement with theoretical solutions or experimental data. Copyright © 2011 John Wiley & Sons, Ltd.     

Electrokinetic flow and energy conversion in a curved microtube

Curved channels are ubiquitous in microfluidic systems. The pressuredriven electrokinetic flow and energy conversion in a curved microtube are investigated analytically by using a perturbation analysis method under the assumptions of the small curvature ratio and the Reynolds number. The results indicate that the curvature of the microtube leads to a skewed pattern in the distribution of the electrical double layer (EDL) potential. The EDL potential at the outer side of the bend is larger than that at the inner side of the bend. The curvature shows an inhibitory effect on the magnitude of the streaming potential field induced by the pressure-driven flow. Since the spanwise pressure gradient is dominant over the inertial force, the resulting axial velocity profile is skewed into the inner region of the curved channel. Furthermore, the flow rate in a curved microtube could be larger than that in a straight one with the same pressure gradient and shape of cross section. The asymptotic solutions of the axial velocity and flow rate in the absence of the electrokinetic effect are in agreement with the classical results for low Reynolds number flows. Remarkably, the curved geometry could be beneficial to improving the electrokinetic energy conversion (EKEC) efficiency.     

Numerical simulation of turbulent free surface flow with two‐equation k–ε eddy‐viscosity models

This paper presents a finite difference technique for solving incompressible turbulent free surface fluid flow problems. The closure of the time‐averaged Navier–Stokes equations is achieved by using the two‐equation eddy‐viscosity model: the high‐Reynolds k–ε (standard) model, with a time scale proposed by Durbin; and a low‐Reynolds number form of the standard k–ε model, similar to that proposed by Yang and Shih. In order to achieve an accurate discretization of the non‐linear terms, a second/third‐order upwinding technique is adopted. The computational method is validated by applying it to the flat plate boundary layer problem and to impinging jet flows. The method is then applied to a turbulent planar jet flow beneath and parallel to a free surface. Computations show that the high‐Reynolds k–ε model yields favourable predictions both of the zero‐pressure‐gradient turbulent boundary layer on a flat plate and jet impingement flows. However, the results using the low‐Reynolds number form of the k–ε model are somewhat unsatisfactory. Copyright © 2004 John Wiley & Sons, Ltd.     

Stochastic equations for continuum and determination of hydraulic drag coefficients for smooth flat plate and smooth round tube with taking into account intensity and scale of turbulent flow

The stochastic equations of continuum are used for determining the hydraulic drag coefficients. As a result, the formulas for the hydraulic drag coefficients dependent on the turbulence intensity and scale instead of only on the Reynolds number are proposed for the classic flows of an incompressible fluid along a smooth flat plate and a round smooth tube. It is shown that the new expressions for the classical drag coefficients, which depend only on the Reynolds number, should be obtained from these new general formulas if to use the well-known experimental data for the initial turbulence. It is found that the limitations of classical empirical and semiempirical formulas for the hydraulic drag coefficients and their deviation from the experimental data depend on different parameters of initial fluctuations in the flow for different experiments in a wide range of Reynolds numbers. On the basis of these new dependencies, it is possible to explain that the differences between the experimental results for the fixed Reynolds number are caused by the difference in the values of flow fluctuations for each experiment instead of only due to the systematic error in the processing of experiments. Accordingly, the obtained general dependencies for the smooth flat plate and the smooth round tube can serve as the basis for clarifying the results of experiments and the experimental formulas, which used for continuum flows in different devices.     

Flow field calculation and its experimental verification for inertia stage of marine gas turbine air intake filter

A numerical method is presented in this paper for the solution of flow field in marine gas turbine air intake filtration channel. The flow field was successfully calculated by this method, and aerodynamic characteristics were obtained for various types of filtration channels. This work is expected to be of practical importance for the design of such filters. Upstream difference was adopted to discretize the non-conservative type N-S equation for two-dimensional, time-dependnet, viscous and incompressible flow, and the stability, convergence, accuracy and artificial viscosity of the resulting difference equation were examined. This equation can be used to calculate viscous flows with Reynolds number up to tens of thousands. Also presented in this paper is a calculation method for treating wall vortex at boundary inflection points. Careful studies show that calculation based on the difference equational and wall vortex treatment proposed here are in good agreement with experimental results.     

Development of highly accurate interpolation method for mesh‐free flow simulations II. Application of CIVA method to incompressible flow simulations

This is the second report on the development of a highly accurate interpolation method, which is called cubic interpolation with volume/area (CIVA) co‐ordinates, for mesh‐free flow simulations. In this paper, the method of determining the c‐parameter of CIVA using a constant curvature condition is first considered for the two‐ and three‐dimensional cases. A computation of a three‐dimensional passive scalar advection problem is performed for accuracy verification and for comparison with widely used methods. Then, an application algorithm of the CIVA method respecting incompressible fluid simulation is presented. As the incompressible condition based on Lagrangian approaches causes problems, in this paper we consider the condition based on the conventional Eulerian approach. The CIVA‐based incompressible flow simulation algorithm enables a highly accurate simulation of many kinds of problems that have complicated geometries and involve complicated phenomena. To confirm the facts, numerical analyzes are executed for some benchmark problems, namely flow in a square cavity, free surface sloshing and moving boundary problems in complex geometries. The results show that the method achieves high accuracy and has high flexibility, even for the flows involving high Reynolds number, complicated geometries, moving boundaries and free surfaces. Copyright © 2000 John Wiley & Sons, Ltd.     

A HYBRID VORTEX METHOD FOR FLOWS OVER A BLUFF BODY

A hybrid vortex method was developed to simulate the two-dimensional viscous incompressible flows over a bluff body numerically. It is based on a combination of the diffusion–vortex method and the vortex-in-cell method by dividing the flow field into two regions. In the region near the body surface the diffusion–vortex method is used to solve the Navier–Stokes equations, while the vortex-in-cell method is used in the exterior domain. Comparison with results obtained by the finite difference method, other vortex methods and experiments shows that the present method is well adapted to calculate two-dimensional external flows at high Reynolds number. It is capable of calculating not only the global characteristics of the separated flow but also the evolution of the fine structure of the flow field with time precisely. The influence of the grid system and region decomposition on the results will also be discussed. © by 1997 John Wiley & Sons, Ltd.     

Modelling of Power-law Fluid Flow Through Porous Media Using Smoothed Particle Hydrodynamics

The flow of non-Newtonian fluids through two-dimensional porous media is analyzed at the pore scale using the smoothed particle hydrodynamics (SPH) method. A fully explicit projection method is used to simulate incompressible flow. This study focuses on a shear-thinning power-law model (n < 1), though the method is sufficiently general to include other stress-shear rate relationships. The capabilities of the proposed method are demonstrated by analyzing a Poiseuille problem at low Reynolds numbers. Two test cases are also solved to evaluate validity of Darcy’s law for power-law fluids and to investigate the effect of anisotropy at the pore scale. Results show that the proposed algorithm can accurately simulate non-Newtonian fluid flows in porous media.     

Boundary element analysis of driven cavity flow for low and moderate Reynolds numbers

A boundary element method for steady two‐dimensional low‐to‐moderate‐Reynolds number flows of incompressible fluids, using primitive variables, is presented. The velocity gradients in the Navier–Stokes equations are evaluated using the alternatives of upwind and central finite difference approximations, and derivatives of finite element shape functions. A direct iterative scheme is used to cope with the non‐linear character of the integral equations. In order to achieve convergence, an underrelaxation technique is employed at relatively high Reynolds numbers. Driven cavity flow in a square domain is considered to validate the proposed method by comparison with other published data. Copyright © 2001 John Wiley & Sons, Ltd.     

Quasihomogeneous model of cavitation flows in diffuser channels

Starting with the experiments carried out by Reynolds in 1894, the flow in Venturi tubes has traditionally been used to study and demonstrate various forms of cavitation. Numerous authors have carried out experimental research on the various flow regimes in diffuser channels [1–7] or have investigated theoretical models of such flows [6, 8]. The occurrence and development of cavitation is closely associated with the phenomenon of turbulent separation complicated by the presence of two-phase flow in the dissipation zone. For a long time these effects were considered separately, until Gogish and Stepanov [9] proposed a single model of cavitation and separation based on the theory of intense interaction of an incompressible potential flow and a turbulent cavitation layer of variable density and embracing the various stages of cavitation. The object of this study is to demonstrate the possibilities of this model with reference to the simple example of flows accompanied by cavitation and separation in plane and axisymmetric diffuser channels of the Venturi tube type with straight and curved walls. The dissipative flow near the walls is described by a quasihomogeneous model of turbulent two-phase flow, in which the presence of two phases is taken into account only by varying the mean density. The potential core of the flow is considered in the one-dimensional formulation. The displacement thickness serves as the flow interaction parameter. The conditions of ocurrence and development of circulatory flows are determined. Examples of symmetrical and nonsymmetrical flows are presented.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 47–54, September–October, 1986.     

Simulation of three-dimensional incompressible flows in generalized curvilinear coordinates using a high-order compact finite-difference lattice Boltzmann method

In the present study, a high-order compact finite-difference lattice Boltzmann method is applied for accurately computing 3-D incompressible flows in the generalized curvilinear coordinates to handle practical and realistic geometries with curved boundaries and nonuniform grids. The incompressible form of the 3-D nineteen discrete velocity lattice Boltzmann method is transformed into the generalized curvilinear coordinates. Herein, a fourth-order compact finite-difference scheme and a fourth-order Runge-Kutta scheme are used for the discretization of the spatial derivatives and the temporal term, respectively, in the resulting 3-D nineteen discrete velocity lattice Boltzmann equation to provide an accurate 3-D incompressible flow solver. A high-order spectral-type low-pass compact filtering technique is applied to have a stable solution. All boundary conditions are implemented based on the solution of the governing equations in the 3-D generalized curvilinear coordinates. Numerical solutions of different 3-D benchmark and practical incompressible flow problems are performed to demonstrate the accuracy and performance of the solution methodology presented. Herein, the 2-D cylindrical Couette flow, the decay of a 3-D double shear wave, the cubic lid-driven cavity flow with nonuniform grids, the flow through a square duct with 90° bend and the flow past a sphere at different flow conditions are considered for validating the present computations. Numerical results obtained show the accuracy and robustness of the present solution methodology based on the implementation of the high-order compact finite-difference lattice Boltzman method in the generalized curvilinear coordinates for solving 3-D incompressible flows over practical and realistic geometries.     

Electromagnetohydrodynamic flows and mass transport in curved rectangular microchannels

Curved microchannels are often encountered in lab-on-chip systems because the effective axial channel lengths of such channels are often larger than those of straight microchannels for a given per unit chip length. In this paper, the effective diffusivity of a neutral solute in an oscillating electromagnetohydrodynamic(EMHD)flow through a curved rectangular microchannel is investigated theoretically. The flow is assumed as a creeping flow due to the extremely low Reynolds number in such microflo...