ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 18 June 2025

Sec. Biomechanics

Volume 13 - 2025 | https://doi.org/10.3389/fbioe.2025.1585345

This article is part of the Research TopicDiagnostic and Predictive Roles of Computational Cardiovascular Hemodynamics in the Management of Cardiovascular DiseasesView all 15 articles

Patient-specific modelling of pulmonary arterial hypertension: wall shear stress correlates with disease severity

  • 1National Heart and Lung Institute, Imperial College London, London, United Kingdom
  • 2Department of Chemical Engineering, Imperial College London, London, United Kingdom
  • 3National Pulmonary Hypertension Service, Imperial College Healthcare NHS Trust, London, United Kingdom
  • 4Medical Research Council, Laboratory of Medical Sciences, Imperial College London, London, United Kingdom

Introduction: Pulmonary arterial hypertension (PAH) requires an invasive right heart catheter (RHC) procedure for diagnosis. Patients can present with initial symptoms and interact with healthcare institutes for up to 3 years before referral for diagnosis. Thus, there is a great need to develop non-invasive tools, to better screen patients and improve early diagnosis rates.

Methods: seven patients diagnosed and treated for PAH were included in this study. Patient-specific computational fluid dynamics (CFD) models were built for all patients, with all model parameters tuned using non-invasive imaging data, including CT, cardiac MR, echocardiogram, and 4D-flow MRI scans–crucially, a 3D inlet velocity profile was derived from 4D-flow MRI.

Results: CFD models were quantitatively and qualitatively well matched with in-vivo 4D-flow hemodynamics. A linear correlation of R2 = 0.84 was found between CFD derived time-averaged wall shear stress (TAWSS) and RHC measured mean pulmonary pressure (key diagnostic value): low TAWSS correlated with high pressure.

Conclusion: This study highlights TAWSS as a potential computational biomarker for PAH. The clinical use of TAWSS to diagnose and stratify PAH patients has the potential to greatly improve patient outcomes. Further work is ongoing to validate these findings in larger cohorts.

1 Introduction

Pulmonary arterial hypertension (PAH) is characterized by high blood pressure (>20 mmHg) in the pulmonary arteries (Simonneau et al., 2019). Current medical guidelines require an invasive right heart catheter (RHC) procedure to measure pulmonary arterial pressure (PAP) and diagnose PAH (Humbert et al., 2023). In the UK and Ireland, an RHC will be conducted at one of nine specialist pulmonary hypertension centres. However, studies have shown that PAH patients can present to local healthcare institutions up to 3 years before diagnosis with initial symptoms such as breathlessness (Kiely et al., 2019a). In the period between first presentation and referral for RHC these patients will continue to interact with the healthcare system through A&E visits, inpatient hospitalisations, outpatient consultations, etc., (Lawrie et al., 2024; Bergemann et al., 2018). Within this time, extensive data, including blood tests, clinical tests, and imaging studies, will be collected (Kiely et al., 2019a; 2019b). Despite a wealth of information available on these patients, and increasing awareness of PAH, there has been little improvement in time from symptom onset to diagnosis. Early diagnosis of PAH is known to improve PAH management and reduce healthcare costs (Lawrie et al., 2024). Thus, there is a great need to develop non-invasive methods to estimate PAP that can utilise the extensive data sets available on patients presenting continuously, to better identify those that are at risk of PAH.

Transthoracic echocardiogram (TTE) is one such non-invasive method from which PAP can be estimated. While TTE is cheap and can easily be performed bedside, several studies have reported its high potential for misdiagnosis - in particular, TTE performs poorly on patients with PAPs close to the diagnostic threshold of 20 mmHg (Testani et al., 2010; Fisher et al., 2009; Rich et al., 2011). 4D-flow MRI is a relatively newer non-invasive technology that captures 3D velocity fields in the scan region of interest at multiple time points within a cardiac cycle. This provides an extensive amount of highly valuable in-vivo hemodynamic data. Studies have identified 4D-flow MRI measured parameters including time-averaged wall shear stress (TAWSS) and vortical flow to correlate with disease severity (Schafer et al., 2016; Cerne et al., 2022; Barker et al., 2015; Reiter et al., 2020; Deux et al., 2022). However, exact magnitudes of velocity derivative parameters are known to be inaccurate when evaluated from 4D-flow directly due to the low spatial and temporal resolution of the scans (Miyazaki et al., 2017; Zhang et al., 2022), with spatial resolution typically in the region of 2.5 mm3 (Zilber et al., 2021), and temporal resolution approximately 0.04 s (equivalent to 25 time points in a cardiac cycle of length 1 s).

