Simulation of single‐ and two‐phase flows on sliding unstructured meshes using finite volume method

The paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.     

A hybrid time stepping scheme combining explicit Runge–Kutta with implicit LU‐SGS for Navier–Stokes equations

A hybrid time stepping scheme is developed and implemented by a combination of explicit Runge–Kutta with implicit LU‐SGS scheme at the level of system matrix. In this method, the explicit scheme is applied to those grid cells of blocks that have large local time steps; meanwhile, the implicit scheme is applied to other grid cells of blocks that have smaller allowable local time steps in the same flow field. As a result, the discretized governing equations can be expressed as a compound of explicit and implicit matrix operator. The proposed method has been used to compute the steady transonic turbulent flow over the RAE 2822 airfoil. The numerical results are found to be in excellent agreement with the experimental data. In the validation case, the present scheme saved at least 50% of the memory resources compared with the fully implicit LU‐SGS. Copyright © 2015 John Wiley & Sons, Ltd.     

A mass‐conserving Level‐Set method for modelling of multi‐phase flows

A mass‐conserving Level‐Set method to model bubbly flows is presented. The method can handle high density‐ratio flows with complex interface topologies, such as flows with simultaneous occurrence of bubbles and droplets. Aspects taken into account are: a sharp front (density changes abruptly), arbitrarily shaped interfaces, surface tension, buoyancy and coalescence of droplets/bubbles. Attention is paid to mass‐conservation and integrity of the interface. The proposed computational method is a Level‐Set method, where a Volume‐of‐Fluid function is used to conserve mass when the interface is advected. The aim of the method is to combine the advantages of the Level‐Set and Volume‐of‐Fluid methods without the disadvantages. The flow is computed with a pressure correction method with the Marker‐and‐Cell scheme. Interface conditions are satisfied by means of the continuous surface force methodology and the jump in the density field is maintained similar to the ghost fluid method for incompressible flows. Copyright © 2005 John Wiley & Sons, Ltd.     

A virtual boundary method with improved computational efficiency using a multi‐grid method

The flow around spherical, solid objects is considered. The boundary conditions on the solid boundaries have been applied by replacing the boundary with a surface force distribution on the surface, such that the required boundary conditions are satisfied. The velocity on the boundary is determined by extrapolation from the flow field. The source terms are determined iteratively, as part of the solution. They are then averaged and are smoothed out to nearby computational grid points. A multi‐grid scheme has been used to enhance the computational efficiency of the solution of the force equations. The method has been evaluated for flow around both moving and stationary spherical objects at very low and intermediate Reynolds numbers. The results shows a second order accuracy of the method both at creeping flow and at Re=100. The multi‐grid scheme is shown to enhance the convergence rate up to a factor 10 as compared to single grid approach. Copyright © 2004 John Wiley & Sons, Ltd.     

An overlapping control volume method for the Navier–Stokes equations on non‐staggered grids

An algorithm, based on the overlapping control volume (OCV) method, for the solution of the steady and unsteady two‐dimensional incompressible Navier–Stokes equations in complex geometry is presented. The primitive variable formulation is solved on a non‐staggered grid arrangement. The problem of pressure–velocity decoupling is circumvented by using momentum interpolation. The accuracy and effectiveness of the method is established by solving five steady state and one unsteady test problems. The numerical solutions obtained using the technique are in good agreement with the analytical and benchmark solutions available in the literature. On uniform grids, the method gives second‐order accuracy for both diffusion‐ and convection‐dominated flows. There is little loss of accuracy on grids that are moderately non‐orthogonal. Copyright © 1999 John Wiley & Sons, Ltd.     

Employment of the second‐moment turbulence closure on arbitrary unstructured grids

