A transformation‐free HOC scheme for incompressible viscous flows on nonuniform polar grids

We recently proposed a transformation‐free higher‐order compact (HOC) scheme for two‐dimensional (2‐D) steady convection–diffusion equations on nonuniform Cartesian grids (Int. J. Numer. Meth. Fluids 2004; 44 :33–53). As the scheme was equipped to handle only constant coefficients for the second‐order derivatives, it could not be extended directly to curvilinear coordinates, where they invariably occur as variables. In this paper, we extend the scheme to cylindrical polar coordinates for the 2‐D convection–diffusion equations and more specifically to the 2‐D incompressible viscous flows governed by the Navier–Stokes (N–S) equations. We first apply the formulation to a problem having analytical solution and demonstrate its fourth‐order spatial accuracy. We then apply it to the flow past an impulsively started circular cylinder problem and finally to the driven polar cavity problem. We present our numerical results and compare them with established numerical and analytical and experimental results whenever available. This new approach is seen to produce excellent comparison in all the cases. Copyright © 2009 John Wiley & Sons, Ltd.     

A finite point method for adaptive three‐dimensional compressible flow calculations

The finite point method (FPM) is a meshless technique, which is based on both, a weighted least‐squares numerical approximation on local clouds of points and a collocation technique which allows obtaining the discrete system of equations. The research work we present is part of a broader investigation into the capabilities of the FPM to deal with 3D applications concerning real compressible fluid flow problems. In the first part of this work, the upwind‐biased scheme employed for solving the flow equations is described. Secondly, with the aim of exploiting the meshless capabilities, an h‐adaptive methodology for 2D and 3D compressible flow calculations is developed. This adaptive technique applies a solution‐based indicator in order to identify local clouds where new points should be inserted in or existing points could be safely removed from the computational domain. The flow solver and the adaptive procedure have been evaluated and the results are encouraging. Several numerical examples are provided in order to illustrate the good performance of the numerical methods presented. Copyright © 2008 John Wiley & Sons, Ltd.     

A local domain‐free discretization method for simulation of incompressible flows over moving bodies

This paper presents a local domain‐free discretization (DFD) method for the simulation of unsteady flows over moving bodies governed by the incompressible Navier–Stokes equations. The discretization strategy of DFD is that the discrete form of partial differential equations at an interior point may involve some points outside the solution domain. All the mesh points are classified as interior points, exterior dependent points and exterior independent points. The functional values at the exterior dependent points are updated at each time step by the approximate form of solution near the boundary. When the body is moving, only the status of points is changed and the mesh can stay fixed. The issue of ‘freshly cleared nodes/cells’ encountered in usual sharp interface methods does not pose any particular difficulty in the presented method. The Galerkin finite‐element approximation is used for spatial discretization, and the discrete equations are integrated in time via a dual‐time‐stepping scheme based on artificial compressibility. In order to validate the present method for moving‐boundary flow problems, two groups of flow phenomena have been simulated: (1) flows over a fixed circular cylinder, a harmonic in‐line oscillating cylinder in fluid at rest and a transversely oscillating cylinder in uniform flow; (2) flows over a pure pitching airfoil, a heaving–pitching airfoil and a deforming airfoil. The predictions show good agreement with the published numerical results or experimental data. Copyright © 2010 John Wiley & Sons, Ltd.     

A high-order compact scheme for solving the 2D steady incompressible Navier-Stokes equations in general curvilinear coordinates

In this paper, a high-order compact finite difference algorithm is established for the stream function-velocity formulation of the two-dimensional steady incompressible Navier-Stokes equations in general curvilinear coordinates. Different from the previous work, not only the stream function and its first-order partial derivatives but also the second-order mixed partial derivative is treated as unknown variable in this work. Numerical examples, including a test problem with an analytical solution, three types of lid-driven cavity flow problems with unusual shapes and steady flow past a circular cylinder as well as an elliptic cylinder with angle of attack, are solved numerically by the newly proposed scheme. For two types of the lid-driven trapezoidal cavity flow, we provide the detailed data using the fine grid sizes, which can be considered the benchmark solutions. The results obtained prove that the present numerical method has the ability to solve the incompressible flow for complex geometry in engineering applications, especially by using a nonorthogonal coordinate transformation, with high accuracy.     

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.     

Improved CVP scheme for laminar incompressible flows

