학술논문

Hierarchically refined and coarsened splines for moving interface problems, with particular application to phase-field models of prostate tumor growth.
Document Type
Journal
Author
Lorenzo, G. (E-CRNA-MM) AMS Author Profile; Scott, M. A. (1-BYU-CEV) AMS Author Profile; Tew, K. (1-BYU-IT) AMS Author Profile; Hughes, T. J. R. (1-TX-CPE) AMS Author Profile; Gomez, H. (1-PURD-SME) AMS Author Profile
Source
Computer Methods in Applied Mechanics and Engineering (Comput. Methods Appl. Mech. Engrg.) (20170101), 319, 515-548. ISSN: 0045-7825 (print).eISSN: 1879-2138.
Subject
65 Numerical analysis -- 65D Numerical approximation and computational geometry
  65D07 Splines

74 Mechanics of deformable solids -- 74L Special subfields of solid mechanics
  74L15 Biomechanical solid mechanics
Language
English
Abstract
Summary: ``Moving interface problems are ubiquitous in science and engineering. To develop an accurate and efficient methodology for this class of problems, we present algorithms for local $h$-adaptivity of hierarchical B-splines to be utilized in isogeometric analysis. We extend Bézier projection, an efficient quadrature-free local projection technique, to the hierarchical setting. In this case, extraction operators may not be invertible. To address this issue we develop a multi-level reconstruction operator which maintains the locality properties of the projection. We also introduce a balance parameter to control the overlap of hierarchical functions leading to improved numerical conditioning. We apply our algorithms to the simulation of localized prostate cancer growth. We model this disease using the phase-field method and a set of diffusion-reaction equations to account for the dynamics of nutrients and a key biomarker termed Prostate Specific Antigen. Our results include examples on simple 2D and 3D domains and a more compelling tissue-scale, patient-specific simulation, which is run over a prostate anatomy extracted from medical images. Our methods for local $h$-adaptivity efficiently capture the evolving interface between the tumor and the neighboring healthy tissue with remarkable accuracy in all cases.''