A diffuse interface method for the Navier–Stokes/Darcy equations: Perfusion profile for a patient-specific human liver based on MRI scans

Stein K.F. Stoter1, Peter Müller1, Luca Cicalese2, Massimiliano Tuveri2, Dominik Schillinger1, Thomas J.R. Hughes3
1Department of Civil, Environmental, and Geo- Engineering, University of Minnesota, USA
2Texas Transplant Center and Hepatobiliary Surgery, University of Texas Medical Branch, USA
3Institute for Computational Engineering and Sciences, The University of Texas at Austin, USA

Tài liệu tham khảo

Beavers, 1967, Boundary conditions at a naturally permeable wall, J. Fluid Mech., 30, 197, 10.1017/S0022112067001375 Saffman, 1971, On the boundary condition at the interface of a porous medium, Stud. Appl. Math., 50, 93, 10.1002/sapm197150293 Jones, 1973, Low Reynolds number flow past a porous spherical shell, 231 Mikelic, 2000, On the interface boundary condition of Beavers, Joseph, and Saffman, SIAM J. Appl. Math., 60, 1111, 10.1137/S003613999833678X Discacciati, 2002, Mathematical and numerical models for coupling surface and groundwater flows, Appl. Numer. Math., 43, 57, 10.1016/S0168-9274(02)00125-3 Burman, 2007, A unified stabilized method for Stokes’ and Darcy’s equations, J. Comput. Appl. Math., 198, 35, 10.1016/j.cam.2005.11.022 Badia, 2009, Unified stabilized finite element formulations for the Stokes and the Darcy problems, SIAM J. Numer. Anal., 47, 1971, 10.1137/08072632X Chidyagwai, 2009, On the solution of the coupled Navier–Stokes and Darcy equations, Comput. Methods Appl. Mech. Engrg., 198, 3806, 10.1016/j.cma.2009.08.012 Layton, 2002, Coupling fluid flow with porous media flow, SIAM J. Numer. Anal., 40, 2195, 10.1137/S0036142901392766 Discacciati, 2007, Robin-Robin domain decomposition methods for the Stokes-Darcy coupling, SIAM J. Numer. Anal., 45, 1246, 10.1137/06065091X Chen, 2011, A parallel Robin-Robin domain decomposition method for the Stokes-Darcy system, SIAM J. Numer. Anal., 49, 1064, 10.1137/080740556 Baber, 2012, Numerical scheme for coupling two-phase compositional porous-media flow and one-phase compositional free flow, IMA J. Appl. Math., 77, 887, 10.1093/imamat/hxs048 Discacciati, 2009, Navier-Stokes/Darcy coupling: modeling, analysis, and numerical approximation, Rev. Mat. Complut., 22, 315, 10.5209/rev_REMA.2009.v22.n2.16263 Hughes, 2005, Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement, Comput. Methods Appl. Mech. Engrg., 194, 4135, 10.1016/j.cma.2004.10.008 Cottrell, 2009 Sethian, 2008, Solving partial differential equations on irregular domains with moving interfaces, with applications to superconformal electrodeposition in semiconductor manufacturing, J. Comput. Phys., 227, 6411, 10.1016/j.jcp.2008.03.001 Moumnassi, 2011, Finite element analysis on implicitly defined domains: An accurate representation based on arbitrary parametric surfaces, Comput. Methods Appl. Mech. Engrg., 200, 774, 10.1016/j.cma.2010.10.002 Joulaian, 2013, Local enrichment of the finite cell method for problems with material interfaces, Comput. Mech., 52, 741, 10.1007/s00466-013-0853-8 Sukumar, 2001, Modeling holes and inclusions by level sets in the extended finite-element method, Comput. Methods Appl. Mech. Engrg., 190, 6183, 10.1016/S0045-7825(01)00215-8 Fries, 2010, The generalized/extended finite element method: An overview of the method and its applications, Internat. J. Numer. Methods Engrg., 84, 253, 10.1002/nme.2914 Legrain, 2012, High order X-FEM and levelsets for complex microstructures: uncoupling geometry and approximation, Comput. Methods Appl. Mech. Engrg., 241, 172, 10.1016/j.cma.2012.06.001 Liehr, 2009, Composite finite elements for 3D image based computing, Comput. Vis. Sci., 12, 171, 10.1007/s00791-008-0093-1 Preusser, 2011, 3D Composite Finite Elements for Elliptic Boundary Value Problems with Discontinuous Coefficients, SIAM J. Sci. Comput., 33, 2115, 10.1137/100791750 Zhang, 2004, Immersed finite element method, Comput. Methods Appl. Mech. Engrg., 193, 2051, 10.1016/j.cma.2003.12.044 Löhner, 2008, Adaptive embedded and immersed unstructured grid techniques, Comput. Methods Appl. Mech. Engrg., 197, 2173, 10.1016/j.cma.2007.09.010 Schillinger, 2012, An isogeometric design-through-analysis methodology based on adaptive hierarchical refinement of NURBS, immersed boundary methods, and T-spline CAD Surfaces, Comput. Methods Appl. Mech. Engrg., 249–250, 116, 10.1016/j.cma.2012.03.017 Rüberg, 2014, A fixed-grid b-spline finite element technique for fluid–structure interaction, Internat. J. Numer. Methods Fluids, 74, 623, 10.1002/fld.3864 Parvizian, 2007, Finite cell method: h- and p- extension for embedded domain methods in solid mechanics, Comput. Mech., 41, 122, 10.1007/s00466-007-0173-y Schillinger, 2012, Small and large deformation analysis with the p- and B-spline versions of the Finite Cell Method, Comput. Mech., 50(4), 445, 10.1007/s00466-012-0684-z Schillinger, 2015, The finite Cell Method: A review in the context of higher-order structural analysis of CAD and image-based geometric models, Arch. Comput. Methods Eng., 22, 391, 10.1007/s11831-014-9115-y Kamensky, 2015, An immersogeometric variational framework for fluid–structure interaction: application to bioprosthetic heart valves, Comput. Methods Appl. Mech. Engrg., 284, 1005, 10.1016/j.cma.2014.10.040 Xu, 2016, The tetrahedral finite cell method for fluids: immersogeometric analysis of turbulent flow around complex geometries, Comput. & Fluids, 10.1016/j.compfluid.2015.08.027 Varduhn, 2016, The tetrahedral finite cell method: Higher-order immersogeometric analysis on adaptive non-boundary-fitted meshes, Int. J. Numer. Methods Engrg., 107, 1054, 10.1002/nme.5207 Stavrev, 2016, Geometrically accurate, efficient, and flexible quadrature techniques for the tetrahedral finite cell method, Comput. Methods Appl. Mech. Engrg., 310, 646, 10.1016/j.cma.2016.07.041 Fries, 2016, Higher-order accurate integration of implicit geometries, Internat. J. Numer. Methods Engrg., 106, 323, 10.1002/nme.5121 Kudela, 2016, Smart octrees: Accurately integrating discontinuous functions in 3D, Comput. Methods Appl. Mech. Engrg., 10.1016/j.cma.2016.04.006 Bueno-Orovio, 2006, Spectral methods for partial differential equations in irregular domains: the spectral smoothed boundary method, SIAM J. Sci. Comput., 28, 886, 10.1137/040607575 Rätz, 2006, PDE’s on surfaces—a diffuse interface approach, Commun. Math. Sci., 4, 575, 10.4310/CMS.2006.v4.n3.a5 Li, 2009, Solving PDEs in complex geometries: a diffuse domain approach, Commun. Math. Sci., 7, 81, 10.4310/CMS.2009.v7.n1.a4 K.Y. Lervåg, J. Lowengrub, Analysis of the diffuse-domain method for solving PDEs in complex geometries, arXiv preprint, 2014. arXiv:1407.7480 Teigen, 2009, A diffuse-interface approach for modeling transport, diffusion and adsorption/desorption of material quantities on a deformable interface, Commun. Math. Sci., 4, 1009 Elliott, 2011, Numerical computation of advection and diffusion on evolving diffuse interfaces, IMA J. Numer. Anal., 31, 786, 10.1093/imanum/drq005 Aland, 2010, Two-phase flow in complex geometries: A diffuse domain approach, Comput. Model. Eng. Sci., 57, 77 Teigen, 2011, A diffuse-interface method for two-phase flows with soluble surfactants, J. Comput. Phys., 230, 375, 10.1016/j.jcp.2010.09.020 Miehe, 2010, A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits, Comput. Methods Appl. Mech. Engrg., 199, 2765, 10.1016/j.cma.2010.04.011 Borden, 2012, A phase-field description of dynamic brittle fracture, Comput. Methods Appl. Mech. Engrg., 217–220, 77, 10.1016/j.cma.2012.01.008 Schillinger, 2015, Isogeometric collocation for phase-field fracture models, Comput. Methods Appl. Mech. Engrg., 284, 583, 10.1016/j.cma.2014.09.032 Mikelic, 2015, A phase-field method for propagating fluid-filled fractures coupled to a surrounding porous medium, SIAM Multiscale Model. Simul., 13, 367, 10.1137/140967118 Vilanova, 2013, Capillary networks in tumor angiogenesis: From discrete endothelial cells to phase-field averaged descriptions via isogeometric analysis, Int. J. Numer. Methods Biomed. Eng., 29, 1015, 10.1002/cnm.2552 Gomez, 2013, Three-dimensional simulation of unstable gravity-driven infiltration of water into a porous medium, J. Comput. Phys., 238, 217, 10.1016/j.jcp.2012.12.018 Anders, 2012, Computational modeling of phase separation and coarsening in solder alloys, Int. J. Solids Struct., 49, 1557, 10.1016/j.ijsolstr.2012.03.018 Liu, 2015, Liquid-vapor phase transition: Thermomechanical theory, entropy stable numerical formulation, and boiling simulations, Comput. Methods Appl. Mech. Engrg., 297, 476, 10.1016/j.cma.2015.09.007 Zhao, 2015, Isogeometric analysis of mechanically coupled Cahn–Hilliard phase segregation in hyperelastic electrodes of Li-ion batteries, Comput. Methods Appl. Mech. Engrg., 297, 325, 10.1016/j.cma.2015.09.008 Peskin, 2002, The immersed boundary method, Acta Numer., 11, 479, 10.1017/S0962492902000077 Griffith, 2005, On the order of accuracy of the immersed boundary method: Higher order convergence rates for sufficiently smooth problems, J. Comput. Phys., 208, 75, 10.1016/j.jcp.2005.02.011 Ramière, 2007, A fictitious domain approach with spread interface for elliptic problems with general boundary conditions, Comput. Methods Appl. Mech. Engrg., 196, 766, 10.1016/j.cma.2006.05.012 Ramière, 2007, A general fictitious domain method with immersed jumps and multilevel nested structured meshes, J. Comput. Phys., 225, 1347, 10.1016/j.jcp.2007.01.026 Maury, 2001, A fat boundary method for the Poisson problem in a domain with holes, J. Sci. Comput., 16, 319, 10.1023/A:1012821728631 Bertoluzza, 2011, Analysis of the fully discrete fat boundary method, Numer. Math., 118, 49, 10.1007/s00211-010-0317-4 Boyer, 1999, Mathematical study of multi-phase flow under shear through order parameter formulation, Asymptot. Anal., 20, 175 Boyer, 2002, A theoretical and numerical model for the study of incompressible mixture flows, Comput. & Fluids, 31, 41, 10.1016/S0045-7930(00)00031-1 Boyer, 2009, A local adaptive refinement method with multigrid preconditionning illustrated by multiphase flows simulations, 15 Karma, 1998, Quantitative phase-field modeling of dendritic growth in two and three dimensions, Phys. Rev. E, 57, 4323, 10.1103/PhysRevE.57.4323 Fenton, 2005, Modeling wave propagation in realistic heart geometries using the phase-field method, Chaos, 15, 013502, 10.1063/1.1840311 Shen, 2010, Numerical approximations of Allen–Cahn and Cahn–Hilliard equations, Discrete Contin. Dyn. Syst., 28, 1669, 10.3934/dcds.2010.28.1669 Donea, 2003 Schwen, 2014, Spatio-temporal simulation of first pass drug perfusion in the liver, PLoS Comput. Biol., 10, e1003499, 10.1371/journal.pcbi.1003499 Masud, 2002, A stabilized mixed finite element method for Darcy flow, Comput. Methods Appl. Mech. Engrg., 191, 4341, 10.1016/S0045-7825(02)00371-7 Malta, 2000, Numerical analysis of a stabilized finite element method for tracer injection simulations, Comput. Methods Appl. Mech. Engrg., 187, 119, 10.1016/S0045-7825(99)00113-9 Michler, 2013, A computationally efficient framework for the simulation of cardiac perfusion using a multi-compartment Darcy porous-media flow model, Int. J. Numer. Methods Biomed. Eng., 29, 217, 10.1002/cnm.2520 Payne, 1998, Analysis of the boundary condition at the interface between a viscous fluid and a porous medium and related modelling questions, J. Math. Pures Appl. (9), 77, 317, 10.1016/S0021-7824(98)80102-5 Jäger, 2001, On the roughness-induced effective boundary conditions for an incompressible viscous flow, J. Differential Equations, 170, 96, 10.1006/jdeq.2000.3814 Martin, 2000, Practical ray tracing of trimmed NURBS surfaces, J. Graphics Tools, 5, 27, 10.1080/10867651.2000.10487519 Netgen Mesh Generator, developed by J. Schoeberl, 2015, http://sourceforge.net/projects/netgen-mesher/ Schöberl, 1997, NETGEN. An advancing front 2D/3D-mesh generator based on abstract rules, Comput. Vis. Sci., 1, 41, 10.1007/s007910050004 Tezduyar, 1991, Stabilized finite element formulations for incompressible flow computations, Adv. Appl. Mech., 28, 1, 10.1016/S0065-2156(08)70153-4 Tezduyar, 2000, Finite element stabilization parameters computed from element matrices and vectors, Comput. Methods Appl. Mech. Engrg., 28, 411, 10.1016/S0045-7825(00)00211-5 Nguyen, 2017, Phase-field boundary conditions for the voxel finite cell method: surface-free stress analysis of CT-based bone structures, Int. J. Numer. Methods Biomed. Eng., 10.1002/cnm.2880 Chen, 2004, Generation, propagation, and annihilation of metastable patterns, J. Differential Equations, 206, 399, 10.1016/j.jde.2004.05.017 Carr, 1989, Metastable patterns in solutions of ut=e2uxx-f(u), Comm. Pure Appl. Math., 42, 523, 10.1002/cpa.3160420502 Ambrosio, 1990, Approximation of functional depending on jumps by elliptic functional via ?-convergence, Comm. Pure Appl. Math., 43, 999, 10.1002/cpa.3160430805 Bourdin, 1999, Image segmentation with a finite element method, ESAIM: Modélisation Mathématique et analyse Numérique, 33, 229, 10.1051/m2an:1999114 Braides, 2002 Bar, 2011, Mumford and Shah model and its applications to image segmentation and image restoration, 1095 Logg, 2012 Evans, 2013, Isogeometric divergence-conforming B-splines for the steady Navier-Stokes equations, Math. Models Methods Appl. Sci., 23, 1421, 10.1142/S0218202513500139 Bazilevs, 2007, Weak imposition of Dirichlet boundary conditions in fluid mechanics, Comput. & Fluids, 36, 12, 10.1016/j.compfluid.2005.07.012 Schillinger, 2016, The non-symmetric Nitsche method for the parameter-free imposition of weak boundary and coupling conditions in immersed finite elements, Comput. Methods Appl. Mech. Engrg., 309, 625, 10.1016/j.cma.2016.06.026 Wahlbin, 2006 Caiazzo, 2014, On classical iterative subdomain methods for the Stokes–Darcy problem, Comput. Geosci., 18, 711, 10.1007/s10596-014-9418-y Schneider, 2004, Preoperative assessment of liver function, Surg. Clin. North Am., 84, 355, 10.1016/S0039-6109(03)00224-X Hossain, 2012, Mathematical modeling of coupled drug and drug-encapsulated nanoparticle transport in patient-specific coronary artery walls, Comput. Mech., 49, 213, 10.1007/s00466-011-0633-2 Liu, 2012, Computational modeling of nanoparticle targeted drug delivery, Rev. Nanosci. Nanotechnology, 1, 66, 10.1166/rnn.2012.1014 Cattaneo, 2014, A computational model of drug delivery through microcirculation to compare different tumor treatments, Int. J. Numer. Methods Biomed. Eng., 30, 1347, 10.1002/cnm.2661 Ricken, 2010, A biphasic model for sinusoidal liver perfusion remodeling after outflow obstruction, Biomech. Model. Mechanobiol., 9, 435, 10.1007/s10237-009-0186-x Birkfellner, 2014 Debbaut, 2013