An improved scheme of the continuity vorticity pressure (CVP) variational equations method is presented. The changes from the original version of the CVP method concern the velocity and the pressure correction equations that are used in the solution procedure and the topology of the grid where the method is applied. The improved CVP scheme is faster, simpler and more stable than the original version of the method. The efficiency and the accuracy of the new scheme are tested and validated through comparison of predictions and of computational time, with numerical results obtained with the SIMPLE method. Moreover, we present extensive comparisons of the results of the improved CVP scheme with numerical and experimental data from various researchers that show excellent agreement for a wide range of benchmark 2D and 3D laminar internal flow problems such as flow over a backward facing step, flow in square, circular and elliptical curved ducts and pulsating flow. Copyright © 2010 John Wiley & Sons, Ltd.     

The (9,5) HOC formulation for the transient Navier–Stokes equations in primitive variable

We recently proposed an improved (9,5) higher order compact (HOC) scheme for the unsteady two‐dimensional (2‐D) convection–diffusion equations. Because of using only five points at the current time level in the discretization procedure, the scheme was seen to be computationally more efficient than its predecessors. It was also seen to capture very accurately the solution of the unsteady 2‐D Navier–Stokes (N–S) equations for incompressible viscous flows in the stream function–vorticity (ψ – ω) formulation. In this paper, we extend the scope of the scheme for solving the unsteady incompressible N–S equations based on primitive variable formulation on a collocated grid. The parabolic momentum equations are solved for the velocity field by a time‐marching strategy and the pressure is obtained by discretizing the elliptic pressure Poisson equation by the steady‐state form of the (9,5) scheme with the Neumann boundary conditions. In particular, for pressure, we adopt a strategy on the collocated grid in conjunction with ideas borrowed from the staggered grid approach in finite volume. We first apply this extension to a problem having analytical solution and then to the famous lid‐driven square cavity problem. We also apply our formulation to the backward‐facing step problem to see how the method performs for external flow problems. The results are presented and are compared with established numerical results. This new approach is seen to produce excellent comparison in all the cases. Copyright © 2007 John Wiley & Sons, Ltd.     

Incompressible Roe‐like scheme for turbulent flows

In the current study, numerical investigation of incompressible turbulent flow is presented. By the artificial compressibility method, momentum and continuity equations are coupled. Considering Reynolds averaged Navier–Stokes equations, the Spalart–Allmaras turbulence model, which has accurate results in two‐dimensional problems, is used to calculate Reynolds stresses. For convective fluxes a Roe‐like scheme is proposed for the steady Reynolds averaged Navier–Stokes equations. Also, Jameson averaging method was implemented. In comparison, the proposed characteristics‐based upwind incompressible turbulent Roe‐like scheme, demonstrated very accurate results, high stability, and fast convergence. The fifth‐order Runge–Kutta scheme is used for time discretization. The local time stepping and implicit residual smoothing were applied as the convergence acceleration techniques. Suitable boundary conditions have been implemented considering flow behavior. The problem has been studied at high Reynolds numbers for cross flow around the horizontal circular cylinder and NACA0012 hydrofoil. Results were compared with those of others and a good agreement has been observed. Copyright © 2012 John Wiley & Sons, Ltd.     

Assessment of regularized delta functions and feedback forcing schemes for an immersed boundary method

We present an improved immersed boundary method for simulating incompressible viscous flow around an arbitrarily moving body on a fixed computational grid. To achieve a large Courant–Friedrichs–Lewy number and to transfer quantities between Eulerian and Lagrangian domains effectively, we combined the feedback forcing scheme of the virtual boundary method with Peskin's regularized delta function approach. Stability analysis of the proposed method was carried out for various types of regularized delta functions. The stability regime of the 4‐point regularized delta function was much wider than that of the 2‐point delta function. An optimum regime of the feedback forcing is suggested on the basis of the analysis of stability limits and feedback forcing gains. The proposed method was implemented in a finite‐difference and fractional‐step context. The proposed method was tested on several flow problems, including the flow past a stationary cylinder, inline oscillation of a cylinder in a quiescent fluid, and transverse oscillation of a circular cylinder in a free‐stream. The findings were in excellent agreement with previous numerical and experimental results. Copyright © 2008 John Wiley & Sons, Ltd.     

An immersed boundary solver for inviscid compressible flows

