- 1BioCardioLab, Bioengineering Unit - Heart Hospital, Fondazione Toscana “G. Monasterio”, Massa, Italy
- 2Philips Research Paris, Suresnes, France
- 3Mines Saint-Étienne, Université Jean Monnet, INSERM, Saint-Étienne, France
- 4Predisurge, Grande Usine Creative 2, Saint-Etienne, France
- 5Department of Radiology, Fondazione Toscana “G. Monasterio”, Massa, Italy
Background: Abdominal Aortic Aneurysm (AAA) is a balloon-like dilatation that can be life-threatening if not treated. Fabricating patient-specific AAA models can be beneficial for in-vitro investigations of hemodynamics, as well as for pre-surgical planning and training, testing the effectiveness of different interventions, or developing new surgical procedures. The current direct additive manufacturing techniques cannot simultaneously ensure the flexibility and transparency of models required by some applications. Therefore, casting techniques are presented to overcome these limitations and make the manufactured models suitable for in-vitro hemodynamic investigations, such as particle image velocimetry (PIV) measurements or medical imaging.
Methods: Two complex patient-specific AAA geometries were considered, and the related 3D models were fabricated through material casting. In particular, two casting approaches, i.e. lost molds and lost core casting, were investigated and tested to manufacture the deformable AAA models. The manufactured models were acquired by magnetic resonance, computed tomography (CT), ultrasound imaging, and PIV. In particular, CT scans were segmented to generate a volumetric reconstruction for each manufactured model that was compared to a reference model to assess the accuracy of the manufacturing process.
Results: Both lost molds and lost core casting techniques were successful in the manufacturing of the models. The lost molds casting allowed a high-level surface finish in the final 3D model. In this first case, the average signed distance between the manufactured model and the reference was () mm. However, this approach was more expensive and time-consuming. On the other hand, the lost core casting was more affordable and allowed the reuse of the external molds to fabricate multiple copies of the same AAA model. In this second case, the average signed distance between the manufactured model and the reference was () mm. However, the final model’s surface finish quality was poorer compared to the model obtained by lost molds casting as the sealing of the outer molds was not as firm as the other casting technique.
Conclusions: Both lost molds and lost core casting techniques can be used for manufacturing patient-specific deformable AAA models suitable for hemodynamic investigations, including medical imaging and PIV.
1. Introduction
An abdominal aortic aneurysm (AAA) is a vascular disease that occurs when the arterial lumen dilates permanently and irreversibly, and the aortic diameter exceeds 30 mm (1). Currently, there are flaws in fully understanding the mechanisms of AAA initiation, development, and rupture, as well as in managing the treatment (2–5). In the biomedical field, extensive research has been performed throughout the last two decades to have insight into the AAA disease by exploiting both in-silico and in-vitro approaches as an effective alternative to animal experiments (6). Indeed, advances in medical imaging have enabled the acquisition of high-resolution images of the anatomy of interest, from which a 3D model can be accurately reconstructed and manufactured (7) for both planning and training (8–15).
In the context of functional investigations, the insertion of 3D models in a mock-circulatory loop (MCL) system provides a useful tool to reproduce and study the cardiovascular system at a patient-specific level in terms of geometry and hemodynamic conditions (16–18). Furthermore, imaging models inserted in MCL systems by means of magnetic resonance imaging (MRI) (19, 20), ultrasound (US) (21–23), and particle imaging velocimetry (PIV) (24) can deepen the knowledge of cardiovascular diseases, providing insights into the pathology conditions and progression, also in combination with in-silico studies (25, 26). However, each imaging modality requires specific features of the manufacturing material. Indeed, MRI applications require nonmagnetic materials, US imaging requires specific values of sound speed, acoustic attenuation, and echogenicity of the material (27–30), whereas models made of optically clear materials with the same refractive index as the blood-mimicking fluid are required for PIV applications (31). Either rigid or deformable models can be manufactured based on the application and purpose of the study. Glass and transparent resins are examples of materials used to manufacture rigid models, whereas rubberlike compounds, polyurethane, and silicone are examples of materials that can be used to fabricate deformable models.
Focusing on AAA 3D models, several works in the literature involved a variety of geometries, materials, manufacturing techniques, and applications. The geometry of AAA phantoms can be either idealized (32–35), which allows performing parametric studies and investigating complex hemodynamics in controlled conditions, or patient-specific (12, 36–38) to observe personalized morphological and hemodynamic features. Mock circulatory loops coupled with MRI systems were used to measure flows and wall shear stress (35), or assess temporal changes in intraluminal pressure and in the AAA sac (39). A model compatible with computerized tomography (CT) and US imaging was also manufactured and its utility was demonstrated during a simulated EVAR procedure (40). Finally, many studies performed PIV measurements to provide accurate quantification of the abnormal velocity fields in AAA by using rigid (32, 41, 42) or deformable silicone (34, 43–47) models.
In order to couple echogenicity, transparency, and deformability of the manufactured models, silicone-based materials have been proven to be the most suitable to manufacture complex anatomical models (30). In this regard, while the material is selected based on the type of application, the fabrication method is selected based on the material. Currently, silicone-based materials are challenging to directly print a model (48). Consequently, material casting techniques are preferable (49).
In material casting, silicone-based materials are poured into a mold, usually 3D-printed, which contains a hollow cavity of the desired shape where poured materials cure and solidify. Once the curing phase is complete, the mold is removed, thus obtaining the final manufactured model. In the case of hollow and thin-walled models, both inner (or core) and outer molds are required. Fabricating molds and deformable models is challenging in the case of complex geometries that require the use of sacrificial molds to fabricate the final model. In these complex cases, the core mold is always sacrificed, whereas the outer molds can be sacrificed (lost molds technique) or preserved (lost core technique) and reused for more castings.
However, despite a quite large production of these types of works, there is a significant lack of accurate information on the manufacturing of deformable phantoms, thus allowing the reproducibility and repeatability of the manufacturing process. This is particularly true in the case of 3D model fabrication through molding techniques, which are user-dependent and require multiple steps (3D printing, mold definition, mold assembly, material pouring, demolding) with different aspects to consider that can influence the success or the failure as well as the final results of the manufacturing process.
In this work, we present a detailed description of all the phases necessary to create complex (when undercuts are present) deformable CT-based AAA models to be inserted in a mock-circulatory loop system for image-based hemodynamic investigations. Both lost molds and lost core techniques are investigated and used to manufacture the AAA models. All design and fabrication phases are described, along with costs and equipment requirements, addressing the main advantages and limitations of the proposed methodologies.
2. Materials and methods
The manufacturing procedure to fabricate deformable models that mimic the compliant behavior of AAAs involved six main steps: (i) model creation; (ii) design of the molds; (iii) 3D-printing of the molds; (iv) surface treatment; (v) material casting; (vi) removal of the molds. Finally, the manufacturing accuracy was assessed. Each of these steps is detailed in the following sections.
2.1. Model creation
2.1.1. Image dataset and segmentation
CT images of the thoracoabdominal anatomy of 85 and 78-year-old patients affected by AAA were used to create the 3D models. The datasets were provided by Rigshospitalet (Copenhagen, DK), which recruited patients from the Copenhagen Aneurysms Cohort (COACH), a cohort of patients with AAAs under surveillance with advanced ultrasound imaging and biochemistry (5).
The images were segmented and the volume of the vessels was reconstructed with the open-source software 3DSLICER (Slicer, NIH) (50) by following the approach previously described in Celi et al. (51). Figures 1A–C show an example of segmentation in the axial, coronal, and longitudinal planes. The segmented AAA geometry is shown in Figure 1D.
Figure 1. Example of image segmentation in axial (A), coronal (B), longitudinal view (C), and segmentation result (D); 3D solid models (E, G) with related cross-sections (Sect) (F, H).
2.1.2. AAA solid model generation
The segmented geometries were exported as stereo-lithography format (stl) files and imported into the solid modeling computer-aided design (CAD) software SpaceClaim (ANSYS, Inc.) where the associated solid models were generated by means of non-uniform rational basis-splines interpolation. The 3D solid models are shown in Figures 1E,G. Both models include the abdominal aortic inlet (AA), aneurysm bulge, Left Iliac artery (LI), and Right Iliac artery (RI), which are the two outlets. The solid models were then imported into the CAD software SOLIDWORKS (Dassault Systèmes S. A., Vélizy-Villacoublay, France). Here, outward offset values of 1.8 and 2 mm were applied to the AAA models, respectively, to fix the thickness as shown by the cross-sections in Figures 1F,H.
Among all those available, we chose these two patient-specific geometries because they are representative of challenges in the casting technique due to their shape complexity. The two models were manufactured by testing two techniques: a lost molds technique for the first AAA model, and a lost core technique for the second one, as described in the following sections. According to the casting technique, the two models were referred to as lost molds model (LMM) and lost core model (LCM).
2.2. Design of the molds
The core and outer molds were designed in SOLIDWORKS for both cases by following the same procedure as the designing of the molds does not depend on the casting technique.
2.2.1. Core molds
Cores define the inner surface of the deformable models. As the first step, a hollow core was designed for each model by extracting the inner wall of the AAA solid models (in Figures 1E,G) and applying an inward offset value of 3 mm to define the thickness of the core. The LI and RI outlets and the AA inlets were closed and extruded 15 mm outward and inward, respectively, to obtain reference pins. Additionally, a 15 mm cylindrical pin was added to the inlets. Those pins allowed a firm alignment and a reference between the core and outer molds, thus maintaining the correct gap for material casting. Figure 2 shows the extrusion and reference pins related to the core of LMM.
2.2.2. Outer molds
The outer molds define the outer surface of the manufacturing deformable models. For both LMM and LCM, outer molds were designed through a multi-steps strategy consisting of:
(1) Creation of a monolithic external mold - An outward offset value of 10 mm was applied to the outer surface of each AAA solid model (Figures 1E,G) after closing inlet and outlets. A monolithic outer mold, solid inside, was obtained. This mold can contain the core and AAA solid model.
(2) Splitting of the outer molds - Due to the complex anatomies, it was necessary to create a specific separation surface for each model to avoid undercuts and to facilitate the assembly. For this reason, a specific separation surface was created for each model and used as the trim tool.The AAA solid models (Figures 3A and 4A) were considered to generate the separation surface for each model. Furthermore, since non-conventional separation surfaces were employed instead of planes, we divided the outer mold appropriately to satisfy the requirement to assemble and dismantle the several components into which the outer mold would be separated, thereby ensuring the right gap between the core and outer molds. First, the pull directions and a draft angle of 4 were selected to create a separation surface using the parting line and parting surfaces SOLIDWORKS tools. Specifically, the pull directions are the directions in which the two halves of a mold separate. In these cases, the pull directions were defined by the planes in Figures 3A and 4A. The so-obtained surfaces are shown in Figures 3B and 4B together with the AAA solid models. These surfaces were further extended outward since a wide enough surface is necessary to split the external mold (Figures 3C and 4C). We hence split the external mold by using those separating surfaces as the trim tools and obtaining the two parts of Figures 3D and 4D. Referring to the LCM only, there is a discontinuity in the separating surface (as shown in Figure 4C) that would have made the use of only two external molds impossible. As a result, this region was omitted from the two outer molds, and a third mold was created (in purple in Figure 4D), which is assembled with the other two orthogonally to the pull direction.
(3) Subtraction of the core and AAA solid model from the outer molds - From solid, the outer molds were made hollow by subtracting the core and AAA solid model from the outer molds. Therefore, the core and the outer molds could be assembled and disassembled successfully while ensuring the correct internal gap between them, which corresponds to 1.8 and 2 mm (thickness of models) for LCM and LMM, respectively.
(4) Designing of sprues, reference pins for the core, and gluing channels. We cut the top of the outer molds by removing material between two concentric circumferences so as to recreate a conical-shaped, 20-mm-high sprue (in Figures 5A,B for the LMM) that would facilitate the material pouring and prevent it from spilling out of the molds. We then added four 90- equispaced pins on the inner surface of the sprue (as shown in Figures 5A,B for the LMM) that could serve as a support and reference for core placement. As a final step, we drew a spline that could internally follow the contour of one of the two outer molds. This spline defined the path to sweep a circular profile to create a groove on one outer mold and a complementary extrusion on the other as shown in detail in Figures 5C,D for the LMM. Radial and axial geometric tolerances were assigned to the outer molds in correspondence with the reference pins to ensure assembly.
Figure 3. Outer mold splitting procedure for LMM: AAA solid model and plane used to define the pull direction (A) and parting surface (B); extrusion of the parting surface outward (C); outer molds after splitting (D).
Figure 4. Outer mold splitting procedure for LCM: AAA solid model and plane used to define the pull direction (A) and parting surface (B); extrusion of the parting surface outward and detail of surface discontinuity (C); outer molds after splitting (D).
Figure 5. Example of outer molds for LMM and related sprue, reference pins (A, B), and gluing channels (C, D).
2.3. 3D printing of the molds
Two casting techniques were tested: (i) lost molds casting, meaning that both core and outer molds were printed in polyvinyl alcohol (PVA), which is a water-soluble material, and dissolved once the casting material has cured; (ii) lost core casting, meaning that only the core was printed in PVA, whereas the outer molds were printed in acrylonitrile–butadiene–styrene (ABS) so that they may be reused to manufacture the same model several times. Core and outer molds were 3D printed with a fused deposition modeling (FDM) printer (A4v3, 3NTR Italy). The G-Code format files were generated with the dedicated software 3NTR SSI (3NTR, Italy).
2.3.1. Lost molds
A PVA filament of 2.85 mm was used as the primary material for printing the molds and as support material. Rectilinear and honeycomb patterns were set for the support and infill, respectively. The support density was set to medium for the outer and to fine for the core molds. The main settings of PVA printing are detailed in Table 1.
2.3.2. Lost core
A PVA filament of 2.85 mm was used as the primary material for printing the core and as support material by setting a fine density. An ABS filament of 2.85 mm was used as the primary material for printing the outer molds, whereas the supports were printed in the high-impact polystyrene-based SSU00 (3NTR, Italy) material. Rectilinear patterns were set for the support as well as the infill. The main settings of ABS printing are detailed in Table 1.
2.4. Surface treatment
Prior to casting, supports were manually separated from the molds. The external surfaces of the cores and the internal surfaces of the outer molds (contact surfaces) were treated to eliminate surface roughness typical of FDM printing. Specifically, these surfaces were first carefully manually sanded with fine-grained sandpaper. Subsequently, we treated cores and outer molds depending on the lost molds or lost core casting technique.
2.4.1. Lost molds
A solution of 15% (by weight) of low molecular weight PVA powder was prepared in distilled water at 120C. This solution was stirred until the solute was completely dissolved, resulting in a transparent solution that was then evenly brushed twice onto the contact surfaces, alternating between brushing and air drying. This process made the contact surfaces smooth, ensuring a good finish to the final model. Finally, those surfaces were treated with a release agent (Mann, Macungie, Pennsylvania) to ensure that the casting material was completely released from the molds during molds’ dissolving. Once the contact surfaces were completely dry, we assembled the core, poured liquid vinyl glue into the gluing channels, and closed the outer mold with its counterpart, ensuring a firm press by means of ties. Acrylic glue was added at the joint edges of the outer molds, and some tape was affixed so as to prevent possible leakage of the casting material. The molds were then left in place for 12 h.
2.4.2. Lost core
The contact surfaces of each mold were treated with clear acrylic paint (Ambro-Sol S.r.l., Italy) to make those surfaces hydrophobic thus avoiding the interaction between PVA and casting material. This interaction could alter the optical properties of the phantom. As for the lost molds casting, the surfaces adhering to the casting material were treated with a release agent (Mann, Macungie, Pennsylvania) to ensure that the casting material was completely released from the molds during external mold disassembling and core dissolving. Once the treated surfaces were completely dry, the molds were assembled and sealed by means of a two-component sealant specific for ABS (presto, European Aerosols, Wolvega, the Netherlands), applied to the joint edges of the outer molds so as to prevent possible leakage of the casting material. This specific type of sealant was employed since it can both provide a firm seal and be easily removed without destroying the molds as an acrylic glue would. Therefore, we waited until the sealant was completely dry before casting the material into the molds.
2.5. Material casting
Patient-specific AAA deformable models were made out of the two-component (silicone and curing agent) silicone material Sylgard 184 (Dow, Wiesbaden, Germany). Properties of the used material can be found in Table 2 (31).
According to the instructions, we prepared the material with a 10:1 (silicone:curing agent) ratio. After stirring, the casting material was degassed in a vacuum chamber for 30 min at 0.8 bar to remove air bubbles resulting from the mixing of the two parts. Degassing cycles were repeated until the total removal of bubbles was achieved. Subsequently, the material was poured into the molds in several batches to prevent air bubbles being trapped trapped in thick layers. For this reason, several cycles of pouring/degassing at 0.8 bar for 30 min were alternated. We concluded the casting phase once the silicone had reached the upper edge of the inner mold, and double-checked that there were no visible air bubbles. The silicone was cured for 48 h at room temperature.
2.6. Removal of the molds
2.6.1. Lost molds
Once the silicone has cured, the PVA components were dissolved by placing the entire assembly of molds in an ultrasonic tank filled with warm water and by securing the assembly to the bottom. 3 to 4 washing cycles were done, the first one about 1 h long at 60C with ultrasound enabled to saturate the PVA faster. The glued parts of the molds were broken up by using small tools. The PVA dissolved leaving both a smooth lumen inside of the silicone and a smooth surface outside of the silicone. The LMM was finally left in cold water for a few hours.
2.6.2. Lost core
Once the silicone has cured, the sealant for ABS was manually removed by sanding the zones of interest and then scratching the sealant using a scalpel. First, the purple component of Figure 4, which is orthogonal to the pull direction was disassembled. One out of two outer molds was also disassembled whereas the remaining outer mold, the core, and the silicone model were placed in an ultrasonic tank filled with warm water with enabled ultrasound. The outer mold was removed after about 20 min and the PVA core dissolved in about 4 h at 70C. The LCM was finally left in cold water for a few hours.
2.7. Accuracy assessment
The accuracy of the manufactured models was assessed by comparing them and the original CAD models in terms of related meshes. The manufactured models were scanned with CT to generate image volumes comparable to the patient ones. The images acquired on each model were then segmented to reconstruct the so-called manufactured mesh (also including the wall thickness) that was compared to the mesh of the original CAD model (the so-called reference mesh). As the first step, the two meshes, which had a comparable number of vertices, were registered through iterative closest point algorithm (52) performing a rigid registration. Once registered, the signed distance was computed for each vertex of the compared mesh with respect to the reference mesh by using CloudCompare (v.2.13.alpha [GPL software], 2022) (53). The generated distance scalar field was hosted by the vertices of the compared mesh. The same procedure was performed for both LMM and LCM.
3. Results
3.1. 3D printing of the molds
Figures 6A–C and 7A–C show the CAD representation of the several components designed to fabricate the LMM and the LCM, respectively. Here, the whole assemblies, which also include the AAA solid models (in red in Figures 6C and 7C), are also given. Symmetrically, Figures 6D–F and 7D–F show the results of the mold printing via FDM for the LMM and LCM, respectively. The FDM printer accurately created outer and core molds, coherently with the original design. In particular, all the complexities of these patient-specific geometries were reproduced. The core and outer molds were easy to assemble, achieving a perfect closure while maintaining the correct inner gap into which the casting material would later creep as shown in Figures 6E and 7E. The sealing of the outer molds using PVA and acrylic glue or the sealant specific for ABS is shown in Figures 6F and 7F for the lost molds and the lost core casting techniques, respectively.
3.2. Material casting and cost analysis
The silicone castings succeeded in both lost molds and lost core casting as shown in Figure 8. The PVA molds dissolved completely without interacting with the casting material. The internal and external surfaces of both silicone models are smooth. Small outgrowths of material are present in the LCM at the joint edges of the three outer molds. This defect is not noticeable in the LMM as the PVA molds allowed for a firmer sealing, and there was no silicone leakage from the outer molds. The main features of the manufactured models are summarized in Table 3, including material costs, necessary equipment, and compatibility with medical imaging modalities.
3.3. Medical imaging applications and accuracy assessment
The manufactured model can be inserted into an MCL system to replicate patient-specific hemodynamic conditions and/or perform PIV investigations. In this case, models are often filled with a blood-mimicking fluid and placed within a purpose-designed box, which may then be filled with water, the same blood-mimicking fluid (in PIV), or gelatin (in US imaging). Figures 9A,B shows LMM immersed in water and a blood-mimicking fluid (mixture of water, glycerol, and urea) that matches the refractive index of Sylgard 184. In the latter case, the models become transparent.
As reported in Table 3, the models are also compatible with the main medical imaging modalities. Examples of 3D rendering of MRI and CT scans acquired on the LMM are shown in Figures 10A,B, respectively. An example of Doppler B-mode acquisition of the same model in longitudinal view is depicted in Figure 10C. Finally, an example of illuminated particles for PIV is shown in Figures 10D,E. Portions proximal to LMM’s inlet were illuminated to reveal no particle distortion and no light scattering at the solid/liquid interface.
Figure 10. Render of MRI (A) and CT (B) scans acquired on LMM; B-mode Doppler acquisition in longitudinal view (C); PIV images (D–E).
CT scans of the two manufactured models were also used to assess the accuracy of the manufacturing with respect to the original design considered as reference. The signed distance between the manufactured and reference meshes for LMM and LCM is depicted in Figures 11A,B, respectively. This distance takes into account differences in thickness, is expressed in [mm], and is depicted as the spatial distribution of the distance scalar field with the related histograms. Values of signed distance expressed as deviation equal to () mm and () were computed for LMM and LCM, respectively. In both cases, the regions with the highest distance to the reference are inlets and outlets. This trend is more evident in the LCM and may be caused by a mild misalignment of the core with respect to the outer molds, resulting in wall thickening (red region of Figure 11B) on one side and wall thinning (blue region of Figure 11B) on the other. Nevertheless, these inaccuracies here are not an issue because the AAA models will be plugged into the experimental rig right in these regions.
Figure 11. Signed distance [mm] computed for LMM (A) and for LCM (B). Histograms of distances are depicted below color bars.
4. Discussion
In the context of the mock-circulatory loop systems, fabricating deformable models represents a crucial step for in-vitro quantitative analysis and device testing. This work presented all the phases to manufacture complex patient-specific AAA models with deformable material through lost molds and lost core casting techniques. In particular, a step-by-step workflow is presented to make the manufacturing of these types of models reproducible and repeatable. In this work, two different geometries have been considered to include two common shapes: a shorter and a longer AAA. However, it is worth highlighting that the presented workflow can be applied to several anatomic configurations, in particular when separation planes are difficult to be defined and/or undercuts are present. Aside from the materials used to print the outer molds and the casting technique, the two models differ in anatomical complexity, size, aneurysmatic bulge shape, and thickness of the final silicone model.
The segmentation of the patient’s images was the first step for manufacturing the patient-specific models. However, this step is user-dependent and can be a source of uncertainty in the resulting 3D-printed model (7). Dealing with patient-specific models, from image segmentation or from statistical shape modelling (54), the CAD phase of the molds necessitated some experience. The complex morphologies presented in this work required the definition of non-conventional separation surfaces, instead of planes, in the mold design phase to ensure correct phantom extraction. Indeed, as many undercuts may occur in such molds, it is not possible to define mold separation from a continuous surface since some regions of the geometry could penetrate the molds during the phantom extraction. In this case, additional separation surfaces must be defined by the user. This aspect and all the phases necessary to identify the separation geometries (including the software commands) and overcome undercuts were described in our work to allow the reproducibility of the fabrication and reduce user dependency.
PVA was chosen due to its hydrosolubility property. However, it is almost four times more expensive than ABS and is highly hygroscopic. For this reason, it was necessary to ensure that molds printed in PVA would not undergo deformation before assembly/material casting caused by the PVA’s sensitivity to humidity. ABS is a plastic material with superior mechanical and thermal properties than conventional plastics, which could also imply greater printing difficulties. Among all the properties, acrylonitrile, butadiene, and styrene provide hardness, rigidity, mechanical strength, impact resistance, and resistance to chemical products to ABS-printed components. ABS can be post-processed (polishing, drilling, gluing, or painting) easily, without special tools or products. Finally, ABS is compatible with silicone. These features determined the choice of ABS to print the outer molds for the LCM. Using ABS instead of PVA for printing outer molds is more advantageous in applications that require many models of the same geometry and is more affordable. In this regard, the FDM printer deposed 623 g of materials for LCM and 227 g for LMM. By normalizing the cost of materials used for FDM for each procedure, the LCM technique costs 43.98 €/kg, including the cost of PVA core, ABS outer molds, and SSU00 supports for printing ABS. The LMM technique costs 121 €/kg as only PVA was used.
A 20% infill and two contours were an acceptable compromise between strength and ease of dissolution of the molds for the LMM. In this case, a fine support density was set for the core to avoid distortions due to PVA hygroscopicity. A 16% infill and 4 contours were instead set to print both outer and core molds for the LCM. Reducing the infill and increasing the contours, silicone had fewer chances to penetrate inside the outer and core, but mold melting was faster now. This parameter variation did not affect the strength and stiffness of the molds. The FDM printer has more nozzles, allowing the simultaneous printing of several materials. For this reason, it was possible to print the supports of the ABS molds in SSU00, allowing for a simpler detachment of the molds from their supports, as well as superior surface quality in the contact regions due to no scratches or roughness following removal. This aspect was important at spots where the molds would come into contact during assembly, such as the gluing channels. The supports of molds in PVA were also printed in PVA through the same nozzle. Using the same material for both molds and supports implied a greater detaching effort.
Manufacturing artifacts and inaccuracies in molds have an impact on the resulting silicone models. In our case, the FDM printer provided an accurate reproduction of the internal and external surfaces of the AAA models (Figures 6 and 7), including a precise replica of the vessel’s curvature at AAA neck and the complex iliac bifurcation in both cases. FDM printers are more affordable than other techniques like laser stereolithography (SLA) printing either for the initial purchase or the used material, printing process costs, and the possibility to print larger volumes. However, SLA would have provided a better surface finish, avoiding the surface roughness typical of the FDM process that required surface treatment. Indeed, surface finish is crucial in determining the outcome of the manufacturing process as surface finish mainly affects the accuracy of measurements or acquisitions in imaging applications.
A careful analysis of casting material was performed to investigate materials able to mimic soft tissues under selected conditions and suitable for imaging applications. The model’s transparency is the most important requirement for PIV, whereas acoustical properties and attenuation coefficient over the diagnostic frequency range (30) are prerequisites of materials in medical ultrasound. Although rigid models could be suitable for flow visualization, the absence of wall compliance could provoke inaccuracies in the hemodynamics and fail in reproducing the actual patient-specific condition. As the result of this analysis, we chose silicone Sylgard 184, a type of silicone rubber commonly used for casting and prototyping applications (49). It is known for its mechanical properties (31, 55), including good tensile strength and elongation. As shown in Figure 10, it is compatible with PIV and other imaging techniques, which can be helpful for pre-surgical planning and research purposes. In addition, Silicone Sylgard 184 has many other useful properties that make it well-suited for use in the fabrication of patient-specific AAA models. These properties include good dimensional stability, low shrinkage, and the ability to replicate fine details. It is also easy to work with and can be easily molded and shaped to create detailed, accurate models with the strength necessary to withstand prolonged physiologically realistic mechanical stresses (27, 55). Additionally, Sylgard 184 can be eventually sealed from exposure to air or water, thus, producing a more stable phantom.
The degassing and casting steps were highly sensitive since any bubble trapped in the material after degassing might harm the outcome. In particular, the presence of bubbles prevents US signal transmission, degrades acquisition quality, and creates scatter artifacts in PIV applications. The use of a controllable vacuum chamber aided in properly degassing the material. It is important to avoid silicone from developing distinguishable layers that would make it difficult to illuminate the various planes for PIV measurements and US acquisitions. Shrinkage can occur during the curing phase, affecting the accuracy of the final model and provoking geometric distortions that affect the model’s flexibility. A major drawback of Sylgard 184 is its tendency to leak out easily from the molds at the junction lines if the molds are not properly sealed. This silicone can also permeate the molds if an adequate number of contours are not set or the contact surfaces are not properly treated. In these regards, the use of the ABS-specific sealant was crucial as was the treatment with PVA glue in the LMM case.
PVA dissolved quickly in warm water (approximately 60/70C) because the printing settings provided for a sufficient trade-off between mold strength and ease of dissolving. Furthermore, the use of an ultrasonic tank sped up the dissolving process, which took three hours for the LMM and up to five hours for the core of the LCM.
The accuracy of the manufactured AAA models, computed in terms of distance to the CAD reference (Figures 11A,B), proved that both casting techniques succeeded in the manufacturing allowing a fine replica of the patient-specific vessels and providing an almost uniform thickness of the final models, with exceptions at inlets and outlets where some wall thickening/thinning occurred especially in the lost core case (Figure 11B).
The presented 3D printing techniques minimized the costs as there was no wastefulness of the casting materials. Considering only the materials (PVA, ABS, SSU00, and Sylgard 184), the two AAA models cost € 38.10 and € 46.78, respectively, to which other costs such as, for example, working hours and energy consumed by the FDM printer should be added. In terms of time, the whole manufacturing process from the image segmentation to the final silicone model can be summarized as a couple of days for segmenting the images and designing the molds, up to three days to have the molds printed, up to two days for the preprocessing (surface treatment and drying time, sealing), half a day for preparing the material, degas and cast it, two days to cure, and one day for the post-processing. Finally, we can estimate the duration of the complete process around 11 days.
In this regard, it is worth noting that, if a laboratory has 3D printers, some equipment, and skilled personnel available, the manufacturing process can be more advantageous than buying a model either in terms of cost or time. Buying the LCM from a specialized 3D printing company would have cost roughly € 800 and taken over a month to receive the model. Furthermore, it would always be possible to reproduce the geometry of a specific patient that is particularly useful to better investigate the anatomy and/or the hemodynamics when the in-vivo analysis is challenging. It is worth noticing that the patients considered in this work did not present intraluminal thrombi (ILT), and no visceral vessels and/or distal vessels were included. However, there are no specific drawbacks to also include these additional features in the phantoms. Regarding the ILT, the related volume can be included in the outer molds of the model, or it can be manufactured through a dedicated casting procedure and then included in the molds before the material casting by following the approach described in Emendi et al. (56). Regarding the inclusion of additional vessels, such as the renal, splenic, and hepatic arteries, the same process can be followed through the definition of appropriate separation surfaces to avoid undercuts and related problems in demolding the phantom. Indeed, being patient-specific geometries, these AAA models faithfully mimic the behavior of the considered vessel, allowing researchers to capture hemodynamic features as in-vivo. In particular, extensive investigations and quantification of hemodynamic features such as velocity, flow patterns, pressures, etc., which play a key role in AAA initiation, development, and potential rupture, can be performed by inserting the AAA models in mock circulatory loop systems (17). Additionally, the AAA models can be used for applications that aim at personalized investigation, such as clinical imaging, surgical training, intervention planning, or therapeutic management that are strictly patient-specific. For instance, a hemodynamic characterization can be performed by processing Color-Doppler images (57) acquired on these models with the noticeable advantage of investigating multiple flow conditions at the inlet, flow splitting at the outlets, and different geometries. However, prospective studies with effective outcomes are beneficial to prove the actual clinical utility of patient-specific AAA models. The accuracy of the manufactured models in representing the patient’s anatomy can be evaluated, possibly by comparing the imaging of the models to the original imaging data. Finally, a comparison between in-vivo and in-vitro hemodynamic measurements should be performed to verify that in-vitro blood flows match in-vivo ones.
5. Conclusion
Two deformable patient-specific AAA models were fabricated using different techniques, both resulting in accurate replicas made of silicone Sylgard 184 with smooth inner and outer surfaces. Our work aimed at describing the entire phantoms’ fabrication workflow in detail and including all the necessary aspects to allow the reproducibility of the process, thus reducing user dependency. The lost molds casting involved printing both the core and outer molds in PVA using an FDM printer and then dissolving the PVA in warm water to create a 1.8 mm thick silicone model. The lost core casting technique involved printing the outer molds in ABS using an FDM printer and the core in PVA, which was then dissolved in warm water to create a 2 mm thick silicone model. At the expense of a less smooth surface, the lost core casting technique was found to be more cost-effective and allowed for the production of multiple replicas of the same model without having to reprint the outer molds. These models can be used in a variety of applications, including biomedical engineering, training, and hemodynamic investigations in MCL systems.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the Danish Scientific Ethics Committees for the Capital Region (De Videnskabsetiske Komiteer for Region Hovedstaden, n° H-18014413). The patients/participants provided their written informed consent to participate in this study.
Author contributions
Conceptualization, MNA, EG, and SC; implementation, MNA, EG, and FB; data resources, AM, AT, and LR; writing – original draft preparation, MNA, EG, and SC; writing – review and editing, MNA, EG, FB, AT, LR, SA, and SC; supervision, LR, SA, and SC. All authors contributed to the article and approved the submitted version.
Funding
This study has received funding from the Marie Sklodowska-Curie grant agreement MeDiTATe project No 859836.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Chaikof EL, Dalman RL, Eskandari MK, Jackson BM, Lee WA, Mansour MA, et al. The society for vascular surgery practice guidelines on the care of patients with an abdominal aortic aneurysm. J Vasc Surg. (2018) 67:2–77. doi: 10.1016/j.jvs.2017.10.044
2. Cyron CJ, Wilson JS, Humphrey J. Mechanobiological stability: a new paradigm to understand the enlargement of aneurysms? J R Soc Interface. (2014) 11:20140680. doi: 10.1098/rsif.2014.0680
3. Bredahl K, Sandholt B, Lönn L, Rouet L, Ardon R, Eiberg J, et al. Three-dimensional ultrasound evaluation of small asymptomatic abdominal aortic aneurysms. Eur J Vasc Endovasc Surg. (2015) 49:289–96. doi: 10.1016/j.ejvs.2014.12.022
4. Yuan Z, Lu Y, Wei J, Wu J, Yang J, Cai Z. Abdominal aortic aneurysm: roles of inflammatory cells. Front Immunol. (2021) 11:609161. doi: 10.3389/fimmu.2020.609161
5. Ghulam QM, Svendsen MBS, Zielinski AH, Eiberg JP. Ultrasound transducer pressure: an unexplored source of abdominal aortic aneurysm measurement error. Ultrasound Med Biol. (2022) 48(9):1778–84. doi: 10.1016/j.ultrasmedbio.2022.05.003
6. Celi S, Cioffi M, Capellini K, Fanni BM, Emanuele G, Vignali E, et al. Advanced non-animal models in biomedical research. JRC Publications Repository (2022).
7. Celi S, Gasparotti E, Capellini K, Vignali E, Fanni BM, Ali LA, et al. 3D printing in modern cardiology. Curr Pharm Des. (2021) 27:1918–30. doi: 10.2174/1381612826666200622132440
8. Vukicevic M, Mosadegh B, Min JK, Little SH. Cardiac 3D printing and its future directions. JACC: Cardiovasc Imaging. (2017) 10:171–84. doi: 10.1016/j.jcmg.2016.12.001
9. Santoro G, Pizzuto A, Rizza A, Cuman M, Federici D, Cantinotti M, et al. Transcatheter treatment of “complex” aortic coarctation guided by printed 3D model. Case Rep. (2021) 3:900–4. doi: 10.1016/j.jaccas.2021.04.036
10. Segaran N, Saini G, Mayer JL, Naidu S, Patel I, Alzubaidi S, et al. Application of 3D printing in preoperative planning. J Clin Med. (2021) 10:917. doi: 10.3390/jcm10050917
11. Gharleghi R, Dessalles CA, Lal R, McCraith S, Sarathy K, Jepson N, et al. 3D printing for cardiovascular applications: From end-to-end processes to emerging developments. Ann Biomed Eng. (2021) 49:1598–618. doi: 10.1007/s10439-021-02784-1
12. Coles-Black J, Barber T, Bolton D, Chuen J. A systematic review of three-dimensional printed template-assisted physician-modified stent grafts for fenestrated endovascular aneurysm repair. J Vasc Surg. (2021) 74:296–306. doi: 10.1016/j.jvs.2020.08.158
13. Pizzuto A, Raimondi F, Celi S, Calabri GB, Spaziani G, Gasparotti E, et al. Transcatheter treatment of native idiopathic multiloculated aortic aneurysm guided by 3D printing technology. JACC: Case Rep. (2022) 8:101662. doi: 10.1016/j.jaccas.2022.09.028
14. Santoro G, Rizza A, Pizzuto A, Berti S, Cuman M, Gasparotti E, et al. Transcatheter treatment of ascending aorta pseudoaneurysm guided by 3D-model technology. Case Rep. (2022) 4:343–7. doi: 10.1016/j.jaccas.2022.01.005
15. Cuman M, Clemente A, Celi S, Santoro G. 3D model-guided transcatheter closure of ascending aorta pseudoaneurysm with the novel Amplatzer Trevisio intravascular delivery system. Catheter Cardiovasc Interv. (2022) 99:140–4. doi: 10.1002/ccd.29937
16. Vignali E, Gasparotti E, Mariotti A, Haxhiademi D, Ait-Ali L, Celi S. High-versatility left ventricle pump, aortic mock circulatory loop development for patient-specific hemodynamic in vitro analysis. ASAIO J. (2022) 68:1272–81. doi: 10.1097/MAT.0000000000001651
17. Bardi F, Gasparotti E, Vignali E, Avril S, Celi S. A hybrid mock circulatory loop for fluid dynamic characterization of 3D anatomical phantoms. IEEE Trans Biomed Eng. (2022) 70(5):1651–61. doi: 10.1109/TBME.2022.3224581
18. Rocchi M, Ingram M, Claus P, D’hooge J, Meyns B, Fresiello L. Use of 3D anatomical models in mock circulatory loops for cardiac medical device testing. Artif Organs. (2022) 47(2):260–72. doi: 10.1111/aor.14433
19. Fanni B, Sauvage E, Celi S, Norman W, Vignali E, Landini L, et al. A proof of concept of a non-invasive image-based material characterization method for enhanced patient-specific computational modeling. Cardiovasc Eng Technol. (2020) 11:532–43. doi: 10.1007/s13239-020-00479-7
20. Zimmermann J, Loecher M, Kolawole FO, Bäumler K, Gifford K, Dual SA, et al. On the impact of vessel wall stiffness on quantitative flow dynamics in a synthetic model of the thoracic aorta. Sci Rep. (2021) 11:1–14. doi: 10.3390/jcm10050917
21. Susanto S, Skote M, Chauhan S. In vitro assessment of combined Doppler ultrasound and CFD modeling in arterial blood flow quantification. Flow Meas Instrum. (2013) 33:218–27. doi: 10.1016/j.flowmeasinst.2013.07.013
22. Asami R, Tanaka T, Shimizu M, Seki Y, Nishiyama T, Sakashita H, et al. Ultrasonic vascular vector flow mapping for 2-D flow estimation. Ultrasound Med Biol. (2019) 45:1663–74. doi: 10.1016/j.ultrasmedbio.2019.02.014
23. Leinan PR, Grønli T, Skjetne P, Wigen MS, Urheim S, Lovstakken L, et al. Comparison of ultrasound vector flow imaging and CFD simulations with PIV measurements of flow in a left ventricular outflow tract phantom. implications for clinical use and in silico studies. Comput Biol Med. (2022) 146:105358. doi: 10.1016/j.compbiomed.2022.105358
24. Saaid H, Voorneveld J, Schinkel C, Westenberg J, Gijsen F, Segers P, et al. Tomographic PIV in a model of the left ventricle: 3D flow past biological and mechanical heart valves. J Biomech. (2019) 90:40–9. doi: 10.1016/j.jbiomech.2019.04.024
25. Bonfanti M, Franzetti G, Homer-Vanniasinkam S, Díaz-Zuccarini V, Balabani S. A combined in vivo, in vitro, in silico approach for patient-specific haemodynamic studies of aortic dissection. Ann Biomed Eng. (2020) 48:2950–64. doi: 10.1007/s10439-020-02603-z
26. Mariotti A, Gasparotti E, Vignali E, Marchese P, Celi S, Salvetti MV. Integrating data in CFD simulations and in experiments of the hemodynamic in healthy and pathologic thoracic aorta. In: International Work-Conference on Bioinformatics and Biomedical Engineering. Springer (2022). p. 208–19.
27. Poepping TL, Nikolov HN, Thorne ML, Holdsworth DW. A thin-walled carotid vessel phantom for Doppler ultrasound flow studies. Ultrasound Med Biol. (2004) 30:1067–78. doi: 10.1016/j.ultrasmedbio.2004.06.003
28. Watts D, Sutcliffe C, Morgan R, Meagher S, Wardlaw J, Connell M, et al. Anatomical flow phantoms of the nonplanar carotid bifurcation, part I: computer-aided design and fabrication. J Ultrasound Biol. (2007) 33:296–302. doi: 10.1016/j.ultrasmedbio.2006.08.003
29. Meagher S, Poepping T, Ramnarine K, Black R, Hoskins P. Anatomical flow phantoms of the nonplanar carotid bifurcation, part II: experimental validation with doppler ultrasound. Ultrasound Med Biol. (2007) 33:303–10. doi: 10.1016/j.ultrasmedbio.2006.08.004
30. Thouvenot A, Poepping T, Peters TM, Chen EC. Characterization of various tissue mimicking materials for medical ultrasound imaging. In: Medical Imaging 2016 Physics of Medical Imaging. Vol. 9783. SPIE (2016). p. 1436–44.
31. Yazdi SG, Geoghegan P, Docherty P, Jermy M, Khanafer A. A review of arterial phantom fabrication methods for flow measurement using PIV techniques. Ann Biomed Eng. (2018) 46:1697–721. doi: 10.1007/s10439-018-2085-8
32. Salsac A-V, Sparks SR, Lasheras JC. Hemodynamic changes occurring during the progressive enlargement of abdominal aortic aneurysms. Ann Vasc Surg. (2004) 18:14–21. doi: 10.1007/s10016-003-0101-3
33. Deplano V, Knapp Y, Bertrand E, Gaillard E. Flow behaviour in an asymmetric compliant experimental model for abdominal aortic aneurysm. J Biomech. (2007) 40:2406–13. doi: 10.1016/j.jbiomech.2006.11.017
34. Wang Y, Joannic D, Patrick J, Keromnes A, Aurélien M, Lalande A, et al. Comparison of flow measurement by 4D flow magnetic resonance imaging and by particles image velocimetry on phantom of abdominal aortic aneurysm. SM Vascular Med. (2016) 1(2):1008.
35. Bauer A, Wegt S, Bopp M, Jakirlic S, Tropea C, Krafft AJ, et al. Comparison of wall shear stress estimates obtained by laser Doppler velocimetry, magnetic resonance imaging and numerical simulations. Exp Fluids. (2019) 60:1–16. doi: 10.1007/s00348-019-2758-6
36. Meess KM, Izzo RL, Dryjski ML, Curl RE, Harris LM, Springer M, et al. 3D printed abdominal aortic aneurysm phantom for image guided surgical planning with a patient specific fenestrated endovascular graft system. In: Medical imaging 2017: imaging informatics for healthcare, research, and applications. Vol. 10138. SPIE (2017). p. 159–72.
37. Kärkkäinen JM, Sandri G, Tenorio ER, Alexander A, Bjellum K, Matsumoto J, et al. Simulation of endovascular aortic repair using 3D printed abdominal aortic aneurysm model and fluid pump. Cardiovasc Intervent Radiol. (2019) 42:1627–34. doi: 10.1007/s00270-019-02257-y
38. Khural M, Gullipalli R, Dubrowski A. Evaluating the use of a generic three-dimensionally (3D) printed abdominal aortic aneurysm model as an adjunct patient education tool. Cureus. (2020) 12(6):e8533. doi: 10.7759/cureus.8533
39. Thirugnanasambandam M, Canchi T, Piskin S, Karmonik C, Kung E, Avril S, et al. Design, development, and temporal evaluation of a magnetic resonance imaging compatible in vitro circulation model using a compliant abdominal aortic aneurysm phantom. J Biomech Eng. (2021) 143:051004. doi: 10.1115/1.4049894
40. Little CD, Mackle EC, Maneas E, Chong D, Nikitichev D, Constantinou J, et al. A patient-specific multi-modality abdominal aortic aneurysm imaging phantom. Int J Comput Assist Radiol Surg. (2022) 17(9):1–7. doi: 10.1007/s11548-022-02612-4
41. Salsac A-V, Sparks S, Chomaz J-M, Lasheras J. Evolution of the wall shear stresses during the progressive enlargement of symmetric abdominal aortic aneurysms. J Fluid Mech. (2006) 560:19–51. doi: 10.1017/S002211200600036X
42. Salman HE, Ramazanli B, Yavuz MM, Yalcin HC. Biomechanical investigation of disturbed hemodynamics-induced tissue degeneration in abdominal aortic aneurysms using computational and experimental techniques. Front Bioeng Biotechnol. (2019) 7:111. doi: 10.3389/fbioe.2019.00111
43. Stamatopoulos C, Mathioulakis D, Papaharilaou Y, Katsamouris A. Experimental unsteady flow study in a patient-specific abdominal aortic aneurysm model. Exp Fluids. (2011) 50:1695–709. doi: 10.1007/s00348-010-1034-6
44. Meyer CA, Bertrand E, Boiron O, Deplano V. Stereoscopically observed deformations of a compliant abdominal aortic aneurysm model. J Biomech Eng. (2011) 133:111004. doi: 10.1115/1.4005416
45. Chen C-Y, Antón R, Hung M-y., Menon P, Finol EA, Pekkan K. Effects of intraluminal thrombus on patient-specific abdominal aortic aneurysm hemodynamics via stereoscopic particle image velocity and computational fluid dynamics modeling. J Biomech Eng. (2014) 136:0310011–19. doi: 10.1115/1.4026160
46. Deplano V, Guivier-Curien C, Bertrand E. 3D analysis of vortical structures in an abdominal aortic aneurysm by stereoscopic PIV. Exp Fluids. (2016) 57:1–11. doi: 10.1007/s00348-016-2263-0
47. Moravia A, Simoëns S, El Hajem M, Bou-Saïd B, Kulisa P, Della-Schiava N, et al. In vitro flow study in a compliant abdominal aorta phantom with a non-newtonian blood-mimicking fluid. J Biomech. (2022) 130:110899. doi: 10.1016/j.jbiomech.2021.110899
48. Kuthe S, Schlothauer A, Bodkhe S, Hulme C, Ermanni P. 3D printed mechanically representative aortic model made of gelatin fiber reinforced silicone composite. Mater Lett. (2022) 320:132396. doi: 10.1016/j.matlet.2022.132396
49. Raffel M, Willert CE, Scarano F, Kähler CJ, Wereley ST, Kompenhans J. Image evaluation methods for PIV. In: Particle image velocimetry. Springer (2018). p. 145–202.
50. Kikinis R, Pieper SD, Vosburgh KG. 3D Slicer: a platform for subject-specific image analysis, visualization, and clinical support. In: Intraoperative imaging and image-guided therapy. Springer (2014). p. 277–89.
51. Celi S, Vignali E, Capellini K, Gasparotti E. On the role and effects of uncertainties in cardiovascular in silico analyses. Front Med Technol. (2021) 3:748908. doi: 10.3389/fmedt.2021.748908
52. Besl PJ, McKay ND. Method for registration of 3-D shapes. In: Sensor fusion IV: control paradigms and data structures. Vol. 1611. SPIE (1992). p. 586–606.
54. Scarpolini MA, Mazzoli M, Celi S. Enabling supra-aortic vessels inclusion in statistical shape models of the aorta: a novel non-rigid registration method. Front Physiol. (2023) 14:1211461. doi: 10.3389/fphys.2023.1211461
55. Vaicekauskaite J, Mazurek P, Vudayagiri S, Skov AL. Mapping the mechanical and electrical properties of commercial silicone elastomer formulations for stretchable transducers. J Mater Chem C. (2020) 8:1273–9. doi: 10.1039/C9TC05072H
56. Emendi M, Støverud KH, Tangen GA, Ulsaker H, Manstad-H F, Di Giovanni P, et al. Prediction of guidewire-induced aortic deformations during EVAR: a finite element and in vitro study. Front Physiol. (2023) 14:1098867. doi: 10.3389/fphys.2023.1098867
Keywords: abdominal aortic aneurysm, 3D printing, fused deposition modeling, lost molds casting, lost core casting
Citation: Antonuccio MN, Gasparotti E, Bardi F, Monteleone A, This A, Rouet L, Avril S and Celi S (2023) Fabrication of deformable patient-specific AAA models by material casting techniques. Front. Cardiovasc. Med. 10:1141623. doi: 10.3389/fcvm.2023.1141623
Received: 10 January 2023; Accepted: 24 August 2023;
Published: 11 September 2023.
Edited by:
Junjie Xiao, Shanghai University, ChinaReviewed by:
Ying He, Dalian University of Technology, ChinaPietro Cefali’, ASST Papa Giovanni XXIII - Vascular Surgery Unit, Italy
© 2023 Antonuccio, Gasparotti, Bardi, Monteleone, This, Rouet, Avril and Celi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Simona Celi cy5jZWxpQGZ0Z20uaXQ=