One method which can overcome the temporal and spatial resolution limitations of 4D-flow is computational fluid dynamics (CFD). CFD has been shown to provide highly detailed hemodynamic information in numerous cardiovascular diseases, with the ability to accurately characterise complex hemodynamic patterns beyond velocity alone (Bordones et al., 2018; Fanni et al., 2023; Sabry et al., 2023; Tsubata et al., 2019; 2023). Within PAH, CFD studies have indicated that TAWSS is lower in PAH patients compared to healthy volunteers (Tang et al., 2012; Pillalamarri et al., 2021). Abnormal TAWSS magnitudes and patterns are known to influence endothelial cell function in the lining of vessel walls (Chiu et al., 2004; 2005; 2007; Hastings et al., 2007; Wang et al., 2013; Arshad et al., 2021; Armour et al., 2024). With the knowledge that PAH leads to arterial stiffening and remodelling, there is a desire and opportunity to better understand how PAH driven hemodynamics contribute to disease progression. Additionally, understanding the connection between TAWSS and PAH may allow for improved diagnosis and treatment development.

While CFD provides more confidence in exact magnitudes of hemodynamics parameters, prior studies that have analysed pulmonary hemodynamics, including those that have specifically identified TAWSS patterns, in PAH patients and the interpretation of their results are limited by the choice of model set up parameters, especially boundary conditions applied with idealised flat or parabolic profiles (Tang et al., 2012; Bordones et al., 2018; Pillalamarri et al., 2021; Fanni et al., 2023). Sabry et al. (2023) significantly improved on these prior studies by implementing 2D inlet velocity profiles, in six subjects, derived from patient-specific 2D-flow MRI scans. Within the region of interest, the main pulmonary artery, being adjacent to the ejection of blood from the right ventricle through the pulmonary valve, flow patterns in the main pulmonary artery are highly non-uniform and individual in each patient. The application of a simplified flat or parabolic inlet velocity profile (as has often been applied in previous CFD models of PAH) will likely be to miss distinct flow features in the pulmonary arteries. While 2D profiles are a step closer towards physiological accuracy, studies in the ascending aorta (anatomically analogous to the main pulmonary artery) have shown that only a patient-specific 3D inlet velocity profile (derived from 4D-flow MRI data) can accurately model hemodynamics (Pirola et al., 2017; Armour et al., 2020). Additionally, Armour et al. (2020) demonstrated the influence of stroke volume on TAWSS, with a 25% reduction in stroke volume leading to a 35% reduction in TAWSS–an important finding to keep in mind when interpreting CFD results based on a population average pulmonary flow profile.

The use of patient-specific 4D-flow MRI based CFD models to derived high resolution hemodynamic information has been demonstrated and validated in various diseased and healthy arterial models (Armour et al., 2022; Yin et al., 2024; Wang et al., 2024). Thus, the aim of this study is to apply a 4D-flow MRI based CFD methodology to a PAH cohort, to accurately evaluate TAWSS and assess any correlation to RHC measured mean PAP. To the best of our knowledge, this is the first study to apply a fully patient-specific hemodynamic workflow by implementing a 3D inlet velocity profile.

2 Materials and methods

Seven patients (P1-7) diagnosed and treated for PAH at Hammersmith Hospital, London were included in this study. Collection and use of all data had approval from the UK National Health Service Health Research Authority. All patients were undergoing regular checkups and monitoring (including, RHC, echocardiogram, CT) as standard of care, with clinical information being collected as part of the TRIPHIC database (IRAS ID: 210461). Additionally, the patients were recruited to undergo a 4D-flow MRI scan alongside their standard scans (IRAS ID: 280472). The resolution of the 4D-flow scans for all patients were in the range of 2 mm × 2.375–2.625 mm × 2.375–2.625 mm. Among the seven patients, the time between the RHC and 4D-flow MRI scan varied between one and 30 days (average of 9 days).

3D geometric models of the pulmonary artery, main bifurcation, and right (RPA) and left pulmonary artery (LPA) branches were reconstructed from the 4D-flow MRI data. Using Mimics (Materialize v22), 2D masks of the pulmonary artery were segmented on each slice of all three velocity phase images as well as the magnitude images. Individually, this would allow for four separate 3D models to be generated, with different sections of the pulmonary artery segmented depending on the dominant direction of flow in each region. Therefore, by combing all layers of the four masks through the inbuilt addition ‘Boolean Operations’, a final complete 3D geometry of the whole pulmonary artery was reconstructed–the geometric models for all seven patients can be seen in Figure 1. Geometries were meshed in ICEM (Ansys, v22), containing a hexahedral core and 10 prism layers at the wall. Mesh sensitivity tests were performed on one exemplar geometry to ensure mesh independent solutions–mesh parameters were deemed acceptable when velocity and wall shear converged below 5%. The chosen mesh parameters were then set for all subsequent geometries to ensure consistency. Final meshes varied between 0.46–1.09 million elements.

Figure 1
www.frontiersin.org

Figure 1. Geometric models of seven PAH patients (P1-7) reconstructed from 4D-flow MRI data.

