Immersed boundary-finite element model of fluid-structure interaction in the aortic root

It has long been recognized that aortic root elasticity helps to ensure efficient aortic valve closure, but our understanding of the functional importance of the elasticity and geometry of the aortic root continues to evolve as increasingly detailed in vivo imaging data become available. Herein, we describe fluid-structure interaction models of the aortic root, including the aortic valve leaflets, the sinuses of Valsalva, the aortic annulus, and the sinotubular junction, that employ a version of Peskin's immersed boundary (IB) method with a finite element (FE) description of the structural elasticity. We develop both an idealized model of the root with three-fold symmetry of the aortic sinuses and valve leaflets, and a more realistic model that accounts for the differences in the sizes of the left, right, and noncoronary sinuses and corresponding valve cusps. As in earlier work, we use fiber-based models of the valve leaflets, but this study extends earlier IB models of the aortic root by employing incompressible hyperelastic models of the mechanics of the sinuses and ascending aorta using a constitutive law fit to experimental data from human aortic root tissue. In vivo pressure loading is accounted for by a backwards displacement method that determines the unloaded configurations of the root models. Our models yield realistic cardiac output at physiological pressures, with low transvalvular pressure differences during forward flow, minimal regurgitation during valve closure, and realistic pressure loads when the valve is closed during diastole. Further, results from high-resolution computations demonstrate that IB models of the aortic valve are able to produce essentially grid-converged dynamics at practical grid spacings for the high-Reynolds number flows of the aortic root

Estimation of wall shear stress using 4D flow cardiovascular MRI and computational fluid dynamics