In this paper, a simple and efficient immersed boundary (IB) method is developed for the numerical simulation of inviscid compressible Euler equations. We propose a method based on coordinate transformation to calculate the unknowns of ghost points. In the present study, the body‐grid intercept points are used to build a complete bilinear (2‐D)/trilinear (3‐D) interpolation. A third‐order weighted essentially nonoscillation scheme with a new reference smoothness indicator is proposed to improve the accuracy at the extrema and discontinuity region. The dynamic blocked structured adaptive mesh is used to enhance the computational efficiency. The parallel computation with loading balance is applied to save the computational cost for 3‐D problems. Numerical tests show that the present method has second‐order overall spatial accuracy. The double Mach reflection test indicates that the present IB method gives almost identical solution as that of the boundary‐fitted method. The accuracy of the solver is further validated by subsonic and transonic flow past NACA2012 airfoil. Finally, the present IB method with adaptive mesh is validated by simulation of transonic flow past 3‐D ONERA M6 Wing. Global agreement with experimental and other numerical results are obtained.     

On the efficiency of linearization schemes and coupled multigrid methods in the simulation of a 3D flow around a cylinder

This paper studies the efficiency of two ways to treat the non‐linear convective term in the time‐dependent incompressible Navier–Stokes equations and of two multigrid approaches for solving the arising linear algebraic saddle point problems. The Navier–Stokes equations are discretized by a second‐order implicit time stepping scheme and by inf–sup stable, higher order finite elements in space. The numerical studies are performed at a 3D flow around a cylinder. Copyright © 2005 John Wiley & Sons, Ltd.     

Three‐dimensional computation for flow‐induced vibrations in in‐line and cross‐flow directions of a circular cylinder

We present numerical results for in‐line and cross‐flow vibrations of a circular cylinder, which is immersed in a uniform flow and is elastically supported by damper‐spring systems to compute vibrations of a rigid cylinder. In the case of a circular cylinder with a low Scruton number, it is well‐known that two types of self‐excited vibrations appear in the in‐line direction in the range of low reduced velocities. On the other hand, a cross‐flow vibration of the circular cylinder can be excited in the range of high reduced velocities. Therefore, we compute the flow‐induced vibrations of the circular cylinder in the wide range of the reduced velocities at low and high Scruton numbers and discuss about excitation mechanisms in the in‐line and cross‐flow directions. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulation of non-linear elastic flows with a general collocated finite-volume method

This paper reports the development and application of a finite-volume based methodology for the calculation of the flow of fluids which follow differential viscoelastic constitutive models. The novelty of the method lies on the use of the non-staggered grid arrangement, in which all dependent variables are located at the center of the control volumes, thus greatly simplifying the adoption of general curvilinear coordinates. The pressure–velocity–stress decoupling was removed by the development of a new interpolation technique inspired on that of Rhie and Chow, AIAA 82 (1982) 998. The differencing schemes are second order accurate and the resulting algebraic equations for each variable are solved in a segregated way (decoupled scheme). The numerical formulation especially designed for the interpolation of the stress field was found to work well and is shown to be indispensable for accurate results. Calculations have been carried out for two problems: the entry flow problem of Eggleton et al., J. Non-Newtonian Fluid Mech. 64 (1996) 269, with orthogonal and non-orthogonal meshes; and the bounded and unbounded flows around a circular cylinder. The results of the simulations compare favourably with those in the literature and iterative convergence has been attained for Deborah and Reynolds numbers similar to, or higher than, those reported for identical flow problems using other numerical methods. The application of the method with non-orthogonal coordinates is demonstrated. The entry flow problem is studied in more detail and for this case differences between Newtonian and viscoelastic fluids are identified and discussed. Viscoelasticity is shown to be responsible for the development of very intense normal stresses, which are tensile in the wall region. As a consequence, the viscoelastic fluid is more intensely decelerated in the wall region than the Newtonian fluid, thus reducing locally the shear rates and the role of viscosity in redeveloping the flow. A layer of high stress-gradients is formed at the wall leading edge and is convected below and away from the wall; its effect is to intensify the aforementioned deviation of elastic fluid from the wall.     

A Cartesian grid technique based on one‐dimensional integrated radial basis function networks for natural convection in concentric annuli

This paper reports a radial basis function (RBF)‐based Cartesian grid technique for the simulation of two‐dimensional buoyancy‐driven flow in concentric annuli. The continuity and momentum equations are represented in the equivalent stream function formulation that reduces the number of equations from three to one, but involves higher‐order derivatives. The present technique uses a Cartesian grid to discretize the problem domain. Along a grid line, one‐dimensional integrated RBF networks (1D‐IRBFNs) are employed to represent the field variables. The capability of 1D‐IRBFNs to handle unstructured points with accuracy is exploited to describe non‐rectangular boundaries in a Cartesian grid, while the method's ability to avoid the reduction of convergence rate caused by differentiation is instrumental in improving the quality of the approximation of higher‐order derivatives. The method is applied to simulate thermally driven flows in annuli between two circular cylinders and between an outer square cylinder and an inner circular cylinder. High Rayleigh number solutions are achieved and they are in good agreement with previously published numerical data. Copyright © 2007 John Wiley & Sons, Ltd.     