Using in-house python codes (Saitta, 2022; Saitta et al., 2023), 3D inlet velocity profiles were extracted from each patient’s 4D-flow MRI scan and applied as the inlet boundary condition. 3-element Windkessel (3EWK) models were tuned for the right and left pulmonary artery in each model following the methodology reported by Pirola et al. (2017). To make the 3EWK parameters patient-specific, a total branch resistance must be estimated from average pressure and flow values for each patient. Average flow to the left and right pulmonary artery was measured from 4D-flow MRI data, while average pressure for each patient was extracted from echocardiogram. Peak pressure was first estimated using the simplified Bernoulli’s equation: peak PAP [mmHg] = 4·v2 + RAP, where v is the peak tricuspid valve velocity measured during the echocardiogram, and RAP is an estimate of the right atrial pressure, set by the clinician at the time of echocardiography depending on inferior vena cava diameter and respiratory variation. Mean PAP was then estimated as mean PAP = 0.61 peak PAP + 2 (Chemla et al., 2024). As the model was developed as a diagnostic tool it was assumed that the RHC invasive pressure readings were not available for tuning the 3EWK parameters. Therefore, the parameters were implemented in the models as first calculated with the input patient-specific average 4D-flow-derived flow and echocardiogram-derived pressure values. All simulations were run for a minimum of three cardiac cycles to ensure the outlet flow and pressure values being calculated by the 3EWK model became periodic, with a variation of less than 1% between the final two cycles. Table 1 reports key flow information of the inlet velocity profile as well as the 3EWK parameters for all patients.

Table 1
www.frontiersin.org

Table 1. Key flow information *extracted from 4D-flow MRI and **measured from right heart catheter, and 3EWK parameters for patients P1-7.

Previous literature has reported that transition to turbulence in the aorta occurs at a critical value of 150α, where α is the Womersley number defined as 2π/period (Kousera et al., 2012). As similar data has not been reported for the pulmonary artery, the same criterion was applied in this study due to the anatomical and functional similarities of the aorta and pulmonary artery. Thus, peak Reynolds number was calculated for all models and was found to be below the critical value for transition to turbulence. Therefore, flow was assumed to be and modelled as laminar in all cases.

The non-Newtonian behaviour of the blood was captured through the Quemada model, with model parameters taken from literature (Mimouni, 2016). All simulations were run with a timestep of 0.005 s in CFX (Ansys, v22) for a minimum of three cardiac cycles to ensure a periodic solution. The final cycle was used for analysis, with post-processing done in EnSight (Ansys, v22). Qualitative and quantitative comparisons of velocity and flow fields between CFD and 4D-MRI data were first conducted for validation. Further analysis of the numerical results was then conducted to evaluate TAWSS values and patterns. The CFD calculated TAWSS values were then correlated against RHC measured mean PAP, as this is the key parameter required for diagnosis of PAH. Any correlation identified between TAWSS and mean PAP would thus indicate a potential non-invasive computational biomarker for PAH.

3 Results

3.1 Validation of CFD against 4D-flow MRI

Figure 2 presents peak systolic velocity fields from 4D-flow MRI data and CFD models. Across all patients the CFD predicted velocities are well matched to the in-vivo 4D-flow data. All patients presented with a clear velocity jet entering the main pulmonary artery, however, the variation in jet distribution and size is clear, further highlighting the importance of a patient-specific inlet velocity profile. P3 observed higher velocities through the pulmonary artery compared to all patients, represented by the increased velocity scale (upper limit of 0.9 m/s vs. 0.5 m/s). This is expected due to P3’s high stroke volume (183 mL). Not only do the CFD models capture well the intricate non-uniform velocity jets in the main pulmonary artery that vary significantly across patients, but also, the unique distribution of flow to the LPA and RPA in each patient is clearly observed. For example, local regions of increased velocities match well between CFD and 4D-flow data in the LPA and RPA of P1, and the RPA of P6 and P7.

Figure 2
www.frontiersin.org

Figure 2. Peak systolic velocity fields for all patients, extracted from 4D-flow MRI data and modelled with CFD. P3 has been separated due to a different velocity scale.

To quantitatively validate the CFD models, the percentage of flow reporting to each branch and the mean peak systolic velocity on a plane in the centra of the main pulmonary artery were measured from both the 4D-flow MRI scans and CFD simulations–Table 2 reports this data. On average CFD velocities in the main pulmonary artery were within 10% of 4D-flow MRI velocities (range 5%–19%), and all CFD models simulated the LPA/RPA distribution within 3.5% of the 4D-flow measurement. The strong quantitative and qualitative agreement in velocity fields between the CFD and 4D-flow MRI data gives confidence in the following step of evaluating derivative hemodynamic parameters from the CFD data.

Table 2
www.frontiersin.org

Table 2. Comparison of hemodynamic parameters measured from 4D-flow MRI and calculated from CFD simulation.

3.2 CFD predicted WSS

Figure 3 shows TAWSS contours for all seven patients. As with the velocity results, varying TAWSS patterns are observed across the cohort. P2 and P4 observed the lowest TAWSS throughout the pulmonary artery; P1 and P5-7 exhibited similar levels of TAWSS, with unique localisations indicating the locations where the pulmonary valve velocity jets were impinging on the wall and diverting to the LPA and RPA. Again, in line with the velocities, P3 observed significantly higher TAWSS than the other patients–this would be expected given the increased velocities in volume of similar size.

