Numerical verification of PolyDG algebraic solvers for the pseudo-stress Stokes problem

72/2025

Numerical verification of PolyDG algebraic solvers for the pseudo-stress Stokes problem

Monday 1st December 2025

Antonietti, P. F.; Cancrini, A.; Ciaramella, G.

This work focuses on the development of efficient solvers for the pseudo-stress formulation of the unsteady Stokes problem, discretised by means of a discontinuous Galerkin method on polytopal grids (PolyDG). The introduction of the pseudo-stress variable is motivated by the growing interest in non-Newtonian flow models and coupled interface problems, where the stress field plays a fundamental role in the physical description. The space-time discretisation of the problem is obtained by combining the PolyDG approach in space with the implicit Euler method for time integration. The resulting linear system, characterised by a symmetric, positive, definite matrix, exhibits deteriorating convergence with standard solvers as the time step decreases. To address this issue, we investigate two tailored strategies: deflated Conjugate Gradient, which mitigates the effect of the most problematic eigenmodes, and collective Block-Jacobi, which exploits the block structure of the system matrix. Numerical experiments show that both approaches yield iteration counts effectively independent of \Delta t, ensuring robust performance with respect to the time step. Future work will focus on extending this robustness to the spatial discretisation parameter h by integrating multigrid strategies with the time-robust solvers developed in this study.

This report, or a modified version of it, has been also submitted to, or published on
Proceedings of Domain Decomposition Methods in Science and Engineering XXIX