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. Stoter; Peter Müller; Luca Cicalese; Massimiliano Tuveri; Dominik Schillinger; Thomas J.R. Hughes

doi:10.1016/j.cma.2017.04.002

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. Stoter, Peter Müller, Luca Cicalese, Massimiliano Tuveri, Dominik Schillinger, Thomas J.R. Hughes

Civil, Environmental, and Geo- Engineering

Research output: Contribution to journal › Article › peer-review

28 Scopus citations

Abstract

We present a diffuse interface method for coupling free and porous-medium-type flows modeled by the Navier–Stokes and Darcy equations. Its essential component is a diffuse geometry model generated from the phase-field solution of a separate initial boundary value problem that is based on the Allen–Cahn equation. Phase-field approximations of the interface and its gradient are then employed to transfer all interface terms in the coupled variational flow formulation into volumetric terms. This eliminates the need for an explicit interface parametrization between the two flow regimes. We illustrate accuracy and convergence for a series of benchmark examples, using standard low-order stabilized finite element discretizations. Our diffuse interface method is particularly attractive for coupled flow analysis on imaging data with complex implicit interfaces, where procedures for deriving explicit surface parametrizations constitute a significant bottleneck. We demonstrate the potential of our method to establish seamless imaging-through-analysis workflows by computing a perfusion profile for a full-scale 3D human liver based on MRI scans.

Original language	English (US)
Pages (from-to)	70-102
Number of pages	33
Journal	Computer Methods in Applied Mechanics and Engineering
Volume	321
DOIs	https://doi.org/10.1016/j.cma.2017.04.002
State	Published - Jul 1 2017

Bibliographical note