Figure 3
www.frontiersin.org

Figure 3. CFD predicted time-average wall shear stress (TAWSS) patterns for all patients. P3 has been separated due to a different TAWSS scale.

Spatially averaged TAWSS values over the whole geometry (main pulmonary artery and initial segment of the LPA and RPA) and in the main pulmonary artery alone (proximal to the bifurcation) are reported in Table 3. TAWSS ranged from 0.63 to 3.62 across the seven patients, with P2 seeing the lowest TAWSS values and P3 the highest. All patients saw higher TAWSS when taking the LPA and RPA into account, however the magnitude of influence of the inclusion of the LPA and RPA varied: P2, P4, and P5 saw less than 0.2 Pa increase, P1, P6, and P7 an increase of 0.3–0.4 Pa, and P3 and increase of 1 Pa.

Table 3
www.frontiersin.org

Table 3. Spatially averaged, over two regions of the geometric model, TAWSS for patients P1-7.

Both sets of TAWSS values (main pulmonary artery alone, and then with the addition of LPA and RPA) were tested for correlations with invasively measured RHC PAP (values reported in Table 1). A linear correlation was observed for mean PAP vs. TAWSS in the main pulmonary artery + LPA + RPA region, with R2 = 0.84 (Figure 4). Considering the main pulmonary artery segment alone, the R2 was 0.77.

Figure 4
www.frontiersin.org

Figure 4. Mean time-averaged wall shear stress (TAWSS) estimated by patient-specific CFD models correlated with mean pulmonary arterial pressure (PAP) measured by right heart catheter.

4 Discussion

PAH is a challenging disease to diagnose, treat, and manage. Symptoms can often be subtle and, as PAH requires a technical and invasive RHC procedure for diagnosis, referral to a specialist centre is required. With this, patients can be interacting with doctors for years before referral (Kiely et al., 2019a). This extended period between symptom onset and diagnosis leads to worse patient outcomes and higher healthcare costs (Lawrie et al., 2024). Yet a wealth of data is available on these patients due to the numerous interactions with their local healthcare system. This provides an opportunity to develop novel non-invasive diagnostic methods based on this data to better identify patients at risk of PAH and improve early diagnosis numbers.

In this study, a methodology to build patient-specific 4D-flow MRI based CFD models, previously developed and validated in aortic dissection (Armour et al., 2022), was applied to a pilot cohort of PAH patients. Pulmonary artery hemodynamics of seven patients were modelled and analysed to identify hemodynamics parameters that are linked to mean PAP (the key diagnostic parameter for PAH).

The velocity fields for all patients shown in Figure 2 were evaluated in detail. The variation in pulmonary valve velocity jet and resulting TAWSS contours demonstrates how individual pulmonary hemodynamics are, and the importance of developing patient-specific CFD models. Within the cohort the size, distribution, and magnitude of the main jet varies, as well as the distribution to the LPA and RPA. From Figure 2, P1 and P6 observe significantly different velocities between their own LPA and RPA, yet P1sees higher TAWSS in the RPA where velocities are higher while P6 sees higher TAWSS in the LPA where velocities are lower. These unbalanced velocities and TAWSS between the LPA and RPA in P1 and P6 correspond to the relatively uneven LPA/RPA flow split measured from the 4D-flow data as reported in Table 1. However, P7 saw the most uneven LPA/RPA flow split, with 69% of inlet flow reporting to the RPA, yet there were less drastic variations in velocities.

The strong match between the CFD simulation velocity fields and 4D-flow MRI data allowed for CFD calculated TAWSS patterns, and importantly magnitudes, to be further interpreted and linked to patient-specific clinical data. This resulted in the identification of a linear correlation between TAWSS and mean PAP of 0.84. Given the sample size of seven patients, this cannot be classed as statistically significant, however, it is a clear indication of a potential computational biomarker of PAH. This aligns with the findings of Schafer et al. who evaluated TAWSS directly from 4D-flow MRI. Importantly, the presented methodology of using CFD to evaluate TAWSS opens up the potential to focus on the difficult area of patients with PAPs near the 20 mmHg diagnostic threshold. With more confidence in exact magnitudes of TAWSS estimated from CFD (as opposed to 4D-flow MRI) larger cohort studies can identify incremental thresholds that better distinguish PAP values and PAH severity.

Interestingly, the correlation between TAWSS and mean PAP slightly reduced to 0.77 when TAWSS was evaluated in the main pulmonary artery alone (excluding the beginning of the LPA and RPA). This may suggest that the bifurcation, both in terms of geometry, alignment, and possibly biomechanical properties is related to high pressure. Exactly how and where in the development of the disease this is related requires further study. In studies that have focused on vortical flow (Cerne et al., 2022; Sabry et al., 2023) vortices were seen to gather and dissipate around the bifurcation area. Such vortical flow would influence TAWSS, further strengthening the hypothesis that the pulmonary artery bifurcation is influential in hemodynamics relating to PAH.

