The simulation of secondary flow effects in turbulent non-circular passage flows

A finite difference method has been developed to predict the overall features of the local mean flow in fully developed turbulent non-circular passage flows. The main transport effects of secondary flow have been identified and simulated with diffusion transport in a simple way which eliminates solution of the cross-plane momentum and continuity equations and produces a compact calculation method. Predictions are presented for four different passage shapes and are discussed in relation to experimental measurements and predictions from other more complex methods. Although some minor details were not predicted, the main effects of secondary flow on the mean flow were found to have been quite well simulated, yielding predictions that are in reasonable overall agreement with experiment.     

A finite element solution of the time-dependent incompressible Navier–Stokes equations using a modified velocity correction method

A finite element solution of the two-dimensional incompressible Navier–Stokes equations has been developed. The present method is a modified velocity correction approach. First an intermediate velocity is calculated, and then this is corrected by the pressure gradient which is the solution of a Poisson equation derived from the continuity equation. The novelty, in this paper, is that a second-order Runge–Kutta method for time integration has been used. Discretization in space is carried out by the Galerkin weighted residual method. The solution is in terms of primitive variables, which are approximated by polynomial basis functions defined on three-noded, isoparametric triangular elements. To demonstrate the present method, two examples are provided. Results from the first example, the driven cavity flow problem, are compared with previous works. Results from the second example, uniform flow past a cylinder, are compared with experimental data.     

A spectral collocation method for compressible,non-similar boundary layers

An efficient and highly accurate algorithm based on a spectral collocation method is developed for numerical solution of the compressible, two-dimensional and axisymmetric boundary layer equations. The numerical method incorporates a fifth-order, fully implicit marching scheme in the streamwise (timelike) dimension and a spectral collocation method based on Chebyshev polynomial expansions in the wall-normal (spacelike) dimension. The discrete governing equations are cast in residual form and the residuals are minimized at each marching step by a preconditioned Richardson iteration scheme which fully couples energy, momentum and continuity equations. Preconditioning on the basis of the finite difference analogues of the governing equations results in a computationally efficient iteration with acceptable convergence properties. A practical application of the algorithm arises in the area of compressible linear stability theory, in the investigation of the effects of transverse curvature on the stability of flows over axisymmetric bodies. The spectral collocation algorithm is used to derive the non-similar mean velocity and temperature profiles in the boundary layer of a ‘fuselage’ (cylinder) in a high-speed (Mach 5) flow parallel to its axis. The stability of the flow is shown to be sensitive to the gradual streamwise evolution of the mean flow and it is concluded that the effects of transverse curvature on stability should not be ignored routinely.     

Performance of a multigrid calculation procedure in three-dimensional sudden expansion flows

The performance of a recently developed calculation procedure for steady incompressible flows is assessed in a variety of three-dimensional sudden expansion type flows representative of those encountered in several types of industrial equipment. The calculation procedure, called here BLIMM (for block-implicit multigrid method), is based on a coupled solution of the three-dimensional momentum and continuity equations in primitive variables, using the multigrid technique. Different Reynolds numbers and finite difference grids are considered for each flow situation. The rates of convergence and the computational times are reported for each case.     

AN EXPLICIT FINITE VOLUME SPATIAL MARCHING METHOD FOR REDUCED NAVIER–STOKES EQUATIONS

This paper develops a spatial marching method for high-speed flows based on a finite volume approach. The method employs the reduced Navier– Stokes equations and a pressure splitting in the streamwise direction based on the Vigneron strategy. For marching from an upstream station to one downstream the modified five-level Runge–Kutta integration scheme due to Jameson and Schmidt is used. In addition, for shock handling and for good convergence properties the method employs a matrix form of the artificial dissipation terms, which has been shown to improve the accuracy of predictions. To achieve a fast rate of convergence, a local time-stepping concept is used. The method retains the time derivative in the governing equations and the solution at every spatial station is obtained in an iterative manner. The developed method is validated against two test cases: (a) supersonic flow past a flat plate; and (b) hypersonic flow past a compression corner involving a strong viscous–inviscid interaction. The computed wall pressure and wall heat transfer coefficients exhibit good general agreement with previous computations by other investigators and with experiments.     

Analysis of lubricated squeezing flow

