A dynamically load-balanced parallel $p$-adaptive implicit high-order flux reconstruction method for under-resolved turbulence simulation

Abstract

We present a dynamically load-balanced parallel $p$-adaptive implicit high-order flux reconstruction method for under-resolved turbulence simulation. The high-order explicit first stage, singly diagonal implicit Runge-Kutta (ESDIRK) method is employed to circumvent the restriction on the time step size. The pseudo transient continuation is coupled with the matrix-free restarted generalized minimal residual (GMRES) method to solve the nonlinear equations at each stage, except the first one, of ESDIRK. We use the spectral decay smoothness indicator as the refinement/coarsening indicator for $p$-adaptation. A dynamic load balancing technique is developed with the aid of the open-source library ParMETIS. The trivial cost, compared to implicit time stepping, of mesh repartitioning and data redistribution enables us to conduct $p$-adaptation and load balancing every time step. An isentropic vortex propagation case is employed to study the impact of element weights used in mesh repartitioning on parallel efficiency. We apply the $p$-adaptive solver for implicit large eddy simulation (ILES) of the transitional flows over a cylinder when Reynolds number (Re) is $3900$ and the SD7003 wing when Re is $60000$. Numerical experiments demonstrate that a significant reduction in the run time (up to $70\%$) and total number of solution points (up to $76\%$) can be achieved with $p$-adaptation.

References in corpus (8)

• Review of Summation-by-parts schemes for initial-boundary-value problems
• Split Form Nodal Discontinuous Galerkin Schemes with Summation-By-Parts Property for the Compressible Euler Equations
• PyFR: An Open Source Framework for Solving Advection-Diffusion Type Problems on Streaming Architectures using the Flux Reconstruction Approach
• A Comparison of Hybridized and Standard DG Methods for Target-Based hp-Adaptive Simulation of Compressible Flow
• Efficient mesh management in Firedrake using PETSc-DMPlex
• A locally p-adaptive approach for Large Eddy Simulation of compressible flows in a DG framework
• An Implicit High-Order Preconditioned Flux Reconstruction Method for Low-Mach-Number Flow Simulation with Dynamic Meshes
• An implicit $P$-multigrid flux reconstruction method for simulation of locally preconditioned unsteady Navier-Stokes equations at low Mach numbers

Cited by in corpus (1)

• Dynamical p-adaptivity for LES of compressible flows in a high order DG framework