The paper presents a finite‐volume calculation procedure using a second‐moment turbulence closure. The proposed method is based on a collocated variable arrangement and especially adopted for unstructured grids consisting of ‘polyhedral’ calculation volumes. An inclusion of 23k in the pressure is analysed and the impact of such an approach on the employment of the constant static pressure boundary is addressed. It is shown that this approach allows a removal of a standard but cumbersome velocity–pressure –Reynolds stress coupling procedure known as an extension of Rhie‐Chow method (AIAA J. 1983; 21 : 1525–1532) for the Reynolds stresses. A novel wall treatment for the Reynolds‐stress equations and ‘polyhedral’ calculation volumes is presented. Important issues related to treatments of diffusion terms in momentum and Reynolds‐stress equations are also discussed and a new approach is proposed. Special interpolation practices implemented in a deferred‐correction fashion and related to all equations, are explained in detail. Computational results are compared with available experimental data for four very different applications: the flow in a two‐dimensional 180o turned U‐bend, the vortex shedding flow around a square cylinder, the flow around Ahmed Body and in‐cylinder engine flow. Additionally, the performance of the methodology is assessed by applying it to different computational grids. For all test cases, predictions with the second‐moment closure are compared to those of the k–εmodel. The second‐moment turbulence closure always achieves closer agreement with the measurements. A moderate increase in computing time is required for the calculations with the second‐moment closure. Copyright © 2004 John Wiley & Sons, Ltd.     

An improved unstructured WENO method for compressible multi‐fluids flow

An improved high‐order accurate WENO finite volume method based on unstructured grids for compressible multi‐fluids flow is proposed in this paper. The third‐order accuracy WENO finite volume method based on triangle cell is used to discretize the governing equations. To have higher order of accuracy, the P1 polynomial is reconstructed firstly. After that, the P2 polynomial is reconstructed from the combination of the P1. The reconstructed coefficients are calculated by analytical form of inverse matrix rather than the numerical inversion. This greatly improved the efficiency and the robustness. Four examples are presented to examine this algorithm. Numerical results show that there is no spurious oscillation of velocity and pressure across the interface and high‐order accurate result can be achieved. Copyright © 2016 John Wiley & Sons, Ltd.     

A Taylor series‐based finite volume method for the Navier–Stokes equations

A Taylor series‐based finite volume formulation has been developed to solve the Navier–Stokes equations. Within each cell, velocity and pressure are obtained from the Taylor expansion at its centre. The derivatives in the expansion are found by applying the Gauss theorem over the cell. The resultant integration over the faces of the cell is calculated from the value at the middle point of the face and its derivatives, which are further obtained from a higher order interpolation based on the values at the centres of two cells sharing this face. The terms up to second order in the velocity and the terms up to first order in pressure in the Taylor expansion are retained throughout the derivation. The test cases for channel flow, flow past a circular cylinder and flow in a collapsible channel have shown that the method is quite accurate and flexible. Copyright © 2008 John Wiley & Sons, Ltd.     

基于非结构/混合网格的高阶精度格式研究进展

尽管以二阶精度格式为基础的计算流体力学(CFD) 方法和软件已经在航空航天飞行器设计中发挥了重要的作用, 但是由于二阶精度格式的耗散和色散较大, 对于湍流、分离等多尺度流动现象的模拟, 现有成熟的CFD 软件仍难以给出满意的结果, 为此CFD 工作者发展了众多的高阶精度计算格式. 如果以适应的计算网格来分类, 一般可以分为基于结构网格的有限差分格式、基于非结构/混合网格的有限体积法和有限元方法,以及各种类型的混合方法. 由于非结构/混合网格具有良好的几何适应性, 基于非结构/混合网格的高阶精度格式近年来备受关注. 本文综述了近年来基于非结构/混合网格的高阶精度格式研究进展, 重点介绍了空间离散方法, 主要包括k-Exact 和ENO/WENO 等有限体积方法, 间断伽辽金(DG) 有限元方法, 有限谱体积(SV) 和有限谱差分(SD) 方法, 以及近来发展的各种DG/FV 混合算法和将各种方法统一在一个框架内的CPR (correctionprocedure via reconstruction) 方法等. 随后简要介绍了高阶精度格式应用于复杂外形流动数值模拟的一些需要关注的问题, 包括曲边界的处理方法、间断侦测和限制器、各种加速收敛技术等. 在综述过程中, 介绍了各种方法的优势与不足, 其间介绍了作者发展的基于"静动态混合重构" 的DG/FV 混合算法. 最后展望了基于非结构/混合网格的高阶精度格式的未来发展趋势及应用前景.     

