A Mach‐uniform preconditioner for incompressible and subsonic flows

In this study, a novel Mach‐uniform preconditioning method is developed for the solution of Euler equations at low subsonic and incompressible flow conditions. In contrast to the methods developed earlier in which the conservation of mass equation is preconditioned, in the present method, the conservation of energy equation is preconditioned, which enforces the divergence free constraint on the velocity field even at the limiting case of incompressible, zero Mach number flows. Despite most preconditioners, the proposed Mach‐uniform preconditioning method does not have a singularity point at zero Mach number. The preconditioned system of equations preserves the strong conservation form of Euler equations for compressible flows and recovers the artificial compressibility equations in the case of zero Mach number. A two‐dimensional Euler solver is developed for validation and performance evaluation of the present formulation for a wide range of Mach number flows. The validation cases studied show the convergence acceleration, stability, and accuracy of the present Mach‐uniform preconditioner in comparison to the non‐preconditioned compressible flow solutions. The convergence acceleration obtained with the present formulation is similar to those of the well‐known preconditioned system of equations for low subsonic flows and to those of the artificial compressibility method for incompressible flows. Copyright © 2013 John Wiley & Sons, Ltd.     

Finite volume scheme for the solution of fluid flow problems on unstructured non‐staggered grids

A recently developed non‐staggered methodology which uses the principle of applying fourth‐order dissipation to the governing pressure‐correction equation is developed so it can be applied to unstructured grids. A finite volume methodology is used for discretization. The fourth‐order dissipation term is found using second‐order gradient operators. This makes it straightforward to incorporate the dissipation term on unstructured grids. The new methodology is compared with solutions from a standard finite volume second‐order flow solver and is also tested for a standard laminar driven‐lid flow problem with grids systems that do not have a uniform structure. Finally, we demonstrate how the new methodology can be used to predict flow over a wavy boundary. Copyright © 2002 John Wiley & Sons, Ltd.     

A pressure‐based unstructured grid method for all‐speed flows

An all‐speed algorithm based on the SIMPLE pressure‐correction scheme and the ‘retarded‐density’ approach has been formulated and implemented within an unstructured grid, finite volume (FV) scheme for both incompressible and compressible flows, the latter involving interaction of shock waves. The collocated storage arrangement for all variables is adopted, and the checkerboard oscillations are eliminated by using a pressure‐weighted interpolation method, similar to that of Rhie and Chow [Numerical study of the turbulent flow past an airfoil with trailing edge separation. AIAA Journal 1983; 21 : 1525]. The solution accuracy is greatly enhanced when a higher‐order convection scheme combined with adaptive mesh refinement (AMR) are used. Copyright © 2000 John Wiley & Sons, Ltd.     

High‐order ENO and WENO schemes for unstructured grids

This work describes the implementation and analysis of high‐order accurate schemes applied to high‐speed flows on unstructured grids. The class of essentially non‐oscillatory schemes (ENO), that includes weighted ENO schemes (WENO), is discussed in the paper with regard to the implementation of third‐ and fourth‐order accurate methods. The entire reconstruction process of ENO and WENO schemes is described with emphasis on the stencil selection algorithms. The stencils can be composed by control volumes with any number of edges, e.g. triangles, quadrilaterals and hybrid meshes. In the paper, ENO and WENO schemes are implemented for the solution of the dimensionless, 2‐D Euler equations in a cell centred finite volume context. High‐order flux integration is achieved using Gaussian quadratures. An approximate Riemann solver is used to evaluate the fluxes on the interfaces of the control volumes and a TVD Runge–Kutta scheme provides the time integration of the equations. Such a coupling of all these numerical tools, together with the high‐order interpolation of primitive variables provided by ENO and WENO schemes, leads to the desired order of accuracy expected in the solutions. An adaptive mesh refinement technique provides better resolution in regions with strong flowfield gradients. Results for high‐speed flow simulations are presented with the objective of assessing the implemented capability. Copyright © 2007 John Wiley & Sons, Ltd.     