As WSS is dependent on velocity gradient, both vessel diameter and stroke volume can significantly impact WSS values due to their influence on velocity. Therefore, correlations of WSS against both vessel diameter and stroke volume (as a surrogate of velocity) were evaluated. Vessel diameter showed a weak R2 value of 0.36 against mean PAP, suggesting that vessel diameter alone cannot predict mean PAP. Stroke volume reported a R2 value of 0.69 against mean PAP. This is a stronger correlation than diameter yet lower still than the R2 value of 0.84 for the CFD derived TAWSS. This suggests that while stroke volume is potentially more closely linked to PAP, there are more factors influencing PH than stroke volume alone–TAWSS also encapsulates variations in the 3D geometry of the pulmonary artery.

Finally, as mentioned in the introduction, TAWSS has been estimated directly from 4D-flow MRI in previous studies. The limitation with these 4D-flow values arises due to the low spatial resolution. For example, a 2D plane cut in the main pulmonary artery of patient P1 contained approximately 4700 and 260 mesh elements for the CFD model and 4D-flow MRI, respectively. The significantly enhanced resolution of the CFD model allows for flow patterns and velocity gradients to be captured in greater detail, particularly near the wall. In the 4D-flow MRI data small deviations in velocities will be lost due to the large voxels, resulting in lower velocity gradients and therefore lower TAWSS values.

As a comparison, TAWSS was estimated directly from the 4D-flow MRI data on planes placed in the MPA, RPA, and LPA using cvi42 (Circle Cardiovascular Imaging Inc.). The average and range of these values among the seven patients are reported in Table 4, along with the corresponding CFD values. The 4D-flow values are substantially lower than the CFD values and cover a smaller range across all three planes (4D-flow MRI: 0.04–0.19 Pa vs. CFD: 0.45–3.59 Pa). Among the seven patients, the 4D-flow MRI values showed no correlation with RHC measured mean PAP (R2 < 0.1). However, this is a small sample set and evaluation in a larger patient cohort will allow for more confident conclusions. Future work will also look to quantitatively compare TAWSS distributions across the entire pulmonary artery wall.

Table 4
www.frontiersin.org

Table 4. Averaged (and range of) TAWSS among the seven patients (P1-7) on planes placed in the main (MPA), right (RPA), and left pulmonary artery (LPA) for the CFD model and 4D-flow MRI data.

A key limitation of this work is the small number of patients included. In this case, this pilot study cohort has been used to identify potential computational biomarkers for PAH diagnosis and severity classification, and further work is ongoing to validate these findings in a significantly larger cohort. In this future work, healthy and non-PAH cardiovascular disease controls will be assessed to overcome the limitation of the current study being limited to PAH only. In terms of the computational model, the assumption of a rigid wall will impact simulated hemodynamics. However, it is known that increased pulmonary artery stiffness comes with PAH (Schafer et al., 2016), thus the impact on WSS values of the rigid wall will likely be reduced, compared to a healthy model. Furthermore, modelling wall motion significantly increase computational complexity, and reduces clinical deployment feasibility. Provided methodology is consistent, suitable rigid model WSS thresholds identified in a larger cohort study can still be used to stratify patients. Finally, the extent of the geometric model was limited to the main pulmonary artery and initial section of the LPA and RPA due to the scan region of the 4D-flow MRI. In future work, segmentation from anatomical CT or MRI will be incorporated to extend the geometric model to the second branch of the pulmonary tree to fully assess TAWSS in the LPA and RPA.

The results of this study demonstrate a clear indication that WSS is linked to pulmonary artery pressure in PAH patients. Specifically, higher pulmonary artery pressure is associated with lower TAWSS. These findings lay the groundwork for future work to develop non-invasive computational models for improved early diagnosis and risk stratification in PAH.

Data availability statement

All processed data, that complies with data governance and ethics, supporting the conclusions of this article will be shared by authors upon reasonable request.

Ethics statement

The studies involving humans were approved by UK National Health Service - Health Research Authority. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

CA: Conceptualization, Formal Analysis, Funding acquisition, Methodology, Writing – original draft, Writing – review and editing. DG: Data curation, Writing – review and editing. BS: Data curation, Writing – review and editing. DO’: Data curation, Writing – review and editing. LH: Data curation, Writing – review and editing. MW: Data curation, Writing – review and editing. XX: Supervision, Writing – review and editing. AL: Supervision, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. CHA is supported by the Imperial College London British Heart Foundation Centre of Research Excellence (RE/18/4/34215) and the British Heart Foundation (FS/IPBSRF/24/27116). DO’R is supported by the Medical Research Council (MC_UP_1605/13); National Institute for Health Research (NIHR) Imperial College Biomedical Research Centre; and the British Heart Foundation (RG/19/6/34387, RE/24/130023, CH/P/23/80008).

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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

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

Armour, C. H., Guo, B., Pirola, S., Saitta, S., Liu, Y., Dong, Z., et al. (2020). The influence of inlet velocity profile on predicted flow in type B aortic dissection. Biomech. Mod. Mechanbiol 20, 481–490. doi:10.1007/s10237-020-01395-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Armour, C. H., Guo, B., Saitta, S., Pirola, S., Liu, Y., Dong, Z., et al. (2022). Evaluation and verification of patient-specific modelling of type B aortic dissection. Comp. Bio Med. 140, 105053. doi:10.1016/j.compbiomed.2021.105053

