Ebrahim M. Kolahdouz

Ebrahim M. Kolahdouz, David Salac

Here a semi-implicit formulation of the gradient augmented level set method is presented. By tracking both the level set and it's gradient accurate subgrid information is provided,leading to highly accurate descriptions of a moving interface. The result is a hybrid Lagrangian-Eulerian method that may be easily applied in two or three dimensions. The new approach allows for the investigation of interfaces evolving by mean curvature and by the intrinsic Laplacian of the curvature. In this work the algorithm, convergence and accuracy results are presented. Several numerical experiments in both two and three dimensions demonstrate the stability of the scheme.

Ali Hasan, Ebrahim M. Kolahdouz, Andinet Enquobahrie et al.

Each year, approximately 300,000 heart valve repair or replacement procedures are performed worldwide, including approximately 70,000 aortic valve replacement surgeries in the United States alone. This paper describes progress in constructing anatomically and physiologically realistic immersed boundary (IB) models of the dynamics of the aortic root and ascending aorta. This work builds on earlier IB models of fluid-structure interaction (FSI) in the aortic root, which previously achieved realistic hemodynamics over multiple cardiac cycles, but which also were limited to simplified aortic geometries and idealized descriptions of the biomechanics of the aortic valve cusps. By contrast, the model described herein uses an anatomical geometry reconstructed from patient-specific computed tomography angiography (CTA) data, and employs a description of the elasticity of the aortic valve leaflets based on a fiber-reinforced constitutive model fit to experimental tensile test data. Numerical tests show that the model is able to resolve the leaflet biomechanics in diastole and early systole at practical grid spacings. The model is also used to examine differences in the mechanics and fluid dynamics yielded by fresh valve leaflets and glutaraldehyde-fixed leaflets similar to those used in bioprosthetic heart valves. Although there are large differences in the leaflet deformations during diastole, the differences in the open configurations of the valve models are relatively small, and nearly identical hemodynamics are obtained in all cases considered.

Guhan Velmurugan, Ebrahim M. Kolahdouz, David Salac

Many interfacial phenomena in physical and biological systems are dominated by high order geometric quantities such as curvature. Here a semi-implicit method is combined with a level set jet scheme to handle stiff nonlinear advection problems. The new method offers an improvement over the semi-implicit gradient augmented level set method previously introduced by requiring only one smoothing step when updating the level set jet function while still preserving the underlying methods higher accuracy. Sample results demonstrate that accuracy is not sacrificed while strict time step restrictions can be avoided.