Shallow flow simulation on dynamically adaptive cut cell quadtree grids

A computationally efficient, high‐resolution numerical model of shallow flow hydrodynamics is described, based on dynamically adaptive quadtree grids. The numerical model solves the two‐dimensional non‐linear shallow water equations by means of an explicit second‐order MUSCL‐Hancock Godunov‐type finite volume scheme. Interface fluxes are evaluated using an HLLC approximate Riemann solver. Cartesian cut cells are used to improve the fit to curved boundaries. A ghost‐cell immersed boundary method is used to update flow information in the smallest cut cells and overcome the time step restriction that would otherwise apply. The numerical model is validated through simulations of reflection of a surge wave at a wall, a low Froude number potential flow past a circular cylinder, and the shock‐like interaction between a bore and a circular cylinder. The computational efficiency is shown to be greatly improved compared with solutions on a uniform structured grid implemented with cut cells. Copyright © 2006 John Wiley & Sons, Ltd.     

A stable second‐order mass‐weighted upwind scheme for unstructured meshes

In this paper, an original second‐order upwind scheme for convection terms is described and implemented in the context of a Control‐Volume Finite‐Element Method (CVFEM). The proposed scheme is a second‐order extension of the first‐order MAss‐Weighted upwind (MAW) scheme proposed by Saabas and Baliga (Numer. Heat Transfer 1994; 26B :381–407). The proposed second‐order scheme inherits the well‐known stability characteristics of the MAW scheme, but exhibits less artificial viscosity and ensures much higher accuracy. Consequently, and in contrast with nearly all second‐order upwind schemes available in the literature, the proposed second‐order MAW scheme does not need limiters. Some test cases including two pure convection problems, the driven cavity and steady and unsteady flows over a circular cylinder, have been undertaken successfully to validate the new scheme. The verification tests show that the proposed scheme exhibits a low level of artificial viscosity in the pure convection problems; exhibits second‐order accuracy for the driven cavity; gives accurate reattachment lengths for low‐Reynolds steady flow over a circular cylinder; and gives constant‐amplitude vortex shedding for the case of high‐Reynolds unsteady flow over a circular cylinder. Copyright © 2005 John Wiley & Sons, Ltd.     

Development of a hybrid finite volume/element solver for incompressible flows

In this paper, we report our development of an implicit hybrid flow solver for the incompressible Navier–Stokes equations. The methodology is based on the pressure correction or projection method. A fractional step approach is used to obtain an intermediate velocity field by solving the original momentum equations with the matrix‐free implicit cell‐centred finite volume method. The Poisson equation derived from the fractional step approach is solved by the node‐based Galerkin finite element method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centres and the auxiliary variable at cell vertices, making the current solver a staggered‐mesh scheme. Numerical examples demonstrate the performance of the resulting hybrid scheme, such as the correct temporal convergence rates for both velocity and pressure, absence of unphysical pressure boundary layer, good convergence in steady‐state simulations and capability in predicting accurate drag, lift and Strouhal number in the flow around a circular cylinder. Copyright © 2007 John Wiley & Sons, Ltd.     

幂律流体圆柱绕流的格子波尔兹曼模拟

基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明:插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中.      

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 modeling of wave interactions with coastal structures by a constrained interpolation profile/immersed boundary method

A high‐order difference method based multiphase model is proposed to simulate nonlinear interactions between water wave and submerged coastal structures. The model is based on the Navier–Stokes equations using a constrained interpolation profile (CIP) method for the flow solver, and employs an immersed boundary method (IBM) for the treatment of wave–structure interactions. A more accurate interface capturing scheme, the volume of fluid/weighed line interface calculation (VOF/WLIC) scheme, is adopted as the interface capturing method. A series of computations are performed to verify the application of the model for simulations of fluid interaction with various structures. These problems include flow over a fixed cylinder, water entry of a circular cylinder and solitary waves passing various submerged coastal structures. Computations are compared with the available analytical, experimental and other numerical results and good agreement is obtained. The results of this study demonstrate the accuracy and applications of the proposed model to simulate the nonlinear flow phenomena and capture the complex free surface flow. Copyright © 2015 John Wiley & Sons, Ltd.