PubMed Abstract | CrossRef Full Text | Google Scholar

Armour, C. H., Salmasi, M. Y., Xu, X. Y., and Pirola, S. (2024). Chapter 19 - modeling of flow-induced mechanosignaling. Biomechanics Living Organs, Biomechanics Aorta, 431–448. doi:10.1016/B978-0-323-95484-6.00004-X

CrossRef Full Text | Google Scholar

Arshad, M., Ghim, M., Mohamied, Y., Sherwin, S. J., and Weinberg, P. D. (2021). Endothelial cells do not align with themean wall shear stress vector. J. R. Soc. Interface 18, 20200772. doi:10.1098/rsif.2020.0772

PubMed Abstract | CrossRef Full Text | Google Scholar

Barker, A. J., Roldán-Alzate, A., Entezari, P., Shah, S. J., Chesler, N. C., Wieben, O., et al. (2015). Four-dimensional flow assessment of pulmonary artery flow and wall shear stress in adult pulmonary arterial hypertension: results from two institutions. Magn. Reson Med. 73, 1904–1913. doi:10.1002/mrm.25326

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergemann, R., Allsopp, J., Jenner, H., Daniels, F. A., Drage, E., Samyshkin, Y., et al. (2018). High levels of healthcare utilization prior to diagnosis in idiopathic pulmonary arterial hypertension support the feasibility of an early diagnosis algorithm: the SPHInX project. Pulm. Circ. 8, 1–9. doi:10.1177/2045894018798613

PubMed Abstract | CrossRef Full Text | Google Scholar

Bordones, A. D., Leroux, M., Kheyfets, V. O., Wu, Y. A., Chen, C. Y., and Finol, E. A. (2018). Computational fluid dynamics modeling of the human pulmonary arteries with experimental validation. Ann. Biomed. Eng. 46 (9), 1309–1324. doi:10.1007/s10439-018-2047-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Cerne, J. W., Pathrose, A., Gordon, D. Z., Sarnari, R., Veer, M., Blaisdell, J., et al. (2022). Evaluation of pulmonary hypertension using 4D flow MRI. J. Magn. Reson Imaging 56, 234–245. doi:10.1002/jmri.27967

PubMed Abstract | CrossRef Full Text | Google Scholar

Chemla, D., Castelain, V., Humbert, M., Hébert, J. L., Simonneau, G., Lecarpentier, Y., et al. (2004). New formula for predicting mean pulmonary artery pressure using systolic pulmonary artery pressure. Chest 126 (4), 1313–1317. doi:10.1378/chest.126.4.1313

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiu, J.-J., Chen, L.-J., Chang, S.-F., Lee, P.-L., Lee, C.-I., Tsai, M. C., et al. (2005). Shear stress inhibits smooth muscle cell induced inflammatory gene expression in endothelial cells role of NF-κB. Arterioscler. Thromb. Vasc. Biol. 25, 963–969. doi:10.1161/01.atv.0000159703.43374.19

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiu, J.-J., Chen, L.-J., Lee, C.-I., Lee, P.-L., Lee, D.-Y., Tsai, M. C., et al. (2007). Mechanisms of induction of endothelial cell E-selectin expression by smooth muscle cells and its inhibition by shear stress. Blood 110, 519–528. doi:10.1182/blood-2006-08-040097

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiu, J.-J., Lee, P.-L., Chen, C.-N., Lee, C.-I., Chang, S.-F., Chen, L.-J, et al. (2004). Shear stress increases ICAM-1 and decreases VCAM-1 and E-selectin expressions induced by tumor necrosis factor-in endothelial cells. Arterioscler. Thromb. Vasc. Biol. 24, 73–79. doi:10.1161/01.atv.0000106321.63667.24

PubMed Abstract | CrossRef Full Text | Google Scholar

Deux, J. F., Crowe, L. A., Genecand, L., Hachulla, A. L., Glessgen, C., Noble, S., et al. (2022). Correlation between pulmonary artery pressure and vortex duration determined by 4D flow MRI in main pulmonary artery in patients with suspicion of chronic thromboembolic pulmonary hypertension (CTEPH). J. Clin. Med. 11, 5237. doi:10.3390/jcm11175237

PubMed Abstract | CrossRef Full Text | Google Scholar

Fanni, B. M., Gasparotti, E., Vignali, E., Capelli, C., Positano, V., and Celi, S. (2023). An integrated in-vitro and in-silico workflow to study the pulmonary bifurcation hemodynamics. Comp. Fluids 260, 105912. doi:10.1016/j.compfluid.2023.105912

CrossRef Full Text | Google Scholar

