Numerical simulation of an acoustically excited plane jet in an incompressible framework and comparison with experimental data.

This paper is concerned with the interaction of a two-dimensional plane jet with transverse plane acoustic waves, which occur, for example, in flue instruments in the vicinity of the nozzle exit. The acoustic excitation is modeled with fluctuating boundary conditions within the framework of an incompressible simulation. This method can be easily implemented in commercial CFD software. The simulations are compared with well-documented measurement data from other authors. It is shown that the model can predict not only the growth rate but also the amplitude of the perturbation, thus it captures the key features of the jet behavior in the receptive and in the linear growth region.

Emerging fractal patterns in a real 3D cerebral aneurysm.

The behaviour of biological fluid flows is often investigated in medical practice to draw conclusions on the physiological or pathological conditions of the considered organs. One area where such investigations are proven to be useful is the flow-related formation and growth of different pathologic malformations of the cerebro-vascular system. In this work, a detailed study is presented on the effect of a cerebral aneurysm on blood transport inside a human brain artery segment. This malformation causes strong flow instabilities that drives the flow system towards chaotic behaviour. The emerging fractal structure and some of its measurable properties have been explored using a method that makes the measurement of these properties feasible even in complicated large three dimensional data sets. We find that, from the investigated chaos parameters, the information dimension turns out to be the most reliable parameter to characterise chaotic advection in the vicinity of the aneurysm sac. We propose that properties of chaotic mixing close to aneurysms might be relevant for the condition of this pathologic malformation.

Central arterial pressure estimation based on two peripheral pressure measurements using one-dimensional blood flow simulation.

Aortic pressure can be estimated using one-dimensional arterial flow simulations. This study demonstrates that two peripheral pressure measurements can be used to acquire the central pressure curve through the patient-specific optimization of a set of system parameters. Radial and carotid pressure measurements and parameter optimization were performed in the case of 62 patients. The two calculated aortic curves were in good agreement, Systolic and Mean Blood Pressures differed on average by 0.5 and -0.5 mmHg, respectively. Good agreement was achieved with the transfer function method as well. The effect of carotid clamping is demonstrated using one resulting patient-specific arterial network.

Comparison of Flow Reduction Efficacy of Nominal and Oversized Flow Diverters Using a Novel Measurement-assisted in Silico Method.

PURPOSE: The high efficacy of flow diverters (FD) in the case of wide-neck aneurysms is well demonstrated, yet new challenges have arisen because of reported posttreatment failures and the growing number of new generation of devices. Our aim is to present a measurement-supported in silico workflow that automates the virtual deployment and subsequent hemodynamic analysis of FDs. In this work, the objective is to analyze the effects of FD deployment variability of two manufacturers on posttreatment flow reduction. METHODS: The virtual deployment procedure is based on detailed mechanical calibration of the flow diverters, while the flow representation is based on hydrodynamic resistance (HR) measurements. Computational fluid dynamic simulations resulted in 5 untreated and 80 virtually treated scenarios, including 2 FD designs in nominal and oversized deployment states. The simulated aneurysmal velocity reduction (AMVR) is correlated with the HR values and deployment scenarios. RESULTS: The linear HR coefficient and AMVR revealed a power-law relationship considering all 80 deployments. In nominal deployment scenarios, a significantly larger average AMVR was obtained (60.3%) for the 64-wire FDs than for 48-wire FDs (51.9%). In oversized deployments, the average AMVR was almost the same for 64-wire and 48-wire device types, 27.5% and 25.7%, respectively. CONCLUSION: The applicability of our numerical workflow was demonstrated, also in large-scale hemodynamic investigations. The study revealed a robust power-law relationship between a HR coefficient and AMVR. Furthermore, the 64 wire configurations in nominal sizing produced a significantly higher posttreatment flow reduction, replicating the results of other in vitro studies.

Fractals and Chaos in the Hemodynamics of Intracranial Aneurysms.

Computing the emerging flow in blood vessel sections by means of computational fluid dynamics is an often applied practice in hemodynamics research. One particular area for such investigations is related to the cerebral aneurysms, since their formation, pathogenesis, and the risk of a potential rupture may be flow-related. We present a study on the behavior of small advected particles in cerebral vessel sections in the presence of aneurysmal malformations. These malformations cause strong flow disturbances driving the system toward chaotic behavior. Within these flows, the particle trajectories can form a fractal structure, the properties of which are measurable by quantitative techniques. The measurable quantities are well established chaotic properties, such as the Lyapunov exponent, escape rate, and information dimension. Based on these findings, we propose that chaotic flow within blood vessels in the vicinity of the aneurysm might be relevant for the pathogenesis and development of this malformation.