An implicit meshless method for application in computational fluid dynamics

An implicit meshless scheme is developed for solving the Euler equations, as well as the laminar and Reynolds‐averaged Navier–Stokes equations. Spatial derivatives are approximated using a least squares method on clouds of points. The system of equations is linearised, and solved implicitly using approximate, analytical Jacobian matrices and a preconditioned Krylov subspace iterative method. The details of the spatial discretisation, linear solver and construction of the Jacobian matrix are discussed; and results that demonstrate the performance of the scheme are presented for steady and unsteady two dimensional fluid flows. Copyright © 2012 John Wiley & Sons, Ltd.     

SIMPLE‐type preconditioners for cell‐centered,colocated finite volume discretization of incompressible Reynolds‐averaged Navier–Stokes equations

This paper contains a comparison of four SIMPLE‐type methods used as solver and as preconditioner for the iterative solution of the (Reynolds‐averaged) Navier–Stokes equations, discretized with a finite volume method for cell‐centered, colocated variables on unstructured grids. A matrix‐free implementation is presented, and special attention is given to the treatment of the stabilization matrix to maintain a compact stencil suitable for unstructured grids. We find SIMPLER preconditioning to be robust and efficient for academic test cases and industrial test cases. Compared with the classical SIMPLE solver, SIMPLER preconditioning reduces the number of nonlinear iterations by a factor 5–20 and the CPU time by a factor 2–5 depending on the case. The flow around a ship hull at Reynolds number 2E9, for example, on a grid with cell aspect ratio up to 1:1E6, can be computed in 3 instead of 15 h.Copyright © 2012 John Wiley & Sons, Ltd.     

One‐shot pseudo‐time method for aerodynamic shape optimization using the Navier–Stokes equations

This paper presents a numerical method for aerodynamic shape optimization problems in compressible viscous flow. It is based on simultaneous pseudo‐time stepping in which stationary states are obtained by solving the pseudo‐stationary system of equations representing the state, costate and design equations. The main advantages of this method are that it blends in nicely with previously existing pseudo‐time‐stepping methods for the state and the costate equations, that it requires no additional globalization in the design space, and that a preconditioner can be used for convergence acceleration which stems from the reduced SQP methods. For design examples of 2D problems, the overall cost of computation can be reduced to less than 2 times the forward simulation runs. Copyright © 2011 John Wiley & Sons, Ltd.     

Parallel‐in‐time simulation of the unsteady Navier–Stokes equations for incompressible flow

In this paper, we present an application of a parallel‐in‐time algorithm for the solution of the unsteady Navier–Stokes model equations that are of parabolic–elliptic type. This method is based on the alternated use of a coarse global sequential solver and a fine local parallel one. A standard finite volume/finite differences first‐order approach is used for discretization of the unsteady two‐dimensional Navier–Stokes equations. The Taylor vortex decay problem and the confined flow around a square cylinder were selected as unsteady flow examples to illustrate and analyse the properties of the parallel‐in‐time method through numerical experiments. The influence of several parameters on the computing time required to perform a parallel‐in‐time calculation on a PC cluster was verified. Among them we have analysed the influence of the number of processors, the number of iterations for convergence, the resolution of the spatial domain and the influence of the time‐step sizes ratio between the coarse and fine grids. Significant computer time saving was achieved when compared with the single processor computing time, particularly when the spatial dimension of the problem is low and the temporal scale is large. Copyright © 2004 John Wiley & Sons, Ltd.     

