Mariano Vázquez

Analysis of hemodynamics and wall mechanics at sites of cerebral aneurysm rupture.

Background It is thought that aneurysms evolve as the result of progressive degradation of the wall in response to abnormal hemodynamics characterized by either high or low wall shear stress (WSS). Objective To investigate the effects of these two different hemodynamic pathways in a series of cerebral aneurysms with known rupture sites. Methods Nine aneurysms in which the rupture site could be identified in three-dimensional images were analyzed. The WSS distribution was obtained from computational fluid dynamics (CFD) simulations. Internal wall stresses were computed using structural wall models under hemodynamic loads determined by the CFD models. Wall properties (thickness and stiffness) were modulated with the WSS distribution (increased or decreased in regions of high or low WSS) to test possible wall degradation pathways. Rupture probability indices (RPI) were calculated to compare different wall models. Results Most rupture sites aligned with the intrasaccular flow stream and downstream of the primary impaction zone. The model that best explained the rupture site (produced higher RPI) in eight of the nine aneurysms (89%) had thinner and stiffer walls in regions of abnormally high WSS. The remaining case (11%) was best explained by a model with thinner and stiffer walls in regions of abnormally low WSS. Conclusions Aneurysm rupture seems to be caused by localized degradation and weakening of the wall in response to abnormal hemodynamics. Image-based computational models assuming wall thinning and stiffening in regions of abnormally high WSS were able to explain most of the observed rupture sites.

Large-scale CFD simulations of the transitional and turbulent regime for the large human airways during rapid inhalation

The dynamics of unsteady flow in the human large airways during a rapid inhalation were investigated using highly detailed large-scale computational fluid dynamics on a subject-specific geometry. The simulations were performed to resolve all the spatial and temporal scales of the flow, thanks to the use of massive computational resources. A highly parallel finite element code was used, running on two supercomputers, solving the transient incompressible Navier-Stokes equations on unstructured meshes. Given that the finest mesh contained 350 million elements, the study sets a precedent for large-scale simulations of the respiratory system, proposing an analysis strategy for mean flow, fluctuations and wall shear stresses on a rapid and short inhalation (a so-called sniff). The geometry used encompasses the exterior face and the airways from the nasal cavity, through the trachea and up to the third lung bifurcation; it was derived from a contrast-enhanced computed tomography (CT) scan of a 48-year-old male. The transient inflow produces complex flows over a wide range of Reynolds numbers (Re). Thanks to the high fidelity simulations, many features involving the flow transition were observed, with the level of turbulence clearly higher in the throat than in the nose. Spectral analysis revealed turbulent characteristics persisting downstream of the glottis, and were captured even with a medium mesh resolution. However a fine mesh resolution was found necessary in the nasal cavity to observe transitional features. This work indicates the potential of large-scale simulations to further understanding of airway physiological mechanics, which is essential to guide clinical diagnosis; better understanding of the flow also has implications for the design of interventions such as aerosol drug delivery.

What a difference in biomechanics cardiac fiber makes

Computational simulations of the heart are a powerful tool for a comprehensive understanding of cardiac function and its intrinsic relationship with its muscular architecture. Cardiac biomechanical models require a vector field representing the orientation of cardiac fibers. A wrong orientation of the fibers can lead to a non-realistic simulation of the heart functionality. In this paper we explore the impact of the fiber information on the simulated biomechanics of cardiac muscular anatomy. We have used the John Hopkins database to perform a biomechanical simulation using both a synthetic benchmark fiber distribution and the data obtained experimentally from DTI. Results illustrate how differences in fiber orientation affect heart deformation along cardiac cycle.

Dynamic load balance applied to particle transport in fluids

This work presents a parallel numerical strategy to transport Lagrangian particles in a fluid using a dynamic load balance strategy. Both fluid and particle solvers are parallel, with two levels of parallelism. The first level is based on a substructuring technique and uses message passing interface (MPI) as the communication library; the second level consists of OpenMP pragmas for loop parallelisation at the node level. When dealing with transient flows, there exist two main alternatives to address the coupling of these solvers. On the one hand, a single-code approach consists in solving the particle equations once the fluid solution has been obtained at the end of a time step, using the same instance of the same code. On the other hand, a multi-code approach enables one to overlap the transport of the particles with the next time-step solution of the fluid equations, and thus obtain asynchronism. In this case, different codes or two instances of the same code can be used. Both approaches will be presented. In addition, a dynamic load balancing library is used on the top of OpenMP pragmas in order to continuously exploit all the resources available at the node level, thus increasing the load balance and the efficiency of the parallelisation and uses the MPI.

Hybrid MPI-OpenMP performance in massively parallel computational fluid dynamics

We analyze in this paper the parallel performance of a computational mechanics code, based on a hybrid MPI-OpenMP parallelization. The specific problems considered are an explicit compressible flow solver and an implicit incompressible flow solver. Their performance using a fully MPI approach and a Hybrid approach are compared. The performance obtained on different architextures are also compared.

Total energy conservation in ALE schemes for compressible flows

