Next Article in Journal
Early Postnatal Neuroinflammation Produces Key Features of Diffuse Brain White Matter Injury in Rats
Next Article in Special Issue
Neurosurgical Microvascular Anastomosis: Systematic Review of the Existing Simulators and Proposal of a New Training Classification System
Previous Article in Journal
White Matter Magnetic Resonance Diffusion Measures in Multiple Sclerosis with Overactive Bladder
Previous Article in Special Issue
Assessment of Interrater Reliability and Accuracy of Cerebral Aneurysm Morphometry Using 3D Virtual Reality, 2D Digital Subtraction Angiography, and 3D Reconstruction: A Randomized Comparative Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fluid–Structure Interaction Simulations of the Initiation Process of Cerebral Aneurysms

1
eulerian-solutions e.U., Leonfeldnerstraße 245, 4040 Linz, Austria
2
Unit Medical Informatics, RISC Software GmbH, Softwarepark 32a, 4232 Hagenberg, Austria
3
Institute of Polymer Product Engineering, Johannes Kepler University Linz, Altenberger Strasse 69, 4040 Linz, Austria
4
Medical Faculty, Johannes Kepler University Linz, Altenberger Strasse 69, 4020 Linz, Austria
5
Department of Neurosurgery, Kepler University Hospital, Wagner-Jauregg-Weg 15, 4020 Linz, Austria
6
Clinical Research Institute for Neuroscience, Johannes Kepler University Linz, 4020 Linz, Austria
*
Author to whom correspondence should be addressed.
Brain Sci. 2024, 14(10), 977; https://doi.org/10.3390/brainsci14100977
Submission received: 25 August 2024 / Revised: 24 September 2024 / Accepted: 25 September 2024 / Published: 27 September 2024

Abstract

:
Background: Hemodynamics during the growth process of cerebral aneurysms are incompletely understood. We developed a novel fluid–structure interaction analysis method for the identification of relevant scenarios of aneurysm onset. Method: This method integrates both fluid dynamics and structural mechanics, as well as their mutual interaction, for a comprehensive analysis. Patients with a single unruptured cerebral aneurysm were included. Results: Overall, three scenarios were identified. In scenario A, wall shear stress (WSS) was low, and the oscillatory shear index (OSI) was high in large areas within the region of aneurysm onset (RAO). In scenario B, the quantities indicated a reversed behavior, where WSS was high and OSI was low. In the last scenario C, a behavior in-between was found, with scenarios A and B coexisting simultaneously in the RAO. Structural mechanics demonstrated a similar but independent trend. Further, we analyzed the change in hemodynamics between the onset and a fully developed aneurysm. While scenarios A and C remained unchanged during aneurysm growth, 47% of aneurysms in scenario B changed into scenario A and 20% into scenario C. Conclusions: In conclusion, these findings suggest that WSS and the OSI are reciprocally regulated, and both low and high WSS/OSI conditions can lead to aneurysm onset.

1. Introduction

Cerebral aneurysms affect 2–5% of the general population [1]. An aneurysm rupture may cause a severe subarachnoid hemorrhage (SAH) with high mortality and morbidity rates [2]. Several studies have demonstrated various clinical, genetic, morphological, or hemodynamic factors contributing to the initiation, growth, and rupture of cerebral aneurysms [3].
Recent computational fluid dynamics (CFD) investigations revealed significantly different hemodynamics in ruptured versus unruptured aneurysms [4,5]. Hemodynamic phenomena can provide mechanical triggers, which are transduced into biological signals, leading to possible aneurysm growth and rupture [6].
Numerous attempts have been made to study the influence of the wall shear stress (WSS), one of the best-characterized fluid-dynamic parameters, on these processes [5,6,7,8]. Meng et al. described a novel view of aneurysm development, proposing that WSS plays a key mechanistic role, as suggested by their high-versus-low WSS theory [9]. In other studies, similar investigations were carried out, with conflicting results [10].
If the exact mechanisms of aneurysm formation could be determined, aneurysm rupture might be predicted and prevented [7]. Additionally, aneurysm growth may occur in up to 18% of affected patients and may be associated with aneurysm destabilization and rupture [2]. Similarly, other studies found a 12-fold increase in rupture risk in growing aneurysms [11].
Currently, several studies have analyzed hemodynamics in aneurysm initiation [7,8,12]. Although many detailed fluid dynamics simulation analyses, including complex flow patterns, were conducted [13,14,15,16,17,18,19,20,21,22], only a few studies have focused on structural dynamics calculations [23,24,25,26]. Most of these studies focused on abdominal aneurysms [24,25,26]. However, a comprehensive evaluation of individual hemodynamics should include the analysis of fluid as well as structural aspects using fluid–structure interaction (FSI) analysis. Therefore, basic research remains important for accurately characterizing and understanding mechanisms of cerebrovascular diseases to enable patient-specific treatment strategy optimization.
The aim of this study is to develop a novel FSI analysis method and investigate the impacts of both hemodynamics and structural mechanics on the initiation and growth of cerebral aneurysms in order to better understand the controversial dynamics of the initiation and growth process of cerebral aneurysms.

2. Methods

2.1. Patient Data

Patients with a single unruptured cerebral aneurysm treated at the Department of Neurosurgery, Kepler University Hospital Linz in 2020 were included for the identification of characteristic scenarios regarding the initiation and growth of cerebral aneurysms. Data were collected retrospectively, and fusiform and dissecting aneurysms were excluded. Additionally, we excluded large aneurysms with a size > 10 mm (n = 4) to minimize the effect of aneurysm size on the FSI analysis, as previously reported [7,27]. Overall, 44 patients were included (26 females, mean age 54 years; ranging from 34–76 years). Microsurgical or endovascular treatments were performed in 19 and 25 patients, respectively. The aneurysms were located at the internal carotid artery (n = 13), middle cerebral artery (MCA; n = 16), anterior communicating artery (n = 9), and posterior circulation (n = 6). This study was approved by the local ethics committee (Ethikkommission der medizinischen Fakultät der Johannes Kepler Universität; EK Nr: 1129/2022), and the requirement for the acquisition of informed consent from patients was waived, owing to the retrospective nature of the research.

2.2. Patient Image Data

