In the present paper, we discuss implementation details of a free and opensource numerical solver based on the finite volume method for numerical simulation of viscoplastic non-Newtonian fluids. In addition to the fact that they are involved in many industrial applications, both their physical properties and their rheological behavior make them challenging for numerical simulation. Viscoplastic fluids are known to behave as solid unless the shear stress reaches a critical level, known as yield-stress, beyond which they behave as liquid. In most cases, both yielded and unyielded regions coexist in the fluid domain. In mathematical model of viscoplastic fluid, the constitutive equation is a non-differentiable function. This is often overcome by using the approximate constitutive equation that has a regularized form, e.g. the Papanastasiou regularization model. Using the same approach, we assess the influence of regularization parameters on simulation convergence and results accuracy. In this study, we give implementation details of viscoplastic fluid models in freeCappuccino open-source Computational Fluid Dynamics code. Moreover, we perform validation on several well known benchmark cases and compare proposed approach with those existing in published literature. We also perform a parametric analysis and show the effect of Reynolds and Bingham numbers on the extent of the yielded regions. Conclusions of the study have relevance in practical application of computational fluid dynamics to viscoplastic fluids in particular and to non-Newtonian fluids in general.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.