3D Printed Patient-Specific Model for In Vitro Hemodynamic Studies and Comparison with In Vivo Findings Using 4D Flow MRI
Rouzbeh R Ahmadian1, Austin P Boyd2, Jeremy D Collins1, James C Carr1, Alex J Barker1, and Michael Markl3

1Radiology, Northwestern University, Chicago, IL, United States, 2Northwestern University, Chicago, IL, United States, 3Radiology & Biomedical Engineering, Northwestern University, Chicago, IL, United States

Synopsis

The advent of 3D printing has opened exciting possibilities for biomedical applications. One of the most intriguing of these possibilities is the ability to use images obtained from radiology scanners (CT or MR) to create 3D models of patient anatomy followed by 3D printing. These models will have all the same geometries of patient anatomy to very high detail. In order to have practical applications, however, these 3D models need to behave similarly to human tissue under standard conditions. Our research utilizes patient-specific 3D printed models in an in vitro fluid dynamic circuit to compare 4D Flow MRI data to that of the patient (in vivo). This direct comparison will allow for validation of the 3D printed model for further biomedical application.

Purpose

Manipulation of vessel geometry via surgical intervention can have unpredictable consequences with regards to blood flow due to variation in specific patient anatomy. It is therefore imperative to devise a methodology to test the efficacy of treatments prior to surgical intervention. Studies have shown that 4D flow MRI can be used for detailed in-vivo assessment of vascular hemodynamics with 3D visualization and quantification of blood flow[1]. Advances in rapid prototyping (3D printing) have provided the ability to convert 3D imaging data into real patient-specific replicas quickly and at low cost[2]. This study was conducted to introduce a novel application of 3D printing to develop an in vitro circuit for blood flow studies. Using 4D flow MRI, hemodynamic characteristics in a 3D printed model of aortic coarctation is demonstrated and compared to in vivo findings.

Methods

3D contrast-enhanced (CE) MRA and 4D flow MRI were obtained from a 41-year-old patient with bicuspid aortic valve (BAV) and coarctation of the aorta (AoCo). Flow visualization and quantification was carried out using velocity mapping (Figure 1A-C).

DICOM files containing the MR results were imported into a segmentation software (Mimics®, Materialise). A digital model of the aorta was generated based on the 3D CE-MRA data and exported as an STL file (Figure 2A). The model was combined with standardized size 1” flange fittings at both ends for integration into the flow circuit (Figure 2B). A 1:1 replica made out of Acrylonitrile-Butadiene-Styrene (ABS) was generated using a 3D printer (Makerbot Replicator® 2X, max. build size 246 x 152 x 155 mm3, layer resolution 100µm) with settings at 0.2mm layer height, 75% infill and five-layer shells (Figure 2C). The model was then combined with vinyl and standard PVC pipe tubing (Figure 2D) to allow for integration into the experimental setup.

Schematic of the flow circuit is illustrated in Figure 3. A clinical Ventricular Assist Device-VAD[4],[5] (MEDOS VAD-Pump-Chamber, MEDOS Medizintechnik AG, Stolberg, Germany; chamber size 60 mL) was used in an MR-compatible flow circuit placed directly inside a 3 Tesla MR system (TRIO, Siemens, Germany). Pulsatile VAD flow was generated using pulsatile pneumatic pressure created by a VAD control unit (MEDOS). A triggering system controlled via National Instruments® Lab View application was used to simultaneously trigger the scanner and the VAD for synchronized gating. Fluid in the circuit was tap water containing 1.08 mmol/L Gadolinium to enhance signal. 4D flow MRI was conducted with 2mm3 voxel size.

Results

In vivo: velocity streamlines showed formation of flow jet patterns with a strong right-handed helix in both the ascending aorta (AAo) and distal to the coarctation (Figure 1B and C). Quantitative analysis revealed highest velocity through the coarctation (2.35 m/s) as well as a high velocity flow jet along the postero-lateral wall of the ascending aorta (DAo) secondary to BAV with fusion of right- and left-coronary leaflets.