Flow simulations on an organ pipe foot model.

The present paper shows numerical simulations of the flow responsible for the sound generation in an organ pipe. Only the foot model of the organ pipe (i.e., with the resonator detached) is investigated by two-dimensional incompressible CFD simulations. It is shown that in spite of the moderately high Reynolds number (Re≈2350) no turbulence modeling is necessary. Free jet simulation (foot model without the upper lip) showed that the jet oscillates due to its natural instability. The velocity profile, the centerline and the width of the jet is determined at different heights above the flue. Edge tone simulations (foot model with the upper lip) were carried out having the upper lip at a constant height but at different x positions. It was found that the strongest and most stable edge tone oscillation occurs if the lower left corner of the upper lip is in the centerline of the jet (optimum position). When the upper lip is far from its optimum position the oscillation of the jet is rather due to the natural instability of the jet than the edge tone phenomenon. The results agree well with the experimental results of Außerlechner et al. [J. Acoust. Soc. Am. 126, 878-886 (2009)] and Außerlechner (Ph.D. thesis, Universität Stuttgart, Stuttgart, Germany) and with former results of the authors [Paál and Vaik, Int. J. Heat Fluid Flow 28, 575-586 (2007); Paál and Vaik, in Conference on Modelling Fluid Flow (CMFF'09), Budapest, Hungary].

Design space exploration of flow diverter hydraulic resistance parameters in sidewall intracranial aneurysms.

Intracranial aneurysms are nowadays treated with endovascular flow diverter devices to avoid sac rupture. This study explores how different linear and quadratic hydrodynamic resistance parameters reduce the flow in the sac for five patient-specific sidewall aneurysms.The 125 performed blood flow simulations included the stents using a Darcy-Forcheimer porous layer approach based on real-life stent characteristics. Time- and space-averaged velocity magnitudes were strongly affected by the linear coefficient with a power-law relationship. Quadratic coefficients alter the flow in a minor way due to the low-velocity levels in the aneurysm sac and neck region.

First blood: An efficient, hybrid one- and zero-dimensional, modular hemodynamic solver.

Low-dimensional (1D or 0D) models can describe the whole human blood circulation, for example, 1D distributed parameter model for the arterial network and 0D concentrated models for the heart or other organs. This paper presents a combined 1D-0D solver, called first_blood, that solves the governing equations of fluid dynamics to model low-dimensional hemodynamic effects. An extended method of characteristics is applied here to solve the momentum, and mass conservation equations and the viscoelastic wall model equation, mimicking the material properties of arterial walls. The heart and the peripheral lumped models are solved with a general zero-dimensional (0D) nonlinear solver. The model topology can be modular, that is, first_blood can solve any 1D-0D hemodynamic model. To demonstrate the applicability of first_blood, the human arterial system, the heart and the peripherals are modelled using the solver. The simulation time of a heartbeat takes around 2 s, that is, first_blood requires only twice the real-time for the simulation using an average PC, which highlights the computational efficiency. The source code is available on GitHub, that is, it is open source. The model parameters are based on the literature suggestions and on the validation of output data to obtain physiologically relevant results.

Simplified coronary flow reserve calculations based on three-dimensional coronary reconstruction and intracoronary pressure data.

BACKGROUND: Measurements of fractional flow reserve (FFR) and/or coronary flow reserve (CFR) are widely used for hemodynamic characterization of coronary lesions. The frequent combination of the epicardial and microvascular disease may indicate a need for complex hemodynamic evaluation of coronary lesions. This study aims at validating the calculation of CFR based on a simple hemodynamic model to detailed computational fluid dynamics (CFD) analysis. METHODS: Three-dimensional (3D) morphological data and pressure values from FFR measurements were used to calculate the target vessel. Nine patients with one intermediate stenosis each, measured by pressure wire, were included in this study. RESULTS: A correlation was found between the determined CFR from simple equations and from a steady flow simulation (r = 0.984, p < 10-5). There was a significant correlation between the CFR values calculated by transient and steady flow simulations (r = 0.94, p < 10-3). CONCLUSIONS: Feasibility was demonstrated of a simple hemodynamic calculation of CFR based on 3D-angiography and intracoronary pressure measurements. A simultaneous determination of both the FFR and CFR values provides the capability to diagnose microvascular dysfunction: the CFR/FFR ratio characterizes the microvascular reserve.

Central arterial pressure and patient-specific model parameter estimation based on radial pressure measurements.

One-dimensional arterial flow simulations are suitable to estimate the aortic pressure from peripheral measurements in a patient-specific arterial network. This study introduces a reduction of the system parameters, and a novel calculation method to estimate the patient-specific set and the aortic curve based on radial applanation tonometry. Peripheral and aortic pressure curves were measured in patients, optimization were carried out. The aortic pressure curves were reproduced well, with an overestimation of the measured Systolic and Mean Blood Pressures on average by 0.6 and 0.5 mmHg respectively, and the Root Mean Square Difference of the curves was 3 mmHg on average.

Impact of Design and Deployment Technique on the Hydrodynamic Resistance of Flow Diverters : An in Vitro Experimental Study.

PURPOSE: Despite the high efficacy of flow diverters (FD) in treating sidewall intracranial aneurysms, failures are reported. One of the physical factors determining efficacy is the flow reducing capacity of the FD that is currently unknown to the operator. Our aim was to measure the flow reducing capacity expressed as the hydrodynamic resistance (HR), the metallic surface area (MSA) and pore density (PD) of two different FD designs and quantitatively investigate the impact of sizing and the deployment technique on these parameters. METHODS: Altogether 38 Pipeline (Medtronic) and P64 (Phenox) FD­s were implanted in holder tubes by a neurointerventionist in nominally sized, oversized and longitudinally compressed or elongated manners. The tubes were placed in a flow model with the flow directed across the FD through a side hole on the tube. HR was expressed by the measured pressure drop as the function of the flow rate. Deployed length, MSA and PD were also measured and correlated with the HR. RESULTS: Both PD and MSA changed with varying deployment length, which correlates well with the change in HR. Oversizing the device by 1 mm in diameter has reduced the HR on average to one fifth of the original value for both manufacturers. CONCLUSION: This study demonstrates experimentally that different FD designs have different flow diverting capacities (HR). Parameters are greatly influenced by radial sizing and longitudinal compression or elongation during implantation. Our results might be useful in procedure planning, predicting clinical outcome, and in patient-specific numerical flow simulations.

Pressure- and 3D-Derived Coronary Flow Reserve with Hydrostatic Pressure Correction: Comparison with Intracoronary Doppler Measurements.

Purpose: To develop a method of coronary flow reserve (CFR) calculation derived from three-dimensional (3D) coronary angiographic parameters and intracoronary pressure data during fractional flow reserve (FFR) measurement. Methods: Altogether 19 coronary arteries of 16 native and 3 stented vessels were reconstructed in 3D. The measured distal intracoronary pressures were corrected to the hydrostatic pressure based on the height differences between the levels of the vessel orifice and the sensor position. Classical fluid dynamic equations were applied to calculate the flow during the resting state and vasodilatation based on morphological data and intracoronary pressure values. 3D-derived coronary flow reserve (CFRp-3D) was defined as the ratio between the calculated hyperemic and the resting flow and was compared to the CFR values simultaneously measured by the Doppler sensor (CFRDoppler). Results: Haemodynamic calculations using the distal coronary pressures corrected for hydrostatic pressures showed a strong correlation between the individual CFRp-3D values and the CFRDoppler measurements (r = 0.89, p < 0.0001). Hydrostatic pressure correction increased the specificity of the method from 46.1% to 92.3% for predicting an abnormal CFRDoppler < 2. Conclusions: CFRp-3D calculation with hydrostatic pressure correction during FFR measurement facilitates a comprehensive hemodynamic assessment, supporting the complex evaluation of macro-and microvascular coronary artery disease.

A novel virtual flow diverter implantation method with realistic deployment mechanics and validated force response.

Virtual flow diverter deployment techniques underwent significant development during the last couple of years. Each existing technique displays advantageous features, as well as significant limitations. One common drawback is the lack of quantitative validation of the mechanics of the device. In the following work, we present a new spring-mass-based method with validated mechanical responses that combines many of the useful capabilities of previous techniques. The structure of the virtual braids naturally incorporates the axial length changes as a function of the local expansion diameter. The force response of the model was calibrated using the measured response of real FDs. The mechanics of the model allows to replicate the expansion process during deployment, including additional effects such as the push-pull technique that is required for the deployment of braided FDs to achieve full opening and proper wall apposition. Furthermore, it is a computationally highly efficient solution that requires little pre-processing and has a run-time of a few seconds on a general laptop and thus allows for exploratory analyses. The model was applied in a patient-specific geometry, where corresponding accurate control measurements in a 3D-printed model were also available. The analysis shows the effects of FD oversizing and push-pull application on the radial expansion, surface density, and on the wall contact pressure.

Hydrodynamic Resistance of Intracranial Flow-Diverter Stents: Measurement Description and Data Evaluation.

PURPOSE: Intracranial aneurysms are malformations forming bulges on the walls of brain arteries. A flow diverter device is a fine braided wire structure used for the endovascular treatment of brain aneurysms. This work presents a rig and a protocol for the measurement of the hydrodynamic resistance of flow diverter stents. Hydrodynamic resistance is interpreted here as the pressure loss versus volumetric flow rate function through the mesh structure. The difficulty of the measurement is the very low flow rate range and the extreme sensitivity to contamination and disturbances. METHODS: Rigorous attention was paid to reproducibility, hence a strict protocol was designed to ensure controlled circumstances and accuracy. Somewhat unusually, the history of the development of the rig, including the pitfalls was included in the paper. In addition to the hydrodynamic resistance measurements, the geometrical properties-metallic surface area, pore density, deployed and unconstrained length and diameter-of the stent deployment were measured. RESULTS: Based on our evaluation method a confidence band can be determined for a given deployment scenario. Collectively analysing the hydrodynamic resistance and the geometric indices, a deeper understanding of an implantation can be obtained. Our results suggest that to correctly interpret the hydrodynamic resistance of a scenario, the deployment length has to be considered. To demonstrate the applicability of the measurement, as a pilot study the results of four intracranial flow diverter stents of two types and sizes have been reported in this work. The results of these measurements even on this small sample size provide valuable information on differences between stent types and deployment scenarios.

Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH)-phase II: rupture risk assessment.

PURPOSE: Assessing the rupture probability of intracranial aneurysms (IAs) remains challenging. Therefore, hemodynamic simulations are increasingly applied toward supporting physicians during treatment planning. However, due to several assumptions, the clinical acceptance of these methods remains limited. METHODS: To provide an overview of state-of-the-art blood flow simulation capabilities, the Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH) was conducted. Seventeen research groups from all over the world performed segmentations and hemodynamic simulations to identify the ruptured aneurysm in a patient harboring five IAs. Although simulation setups revealed good similarity, clear differences exist with respect to the analysis of aneurysm shape and blood flow results. Most groups (12/71%) included morphological and hemodynamic parameters in their analysis, with aspect ratio and wall shear stress as the most popular candidates, respectively. RESULTS: The majority of groups (7/41%) selected the largest aneurysm as being the ruptured one. Four (24%) of the participating groups were able to correctly select the ruptured aneurysm, while three groups (18%) ranked the ruptured aneurysm as the second most probable. Successful selections were based on the integration of clinically relevant information such as the aneurysm site, as well as advanced rupture probability models considering multiple parameters. Additionally, flow characteristics such as the quantification of inflow jets and the identification of multiple vortices led to correct predictions. CONCLUSIONS: MATCH compares state-of-the-art image-based blood flow simulation approaches to assess the rupture risk of IAs. Furthermore, this challenge highlights the importance of multivariate analyses by combining clinically relevant metadata with advanced morphological and hemodynamic quantification.

