Sergio González-Andrade

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.

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 the divergence conforming and discontinuous Galerkin formulations to incorporate upwind discretizations to stabilize the formulation. The stability of the continuous problem and the full-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.

This paper aims to simulate viscoplastic flow in a shallow-water regime. We specifically use the Bingham model in which the material behaves as a solid if the stress is below a certain threshold, otherwise, it moves as a fluid. The main difficulty of this problem is the coupling of the shallow-water equations with the viscoplastic constitutive laws and the high computational effort needed in its solution. Although there have been many studies of this problem, most of these works use explicit methods with simplified empirical models. In our work, to accommodate non-uniform grids and complicated geometries, we use the discontinuous Galerkin method to solve shallow viscoplastic flows. This method is attractive due to its high parallelization, h- and p-adaptivity, and ability to capture shocks. Additionally, we treat the discontinuities in the interfaces between elements with numerical fluxes that ensure a stable solution of the nonlinear hyperbolic equations. To couple the Bingham model with the shallow-water equations, we regularize the problem with three alternatives. Finally, in order to show the effectiveness of our approach, we perform numerical examples for the usual benchmarks of the shallow-water equations.

In this paper, we propose a dual-mixed formulation for stationary viscoplastic flows with yield, such as the Bingham or the Herschel-Bulkley flow. The approach is based on a Huber regularization of the viscosity term and a two-fold saddle point nonlinear operator equation for the resulting weak formulation. We provide the uniqueness of solutions for the continuous formulation and propose a discrete scheme based on Arnold-Falk-Winther finite elements. The discretization scheme yields a system of slantly differentiable nonlinear equations, for which a semismooth Newton algorithm is proposed and implemented. Local superlinear convergence of the method is also proved. Finally, we perform several numerical experiments in two and three dimensions to investigate the behavior and efficiency of the method.