A thin film of low-viscosity lubricating liquid between a solid wall and a viscous material reduces shear stress on the latter and tends to make it flow as though it were slipping along the wall. The result when the lubricated material is being squeezed out of the gap between approaching parallel plates is flow more nearly irrotational, or extensional, the more effective the lubricating film on the plates. Two Newtonian analyses of this flow situation are reported. One is an approximate, asymptotic analytical solution for Newtonian lubricating flow in the films and combined mixed flow, shear and extension, in the viscous layer. The second is a full two-dimensional axisymmetric solution of the momentum and continuity equations along with the kinematic condition which governs the motion of the interface. Both analyses indicate that there are two limiting flow regimes, depending on the ratio of the thickness of each of the two phases to radius and on the viscosity ratio of the two liquids. In one limit the flow is parallel squeezing and the lubricant layer slowly thins and persists a long time. In the other the lubricant is expelled preferentially. Implications of the results are discussed for rheological characterization of viscoelastic liquids and for prediction of lubricated or autolubricated flows in processing situations.     

A stable fast marching scheme for computational fluid mechanics

A stable fast marching scheme has been developed for the solution of coupled parabolic partial differential equations such as the Navier–Stokes equations. The scheme was developed with the aid of a stability analysis. The implementation of the method in standard alternating direction implicit algorithms is straightforward. The scheme was tested on the problem of natural convection in a square cavity. The number of iterations required for convergence was significantly reduced compared to conventional methods.     

Three-dimensional reduced Navier–Stokes solutions for subsonic separated and non-separated flows,using a global pressure relaxation procedure

The reduced Navier–Stokes and thin layer approximations to the Navier–Stokes equations are used to obtain solutions for viscous subsonic three-dimensional flows. A spatial marching method is combined with a direct sparse matrix solver to obtain successive solutions in a global relaxation process. Results have been obtained for flow fields with and without regions of flow reversal.     

Finite element,stream function–vorticity solution of steady laminar natural convection

Stream function–vorticity finite element solution of two-dimensional incompressible viscous flow and natural convection is considered. Steady state solutions of the natural convection problem have been obtained for a wide range of the two independent parameters. Use of boundary vorticity formulae or iterative satisfaction of the no-slip boundary condition is avoided by application of the finite element discretization and a displacement of the appropriate discrete equations. Solution is obtained by Newton–Raphson iteration of all equations simultaneously. The method then appears to give a steady solution whenever the flow is physically steady, but it does not give a steady solution when the flow is physically unsteady. In particular, no form of asymmetric differencing is required. The method offers a degree of economy over primitive variable formulations. Physical results are given for the square cavity convection problem. The paper also reports on earlier work in which the most commonly used boundary vorticity formula was found not to satisfy the no-slip condition, and in which segregated solution procedures were attempted with very minimal success.     

Three-dimensional surface water modelling using a mesh adaptive technique

Numerical simulation of open water flow in natural courses seems to be doomed to one- or two-dimensional numerical simulations. Investigations of flow hydrodynamics through the application of three-dimensional models actually have very few appearances in the literature. This paper discusses the development and the initial implementation of a general three-dimensional and time-dependent finite volume approach to simulate the hydrodynamics of surface water flow in rivers and lakes. The slightly modified Navier-Stokes equations, together with the continuity and the water depth equations, form the theoretical basis of the model. A body-fitted time-dependent co-ordinate system has been used in the solution process, in order to accommodate the commonly complex and irregular boundary and bathymetry of natural water courses. The proposed adaptive technique allows the mesh to follow the movement of the water boundaries, including the unsteady free-water surface. The primitive variable equations are written in conservative form in the Cartesian co-ordinate system, and the computational procedure is executed in the moveable curvilinear co-ordinate system. Special stabilizing techniques are introduced in order to eliminate the oscillating behaviour associated with the finite volume formulation. Also, a new and comprehensive approximation for the pressure forces at the faces of a control volume is presented. Finally, results of several tests demonstrate the performance of the finite volume approach coupled with the adaptive technique employed in the three-dimensional time-dependent mesh system.     

Coupled Groundwater Flow and Transport in Porous Media. A Conservative or Non-conservative Form?