A new hypothesis on the role of vessel topology in cerebral aneurysm initiation.

Aneurysm pathogenesis is thought to be strongly linked with hemodynamical effects. According to our current knowledge, the formation process is initiated by locally disturbed flow conditions. The aim of the current work is to provide a numerical investigation on the role of the flow field at the stage of the initiation, before the aneurysm formation. Digitally reconstructed pre-aneurysmal geometries are used to examine correlations of the flow patterns to the location and direction of the aneurysms formed later. We argue that a very specific rotational flow pattern is present in all the investigated cases marking the location of the later aneurysm and that these flow patterns provide the mechanical load on the wall that can lead to a destructive remodelling in the vessel wall. Furthermore, these patterns induce elevated vessel surface related variables (e.g. wall shear stress (WSS), wall shear stress gradient (WSSG) and oscillatory shear index (OSI)), in agreement with the previous findings. We emphasise that the analysis of the flow patterns provides a deeper insight and a more robust numerical methodology compared to the sole examination of the aforementioned surface quantities.

Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH): Phase I: Segmentation.

PURPOSE: Advanced morphology analysis and image-based hemodynamic simulations are increasingly used to assess the rupture risk of intracranial aneurysms (IAs). However, the accuracy of those results strongly depends on the quality of the vessel wall segmentation. METHODS: To evaluate state-of-the-art segmentation approaches, the Multiple Aneurysms AnaTomy CHallenge (MATCH) was announced. Participants carried out segmentation in three anonymized 3D DSA datasets (left and right anterior, posterior circulation) of a patient harboring five IAs. Qualitative and quantitative inter-group comparisons were carried out with respect to aneurysm volumes and ostia. Further, over- and undersegmentation were evaluated based on highly resolved 2D images. Finally, clinically relevant morphological parameters were calculated. RESULTS: Based on the contributions of 26 participating groups, the findings reveal that no consensus regarding segmentation software or underlying algorithms exists. Qualitative similarity of the aneurysm representations was obtained. However, inter-group differences occurred regarding the luminal surface quality, number of vessel branches considered, aneurysm volumes (up to 20%) and ostium surface areas (up to 30%). Further, a systematic oversegmentation of the 3D surfaces was observed with a difference of approximately 10% to the highly resolved 2D reference image. Particularly, the neck of the ruptured aneurysm was overrepresented by all groups except for one. Finally, morphology parameters (e.g., undulation and non-sphericity) varied up to 25%. CONCLUSIONS: MATCH provides an overview of segmentation methodologies for IAs and highlights the variability of surface reconstruction. Further, the study emphasizes the need for careful processing of initial segmentation results for a realistic assessment of clinically relevant morphological parameters.