We extracted the aneurysm geometries from medical image data obtained via digital subtraction angiography (DSA). Since these images showed only the cerebral blood vessels with high contrast and no other anatomical structures, segmentation was relatively straightforward and could be performed with intensity thresholding and minor editing of the voxel volume. A hydrodynamic entrance length of 5–10 times the inlet and outlet diameters was utilized for the inflow and outflow vessels, respectively, in order to define the cutoff positions during the reconstruction of the vascular geometry [28,29]. In case the quality of the DSA images did not allow such a length, the maximum possible length was taken. The segmentation was then converted into a surface mesh (STL format), the centerlines of the mesh were calculated, and inlet as well as outlet planes were placed perpendicular to these lines.
There are multiple alternatives to investigating the aneurysms onset. Possibly the most ideal approach is the investigation of de novo aneurysms, as shown with 10 aneurysms in [7]. Often, it is difficult to obtain patient data of a large number of de novo aneurysms; thus, in literature, a retrospective approach is used alternatively, where aneurysms are manually removed from the parent vessel. While this approach may have some limitations compared to the first approach, it is still mostly utilized in literature [8,27,30,31,32,33,34].
In this work, the retrospective approach was utilized, due to the available patient images. The CAD models in the form of the generated STL files of the aneurysms (see Figure 1 left) were modified in a subsequent step by removing the aneurysm sac with the help of a Laplacian smoother algorithm [35,36,37,38] in order to obtain the geometry of the parent vessel without the aneurysm sac (see Figure 1 center).
In Figure 1, an estimated black line is shown, which indicates the approximate region of aneurysm onset (RAO). This region is important for understanding the initiation process of aneurysm growth.
The generated STL files (both with and without the aneurysm sac) were used to generate the volumetric calculation grid for the FSI simulations. In the first step, the fluid domain was established, and in the second subsequent step, the wall domain was created, with a constant wall thickness of 0.3 mm [22], by extruding and creating a separate vessel wall domain with the given thickness.

2.3. Hemodynamic and Structural Mechanical Modeling

A finite volume solver [39], based on the computational fluid dynamics framework OpenFOAM v2012 [40,41], alongside the fluid–structure interaction library solids4Foam [42], was utilized to numerically solve the unsteady equations of hemodynamics in the simulation process. The models applied the principles of mass and momentum conservation through both the continuity equation and the Navier–Stokes equations. For structural mechanics, Newton’s second law was employed. These equations were solved in a sequential manner and were coupled with a fluid–structure interaction (FSI) boundary condition. The details, as well as the experimental validation of the simulation methods, were described in [39].
The fluid dynamic boundary conditions at the inflow were defined by a pulsatile flow profile, characterized by a temporal velocity curve derived from published data [43]. The outflow conditions were set to a time-dependent value defined by experimental pressure measurements [44]. Boundary conditions along the interior vessel walls assumed a no-slip condition. Blood was represented in the model as a Newtonian fluid, with viscosity set at 0.04 poise and a constant density of 1.06 g/cm3. The simulation lasted for a single cardiac cycle of 1 s, corresponding to a heart rate of 60 beats per minute, discretized into 100 temporal increments (Δt = 0.01 s).
The vessel walls at the in- and outflow faces were fixed on the solid mechanical side. The vessel wall was assumed to be linear elastic, with a Young’s modulus value E of 2490 kPa [45] and a Poisson’s ratio ν of 0.49, which assumed an almost incompressible material, as described in [25,46]. Results for this Young’s modulus are shown in Figure 2, Figure 3 and Figure 4. We changed the Young’s modulus from the defined value of 2490 kPa from [45] to 5700 and 500 kPa [23], respectively, to assess the influence of vessel wall stiffness.
The investigated parameter of the oscillatory shear index OSI was defined as:
O S I = 0.5 1 0 T τ d t 0 T τ d t ,
where τ is the wall shear stress.
Investigated cohorts did not show normal distribution; thus, standard deviation was not given, and the significance of the difference between the means was determined with a Mann–Whitney U test (p-values < 0.05).

2.4. Convergence Analysis

In order to guarantee the accuracy of the simulations while maintaining a reasonable simulation runtime, the results of a convergence analysis are presented with the exemplary aneurysm in Figure 1. For this, the following numeric parameters were analyzed:
  • Time step independence;
  • Calculation mesh convergence;
  • Newtonian vs. non-Newtonian fluid;
  • One cardiac cycle vs. five cardiac cycles.
Table 1 shows the convergence results, including the required runtime for a simulation. A good time step convergence can be found, and a time step size of 0.01 s was utilized in all simulations, guaranteeing accuracy and also a reasonable simulation runtime.
The change in the fluid volume element edge length delivered a good convergence at 0.2 mm. The given values were target values for the meshing algorithm, and due to the irregularity of the geometry, small local deviations were possible.
The change in the solid volume element size was defined by the utilized number of cells over the thickness of the vessel wall; as in all simulations, a constant thickness of 0.3 mm [22] was utilized. Three cells over the thickness offered the best compromise between accuracy and simulation runtime.
The influence of a non-Newtonian viscosity model was also compared to the constant-viscosity approach by utilizing the Carreau–Yasuda model, as presented in [47]. Since the investigated aneurysms were located in arteries, the arterial blood flow with the assumed volume flux from [43] seemed to be dominated by high shear rates, thus resulting in a negligible change in fluid parameters, due to non-Newtonian effects, and no considerable change in structural mechanical parameters.
Since neither fluid nor solid mechanical properties were assumed to be viscoelastic, no time delay was introduced into the simulations. Due to the almost incompressibility of both the fluid as well as the vessel wall material [25,46], no significant differences can be seen when comparing one single cycle to five cycles.
Thus, all simulations were run with the base simulation settings of Δt = 0.01 s, Δx = 0.2 mm, 3 solid cells, constant viscosity, and one single cardiac cycle.

3. Results

3.1. Hemodynamic Scenarios in the Parent Artery and Cerebral Aneurysm

Three different scenarios can be identified (see first and third columns in Figure 2 for behavior without the aneurysm sac).
  • In Scenario A, a large area with low WSS (<~0.1 Pa) and an increased OSI (>~0.1) was clearly visible.
  • In Scenario B, high WSS values (i.e., large area with WSS > ~0.5 Pa) were observed, while OSI values remained low (<~0.05).
  • In Scenario C, an intermediate behavior was observed, where certain parts of the region of aneurysm growth clearly showed low WSS, while other regions displayed high WSS values.
