A Scalable Monolithic Modified Newton Multigrid Framework for Time-Dependent $p$-Navier-Stokes Flow

arXiv:2603.2870638.3h-index: 9
AI Analysis

This work addresses computational bottlenecks in simulating non-Newtonian fluid flows for researchers in computational fluid dynamics, though it is incremental as it builds on existing methods with specific improvements.

The paper tackles the challenge of solving large nonlinear monolithic saddle-point systems from time-dependent $(p,\\delta)$-Navier-Stokes models in the shear-thinning regime by developing a scalable monolithic modified Newton framework that replaces the ill-conditioned exact constitutive tangent with a better-conditioned surrogate, achieving robustness in numerical tests with scalable parallel performance.

Fully implicit tensor-product space-time discretizations of time-dependent $(p,δ)$-Navier-Stokes models yield, on each time step, large nonlinear monolithic saddle-point systems. In the shear-thinning regime $1<p<2$, especially as $p\downarrow 1$ and $δ\downarrow 0$, the decisive difficulty is the constitutive tangent: its ill-conditioning impairs Newton globalization and the preconditioning of the arising linear systems. We therefore develop a scalable monolithic modified Newton framework for tensor-product space-time finite elements in which the exact constitutive tangent in the Jacobian action is replaced by a better-conditioned surrogate. Picard and exact Newton serve as reference linearizations within the same algebraic framework. Scalability is achieved through matrix-free operator evaluation, a monolithic multigrid V-cycle preconditioner, order-preserving reduced Gauss-Radau time quadrature, and an inexact space-time Vanka smoother with single-time-point coefficient freezing in local patch matrices. We prove coercivity of the linearized viscous-Nitsche term in the uniformly elliptic regime $ν_\infty>0$ and consistency of the reduced time quadrature. Numerical tests demonstrate robustness with respect to model parameters, nonlinear and linear iteration counts, and scalable parallel performance.

Foundations

The foundational work for this paper's niche, ranked by how specifically the neighbourhood builds on it — not by global fame.

Your Notes