Fisher, M. R., Forfia, P. R., Chamera, E., Housten-Harris, T., Champion, H. C., Girgis, R. E., et al. (2009). Accuracy of Doppler echocardiography in the hemodynamic assessment of pulmonary hypertension. Am. J. Respir. Crit. Care Med. 179, 615–621. doi:10.1164/rccm.200811-1691oc

PubMed Abstract | CrossRef Full Text | Google Scholar

Hastings, N. E., Simmers, M. B., McDonald, O. G., Wamhoff, B. R., and Blackman, B. R. (2007). Atherosclerosis-prone hemodynamics differentially regulates endothelialand smooth muscle cell phenotypes and promotes pro-inflammatory priming. Am. J. Physiol. Cell Physiol. 293, C1824–C1833. doi:10.1152/ajpcell.00385.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

Humbert, M., Kovacs, G., Hoeper, M. M., Badagliacca, R., Berger, R. M., Brida, M., et al. (2023). 2022 ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension. Eur. Respir. J. 61 (1), 2200879. doi:10.1183/13993003.00879-2022

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiely, D. G., Doyle, O., Drage, E., Jenner, H., Salvatelli, V., Daniels, F. A., et al. (2019a). Utilising artificial intelligence to determine patients at risk of a rare disease: idiopathic pulmonary arterial hypertension. Pulm. Circ. 9, 1–9. doi:10.1177/2045894019890549

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiely, D. G., Lawrie, A., and Humbert, M. (2019b). Screening strategies for pulmonary arterial hypertension. Eur. Heart J. Suppl. 21, K9–K20. doi:10.1093/eurheartj/suz204

PubMed Abstract | CrossRef Full Text | Google Scholar

Kousera, C. A., Wood, N. B., Seed, W. A., Torii, R., O'Regan, D., and Xu, X. Y. (2012). A numerical study of aortic flow stability and comparison with in vivo flow measurements. J. Biomech. Eng. 135 (1), 011003. doi:10.1115/1.4023132

PubMed Abstract | CrossRef Full Text | Google Scholar

Lawrie, A., Hamilton, N., Wood, S., Exposto, F., Muzwidzwa, R., Raiteri, L., et al. (2024). Association of risk assessment at diagnosis with healthcare resource utilization and health-related quality of life outcomes in pulmonary arterial hypertension. Pulm. Circ. 14, e12399. doi:10.1002/pul2.12399

PubMed Abstract | CrossRef Full Text | Google Scholar

Mimouni, Z. (2016). The rheological behavior of human blood—comparison of two models. Open J. Biophys. 6 (2), 29–33. doi:10.4236/ojbiphy.2016.62004

CrossRef Full Text | Google Scholar

Miyazaki, S., Itatani, K., Furusawa, T., Nishino, T., Sugiyama, M., Takehara, Y., et al. (2017). Validation of numerical simulation methods in aortic arch using 4D Flow MRI. Heart Vessels 32, 1032–1044. doi:10.1007/s00380-017-0979-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Pillalamarri, N. R., Piskin, S., Patnaik, S. S., Murali, S., and Finol, E. A. (2021). Patient-specific computational analysis of hemodynamics in adult pulmonary hypertension. Ann. Biomed. Eng. 49, 3465–3480. doi:10.1007/s10439-021-02884-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Pirola, S., Jarral, O. A., O’Regan, D. P., Asimakopoulos, G., Anderson, J. R., Pepper, J. R., et al. (2017). Computational study of aortic hemodynamics for patients with an abnormal aortic valve: the importance of secondary flow at the ascending aorta inlet. Apl. Bioeng. 2 (2), 026101. doi:10.1063/1.5011960

CrossRef Full Text | Google Scholar

Reiter, U., Kovacs, G., Reiter, C., Kräuter, C., Nizhnikava, V., Fuchsjäger, M., et al. (2020). MR 4D flow-based mean pulmonary arterial pressure tracking in pulmonary hypertension. Eur. Radiol. 31, 1883–1893. doi:10.1007/s00330-020-07287-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Rich, J. D., Shah, S. J., Swamy, R. S., Kamp, A., and Rich, S. (2011). Inaccuracy of Doppler echocardiographic estimates of pulmonary artery pressures in patients with pulmonary hypertension. Chest 139 (5), 988–993. doi:10.1378/chest.10-1269

PubMed Abstract | CrossRef Full Text | Google Scholar

Sabry, M., Lamata, P., Sigfridsson, A., Kermati, H., Fyrdahl, A., Ugander, M., et al. (2023). Vortex duration time to infer pulmonary hypertension: in-silico emulation and dependence on quantification technique functional imaging and modeling of the heart, FIMH 2023. LNCS 13958, 425–434. doi:10.1007/978-3-031-35302-4_44

CrossRef Full Text | Google Scholar

Saitta, S. (2022). saitta-s/flow4D: flow4D. Zenodo. doi:10.5281/zenodo.7236015

CrossRef Full Text | Google Scholar

Saitta, S., Maga, L., Armour, C., Votta, E., O'Regan, D. P., Salmasi, M. Y., et al. (2023). Data-driven generation of 4D velocity profiles in the aneurysmal ascending aorta. Comp. Meth. Prog. Biomed. 233, 107468. doi:10.1016/j.cmpb.2023.10746

