In vitro fluid dynamics of the Ahmed glaucoma valve modified with expanded polytetrafluoroethylene.

PURPOSE: Long-term intraocular pressure reduction by glaucoma drainage devices (GDDs) is often limited by the fibrotic capsule that forms around them. Prior work demonstrates that modifying a GDD with a porous membrane promotes a vascularized and more permeable capsule. This work examines the in vitro fluid dynamics of the Ahmed valve after enclosing the outflow tract with a porous membrane of expanded polytetrafluoroethylene (ePTFE). MATERIALS AND METHODS: The control and modified Ahmed implants (termed porous retrofitted implant with modified enclosure or PRIME-Ahmed) were submerged in saline and gelatin and perfused in a system that monitored flow (Q) and pressure (P). Flow rates of 1-50 μl/min were applied and steady state pressure recorded. Resistance was calculated by dividing pressure by flow. RESULTS: Modifying the Ahmed valve implant outflow with expanded ePTFE increased pressure and resistance. Pressure at a flow of 2 μl/min was increased in the PRIME-Ahmed (11.6 ± 1.5 mm Hg) relative to the control implant (6.5 ± 1.2 mm Hg). Resistance at a flow of 2 μl/min was increased in the PRIME-Ahmed (5.8 ± 0.8 mm Hg/μl/min) when compared to the control implant (3.2 ± 0.6 mm Hg/μl/min). CONCLUSIONS: Modifying the outflow tract of the Ahmed valve with a porous membrane adds resistance that decreases with increasing flow. The Ahmed valve implant behaves as a variable resistor. It is partially open at low pressures and provides reduced resistance at physiologic flow rates.

A New Method for Modeling Free Surface Flows and Fluid-structure Interaction with Ocean Applications

The computational modeling of ocean waves and ocean-faring devices poses numerous challenges. Among these are the need to stably and accurately represent both the fluid-fluid interface between water and air as well as the fluid-structure interfaces arising between solid devices and one or more fluids. As techniques are developed to stably and accurately balance the interactions between fluid and structural solvers at these boundaries, a similarly pressing challenge is the development of algorithms that are massively scalable and capable of performing large-scale three-dimensional simulations on reasonable time scales. This dissertation introduces two separate methods for approaching this problem, with the first focusing on the development of sophisticated fluid-fluid interface representations and the second focusing primarily on scalability and extensibility to higher-order methods.

We begin by introducing the narrow-band gradient-augmented level set method (GALSM) for incompressible multiphase Navier-Stokes flow. This is the first use of the high-order GALSM for a fluid flow application, and its reliability and accuracy in modeling ocean environments is tested extensively. The method demonstrates numerous advantages over the traditional level set method, among these a heightened conservation of fluid volume and the representation of subgrid structures.

Next, we present a finite-volume algorithm for solving the incompressible Euler equations in two and three dimensions in the presence of a flow-driven free surface and a dynamic rigid body. In this development, the chief concerns are efficiency, scalability, and extensibility (to higher-order and truly conservative methods). These priorities informed a number of important choices: The air phase is substituted by a pressure boundary condition in order to greatly reduce the size of the computational domain, a cut-cell finite-volume approach is chosen in order to minimize fluid volume loss and open the door to higher-order methods, and adaptive mesh refinement (AMR) is employed to focus computational effort and make large-scale 3D simulations possible. This algorithm is shown to produce robust and accurate results that are well-suited for the study of ocean waves and the development of wave energy conversion (WEC) devices.

Influence of Outlet Boundary Conditions on Cerebrovascular Aneurysm Hemodynamics

Computational fluid dynamic (CFD) studies of blood flow in cerebrovascular aneurysms have potential to improve patient treatment planning by enabling clinicians and engineers to model patient-specific geometries and compute predictors and risks prior to neurovascular intervention. However, the use of patient-specific computational models in clinical settings is unfeasible due to their complexity, computationally intensive and time-consuming nature. An important factor contributing to this challenge is the choice of outlet boundary conditions, which often involves a trade-off between physiological accuracy, patient-specificity, simplicity and speed. In this study, we analyze how resistance and impedance outlet boundary conditions affect blood flow velocities, wall shear stresses and pressure distributions in a patient-specific model of a cerebrovascular aneurysm. We also use geometrical manipulation techniques to obtain a model of the patient’s vasculature prior to aneurysm development, and study how forces and stresses may have been involved in the initiation of aneurysm growth. Our CFD results show that the nature of the prescribed outlet boundary conditions is not as important as the relative distributions of blood flow through each outlet branch. As long as the appropriate parameters are chosen to keep these flow distributions consistent with physiology, resistance boundary conditions, which are simpler, easier to use and more practical than their impedance counterparts, are sufficient to study aneurysm pathophysiology, since they predict very similar wall shear stresses, time-averaged wall shear stresses, time-averaged pressures, and blood flow patterns and velocities. The only situations where the use of impedance boundary conditions should be prioritized is if pressure waveforms are being analyzed, or if local pressure distributions are being evaluated at specific time points, especially at peak systole, where the use of resistance boundary conditions leads to unnaturally large pressure pulses. In addition, we show that in this specific patient, the region of the blood vessel where the neck of the aneurysm developed was subject to abnormally high wall shear stresses, and that regions surrounding blebs on the aneurysmal surface were subject to low, oscillatory wall shear stresses. Computational models using resistance outlet boundary conditions may be suitable to study patient-specific aneurysm progression in a clinical setting, although several other challenges must be addressed before these tools can be applied clinically.