Measurement of flow diverter hydraulic resistance to model flow modification in and around intracranial aneurysms.

Flow diverters (FDs) have been successfully applied in the recent decade to the treatment of intracranial aneurysms by impairing the communication between the flows in the parent artery and the aneurysm and, thus, the blood within the aneurysm sac. It would be desirable to have a simple and accurate computational method to follow the changes in the peri- and intraaneurysmal flow caused by the presence of FDs. The detailed flow simulation around the intricate wire structure of the FDs has three disadvantages: need for high amount of computational resources and highly skilled professionals to prepare the computational grid, and also the lack of validation that makes the invested effort questionable. In this paper, we propose a porous layer method to model the hydraulic resistance (HR) of one or several layers of the FDs. The basis of this proposal is twofold: first, from an application point of view, the only interesting parameter regarding the function of the FD is its HR; second, we have developed a method to measure the HR with a simple apparatus. We present the results of these measurements and demonstrate their utility in numerical simulations of patient-specific aneurysm simulations.

Impact of aneurysmal geometry on intraaneurysmal flow: a computerized flow simulation study.

INTRODUCTION: This study was performed to assess the effect of aneurysm geometry on parameters that may have an impact on the natural history of intracranial aneurysms, such as intraaneurysmal flow pressure and shear stress. METHODS: Flow was simulated in 21 randomly selected aneurysms using finite volume modeling. Ten aneurysms were classified as side-wall aneurysms, with either single-sided or circumferential involvement of the parent artery wall, and 11 as bifurcation aneurysms (symmetric or asymmetric), with an axis either perpendicular or parallel to the parent artery. The flow patterns were classified as either jet or vortex types (with regular or irregular vortex flow). Pressures and shear stresses were characterized as evenly or unevenly distributed over the aneurysm wall and neck. RESULTS: All side-wall and four of the bifurcation aneurysms with a perpendicular axis had a vortex type flow pattern and seven bifurcation aneurysms with a parallel axis (four symmetric and two asymmetric) had a jet flow pattern. Jet type flow was associated with an uneven pressure distribution in seven out of seven aneurysms. Vortex type flow resulted in an even pressure distribution in five out of six aneurysms with an irregular flow pattern and six out of eight with a regular flow pattern. No firm relationship could be established between any geometrical type and shear stress distribution. Only 1 of 14 aneurysms with a perpendicular axis, but 4 of 7 aneurysms with a parallel axis, had ruptured. CONCLUSION: Aneurysm geometry does have an impact on flow conditions. Aneurysms with a main axis parallel to the parent artery have a tendency to have a jet flow pattern and uneven distribution of unsteady pressure. These aneurysms may have a higher rate of rupture as than those with a main axis perpendicular to the parent artery.