A 3D second‐order accurate projection‐based Finite Volume code on non‐staggered,non‐uniform structured grids with continuity preserving properties: application to buoyancy‐driven flows

It is well known that exact projection methods (EPM) on non‐staggered grids suffer for the presence of non‐solenoidal spurious modes. Hence, a formulation for simulating time‐dependent incompressible flows while allowing the discrete continuity equation to be satisfied up to machine‐accuracy, by using a Finite Volume‐based second‐order accurate projection method on non‐staggered and non‐uniform 3D grids, is illustrated. The procedure exploits the Helmholtz–Hodge decomposition theorem for deriving an additional velocity field that enforces the discrete continuity without altering the vorticity field. This is accomplished by first solving an elliptic equation on a compact stencil that is by performing a standard approximate projection method (APM). In such a way, three sets of divergence‐free normal‐to‐face velocities can be computed. Then, a second elliptic equation for a scalar field is derived by prescribing that its additional discrete gradient ensures the continuity constraint based on the adopted linear interpolation of the velocity. Characteristics of the double projection method (DPM) are illustrated in details and stability and accuracy of the method are addressed. The resulting numerical scheme is then applied to laminar buoyancy‐driven flows and is proved to be stable and efficient. Copyright © 2006 John Wiley & Sons, Ltd.     

A non‐diffusive,divergence‐free,finite volume‐based double projection method on non‐staggered grids

Second‐order accurate projection methods for simulating time‐dependent incompressible flows on cell‐centred grids substantially belong to the class either of exact or approximate projections. In the exact method, the continuity constraint can be satisfied to machine‐accuracy but the divergence and Laplacian operators show a four‐dimension nullspace therefore spurious oscillating solutions can be introduced. In the approximate method, the continuity constraint is relaxed, the continuity equation being satisfied up to the magnitude of the local truncation error, but the compact Laplacian operator has only the constant mode. An original formulation for allowing the discrete continuity equation to be satisfied to machine‐accuracy, while using a finite volume based projection method, is illustrated. The procedure exploits the Helmholtz–Hodge decomposition theorem for deriving an additional velocity field that enforces the discrete continuity without altering the vorticity field. This is accomplished by solving a second elliptic field for a scalar field obtained by prescribing that its additional discrete gradients ensure discrete continuity based on the previously adopted linear interpolation of the velocity. The resulting numerical scheme is applied to several flow problems and is proved to be accurate, stable and efficient. This paper has to be considered as the companion of: 'F. M. Denaro, A 3D second‐order accurate projection‐based finite volume code on non‐staggered, non‐uniform structured grids with continuity preserving properties: application to buoyancy‐driven flows. IJNMF 2006; 52 (4):393–432. Now, we illustrate the details and the rigorous theoretical framework. Copyright © 2006 John Wiley & Sons, Ltd.     

A deconvolution‐based fourth‐order finite volume method for incompressible flows on non‐uniform grids

This paper is concerned with the development of a new high‐order finite volume method for the numerical simulation of highly convective unsteady incompressible flows on non‐uniform grids. Specifically, both a high‐order fluxes integration and the implicit deconvolution of the volume‐averaged field are considered. This way, the numerical solution effectively stands for a fourth‐order approximation of the point‐wise one. Moreover, the procedure is developed in the framework of a projection method for the pressure–velocity decoupling, while originally deriving proper high‐order intermediate boundary conditions. The entire numerical procedure is discussed in detail, giving particular attention to the consistent discretization of the deconvolution operation. The present method is also cast in the framework of approximate deconvolution modelling for large‐eddy simulation. The overall high accuracy of the method, both in time and space, is demonstrated. Finally, as a model of real flow computation, a two‐dimensional time‐evolving mixing layer is simulated, with and without sub‐grid scales modelling. Copyright © 2003 John Wiley & Sons, Ltd.     

A staggered control volume scheme for unstructured triangular grids

