An Implicit High-Order Preconditioned Flux Reconstruction Method for Low-Mach-Number Flow Simulation with Dynamic Meshes

Files

1903.06304v2.pdf (2.34 MB)

Links to Files

https://onlinelibrary.wiley.com/doi/10.1002/fld.4759?msockid=141ccadc7ebf6cb12145ddfe7fc76dc0

Permanent Link

http://hdl.handle.net/11603/13339
 https://doi.org/10.1002/fld.4759

Collections

UMBC Mechanical Engineering Department
UMBC Faculty Collection
UMBC Student Collection

Author/Creator

Author/Creator ORCID

Date

2019-06-22

Type of Work

journal articles
preprints
Text

Department

Program

Citation of Original Publication

Wang, Lai, and Meilin Yu. “An Implicit High-Order Preconditioned Flux Reconstruction Method for Low-Mach-Number Flow Simulation with Dynamic Meshes.” International Journal for Numerical Methods in Fluids 91, no. 7 (2019): 348–66. https://doi.org/10.1002/fld.4759.

Rights

This is the pre-peer reviewed version of the following article: Wang, Lai, and Meilin Yu. “An Implicit High-Order Preconditioned Flux Reconstruction Method for Low-Mach-Number Flow Simulation with Dynamic Meshes.” International Journal for Numerical Methods in Fluids 91, no. 7 (2019): 348–66. https://doi.org/10.1002/fld.4759., which has been published in final form at https://doi.org/10.1002/fld.4759. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Use of Self-Archived Versions.

Subjects

Flux Reconstruction/Correction Procedure via Reconstruction
Low-Mach- Number Preconditioning
High-Order Methods
Dual-Time Stepping
Unsteady Flows
Dynamic Meshes
Weak Compressibility

Abstract

A fully implicit high-order preconditioned ux reconstruction/correction procedure via reconstruction (FR/CPR) method is developed to solve the compressible Navier{Stokes equations at low Mach numbers. A dual-time stepping approach with the second-order backward di erentiation formula (BDF2) is employed to ensure temporal accuracy for unsteady ow simulation. When dynamic meshes are used to handle moving/deforming domains, the geometric conservation law (GCL) is implicitly enforced to eliminate errors due to the resolution discrepancy between BDF2 and the spatial FR/CPR discretization. The large linear system resulting from the spatial and temporal discretizations is tackled with the restarted Generalized Minimal Residual (GMRES) solver in the PETSc (Portable, Extensible Toolkit for Scienti c Computation) library. Through several benchmark steady and unsteady numerical tests, the preconditioned FR/CPR methods have demonstrated good convergence and accuracy for simulating ows at low Mach numbers. The new ow solver is then used to study the e ects of Mach number on unsteady force generation over a plunging airfoil when operating in low-Mach-number ows. It is observed that weak compressibility has a signi cant impact on thrust generation but a negligible e ect on lift generation of an oscillating airfoil.