We did not find low WSS/low OSI or high WSS/high OSI scenarios in this study.
Figure 2. Hemodynamic scenarios. First two columns show the wall shear stress, the first column without the aneurysm sac and the second with the aneurysm sac; third and fourth columns show the oscillatory shear index OSI, the third column without the aneurysm sac and the fourth column with the aneurysm sac.
Figure 2. Hemodynamic scenarios. First two columns show the wall shear stress, the first column without the aneurysm sac and the second with the aneurysm sac; third and fourth columns show the oscillatory shear index OSI, the third column without the aneurysm sac and the fourth column with the aneurysm sac.
Brainsci 14 00977 g002
The results indicated that WSS and the OSI were reciprocally regulated. Table 2 shows the mean values of the three scenarios without an aneurysm sac. All three scenarios exhibited a significant difference (p < 0.05 in U-test), with scenario A clearly separating itself from the other two. A smaller yet still significant difference was found between scenarios B and C, with scenario C being in-between the two.
The change in Young’s modulus in the structural mechanical part of the simulations from 2490 to 5700 kPa did not affect the hemodynamic part considerably in the geometries without an aneurysm sac. Lowering the value to 500 kPa slightly modified it. Statistical differences between the scenarios remained significant.
In the geometries with aneurysm sacs (second and fourth column in Figure 2), a similar, reciprocally regulated behavior of WSS and the OSI was observed. Additionally, the WSS was lowered, compared to the geometries without an aneurysm sac, approximately by a factor of 5–15, due to the growth of the sac, with an increase in the OSI by a factor of approximately 2–3 in those regions. In Scenario C, the aneurysm sac could be clearly divided into two separate regions, where either WSS was high and the OSI value was low, or vice versa. Aneurysm growth was directed primarily in the direction of low WSS and high OSI in all cases with an aneurysm sac.
Table 3 shows the mean values of aneurysms belonging to the three scenarios for the geometries with an aneurysm sac. A significant difference between the scenarios was evident (all p-values < 0.05). Also, here, Scenario A showed a strong difference compared to scenario B, with scenario C being in-between the two.
The change in Young’s modulus from 2490 kPa to the higher value of 5700 kPa increased the average WSS in aneurysms and reduced the OSI. Introducing the lower value of 500 kPa lowered WSS and increased the OSI to a certain degree. Statistical differences between the scenarios remained significant.
Table 4 shows the distribution of scenarios in percentages amongst the patients with and without aneurysms. In geometries without aneurysms, scenarios B and C were more prominent with similar percentages, whereas scenario A was less likely. Strong contrast geometries with aneurysms showed a tendency towards scenario C. The percentage of scenario A was significantly increased, whereas scenario B dropped to approximately 11% of all cases. Despite the change in Young’s modulus, all scenarios remained unaltered.
Scenarios A and C did not change due to growth in any of the aneurysms. A very different behavior can be seen in geometries in scenario B. Although this scenario can be maintained in 33.42% of the cases, it can change to scenario A (46.65%), as well as to scenario C (19.93%). (Figure 3 shows illustrative cases).
Figure 3. Change in hemodynamic features for scenario B. Individual rows show individual aneurysm geometries (row 1—scenario B to scenario A, row 2—scenario B remains scenario B, row 3—scenario B to scenario C); first two columns show the wall shear stress without the aneurysm sac (first column) and with the aneurysm sac (second column); third and fourth columns show the oscillatory shear index OSI without the aneurysm sac (third column) and with the aneurysm sac (fourth column).
Figure 3. Change in hemodynamic features for scenario B. Individual rows show individual aneurysm geometries (row 1—scenario B to scenario A, row 2—scenario B remains scenario B, row 3—scenario B to scenario C); first two columns show the wall shear stress without the aneurysm sac (first column) and with the aneurysm sac (second column); third and fourth columns show the oscillatory shear index OSI without the aneurysm sac (third column) and with the aneurysm sac (fourth column).
Brainsci 14 00977 g003

3.2. Structural Mechanical Scenarios in the Parent Artery and Cerebral Aneurysm

Two scenarios (D and E) can be identified in the structural mechanical results; however, these scenarios did not fully correspond to the hemodynamic ones (A, B, and C). Thus, structural mechanics and hemodynamics were not directly correlated.
  • Scenario D was characterized by regions with elevated equivalent wall stresses (MISES) (i.e., areas > ~70 kPa visible; see Figure 4, first row).
  • Scenario E indicated regions with low equivalent wall stresses (i.e., areas mostly < ~50 kPa) in both the region of aneurysm onset and the aneurysm wall (see Figure 4, second row).
The aneurysm in the third row of Figure 4 showed a mixed behavior, where the parent vessel without a sac was characterized mostly by scenario E and the configuration with a sac by scenario D.
Figure 4. Structural dynamic features. The first two columns show the equivalent stress, first column without the aneurysm sac and the second column with the aneurysm sac; the third and fourth columns show the equivalent strain, the third column without the aneurysm sac and the fourth column with the aneurysm sac.
Figure 4. Structural dynamic features. The first two columns show the equivalent stress, first column without the aneurysm sac and the second column with the aneurysm sac; the third and fourth columns show the equivalent strain, the third column without the aneurysm sac and the fourth column with the aneurysm sac.
Brainsci 14 00977 g004
Equivalent strain values were relatively low (<0.1 [23]) in all geometries (except the simulation with low Young’s modulus; see model limitations for more details). Table 5 and Table 6 show the mean values of geometries without and with an aneurysm sac, respectively. The statistical difference in the conducted U-test was significant. There was increased stress, as well as strain, due to aneurysm growth.
The change in Young’s modulus (see Table 5 and Table 6) did not change the mean stress considerably, as due to the workflow in FSI, the force of the blood exerted by the blood pressure upon the vessel wall remained constant; however, the strain changed according to the linear material low in a mostly linear fashion.
The results suggested a certain dynamic between the two scenarios. In geometries without aneurysms, scenario E with lowered stress values was more prominent (see Table 7), whereas the initiation of an aneurysm clearly showed an increase in stress towards scenario D. Some Scenario E aneurysms (11.4%) remained in this scenario; however, most changed to scenario D. None of the geometries showed a decrease in stress after aneurysm growth (i.e., no change from scenario D towards E).

4. Discussion

4.1. Hemodynamics

