A novel formulation for Neumann inflow boundary conditions in biomechanics

Neumann boundary conditions prescribing the total momentum flux at inflow boundaries of biomechanical problems are proposed in this study. This approach enables the simultaneous application of velocity/flow rate and pressure curves at inflow boundaries. As the basic numerical method, a residual-based variational multiscale (or stabilized) finite element method is presented. The focus of the numerical examples in this work is on respiratory flows with complete flow reversals. However, the proposed formulation is just as well suited for cardiovascular flow problems with partial retrograde flow. Instabilities, which were reported for such problems in the literature, are resolved by the present approach without requiring the additional consideration of a Lagrange multiplier technique. The suitability of the approach is demonstrated for two respiratory flow examples, a rather simple tube and complex tracheobronchial airways (up to the fourth generation, segmented from end-expiratory CT images). For the latter example, the boundary conditions are generated from mechanical ventilation data obtained from an intensive care unit patient suffering from acute lung injury. For the tube, analytical pressure profiles can be replicated, and for the tracheobronchial airways, a correct distribution of the prescribed total momentum flux at the inflow boundary into velocity and pressure part is observed.

[1]  Christian Vergara,et al.  Nitsche’s Method for Defective Boundary Value Problems in Incompressibile Fluid-dynamics , 2011, J. Sci. Comput..

[2]  A. Green,et al.  Modelling of peak-flow wall shear stress in major airways of the lung. , 2004, Journal of biomechanics.

[3]  G. Hulbert,et al.  A generalized-α method for integrating the filtered Navier–Stokes equations with a stabilized finite element method , 2000 .

[4]  Wolfgang A. Wall,et al.  Towards a comprehensive computational model for the respiratory system , 2010 .

[5]  Charles A. Taylor,et al.  Effect of Curvature on Displacement Forces Acting on Aortic Endografts: A 3-Dimensional Computational Analysis , 2009, Journal of endovascular therapy : an official journal of the International Society of Endovascular Specialists.

[6]  T. Hughes,et al.  Isogeometric variational multiscale modeling of wall-bounded turbulent flows with weakly enforced boundary conditions on unstretched meshes , 2010 .

[7]  Wolfgang A. Wall,et al.  An algebraic variational multiscale-multigrid method for large-eddy simulation of turbulent variable-density flow at low Mach number , 2010, J. Comput. Phys..

[8]  Wolfgang A. Wall,et al.  An algebraic variational multiscale-multigrid method based on plain aggregation for convection-diffusion problems , 2009 .

[9]  H. M. Bücker,et al.  Sensitivity of optimal shapes of artificial grafts with respect to flow parameters , 2010 .

[10]  M Ismail,et al.  Coupled and reduced dimensional modeling of respiratory mechanics during spontaneous breathing , 2013, International journal for numerical methods in biomedical engineering.

[11]  Charles A. Taylor,et al.  Augmented Lagrangian method for constraining the shape of velocity profiles at outlet boundaries for three-dimensional finite element simulations of blood flow , 2009 .

[12]  John Kim,et al.  DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOWS UP TO RE=590 , 1999 .

[13]  P. Blanco,et al.  A unified variational approach for coupling 3D-1D models and its blood flow applications , 2007 .

[14]  M. Kronbichler,et al.  An algebraic variational multiscale-multigrid method for large eddy simulation of turbulent flow , 2010 .

[15]  Frédéric Valentin,et al.  An unusual stabilized finite element method for a generalized Stokes problem , 2002, Numerische Mathematik.

[16]  Clement Kleinstreuer,et al.  Transient airflow structures and particle transport in a sequentially branching lung airway model , 2002 .

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

[18]  Alfio Quarteroni,et al.  Numerical Treatment of Defective Boundary Conditions for the Navier-Stokes Equations , 2002, SIAM J. Numer. Anal..

[19]  L. Franca,et al.  On an Improved Unusual Stabilized Finite Element Method for theAdvective-Reactive-Diffusive Equation , 1999 .

[20]  Wolfgang A. Wall,et al.  An algebraic variational multiscale–multigrid method for large‐eddy simulation of turbulent pulsatile flows in complex geometries with detailed insight into pulmonary airway flow , 2013 .

[21]  Wolfgang Schröder,et al.  Numerical investigation of the three-dimensional flow in a human lung model. , 2008, Journal of biomechanics.

[22]  Wolfgang A Wall,et al.  Structured Tree Impedance Outflow Boundary Conditions for 3D Lung Simulations , 2010, Journal of biomechanical engineering.

[23]  Kenneth E. Jansen,et al.  A better consistency for low-order stabilized finite element methods , 1999 .

[24]  Erik Burman,et al.  Stabilized finite element methods for the generalized Oseen problem , 2007 .

[25]  Eric A. Hoffman,et al.  Airway Wall Stiffening Increases Peak Wall Shear Stress: A Fluid–Structure Interaction Study in Rigid and Compliant Airways , 2010, Annals of Biomedical Engineering.

[26]  Wolfgang A. Wall,et al.  Residual‐based variational multiscale methods for laminar, transitional and turbulent variable‐density flow at low Mach number , 2011 .

[27]  Wolfgang A. Wall,et al.  An algebraic variational multiscale-multigrid method for large-eddy simulation: generalized-α time integration, Fourier analysis and application to turbulent flow past a square-section cylinder , 2011 .

[28]  D. Saloner,et al.  Numerical analysis of flow through a severely stenotic carotid artery bifurcation. , 2002, Journal of biomechanical engineering.

[29]  W. Wall,et al.  Nanoparticle transport in a realistic model of the tracheobronchial region , 2010 .

[30]  Arne Erik Holdø,et al.  Numerical simulation of respiratory flow patterns within human lung , 2002, Respiratory Physiology & Neurobiology.

[31]  Thomas J. R. Hughes,et al.  Conservation properties for the Galerkin and stabilised forms of the advection–diffusion and incompressible Navier–Stokes equations , 2005 .

[32]  Charles A. Taylor,et al.  Outflow boundary conditions for three-dimensional finite element modeling of blood flow and pressure in arteries , 2006 .

[33]  C. Taylor,et al.  Predictive medicine: computational techniques in therapeutic decision-making. , 1999, Computer aided surgery : official journal of the International Society for Computer Aided Surgery.

[34]  L. Formaggia,et al.  On the stability of the coupling of 3D and 1D fluid-structure interaction models for blood flow simulations , 2007 .

[35]  T. Hughes,et al.  The Galerkin/least-squares method for advective-diffusive equations , 1988 .

[36]  Kenneth R. Lutchen,et al.  An Anatomically Based Hybrid Computational Model of the Human Lung and its Application to Low Frequency Oscillatory Mechanics , 2006, Annals of Biomedical Engineering.

[37]  Ananth V. Annapragada,et al.  Computational Fluid Dynamics Simulation of Airflow and Aerosol Deposition in Human Lungs , 2004, Annals of Biomedical Engineering.

[38]  L. Franca,et al.  Stabilized finite element methods. II: The incompressible Navier-Stokes equations , 1992 .

[39]  William L. Roberts,et al.  A new methodology for targeting drug-aerosols in the human respiratory system , 2008 .

[40]  Eric A Hoffman,et al.  Characteristics of airflow in a CT-based ovine lung: a numerical study. , 2007, Journal of applied physiology.

[41]  Charles A. Taylor,et al.  On Coupling a Lumped Parameter Heart Model and a Three-Dimensional Finite Element Aorta Model , 2009, Annals of Biomedical Engineering.