A numerical study of isotropic and anisotropic constitutive models with relevance to healthy and unhealthy cerebral arterial tissues

This paper presents an analysis of hyperelastic constitutive models for continuous bodies both from a modeling and numerical point of view. Contributions are made within the context of finite element numerical simulations. Numerical results with relevance to flows in the cardiovascular system are outlined in the case of a sophisticated fluid–structure interaction problem, in specific complex geometries of anatomically accurate cerebral arteries in diseased state. In this regard, the work carefully outlines the numerical validation of constitutive models for healthy and unhealthy cerebral arterial tissues by means of simulations of static inflation tests on an idealized specimen of anterior cerebral artery (ACA). The healthy tissue is described by means of isotropic and anisotropic models that, are fitted with respect to experimental data describing the mechanical behavior of the ACA; the numerical results are presented highlighting the most important numerical aspects influencing the correct and efficient simulation of the mechanics of continuous bodies such as, for instance, the arterial wall. We further consider numerical simulations of unhealthy conditions of the tissue by taking into account different levels of weakening of its mechanical properties. Taking the cerebral cardiovascular system as a challenging test problem, we focus on the study of the effects of the imposed mechanical levels of degradation on kinematic quantities of interest by simulating static inflation tests for the different models. This work does not aim to propose a new mathematical model for the mechanical damage occurring at the onset of cardiovascular diseases such as cerebral aneurysms. The modeling and numerical techniques presented may be applied to a wide range of problems, equally challenging to that of the cardiovascular system with complex structural models and fluid-structure coupling.

Simulation of shear-driven flows

Roland Bouffanais

The research work reported in the present dissertation is aimed at the analysis of complex physical phenomena involving instabilities and nonlinearities occurring in fluids through state-of-the-art numerical modeling. Solutions of intricate fluid physics problems are devised in two particularly arduous situations: fluid domains with moving boundaries and the high-Reynolds-number regime dominated by nonlinear convective effects. Shear-driven flows of incompressible Newtonian fluids enclosed in cavities of varying geometries are thoroughly investigated in the two following frameworks: transition with a free surface and confined turbulence. The physical system we consider is made of an incompressible Newtonian fluid filling a bounded, or partially bounded cavity. A series of shear-driven flows are easily generated by setting in motion some part of the container boundary. These driven-cavity flows are not only technologically important, they are of great scientific interest because they display almost all physical fluid phenomena that can possibly occur in incompressible flows, and this in the simplest geometrical settings. Thus corner eddies, secondary flows, longitudinal vortices, complex three-dimensional patterns, chaotic particle motions, nonuniqueness, transition, and turbulence all occur naturally and can be studied in the same geometry. This facilitates the comparison of results from experiments, analysis, and computation over the whole range of Reynolds numbers. The flows under consideration are part of a larger class of confined flows driven by linear or angular momentum gradients. This dissertation reports a detailed study of a novel numerical method developed for the simulation of an unsteady free-surface flow in three-space-dimensions. This method relies on a moving-grid technique to solve the Navier-Stokes equations expressed in the arbitrary Lagrangian-Eulerian (ALE) kinematics and discretized by the spectral element method. A comprehensive analysis of the continuous and discretized formulations of the general problem in the ALE frame, with nonlinear, non-homogeneous and unsteady boundary conditions is presented. In this dissertation, we first consider in the internal turbulent flow of a fluid enclosed in a bounded cubical cavity driven by the constant translation of its lid. The solution of this flow relied on large-eddy simulations, which served to improve our physical understanding of this complex flow dynamics. Subsequently, a novel subgrid model based on approximate deconvolution methods coupled with a dynamic mixed scale model was devised. The large-eddy simulation of the lid-driven cubical cavity flow based on this novel subgrid model has shown improvements over traditional subgrid-viscosity type of models. Finally a new interpretation of approximate deconvolution models when used with implicit filtering as a way to approximate the projective grid filter was given. This led to the introduction of the grid filter models. Through the use of a newly-developed method of numerical simulation, in this dissertation we solve unsteady flows with a flat and moving free-surface in the transitional regime. These flows are the incompressible flow of a viscous fluid enclosed in a cylindrical container with an open top surface and driven by the steady rotation of the bottom wall. New flow states are investigated based on the fully three-dimensional solution of the Navier-Stokes equations for these free-surface cylindrical swirling flows, without resorting to any symmetry properties unlike all other results available in the literature. To our knowledge, this study delivers the most general available results for this free-surface problem due to its original mathematical treatment. This second part of the dissertation is a basic research task directed at increasing our understanding of the influence of the presence of a free surface on the intricate transitional flow dynamics of shear-driven flows.

EPFL2007

One dimensional and multiscale models for blood flow circulation

The aim of this work is to provide mathematically sound and computationally effective tools for the numerical simulation of the interaction between fluid and structures as occurring, for instance, in the simulation of the human cardiovascular system. This problem is global, in the sense that local changes can modify the solution far away. From the point of view of computing and modelling this calls for the use of multiscale methods, where simplified models are used to treat the global problem leaving to more accurate models the local description. Moreover it is characterised by the appearance of pressure waves inside the vessels. In large arteries the vessel wall dynamics can be described by a thin elastic membrane model (Navier equation) while the fluid motion can be represented by the Navier-Stokes equations for incompressible Newtonian fluids. Unfortunately, given the high levels of details furnished by this model, its computational complexity is dramatically high. Therefore reduced models have been developed. In particular, one-dimensional models, originally introduced by Euler, seem to be appropriate for the study in the time-space domain of pressure wave propagation induced by the interaction between the fluid and the vessel wall in the arterial tree. These reduced models are obtained after integrating the Navier-Stokes equations over a vessel section, supposed to be circular, and assuming an algebraic wall law to describe the relationship between pressure and wall deformation. They can be used in place of the more complex three dimensional fluid-structure models or in cooperation with them (multiscale approach). The first part of this work deals with one dimensional models. A reduced 1D model taken from literature is presented and analysed. Some extensions of the basic model, in particular with respect to vessel wall law (generalised string model) and more complex geometries (bifurcated and curved arteries), are also considered. Numerical schemes are proposed and some numerical results are presented. In the second part of this thesis we focus on a multiscale model. We consider a 55 arterial tree, described by the 1D model, coupled with lumped parameter models for heart and capillaries. In particular, specific attention has been devoted to the coupling between the left ventricle and the arterial system, whose physiopathological relevance is well known. This mathematical model gives good results in numerical tests and is able to describe the relevant features of the pressure wave propagation and reflections within the arterial system.

EPFL2004

One dimensional and multiscale models for blood flow circulation

EPFL2004

Simulation of shear-driven flows

Roland Bouffanais

EPFL2007