Abstract
Purpose:
To compare the flow patterns and hemodynamics of the AFX stent-graft and the covered endovascular reconstruction of aortic bifurcation (CERAB) configuration using laser particle image velocimetry (PIV) experiments.
Materials and Methods:
Two anatomically realistic aortoiliac phantoms were constructed using polydimethylsiloxane polymer. An AFX stent-graft with a transparent cover made with a new method was inserted into one phantom. A CERAB configuration using Atrium’s Avanta V12 with transparent covers made with a previously established method was inserted into the other phantom, both modified stent-grafts were suitable for laser PIV, enabling visualization of the flow fields and quantification of time average wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT).
Results:
Disturbed flow was observed at the bifurcation region of the AFX, especially at the end systolic velocity (ESV) time-point where recirculation was noticeable due to vortical flow. In contrast, predominantly unidirectional flow was observed at the CERAB bifurcation. These observations were confirmed by the quantified hemodynamic results from PIV analysis where mean TAWSS of 0.078 Pa (range: 0.009–0.242 Pa) was significantly lower in AFX as compared with 0.229 Pa (range: 0.013–0.906 Pa) for CERAB (p<0.001). Mean OSI of 0.318 (range: 0.123–0.496) in AFX was significantly higher than 0.252 (range: 0.055–0.472) in CERAB (p<0.001). Likewise, mean RRT of 180 Pa−1 (range: 9–3603 Pa−1) in AFX was also significantly higher than 88 Pa−1 (range: 2–840 Pa−1) in CERAB (p=0.0086).
Conclusion:
In this in vitro study, the flow pattern of a modified AFX stent-graft was found to be more disturbed especially at the end systolic phase, its hemodynamic outcomes less desirable than CERAB configuration.
Clinical Relevance:
While the AFX stent-graft has an advantage over the CERAB configuration in eliminating radial mismatch, and maintaining the anatomical bifurcation for future endovascular intervention, this in vitro study revealed that the associated lower TAWSS, higher OSI and RRT may predispose to thrombosis and are, thus, less desirable as compared to a CERAB configuration. Further investigation is warranted to confirm whether these findings translate into the clinical setting.
Keywords
Introduction
The covered endovascular reconstruction of the aortic bifurcation (CERAB) technique is a relatively new approach in treating extensive aortoiliac occlusive disease.1,2 One of the advantages of the CERAB configuration as compared to the traditional kissing stent procedure is that it reduces radial mismatch, which is defined as the discrepancy between the stented lumen and the vessel lumen after stent placement. 3 This leads to more favorable flow conditions. Initial clinical data for CERAB are encouraging, with 87.3% and 82.3% primary patency rate at 1- and 2-year follow-up in a group of 103 patients with 88% TransAtlantic Inter-Society Consensus II D lesions 3 as compared with 91.5% primary patency rate at 2-year follow-up for open surgery. 4 Taeymans et al 5 further updated the CERAB results with a total of 130 patients showing that a primary patency of 82.0% was achieved at 3-year follow-up as compared with 90.2% for open surgery. 6
More recently, it was suggested that the use of a unibody stent-graft, designed for aneurysm treatment, could have several advantages for the treatment of aortoiliac occlusive disease over CERAB. It preserves the aortic bifurcation and thus allows for future cross-over endovascular interventions and avoids limb competition in the distal aorta.7 –10 Like CERAB, it also protects against potentially fatal aortoiliac rupture during dilatation of heavily calcified lesions. The AFX unibody stent-graft (Endologix, Inc, Irvine, CA, USA) is a thin-walled expanded polytetrafluoroethylene (ePTFE) sheath with an endoskeleton made of Elgiloy (cobalt-chromium alloy). The benefit of the AFX stent-graft lies in the fact that radial mismatch is completely eliminated because of the one-piece bifurcation design. 7 The safety, efficacy, and early patency rates of the AFX stent-graft for treatment of aortoiliac occlusive disease were examined by several groups, totaling 128 patients.7 –10 Primary patency rate for these studies were reported in the range of 80% to 100% at 1-year follow-up. The primary and secondary patency for the largest study with 90 patients were 78.8% and 100% at 3-year follow-up.
Even though both the off-label use of Advanta V12 in CERAB and AFX endograft are effective in treating extensive aortoiliac occlusive disease, the hemodynamic differences between them are currently unknown. As such, an in vitro study would provide new and valuable information for vascular surgeons and existing stent manufacturers to make further improvement in the treatment of aortoiliac occlusive disease. Particle image velocimetry (PIV) is a powerful experimental method to analyze hemodynamic parameters such as flow pattern, wall shear stress (WSS), oscillatory flow index (OSI) and relative residence time (RRT) in an in vitro setting. Laser PIV is the current gold standard for flow imaging,11 –13 its use will help to better understand the influence of hemodynamic factors on blood flow in general and their impact on atherosclerosis in particular. A new method in making transparent cover for AFX stent-graft was introduced for the laser PIV study. The aim of this study is to test the hypothesis that the AFX stent-graft is superior to CERAB configuration where undisturbed flow patterns, physiological WSS 0.5 to 1.2 Pa, lower OSI and RRT are considered as better hemodynamic outcomes. 14 Figure 1 illustrates the 2 approaches with Figure 1A and D showing the anterior view of AFX and CERAB, respectively.

