Boundary conditions by Schwarz-Christoffel mapping in anatomically accurate hemodynamics

Appropriate velocity boundary conditions are a prerequisite in computational hemodynamics. A method for mapping analytical or experimental velocity profiles on anatomically realistic boundary cross-sections is presented. Interpolation is required because the computational and experimental domains are seldom aligned. In the absence of velocity information one alternative is the adaptation of analytical profiles based on volumetric flux constraints. The presented algorithms are based on the Schwarz-Christoffel (S-C) mapping of singly or doubly connected polygons to the unit circle or an annulus with unary external radius. S-C transformations are combined to construct a one-to-one invertible map between the target surface and the measurement domain or the support of the source analytical profile. The proposed technique permits us to segment each space separately and map one onto the other in its entirety. Tests are performed with normal velocity boundary conditions for computational simulations of blood flow in the ascending aorta and cerebrospinal fluid flow in the spinal cavity. Mappings of axisymmetric velocity profiles of the Womersley type through a simply connected circular pipe as well as through a doubly connected circular annulus, and interpolations from in-vivo phase-contrast magnetic resonance imaging velocity measurements under instantaneous volumetric flux constraints are considered.