Implicit-explicit Crank-Nicolson scheme for Oseen's equation at high Reynolds number

Erik Burman,Deepika Garg,Johnny Guzman
2024-05-21
Abstract:In this paper we continue the work on implicit-explicit (IMEX) time discretizations for the incompressible Oseen equations that we started in \cite{BGG23} (E. Burman, D. Garg, J. Guzmàn, {\emph{Implicit-explicit time discretization for Oseen's equation at high Reynolds number with application to fractional step methods}}, SIAM J. Numer. Anal., 61, 2859--2886, 2023). The pressure velocity coupling and the viscous terms are treated implicitly, while the convection term is treated explicitly using extrapolation. Herein we focus on the implicit-explicit Crank-Nicolson method for time discretization. For the discretization in space we consider finite element methods with stabilization on the gradient jumps. The stabilizing terms ensures inf-sup stability for equal order interpolation and robustness at high Reynolds number. Under suitable Courant conditions we prove stability of the implicit-explicit Crank-Nicolson scheme in this regime. The stabilization allows us to prove error estimates of order $O(h^{k+\frac12} + \tau^2)$. Here $h$ is the mesh parameter, $k$ the polynomial order and $\tau$ the time step. Finally we discuss some fractional step methods that are implied by the IMEX scheme. Numerical examples are reported comparing the different methods when applied to the Navier-Stokes' equations.
Numerical Analysis
What problem does this paper attempt to address?
This paper attempts to solve the problems of stability and error estimation of the implicit - explicit (IMEX) time discretization method for the Oseen equation at high Reynolds numbers. Specifically, the paper focuses on using the Crank - Nicolson method for time discretization under high Reynolds number conditions, combined with finite - element space discretization and explicit treatment of the convection term (by extrapolation). The main objectives of the paper include: 1. **Stability analysis**: Prove the stability of the implicit - explicit Crank - Nicolson scheme at high Reynolds numbers under appropriate Courant conditions (such as the standard hyperbolic CFL condition or the slightly stronger 4/3 - CFL condition). 2. **Error estimation**: Provide error estimates in the form of \(O(h^{k + 1/2}+\tau^{2})\) at high Reynolds numbers, where \(h\) is the mesh parameter, \(k\) is the polynomial order, and \(\tau\) is the time step. 3. **Numerical experiments**: Compare the performance of different methods when applied to the Navier - Stokes equations through numerical experiments, especially discuss some step - by - step methods implied by the IMEX scheme. ### Key technical points - **Implicit - explicit (IMEX) time discretization**: The pressure - velocity coupling terms and viscous terms are treated implicitly, and the convection term is treated explicitly and implemented by extrapolation. - **Space discretization**: Use the continuous finite - element method with gradient - jump stabilization to ensure stability and robustness at high Reynolds numbers. - **Stabilization techniques**: Introduce the Continuous Interior Penalty (CIP) method to impose penalty terms on the gradient jumps of the cell faces to enhance stability. - **Error analysis**: Prove the stability and error estimates at high Reynolds numbers through a series of technical lemmas and inequalities. ### Main contributions - **Theoretical results**: Provide the stability analysis and error estimates of the IMEX Crank - Nicolson scheme at high Reynolds numbers for the first time in the literature. - **Numerical verification**: Verify the validity of the theoretical results through numerical experiments and show the performance of different methods in practical applications. ### Conclusion Through rigorous mathematical analysis and numerical experiments, the paper proves the stability and accuracy of the IMEX Crank - Nicolson scheme at high Reynolds numbers, providing an effective numerical solution method for high - Reynolds - number fluid - dynamics problems.