Computed tomography (CT) images of the flow phantoms. (A) Anterior view of AFX. (B) Lateral view of AFX. (C) Superior view of AFX. (D) Anterior view of covered endovascular reconstruction of aortic bifurcation (CERAB). (E) Lateral view of CERAB. (F) Superior view of CERAB. The red lines indicate the location of the lase sheet.
Materials and Methods
Model Designs and Stent Deployment
Two identical aortoiliac bifurcation anatomies in which physiologic flow was reproduced, were used in this experimental study. A CERAB configuration with modified covered stents was deployed in one flow model and a modified AFX unibody stent in the other. The detailed production process for making these flow models has been previously described. 15 The CERAB model was configured using 3 modified covered stents comprising two 8-mm and one 12-mm balloon expandable stents. These stents were derived from Advanta V12 stents (Getinge Group, Wayne, MA, USA) where the ePTFE graft covers were manually removed and subsequently individually coated with a transparent polyurethane cover (Tecoflex TPU clear; The Lubrizol Corporation, Wickliffe, OH, USA) using a dip coating process, which facilitated the visualization of flow patterns inside the flow model.3,16 The polyurethane dipping process was not suitable to replace the ePTFE cover of the AFX stent because of the large spaces between the struts in its expanded configuration. The AFX stent cage was supplied by Endologix with 22-mm main body and 13-mm limbs. The cover with 18-mm body and 10-mm limb was 1 mm oversized compared with the wall-less, rigid polydimethylsiloxane (PDMS) phantom with 16-mm main body and 8-mm limb. The AFX ePTFE cover was replaced by low-density polyethylene (LDPE) film in the shape of 2 cylindrical-shaped grafts, which were stitched together to form the unibody. The covering material was on the outer layer, as with the original stent. Details of the AFX transparent cover-making process are found in Appendix A. A transparent LDPE cover was stitched onto the proximal frame the same way as the original, opaque ePTFE cover was secured on the frame of the AFX endograft. By replicating this feature, it enables the cover to move independently of the frame and we expected it to have the same sealing effect as the original AFX and allows for expansion throughout the cardiac cycle.
Experimental Flow Circuit
An experimental flow circuit as shown in Figure 2 was used to study and compare the hemodynamics in the AFX and CERAB. The flow loop was driven by a programmable hydraulic piston-driven pump (SuperPump) and its accompanying software Vivitest (Vivitro Labs Inc, Victoria, CA). A blood-mimicking fluid comprising water (44.07%), glycerol (34.52%), and urea (21.41%) was used to obtain a fluid with a refractive index optically matched to that of the PDMS model and physiologically matched to that of human blood (dynamic viscosity 4.2 mPa·s, density 1114 kg/m3). 17 The blood-mimicking fluid was circulated through the flow loop by the SuperPump and a suprarenal flow profile mimicking human resting conditions with a mean flow rate of 1.6 L/min (or a stroke volume of 26.7 mL) and 60 beats per minute heart rate was generated by the Vivitest. A CORRI_FLOW M55 Coriolis flow meter (Bronkhorst High-Tech BV, Ruurlo, the Netherlands) was used to provide a real time display of the targeted suprarenal flow profile. Three UF8B ultrasound flow sensors (Cynergy3 Components Ltd, Wimborne, Dorset, UK) and 3 needle valves article number 136230 (Riegler Medical, Nieder-Ramstadt, Germany) were connected to left renal outlet, right renal outlet and the combined iliac outlets, respectively. A display unit provided real-time display of the flow rates of blood-mimicking fluid flowing through the ultrasound flow sensors. The needle valves were manually adjusted to achieve a mean pressure of 100 mm Hg and approximately 25% outflow, or 0.4 L/min, at the left and right renal outlet and approximately 50%, or 0.8 L/min, outflow at the combined iliac outlet. In addition, the 2 iliac outlets were connected to a single compliance chamber. A DS44 handheld sphygmomanometer (Welch Allyn, Skaneateles Falls, NY, USA) was used to set the distal peripheral pressure between 80 and 120 mm Hg.

