A high order Discontinuous Galerkin – Fourier incompressible 3D Navier–Stokes solver with rotating sliding meshes

• Published in: Journal of computational physics 2012-08, Vol.231 (21), p.7037-7056
• Main Authors: Ferrer, Esteban, Willden, Richard H.J.
• Format: Article
• Language: English
• Subjects: Airfoils; Arbitrary Lagrangian–Eulerian; Computational techniques; Curved boundaries; Exact sciences and technology; Finite element method; Fourier extension; Galerkin methods; High order Discontinuous Galerkin; Incompressible Navier–Stokes; Mathematical analysis; Mathematical methods in physics; Modal basis functions; Navier-Stokes equations; Physics; Rotating; Rotating sliding mesh; Sliding; Solvers; Symmetric Interior Penalty Galerkin
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2012.04.039

We present the development of a sliding mesh capability for an unsteady high order (order⩾3) h/p Discontinuous Galerkin solver for the three-dimensional incompressible Navier–Stokes equations. A high order sliding mesh method is developed and implemented for flow simulation with relative rotational motion of an inner mesh with respect to an outer static mesh, through the use of curved boundary elements and mixed triangular–quadrilateral meshes. A second order stiffly stable method is used to discretise in time the Arbitrary Lagrangian–Eulerian form of the incompressible Navier–Stokes equations. Spatial discretisation is provided by the Symmetric Interior Penalty Galerkin formulation with modal basis functions in the x–y plane, allowing hanging nodes and sliding meshes without the requirement to use mortar type techniques. Spatial discretisation in the z-direction is provided by a purely spectral method that uses Fourier series and allows computation of spanwise periodic three-dimensional flows. The developed solver is shown to provide high order solutions, second order in time convergence rates and spectral convergence when solving the incompressible Navier–Stokes equations on meshes where fixed and rotating elements coexist. In addition, an exact implementation of the no-slip boundary condition is included for curved edges; circular arcs and NACA 4-digit airfoils, where analytic expressions for the geometry are used to compute the required metrics. The solver capabilities are tested for a number of two dimensional problems governed by the incompressible Navier–Stokes equations on static and rotating meshes: the Taylor vortex problem, a static and rotating symmetric NACA0015 airfoil and flows through three bladed cross-flow turbines. In addition, three dimensional flow solutions are demonstrated for a three bladed cross-flow turbine and a circular cylinder shadowed by a pitching NACA0012 airfoil.