A Very High-Order Accurate Staggered Finite Volume Scheme for the Stationary Incompressible Navier–Stokes and Euler Equations on Unstructured Meshes

International audienceWe propose a sixth-order staggered finite volume scheme based on polynomial reconstructions to achieve high accurate numerical solutions for the incompressible Navier-Stokes and Euler equations. The scheme is equipped with a fixed-point algorithm with solution relaxation to speed-up the convergence and reduce the computation time. Numerical tests are provided to assess the effectiveness of the method to achieve up to sixth-order con-2 Ricardo Costa et al. vergence rates. Simulations for the benchmark lid-driven cavity problem are also provided to highlight the benefit of the proposed high-order scheme

6th-order finite volume approximations for the stokes equations with a curved boundary

A new solver for the Stokes equations based on the finite volume method is proposed using very accurate polynomial reconstruction to provide a 6th-order scheme. We face two main difficulties: the gradient-divergence duality where the divergence free condition will impose the pressure gradient, and on the other hand, we assume that the domain has a regular curved boundary. The last point implies that a simple approximation of the boundary using piecewise segment lines dramatically reduces the scheme accuracy to at most a second-order one. We propose a new and simple technology which enables to restore the full scheme accuracy based on a specific polynomial reconstruction only using the Gauss points of the curved boundary and does not require any geometrical transformation.Fundação para a Ciência e a Tecnologia (FCT)This research was financed by FEDER Funds through Programa Operational Fatores de Competitividade — COMPETE and by Portuguese Funds FCT — Fundação para a Ciência e a Tecnologia, within the Projects PEst-C/MAT/UI0013/2014, PTDC/MAT/121185/2010 and FCT-ANR/MAT-NAN/0122/201

Convergence analysis of a colocated finite volume scheme for the incompressible Navier-Stokes equations on general 2 or 3D meshes

We study a colocated cell centered finite volume method for the approximation of the incompressible Navier-Stokes equations posed on a 2D or 3D finite domain. The discrete unknowns are the components of the velocity and the pressures, all of them colocated at the center of the cells of a unique mesh; hence the need for a stabilization technique, which we choose of the Brezzi-Pitk\"aranta type. The scheme features two essential properties: the discrete gradient is the transposed of the divergence terms and the discrete trilinear form associated to nonlinear advective terms vanishes on discrete divergence free velocity fields. As a consequence, the scheme is proved to be unconditionally stable and convergent for the Stokes problem, the steady and the transient Navier-Stokes equations. In this latter case, for a given sequence of approximate solutions computed on meshes the size of which tends to zero, we prove, up to a subsequence, the $L^2$-convergence of the components of the velocity, and, in the steady case, the weak $L^2$-convergence of the pressure. The proof relies on the study of space and time translates of approximate solutions, which allows the application of Kolmogorov's theorem. The limit of this subsequence is then shown to be a weak solution of the Navier-Stokes equations. Numerical examples are performed to obtain numerical convergence rates in both the linear and the nonlinear case.Comment: submitted September 0

Arbitrary high order discontinuous Galerkin methods for the shallow water and incompressible Navier-Stokes equations on unstructured staggered meshes

