Stability analysis of a singlerate and multirate predictor-corrector scheme for overlapping grids

by   Ketan Mittal, et al.

We use matrix stability analysis for a singlerate and multirate predictor-corrector scheme (PC) used to solve the incompressible Navier-Stokes equations (INSE) in overlapping grids. By simplifying the stability analysis with the unsteady heat equation in 1D, we demonstrate that, as expected, the stability of the PC scheme increases with increase in the resolution and overlap of subdomains. For singlerate timestepping, we also find that the high-order PC scheme is stable when the number of corrector iterations (Q) is odd. This difference in the stability of odd- and even-Q is novel and has not been demonstrated in the literature for overlapping grid-based methods. We address the odd-even behavior in the stability of the PC scheme by modifying the last corrector iterate, which leads to a scheme whose stability increases monotonically with Q. For multirate timestepping, we observe that the stability of the PC scheme depends on the timestep ratio (η). For η=2, even-Q is more stable than odd-Q. For η≥3, even-Q is more stable than odd-Q for a small nondimensional timestep size and the odd-even behavior vanishes as the timestep size increases. The stability analysis presented in this work gives novel insight into a high-order temporal discretization for ODEs and PDEs, and has helped us develop an improved PC scheme for solving the incompressible Navier-Stokes equations.



There are no comments yet.


page 29


First-order positivity-preserving entropy stable spectral collocation scheme for the 3-D compressible Navier-Stokes equations

In this paper, we extend the positivity-preserving, entropy stable first...

Non-overlapping block smoothers for the Stokes equations

Overlapping block smoothers efficiently damp the error contributions fro...

Defect-Deferred Correction Method Based on a Subgrid Artificial Viscosity Modeling

An alternative first step approximation based on subgrid artificial visc...

Multirate Timestepping for the Incompressible Navier-Stokes Equations in Overlapping Grids

We develop a multirate timestepper for semi-implicit solutions of the un...

A Dissipation Theory for Three-Dimensional FDTD with Application to Stability Analysis and Subgridding

The finite-difference time-domain (FDTD) algorithm is a popular numerica...

Energy stability analysis of turbulent incompressible flow based on the triple decomposition of the velocity gradient tensor

In the context of flow visualization a triple decomposition of the veloc...

Energy conservative SBP discretizations of the acoustic wave equation in covariant form on staggered curvilinear grids

We develop a numerical method for solving the acoustic wave equation in ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.