Gábor Janiga

Optimization and computational fluid dynamics

Generalities and methods.- A Few Illustrative Examples of CFD-based Optimization.- Mathematical Aspects of CFD-based Optimization.- Adjoint Methods for Shape Optimization.- Specific Applications of CFD-based Optimization to Engineering Problems.- Efficient Deterministic Approaches for Aerodynamic Shape Optimization.- Numerical Optimization for Advanced Turbomachinery Design.- CFD-based Optimization for Automotive Aerodynamics.- Multi-objective Optimization for Problems Involving Convective Heat Transfer.- CFD-based Optimization for a Complete Industrial Process: Papermaking.

Variability of Computational Fluid Dynamics Solutions for Pressure and Flow in a Giant Aneurysm: The ASME 2012 Summer Bioengineering Conference CFD Challenge

Stimulated by a recent controversy regarding pressure drops predicted in a giant aneurysm with a proximal stenosis, the present study sought to assess variability in the prediction of pressures and flow by a wide variety of research groups. In phase I, lumen geometry, flow rates, and fluid properties were specified, leaving each research group to choose their solver, discretization, and solution strategies. Variability was assessed by having each group interpolate their results onto a standardized mesh and centerline. For phase II, a physical model of the geometry was constructed, from which pressure and flow rates were measured. Groups repeated their simulations using a geometry reconstructed from a micro-computed tomography (CT) scan of the physical model with the measured flow rates and fluid properties. Phase I results from 25 groups demonstrated remarkable consistency in the pressure patterns, with the majority predicting peak systolic pressure drops within 8% of each other. Aneurysm sac flow patterns were more variable with only a few groups reporting peak systolic flow instabilities owing to their use of high temporal resolutions. Variability for phase II was comparable, and the median predicted pressure drops were within a few millimeters of mercury of the measured values but only after accounting for submillimeter errors in the reconstruction of the life-sized flow model from micro-CT. In summary, pressure can be predicted with consistency by CFD across a wide range of solvers and solution strategies, but this may not hold true for specific flow patterns or derived quantities. Future challenges are needed and should focus on hemodynamic quantities thought to be of clinical interest.

Impact of Stents and Flow Diverters on Hemodynamics in Idealized Aneurysm Models

Cerebral aneurysms constitute a major medical challenge as treatment options are limited and often associated with high risks. Statistically, up to 3% of patients with a brain aneurysm may suffer from bleeding for each year of life. Eight percent of all strokes are caused by ruptured aneurysms. In order to prevent this rupture, endovascular stenting using so called flow diverters is increasingly being regarded as an alternative to the established coil occlusion method in minimally invasive treatment. Covering the neck of an aneurysm with a flow diverter has the potential to alter the hemodynamics in such a way as to induce thrombosis within the aneurysm sac, stopping its further growth, preventing its rupture and possibly leading to complete resorption. In the present study the influence of different flow diverters is quantified considering idealized patient configurations, with a spherical sidewall aneurysm placed on either a straight or a curved parent vessel. All important hemodynamic parameters (exchange flow rate, velocity, and wall shear stress) are determined in a quantitative and accurate manner using computational fluid dynamics when varying the key geometrical properties of the aneurysm. All simulations are carried out using an incompressible, Newtonian fluid with steady conditions. As a whole, 72 different cases have been considered in this systematic study. In this manner, it becomes possible to compare the efficiency of different stents and flow diverters as a function of wire density and thickness. The results show that the intra-aneurysmal flow velocity, wall shear stress, mean velocity, and vortex topology can be considerably modified thanks to insertion of a suitable implant. Intra-aneurysmal residence time is found to increase rapidly with decreasing stent porosity. Of the three different implants considered in this study, the one with the highest wire density shows the highest increase of intra-aneurysmal residence time for both the straight and the curved parent vessels. The best hemodynamic modifications are always obtained for a small aneurysm diameter.

Realistic virtual intracranial stenting and computational fluid dynamics for treatment analysis

In order to support the decisions of medical experts and to develop better stent designs, the availability of a simulation tool for virtual stenting would be extremely useful. An innovative virtual stenting technique is described in this work, which is directly applicable for complex patient-specific geometries. A basilar tip aneurysm provided for the Virtual Intracranial Stenting Challenge 2010 is considered to demonstrate the advantages of this approach. A free-form deformation is introduced for a wall-tight stent deployment. Numerical flow simulations on sufficiently fine computational meshes are performed for different configurations in order to characterize the inflow rate into the aneurysm and the corresponding residence time in the aneurysm sac. A Neuroform and a SILK stent have been deployed at various locations and the computed residence times have been evaluated and compared, demonstrating the advantage associated with a lower stent porosity. It has been found that the SILK stent leads to a large increase in the residence time and to a significant reduction in the maximum wall shear stress in the aneurysm sac. This is only observed when placing the stent in the appropriate position, showing that virtual stenting might be employed for operation support.