The purpose of this work is to introduce and validate a new staggered control volume method for the simulation of 2D/axisymmetric incompressible flows. The present study introduces a numerical procedure for solving the Navier–Stokes equations using the primitive variable formulation. The proposed method is an extension of the staggered grid methodology to unstructured triangular meshes for a control volume approach which features ease of handling of irregularly shaped domains. Two alternative elements are studied: transported scalars are stored either at the sides of an element or at its vertices, while the pressure is always stored at the centre of an element. Two interpolation functions were investigated for the integration of the momentum equations: a skewed mass-weighted upwind function and a flow-oriented exponential shape function. The momentum equations are solved over the covolume of a side or of a vertex and the pressure–velocity coupling makes use of a localized linear reconstruction of the discontinuous pressure field surrounding an element in order to obtain the pressure gradient terms. The pressure equation is obtained through a discretization of the continuity equation which uses the triangular element itself as the control volume. The method is applied to the simulation of the following test cases: backward-facing step flow, flow over a two-dimensional obstacle and flow in a pipe with sudden contraction of cross-sectional area. All numerical investigations are compared with experimental data from the literature. A grid convergence and error analysis study is also carried out for flow in a driven cavity. Results compared favourably with experimental data and so the new control volume scheme is deemed well suited for the prediction of incompressible flows in complex geometries. © 1997 John Wiley & Sons, Ltd.     

Oscillations in high‐order finite difference solutions of stiff problems on non‐uniform grids

This work investigates the mitigation and elimination of scheme‐related oscillations generated in compact and classical fourth‐order finite difference solutions of stiff problems, represented here by the Burgers and Reynolds equations. The regions where severe gradients are anticipated are refined by the use of subdomains where the grid is distributed according to a geometric progression. It is observed that, for multi‐domain solutions, both the classical and compact fourth‐order finite difference schemes can exhibit spurious oscillations. When present, the oscillations are initially generated around the interface between the uniform and non‐uniform grid subdomains. Based on a thorough study of the grid distribution effects, it is shown that the numerical oscillations are caused by inadequate geometric progression ratios within the non‐uniformly discretized subdomains. Indeed, accurate solutions are obtainable if and only if the grid ratios in the non‐uniform subdomains are greater than a critical threshold ratio. It is concluded that high‐order classical and compact schemes can be used with confidence to efficiently solve one‐ or two‐dimensional problems whose solutions exhibit sharp gradients in very thin regions, provided that the numerically generated oscillations are eliminated by an appropriate choice of grid distribution within the non‐uniformly discretized subdomains. Copyright © 1999 John Wiley & Sons, Ltd.     

A 3D staggered Lagrangian scheme for ideal magnetohydrodynamics on unstructured meshes

In this paper, we propose a 3D staggered Lagrangian scheme for the ideal magnetohydrodynamics (MHD) on unstructured meshes. All the thermal variables and the magnetic induction are defined in the cell centers while the fluid velocity is located at the nodes. The meshes are compatibly discretized to ensure the geometric conservation laws in Lagrangian computation by the classical subcell method, then the momentum equation is discretized using the subcell forces and the specific internal energy equation is obtained by the total energy conservation. Invoking the Galilean invariance, magnetic flux conservation, and the thermodynamic consistency, the expressions of subcell force as well as the cell-centered velocity are derived. Besides, the magnetic divergence-free constraint is fulfilled by a projection method after each time step. Various numerical tests are presented to assert the robustness and accuracy of our scheme.     

An eigenvector‐based linear reconstruction scheme for the shallow‐water equations on two‐dimensional unstructured meshes

This paper presents a new approach to MUSCL reconstruction for solving the shallow‐water equations on two‐dimensional unstructured meshes. The approach takes advantage of the particular structure of the shallow‐water equations. Indeed, their hyperbolic nature allows the flow variables to be expressed as a linear combination of the eigenvectors of the system. The particularity of the shallow‐water equations is that the coefficients of this combination only depend upon the water depth. Reconstructing only the water depth with second‐order accuracy and using only a first‐order reconstruction for the flow velocity proves to be as accurate as the classical MUSCL approach. The method also appears to be more robust in cases with very strong depth gradients such as the propagation of a wave on a dry bed. Since only one reconstruction is needed (against three reconstructions in the MUSCL approach) the EVR method is shown to be 1.4–5 times as fast as the classical MUSCL scheme, depending on the computational application. Copyright © 2006 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 Godunov‐type fractional semi‐implicit method based on staggered grid for dam‐break modeling

A semi‐implicit finite volume model based upon staggered grid is presented for solving shallow water equation. The model employs a time‐splitting scheme that uses a predictor–corrector method for the advection term. The fluxes are calculated based on a Riemann solver in the prediction step and a downwind scheme in the correction step. A simple TVD scheme is employed for shock capturing purposes in which the Minmond limiter is used for flux functions. As a consequence of using staggered grid, an ADI method is adopted for solving the discretized equations for 2‐D problems. Several 1‐D and 2‐D flows have been modeled with satisfactory results when compared with analytical and experimental test cases. The model is also capable of simulating supercritical as well as subcritical flow. Copyright © 2010 John Wiley & Sons, Ltd.     

Running unstructured grid‐based CFD solvers on modern graphics hardware

Techniques used to implement an unstructured grid solver on modern graphics hardware are described. The three‐dimensional Euler equations for inviscid, compressible flow are considered. Effective memory bandwidth is improved by reducing total global memory access and overlapping redundant computation, as well as using an appropriate numbering scheme and data layout. The applicability of per‐block shared memory is also considered. The performance of the solver is demonstrated on two benchmark cases: a NACA0012 wing and a missile. For a variety of mesh sizes, an average speed‐up factor of roughly 9.5 × is observed over the equivalent parallelized OpenMP code running on a quad‐core CPU, and roughly 33 × over the equivalent code running in serial. Copyright © 2010 John Wiley & Sons, Ltd.     

非结构重叠网格方法研究及应用

航空、航天和兵器技术等领域的研究中存在大量包含运动边界的流场。非结构重叠网格方法是一种高效的处理动边界问题的新方法。围绕相对运动的每个物体单独生成非结构网格,在网格重叠区域通过搜索和插值完成网格系之间的信息传递,提出了动态八叉树搜索算法,发展了绝对坐标系和相对坐标系相结合的流场求解方式,采用二阶精度Van Leer/Hanel格式和四阶Runge-Kutta法分别进行空间和时间离散,形成了一种新的非结构重叠网格算法。对三维Riemann问题的求解结果与精确解能很好吻合,证明了本文的重叠网格算法具有较好的时空离散精度和插值精度。对7.62mm步枪射击过程进行了数值模拟,描述了弹丸离开膛口后膛口流场的发展过程,与实验结果体现的发展过程较为吻合,验证了本文提出的非结构网格算法体系具有较好的计算性能,是研究含动边界复杂流场的一种有效手段。     

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.     

A sliding interface method for unsteady unstructured flow simulations

The objective of this work is to develop a sliding interface method for simulations involving relative grid motion that is fast and efficient and involves no grid deformation, remeshing, or hole cutting. The method is implemented into a parallel, node‐centred finite volume, unstructured viscous flow solver. The rotational motion is accomplished by rigidly rotating the subdomain representing the moving component. At the subdomain interface boundary, the faces along the interface are extruded into the adjacent subdomain to create new volume elements forming a one‐cell overlap. These new volume elements are used to compute a flux across the subdomain interface. An interface flux is computed independently for each subdomain. The values of the solution variables and other quantities for the nodes created by the extrusion process are determined by linear interpolation. The extrusion is done so that the interpolation will maintain information as localized as possible. The grid on the interface surface is arbitrary. The boundary between the two subdomains is completely independent from one another; meaning that they do not have to connect in a one‐to‐one manner and no symmetry or pattern restrictions are placed on the grid. A variety of numerical simulations were performed on model problems and large‐scale applications to examine conservation of the interface flux. Overall solution errors were found to be comparable to that for fully connected and fully conservative simulations. Excellent agreement is obtained with theoretical results and results from other solution methodologies. Copyright © 2006 John Wiley & Sons, Ltd.     

