An implicit gas-kinetic scheme for turbulent flow on unstructured hybrid mesh

In this study, an implicit scheme for the gas-kinetic scheme (GKS) on the unstructured hybrid mesh is proposed. The Spalart–Allmaras (SA) one equation turbulence model is incorporated into the implicit gas-kinetic scheme (IGKS) to predict the effects of turbulence. The implicit macroscopic governing equations are constructed and solved by the matrix-free lower-upper symmetric-Gauss–Seidel (LU-SGS) method. To reduce the number of cells and computational cost, the hybrid mesh is applied. A modified non-manifold hybrid mesh data(NHMD) is used for both unstructured hybrid mesh and uniform grid. Numerical investigations are performed on different 2D laminar and turbulent flows. The convergence property and the computational efficiency of the present IGKS method are investigated. Much better performance is obtained compared with the standard explicit gas-kinetic scheme. Also, our numerical results are found to be in good agreement with experiment data and other numerical solutions, demonstrating the good applicability and high efficiency of the present IGKS for the simulations of laminar and turbulent flows.     

An improved multislope MUSCL scheme for solving shallow water equations on unstructured grids

This paper describes an improved vector manipulation multislope monotone upstream-centred scheme for conservation laws (MUSCL) reconstruction for solving the shallow water equations on unstructured grids. This improved MUSCL reconstruction method includes a bigger stencil for the interpolation and saves time for determining the geometric relations compared to the original vector manipulation method, so it is computationally more efficient and straightforward to implement. Four examples involving an analytical solution, laboratory experiments and field-scale measurements are used to test the performance of the proposed scheme. It has been proven that the proposed scheme can provide comparable accuracy and higher efficiency compared to the original vector manipulation method. With the increasing of the number of cells, the advantage of the proposed scheme becomes more apparent.     

One- and two-equation turbulence models for the prediction of complex cascade flows using unstructured grids

One- and two-equation, low-Reynolds eddy-viscosity turbulence models are employed in the context of a primitive variable, finite volume, Navier-Stokes solver for unstructured grids. Through the study of the complex flow in a controlled-diffusion compressor cascade at off-design conditions, the ability of the models under consideration to predict the laminar separation bubble close to the leading edge and the boundary layer development is investigated. In order to control the unphysical growth of turbulent kinetic energy near the leading edge stagnation point, appropriate modifications to the conventional models are employed and tested. All of them improve the leading edge flow patterns and significantly affect the size of the predicted laminar separation bubble. The use of an adequately refined mesh around the airfoil, that is formed by triangles placed in a quasi-structured way, allows for the generation of grid elements of moderate aspect ratios. This helps to readily overcome any relevant problems of accuracy; a second-order upwind scheme without flux limiters or least squares approximations is successfully employed for the gradients. The test case includes quasi-3D effects by considering the streamtube thickness variation in the governing equations.     

An implicit matrix-free Discontinuous Galerkin solver for viscous and turbulent aerodynamic simulations

This paper presents some recent advancements of the computational efficiency of a Discontinuous Galerkin (DG) solver for the Navier–Stokes (NS) and Reynolds Averaged Navier Stokes (RANS) equations. The implementation and the performance of a Newton–Krylov matrix-free (MF) method is presented and compared with the matrix based (MB) counterpart. Moreover two solution strategies, developed in order to increase the solver efficiency, are discussed and experimented. Numerical results of some test cases proposed within the EU ADIGMA (Adaptive Higher-Order Variational Methods for Aerodynamic Applications in Industry) project demonstrate the capabilities of the method.     

Parallel computation of unsteady three-dimensional incompressible viscous flow using an unstructured multigrid method

The development and validation of a parallel unstructured tetrahedral non-nested multigrid (MG) method for simulation of unsteady 3D incompressible viscous flow is presented. The Navier-Stokes solver is based on the artificial compressibility method (ACM) and a higher-order characteristics-based finite-volume scheme on unstructured MG. Unsteady flow is calculated with an implicit dual time stepping scheme. The parallelization of the solver is achieved by a MG domain decomposition approach (MG-DD), using the Single Program Multiple Data (SPMD) programming paradigm. The Message-Passing Interface (MPI) Library is used for communication of data and loop arrays are decomposed using the OpenMP standard. The parallel codes using single grid and MG are used to simulate steady and unsteady incompressible viscous flows for a 3D lid-driven cavity flow for validation and performance evaluation purposes. The speedups and efficiencies obtained by both the parallel single grid and MG solvers are reasonably good for all test cases, using up to 32 processors on the SGI Origin 3400. The parallel results obtained agree well with those of serial solvers and with numerical solutions obtained by other researchers, as well as experimental measurements.     

Classical and variational multiscale LES of the flow around a circular cylinder on unstructured grids

The effects of numerical viscosity, subgrid scale (SGS) viscosity and grid resolution are investigated in LES and VMS-LES simulations of the flow around a circular cylinder at Re=3900 on unstructured grids. The separation between the largest and the smallest resolved scales in the VMS formulation is obtained through a variational projection operator and finite-volume cell agglomeration. Three different non-dynamic eddy-viscosity SGS models are used both in classical and in VMS-LES. The so-called small-small formulation is used in VMS-LES, i.e. the SGS viscosity is computed as a function of the smallest resolved scales. Two different grid resolutions are considered. It is found that, for each considered SGS model, the amount of SGS viscosity introduced in the VMS-LES formulation is significantly lower than in classical LES. This, together with the fact that in the VMS formulation the SGS viscosity only acts on the smallest resolved scales, has a strong impact on the results. However, a significant sensitivity of the results to the considered SGS model remains also in the VMS-LES formulation. Moreover, passing from classical LES to VMS-LES does not systematically lead to an improvement of the quality of the numerical predictions.     

High resolution finite volume computations on unstructured grids using solution dependent weighted least squares gradients

An improved high resolution finite volume method based on linear and quadratic variable reconstructions using solution dependent weighted least squares (SDWLS) gradients has been presented here. An extended stencil consisting of vertex-based neighbours of a cell is used in the higher order reconstructions for inviscid flux computations. A QR algorithm with Householder transformation is used to solve the weighted least squares problem. In case of Navier–Stokes equations, viscous fluxes are discretized in a central differencing manner based on the Coirier’s diamond path. A few inviscid and viscous test cases are solved in order to demonstrate the efficacy of the present method. Progressive improvements in solution accuracy are observed with the increase in the order of variable reconstructions. In most cases, results of quadratic reconstruction show significant improvements over that of linear reconstruction.     

Numerical solutions of a multi-class traffic flow model on an inhomogeneous highway using a high-resolution relaxed scheme

A high-resolution relaxed scheme which requires little information of the eigenstructure is presented for the multi-class Lighthill-Whitham-Richards(LWR) model on an inhomogeneous highway.The scheme needs only an estimate of the upper boundary of the maximum of absolute eigenvalues.It is based on incorporating an improved fifth-order weighted essentially non-oscillatory(WENO) reconstruction with relaxation approximation.The scheme benefits from the simplicity of relaxed schemes in that it requires no exact or approximate Riemann solvers and no projection along characteristic directions.The effectiveness of our method is demonstrated in several numerical examples.