A new formulation for the modeling of density coupled flow and transport in porous media is presented. This formulation is based on the development of the mass balance equation by using the conservative form. The system of equations obtained by coupling the flow and transport equations using a state equation is solved by a combination of the mixed hybrid finite element method (MHFEM) and the discontinuous finite element method (DFEM). The former is applied in order to solve the flow equation and the dispersive part of the transport equation, whilst the latter is used to solve the advective part of the transport equation. Although the advantages of the MHFEM are known (efficiency calculation of velocity field and continuity of fluxes from one element to an adjacent one), its application in a classical development form (volumetric fluxes as unknowns) leads to the non-conservative version of the mass balance equation. The associated matrix of the system of equations obtained by hybridization is positive definite but non-symmetrical. By using a new approach (mass fluxes as unknowns) the conservative form of the continuity equation is preserved and the associated matrix of the system of equations obtained by hybridization becomes symmetrical. When applied to Elder's problem involving a strong density contrast, this new approach, with a lower calculation cost, leads to similar or identical results to those found in the specialized literature. The comparison between the conservative and non-conservative formulations solved with the same MHFEM and DFEM combination emphasizes the rigor and the pertinence of this new approach. Furthermore, we show the existence of a limit refinement defining the stability of the numerical solution for Elder's problem.     

Prediction of incompressible flow separation with the approximate factorization technique

This paper describes one application of the approximate factorization technique to the solution of incompressible steady viscous flow problems in two dimensions. The velocity-pressure formulation of the Navier-Stokes equations written in curvilinear non-orthogonal co-ordinates is adopted. The continuity equation is replaced with one equation for the pressure by means of the artificial compressibility concept to obtain a system parabolic in time. The resulting equations are discretized in space with centred finite differences, and the steady state solution obtained by a time-marching ADI method requiring to solve 3 x 3 block tridiagonal linear systems. An optimized fourth-order artificial dissipation is introduced to damp the numerical instabilities of the artificial compressibility equation and ensure convergence. The resulting solver is applied to the prediction of a wide variety of internal flows, including both streamlined boundaries and sharp corners, and fast convergence and good results obtained for all the configurations investigated.     

时间推进方法在叶轮机械内部流场计算中的进展

时间推进方法目前广泛应用于叶轮机械内部流场的计算,这种方法可以划分成显式时间推进方法和隐式时间推进方法．本文简要介绍了时间推进方法用于叶轮机械内部流场计算的起源及这种方法的优点和缺点,重点记述了国内外叶轮机界应用显式时间推进方法和隐式时间推进方法计算叶轮机内部流场的最新进展,并简要介绍世界上一些著名的研究机构开发的时间推进方法计算程序所采用的数值计算方法      

A numerical‐variational procedure for laminar flow in curved square ducts

A new numerical method is presented for the solution of the Navier–Stokes and continuity equations governing the internal incompressible flows. The method denoted as the CVP method consists in the numerical solution of these equations in conjunction with three additional variational equations for the continuity, the vorticity and the pressure field, using a non‐staggered grid. The method is used for the study of the characteristics of the laminar fully developed flows in curved square ducts. Numerical results are presented for the effects of the flow parameters like the curvature, the Dean number and the stream pressure gradient on the velocity distributions, the friction factor and the appearance of a pair of vortices in addition to those of the familiar secondary flow. The accuracy of the method is discussed and the results are compared with those obtained by us, using a variation of the velocity–pressure linked equation methods denoted as the PLEM method and the results obtained by other methods. Copyright © 2004 John Wiley & Sons, Ltd.     

微可压缩模型预处理技术研究

微可压缩模型(slightly compressible model,SCM)是求解低马赫数流动的一种有效模型, SCM在求解过程中不必满足速度散度为零的条件,可直接应用时间推进方法求解得到不可压缩N-S方程的解. 深入研究了该模型的效率和精度; 为了提高收敛速度将预处理技术引入到该模型中, 推导了预处理后的控制方程和特征系统, 并构造了预处理后的通量. 通过对圆柱绕流、方腔流动、NACA0012翼型和6:1椭球的数值模拟, 一方面, 进一步展示了SCM的可行性与健壮性, 表明SCM适合于低马赫流动的数值模拟; 另一方面, 充分验证了预处理技术在微可压缩模型中的作用, 一定程度上解决了低马赫数流动的收敛问题, 并提高了求解的准确性和精度. 这为SCM应用于工程实际创造了一定的条件.      

Finite volume multigrid method of the planar contraction flow of a viscoelastic fluid