The role of WSS in the aneurysm initiation process has recently been discussed, with divergent results [9]. Kulcsar and Meng suggested that the combination of high WSS and high positive spatial WSS gradient (WSSG) might play a role in aneurysm formation [9,12]. According to Fujimura, the aneurysmal initiation area may correspond to the highest WSS divergence (WSSD point) [7]. Zimny et al. analyzed hemodynamic changes in patients with MCA aneurysms and postulated that a higher WSS may impact aneurysm formation, whereas WSSG may promote this process [8].
In this study, using the distribution of WSS, as well as the OSI, the main hemodynamic processes of aneurysm initiation could be identified. Scenario B resembled patterns identified in previously published studies [9]. Scenario A clearly showed low WSS values in the RAO, along with high OSI values. Interestingly, in this work, additional aneurysms with a complex hemodynamic behavior were presented in-between, as illustrated in Scenario C; certain parts of the region of aneurysm growth clearly showed low WSS values, and other regions showed high WSS values. However, we did not observe a specific low WSS/low OSI or high WSS/high OSI scenario.
Aneurysm initiation has been considered to be partially dependent on chronic inflammation in the vascular wall [7]. Our results may indicate that these processes might be triggered by high WSS conditions, as well as high OSI conditions. However, WSS and the OSI seem to be reciprocally regulated in this stage of the disease. Similarly, for aneurysm growth, Meng et al. postulated the “high-versus-low wall shear stress” controversy, indicating the complexity of this disease. Our results seem to suggest that one of these alternatives in one aneurysm does not rule the other out in another one. While low WSS may be associated with an inflammatory cell-mediated pathway, high WSS may lead to a mural cell-mediated pathology. Both pathways may trigger aneurysm growth and rupture [9]. However, our results suggest that these two pathways may play a role in aneurysm initiation and can coexist.
The investigation revealed a significant dynamic between identified scenarios. WSS did not increase in any of the investigated geometries but decreased after aneurysm initiation. The only change occurred in scenario B, where WSS levels were either maintained (scenario BB in Figure 3), lowered (scenario BA in Figure 3), or partially lowered (scenario BC in Figure 3). Again, the OSI was regulated in an inverse way. Aneurysms with low WSS and a high OSI (scenario A) did not change this behavior due to aneurysm initiation, and a mixed distribution of WSS and OSI (scenario C) was maintained after the start of growth. Our results provide evidence that the inflammatory cell-mediated pathway with low WSS and a high OSI might be the more dominant one during aneurysm initiation.
The analysis of vessel wall stiffness (three different Young’s moduli) strengthened our analysis. Changing the value from 2490 kPa [45] to the maximum value of 5700 kPa in [23] did not influence the hemodynamic behavior. Lowering the value to the lowest reported stiffness in [23] resulted in small changes; however, the identified behavior and dynamics of scenarios did not change.

4.2. Structural Mechanics

The structural mechanical behavior of the vessel and aneurysm wall appeared to be independent of the hemodynamics, showing two different scenarios (D and E). Although scenario D showed elevated stress values (>70–100 kPa), the stress level was still far away from the critical value of 421 kPa, as determined in [38], where rupture could be estimated. Similarly, in the other scenario with even lower stress values, rupture may not be expected.
Aneurysm initiation induced stress in the vessel wall in more than half of the investigated geometries. With the addition of cases with already elevated stress from the start, aneurysm initiation and growth clearly placed additional stress onto the vessel tissue rather than relieving it.
The ES in the geometries without and with an aneurysm sac were below 0.1. Cho et al. [23] provided a critical ES threshold value of 0.3, which was not reached in any of the geometries (with the exception of the simulations at a low Young’s modulus and with an aneurysm sac; see model limitations). In hemodynamics, while the decreased pathway of wall shear stress seemed to be the more dominant one, in the vessel wall, stresses and strains seemed to increase due to aneurysm growth. Since both stress and strain values were significantly below the critical thresholds, as identified in [23,38], respectively, increasing the structural mechanical quantities could be considered a compromise in favor of the preferable reduction in WSS and increase in the OSI in hemodynamics.
The influence of stiffness can be clearly seen in the change in the equivalent strain. By changing the Young’s modulus, the strain changed accordingly in an almost linear fashion. Stress values did not change, as they were defined by the FSI workflow via the force exerted by the blood pressure onto the vessel wall, which remained constant. The small change in strain for the geometries without an aneurysm sac resulted in a negligible or very small change in hemodynamics (see Table 2 and Table 5). The increased change in strain for the geometries with aneurysms resulted in a corresponding increased modification in hemodynamics (see Table 3 and Table 6). However, these changes could only be seen in the mean values, and the scenarios and their dynamics still remained unaltered.
Both the behaviors in stress and strain underline the fact that none of the aneurysms ruptured in reality. However, our results suggest that aneurysm initiation and growth are mainly associated with a reduction in WSS, along with an increase in stress in the vessel wall.

4.3. Model Limitations

Our investigation had certain limitations due to uncertainties in input data, which had to be approximated with certain assumptions. The assumption was made that during aneurysm development, the main shape of the artery does not change considerably (e.g., angle of outflow arteries, etc.). In the FSI simulations, an assumption regarding fixing the geometry must be made. Here, the blood vessel was fixed at the ends of the inflow vessel, as well as the outflow vessels. This does not correspond to real life, as the blood vessels are connected to other vessels. Additionally, the investigated parts may be in contact with parts of brain tissue as well. This contact was neglected in the simulations.
Parent vessels without aneurysms were created by manually removing the sac. Thus, this manual removal introduced a degree of subjective error. Therefore, extra attention was paid during aneurysm removal to minimize such subjective user error.
For the generation of the investigated scenarios, only unruptured but treated aneurysms were considered to avoid the influence of the change in shape due to rupture. In the next step, the influence of ruptured aneurysms will be investigated.
In this work, a constant vessel wall thickness of 0.3 mm was assumed [22]. This topic is currently under investigation, and further studies on wall thickness are needed in future steps to advance toward a more reliable description of aneurysm behavior.
Simulations with an aneurysm sac with the lowest reported value of Young’s modulus in [23] showed increased strain values of >0.5, where the limits of the utilized linear elastic model were reached. This stiffness was considered to be a worst-case scenario for blood vessels. In future investigations, a non-linear approach will be utilized to better describe the material behavior at high strain values.
The data for the pulsatile blood flow was taken from experimental data found in the literature [36]. In experiments, systematic errors occur, which add additional uncertainty into the simulations. Since there was no information on the material properties for the tissues of the individual investigated aneurysms in this work, an assumption from literature was made for the Young’s modulus of the material.

5. Conclusions