Cerebral Blood Flow in a Healthy Circle of Willis and Two Intracranial Aneurysms: Computational Fluid Dynamics Versus Four-Dimensional Phase-Contrast Magnetic Resonance Imaging

Computational fluid dynamics (CFD) opens up multiple opportunities to investigate the hemodynamics of the human vascular system. However, due to numerous assumptions the acceptance of CFD among physicians is still limited in practice and validation through comparison is mandatory. Time-dependent quantitative phase-contrast magnetic resonance imaging PC-MRI measurements in a healthy volunteer and two intracranial aneurysms were carried out at 3 and 7 Tesla. Based on the acquired images, three-dimensional (3D) models of the aneurysms were reconstructed and used for the numerical simulations. Flow information from the MR measurements were applied as boundary conditions. The four-dimensional (4D) velocity fields obtained by CFD and MRI were qualitatively as well as quantitatively compared including cut planes and vector analyses. For all cases a high similarity of the velocity patterns was observed. Additionally, the quantitative analysis revealed a good agreement between CFD and MRI. Deviations were caused by minor differences between the reconstructed vessel models and the actual lumen. The comparisons between diastole and systole indicate that relative differences between MRI and CFD are intensified with increasing velocity. The findings of this study lead to the conclusion that CFD and MRI agree well in predicting intracranial velocities when realistic geometries and boundary conditions are provided. Due to the considerably higher temporal and spatial resolution of CFD compared to MRI, complex flow patterns can be further investigated in order to evaluate their role with respect to aneurysm formation or rupture. Nevertheless, special care is required regarding the vessel reconstruction since the geometry has a major impact on the subsequent numerical results.

The Computational Fluid Dynamics Rupture Challenge 2013—Phase I: Prediction of Rupture Status in Intracranial Aneurysms

BACKGROUND AND PURPOSE: Rupture risk assessment for intracranial aneurysms remains challenging, and risk factors, including wall shear stress, are discussed controversially. The primary purpose of the presented challenge was to determine how consistently aneurysm rupture status and rupture site could be identified on the basis of computational fluid dynamics. MATERIALS AND METHODS: Two geometrically similar MCA aneurysms were selected, 1 ruptured, 1 unruptured. Participating computational fluid dynamics groups were blinded as to which case was ruptured. Participants were provided with digitally segmented lumen geometries and, for this phase of the challenge, were free to choose their own flow rates, blood rheologies, and so forth. Participants were asked to report which case had ruptured and the likely site of rupture. In parallel, lumen geometries were provided to a group of neurosurgeons for their predictions of rupture status and site. RESULTS: Of 26 participating computational fluid dynamics groups, 21 (81%) correctly identified the ruptured case. Although the known rupture site was associated with low and oscillatory wall shear stress, most groups identified other sites, some of which also experienced low and oscillatory shear. Of the 43 participating neurosurgeons, 39 (91%) identified the ruptured case. None correctly identified the rupture site. CONCLUSIONS: Geometric or hemodynamic considerations favor identification of rupture status; however, retrospective identification of the rupture site remains a challenge for both engineers and clinicians. A more precise understanding of the hemodynamic factors involved in aneurysm wall pathology is likely required for computational fluid dynamics to add value to current clinical decision-making regarding rupture risk.

Endothelialization of over- and undersized flow-diverter stents at covered vessel side branches: An in vivo and in silico study.

Although flow-diverting devices are promising treatment options for intracranial aneurysms, jailed side branches might occlude leading to insufficient blood supply. Especially differences in the local stent strut compression may have a drastic influence on subsequent endothelialization. To investigate the outcome of different treatment scenarios, over- and undersized stent deployments were realized experimentally and computationally. Two Pipeline Embolization Devices were placed in the right common carotid artery of large white swine, crossing the right ascending pharyngeal artery. DSA and PC-MRI measurements were acquired pre- and post-stenting and after three months. To evaluate the stent strut endothelialization and the corresponding ostium patency, the swine were sacrificed and scanning electron microscopy measurements were carried out. A more detailed analysis of the near-stent hemodynamics was enabled by a realistic virtual stenting in combination with highly resolved Computational Fluid Dynamics simulations using case-specific boundary conditions. The oversizing resulted in an elongated stent deployment with more open stent pores, while for the undersized case a shorter deployment with more condensed pores was present. In consequence, the side branch of the first case remained patent after three months and the latter almost fully occluded. The virtual investigation confirmed the experimental findings by identifying differences between the individual velocities as well as stent shear stresses at the distal part of the ostia. The choice of flow-diverting device and the subsequent deployment strategy strongly influences the patency of jailed side branches. Therefore, careful treatment planning is required, to guarantee sufficient blood supply in the brain territories supplied those branches.