This paper reports on a numerical algorithm for the steady flow of viscoelastic fluid. The conservative and constitutive equations are solved using the finite volume method (FVM) with a hybrid scheme for the velocities and first‐order upwind approximation for the viscoelastic stress. A non‐uniform staggered grid system is used. The iterative SIMPLE algorithm is employed to relax the coupled momentum and continuity equations. The non‐linear algebraic equations over the flow domain are solved iteratively by the symmetrical coupled Gauss–Seidel (SCGS) method. In both, the full approximation storage (FAS) multigrid algorithm is used. An Oldroyd‐B fluid model was selected for the calculation. Results are reported for planar 4:1 abrupt contraction at various Weissenberg numbers. The solutions are found to be stable and smooth. The solutions show that at high Weissenberg number the domain must be long enough. The convergence of the method has been verified with grid refinement. All the calculations have been performed on a PC equipped with a Pentium III processor at 550 MHz. Copyright © 2001 John Wiley & Sons, Ltd.     

Shock capturing for steady,supersonic, two-dimensional isentropic flow

A finite difference scheme is presented for the solution of the two-dimensional equations of steady, supersonic, isentropic flow. The scheme incorporates numerical characteristic decomposition, is shock-capturing by design and incorporates space marching as a result of the assumption that the flow is wholly supersonic in at least one space dimension. Results are shown for problems involving oblique hydraulic jumps and reflection from a wall.     

AN ARTIFICIAL DISSIPATION FORMULATION FOR A SEMI-IMPLICIT,PRESSURE BASED SOLUTION SCHEME FOR VISCOUS AND INVISCID FLOWS

The formulation of artificial dissipation terms for a semi-implicit, pressure based flow solver, similar to SIMPLE type formulations, is presented and is applied to both the Euler and the Navier-Stokes equations. The formulation uses generalized coordinates and a non-staggered grid. This formulation is compared to some SIMPLE and time marching formulations. The relationship between SIMPLE and time marching formulations is discussed briefly. The artificial dissipation inherent in some commonly used semi-implicit formulations, e.g. upwind differencing, powerlaw, QUICK and pressure weighting, is investigated. The scheme used here includes these dissipation terms directly, but retains the ability to mimic previous schemes. The potential for errors introduced by the simultaneous use of artificial dissipation in the continuity equation and central differencing of convective terms, is revealed. The effect of the amount of dissipation on the accuracy of the solution and the convergence rate is quantitatively demonstrated for two-dimensional inviscid flow in a mildly curved duct, three-dimensional laminar flow in a square cross section elbow with strong secondary flows, and two-dimensional turbulent flow through a turbine nozzle. The spurious effects of artificial dissipation, particularly second order dissipation, inherent in some commonly used algorithms, is clearly shown. The effect of artificial dissipation on the convergence rate is also demonstrated. The main conclusion drawn from the results is that the minimum amount of artificial dissipation that gives the required accuracy, but also an adequate convergence rate for a particular case, has to be used. This amount of dissipation is case dependent. The direct inclusion of artificial dissipation terms provides control over the amount of dissipation used.     

Lagrangian finite element analysis applied to viscous free surface fluid flow

A new Lagrangian finite element formulation is presented for time-dependent incompressible free surface fluid flow problems described by the Navier-Stokes equations. The partial differential equations describing the continuum motion of the fluid are discretized using a Galerkin procedure in conjunction with the finite element approximation. Triangular finite elements are used to represent the dependent variables of the problem. An effective time integration procedure is introduced and provides a viable computational method for solving problems with equality of representation of the pressure and velocity fields. Its success has been attributed to the strict enforcement of the continuity constraint at every stage of the iterative process. The capabilities of the analysis procedure and the computer programs are demonstrated through the solution of several problems in viscous free surface fluid flow. Comparisons of results are presented with previous theoretical, numerical and experimental results.     

Development of a computer program for the prediction of flow in a rotating cavity

A computer program has been developed to predict laminar source-sink flow in a rotating cylindrical cavity. Although the program is based on a standard finite difference technique for recirculating flow, it incorporates two novel features. Step changes in grid size are employed to obtain sufficient resolution in the boundary layers and special treatment is given to the solution of the pressure correction equations, in the ‘SIMPLE’ algorithm, in order to improve the convergence properties of the method. Results are presented both for the flow in an infinite rotating cylindrical annulus and a finite rotating cylindrical cavity, with the inner cylindrical surface acting as a uniform source and the outer cylinder as a sink. These show good agreement with existing analytical solutions and illustrate some of the problems associated with the computation of rapidly rotating flows.