A polytopal discontinuous Galerkin method for the pseudo-stress formulation of the unsteady Stokes problem

Paola F. Antonietti,Michele Botti,Alessandra Cancrini,Ilario Mazzieri

2024-08-16

Abstract:This work aims to construct and analyze a discontinuous Galerkin method on polytopal grids (PolydG) to solve the pseudo-stress formulation of the unsteady Stokes problem. The pseudo-stress variable is introduced due to the growing interest in non-Newtonian flows and coupled interface problems, where stress assumes a fundamental role. The space-time discretization of the problem is achieved by combining the PolydG approach with the implicit theta-method time integration scheme. For both the semi- and fully-discrete problems we present a detailed stability analysis. Moreover, we derive convergence estimates for the fully discrete space-time discretization. A set of verification tests is presented to verify the theoretical estimates and the application of the method to cases of engineering interest.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop and analyze a method for solving the pseudo - stress form of the unsteady Stokes problem using the discontinuous Galerkin method (PolydG) on a polyhedral mesh. Specifically, the author introduced the pseudo - stress variable \(\sigma(u, p)=\mu \nabla u - pI_d\), where \(u\) is the fluid velocity, \(p\) is the pressure, \(\mu\) is the fluid viscosity, and \(I_d\) is the \(d\)-dimensional identity matrix.

### Main research objectives
1. **Construct and analyze the discontinuous Galerkin method**: Propose a discontinuous Galerkin method based on a polyhedral mesh for the unsteady Stokes problem in the pseudo - stress form.
2. **Space - time discretization**: Combine PolydG space discretization with the implicit \(\theta\)-method time integration scheme to achieve space - time discretization.
3. **Stability analysis**: Conduct a detailed stability analysis of the semi - discrete and fully - discrete problems.
4. **Convergence estimate**: Derive the convergence estimate of the fully - discrete space - time discretization.
5. **Numerical verification**: Verify the theoretical estimates through a series of verification tests and demonstrate the effectiveness of this method in engineering applications, such as the flow around a cylinder problem.

### Significance of the pseudo - stress variable
The introduction of the pseudo - stress variable is mainly to meet the requirements of non - Newtonian fluid flow and coupled interface problems, in which stress plays a crucial role. In addition, the pseudo - stress variable can simplify the design of approximation methods for complex nonlinear flow problems and their numerical solutions.

### Space - time discretization
- **Space discretization**: The PolydG method is adopted.
- **Time discretization**: The implicit \(\theta\)-method is adopted. In particular, the Crank - Nicolson method is used for full discretization.

### Stability and convergence analysis
- **Stability**: Prove the stability of the semi - discrete and fully - discrete problems.
- **Convergence**: Establish the error estimate of the fully - discrete problem under an appropriate discrete norm.

### Numerical experiments
- **Verification tests**: Verify the correctness of the theoretical estimates.
- **Engineering applications**: Such as the application in the flow around a cylinder problem.

Through these works, the author provides a new and effective numerical method for solving complex unsteady Stokes problems and conducts a strict theoretical analysis of its stability and convergence.