This paper presents a stabilized formulation for the generalized Navier-Stokes equations for weak enforcement of essential boundary conditions. The non-Newtonian behavior of blood is modeled via shear-rate dependent constitutive equations. The boundary terms for weak enforcement of Dirichlet boundary conditions are derived via locally resolving the fine-scale variational equation facilitated by the Variational Multiscale (VMS) framework. The proposed method reproduces the consistency and stabilization terms that are present in the Nitsche type approaches. In addition, for the shear-rate fluids, two more boundary terms appear. One of these terms is the viscosity-derivative term and is a function of the shear-rate, while the other term is a zeroth-order term. These terms play an important role in attaining optimal convergence rates for the velocity and pressure fields in the norms considered. A most significant contribution is the form of the stabilization tensors that are also variationally derived. Employing edge functions the edge stabilization tensor is numerically evaluated, and it adaptively adjusts itself to the magnitude of the boundary residual. The resulting formulation is variationally consistent and the weakly imposed no-slip boundary condition leads to higher accuracy of the spatial gradients for coarse boundary-layer meshes when compared with the traditional strongly imposed boundary conditions. This feature of the present approach will be of significance in imposing interfacial continuity conditions across non-matching discretizations in blood-artery interaction problems. A set of test cases is presented to investigate the mathematical attributes of the method and a patient-specific case is presented to show its clinical relevance.
Citation: Soonpil Kang, Sharbel Nashar, Elizabeth R. Livingston, Arif Masud. Weakly imposed boundary conditions for shear-rate dependent non-Newtonian fluids: application to cardiovascular flows[J]. Mathematical Biosciences and Engineering, 2021, 18(4): 3855-3886. doi: 10.3934/mbe.2021193
This paper presents a stabilized formulation for the generalized Navier-Stokes equations for weak enforcement of essential boundary conditions. The non-Newtonian behavior of blood is modeled via shear-rate dependent constitutive equations. The boundary terms for weak enforcement of Dirichlet boundary conditions are derived via locally resolving the fine-scale variational equation facilitated by the Variational Multiscale (VMS) framework. The proposed method reproduces the consistency and stabilization terms that are present in the Nitsche type approaches. In addition, for the shear-rate fluids, two more boundary terms appear. One of these terms is the viscosity-derivative term and is a function of the shear-rate, while the other term is a zeroth-order term. These terms play an important role in attaining optimal convergence rates for the velocity and pressure fields in the norms considered. A most significant contribution is the form of the stabilization tensors that are also variationally derived. Employing edge functions the edge stabilization tensor is numerically evaluated, and it adaptively adjusts itself to the magnitude of the boundary residual. The resulting formulation is variationally consistent and the weakly imposed no-slip boundary condition leads to higher accuracy of the spatial gradients for coarse boundary-layer meshes when compared with the traditional strongly imposed boundary conditions. This feature of the present approach will be of significance in imposing interfacial continuity conditions across non-matching discretizations in blood-artery interaction problems. A set of test cases is presented to investigate the mathematical attributes of the method and a patient-specific case is presented to show its clinical relevance.
| [1] |
L. Botti, M. Piccinelli, B. Ene-Iordache, A. Remuzzi, L. Antiga, An adaptive mesh refinement solver for large-scale simulation of biological flows, Int. J. Numer. Method. Biomed. Eng., 26 (2010), 86-100. doi: 10.1002/cnm.1257
|
| [2] |
G. De Santis, P. Mortier, M. De Beule, P. Segers, P. Verdonck, B. Verhegghe, Patient-specific computational fluid dynamics: structured mesh generation from coronary angiography, Med. Biol. Eng. Comput., 48 (2010), 371-380. doi: 10.1007/s11517-010-0583-4
|
| [3] | J. Müller, O. Sahni, X. Li, K. E. Jansen, M. S. Shephard, C. A. Taylor, Anisotropic adaptive finite element method for modelling blood flow, Comput. Methods Biomech. Biomed. Eng., 8 (2005), 295-305. |
| [4] |
Y. Bazilevs, T. J. R. Hughes, Weak imposition of Dirichlet boundary conditions in fluid mechanics, Comput. Fluids, 36 (2007), 12-26. doi: 10.1016/j.compfluid.2005.07.012
|
| [5] |
Y. Bazilevs, C. Michler, V. M. Calo, T. J. R. Hughes, Weak Dirichlet boundary conditions for wall-bounded turbulent flows, Comput. Methods Appl. Mech. Eng., 196 (2007), 4853-4862. doi: 10.1016/j.cma.2007.06.026
|
| [6] | R. J. Cier, S. Rojas, V. M. Calo, A nonlinear weak constraint enforcement method for advection-dominated diffusion problems, Mech. Res. Commun., In press. |
| [7] |
P. Hansbo, M. Juntunen, Weakly imposed Dirichlet boundary conditions for the Brinkman model of porous media flow, Appl. Numer. Math., 59 (2009), 1274-1289. doi: 10.1016/j.apnum.2008.07.003
|
| [8] |
S. Xu, B. Gao, M. C. Hsu, B. Ganapathysubramanian, A residual-based variational multiscale method with weak imposition of boundary conditions for buoyancy-driven flows, Comput. Methods Appl. Mech. Eng., 352 (2019), 345-368. doi: 10.1016/j.cma.2019.03.057
|
| [9] |
D. Kamensky, M. C. Hsu, D. Schillinger, J. A. Evans, A. Aggarwal, Y. Bazilevs et al., An immersogeometric variational framework for fluid-structure interaction: Application to bioprosthetic heart valves, Comput. Methods Appl. Mech. Eng., 284 (2015), 1005-1053. doi: 10.1016/j.cma.2014.10.040
|
| [10] | S. Kang, A. Masud, A Variational Multiscale method with immersed boundary conditions for incompressible flows, Meccanica, 56 (2021) 1397–1422 |
| [11] | F. Xu, E. L. Johnson, C. Wang, A. Jafari, C. H. Yang, M. S. Sacks et al., Computational investigation of left ventricular hemodynamics following bioprosthetic aortic and mitral valve replacement, Mech. Res. Commun., (2021), In press. |
| [12] |
E. Burman, A penalty-free nonsymmetric Nitsche-type method for the weak imposition of boundary conditions, SIAM J. Numer. Anal., 50 (2012), 1959-1981. doi: 10.1137/10081784X
|
| [13] |
Y. I. Cho, K. R. Kensey, Effects of the non-Newtonian viscosity of blood on flows in a diseased arterial vessel. Part 1: Steady flows, Biorheology, 28 (1991), 241-262. doi: 10.3233/BIR-1991-283-415
|
| [14] |
Elhanafy, A. Guaily, A. Elsaid, Numerical simulation of blood flow in abdominal aortic aneurysms: Effects of blood shear-thinning and viscoelastic properties, Math. Comput. Simul., 160 (2019), 55-71. doi: 10.1016/j.matcom.2018.12.002
|
| [15] |
J. C. Weddell, J. Kwack, P. I. Imoukhuede, A. Masud, Hemodynamic analysis in an idealized artery tree: differences in wall shear stress between Newtonian and non-Newtonian blood models, PLoS One. DOI:10.1371/journal. pone. 0124575. doi: 10.1371/journal.pone.0124575
|
| [16] |
Aguirre, E. Castillo, M. Cruchaga, R. Codina, J. Baiges, Stationary and time-dependent numerical approximation of the lid-driven cavity problem for power-law fluid flows at high Reynolds numbers using a stabilized finite element formulation of the VMS type, J. Non-Newton. Fluid Mech., 257 (2018), 22-43. doi: 10.1016/j.jnnfm.2018.03.014
|
| [17] | J. Kwack, S. Kang, G. Bhat, A. Masud, Time-dependent outflow boundary conditions for blood flow in the arterial system, In: Bazilevs Y., Takizawa K. (eds) Advances in Computational Fluid-Structure Interaction and Flow Simulation, Modeling and Simulation in Science, Engineering and Technology, Springer, DOI10.1007/978-3-319-40827-9_28, (2016), 359-378. |
| [18] |
J. Kwack, A. Masud, A stabilized mixed finite element method for shear-rate dependent non-Newtonian fluids: 3D benchmark problems and application to blood flow in bifurcating arteries, Comput. Mech., 53 (2014), 751-776. doi: 10.1007/s00466-013-0928-6
|
| [19] | J. Kwack, A. Masud, K. R. Rajagopal, Stabilized mixed three-field formulation for a generalized incompressible Oldroyd-B model, Int. J. Numer. Methods Fluids, 83 (2016), 704-734. |
| [20] |
Masud, J. Kwack, A stabilized mixed finite element method for the incompressible shear-rate dependent non-Newtonian fluids: Variational Multiscale framework and consistent linearization, Comput. Methods Appl. Mech. Eng., 200 (2011), 577-596. doi: 10.1016/j.cma.2010.08.012
|
| [21] | Masud, M. Masud, Blood flow analysis in arteries with atherosclerosis, in Proceedings of the 1995 Bioengineering Conference, Beever Creek, Colorado, USA, Jun. 28-Jul. 2, (1995), 525-526. |
| [22] |
T. J. Truster, A. Masud, Primal interface formulation for coupling multiple PDEs: A consistent derivation via the Variational Multiscale method, Comput. Methods Appl. Mech. Eng., 268 (2014), 194-224. doi: 10.1016/j.cma.2013.08.005
|
| [23] |
P. Chen, T. J. Truster, A. Masud, Interfacial stabilization at finite strains for weak and strong discontinuities in multi-constituent materials, Comput. Methods Appl. Mech. Eng., 328 (2018), 717-751. doi: 10.1016/j.cma.2017.06.020
|
| [24] |
J. Boyd, J. M. Buick, S. Green, Analysis of the Casson and Carreau-Yasuda non-Newtonian blood models in steady and oscillatory flows using the lattice Boltzmann method, Phys. Fluids., 19 (2007), 093103. doi: 10.1063/1.2772250
|
| [25] |
T. J. Truster, P. Chen, A. Masud, Finite strain primal interface formulation with consistently evolving stabilization, Int. J. Numer. Methods Eng., 102 (2015), 278-315. doi: 10.1002/nme.4763
|
| [26] |
U. Ghia, K. N. Ghia, C. T. Shin, High-Re solutions for incompressible flow using the Navier-stokes equations and a multigrid method, J. Comput. Phys., 48 (1982), 387-411. doi: 10.1016/0021-9991(82)90058-4
|
| [27] |
Masud, R. Calderer, A variational multiscale stabilized formulation for the incompressible Navier-Stokes, Comput. Mech., 44 (2009), 145-160. doi: 10.1007/s00466-008-0362-3
|
| [28] |
R. Mahmood, S. Bilal, I. Khan, N. Kousar, A. H. Seikh, E. M. Sherif, A comprehensive finite element examination of Carreau Yasuda fluid model in a lid driven cavity and channel with obstacle by way of kinetic energy and drag and lift coefficient measurements, J. Mater. Res. Technol., 9 (2020), 1785-1800. doi: 10.1016/j.jmrt.2019.12.010
|
| [29] |
P. H. Bovendeerd, A. A. Van Steenhoven, F. N. V. de Vosse, G. Vossers, Steady entry flow in a curved pipe, J. Fluid Mech., 177 (1987), 233-246. doi: 10.1017/S0022112087000934
|
| [30] |
E. Vignon-Clementel, C. A. Figueroa, K. E. Jansen, C. A. Taylor, Outflow boundary conditions for three-dimensional finite element modeling of blood flow and pressure inarteries, Comput. Methods Appl. Mech. Eng., 195 (2006) 3776-3796. doi: 10.1016/j.cma.2005.04.014
|
| [31] |
M. E. Moghadam, Y. Bazilevs, T. Y. Hsia, I. E. Vignon-Clementel, A. L. Marsden, A comparison of outlet boundary treatments for prevention of backflow divergence with relevance to blood flow simulations, Comput. Mech., 48 (2011), 277-291. doi: 10.1007/s00466-011-0599-0
|
| [32] |
J. Liu, W. Yang, I. S. Lan, A. L. Marsden, Fluid-structure interaction modeling of blood flow in the pulmonary arteries using the unified continuum and variational multiscale formulation, Mech. Res. Commun., 107 (2020), 103556. doi: 10.1016/j.mechrescom.2020.103556
|
Figures(27) / Tables(5)
Soonpil Kang, Sharbel Nashar, Elizabeth R. Livingston, Arif Masud. Weakly imposed boundary conditions for shear-rate dependent non-Newtonian fluids: application to cardiovascular flows[J]. Mathematical Biosciences and Engineering, 2021, 18(4): 3855-3886. doi: 10.3934/mbe.2021193
DownLoad: