Parallel fast isogeometric L2 projection solver with GALOIS system for 3D tumor growth simulations

Abstract We present a parallel solver for isogeometric L2 projections implemented in GALOIS system. We use the solver for isogeometric finite element method simulations of tumor growth. First, we describe the dependencies between the quantities describing the melanoma tumor growth in three-dimensions. Later, we derive the system of partial differential equations modeling the tumor dynamics, including tumor cell density, tumor angiogenic factor, flux, pressure, extracellular and degraded extracellular matrices. The continuous model is coupled every ten-time steps with graph grammar model describing the vasculature network. The resulting system of PDEs is solved with our parallel solver over the shared memory Linux cluster node. We show high computational efficiency and accuracy of our model and we demonstrate its perfect parallel scalability with the number of CPU cores. We discuss the application of our solver in predictive oncology.

[1]  Keshav Pingali,et al.  IGA-ADS: Isogeometric analysis FEM using ADS solver , 2017, Comput. Phys. Commun..

[2]  T. Belytschko,et al.  A generalized finite element formulation for arbitrary basis functions: From isogeometric analysis to XFEM , 2010 .

[3]  David A. Yuen,et al.  A 3-D model of tumor progression based on complex automata driven by particle dynamics , 2009, Journal of Molecular Modeling.

[4]  M. Chaplain,et al.  Continuous and discrete mathematical models of tumor-induced angiogenesis , 1998, Bulletin of mathematical biology.

[5]  Witold Dzwinel,et al.  A concept of a prognostic system for personalized anti-tumor therapy based on supermodeling , 2017, ICCS.

[6]  M A Weinstock,et al.  The global burden of melanoma: results from the Global Burden of Disease Study 2015 , 2017, The British journal of dermatology.

[7]  H Rieger,et al.  Physical determinants of vascular network remodeling during tumor growth , 2010, The European physical journal. E, Soft matter.

[8]  Victor M. Calo,et al.  Isogeometric Variational Multiscale Large-Eddy Simulation of Fully-developed Turbulent Flow over a Wavy Wall , 2012 .

[9]  Yuri Bazilevs,et al.  High-performance computing of wind turbine aerodynamics using isogeometric analysis , 2011 .

[10]  Victor M. Calo,et al.  Computational cost estimates for parallel shared memory isogeometric multi-frontal solvers , 2014, Comput. Math. Appl..

[11]  T. Hughes,et al.  Isogeometric analysis of the isothermal Navier-Stokes-Korteweg equations , 2010 .

[12]  Oliver Weeger,et al.  Isogeometric Finite Element Analysis of Nonlinear Structural Vibrations , 2015 .

[13]  Thomas J. R. Hughes,et al.  Isogeometric Analysis for Topology Optimization with a Phase Field Model , 2012 .

[14]  Alessandro Reali,et al.  Explicit higher-order accurate isogeometric collocation methods for structural dynamics , 2018, Computer Methods in Applied Mechanics and Engineering.

[15]  S. McDougall,et al.  Mathematical modeling of tumor-induced angiogenesis. , 2006, Annual review of biomedical engineering.

[16]  M A Rupnick,et al.  Quantitative analysis of random motility of human microvessel endothelial cells using a linear under-agarose assay. , 1988, Laboratory investigation; a journal of technical methods and pathology.

[17]  Michael T. Heath.  A tale of two laws , 2015, Int. J. High Perform. Comput. Appl..

[18]  T. Hughes,et al.  Isogeometric analysis of the Cahn–Hilliard phase-field model , 2008 .

[19]  D. N. Herting,et al.  Finite elements: Computational aspects: Vol. III, by G.F. Carey and J. Tinsley Oden, Prentice-Hall, Englewood Cliffs, NJ, 1984 , 1985 .

[20]  Victor M. Calo,et al.  A finite strain Eulerian formulation for compressible and nearly incompressible hyperelasticity using high‐order B‐spline finite elements , 2012 .

[21]  Mark A. J. Chaplain,et al.  Mathematical modeling of cancer cell invasion of tissue: biological insight from mathematical analysis and computational simulation , 2011, Journal of mathematical biology.

[22]  Giancarlo Sangalli,et al.  Mathematical analysis of variational isogeometric methods* , 2014, Acta Numerica.

[23]  Witold Dzwinel,et al.  Application of fast isogeometric L2 projection solver for tumor growth simulations , 2017 .

[24]  Victor M. Calo,et al.  Mathematical modeling of coupled drug and drug-encapsulated nanoparticle transport in patient-specific coronary artery walls , 2012 .

[25]  Witold Dzwinel,et al.  Supermodeling in Simulation of Melanoma Progression , 2016, ICCS.

[26]  Victor M. Calo,et al.  Multiphysics model for blood flow and drug transport with application to patient-specific coronary artery flow , 2008 .