Electronic version of an article published as Journal of mechanics in medicine and biology, 0, 1750046 (2016), 16 pages. DOI:10.1142/S0219519417500464 © World Scientific Publishing CompanyIn the last few years, wall shear stress (WSS) has arisen as a new diagnostic indicator in patients with arterial disease. There is a substantial evidence that the WSS plays a significant role, together with hemodynamic indicators, in initiation and progression of the vascular diseases. Estimation of WSS values, therefore, may be of clinical significance and the methods employed for its measurement are crucial for clinical community. Recently, four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) has been widely used in a number of applications for visualization and quantification of blood flow, and although the sensitivity to blood flow measurement has increased, it is not yet able to provide an accurate three-dimensional (3D) WSS distribution. The aim of this work is to evaluate the aortic blood flow features and the associated WSS by the combination of 4D flow cardiovascular magnetic resonance (4D CMR) and computational fluid dynamics technique. In particular, in this work, we used the 4D CMR to obtain the spatial domain and the boundary conditions needed to estimate the WSS within the entire thoracic aorta using computational fluid dynamics. Similar WSS distributions were found for cases simulated. A sensitivity analysis was done to check the accuracy of the method. 4D CMR begins to be a reliable tool to estimate the WSS within the entire thoracic aorta using computational fluid dynamics. The combination of both techniques may provide the ideal tool to help tackle these and other problems related to wall shear estimation.Peer ReviewedPostprint (author's final draft

Stress-strain analysis of aortic aneurysms

Tato práce se zabývá problematikou aneurysmat břišní aorty a možností využít konečnoprvkovou deformačně-napěťovou analýzu těchto aneurysmat ke stanovení rizika ruptury. První část práce je věnována úvodu do problematiky, popisu kardiovaskulární soustavy člověka s důrazem na abdominální aortu, anatomii, fyziologii a patologii stěny tepny s důrazem na procesy vedoucí ke vzniku aneurysmatu. Dále se práce věnuje rizikovým faktorům přispívajících ke vzniku aneurysmat spolu s analýzou současných klinických postupů ke stanovení rizika ruptury spolu se srovnáním navrhovaného kritéria maximálního napětí. Dominantní část této disertace je věnována identifikaci faktorů ovlivňujících napjatost a deformaci stěny aneurysmatu spolu s návrhem nových postupů, prezentací vlastních poznatků vedoucích ke zpřesnění určení rizika ruptury pomocí deformačně- napěťové analýzy a metody konečných prvků. Nejprve je analyzován vliv geometrie, vedoucí k závěru, že je nezbytné používání individuálních geometrií pacienta. Dále je pozornost zaměřena na odbočující tepny, které ve stěně působí jako koncentrátor napětí a mohou tedy ovlivňovat napjatost v ní. Jako další podstatný faktor byl identifikován vliv nezatížené geometrie a bylo napsáno makro pro její nalezení, které bylo opět zahrnuto jako standardní součást do výpočtového modelu. Mechanické vlastnosti jak stěny aneurysmatu, tak intraluminálního trombu jsou experimentálně testovány pomocí dvouosých zkoušek. Také je zde analyzován vliv modelu materiálu, kde je ukázáno, že srovnávání maximálních napětí u jednotlivých modelů materiálu není vhodné díky zcela rozdílným gradientům napětí ve stěně aneurysmatu. Dále je zdůrazněna potřeba znalosti distribuce kolagenních vláken ve stěně a navržen program k jejímu získání. Intraluminální trombus je analyzován ve dvou souvislostech. Jednak je ukázán vliv jeho ruptury na napětí ve stěně a jednak je analyzován vliv jeho poroelastické struktury na totéž. Posledním identifikovaným podstatným faktorem je zbytková napjatost ve stěně. Její významnost je demonstrována na několika aneurysmatech a i tato je zahrnuta jako integrální součást do našeho výpočtového modelu.Na závěr jsou pak navrženy další možné směry výzkumu.This thesis deals with abdominal aortic aneurysms and the possibility of using finite element method in assessment of their rupture risk. First part of the thesis is dedicated to an introduction into the problem, description of human cardiovascular system where the abdominal aorta, its anatomy, physiology and pathology is emphasized. There Processes leading to formationing of abdominal aortic aneurysms are also discussed. Risk factors contributing to creation of aneurysms are discussed next. Finally, an analysis of current clinical criteria which determine rupture risk of an abdominal aortic aneurysm is presented and compared with the new maximum stress criterion being currently in development. Main part of the thesis deals with the identification of relevant factors which affect stress and deformation of aneurysmal wall. This is connected with proposals of new approaches leading to predicting the rupture risk more accurately by using finite element stress-strain analysis. The impact of geometry is analyzed first with the conclusion that patient-specific geometry is a crucial input in the computational model. Therefore its routine reconstruction has been managed. Attention is then paid to the branching arteries which were neglected so far although they cause a stress concentration in arterial wall. The necessity of knowing the unloaded geometry of aneurysm is then emphasized. Therefore a macro has been written in order to be able to find the unloaded geometry for any patient-specific geometry of aneurysm. Mechanical properties of both aneurysmal wall and intraluminal thrombus were also experimentally tested and their results were fitted by an isotropic material model. The effect of the material model itself has been also investigated by comparing whole stress fields of several aneurysms. It has been shown that different models predict completely different stresses due to different stress gradients in the aneurysmal wall. The necessity of known collagen fiber distribution in arterial wall is also emphasized. A special program is then presented enabling us to obtain this information. Effect of intraluminal thrombus on the computed wall stress is analyzed in two perspectives. First the effect of its failure on wall stress is shown and also the impact of its poroelastic structure is analyzed. Finally the residual stresses were identified as an important factor influencing the computed wall stress in aneurysmal wall and they were included into patient-specific finite element analysis of aneurysms. Further possible regions of investigation are mentioned as the last part of the thesis.

Comparative finite element modelling of aneurysm formation and physiologic inflation in the descending aorta

Despite the general interest in aneurysm rupture prediction, the aneurysm formation has received limited attention. The goal of this study is to assess whether an aneurysm may be instigated in a healthy model of an aorta inflated by a supra-physiological pressure. The effect of two main aspects on numerical predictions has been explored: i) the geometric design and ii) the constitutive law adopted to represent the material properties. Firstly, higher values of wall stress and displacement magnitude were generated in the physiologic model compared to the cylindrical one when assigning the same material properties. Secondly, greater deformations are observed in the anisotropic model compared to the isotropic one

On Coupling a Lumped Parameter Heart Model and a Three-Dimensional Finite Element Aorta Model

Aortic flow and pressure result from the interactions between the heart and arterial system. In this work, we considered these interactions by utilizing a lumped parameter heart model as an inflow boundary condition for three-dimensional finite element simulations of aortic blood flow and vessel wall dynamics. The ventricular pressure–volume behavior of the lumped parameter heart model is approximated using a time varying elastance function scaled from a normalized elastance function. When the aortic valve is open, the coupled multidomain method is used to strongly couple the lumped parameter heart model and three-dimensional arterial models and compute ventricular volume, ventricular pressure, aortic flow, and aortic pressure. The shape of the velocity profiles of the inlet boundary and the outlet boundaries that experience retrograde flow are constrained to achieve a robust algorithm. When the aortic valve is closed, the inflow boundary condition is switched to a zero velocity Dirichlet condition. With this method, we obtain physiologically realistic aortic flow and pressure waveforms. We demonstrate this method in a patient-specific model of a normal human thoracic aorta under rest and exercise conditions and an aortic coarctation model under pre- and post-interventions