In this work we present a new class of well-balanced, arbitrary high order accurate semi-implicit discontinuous Galerkin methods for the solution of the shallow water and incompressible Navier-Stokes equations on staggered unstructured curved meshes. Isoparametric finite elements are used to take into account curved domain boundaries. Regarding two-dimensional shallow water equations, the discrete free surface elevation is defined on a primal triangular grid, while the discrete total height and the discrete velocity field are defined on an edge-based staggered dual grid. Similarly, for the two-dimensional incompressible Navier-Stokes case, the discrete pressure is defined on the main triangular grid and the velocity field is defined on the edge-based staggered grid. While staggered meshes are state of the art in classical finite difference approximations of the incompressible Navier-Stokes equations, their use in the context of high order DG schemes is novel and still quite rare. High order (better than second order) in time can be achieved by using a space-time finite element framework, where the basis and test functions are piecewise polynomials in both space and time. Formal substitution of the discrete momentum equation on the dual grid into the discrete continuity equation on the primary grid yields a very sparse system for the scalar pressure involving only the direct neighbor elements, so that it becomes a block four-point system in 2D and a block five-point system for 3D tetrahedral meshes. The resulting linear system is conveniently solved with a matrix-free GMRES algorithm. Note that the same space-time DG scheme on a collocated grid would lead to ten non-zero blocks per element in 2D and seventeen non-zero blocks in 3D, since substituting the discrete velocity into the discrete continuity equation on a collocated mesh would involve also neighbors of neighbors. From numerical experiments we find that our linear system is well-behaved and that the GMRES method converges quickly even without the use of any preconditioner, which is a unique feature in the context of high order implicit DG schemes. A very simple and efficient Picard iteration is then used in order to derive a space-time pressure correction algorithm that achieves also high order of accuracy in time, which is in general a non-trivial task in the context of high order discretizations for the incompressible Navier-Stokes equations. The special case of high order in space low order in time allows us to recover further regularity about the main linear system for the pressure, such as the symmetry and the positive semi-definiteness in the general case. This allows us to use a very fast linear solver such as the conjugate gradient (CG) method. The flexibility and accuracy of high order space-time DG methods on curved unstructured meshes allows to discretize even complex physical domains with very coarse grids in both space and time. We will further extend the previous method to three-dimensional incompressible Navier-Stokes system using a tetrahedral main grid and a corresponding face-based hexaxedral dual grid. The resulting dual mesh consists in non-standard 5-vertex hexahedral elements that cannot be represented using tensor products of one dimensional basis functions. Indeed a modal polynomial basis will be used for the dual mesh. This new family of numerical schemes is verified by solving a series of typical numerical test problems and by comparing the obtained numerical results with available exact analytical solutions or other numerical reference data. Furthermore, the comparison with available experimental results will be presented for incompressible Navier-Stokes equations

Development of CFD codes for the numerical resolution of potential flow and the incompressible form of the Navier-Stokes equations.

The aim of this project has been to consolidate and amplify the knowledge on fluid dynamics and CFD. It will be developed and verified different numerical codes for each physical case (potential flow, convection-diffusion equation and Navier-Stokes equation).The first chapter will consist on an introduction, where it will be explained the aimand the requirements of this study, a background of numerical methods, the reportor ganisation and it will be also mentioned the previous knowledge on numerical methods.On the following three chapters it will be solved and analysed three different cases of numerical approaches to fluid mechanics (non-viscid flows, convection-diffusionequations and, finally, the incompressible flow of Navier-Stokes equations). In each chapter it will be an introduction, a problem definition, a methodology of resolution (where it will be explained the procedure and the algorithm of the code), an analysis ofthe results (numerical and physical results) and finally a conclusion.On the fifth chapter it will be presented the budget, the task planning and the environmental impact and on the sixth chapter will consist on a conclusion and recommendations for a future work

Institute for Computational Mechanics in Propulsion (ICOMP)

The Institute for Computational Mechanics in Propulsion (ICOMP) is a combined activity of Case Western Reserve University, Ohio Aerospace Institute (OAI) and NASA Lewis. The purpose of ICOMP is to develop techniques to improve problem solving capabilities in all aspects of computational mechanics related to propulsion. The activities at ICOMP during 1991 are described

Study for the computational resolution of conservation equations of mass, momentum and energy. Numerical analysis and turbulence models

his undergraduate thesis presents a study on the numerical solution of the Navier-Stokes equations under various assumptions. To achieve this, several codes have been developed in C++ language to evaluate and verify the cases proposed by the Heat Transfer Technological Center (CTTC). The thesis is divided into different chapters. Firstly, there is an introduction to numerical methods, where spatial and temporal discretizations are explained, along with algorithms for solving systems of equations. Then, the main body of the work consists of four chapters. The first chapter addresses the heat conduction phenomenon and solves a transient two-dimensional conduction case. The second chapter deals with the general convection-diffusion equation and applies it to four problems to validate the code. The third chapter implements the Fractional Step Method (FSM) and solves the Lid-Driven Cavity, Differential Heated Cavity, and Square Cylinder cases. The fourth chapter solves the Burgers’ equation in the Fourier space. All four chapters include theoretical development, the algorithm’s structure for solving the case, and one or several verification cases, along with their respective reference results. Finally, there is a last chapter about the Turbulent Flow problem, where the extension to three dimensions of the numerical algorithm is explained