A Discontinuous Galerkin and Semismooth Newton Approach for the Numerical Solution of Bingham Flow with Variable Density

This paper is devoted to the study of Bingham flow with variable density. We propose a local bi-viscosity regularization of the stress tensor based on a Huber smoothing step. Next, our computational approach is based on a second-order, divergence-conforming discretization of the Huber regularized Bingham constitutive equations, coupled with a discontinuous Galerkin scheme for the mass density. We take advantage of the properties of divergence-conforming and discontinuous Galerkin formulations to effectively incorporate upwind discretizations, thereby ensuring the stability of the formulation. The stability of the continuous problem and the fully discrete scheme are analyzed. Further, a semismooth Newton method is proposed for solving the obtained fully discretized system of equations at each time step. Finally, several numerical examples that illustrate the main features of the problem and the properties of the numerical scheme are presented.

Automated summary

This paper focuses on the study of Bingham flow with variable density and proposes a local bi-viscosity regularization of the stress tensor. The computational approach is based on a second-order, divergence-conforming discretization and a discontinuous Galerkin scheme. The stability of the formulation is ensured by incorporating upwind discretizations. A semismooth Newton method is used to solve the fully discretized system of equations at each time step. Numerical examples are presented to illustrate the main features of the problem and the properties of the numerical scheme.