Funding Information:
D. Schillinger gratefully acknowledges support from the National Science Foundation through grant CISE-1565997 and the NSF CAREER Award No. 1651577. We also acknowledge the Minnesota Supercomputing Institute (MSI) of the University of Minnesota for providing computing resources that have contributed to the research results reported within this paper (https://www.msi.umn.edu/). We thank Dr. André Massing for his help with implementing conforming interfaces in FEniCS, and Dr. David Kamensky for insightful discussions on stabilization techniques.

Publisher Copyright:
© 2017 Elsevier B.V.

Keywords

Allen–Cahn equation
Diffuse interface method
Imaging data
Liver
Navier–Stokes/Darcy coupling
Phase-field approximation

Publisher link

10.1016/j.cma.2017.04.002

Find It

Find It at U of M Libraries

Cite this

@article{3f92e97306164f3e9b625bea43aa182f,

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

abstract = "We present a diffuse interface method for coupling free and porous-medium-type flows modeled by the Navier–Stokes and Darcy equations. Its essential component is a diffuse geometry model generated from the phase-field solution of a separate initial boundary value problem that is based on the Allen–Cahn equation. Phase-field approximations of the interface and its gradient are then employed to transfer all interface terms in the coupled variational flow formulation into volumetric terms. This eliminates the need for an explicit interface parametrization between the two flow regimes. We illustrate accuracy and convergence for a series of benchmark examples, using standard low-order stabilized finite element discretizations. Our diffuse interface method is particularly attractive for coupled flow analysis on imaging data with complex implicit interfaces, where procedures for deriving explicit surface parametrizations constitute a significant bottleneck. We demonstrate the potential of our method to establish seamless imaging-through-analysis workflows by computing a perfusion profile for a full-scale 3D human liver based on MRI scans.",

keywords = "Allen–Cahn equation, Diffuse interface method, Imaging data, Liver, Navier–Stokes/Darcy coupling, Phase-field approximation",

author = "Stoter, {Stein K.F.} and Peter M{\"u}ller and Luca Cicalese and Massimiliano Tuveri and Dominik Schillinger and Hughes, {Thomas J.R.}",

note = "Funding Information: D. Schillinger gratefully acknowledges support from the National Science Foundation through grant CISE-1565997 and the NSF CAREER Award No. 1651577. We also acknowledge the Minnesota Supercomputing Institute (MSI) of the University of Minnesota for providing computing resources that have contributed to the research results reported within this paper (https://www.msi.umn.edu/). We thank Dr. Andr{\'e} Massing for his help with implementing conforming interfaces in FEniCS, and Dr. David Kamensky for insightful discussions on stabilization techniques. Publisher Copyright: {\textcopyright} 2017 Elsevier B.V.",

year = "2017",

month = jul,

day = "1",

doi = "10.1016/j.cma.2017.04.002",

language = "English (US)",

volume = "321",

pages = "70--102",

journal = "Computer Methods in Applied Mechanics and Engineering",

issn = "0374-2830",

publisher = "Elsevier",

}

TY - JOUR

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

AU - Stoter, Stein K.F.

AU - Müller, Peter

AU - Cicalese, Luca

AU - Tuveri, Massimiliano

AU - Schillinger, Dominik

AU - Hughes, Thomas J.R.

N1 - Funding Information: D. Schillinger gratefully acknowledges support from the National Science Foundation through grant CISE-1565997 and the NSF CAREER Award No. 1651577. We also acknowledge the Minnesota Supercomputing Institute (MSI) of the University of Minnesota for providing computing resources that have contributed to the research results reported within this paper (https://www.msi.umn.edu/). We thank Dr. André Massing for his help with implementing conforming interfaces in FEniCS, and Dr. David Kamensky for insightful discussions on stabilization techniques. Publisher Copyright: © 2017 Elsevier B.V.

PY - 2017/7/1

Y1 - 2017/7/1

N2 - We present a diffuse interface method for coupling free and porous-medium-type flows modeled by the Navier–Stokes and Darcy equations. Its essential component is a diffuse geometry model generated from the phase-field solution of a separate initial boundary value problem that is based on the Allen–Cahn equation. Phase-field approximations of the interface and its gradient are then employed to transfer all interface terms in the coupled variational flow formulation into volumetric terms. This eliminates the need for an explicit interface parametrization between the two flow regimes. We illustrate accuracy and convergence for a series of benchmark examples, using standard low-order stabilized finite element discretizations. Our diffuse interface method is particularly attractive for coupled flow analysis on imaging data with complex implicit interfaces, where procedures for deriving explicit surface parametrizations constitute a significant bottleneck. We demonstrate the potential of our method to establish seamless imaging-through-analysis workflows by computing a perfusion profile for a full-scale 3D human liver based on MRI scans.

AB - We present a diffuse interface method for coupling free and porous-medium-type flows modeled by the Navier–Stokes and Darcy equations. Its essential component is a diffuse geometry model generated from the phase-field solution of a separate initial boundary value problem that is based on the Allen–Cahn equation. Phase-field approximations of the interface and its gradient are then employed to transfer all interface terms in the coupled variational flow formulation into volumetric terms. This eliminates the need for an explicit interface parametrization between the two flow regimes. We illustrate accuracy and convergence for a series of benchmark examples, using standard low-order stabilized finite element discretizations. Our diffuse interface method is particularly attractive for coupled flow analysis on imaging data with complex implicit interfaces, where procedures for deriving explicit surface parametrizations constitute a significant bottleneck. We demonstrate the potential of our method to establish seamless imaging-through-analysis workflows by computing a perfusion profile for a full-scale 3D human liver based on MRI scans.

KW - Allen–Cahn equation

KW - Diffuse interface method

KW - Imaging data

KW - Liver

KW - Navier–Stokes/Darcy coupling

KW - Phase-field approximation

UR - http://www.scopus.com/inward/record.url?scp=85018410925&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85018410925&partnerID=8YFLogxK

U2 - 10.1016/j.cma.2017.04.002

DO - 10.1016/j.cma.2017.04.002

M3 - Article

AN - SCOPUS:85018410925

SN - 0374-2830

VL - 321

SP - 70

EP - 102

JO - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

ER -