Automatic Detection and Visualization of Qualitative Hemodynamic Characteristics in Cerebral Aneurysms

Cerebral aneurysms are a pathological vessel dilatation that bear a high risk of rupture. For the understanding and evaluation of the risk of rupture, the analysis of hemodynamic information plays an important role. Besides quantitative hemodynamic information, also qualitative flow characteristics, e.g., the inflow jet and impingement zone are correlated with the risk of rupture. However, the assessment of these two characteristics is currently based on an interactive visual investigation of the flow field, obtained by computational fluid dynamics (CFD) or blood flow measurements. We present an automatic and robust detection as well as an expressive visualization of these characteristics. The detection can be used to support a comparison, e.g., of simulation results reflecting different treatment options. Our approach utilizes local streamline properties to formalize the inflow jet and impingement zone. We extract a characteristic seeding curve on the ostium, on which an inflow jet boundary contour is constructed. Based on this boundary contour we identify the impingement zone. Furthermore, we present several visualization techniques to depict both characteristics expressively. Thereby, we consider accuracy and robustness of the extracted characteristics, minimal visual clutter and occlusions. An evaluation with six domain experts confirms that our approach detects both hemodynamic characteristics reasonably.

Blood Flow Clustering and Applications inVirtual Stenting of Intracranial Aneurysms

Understanding the hemodynamics of blood flow in vascular pathologies such as intracranial aneurysms is essential for both their diagnosis and treatment. Computational fluid dynamics (CFD) simulations of blood flow based on patient-individual data are performed to better understand aneurysm initiation and progression and more recently, for predicting treatment success. In virtual stenting, a flow-diverting mesh tube (stent) is modeled inside the reconstructed vasculature and integrated in the simulation. We focus on steady-state simulation and the resulting complex multiparameter data. The blood flow pattern captured therein is assumed to be related to the success of stenting. It is often visualized by a dense and cluttered set of streamlines.We present a fully automatic approach for reducing visual clutter and exposing characteristic flow structures by clustering streamlines and computing cluster representatives. While individual clustering techniques have been applied before to streamlines in 3D flow fields, we contribute a general quantitative and a domain-specific qualitative evaluation of three state-of-the-art techniques. We show that clustering based on streamline geometry as well as on domain-specific streamline attributes contributes to comparing and evaluating different virtual stenting strategies. With our work, we aim at supporting CFD engineers and interventional neuroradiologists.

An automatic CFD-based flow diverter optimization principle for patient-specific intracranial aneurysms

The optimal treatment of intracranial aneurysms using flow diverting devices is a fundamental issue for neuroradiologists as well as neurosurgeons. Due to highly irregular manifold aneurysm shapes and locations, the choice of the stent and the patient-specific deployment strategy can be a very difficult decision. To support the therapy planning, a new method is introduced that combines a three-dimensional CFD-based optimization with a realistic deployment of a virtual flow diverting stent for a given aneurysm. To demonstrate the feasibility of this method, it was applied to a patient-specific intracranial giant aneurysm that was successfully treated using a commercial flow diverter. Eight treatment scenarios with different local compressions were considered in a fully automated simulation loop. The impact on the corresponding blood flow behavior was evaluated qualitatively as well as quantitatively, and the optimal configuration for this specific case was identified. The virtual deployment of an uncompressed flow diverter reduced the inflow into the aneurysm by 24.4% compared to the untreated case. Depending on the positioning of the local stent compression below the ostium, blood flow reduction could vary between 27.3% and 33.4%. Therefore, a broad range of potential treatment outcomes was identified, illustrating the variability of a given flow diverter deployment in general. This method represents a proof of concept to automatically identify the optimal treatment for a patient in a virtual study under certain assumptions. Hence, it contributes to the improvement of virtual stenting for intracranial aneurysms and can support physicians during therapy planning in the future.