Our simulations suggest that hemodynamics, through the reciprocal regulation of WSS and OSI, are closely associated with aneurysm initiation and growth. Hemodynamics and structural mechanics appear to present separate and independent scenarios, with two hemodynamic scenarios remaining stable during growth, while the third showed a tendency towards change. Structural mechanics indicate increases in stress and strain as a result of aneurysm growth. It is important to note that our analysis was limited to unruptured aneurysms, which presents a limitation in fully understanding the rupture risk. Further, although 44 patients were carefully analyzed, we need to confirm these results in a prospective multi-center investigation. However, we believe our findings provide a solid foundation for future research aimed at predicting aneurysm rupture risk by considering their morphology, hemodynamics, and structural mechanics.

Author Contributions

Conceptualization: J.N. and M.G.; methodology: J.N., M.G., W.F., V.M.M., V.M., S.T., J.M., H.S., J.O., N.S. and P.-R.R.; data collection: M.G., W.F., V.M.M., V.M., S.T., H.S., J.O., N.S. and P.-R.R.; formal analysis and investigation: J.N. and M.G.; writing: J.N. and M.G.; funding acquisition: J.N., M.G., Z.M. and A.G.; resources: J.N., M.G., W.F., V.M.M., V.M., S.T., J.M., H.S., J.O., N.S. and P.-R.R.; supervision: J.N., M.G., Z.M. and A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by research subsidies granted by the government of Upper Austria via the Forschungsförderungsgesellschaft (Austrian Research Promotion Agency) [grant number 872604 (Project MEDUSA) and FO999895610 (Project ARES)]. RISC Software GmbH is a member of UAR (Upper Austrian Research) Innovation Network.

Institutional Review Board Statement

This study was approved by the local ethics committee (Ethikkommission der medizinischen Fakultät der Johannes Kepler Universität; EK Nr: 1129/2022), and the requirement for the acquisition of informed consent from patients was waived, owing to the retrospective nature of the research.

Informed Consent Statement

The requirement for the acquisition of informed consent from patients was waived, owing to the retrospective nature of the research.

Data Availability Statement

The data are not publicly available due to the fact that authors want to keep patient data private, however data can be made available upon request.

Conflicts of Interest

Author Jozsef Nagy was employed by the company eulerian-solutions e.U., Authors Wolfgang Fenz, Stefan Thumfart were employed by the company RISC Software GmbH. All authors certify that they have no affiliations with or involvement in any organization or entity with any financial interest (such as honoraria; educational grants; participation in speakers’ bureaus; membership, employment, consultancies, stock ownership, or other equity interest; and expert testimony or patent-licensing arrangements) or non-financial interest (such as personal or professional relationships, affiliations, knowledge, or beliefs) in the subject matter or materials discussed in this manuscript. Generative AI was not used in the process of scientific writing.

