^{1}

^{2}

In this work, we first investigated the hemodynamic parameters in the case of a normal aortic arch anatomy and in the case of aortic coarctation anatomy, both generated by using CFX-ICEM-ANSYS simulations. Then, we compared these results with those obtained for a proposed model without and with aortic coarctation, while introducing a real tridimensional magnetic resonance imaging geometry in the simulation process. The conclusion is that our proposed model reproduces, with a high agreement, the real case obtained from imaging data.

One of the main goals for numerical investigations of the blood flow in the thoracic aorta, and in particular in the aortic arch, is to well understand changes flow patterns and associated pressures in the presence of aortic pathologies, such as aortic coarctation or atherosclerotic diseases. Indeed, aortic coarctation induces an internal geometry deformation of the aortic wall, which induces at its turn significant changes in the hemodynamic characteristics of the flow. These phenomena would affect the normal blood circulation, and strong turbulence would occur, and then affect the aortic wall, leading ultimately to the need of surgical intervention to avoid fatal events.

More precisely, the aortic coarctation is a narrowing in the aorta (

This aortic disease appears as a good physiological and physical problem which needs to be treated theoretically using the fluid mechanics tools, experimentally by using medical imaging equipments and numerically by using computational simulations performed by means of specific software of fluid dynamics.

In literature, we found many experimental and computational investigations of blood flow in the human aortic arch, [

In order to enhance the performance of these models and to reflect the real situation, they have to correctly take into account several features like the non-Newtonian behavior of the blood, the unsteady and pulsatile aspects of the flow and the distension of the aortic wall. Hence, the role of the non-Newtonian behavior of the blood, which may be predominant under certain flow conditions [

The main goal of this work is double. Firstly, we were interested in studying the simulation of the blood flow in thoracic aorta, by using some imported recent tridimensional magnetic resonance imaging data of a real healthy aorta and a real aorta in presence of coarctation, provided by “Laboratoire d’Imagerie Médicale-Un- iversité de Pierre et Marie Curie” (LIB-UPMC), in the context of a specific CFD software (ANSYS-ICEM-CFX), to build a meshing and elaborate an adapted geometry, and to examine the obtained results of the realized computational simulation. Secondly, we compared these results with those obtained from a computational simulation using the same software on a geometrical model of a healthy aorta proposed in [

Firstly, the blood is considered as an incompressible Newtonian fluid, which is generally accepted as a reasonable first approximation to the actual non-Newtonian behavior of blood at shear rates observed in large arteries [

Furthermore, the blood fluid is assumed to have a density of 1.0595 gr/cm^{3}, dynamical viscosity of 3.45 cp, and molar mass of 3527.89 g/mol.

In all our simulations, we apply the k-є turbulence model with a Reynolds number Re = 13,148 for the first case (proposed aorta model without coarctation), Re = 12,835 for the second case (proposed aorta model with coarctation), Re = 13,810 in the third case (real aorta model without coarctation), and Re = 27,046 for the last case (real aorta model with coarctation). In order to ensure the convergence of each time step, we set the convergence criterion for the relative residual for all dependent variables equal to 1 × 10^{−}^{4}. We also adopt the non-buoyant model, with subsonic flow regimen, normal speed of 150 cm/s, pressure profile bend of 0.05, relative pressure of 60 mmHg and reference pressure of 760 mmHg (i.e. 1 atm). Our computational simulation results using ANSYS-ICEM-CFX 13.0 concerns the profiles of pressure, velocity and WSS in the four cases.

Furthermore, the MRI technique allows measurement the blood flow. At the cardiovascular level, it permits to acquire not only dynamical anatomic images of heart and large vessels, but also functional images of blood flow by means of phase contrast. Applied to the thoracic aorta, this technique constitutes a reliable analysis method. Because of its noninvasive nature, it is the reference for the evolutionary track. For the diagnosis of diseases of the aortic wall, its use varies among centers. In fact, Velocity ENCoded gradient echo imaging (VENC), also known as phase contrast imaging, is an MR technique for quantifying blood flow velocities and volumes. Measuring the phase shift that occurs as protons in the blood move through a magnetic field, the velocity and direction of the blood can be obtained. In order to obtain accurate measurements from a phase contrast image, a plane must be selected perpendicular to the path of the flowing blood. The technique creates a velocity encoded image for multiple phases of the heart cycle, thus creating a cine dataset. In the cine produced, stationary tissue appears gray with tissue moving through the plane appearing as shades of either white or black, depending on the direction. The bright or dark the tissue is, the faster it is moving through the imaging plane. Quantification of blood flow is performed by software that requires the user to outline the vessel of interest on a single phase of the cardiac cycle and the remaining phases are automatically processed. The software produces thus time-veloc- ity and time-flow curves. This technique is valuable for evaluating stenosis, in our case the coarctations. A pressure gradient across a lesion can be calculated from blood velocities. Blood flow is quantified from the individual blood velocities at each point over the cross section of the blood vessel [

In this work, we used the LIB-UPMC imaging data concerning a healthy aorta and an aorta with a coarctation as well as data concerning the blood flow (aorta section, velocity and flow-rate). Using these imaging data, we constructed geometries of these two real aortas and performed on them the same computational simulation as used on our models, in order to compare results. We use also the blood flow imaging data to characterize aorta section, velocity and flow rate as well as the pressure field.

Firstly, in the proposed model, [

In this first case, the computational simulation gives the following results (Figures 3-5).

Concerning, our model with coarctation, [

In this second case, the computational simulation gives the following results (Figures 7-9).

Importing imaging data concerning a real healthy human aorta, we have constructed its geometry using ICEM-CFX code, and we obtained the following meshing (148,319 nodes and 870,374 tetrahedral elements) and geometry (

In this third case, the computational simulation gives the following results (Figures 11-13).

Importing imaging data concerning a real human aorta with coarctation, we have constructed its geometry using ICEM-CFX code, and we obtained the following meshing (154,331 nodes and 901,105 tetrahedral elements) and

geometry (

In this fourth case, the computational simulation gives the following results (Figures 15-17).

Comparing the results obtained for the case of a healthy aorta, the computational simulation gives practically the same profiles of pressure velocity and WSS for the real geometry and the proposed model. This permits us to say that our model without coarctation is a good model.

In the case of an aorta with a coarctation, the confrontation between the results obtained from the computational simulation for the case of the real model and the proposed model shows that there is a good agreement between them for velocity profiles as well as the minimal values of pressure and WSS. Nevertheless, there is a

notable difference between the maximal value of pressure and WSS between the two cases. In fact, the maximal value of pressure in the case of real model is equal to 21,250 Pa, nearly the double of that for the case of proposed model (10,830 Pa), and the maximal value of WSS in the case of the proposed model is equal to 37.08 Pa whereas it is equal to 109.2 Pa in the case of real model. These disparities may be due essentially to the differences between the geometrical sizes of the real model and the proposed model and certainly also to the difference between the degree of reduction at the level of the coarctation in the two cases.

Despite these differences, our computational simulation has however allowed the partial validation of our model of aorta without and with coarctation. It remains to improve slightly this model in order to approach the most possible the real case.

We plan to treat the generalization of this work for the case of a pulsatile unsteady blood flow with taking into account the elasticity of the walls. In this case, the obtained results from computational simulation will be compared with the real measurements obtained in vivo by means of medical imaging tools.

As a continuity of this work, we plan also to treat in a future work the blood flow imaging data characterizing the aorta section, the velocity and the flow rate as well as the pressure field.

Furthermore, it is axiomatic that humans would have aortic arches with different sizes and with slightly different forms, depending on age and gender [

The authors are very grateful to H. Benali (Head of LIB-UPMC-Paris) and especially to N. Kachenoura (Head of cardiovascular imaging group at LIB-UPMC-Paris) for their crucial help, useful discussions and for providing them recent tridimensional magnetic resonance imaging data of a real healthy aorta and a real aorta in presence of coarctation, during their scientific visit at LIB-UPMC-Paris.

Fatima Moumen,Abed-El-Farid Djemaï, (2016) Computational Fluid Dynamics Analysis of the Aortic Coarctation. Natural Science,08,271-283. doi: 10.4236/ns.2016.86032