Experimental flow circuit and laser particle image velocimetry (PIV) setup.
Flow Visualization
In order to facilitate the visualization of the flow pattern inside the phantom, a continuous laser sheet was projected horizontally into the center of the stented lumen region of interest (ROI). The position of the phantom was manually adjusted in such a way that the laser sheet was projected onto the center of inflow and bifurcation ROIs as shown in Figure 3. Even though the ROIs for the AFX and CERAB are not identical due to minor differences in their configurations however, the scientific method of determining the velocity vector and wall shear stress using PIV method is the same in both cases. Therefore, it is not expected to have any direct influence on comparing the results. The laser sheet was slightly obscured by the renal branch at the inflow proximal section of AFX therefore its inflow ROI was shifted distally slightly as compared with that of the CERAB inflow. In addition, the bifurcation ROI for the AFX stent-graft was close to matching the native bifurcation because the configuration of the unibody AFX stent-graft was designed to mimic just that. However, the bifurcation ROI for the CERAB was located about 20 mm upstream. Rhodamine-coated fluorescent polymethyl methacrylate particles (size, 1–20 μm; density 1190 kg/m3) were suspended in the blood mimicking fluid. Digitized frames were collected by illuminating with a continuous-wave laser (5 W DPSS laser, 532 nm; Cohlibri, Lightline, Germany) and by recording using a high-speed camera (FASTCAM SA-Z; Photron Inc, West Wycombe, Buckinghamshire, UK). The images with 1024×1024 pixels resolution and 8 bit/pixel grayscale intensity were captured at a frequency of 2000 Hz. The camera was mounted with its optical axis perpendicular to the laser sheet. Two sets of images were captured for the CERAB-stented phantom at the proximal inflow and aortoiliac bifurcation and then repeated for the AFX-stented phantom.

Images of inflow region illuminated by laser during particle image velocimetry (PIV) experiment. (A) Covered endovascular reconstruction of aortic bifurcation (CERAB) configuration. (B) AFX stent-graft. The red dots are the origins and the redlines are the segment along vessel wall where wall shear data were analyzed. (a) CERAB inflow, (b) CERAB bifurcation right edge, (c) CERAB bifurcation left edge, (d) AFX inflow, (e) AFX bifurcation right edge, (f) AFX bifurcation left edge. x, position in vertical axis; y, position in horizontal axis; u, velocity in vertical axis; v, velocity in horizontal axis.
Data Processing
Data processing was described previously. 18 Since there are 2 regions of interest per stent configuration, data processing was performed on a total of 4 data sets. For each data set, 10 cardiac cycles were obtained and averaged to provide one averaged cycle per data set. Using an in-house developed MATLAB (version 2016b; The MathWorks Inc, Natick, MA, USA) script, the following hemodynamic parameters were calculated along the vessel wall: inflow, bifurcation right edge and bifurcation left edge were as shown in Figure 3. TAWSS, OSI, and RRT were computed according to the equations commonly used in hemodynamic studies. 19
TAWSS is computed as the average WSS over a cardiac cycle where only the magnitude of the shear stress is taken into consideration. It is the average shear stress axially exerted on the vessel wall due to the blood flow motion. A value of 0.5 to 1.2 Pa is considered as physiological WSS in human arterial blood flow and values below and above this range are considered low and high WSS, respectively.20,21 On the other hand, OSI takes into account the WSS direction with respect to the mean blood flow direction during the cardiac cycle. The OSI can vary from 0 to 0.5, 0 meaning no oscillation and 0.5 implying maximum oscillation. OSI has been studied as a relevant biomarker, for the assessment of vessel areas that are predisposed to develop atherosclerotic plaque.22,23 RRT indicates the residence time of fluid particles spent at the vessel wall. Wall region with low WSS and elevated RRT may be prone to thrombus deposition.24,25
Figure 4 shows the difference between flow profiles measured by flow sensor as compared with that determined with PIV method at the inflow over twenty cardiac cycles. The variation of flow rate is one magnitude larger using flow sensor measurement as compared with PIV method based on the statistical analysis of standard deviations.