References

  1. Dhar, S.; Tremmel, M.; Mocco, J.; Kim, M.; Yamamoto, J.; Siddiqui, A.H.; Hopkins, L.N.; Meng, H. Morphology parameters for intracranial aneurysm rupture risk assessment. Neurosurgery 2008, 63, 185–197. [Google Scholar] [CrossRef] [PubMed]
  2. Hackenberg, K.; Hänggi, D.; Etminan, N. Unruptured intracranial aneurysms contemporary data and management. Stroke 2018, 49, 2268–2275. [Google Scholar] [CrossRef] [PubMed]
  3. Jirjees, S.; Htun, Z.M.; Aldawudi, I.; Katwal, P.C.; Khan, S. Role of Morphological and Hemodynamic Factors in Predicting Intracranial Aneurysm Rupture: A Review. Cureus 2020, 12, e9178. [Google Scholar] [CrossRef] [PubMed]
  4. Detmer, F.J.; Chung, B.J.; Mut, F.; Pritz, M.; Slawski, M.; Hamzei-Sichani, F.; Kallmes, D.; Putman, C.; Jimenez, C.; Cebral, J.R. Development of a statistical model for discrimination of rupture status in posterior communicating artery aneurysms. Acta Neurochir. 2018, 160, 1643–1652. [Google Scholar] [CrossRef]
  5. Jiang, P.; Liu, Q.; Wu, J.; Chen, X.; Li, M.; Li, Z.; Yang, S.; Guo, R.; Gao, B.; Cao, Y.; et al. A Novel Scoring System for Rupture Risk Stratification of Intracranial Aneurysms: A Hemodynamic and Morphological Study. Front Neurosci. 2018, 12, 596. [Google Scholar] [CrossRef]
  6. Soldozy, S.; Norat, P.; Elsarrag, M.; Chatrath, A.; Costello, J.S.; Sokolowski, J.D.; Tvrdik, P.; Kalani, M.Y.S.; Park, M.S. The biophysical role of hemodynamics in the pathogenesis of cerebral aneurysm formation and rupture. Neurosurg. Focus. 2019, 47, E11. [Google Scholar] [CrossRef]
  7. Fujimura, S.; Tanaka, K.; Takao, H.; Okudaira, T.; Koseki, H.; Hasebe, A.; Suzuki, T.; Uchiyama, Y.; Ishibashi, T.; Otani, K.; et al. Computational fluid dynamic analysis of the initiation of cerebral aneurysms. J. Neurosurg. 2022, 137, 335–343. [Google Scholar] [CrossRef]
  8. Zimny, M.; Kawlewska, E.; Hebda, A.; Wolański, W.; Ładziński, P.; Kaspera, W. Wall shear stress gradient is independently associated with middle cerebral artery aneurysm development: A case-control CFD patient-specific study based on 77 patients. BMC Neurol. 2021, 21, 281. [Google Scholar] [CrossRef] [PubMed]
  9. Meng, H.; Tutino, V.M.; Xiang, J.; Siddiqui, A. High WSS or Low WSS? Complex Interactions of Hemodynamics with Intracranial Aneurysm Initiation, Growth, and Rupture: Toward a Unifying Hypothesis. AJNR Am. J. Neuroradiol. 2014, 35, 1254–1262. [Google Scholar] [CrossRef]
  10. Rajabzadeh-Oghaz, H.; Siddiqui, A.H.; Asadollahi, A.; Kolega, J.; Tutino, V.M. The association between hemodynamics and wall characteristics in human intracranial aneurysms: A review. Neurosurg. Rev. 2022, 45, 49–61. [Google Scholar] [CrossRef]
  11. Villablanca, J.P.; Duckwiler, G.R.; Jahan, R.; Tateshima, S.; Martin, N.A.; Frazee, J.; Gonzalez, N.R.; Sayre, J.; Vinuela, F.V. Natural history of asymptomatic unruptured cerebral aneurysms evaluated at CT angiography: Growth and rupture incidence and correlation with epidemiologic risk factors. Neuroradiology 2013, 269, 258–265. [Google Scholar] [CrossRef] [PubMed]
  12. Kulcsár, Z.; Ugron, A.; Marosfői, M.; Berentei, Z.; Paál, G.; Szikora, I. Hemodynamics of cerebral aneurysm initiation: The role of wall shear stress and spatial wall shear stress gradient. AJNR Am. J. Neuroradiol. 2011, 32, 587–594. [Google Scholar] [CrossRef]
  13. Xiang, J.; Natarajan, S.K.; Tremmel, M.; Ma, D.; Mocco, J.; Hopkins, L.N.; Siddiqui, A.H.; Levy, E.I.; Meng, H. Hemodynamic-morphologic discriminants for intracranial aneurysm rupture. Stroke 2011, 42, 144–152. [Google Scholar] [CrossRef] [PubMed]
  14. Jou, L.D.; Lee, D.H.; Morsi, H.; Mawad, M. Wall shear stress on ruptured and unruptured intracranial aneurysms at the internal carotid artery. AJNR Am. J. Neuroradiol. 2008, 29, 1761–1767. [Google Scholar] [CrossRef] [PubMed]
  15. Shojima, M.; Oshima, M.; Takagi, K.; Torii, R.; Hayakawa, M.; Katada, K.; Morita, A.; Kirino, T. Magnitude and role of wall shear stress on cerebral aneurysm: Computational fluid dynamic study of 20 middle cerebral artery aneurysms. Stroke 2004, 35, 2500–2505. [Google Scholar] [CrossRef]
  16. Boussel, L.; Rayz, V.; McCulloch, C.; Martin, A.; Acevedo-Bolton, G.; Lawton, M.; Higashida, R.; Smith, W.S.; Young, W.L.; Saloner, D. Aneurysm growth occurs at region of low wall shear stress: Patient-specific correlation of hemodynamics and growth in a longitudinal study. Stroke 2008, 39, 2997–3002. [Google Scholar] [CrossRef]
  17. Acevedo-Bolton, G.; Jou, L.-D.; Dispensa, B.P.; Lawton, M.T.; Higashida, R.T.; Martin, A.J.; Young, W.L.; Saloner, D. Estimating the hemodynamic impact of interventional treatments of aneurysms: Numerical simulation with experimental validation: Technical case report. Neurosurgery 2006, 59, E429–E430. [Google Scholar] [CrossRef] [PubMed]
  18. Castro, M.A.; Putman, C.M.; Sheridan, M.J.; Cebral, J.R. Hemodynamic patterns of anterior communicating artery aneurysms: A possible association with rupture. AJNR Am. J. Neuroradiol. 2009, 30, 297–302. [Google Scholar] [CrossRef]
  19. Qian, Y.; Takao, H.; Umezu, M.; Murayama, Y. Risk analysis of unruptured aneurysms using computational fluid dynamics technology: Preliminary results. AJNR Am. J. Neuroradiol. 2011, 32, 1948–1955. [Google Scholar] [CrossRef]
  20. Sugiyama, S.-I.; Meng, H.; Funamoto, K.; Inoue, T.; Fujimura, M.; Nakayama, T.; Omodaka, S.; Shimizu, H.; Takahashi, A.; Tominaga, T. Hemodynamic analysis of growing intracranial aneurysms arising from a posterior inferior cerebellar artery. World Neurosurg. 2012, 78, 462–468. [Google Scholar] [CrossRef]
  21. Takao, H.; Murayama, Y.; Otsuka, S.; Qian, Y.; Mohamed, A.; Masuda, S.; Yamamoto, M.; Abe, T. Hemodynamic differences between unruptured and ruptured intracranial aneurysms during observation. Stroke 2012, 43, 1436–1439. [Google Scholar] [CrossRef] [PubMed]
  22. Acosta, J.M.; Cayron, A.F.; Dupuy, N.; Pelli, G.; Foglia, B.; Haemmerli, J.; Allémann, E.; Bijlenga, P.; Kwak, B.R.; Morel, S. Effect of Aneurysm and Patient Characteristics on Intracranial Aneurysm Wall Thickness. Front. Cardiovasc. Med. 2021, 8, 775307. [Google Scholar] [CrossRef] [PubMed]
  23. Cho, C.K.; Yang, H.; Kim, J.J.; Oh, J.H.; Kim, Y.B. Prediction of rupture risk in cerebral aneurysms by comparing clinical cases with fluid–structure interaction analyses. Sci. Rep. 2020, 10, 18237. [Google Scholar] [CrossRef] [PubMed]
  24. Gasser, T.C.; Auer, M.; Labruto, F.; Swedenborg, J.; Roy, J. Biomechanical rupture risk assessment of abdominal aortic aneurysms: Model complexity versus predictability of finite element simulations. Eur. J. Vasc. Endovasc. Surg. 2010, 40, 176–185. [Google Scholar] [CrossRef] [PubMed]
  25. Reeps, C.; Gee, M.; Maier, A.; Gurdan, M.; Eckstein, H.H.; Wall, W.A. The impact of model assumptions on results of computational mechanics in abdominal aortic aneurysm. J. Vasc. Surg. 2010, 51, 679–688. [Google Scholar] [CrossRef]
  26. Rissland, P.; Alemu, Y.; Einav, S.; Ricotta, J.; Bluestein, D. Abdominal aortic aneurysm risk of rupture: Patient-specific FSI simulations using anisotropic model. J. Biomech. Eng. 2009, 131, 031001. [Google Scholar] [CrossRef]
  27. Huang, Z.; Zeng, M.; Tao, W.G.; Zeng, F.Y.; Chen, C.Q.; Zhang, L.B.; Chen, F.H. A Hemodynamic Mechanism Correlating with the Initiation of MCA Bifurcation Aneurysms. Am. J. Neuroradiol. 2020, 41, 1217–1224. [Google Scholar] [CrossRef]
  28. Bergman, T.L.; Lavine, A.S.; Incropera, F.P.; De Witt, D.P. Fundamentals of Heat and Mass Transfer; John Wiley & Sons: Hoboken, NJ, USA, 2018; ISBN 978-1-119-35388-1. [Google Scholar]
  29. Çengel, Y.A.; Cimbala, J.M. Fluid Mechanics: Fundamentals and Applications, 4th ed.; McGraw-Hill Higher Education: Boston, MA, USA, 2019; ISBN 978-9-813-15788-0. [Google Scholar]
  30. Guo, H.; Liu, J.F.; Li, C.H.; Wang, J.W.; Li, H.; Gao, B.L. Greater hemodynamic stresses initiate aneurysms on major cerebral arterial bifurcations. Front. Neurol. 2023, 14, 1265484. [Google Scholar] [CrossRef]
  31. Guo, H.; Yang, S.T.; Wang, J.W.; Li, H.; Gao, B.L.; Li, C.H. High hemodynamic stresses induce aneurysms at internal carotid artery bends. Medicine 2023, 102, e34587. [Google Scholar] [CrossRef]
  32. Geers, A.J.; Morales, H.G.; Larrabide, I.; Butakoff, C.; Bijlenga, P.; Frangi, A.F. Wall shear stress at the initiation site of cerebral aneurysms. Biomech. Model. Mechanobiol. 2017, 16, 97–115. [Google Scholar] [CrossRef]
  33. Mantha, A.R.; Karmonik, C.; Benndorf, G.; Strother, C.; Metcalfe, R.W. Hemodynamics in a cerebral artery before and after the formation of an aneurysm. Am. J. Neuroradiol. 2006, 27, 1113–1118. [Google Scholar] [PubMed]
  34. Baek, H.; Jayaraman, M.V.; Karniadakis, G.E. Wall shear stress and pressure distribution on aneurysms and infundibulae in the posterior communicating artery bifurcation. Ann. Biomed. Eng. 2009, 37, 2469–2487. [Google Scholar] [CrossRef] [PubMed]
  35. Castro, M.A.; Putman, C.M.; Cebral, J.R. Computational analysis of anterior communicating artery aneurysm shear stress before and after aneurysm formation. J. Phys. Conf. Ser. 2011, 332, 012001. [Google Scholar] [CrossRef]
  36. Lauric, A.; Hippelheuser, J.; Safain, M.G.; Malek, A.M. Curvature effect on hemodynamic conditions at the inner bend of the carotid siphon and its relation to aneurysm formation. J. Biomech. 2014, 47, 3018–3027. [Google Scholar] [CrossRef]
  37. Shimogonya, Y.; Ishikawa, T.; Imai, Y.; Matsuki, N.; Yamaguchi, T. Can temporal fluctuation in spatial wall shear stress gradient initiate a cerebral aneurysm? A proposed novel hemodynamic index, the gradient oscillatory number (GON). J. Biomech. 2009, 42, 550–554. [Google Scholar] [CrossRef]
  38. Cebral, J.R.; Sheridan, M.; Putman, C.M. Hemodynamics and bleb formation in intracranial aneurysms. AJNR Am. J. Neuroradiol. 2010, 31, 304–310. [Google Scholar] [CrossRef]
  39. Nagy, J.; Maier, J.; Miron, V.; Fenz, W.; Major, Z.; Gruber, A.; Gmeiner, M. Methods, Validation and Clinical Implementation of a Simulation Method of Cerebral Aneurysms. J. Biomed. Eng. Biosci. 2023, 10, 10–19. [Google Scholar] [CrossRef]
  40. Weller, H.G.; Tabor, G.; Jasak, H.; Fureby, C. A tensorial approach to computational continuum mechanics using object orientated techniques. Comput. Phys. 1998, 12, 620–631. [Google Scholar] [CrossRef]
  41. Greenshields, C.J.; Weller, H.G. Notes on Computational Fluid Dynamics: General Principles; CFD Direct: Reading, UK, 2022; ISBN 978-1-3999-2078-0. [Google Scholar]
  42. Cardiff, P.; Karač, A.; Jaeger, P.D.; Jasak, H.; Nagy, J.; Ivanković, A.; Tuković, Ž. An open-source finite volume toolbox for solid mechanics and fluid-solid interaction simulations. arXiv 2018, arXiv:1808.10736v2. Available online: https://arxiv.org/abs/1808.10736 (accessed on 20 September 2024).
  43. Blanco, P.; Müller, L.; Spence, J.D. Blood pressure gradients in cerebral arteries: A clue to pathogenesis of cerebral small vessel disease. Stroke Vasc. Neurol. 2017, 2, 108–117. [Google Scholar] [CrossRef]
  44. Lorenzetti, F.; Suominen, S.; Tukianen, E. Evaluation of Blood Flow in Free Microvascular Flaps. J. Reconstr. Microsurg. 2001, 17, 163–167. [Google Scholar] [CrossRef] [PubMed]
  45. Tóth, B.K. The Mechanical Interaction between the Red Blood Cells and the Blood Vessels. Ph.D. Thesis, University of Budapest, Budapest, Hungary, 2011. [Google Scholar]
  46. Di Martino, E.S.; Guadagni, G.; Fumero, A.; Ballerini, G.; Spirito, R.; Biglioli, P.; Redaelli, A. Fluid-structure interaction within realistic three-dimensional models of the aneurysmatic aorta as a guidance to assess the risk of rupture of the aneurysm. Med. Eng. Phys. 2001, 23, 647–655. [Google Scholar] [CrossRef] [PubMed]
  47. Lynch, S.; Nama, N.; Figueroa, C.A. Effects of non-Newtonian viscosity on arterial and venous flow and transport. Sci. Rep. 2022, 12, 20568. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Typical CAD geometry of a cerebral artery, with aneurysm sac (left), aneurysm sac removed (center), and indicator line of the region of aneurysm onset (RAO, right).
