Unconditional stability in 3d of a sparse grad-div approximation of the Navier-Stokes equations

12/13/2021
by   William Layton, et al.
0

Inclusion of a term -γ∇∇· u, forcing ∇· u to be pointwise small, is an effective tool for improving mass conservation in discretizations of incompressible flows. However, the added grad-div term couples all velocity components, decreases sparsity and increases the condition number in the linear systems that must be solved every time step. To address these three issues various sparse grad-div regularizations and a modular grad-div method have been developed. We develop and analyze herein a synthesis of a fully decoupled, parallel sparse grad-div method of Guermond and Minev with the modular grad-div method. Let G^∗=-diag(∂ _x^2,∂ _y^2,∂ _z^2) denote the diagonal of G=-∇∇·, and α≥ 0 an adjustable parameter. The 2-step method considered is 1 : u^n+1-u^n/k+u^n·∇u^n+1+∇ p^n+1-νΔu^n+1=f ∇·u^n+1=0, 2 : [ 1/kI+(γ +α )G^∗] u^n+1=1/k u^n+1+[ (γ +α )G^∗-γ G] u^n. The analysis also establishes that the method controls the persistent size of ‖∇· u ‖ in general and controls the transients in ‖∇· u‖ for a cold start when α >0.5γ. Consistent numerical tests are presented.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
01/11/2022

A Doubly Adaptive Penalty Method for the Navier Stokes Equations

We develop, analyze and test adaptive penalty parameter methods. We prov...
research
09/05/2023

Variable Time Step Method of DAHLQUIST, LINIGER and NEVANLINNA (DLN) for a Corrected Smagorinsky Model

Turbulent flows strain resources, both memory and CPU speed. The DLN met...
research
10/31/2019

An analysis of a linearly extrapolated BDF2 subgrid artificial viscosity method for incompressible flows

This report extends the mathematical support of a subgrid artificial vis...
research
01/27/2020

Modular grad-div stabilization for multiphysics flow problems

This paper considers a modular grad-div stabilization method for approxi...
research
02/10/2023

Approximation and Structured Prediction with Sparse Wasserstein Barycenters

We develop a general theoretical and algorithmic framework for sparse ap...

Please sign up or login with your details

Forgot password? Click here to reset