A fourth‐order Cartesian local mesh refinement method for the computational fluid dynamics of physiological flow in multi‐generation branched vessels

Since abnormal fluid states in our body cause critical diseases, patient‐specific computational fluid dynamics (CFD) probably become a standard diagnosis tool in the near future. The vessels in our body are multiple‐branched tubes, which makes it difficult to obtain accurate solutions from conventional CFD methods. In this report, we propose a fourth‐order local mesh refinement (LMR) method based on an interpolated differential operator scheme for simulating flow in multi‐generation branched vessels. The proposed LMR method has the accuracy of fourth‐order for three‐dimensional advection and diffusion equations, respectively. We describe how to apply the LMR method to patient‐specific pulmonary airflow simulations. In our method, the computational mesh size is determined locally by geometrical parameters: the diameter of airways and the distance from the airway wall. To demonstrate our method, an LMR model and a fine mesh model were compared for flow in the central airway, and there was no significant difference between results. We also show the applicability of the method to a maximum eleventh‐generation airway model, where the number of computational nodes was reduced by 85% compared with the case using uniform fine meshes. Copyright © 2010 John Wiley & Sons, Ltd.

[1]  Kenichi Funamoto,et al.  Numerical Validation of MR-Measurement-Integrated Simulation of Blood Flow in a Cerebral Aneurysm , 2009, Annals of Biomedical Engineering.

[2]  C. Caro Discovery of the Role of Wall Shear in Atherosclerosis , 2009, Arteriosclerosis, thrombosis, and vascular biology.

[3]  G E Karniadakis,et al.  LARGE‐SCALE SIMULATION OF THE HUMAN ARTERIAL TREE , 2009, Clinical and experimental pharmacology & physiology.

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

[5]  Phillip Colella,et al.  A cell-centered adaptive projection method for the incompressible Navier-Stokes equations in three dimensions , 2007, J. Comput. Phys..

[6]  E. Hoffman,et al.  Characteristics of the turbulent laryngeal jet and its effect on airflow in the human intra-thoracic airways , 2007, Respiratory Physiology & Neurobiology.

[7]  Yohsuke Imai,et al.  Accuracy study of the IDO scheme by Fourier analysis , 2006, J. Comput. Phys..

[8]  Yohsuke Imai,et al.  A higher-order implicit IDO scheme and its CFD application to local mesh refinement method , 2006 .

[9]  Yohsuke Imai,et al.  Stable coupling between vector and scalar variables for the IDO scheme on collocated grids , 2006, J. Comput. Phys..

[10]  P. Colella,et al.  A node-centered local refinement algorithm for Poisson's equation in complex geometries , 2004 .

[11]  S. Sherwin,et al.  One-dimensional modelling of a vascular network in space-time variables , 2003 .

[12]  David A. Steinman,et al.  Image-Based Computational Fluid Dynamics Modeling in Realistic Arterial Geometries , 2002, Annals of Biomedical Engineering.

[13]  T. Yabe,et al.  The constrained interpolation profile method for multiphase analysis , 2001 .

[14]  Takayuki Aoki,et al.  Interpolated differential operator (IDO) scheme for solving partial differential equations , 1997 .

[15]  S. Lele Compact finite difference schemes with spectral-like resolution , 1992 .

[16]  Takashi Yabe,et al.  A universal solver for hyperbolic equations by cubic-polynomial interpolation I. One-dimensional solver , 1991 .

[17]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[18]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[19]  C. Carrington,et al.  Morphometry of the Human Lung , 1965, The Yale Journal of Biology and Medicine.

[20]  T Gemci,et al.  Computational model of airflow in upper 17 generations of human respiratory tract. , 2008, Journal of biomechanics.

[21]  F. Xiao,et al.  Three-dimensional numerical simulation of flows with complex geometries in a regular Cartesian grid and its application to blood flow in cerebral artery with multiple aneurysms , 2005 .