Flow waveform: blue line represents flow waveform obtained with flow sensor at the inlet of the system, red line represents the flow profile obtained by particle image velocimetry (PIV) analysis in the inflow region of the covered endovascular reconstruction of aortic bifurcation (CERAB) configuration inside the phantom. Both flow profiles are averaged over ten cardiac cycles; hence, shadow bars represent 2 standard deviations (SD) over 10 cycles. The green, orange, and purple dots indicate peak systolic velocity (PSV), end systolic velocity (ESV), and peak diastolic velocity (PDV) time-points, respectively.
Results
Inflow Pattern
Figure 5 shows the inflow patterns for AFX and CERAB at PSV, ESV and PDV time-points. At the PSV time-point, the maximum velocity was achieved as the blood-mimicking fluid accelerated into the phantom. A disturbed flow pattern was observed at the AFX inflow especially at the ESV time-point as shown in Figure 5B. Flow reversal was observed during the end-systolic phase and recirculation was noticeable in the AFX phantom (Video 1, online only; B). In contrast, unidirectional flow was observed in the CERAB inflow without noticeable disturbances (Video 2, online only; E).

Inflow patterns: (A) AFX at peak systolic velocity (PSV) time-point, (B) AFX at end systolic velocity (ESV) time-point, (C) AFX at peak diastolic velocity (PDV) time-point, (D) covered endovascular reconstruction of aortic bifurcation (CERAB) at PSV time-point, (E) CERAB at ESV time-point, and (F) CERAB at PDV time-point.
Bifurcation Flow Pattern
Bifurcation flow pattern for AFX and CERAB at PSV, ESV and PDV time-points are shown in Figure 6. Disturbed flow was observed in the AFX stent in all 3 phases of the triphasic flow especially more pronounced during the end-systolic phase (Video 3, online only; A-6C). In contrast, only slightly disturbed flow was observed in the CERAB configuration during the end-systolic phase (Video 4, online only; E).