The numerical prediction of interaction phenomena between a compressible flow model with a moving domain and other physical models requires that the work performed on the fluid is properly translated into total fluid energy variation. We present a numerical model relying on an Arbitrary Lagrangian-Eulerian (ALE) unstructured vertex-centered finite volume that satisfies this condition together with the Geometric Conservation Law. We apply this numerical scheme to the solution of a 3D fluid-structure interaction problem. The results are contrasted with those obtained by the energy non-conservative counterpart.

Implications of bipolar voltage mapping and magnetic resonance imaging resolution in biventricular scar characterization after myocardial infarction

Aims We aimed to study the differences in biventricular scar characterization using bipolar voltage mapping compared with state-of-the-art in vivo delayed gadolinium-enhanced cardiac magnetic resonance (LGE-CMR) imaging and ex vivo T1 mapping. Methods and results Ten pigs with established myocardial infarction (MI) underwent in vivo scar characterization using LGE-CMR imaging and high-density voltage mapping of both ventricles using a 3.5-mm tip catheter. Ex vivo post-contrast T1 mapping provided a high-resolution reference. Voltage maps were registered onto the left and right ventricular (LV and RV) endocardium, and epicardium of CMR-based geometries to compare voltage-derived scars with surface-projected 3D scars. Voltage-derived scar tissue of the LV endocardium and the epicardium resembled surface projections of 3D in vivo and ex vivo CMR-derived scars using 1-mm of surface projection distance. The thinner wall of the RV was especially sensitive to lower resolution in vivo LGE-CMR images, in which differences between normalized low bipolar voltage areas and CMR-derived scar areas did not decrease below a median of 8.84% [interquartile range (IQR) (3.58, 12.70%)]. Overall, voltage-derived scars and surface scar projections from in vivo LGE-CMR sequences showed larger normalized scar areas than high-resolution ex vivo images [12.87% (4.59, 27.15%), 18.51% (11.25, 24.61%), and 9.30% (3.84, 19.59%), respectively], despite having used optimized surface projection distances. Importantly, 43.02% (36.54, 48.72%) of voltage-derived scar areas from the LV endocardium were classified as non-enhanced healthy myocardium using ex vivo CMR imaging. Conclusion In vivo LGE-CMR sequences and high-density voltage mapping using a conventional linear catheter fail to provide accurate characterization of post-MI scar, limiting the specificity of voltage-based strategies and imaging-guided procedures.

A Partitioned Methodology for Conjugate Heat Transfer on Dynamic Structures

A partitioned coupling approach for conjugate heat transfer applications is presented. The coupling scheme is based on the extension of the parallel algebraic domain composition method already validated in fluid-structure interactions problems for thermal coupling. The method alters the original Dirichlet-Neumann approach enforcing the boundary conditions over the subdomains through matrix operations. The algorithm is tested on two benchmark cases with conjugate heat transfer: flow over a heated cylinder and flow over a flat-plate. The results indicate good agreement with previous research and encourages its application for large-scale problems.

Study of the Wall Thermal Condition Effect in a Lean-Premixed Downscaled Can Combustor Using Large-Eddy Simulation

The primary purpose of this study is to evaluate the ability of LES, with a turbulent combustion model based on steady flamelets, to predict the flame stabilization mechanisms in an industrial can combustor at full load conditions. The test case corresponds to the downscaled Siemens can combustor tested in the high pressure rig at the DLR. The effects of the wall temperature on the prediction capabilities of the codes is investigated by imposing several heat transfer conditions at the pilot and chamber walls. The codes used for this work are Alya and OpenFOAM, which are well established CFD codes in the fluid mechanics community. Prior to the simulation, results for 1-D laminar flames at the operating conditions of the combustor are compared with the detailed solutions. Subsequently, results from both codes at the mid-plane are compared against the experimental data available. Acceptable results are obtained for the axial velocity, while discrepancies are more evident for the mixture fraction and the temperature, particularly with Alya. However, both codes showed that the heat losses influence the size and length of the pilot and main flame.

Fully-Coupled Electromechanical Simulations of the LV Dog Anatomy Using HPC: Model Testing and Verification

Verification of electro-mechanic models of the heart require a good amount of reliable, high resolution, thorough in-vivo measurements. The detail of the mathematical models used to create simulations of the heart beat vary greatly. Generally, the objective of the simulation determines the modeling approach. However, it is important to exactly quantify the amount of error between the various approaches that can be used to simulate a heart beat by comparing them to ground truth data. The more detailed the model is, the more computing power it requires, we therefore employ a high-performance computing solver throughout this study. We aim to compare models to data measured experimentally to identify the effect of using a mathematical model of fibre orientation versus the measured fibre orientations using DT-MRI. We also use simultaneous endocardial stimuli vs an instantaneous myocardial stimulation to trigger the mechanic contraction. Our results show that synchronisation of the electrical and mechanical events in the heart beat are necessary to create a physiological timing of hemodynamic events. Synchronous activation of all of the myocardium provides an unrealistic timing of hemodynamic events in the cardiac cycle. Results also show the need of establishing a protocol to quantify the zero-pressure configuration of the left ventricular geometry to initiate the simulation protocol; however, the predicted zero-pressure configuration of the same geometry was different, depending on the origin of the fibre field employed.