Figure 1. Typical CAD geometry of a cerebral artery, with aneurysm sac (left), aneurysm sac removed (center), and indicator line of the region of aneurysm onset (RAO, right).
Brainsci 14 00977 g001
Table 1. Convergence analysis of time step size and mesh size independence, as well as the influence of the non-Newtonian viscosity model and multiple cardiac cycles. (Base simulation settings consisted of Δt = 0.01 s, Δx = 0.2 mm, 3 solid cells, Newtonian model, and one single cardiac cycle).
Table 1. Convergence analysis of time step size and mesh size independence, as well as the influence of the non-Newtonian viscosity model and multiple cardiac cycles. (Base simulation settings consisted of Δt = 0.01 s, Δx = 0.2 mm, 3 solid cells, Newtonian model, and one single cardiac cycle).
ParameterWSS [Pa]OSI [-]ES [-]MISES [kPa]Runtime [s]
Time independenceΔt = 0.1 s0.0830.1480.03373.54347
Δt = 0.01 s0.0850.1500.03475.41715
Δt = 0.001 s0.0860.1510.03575.431872
Fluid mesh convergenceΔx = 0.3 mm0.0800.1390.03169.84262
Δx = 0.25 mm0.0830.1450.03372.14426
Δx = 0.2 mm0.0860.1500.03475.41715
Δx = 0.15 mm0.0880.1520.03576.111592
Solid mesh convergence2 cells0.0860.1500.03273.94578
3 cells0.0860.1500.03475.41715
4 cells0.0860.1510.03475.421643
Viscosity modelNewtonian0.0850.1500.03475.41715
Non-Newtonian0.0860.1510.03475.41727
Number of cardiac cycles10.0850.1500.03475.41715
50.0850.1500.03475.413512
Table 2. Mean values of scenarios A, B, and C, as well as p-values between scenarios without aneurysms.
Table 2. Mean values of scenarios A, B, and C, as well as p-values between scenarios without aneurysms.
Young’s Modulus Mean Valuesp-Values
ABCA vs. BA vs. CB vs. C
2490 kPaWSS [Pa]0.2611.7721.1761.90 × 10−48.86 × 10−43.82 × 10−2
OSI [-]0.0750.0180.037.52 × 10−42.88 × 10−34.22 × 10−2
5700 kPaWSS [Pa]0.2591.7681.1641.86 × 10−48.74 × 10−43.73 × 10−2
OSI [-]0.0770.0200.0327.59 × 10−42.92 × 10−34.29 × 10−2
500 kPaWSS [Pa]0.2311.5681.0401.88 × 10−48.81 × 10−43.83 × 10−2
OSI [-]0.0840.0200.0347.56 × 10−42.89 × 10−34.25 × 10−2
Table 3. Mean values of scenarios A, B, and C, as well as p-values between scenarios with aneurysms.
Table 3. Mean values of scenarios A, B, and C, as well as p-values between scenarios with aneurysms.
Young’s Modulus Mean Valuesp-Values
ABCA vs. BA vs. CB vs. C
2490 kPaWSS [Pa]0.0250.390.1933.00 × 10−38.54 × 10−41.83 × 10−2
OSI [-]0.1540.0370.0561.23 × 10−23.72 × 10−33.22 × 10−2
5700 kPaWSS [Pa]0.0340.4410.2431.98 × 10−36.37 × 10−41.13 × 10−2
OSI [-]0.1130.0330.0449.54 × 10−32.27 × 10−31.92 × 10−2
500 kPaWSS [Pa]0.0160.2860.1372.87 × 10−37.23 × 10−41.95 × 10−2
OSI [-]0.2410.0500.0781.33 × 10−23.89 × 10−33.12 ×10−2
Table 4. Percentage of hemodynamic scenarios with and without aneurysms.
Table 4. Percentage of hemodynamic scenarios with and without aneurysms.
ABC
without aneurysm22.7%34.1%43.2%
with aneurysm38.6%11.4%50.0%
Table 5. Mean values of scenarios D and E, as well as p-values between scenarios without aneurysms.
Table 5. Mean values of scenarios D and E, as well as p-values between scenarios without aneurysms.
Young’s Modulus Mean Valuesp-Value
DED vs. E
2490 kPaMISES [kPa]76.545.81.27 × 10−3
ES [-]0.0330.0194.14 × 10−3
5700 kPaMISES [kPa]77.346.31.65 × 10−3
ES [-]0.0140.0084.89 × 10−3
500 kPaMISES [kPa]77.046.21.33 × 10−3
ES [-]0.1670.0954.43 × 10−3
Table 6. Mean values of scenarios D and E, as well as p-values between scenarios with aneurysms.
Table 6. Mean values of scenarios D and E, as well as p-values between scenarios with aneurysms.
Young’s Modulus Mean Valuesp-Value
DED vs. E
2490 kPaMISES [kPa]83.649.01.55 × 10−4
ES [-]0.1030.0461.38 × 10−3
5700 kPaMISES [kPa]84.249.32.67 × 10−4
ES [-]0.0430.0193.39 × 10−3
500 kPaMISES [kPa]83.949.11.94 × 10−4
ES [-]0.5150.231.77 × 10−3
Table 7. Percentage of structural mechanic scenarios with and without aneurysms.
Table 7. Percentage of structural mechanic scenarios with and without aneurysms.
DE
without aneurysm43.2%56.8%
with aneurysm88.6%11.4%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nagy, J.; Fenz, W.; Miron, V.M.; Thumfart, S.; Maier, J.; Major, Z.; Stefanits, H.; Oberndorfer, J.; Stroh, N.; Mazanec, V.; et al. Fluid–Structure Interaction Simulations of the Initiation Process of Cerebral Aneurysms. Brain Sci. 2024, 14, 977. https://doi.org/10.3390/brainsci14100977