In Vitro: 3D printing took approximately 5 hours for a total cost of $10US. Computer time for segmentation, 3D modeling and fitting design was about 90 minutes. 4D flow analysis of the in vitro study again included 3D flow visualization using 3D streamlines. In addition, flow quantification was performed using planes at various sections of the aortic model (Figure 4A). As demonstrated in Figures 4A-D flow characteristics inside the 3D printed model showed vortex formation in the arch and helical flow in the Dao, but no vortex in the AAo. Velocity and flow quantification showed peak velocity of 0.7 m/s and 155 ml/s respectively.

Discussion

Flow through the 3D printed patient-specific model seems to strongly correlate with the in vivo 4D flow MRI findings. This is indicative of how much blood vessel geometry can affect hemodynamics. As mentioned before, flow characteristics inside the 3D printed model had vortex formation in the arch and helical flow in the Dao, but no vortex in the AAo. We believe the absence of the AAo vortex is due to lack of BAV in the 3D model. Moreover, the velocities in the 3D printed model were drastically lower than in vivo findings. We believe this is due to low stroke volume of the VAD (60ml) compared to an adult human left ventricle (80-100ml). Lastly, the model itself can be enhanced by substituting flexible (compliant) plastics for rigid ABS as well as integrating the branches (subclavian and carotids). The results were encouraging that combined with 4D flow MRI, 3D printing might have the potential to enable better in-vitro simulation for intervention planning such as graft repair and testing of post-operative blood flow characteristics, which may facilitate procedure planning and surgical simulation.

Acknowledgements

Grant support by NIH K25HL119608, NIH R01HL115828, and the Northwestern University McCormick Engineering Research Catalyst Award

References

[1] Markl, M. Schnell, S. Barker, AJ. 4D Flow Imaging: Current Status to Future Clinical Applications. Curr Cardiol Rep. 2014; 16(5): 481.

[2] Marro, A. Bandukwala, T. Mak W. Three-Dimensional Printing and Medical Imaging: A Review of the Methods and Applications. Curr Probl Diagn Radiol. 2015; S0363-0188(15)00112-7.

[3] Allen BD, Baker AJ, Carr JC, Silverberg, RA, Markl M. Time-resolved three-dimensional phase contrast MRI evaluation of bicuspid aortic valve and coarctation of the aorta. Eur Heart J Cardiovasc Imaging. 2013 Apr; 14(4):3999.

[4] Lorenz R, Benk C, Bock J, Stalder AF, Kornik JG, Henning J, Markl M. Closed Circuit MR Compatible Pulsatile Pump System Using a Ventricular Assist Device and Pressure Control Unit. Magn Reson Med. 2012 Jan; 67(1):258-68.

[5] Benk C, Lorenz, R, Beyersdof F, Bock J, Klemm R, Korvink, JG Markl M. J Thorac Cardiovasc Surg 2011 Nov: 142(5):1019-26.

Figures

In Vivo 4D Flow MRI Findings. A. MRA of patient showing coarctation B. Right posterior view of the 4D flow MR streamlines C. LAO View. Note the vortex formation in both AAo and arch, as well as helical blow flow down the DAo after the coarctation.

Figure 2a. 3D Segmentation of aorta. Note the wide AAo opening and narrowing at the AoCo.

Figure 2b. 3D computer model with fittings. Note the shortening of DAo limb. The aortic branches are also fully sealed, leaving a complete aortic arch as the internal geometry.

Figure 2c. 3D printed model. Note the detailed fitting and small ridges on the surface indicating the extrusion methodology of the additive manufacturing.

Figure 2d. 3D printed model with fittings in situ of the in vitro model. The clear tubes are vinyl and used to connect the PVC piping in the rest of the circuit to the 3D printed model.

Figure 3. In Vitro Circuit in the MRI scanner. The triggering mechanism as well as control systems are shown. Note the connections to the water reservoir in order to complete the circuit for flow studies.

Figure 4A. Still image from 4D flow video indicating the

4B. Flow Characteristics in In Vitro Model. Side view (AAo on Left)

4C. Flow Characteristics in In Vitro Model. Anterior view (AAo on Left)

4D. Flow Characteristics in In Vitro Model. Posterior view (AAo on Right)



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
2718