TY - JOUR
T1 - Anatomically realistic lumen motion representation in patient-specific space–time isogeometric flow analysis of coronary arteries with time-dependent medical-image data
AU - Yu, Yuxuan
AU - Zhang, Yongjie Jessica
AU - Takizawa, Kenji
AU - Tezduyar, Tayfun E.
AU - Sasaki, Takafumi
N1 - Funding Information:
Y. Yu and Y.J. Zhang were supported in part by NSF CAREER Award OCI-1149591. This work was also supported (K. Takizawa and T. Sasaki) in part by JST-CREST; Grant-in-Aid for Scientific Research (A) 18H04100 from Japan Society for the Promotion of Science; and Rice–Waseda research agreement. It was also supported (T.E. Tezduyar) in part by Top Global University Project of Waseda University. Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Funding Information:
Y.?Yu and Y.J.?Zhang were supported in part by NSF CAREER Award OCI-1149591. This work was also supported (K.?Takizawa and T.?Sasaki) in part by JST-CREST; Grant-in-Aid for Scientific Research (A) 18H04100 from Japan Society for the Promotion of Science; and Rice?Waseda research agreement. It was also supported (T.E.?Tezduyar) in part by Top Global University Project of Waseda University.
Publisher Copyright:
© 2019, Springer-Verlag GmbH Germany, part of Springer Nature.
PY - 2020/2/1
Y1 - 2020/2/1
N2 - Patient-specific computational flow analysis of coronary arteries with time-dependent medical-image data can provide valuable information to doctors making treatment decisions. Reliable computational analysis requires a good core method, high-fidelity space and time discretizations, and an anatomically realistic representation of the lumen motion. The space–time variational multiscale (ST-VMS) method has a good track record as a core method. The ST framework, in a general context, provides higher-order accuracy. The VMS feature of the ST-VMS addresses the computational challenges associated with the multiscale nature of the unsteady flow in the artery. The moving-mesh feature of the ST framework enables high-resolution flow computation near the moving fluid–solid interfaces. The ST isogeometric analysis is a superior discretization method. With IGA basis functions in space, it enables more accurate representation of the lumen geometry and increased accuracy in the flow solution. With IGA basis functions in time, it enables a smoother representation of the lumen motion and a mesh motion consistent with that. With cubic NURBS in time, we obtain a continuous acceleration from the lumen-motion representation. Here we focus on making the lumen-motion representation anatomically realistic. We present a method to obtain from medical-image data in discrete form an anatomically realistic NURBS representation of the lumen motion, without sudden, unrealistic changes introduced by the higher-order representation. In the discrete projection from the medical-image data to the NURBS representation, we supplement the least-squares terms with two penalty terms, corresponding to the first and second time derivatives of the control-point trajectories. The penalty terms help us avoid the sudden unrealistic changes. The computation we present demonstrates the effectiveness of the method.
AB - Patient-specific computational flow analysis of coronary arteries with time-dependent medical-image data can provide valuable information to doctors making treatment decisions. Reliable computational analysis requires a good core method, high-fidelity space and time discretizations, and an anatomically realistic representation of the lumen motion. The space–time variational multiscale (ST-VMS) method has a good track record as a core method. The ST framework, in a general context, provides higher-order accuracy. The VMS feature of the ST-VMS addresses the computational challenges associated with the multiscale nature of the unsteady flow in the artery. The moving-mesh feature of the ST framework enables high-resolution flow computation near the moving fluid–solid interfaces. The ST isogeometric analysis is a superior discretization method. With IGA basis functions in space, it enables more accurate representation of the lumen geometry and increased accuracy in the flow solution. With IGA basis functions in time, it enables a smoother representation of the lumen motion and a mesh motion consistent with that. With cubic NURBS in time, we obtain a continuous acceleration from the lumen-motion representation. Here we focus on making the lumen-motion representation anatomically realistic. We present a method to obtain from medical-image data in discrete form an anatomically realistic NURBS representation of the lumen motion, without sudden, unrealistic changes introduced by the higher-order representation. In the discrete projection from the medical-image data to the NURBS representation, we supplement the least-squares terms with two penalty terms, corresponding to the first and second time derivatives of the control-point trajectories. The penalty terms help us avoid the sudden unrealistic changes. The computation we present demonstrates the effectiveness of the method.
KW - Anatomically realistic lumen motion
KW - Coronary arteries
KW - Patient-specific computational flow analysis
KW - Penalty spline
KW - Space–time isogeometric analysis
KW - Space–time variational multiscale method
KW - Time-dependent medical-image data
UR - http://www.scopus.com/inward/record.url?scp=85074482178&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85074482178&partnerID=8YFLogxK
U2 - 10.1007/s00466-019-01774-4
DO - 10.1007/s00466-019-01774-4
M3 - Article
AN - SCOPUS:85074482178
SN - 0178-7675
VL - 65
SP - 395
EP - 404
JO - Computational Mechanics
JF - Computational Mechanics
IS - 2
ER -