H(div)-conforming and discontinuous Galerkin approach for Herschel–Bulkley flow with density-dependent viscosity and yield stress

This paper presents a comprehensive study on Herschel–Bulkley flow, where the flow parameters are dependent on the density. The Herschel–Bulkley model is a generalized power-law model used to simulate viscoplastic fluids defined by a plasticity threshold. We consider the case where the plasticity threshold and the viscosity depend on the shear rate and fluid density. To analyze this model, we use a Huber regularization of the stress and propose an H(div)-conforming and discontinuous Galerkin (DG) numerical approximation for the coupled equations governing the flow. We discuss the stability and existence of discrete solutions and propose a semismooth Newton linearization for the numerical solution of the discretized system. Our numerical scheme is validated through several experiments that explore the behavior of Herschel–Bulkley flow under different conditions. The results demonstrate the robustness of our numerical method.