Bifurcation flow pattern: (A) AFX at peak systolic velocity (PSV) time-point, (B) AFX at end systolic velocity (ESV) time-point, (C) AFX at peak diastolic velocity (PDV) time-point, (D) covered endovascular reconstruction of aortic bifurcation (CERAB) at PSV time-point, (E) CERAB at ESV time-point and (F) CERAB at PDV time-point.
Hemodynamic Comparison
Table 1 summarizes the hemodynamic comparison (TAWSS, OSI, and RRT) of AFX and CERAB at the inflow section and the bifurcation. The mean and range for each set of TAWSS, OSI and RRT data were determined and presented in the table. In addition, the p-values for each pairwise Student’s t test between AFX stent-graft and CERAB configuration are also included in the table. The quantified hemodynamic parameters are in agreement with the observation of more disturbed flow in the bifurcation of the AFX stent-graft as compared with CERAB configuration. More details of hemodynamic plots of TAWSS, OSI, and RRT along the vessel wall can be found in Appendix B.
Time-Averaged Wall Shear Stress (TAWSS), Oscillatory Shear Index (OSI), and Relative Residence Time (RRT) Results at Inflow and Bifurcation for AFX Stent-Graft and CERAB Configuration a .
Abbreviation: CERAB, covered endovascular reconstruction of aortic bifurcation.
“n” is the number of spatial points available at the regions of interest.
Discussion
In the current study, a method in making transparent cover for AFX stent-graft was introduced. Laser PIV experiments were performed in AFX and CERAB flow phantoms. The images obtained from the experiments were processed to analyze flow patterns and compare TAWSS, OSI, and RRT. A more disturbed and oscillatory flow was observed in the flow pattern of the bifurcation region in the AFX stent-graft as compared with the CERAB covered stent configuration that showed an undisturbed flow pattern. This result was supported by the quantified hemodynamic parameters where TAWSS was lower in AFX and RRT was higher in AFX stent-graft as compared with the CERAB configuration. This was contrary to the hypothesis of our study. The AFX stent-graft designed for treatment of aneurysmal disease, but used by some surgeons for treatment of occlusive disease, has the advantage over a CERAB configuration because radial mismatch between the stented lumen and vessel lumen is eliminated due to its unibody design. Moreover, the aortic bifurcation is preserved permitting repeated long-term endovascular access. However, from flow pattern and hemodynamic perspective, this study showed that CERAB configuration potentially has the advantage of preventing thrombosis due to less disturbed flow, higher TAWSS, lower OSI and RRT as compared with AFX. WSS below physiological level has been shown to be predictive of plaque development based on the hypothesis that low shear stress promote modification of mass transport of atherogenic substances between the lumen and the vessel wall to cause atherosclerosis. Low and oscillating WSS is associated with the development of thrombosis. 14
TAWSS and OSI in this study ranged from 0.009 to 0.906 Pa and 0.1 to 0.496, respectively as compared with mean inflow WSS at PSV, ESV, and PDV time-points of 1 Pa, −0.3 Pa, and 0.4 Pa; and OSI ranged from 01 to 0.2 at the bifurcation in the previous similar hemodynamic study, 3 which is in line with this current study. In addition, WSS values ranging from 0 to 0.80 Pa and OSI values ranging from 0.01 to 0.49 were obtained in the CERAB results involving WSS study, 26 also in agreement with our current study. There was a slightly more disturbed flow at ESV at the CERAB bifurcation in the current study compared with what was observed in the control model of the previous hemodynamic comparison study. 3 The higher WSS values in the CERAB configuration in the bifurcation could also be attributed to the fact that the flow transition from the bifurcation in the distal aorta to the iliac branches is more orderly with higher velocity near the boundary layer along the vessel wall, which might provide protection from thrombosis. It is important to note that the TAWSS observed in both the AFX stent-graft and the CERAB configuration at the inflow and at the bifurcation in this in vitro study were lower than the physiological WSS value of 0.5 to 1.2 Pa in human arteries.20,22,23 In addition, the maximum value of TAWSS for CERAB at inflow was 0.906 Pa at the bifurcation which was significantly higher than 0.242 Pa for AFX therefore more protective of thrombosis for CERAB in the bifurcation region (p=5.53×10−6). Since blood is not in direct contact with endothelial cells, due to the ePTFE graft material used in the covers of both AFX stent-graft and Atrium V12 stents, the graft material provides some protection even though TAWSS might be below physiological values of 0.5 Pa. 20 However, platelet deposition and neointimal hyperplasia is a risk factor for synthetic graft. Furthermore, graft thrombosis is possible even with use of cover. 22
There are several limitations to our study: First, this in vitro study used polyurethane and LDPE for the modified transparent covers with coefficient of friction (0.12–0.70) and (0.20–2.5), respectively, which are slightly higher than that of the original opaque ePTFE covers (0.05–0.10) therefore might influence the WSS. Second, the PIV method is a two-dimensional analysis of a 3-dimensional flow phenomenon therefore only one slice of the 2-dimensional flow pattern centered on the lumen could be recorded for each PIV analysis. Multiple slices would be required to get a more complete picture of 3-dimensional flow across the lumen as a whole, which could reveal areas of stagnant flow with low TAWSS in the helical CERAB configuration. Third, the laser sheet was obscured by the renal branch of the AFX phantom therefore unable to capture PIV images upstream of the AFX stent-graft. In the existing in vitro setup, a suprarenal flow profile is prescribed at the inlet and a triphasic flow profile at the infrarenal aorta is achieved by adjusting resistances of renal and iliac arteries in addition to the adjustment of the compliance of the iliac arteries. The renal arteries can be eliminated if the pump can provide an infrarenal flow profile. An alternative approach would be to apply an infrarenal flow profile at the inlet directly using a novel roller pump, 27 or a similar pump, in order to eliminate the renal branches in the flow model as to allow an unobstructed laser sheet entry into the model to produce complete AFX inflow. Finally, the proximity of the AFX to the renal arteries could explain part of the recirculations in the inflow, 26 which could generate strong flow disturbances that are less normalized at the AFX relative to the CERAB. However, clinically AFX is likely to be deployed nearer to the renal arteries than CERAB, which remains low.
The aortoiliac model used is considered rigid, which in many cases can be considered appropriate given the often highly calcified lesions in patients with aortoiliac occlusive disease. The AFX graft and CERAB configuration have different properties in terms of flexibility, which may also influence patency results. However, this was outside the scope of the current study. The disturbed and oscillatory flow might have been caused by the inner skeleton of the AFX stent-graft, protruding into the vessel lumen. In addition, the infolding of the graft material could also be a factor as the stent-graft was slightly oversized. We followed the manufacturer recommendation of 10% to 20% oversize where the stent graft was 1 mm oversized as compared with the lumen diameter. In the current study a transparent AFX device was designed, in order to enable laser PIV, considered the gold standard of flow imaging. Future studies using the original AFX, in combination with ultrasound-PIV might confirm whether infolding of LDPE graft material contributed to the disturbed flow. Another future study using computational fluid dynamics (CFD) methods to simulate the experimental setup can be performed to predict the hemodynamics of the 3-dimensional flow model. When the results of this study are combined with previous literature, 3 the hemodynamic consequences of bare metal stent, covered kissing stents, CERAB configuration and the AFX stent-graft for treating aortoiliac occlusive disease can be compared. From a hemodynamic perspective, CERAB seems to have a slightly better performance compared with the other 3 types of stents. The AFX stent-graft has the geometric advantage of eliminating radial mismatch but this is counteracted by disturbed flow in the bifurcation region. Therefore, it might be beneficial to investigate design improvement such as an inner ePTFE graft to smoothen the protruding strut links. This approach potentially reduces disturbed flow if it is to be considered a definitive treatment for occlusive disease rather than aneurysm for which it was designed. In the present study, we were only interested in comparing the hemodynamics of the 2 stent configurations. As such a negative control does not add value and we did not include a negative control without a stent graft.
Conclusion
In this in vitro study, it was found that the AFX may be less desirable as compared to the CERAB configuration with regards to flow changes, and when used for treatment of aortoiliac occlusive disease. Further investigations are warranted to confirm whether the in vitro hemodynamic results identified in this study are translated into the clinical setting and correlate with an increased incidence of thrombosis.
Footnotes
Appendix A
Fashioning Process for AFX Transparent Cover.
| Step | Description | Process |
|---|---|---|
| 1 | Remove ePTFE cover from the original AFX stent-graft. | |
| 2 | Cut LDPE film to the desired shape. | |
| 3 | Heat-weld the LDPE films to form 2 cylindrical shaped tubes with 1 mm oversized compared with model, ie, 17-mm diameter cover for the 16-mm aortic diameter in the model and 9-mm diameter cover for the 8-mm iliac diameter in the model. | |
| 4 | Suture the 2 parts to form the LDPE cover | |
| 5 | Insert the collapsed stents into silicone tubes. | |
| 6 | Insert collapsed stent into LDPE cover | |
| 7 | Redeploy self-expanding stent suture ends of cover to stents after trimming. | |
| 8 | Insert collapsed LDPE covered stent into silicone tubes. | |
| 9 | Insert collapsed stent with LDPE cover into flow phantom | |
| 10 | Redeploy the neo AFX self-expanding stent with LDPE cover in the flow phantom |
Abbreviations: ePTFE, expanded polytetrafluoroethylene; LDPE, low-density polyethylene.
Appendix B
Author Contributions
Conception and design: Chong, Groot Jebbink, Jansen, and Sun. Analysis and interpretations: Chong, Groot Jebbink, and Mirgolbabaee. Data collection: Chong, Groot Jebbink, and Mirgolbabaee. Writing the article: Chong and Groot Jebbink. Critical revision of the article: Chong, Doyle, Groot Jebbink, Mirgolbabaee, Jansen, Reijnen, Sun, van de Velde, and Versluis. Statistical analysis: Chong. Overall responsibility: Chong, Doyle, Jansen, and Sun
Declaration of Conflicting Interests
The author(s) declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: Michel Reijnen is a consultant with Endologix and Bentley Innomed.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Curtin University (Grant/Award Number: FieldTrip ST3483).
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
