A continuous finite element framework for the pressure Poisson equation allowing non‐Newtonian and compressible flow behavior

Computing pressure fields from given flow velocities is a task arising in various engineering, biomedical and scientific computing applications. The so-called pressure Poisson equation (PPE) derived from the balance of linear momentum provides an attractive framework for such a task. However, the PPE increases the regularity requirements on the pressure and velocity spaces, thereby imposing theoretical and practical challenges for its application. In order to stay within a Lagrangian finite element framework, it is common practice to completely neglect the influence of viscosity and compressibility when computing the pressure, which limits the practical applicability of the pressure Poisson method. In this context, we present a mixed finite element framework which enables the use of this popular technique with generalised Newtonian fluids (e.g., blood) and compressible flows, while allowing standard finite element spaces to be employed for the given data and unknowns. This is accomplished through the use of appropriate vector calculus identities and simple projections of certain flow quantities. In the compressible case, the mixed formulation also includes an additional equation for retrieving the density field from the given velocities so that the pressure can be accurately determined. The potential of this new approach is showcased through numerical examples.

[1]  Hans Johnston,et al.  Accurate, stable and efficient Navier-Stokes solvers based on explicit treatment of the pressure term , 2004 .

[2]  Y. Cho,et al.  Effects of the non-Newtonian viscosity of blood on flows in a diseased arterial vessel. Part 1: Steady flows. , 1991, Biorheology.

[3]  David Nordsletten,et al.  Relative pressure estimation from velocity measurements in blood flows: State‐of‐the‐art and new approaches , 2018, International journal for numerical methods in biomedical engineering.

[4]  H. Binous Introducing Non-Newtonian Fluid Mechanics Computations with Mathematica in the Undergraduate Curriculum. , 2007 .

[5]  S. Frey,et al.  Finite element approximations for quasi-Newtonian flows employing a multi-field GLS method , 2011 .

[6]  T. Sochi Analytical solutions for the flow of Carreau and Cross fluids in circular pipes and thin slits , 2015, Rheologica Acta.

[7]  Cameron V. King,et al.  Assessment of pressure field calculations from particle image velocimetry measurements , 2010 .

[8]  Jie Liu,et al.  Open and traction boundary conditions for the incompressible Navier-Stokes equations , 2009, J. Comput. Phys..

[9]  Rodolfo Ruben Rosales,et al.  An efficient method for the incompressible Navier-Stokes equations on irregular domains with no-slip boundary conditions, high order up to the boundary , 2010, J. Comput. Phys..

[10]  B. W. Oudheusden,et al.  PIV-based pressure measurement , 2013 .

[11]  Olaf Steinbach,et al.  On the influence of the wall shear stress vector form on hemodynamic indicators , 2017, Comput. Vis. Sci..

[12]  Pavel B. Bochev,et al.  An Absolutely Stable Pressure-Poisson Stabilized Finite Element Method for the Stokes Equations , 2004, SIAM J. Numer. Anal..

[13]  Jaroslav Hron,et al.  Determination of pressure data from velocity data with a view toward its application in cardiovascular mechanics. Part 1. Theoretical considerations , 2016 .

[14]  Rolf Rannacher,et al.  ARTIFICIAL BOUNDARIES AND FLUX AND PRESSURE CONDITIONS FOR THE INCOMPRESSIBLE NAVIER–STOKES EQUATIONS , 1996 .

[15]  Jian-Guo Liu,et al.  Stable and accurate pressure approximation for unsteady incompressible viscous flow , 2010, J. Comput. Phys..

[16]  Vincent J. Ervin,et al.  Numerical Approximation of a Quasi-Newtonian Stokes Flow Problem with Defective Boundary Conditions , 2007, SIAM J. Numer. Anal..

[17]  L G Raguin,et al.  Methodology to estimate the relative pressure field from noisy experimental velocity data , 2008 .

[18]  Hyo Won Choi,et al.  Numerical study of the impact of non-Newtonian blood behavior on flow over a two-dimensional backward facing step. , 2005, Biorheology.

[19]  Violeta Karyofylli,et al.  Simplex space‐time meshes in compressible flow simulations , 2018, International Journal for Numerical Methods in Fluids.

[20]  Thomas Richter,et al.  Fluid-structure Interactions , 2017 .

[21]  Richard M. Leahy,et al.  Determining cardiac velocity fields and intraventricular pressure distribution from a sequence of ultrafast CT cardiac images , 1994, IEEE Trans. Medical Imaging.

[22]  Zhiqiang Sheng,et al.  A high‐order scheme for the incompressible Navier–Stokes equations with open boundary condition , 2013 .

[23]  Pablo Lamata,et al.  Aortic Relative Pressure Components Derived from Four-Dimensional Flow Cardiovascular Magnetic Resonance , 2013, Magnetic resonance in medicine.

[24]  Thomas A. Manteuffel,et al.  Weighted least-squares finite elements based on particle imaging velocimetry data , 2010, J. Comput. Phys..

[25]  Graham F. Carey,et al.  Performance of iterative methods for Newtonian and generalized Newtonian flows , 1989 .

[26]  A. Huerta,et al.  Finite Element Methods for Flow Problems , 2003 .

[27]  Jean-Frédéric Gerbeau,et al.  Modified Navier-Stokes equations for the outflow boundary conditions in hemodynamics , 2016 .

[28]  Numerical pressure retrieval from velocity measurement of a turbulent tornado-like vortex , 2018 .

[29]  John O Dabiri,et al.  An algorithm to estimate unsteady and quasi-steady pressure fields from velocity field measurements , 2013, Journal of Experimental Biology.

[30]  D N Firmin,et al.  Computation of flow pressure fields from magnetic resonance velocity mapping , 1996, Magnetic resonance in medicine.

[31]  Arnold Reusken,et al.  An extended pressure finite element space for two-phase incompressible flows with surface tension , 2007, J. Comput. Phys..

[32]  Wei Gao,et al.  A hybrid finite volume/finite element method for incompressible generalized Newtonian fluid flows on unstructured triangular meshes , 2009 .

[33]  Pablo Lamata,et al.  A finite-element approach to the direct computation of relative cardiovascular pressure from time-resolved MR velocity data , 2012, Medical Image Anal..

[34]  R. Sani,et al.  On pressure boundary conditions for the incompressible Navier‐Stokes equations , 1987 .