A fourth-order accurate adaptive solver for incompressible flow problems

We present a numerical solver for the incompressible Navier-Stokes equations that combines fourth-order-accurate discrete approximations and an adaptive tree grid (i.e. h-refinement). The scheme employs a novel compact-upwind advection scheme and a 4th-order accurate projection algorithm whereby the numerical solution exactly satisfies the incompressibility constraint. Further, we introduce a new refinement indicator that is tailored to this solver. We show tests and examples to illustrate the consistency, convergence rate and the application for the adaptive solver. The combination of the solver scheme and the proposed grid adaptation algorithm result in fourth-order convergence rates whilst only tuning a single grid-refinement parameter. The speed performance is benchmarked against a well-established second-order accurate adaptive solver alternative. We conclude that the present 4th order solver is an efficient design for problems with strong localization in the spatial-temporal domain and where a high degree of convergence for the solution statistics is desired. (C) 2022 The Author(s). Published by Elsevier Inc.

Automated summary

This article presents a numerical solver for the incompressible Navier-Stokes equations. The solver combines fourth-order-accurate discrete approximations and an adaptive tree grid to achieve high accuracy and efficiency. The solver employs a novel compact-upwind advection scheme and a 4th-order accurate projection algorithm to satisfy the incompressibility constraint. The paper also introduces a new refinement indicator tailored to this solver and demonstrates the consistency and convergence rate of the adaptive solver through tests and examples.