AMA Style

Nagy J, Fenz W, Miron VM, Thumfart S, Maier J, Major Z, Stefanits H, Oberndorfer J, Stroh N, Mazanec V, et al. Fluid–Structure Interaction Simulations of the Initiation Process of Cerebral Aneurysms. Brain Sciences. 2024; 14(10):977. https://doi.org/10.3390/brainsci14100977

Chicago/Turabian Style

Nagy, Jozsef, Wolfgang Fenz, Veronika M. Miron, Stefan Thumfart, Julia Maier, Zoltan Major, Harald Stefanits, Johannes Oberndorfer, Nico Stroh, Vanessa Mazanec, and et al. 2024. "Fluid–Structure Interaction Simulations of the Initiation Process of Cerebral Aneurysms" Brain Sciences 14, no. 10: 977. https://doi.org/10.3390/brainsci14100977

APA Style

Nagy, J., Fenz, W., Miron, V. M., Thumfart, S., Maier, J., Major, Z., Stefanits, H., Oberndorfer, J., Stroh, N., Mazanec, V., Rauch, P. -R., Gruber, A., & Gmeiner, M. (2024). Fluid–Structure Interaction Simulations of the Initiation Process of Cerebral Aneurysms. Brain Sciences, 14(10), 977. https://doi.org/10.3390/brainsci14100977

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop