A two‐step monolithic method for the efficient simulation of incompressible flows

We propose a simple technique for improving computationally the efficiency of monolithic velocity–pressure solvers for incompressible flow problems. The idea consists in solving the discrete nonlinear system of governing equations in two steps: introducing ‘artificial’ compressibility first and afterwards correcting the solution by solving the original incompressible system. The speed‐up is obtained because of a better conditioning of the modified discrete system solved at the prediction step. The formulation can be easily implemented into existing monolithic codes requiring minor modification only. The paper concludes with two examples validating the formulation and facilitating the estimation of the obtained speed‐up. For the tests chosen, an average speed‐up is approximately double, suggesting that the method is a feasible approach for incompressible flows' simulation. Copyright © 2014 John Wiley & Sons, Ltd.     

A gas‐kinetic scheme for the simulation of turbulent flows on unstructured meshes

This paper presents a numerical method for simulating turbulent flows via coupling the Boltzmann BGK equation with Spalart–Allmaras one equation turbulence model. Both the Boltzmann BGK equation and the turbulence model equation are carried out using the finite volume method on unstructured meshes, which is different from previous works on structured grid. The application of the gas‐kinetic scheme is extended to the simulation of turbulent flows with arbitrary geometries. The adaptive mesh refinement technique is also adopted to reduce the computational cost and improve the efficiency of meshes. To organize the unstructured mesh data structure efficiently, a non‐manifold hybrid mesh data structure is extended for polygonal cells. Numerical experiments are performed on incompressible flow over a smooth flat plate and compressible turbulent flows around a NACA 0012 airfoil using unstructured hybrid meshes. These numerical results are found to be in good agreement with experimental data and/or other numerical solutions, demonstrating the applicability of the proposed method to simulate both subsonic and transonic turbulent flows. Copyright © 2016 John Wiley & Sons, Ltd.     

An arbitrary Lagrangian–Eulerian method based on the adaptive Riemann solvers for general equations of state

Approximate or exact Riemann solvers play a key role in Godunov‐type methods. In this paper, three approximate Riemann solvers, the MFCAV, DKWZ and weak wave approximation method schemes, are investigated through numerical experiments, and their numerical features, such as the resolution for shock and contact waves, are analyzed and compared. Based on the analysis, two new adaptive Riemann solvers for general equations of state are proposed, which can resolve both shock and contact waves well. As a result, an ALE method based on the adaptive Riemann solvers is formulated. A number of numerical experiments show good performance of the adaptive solvers in resolving both shock waves and contact discontinuities. Copyright © 2008 John Wiley & Sons, Ltd.