PubMed Abstract | CrossRef Full Text | Google Scholar

Schafer, M., Kheyfets, V. O., Schroeder, J. D., Dunning, J., Shandas, R., Buckner, J. K., et al. (2016). Main pulmonary arterial wall shear stress correlates with invasive hemodynamics and stiffness in pulmonary hypertension. Pulm. Circ. 6 (1), 37–45. doi:10.1086/685024

PubMed Abstract | CrossRef Full Text | Google Scholar

Simonneau, G., Montani, D., Celermajer, D. S., Denton, C. P., Gatzoulis, M. A., Krowka, M., et al. (2019). Haemodynamic definitions and updated clinical classification of pulmonary hypertension. Eur. Respir. J. 53, 1801913. doi:10.1183/13993003.01913-2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, B. T., Pickard, S. S., Chan, F. P., Tsao, P. S., Taylor, C. A., and Feinstein, J. A. (2012). Wall shear stress is decreased in the pulmonary arteries of patients with pulmonary arterial hypertension: an image-based, computational fluid dynamics study. Pulm. Circ. 2 (4), 470–476. doi:10.4103/2045-8932.105035

PubMed Abstract | CrossRef Full Text | Google Scholar

Testani, J. M., St. John Sutton, M. G., and Wiegers, S. E. (2010). Accuracy of noninvasively determine pulmonary artery systolic pressure. Am. J. Cardiol. 105, 1192–1197. doi:10.1016/j.amjcard.2009.11.048

CrossRef Full Text | Google Scholar

Tsuabata, H., Nakanishi, N., Itatani, K., Takigami, M., Matsubara, Y., Ogo, T., et al. (2019). Pulmonary blood flow dynamics analyzed by computational fluid dynamics in the patients with chronic thromboembolic pulmonary hypertension before and after balloon pulmonary angioplasty. Circ 138, A11901. doi:10.1038/s41598-023-33727-6

CrossRef Full Text | Google Scholar

Tsuabata, H., Nakanishi, N., Itatani, K., Takigami, M., Matsubara, Y., Ogo, T., et al. (2023). Pulmonary artery blood flow dynamics in chronic thromboembolic pulmonary hypertension. Sci. Rep. 13, 6490. doi:10.1038/s41598-023-33727-6

CrossRef Full Text | Google Scholar

Wang, C., Baker, B. M., Chen, C. S., and Schwartz, M. A. (2013). Endothelial cell sensing of flow direction. Arterioscler. Thromb. Vasc. Biol. 33, 2130–2136. doi:10.1161/atvbaha.113.301826

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, K., Armour, C. H., Guo, B., Dong, Z., and Xu, X. Y. (2024). A new method for scaling inlet flow waveform in hemodynamic analysis of aortic dissection. Num Meth Biomed. Eng. 40 (9), e3855. doi:10.1002/cnm.3855

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, Z., Armour, C. H., Kandail, H., O’Regan, D. P., Bahrami, T., Mirsadraee, S., et al. (2024). Fluid–structure interaction analysis of a healthy aortic valve and its surrounding haemodynamics. Num Meth Biomed. Eng. 40 (11), e3865. doi:10.1002/cnm.3865

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Rothenberger, S. M., Brindise, M. C., Markl, M., Rayz, V. L., and Vlachos, P. P. (2022). Wall shear stress estimation for 4D flow MRI using Navier–Stokes equation correction. Ann. Biomed. Eng. 50 (12), 1810–1825. doi:10.1007/s10439-022-02993-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zilber, Z. A., Boddu, A., Malaisrie, S. C., Hoel, A. W., Mehta, C. K., Vassallo, P., et al. (2021). Noninvasive morphologic and hemodynamic evaluation of type B aortic dissection: state of the art and future perspectives. Radiol. Cardiothorac. Imag. 3 (3), e200456. doi:10.1148/ryct.2021200456

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: pulmonary arterial hypertension, computational fluid dynamics, 4D-flow MRI, wall shear stress, non-invasive assessment, computational biomarker

Citation: Armour CH, Gopalan D, Statton B, O’Regan DP, Howard L, Wilkins MR, Xu XY and Lawrie A (2025) Patient-specific modelling of pulmonary arterial hypertension: wall shear stress correlates with disease severity. Front. Bioeng. Biotechnol. 13:1585345. doi: 10.3389/fbioe.2025.1585345

Received: 28 February 2025; Accepted: 03 June 2025;
Published: 18 June 2025.

Edited by:

Harvey Ho, University of Auckland, New Zealand

Reviewed by:

Shahin Moledina, University College London, United Kingdom
Peshala Thibbotuwawa Gamage, Florida Institute of Technology, United States

Copyright © 2025 Armour, Gopalan, Statton, O’Regan, Howard, Wilkins, Xu and Lawrie. 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: C. H. Armour, Y2hsb2UuYXJtb3VyQGltcGVyaWFsLmFjLnVr

These authors share senior authorship

Disclaimer: 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.