TY - JOUR
T1 - YZβ discontinuity capturing for advection-dominated processes with application to arterial drug delivery
AU - Bazilevs, Y.
AU - Calo, V. M.
AU - Tezduyar, T. E.
AU - Hughes, T. J.R.
PY - 2007/7/20
Y1 - 2007/7/20
N2 - The YZβ discontinuity-capturing operator, recently introduced in (Encyclopedia of Computational Mechanics, Vol. 3, Fluids. Wiley: New York, 2004) in the context of compressible flows, is applied to a time-dependent, scalar advection-diffusion equation with the purpose of modelling drug delivery processes in blood vessels. The formulation is recast in a residual-based form, which reduces to the previously proposed formulation in the limit of zero diffusion and source term. The NURBS-based isogeometric analysis method, proposed by Hughes et al. (Comput. Methods Appl. Mech. Eng. 2005; 194:4135-4195), was used for the numerical tests. Effects of various parameters in the definition of the YZβ operator are examined on a model problem and the better performer is singled out. While for low-order B-spline functions discontinuity capturing is necessary to improve solution quality, we find that high-order, high-continuity B-spline discretizations produce sharp, nearly monotone layers without the aid of discontinuity capturing. Finally, we successfully apply the YZβ approach to the simulation of drug delivery in patient-specific coronary arteries.
AB - The YZβ discontinuity-capturing operator, recently introduced in (Encyclopedia of Computational Mechanics, Vol. 3, Fluids. Wiley: New York, 2004) in the context of compressible flows, is applied to a time-dependent, scalar advection-diffusion equation with the purpose of modelling drug delivery processes in blood vessels. The formulation is recast in a residual-based form, which reduces to the previously proposed formulation in the limit of zero diffusion and source term. The NURBS-based isogeometric analysis method, proposed by Hughes et al. (Comput. Methods Appl. Mech. Eng. 2005; 194:4135-4195), was used for the numerical tests. Effects of various parameters in the definition of the YZβ operator are examined on a model problem and the better performer is singled out. While for low-order B-spline functions discontinuity capturing is necessary to improve solution quality, we find that high-order, high-continuity B-spline discretizations produce sharp, nearly monotone layers without the aid of discontinuity capturing. Finally, we successfully apply the YZβ approach to the simulation of drug delivery in patient-specific coronary arteries.
KW - Advection-diffusion equation
KW - Discontinuity capturing
KW - Drug delivery
KW - Fluids
KW - Interior layers
KW - Isogeometric analysis
KW - Navier-Stokes equations
UR - http://www.scopus.com/inward/record.url?scp=34347362891&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=34347362891&partnerID=8YFLogxK
U2 - 10.1002/fld.1484
DO - 10.1002/fld.1484
M3 - Article
AN - SCOPUS:34347362891
VL - 54
SP - 593
EP - 608
JO - International Journal for Numerical Methods in Fluids
JF - International Journal for Numerical Methods in Fluids
SN - 0271-2091
IS - 6-8
ER -