A robust multi‐grid pressure‐based algorithm for multi‐fluid flow at all speeds

This paper reports on the implementation and testing, within a full non‐linear multi‐grid environment, of a new pressure‐based algorithm for the prediction of multi‐fluid flow at all speeds. The algorithm is part of the mass conservation‐based algorithms (MCBA) group in which the pressure correction equation is derived from overall mass conservation. The performance of the new method is assessed by solving a series of two‐dimensional two‐fluid flow test problems varying from turbulent low Mach number to supersonic flows, and from very low to high fluid density ratios. Solutions are generated for several grid sizes using the single grid (SG), the prolongation grid (PG), and the full non‐linear multi‐grid (FMG) methods. The main outcomes of this study are: (i) a clear demonstration of the ability of the FMG method to tackle the added non‐linearity of multi‐fluid flows, which is manifested through the performance jump observed when using the non‐linear multi‐grid approach as compared to the SG and PG methods; (ii) the extension of the FMG method to predict turbulent multi‐fluid flows at all speeds. The convergence history plots and CPU‐times presented indicate that the FMG method is far more efficient than the PG method and accelerates the convergence rate over the SG method, for the problems solved and the grids used, by a factor reaching a value as high as 15. Copyright © 2003 John Wiley & Sons, Ltd.     

Calculation of three‐dimensional turbulent flow with a finite volume multigrid method

A numerical method for the efficient calculation of three‐dimensional incompressible turbulent flow in curvilinear co‐ordinates is presented. The mathematical model consists of the Reynolds averaged Navier–Stokes equations and the k–ε turbulence model. The numerical method is based on the SIMPLE pressure‐correction algorithm with finite volume discretization in curvilinear co‐ordinates. To accelerate the convergence of the solution method a full approximation scheme‐full multigrid (FAS‐FMG) method is utilized. The solution of the k–ε transport equations is embedded in the multigrid iteration. The improved convergence characteristic of the multigrid method is demonstrated by means of several calculations of three‐dimensional flow cases. Copyright © 1999 John Wiley & Sons, Ltd.     

Two‐level finite element method with a stabilizing subgrid for the incompressible Navier–Stokes equations

We consider the Galerkin finite element method for the incompressible Navier–Stokes equations in two dimensions. The domain is discretized into a set of regular triangular elements and the finite‐dimensional spaces employed consist of piecewise continuous linear interpolants enriched with the residual‐free bubble functions. To find the bubble part of the solution, a two‐level finite element method with a stabilizing subgrid of a single node is described, and its application to the Navier–Stokes equation is displayed. Numerical approximations employing the proposed algorithm are presented for three benchmark problems. The results show that the proper choice of the subgrid node is crucial in obtaining stable and accurate numerical approximations consistent with the physical configuration of the problem at a cheap computational cost. Copyright © 2008 John Wiley & Sons, Ltd.     

A new hybrid turbulence modelling strategy for industrial CFD

This paper presents a new strategy for turbulence model employment with emphasis on the model's applicability for industrial computational fluid dynamics (CFD). In the hybrid modelling strategy proposed here, the Reynolds stress and mean rate of strain tensors are coupled via Boussinesq's formula as in the standard k–εmodel. However, the turbulent kinetic energy is calculated as the sum of the normal Reynolds‐stress components, representing the solutions of the appropriate transport equations. The equations governing the Reynolds‐stress tensor and dissipation rate have been solved in the framework of a ‘background’ second‐moment closure model. Furthermore, the structure parameter C‐µ has been re‐calculated from a newly proposed functional dependency rather than kept constant. This new definition of C‐µ has been assessed by using direct numerical simulation (DNS) results of several generic flow configurations featuring different phenomena such as separation, reattachment and rotation. Comparisons show a large departure of C‐µ from the commonly used value of 0.09. The model proposed is computationally validated in a number of well‐proven fluid flow benchmarks, e.g. backward‐facing step, 180° turn‐around duct, rotating pipe, impinging jet and three‐dimensional (3D) Ahmed body. The obtained results confirm that the present hybrid model delivers a robust solution procedure while preserving most of the physical advantages of the Reynolds‐stress model over simple k–εmodels. A low Reynolds number version of the hybrid model is also proposed and discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