A new r‐ratio formulation for TVD schemes for vertex‐centered FVM on an unstructured mesh

The r‐ratio is a parameter that measures the local monotonicity, by which a number of high‐resolution and TVD schemes can be formed. A number of r‐ratio formulations for TVD schemes have been presented over the last few decades to solve the transport equation in shallow waters based on the finite volume method (FVM). However, unlike structured meshes, the coordinate directions are not clearly defined on an unstructured mesh; therefore, some r‐ratio formulations have been established by approximating the solute concentration at virtual nodes, which may be estimated from different assumptions. However, some formulations may introduce either oscillation or diffusion behavior within the vertex‐centered (VC) framework. In this paper, a new r‐ratio formulation, applied to an unstructured grid in the VC framework, is proposed and compared with the traditional r‐ratio formulations. Through seven commonly used benchmark tests, it is shown that the newly proposed r‐ratio formulation obtains better results than the traditional ones with less numerical diffusion and spurious oscillation. Moreover, three commonly used TVD schemes—SUPERBEE, MINMOD, and MUSCL—and two high‐order schemes—SOU and QUICK—are implemented and compared using the new r‐ratio formulation. The new r‐ratio formulation is shown to be sufficiently comprehensive to permit the general implementation of a high‐resolution scheme within the VC framework. Finally, the sensitivity test for different grid types demonstrates the good adaptability of this new r‐ratio formulation. Copyright © 2015 John Wiley & Sons, Ltd.     

A robust algorithm for computing fluid flows on highly non‐smooth staggered grids

A new method for computing the fluid flow in complex geometries using highly non‐smooth and non‐orthogonal staggered grid is presented. In a context of the SIMPLE algorithm, pressure and physical tangential velocity components are used as dependent variables in momentum equations. To reduce the sensitivity of the curvature terms in response to coordinate line orientation change, these terms are exclusively computed using Cartesian velocity components in momentum equations. The method is then used to solve some fairly complicated 2‐D and 3‐D flow field using highly non‐smooth grids. The accuracy of results on rough grids (with sharp grid line orientation change and non‐uniformity) was found to be high and the agreement with previous experimental and numerical results was quite good. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical analysis of conservative unstructured discretisations for low Mach flows

Unstructured meshes allow easily representing complex geometries and to refine in regions of interest without adding control volumes in unnecessary regions. However, numerical schemes used on unstructured grids have to be properly defined in order to minimise numerical errors. An assessment of a low Mach algorithm for laminar and turbulent flows on unstructured meshes using collocated and staggered formulations is presented. For staggered formulations using cell‐centred velocity reconstructions, the standard first‐order method is shown to be inaccurate in low Mach flows on unstructured grids. A recently proposed least squares procedure for incompressible flows is extended to the low Mach regime and shown to significantly improve the behaviour of the algorithm. Regarding collocated discretisations, the odd–even pressure decoupling is handled through a kinetic energy conserving flux interpolation scheme. This approach is shown to efficiently handle variable‐density flows. Besides, different face interpolations schemes for unstructured meshes are analysed. A kinetic energy‐preserving scheme is applied to the momentum equations, namely, the symmetry‐preserving scheme. Furthermore, a new approach to define the far‐neighbouring nodes of the quadratic upstream interpolation for convective kinematics scheme is presented and analysed. The method is suitable for both structured and unstructured grids, either uniform or not. The proposed algorithm and the spatial schemes are assessed against a function reconstruction, a differentially heated cavity and a turbulent self‐igniting diffusion flame. It is shown that the proposed algorithm accurately represents unsteady variable‐density flows. Furthermore, the quadratic upstream interpolation for convective kinematics scheme shows close to second‐order behaviour on unstructured meshes, and the symmetry‐preserving is reliably used in all computations. Copyright © 2016 John Wiley & Sons, Ltd.