A new characteristic approach for incompressible thermo‐flow in Cartesian and non‐Cartesian grids

A virtual‐characteristic approach is developed for thermo‐flow with finite‐volume methodology in which a multidimensional characteristic (MC) scheme is applied along with artificial compressibility. To obtain compatibility equations and pseudo‐characteristics, energy equation is taken into account in the MC scheme. With this inherent upwinding of convective fluxes, no artificial viscosity is required even at high Reynolds numbers. Another remarkable advantage of the MC scheme lies in its faster convergence rate with respect to the averaging scheme that is found to exhibit substantial delays in convergence. As benchmarks, forced and mixed convections in a cavity and in flow over cylinder and between parallel plates are examined for a wide range of Reynolds, Grashof, and Prandtl numbers. The MC and averaging schemes are applied for simulation purposes. Results show the better performance of the MC scheme in forced and mixed convections. Results confirm the robustness of the MC scheme in terms of accuracy and convergence. Copyright © 2015 John Wiley & Sons, Ltd.     

Prediction of unobstructed flow for co‐rotating multi‐disk drive in an enclosure

The flow structure and isotherms of hard disk drives (HDD) are investigated using a finite element method (FEM). The governing equations are based on the three‐dimensional axisymmetric Navier–Stokes partial differential equations (PDEs) with Galerkin FE formulation. Co‐rotating models are selected that include the non‐ventilated configuration within an enclosure. With various operating conditions for the disk system, the following important parameters are considered: disk number (n), rotational speed (Ro), and wall temperature. The flow structure changes rapidly when the rotational Reynolds number (Re_ϕ) is increased. The flow has a greater tendency to flow radially outwards and the swirling velocity tends to be more vertically orientated, especially for high Re_ϕ values. The isotherms only have small varying regions near the rotating axis, forming outward arcs near the wall and inward arcs near the end gap of the disk. Different from the case without the enclosure, the vorticities exist along the outer disk ends. Both the swirling velocity and isotherms indicate nearly symmetrical characteristics, as expected. A higher temperature gradient occurs near the right outer disk ends, which implies the characteristic of higher heat flux. A commercial computational fluid dynamic (CFD) code, CFX‐5, was chosen to validate the results. Copyright © 2001 John Wiley & Sons, Ltd.     

A multi‐moment transport model on cubed‐sphere grid

A conservative, single‐cell‐based semi‐Lagrangian transport model is proposed in this paper. Using multi‐moment concept, an additional moment, i.e. volume‐integrated average (VIA), is treated as the model variable besides the point value (PV) updated in the traditional semi‐Lagrangian schemes. A quadratic interpolation function is constructed based on local degrees of freedom defined within each single cell. The PV moment is advanced by the semi‐Lagrangian formulation, whereas the VIA moment is updated by a finite volume formulation to rigorously ensure the numerical conservation. The numerical fluxes are computed from the PV moments defined along the boundary edges of the control volume. The scheme is extended to the spherical geometry through the application of the cubed‐sphere grid that eliminates the polar singularity in the conventional longitude/latitude coordinates by using the quasi‐uniform grid spacing covering the whole sphere. The single‐cell‐based scheme is well suited for the treatment of the connections between different patches. A simple quasi‐monotone limiter to the PV moment is applied to suppress non‐physical oscillations. The proposed scheme has been validated via representative benchmark tests and the performance is competitive to other existing transport schemes. Copyright © 2010 John Wiley & Sons, Ltd.