Medical Decision Making for Cardiac MRI with CFD “Detection of Severe Stenosis Using a 5D Model of the Descending Aorta”

The aim of this study is to develop a reliable 5D (x, y, z, time, flow dimension) model for medical decision making. Sophisticated techniques for the assessment of serious stenosis were developed using time-dependent instantaneous pressure gradients through the aorta (flow rate, Reynolds number, velocity, etc.). A 74 cardiac MRI scan and 3057 scans were performed on a 10-year-old patient with congenital valve and valvular aortic stenosis on sensitive MRI and coarctation (operated and then dilated) in the sense of shone syndrome. The occlusion rate was estimated to be 80.5%. The stenosis area was approximately 15 mm long and 10 mm high. The fluid solver (NS) exhibited a significant shear stress of −3.735 × 10−5 Pa within the first 10 iterations. There was a significant drop in the flux mass of −0.0050 (kg/s), as well as high blood turbulence in vortex field lines and low geometry Reynolds cells. The fifth dimension was used for negative velocity prediction (−81.4 cm/s). The discoveries of the 5D aortic simulation are convincing based on the evaluation of its physical and biomedical features.


Introduction
Visualization is crucial for the display and understanding of medical image data. For diagnostic and surgical planning, radiologists and surgeons must be able to evaluate the data appropriately. Many imaging systems' data can incorporate both functional and structural information, resulting in 4D datasets. When the image contains spectral information, it can be extended to 5D in some circumstances. Overall, 5D imaging reveals more information than 4D imaging. However, there are various ways to visualize 4D medical data, visualizing 5D medical data. The inability to properly portray 5D medical datasets on a 2D screen has drawbacks. Currently, there are five generations of visualization techniques [1] in medicine: 1D waveform display, 2D sequence display, 3D dataset visualization, multidimensional dataset visualization, and virtual reality visualization. Dynamic volume datasets, such as functional magnetic resonance imaging (fMRI) and positron emission tomography (PET) images, are multidimensional medical images that incorporate four dimensions: three for space and one for time. Five-dimensional medical images were created as a result of advancements in medical imaging methods; incorporating spectral information into an image collection is a good example. The representation of 5D medical data on a 2D screen poses significant hurdles. This study [2] used 5D brain electrical impedance tomography (EIT) images as a case study. However, the approach presented in this study can be applied to different imaging technologies. EIT is a relatively recent approach to medical imaging. EIT is based on the fact that different tissues have varied impedances. By injecting a tiny current via sensors encircling the region to be photographed, EIT imaging uses this attribute (i.e., impedance). EIT imaging is a low-cost, safe, and portable option. EIT offers great temporal resolution, but poor spatial resolution In this study, we developed the concept and modeling of the 5D cardiac system (3D model + temporal dimension functional dimension of the flow) detailed in our research work mentioned in [13,14]. This strategy consists of reconstructing a 3D geometry of the descending aorta and diffusing a 2D viscous laminar fluid to detect the zones of narrowing. Blood flow is obstructed at the location of the stenosis in the restricted aorta, restricting the movement of blood cells and creating turbulence in the internal aortic wall. However, blood movement is affected by several other factors, including the thermophysical properties, including viscosity, surface tension, and wettability [15]. Over the last few decades, numerous statistical models of blood rheology have been researched, and various secret features, such as fluid behavior and flow behavior, have been uncovered as a result of these experiments [16]. Consequently, assessing the significance of wall pressure and wall strain stress in arteries is crucial for medical researchers. Pulsed flow behavior was documented by several scholars, and its variations are constantly damped, which can be due to the elasticity of blood vessels. Hemodynamics are affected by blood vessel disease because they disrupt the flow cycle, resulting in a reduction in the wall pressure and shear tension in the arteries [17]. The biochemical mechanisms of blood-related illnesses are the subject of several studies. However, understanding the fundamental physics of the disease in order to understand the mechanism and hence pave the way for less invasive and more sustainable strategies for their prevention is an enormous development. The parietal strain and parietal tension exerted by blood on the inner periphery of the artery are not recognized by doctors using current imaging techniques [18]. The breakthrough in anatomy and cardiovascular physiology physiological simulation has created an opportunity to bridge this gap. Recent advances in the field of computational fluid dynamics (CFD) have made it possible to model blood flow in heart structure geometry. It is a less invasive method, and a device may interpret the blood flow pattern of a disease-related artery. As a consequence, in the areas of congenital heart valve, coronary, myocardial, and peripheral artery disorders, CFD has become a clinical testing instrument for medical practitioners. The CFD module, on either hand, includes non-clinical experience, technical software, efficient computer systems, and a vast number of calculations because it depends on the exact requirements of geometric and flow boundary conditions [19]. It is necessary to implement CFDs in standard clinical procedures, owing to these conditions and restrictions. Viscous dissipation can also be measured using the viscous term of the Navier-Stokes equation, which removes the need for friction and relies solely on internal blood flow velocities, which can be measured noninvasively using 4D flow MRI [20].
Our main contribution is to present a numerical simulation of laminar blood flow in 3D aortic modeling in the presence of a left subclavian aortic coarctation, and an analytical study is conducted to study the impact of a solver of dynamic fluid on the detection of aortic stenosis.

Dynamics of Laminar Viscous: Navier-Stokes Equations
Fluid flow simulations are based on Newtonian and fluid property physics. Physical rules, such as the maintenance of kinetic energy, which contribute to the equation of motion, must be satisfied by these features. The viscous stress tensor is linearly related to the rate of the strain tensor in a Newtonian fluid. Assume Stokes flows (low Reynolds number) and only slight spatial variations in hydrostatic pressure [21].
Let Ω ⊂ R2 be an open domain bounded and connected to the Lipschitz boundary Γ. Consider the Navier-Stokes equation With the initial condition v (0) = v0, where v is the velocity, p is the pressure, ν is the viscosity (or the inverse of the Reynolds number, i.e., ν = 1/Re), it is an external timedependent body force when the Reynolds number approaches a critical value or minor fluctuations are implemented into the flow. It is well understood that the flow transitions from laminar to turbulent [22].
The fluid was intended to be incompressible in this simulation, which meant that the density had to be steady and Newtonian. The vessel wall had to be impervious to deformation. A volumetric flow waveform was defined at the inlet of the pulsatile flow. For the entry of this model, a fully defined laminar velocity profile was used. The blood fluid was set to have a density of 1056 (kg/m 3 ) and a viscosity of 0.06 (Pa). These fluid properties set the Womersley parameter (α), which measures the frequency of the pulsation, defined as [23]: where ω is the angular frequency, and υ is the kinematic viscosity of the fluid. Properties of the fluids are mentioned for this model as well as the boundary conditions. In laminar flow regimes, the related term of the Navier-Stokes equation can be used to measure the volume of the viscous dissipation process [20] where Φ vD is the viscous dissipation per unit volume based on the viscous dissipation term, and µ is the dynamic viscosity. δij = 1 for i = j and δij = 0 for i = j, where i and j are the principal directions (x, y, z) in Equation (4), which consists of the elastic viscosity of the velocity field and the spatial derivatives. If the speed field is defined, such as from 4D flux MRI measurements, Equation (4) can be used to quantify the viscous dissipation per unit volume. The integral of the viscous dissipation of the unit was used to measure the overall viscous dissipation.

Boundary Conditions
The condition was extended to the fluid model, and the vessel walls were considered to have a slip-resistant boundary. The porosity of vessels is often underestimated. As a result, the fluid is not supposed to move through the vessel walls as it passes through it. The boundary conditions applied to the vessel walls are [23]: There must be no increase in velocity in the steering radial along the line's equatorial plane, essentially reducing the radial portion of the velocity to zero. Consequently, the applied boundary condition is The incompressible Newton flow rate through the established geometry is depicted by a parabolic velocity profile. As the Dirichlet boundary state at the output, a variable and spatially uniform pressure boundary condition is applied, and for the limit condition, the velocity Neumann (zero gradient of velocity in the axial direction at the output) is applied [24].
The wall of the aorta was rigid. Consequently, the walls of the flow domains were static and stationary. Boundary conditions for speed were imposed on the wall with no slip or flow. The non-flow condition means that the velocity in the normal direction of the wall is zero because the walls are impervious. The pressure gradient perpendicular to the wall was considered negligible using the Navier-Stokes equations [25].
Blood is a vigorous medium because its viscosity varies as a feature of the shear deformation rate, making it behave like a non-Newtonian fluid. Blood is biologically smoother in the systolic crest than when it moves slowly, as in the diastolic crest. This occurs as red blood cells clump together. The impact of thinning blood shear is as follows [16,23]: In this equation, µ 0 is the viscosity of blood at zero strain rates, and µ ∞ is the viscosity of the blood at an infinite or very high strain rate. The constants m, n, and a were experimentally determined.

Characteristics of Blood Flow Fluids
The sensitivity of a fluid to deformation under shear stress is known as viscosity. This is a type of fluid "friction", which explains the internal resistance of the fluid to the flow. The viscosity of a fluid is strongly influenced by the bonds between the molecules. The viscosity is mathematically defined as the ratio of the shear stress to the velocity gradient. The majority of fluids are Newtonian fluids, which have a steady viscosity. Plasma, blood cells, and other substances carried in the blood make up the blood. The number of particles in the plasma induces non-Newtonian activity in the blood, which means that the viscosity varies with the flow shear rate. The blood flow exhibits Newtonian flow behavior when the shear rate is sufficiently high. In normal situations, however, it is not possible to disregard the fluid's non-Newtonian behavior [26].

Geometry Design
In the "no fall" condition placed at the boundary, fluid in contact with the vascular wall can pass at the same pace as the wall. This is because the shear forces imposed by the wall on the fluid will ultimately cause the boundary flow to have the same velocity as the downstream surface, which is a rational statement in this analysis. Meanwhile, to preserve mass equilibrium, the flow away from the vessel wall continues to speed up until it reaches an established profile, where no further velocity profile changes occur in the flow direction. This causes a non-zero speed gradient, ∂u ∂y , where u and y are the velocity component in the direction of flow and the coordinate of the space perpendicular to the direction of flow, as shown in Figure 1. Mathematical relationship between shear stress of the wall, τ, and viscosity, µ blood, as explained in Section 2.3.1, is indicated in Equation (8) [23]: It was shown that achieving this quantity using experimental methods and relying on the experience of the fluid's viscosity and velocity profile near the vessel wall is challenging. This sum is easier to approximate using CFD techniques, but it is based on mesh consistency.
The governing equations for simulating the hydrodynamic flow of blood through a stenosing aorta with sufficient boundary conditions are as follows [16]: Laminar pulsed flows describe the boundary conditions at the entrance of the stenosing aorta. Because all of the geometries' inputs are circular, a speed limit is enforced to introduce spatial and temporal differences in the pulsed flow. Patients with aortic stenosis are most likely to have valvular stenosis. This study aimed to determine the effect of physiological conditions on the rupture wall of the stenosing artery.

Reynold Number
Laminar and turbulent flows are the two main types of flows that arise. Turbulence is a measure of the degree of oscillation of the surrounding fluid in the direction normal to the fluid flow, and it can be made on the basis of both the energy (kinetic energy of turbulence) and energy dissipation (dissipation turbulence rate) included in this motion. There is also little distinction between laminar and turbulent flows in terms of fluid particles. In chaotic flows, the particles appear to behave spontaneously, which is similar to the initial flow direction on average. The average fluid velocity, density, viscosity, and vessel diameter all play a role in determining when the fluids become turbulent (for internal flows). The Reynolds number Re is the name given to this characteristic value, as shown in Equation (12): where U is the average velocity of the fluid, D is the diameter of the vessel, ρ is the density of the fluid, µ is the fluid viscosity (dynamic), and ν is the kinematic viscosity of the fluid. The fluid tends to become turbulent at a Reynolds number of approximately 2000 for internal flows [27]. The laminar and turbulent flows have different profiles, as shown in Figure 2. It is extremely unlikely that blood can enter turbulent flow under natural circulation conditions. Only severely stenosed vessels, extremely abnormal flow paths, and even certain pathological disorders that impair blood flow properties cause this. Non-Newtonian fluids have a flatter profile toward the axes, so the formed profiles are different. Blood appears to have a uniform velocity profile rather than a developed profile because of the diverse flow paths and irregular sizes of arteries [28,29].

Assessment of the Rate of Severe Stenosis
The amount of cross section that is "blocked" by occlusion determines the degree of stenosis. The narrowest segment (called "groove") of the blockage is usually considered the cross section for stenoses with variable geometry. Figure 3 depicts a common model of a tube with a single stenosis in which the original radius of the vessel R0 is obstructed by a stenosis with a normalized maximum height of value = h/R 0 , where h is the maximum height of the stenosis. To obtain the unit length, the available radius, R 0 = 1. Because the standardized height of the stenosis must be 0 ≤ δ <1, we can deduce the following [23]. The normalized cross section of light through CD = π The normalized cross-section of light on AB = π (1 − δ) 2 Thus, the percentage of blocked area can be derived; therefore, the severity of the stenosis is obtained as follows:

Design and Description of the Solution Methodology
In this section, we apply our solution's architecture and general technique, which includes the use of the Navier-Stokes equation in 2D to determine the action of blood fluid in a 3D model. Figure 4 shows the steps of each component of the solution, as well as an overview of the advantages of using CFD modules and the methods used. Our benchmark is composed of three parts, as indicated in Figure 4: 1-3D Medical Image Processing Software and Mimics Innovation Suite and Materialize for export, import, and management of 3D modeling of the stenosing aorta format; 2-ANSYS Fluent for computational fluid dynamics (CFD) manipulation and fifth dimensional "flux" modeling; 3-MR 4D Flow data for the visualization of aortic blood flow as streamlines, pathlines, and color-coded vectors with Pie Medical Imaging stands for expertise in cardiovascular quantitative analysis software (confirmed to have stenosis by human experts).
The main purpose of CFD is to solve the problem by discretizing the paradigm into small cells. The CFD process is divided into three stages. The first phase is pre-processed. The domain of the study is defined and generated during this step, and the flux regions are discretized to form a mesh of cells. Better findings would be obtained with a finer mesh, but this would require more processing capacity. Several codes now have methods for adapting meshes after multiple iterations to obtain more precise meshes, if desired. The properties of the fluid, as well as any other phenomenon to be investigated, are evaluated at this level. It is essential to develop the required boundary conditions for the model until the medium begins, and the equations that need to be solved are determined. This is achieved by adding known or controllable values to the nodes or cells along the domain boundary. These limit values can be used to evaluate the solution for the remainder of the fluid velocity. It is projected that defining the domain geometry and creating the grid would take up more than half of the time spent on a CFD project. Although current software has made this move simpler, rational skills and knowledge are needed to build real-world meshes and provide correct answers without using more computing power than is appropriate. In the second step, the mesh generated in the previous step is solved. Assigning an initial value to cells first initializes the model and performs preliminary assessments of solution techniques. The final stage is post-processing, which involves extracting and analyzing the data produced in the previous step. Rather than exporting raw data, today's commercially available applications aim to improve the accessibility of the findings by optimizing visual appearance and data manipulation. This provides a graphical user interface that can depict geometry and grid, vector displacements, contour and surface tracing, particle detection, and other features based on the needs of the user. This move often requires fundamental knowledge to understand these findings.

Tools for 5D Imaging
To meet the increasing need for CFD in industrial applications, several businesses have created a variety of applications. To satisfy the existing demands of search and processing resources, the simulation capabilities as well as the user interface have been greatly enhanced. ANSYS (https://www.ansys.com/products/fluids/ansys-fluent, accessed on 21 November 2021 (Academic and Professional version of Enterprise)), which is still one of the most widely used commercial codes for structural research, has also played a key role in the establishment of a package in the market. The benefit of this program is that it allows multiple supported programs to run on a single platform and be connected to one another. This was initially designed to provide developers with a single board on which they could construct a product from the ground up, validate it, and refine it. This includes constructing the target model, discretizing it into finite elements, running the requisite models, supplying the required data, and optimizing the model based on the available data. The number of applications capable of running in this area grew in tandem with the market demand. Finally, in addition to basic structural models, a wide variety of mechanics are included, including fluid flow, electromagnetic, and thermal [30][31][32].
We used the product of Pie Medical Imaging "Caas MR 4D Flow" (https://www.piemed icalimaging.com/product/mr-solutions/caas-mr-4d-flow/, accessed on 21 November 2021 (Private access to the user with a procuration of agreement for the use of the product)) to analyze the valve behavior, considering the difficulty of the reconstruction of the valve structure from the flow sections. This cardiovascular solution aims to analyze blood flow in cardiovascular structures. Blood flow can be evaluated using three-way phase-contrast MRI images and the corresponding MRI images. In addition, blood flow can be analyzed by retrospective reformatting of flow planes. Caas MR 4D Flow is software that allows a cardiologist or radiologist to imagine and analyze blood flow in cardiovascular systems using multi-slice and multi-MRI images. It uses MRI and speed-coded images to examine blood flow in the heart and large arteries, and it contains the following features: analysis of quantitative cardiovascular results, segmentation of cardiovascular structures, and simulation of the intensity and direction of blood flow.
Whenever Caas MR 4D Flow findings are used in a clinical setting to sustain diagnosis, they should not be considered as the only compelling basis for clinical decision making.

Methods and Materials
For a 10-year-old patient with a history of congenital valve and valvular aortic stenosis on close MRI and coarctation (operated and then dilated) in the sense of shone syndrome, 74 cardiac MRI scans and 3057 images were obtained. Technically, reconstruction and segmentation of the descending aorta were performed using the following: 44 TRICKS angiographic slices in dynamic acquisition on the thoracic aorta; Injection perfusion sequences after injection; Sequences ciné-fiesta T2 short-axis 4 cavities; Subsequent infusion sequences short-axis.

Clinical Diagnostic Report of LV:
For the ascending thoracic aorta without parietal abnormality, its main measurements are as follows: Tricuspid valve. Aortic ring with 8.5 mm diameter; Aortic stenosis at 0.42 cm 2 with reduction of sigmoid opening at 5 mm; At the sino-tubular junction: 25 mm; 1/3 medium of the ascending aorta: 18 mm; Horizontal aorta: 16 mm; Size disparity with aortic stenosis at the isthmic level extended over 10 mm, reducing approximately 65% of its lumen by 6 mm in diameter; Mitral valve of normal diameter 2.5-4.3 cm 2 ; The mass of VG tele-diastolic 90 g and tele-diastolic 70 g; The systolic ejection function, estimated according to the 82% contour method.

Features of TRICKS cuts and acquisition protocols:
Magnetic field: 1.

Model of the Geometry of the Descending Aorta in 3D
In this section, we proceed to the first phase of our 5D model, which consists of defining the geometry of the descending aorta in 3D. This step states the reconstruction of the aortic model from the TRICKS sections, as shown in Figure 5. The multiplanar reconstruction in Figure 6 is a crucial step in the 3D reconstruction of the descending aorta presented in the axial plane in Figure 7.  We applied a threshold to the areas of interest to define the active mask. It was applied based on a low limit and a high limit. The mask contained pixels with a value between two threshold limits. The upper and lower thresholds were limited to the maximum and minimum intensities, respectively. We present the concept of the generic algorithm in Table 2 as follows.  Based on the detailed mathematical concept of [33], the thresholding is typically used in scans where a certain anatomy has a very distinct set of values (high contrast). Keeping only certain values of the medical image allows us to visualize the border of the aorta; in our case, the minimum and maximum threshold values were determined automatically. The threshold value for Mimi was approximately 270 (Hounsfield scale). A green mask was created after thresholding. The minimum threshold value for this case was 1235 (grey values) or 211 (Hounsfield values). Thresholding needs to be done before region growth since all previous work is lost after changing the threshold value.
The updates according to the defined limits give a reconstruction of the model of the aorta in 3D based on the multiplanar reconstruction of the TRICKS cuts in Figure 8.
The second step was to divide the single mask for the entire model into two separate masks. This tool allows easy and fast separation of anatomical parts, such as the anatomical model of the heart and aorta. Two masks with the region of interest must be created and used as inputs to divide the selected mask, as shown in Table 3.
We chose for the geometry calculation the "optimal" quality, partly corresponding to the input mask and custom parameters. The last step of this section was to clean the waste generated by the 3D reconstruction, as well as the smoothing and refining of the aorta model presented in Figure 9. The best is that the accuracy of the three-dimensional reconstructions obtained depends essentially on the conditions and quality of acquisition of cuts, thickness of the cuts, distance between two successive acquisition slices, and the quality of the segmentation and the reconstruction algorithm used. The final geometry is shown in Figure 10. Based on a few statistical measures in Table 4, it can be seen that the thickness analysis presented in Figure 11 of the aorta can identify the stenosing portion of the aorta.  Note that the local curvature, as shown in Figure 12, of the segment that contains the shrinkage takes negative values. At this stage, the surface deformation as well as the measure of local curvature can be useful indices for aid in diagnosis. Based on the study in Section 2.3.4, the height of the obstacle, h, is calculated in Table 5 with respect to the unobstructed radius, which gives δ in Equation (13) a value of the severity of the stenosis to each part of the segment of the aorta.

Generation of the Mesh Geometry
In this section, it is crucial to prepare a proper format for the aortic geometry to proceed to the mesh operation. The 3D model was transformed into a geometric smoothing (IGS) format to divide the entire volume into facets that contain triangles, edges, and intersection nodes. This IGS export approach is used to process surface-based functional data, primarily to recreate functions from noisy observations. The proposed IGS method has a broad range of applications, especially in industrial settings. Based on the simulations, the IGS appears to be similar to other commonly used CFD processes. The mathematical modeling and simulation environment presentation is detailed by [34] with precision analysis for IGS models described by [35]. The second phase consists of decomposing the IGS model into three elements, the body, inlet, and outlet, as shown in Figure 13.

Setting of the Solution
To solve the 2D Navier-Stokes equation, the configuration of this fluid solver is based on the pressure with absolute velocity formation in a stable time. Our fluid is blood with a density of 1056 (kg/m 3 ) and a viscosity of 6 (kg/m 3 ). The boundary condition must be zero for the aortic model. The execution of the laminar flow model defines the velocities and pressure in the viscous flow field (the Navier-Stokes equation) to solve the displacement of the internal facet of the 3D model of the aorta. During the resolution during the first 10 iterations, the velocity residuals along the three axes (x, y, z) are estimated in Figure 15. The dimensions of velocity and pressure are combined in the momentum that governs the fluid equations. The coupled and independent methods are the two primary methods for solving discrete time equation algebraic equations. The simultaneous solution of the velocity and pressure parameters characterizes a coupled system. However, it is not commonly used in technical problems, owing to its poor processing power and high memory requirements. Coupled methods are commonly used for calculating compressible flows, whereas separate methods are favored for calculating incompressible flows. Unlike a coupled solution, a single approach solves the velocity and pressure fields independently or sequentially. It has the benefit of reducing computer memory and processing time, making it more effective for studying incompressible environment simulation fluids, as in our case for aortic stenosis modeling [36]. The pressure analysis of the proposed model is shown in Figure 16.
From this illustration, it can be concluded that the pressure at the entrance and exit is negligible compared to the body outside and inside the model of the stenotic aorta. The pressure intersection (Pa) between these two zones (inside and outside) of the wall at a critical point at a position of 15 mm shows a remarkable overlap between the red and blue parts, which explains why turbulence can occur because of the strength of the parietal wall exerted by the external and internal body. The location of blood flow in a case of aortic dissection with a complicated geometric feature, qualitatively and quantitatively, based on the evolution of vortex structures and their interaction in the narrowing region throughout a cardiac cycle can provide an index of the presence of stenosis from a 15 mm length of the aortic segment. A vortex field simulation in the wall of the internal aorta is shown in Figure 17 as well as the flow rate mass in Figure 18 to evaluate the vortex behavior. We deduced from these two results that at the position of 15 mm, there is a very high vortex magnitude that reaches (800,000 (1/s)) as well as a mass of negative fluxes during the first six iterations during the solver calculation. This reflux shows that it occludes a segment of the aorta. These indices support the results reported in the clinical assessment with isthmic extensor aortic stenosis extending over 10 mm in length, reducing approximately 65% of its lumen by 6 mm in diameter, compared to what we estimated to be a 15 mm wide aortic stenosis; a 5 mm error rate was detected.
Owing to the very small size of the geometry, the Reynolds number is small enough for the laminar flow that appears throughout the aorta model, as shown in Figure 19. In our case of aortic coarctation, the highly disturbed flow areas exhibited Stokes flow characteristics. The friction factor, vortex length, shear stress, leakage flow strain, and turbulence results obtained with the low-Re turbulence model were compared with experimental data and findings achieved in terms of velocity profile, vortex length, shear stress, and turbulence.
Because of the unfavorable pressure gradient produced in the expanding segment of the stenosis tube, downstream flow separation was observed, and its size increased with the Reynolds number. The wall shear stress along the stenosis aorta peaked at the stenotic throat, and the tip was notably lower for a longer stenosis. The critical Reynolds number at which the blood flow becomes transient or turbulent distal stenosis was accurately predicted using this model. More surprisingly, the vortex length measured with the formula (Low-Re) nearly matched the vortex length expected by laminar flow modeling on the Re spectrum of the laminar flow. The constant flux was caused by the low Reynolds number.
The flow becomes unstable at a relatively high Reynolds number. Finally, the analysis clearly shows that the proposed model is appropriate for studying blood flow in specific occlusion zones [37].

Solution with ANSYS-Fluent
The average size of the stagnation and shunt zones was predicted by CFD, but the length of the current line and the variations in speed due to aorta flow were underestimated. However, for quantitative confirmation of CFD findings and the quest for flow effects, such as tortuosity and laminar flow behavior, the measurement accuracy must be increased. At almost the same time, the CFD simulation effects are represented not only in terms of pressure drop, but also in terms of vector orientation and velocity, current lines, and regions of blood flow stagnation [38]. This time-based flow concept is used in the 3D modeling of the stenosing aorta to generate the 5D format shown in Figure 20. Owing to the fiber of this model, the fluid flow is deflected, resulting in a sinuous orientation of the flow lines, which shows a three-dimensional flow around the wall of the aorta, mainly directed from the inlet to the outlet. The area of stenosis in the middle, indicated by lower speeds, is shown in Figure 21. Focusing on the velocity in the stenosing part, we noticed that the flux field lines were delayed and decreased from 6.034 to 1.207 ×10 1 to 1.810 × 10 1 ms −1 . Interestingly, the detection of negative pressure in the value narrowing zone was −3.735×10 5 [Pa], as shown in Figure 22.

Cardiovascular Solution with Pie Medical Imaging "Caas 4D Flow"
Our clinical case of aortic coarctation led to narrowing of the aortic valve with shone syndrome (two single-pillar valvular leaflet). In this case, the 5D model for the valve structure is corrosive, given the morphological dissection of the valve leaflet. The definition of the inputs and outputs, as well as the conditions of the limits of the aortic valve, are complicated because the valve acts as a blood pump. For this reason, the use of this research product that enhances the fifth dimension of 2D blood flow is used to predict and simulate the occlusion rate at the aortic valve level. Caas MR 4D Flow allows examining the 2D blood flow by reformatting planes in a 3D volume retrospectively. Additional plans can be placed on the 3D volume's middle side. Standard 2D flow parameters, such as the flow rate, forward flow, and back flow, can be measured using 2D flow analysis. In addition, an extreme eccentric flow measurement in a plane was calculated based on the method of moving the flow [39]. The contour is redefined on its initial contour, initially shown by the Caas MR 4D flow in Figure 23. The displacement of the flow is defined as the distance between the center of the light and the "center of velocity" of the flow, normalized with respect to the diameter of the light. The center of velocity (C vel ) was calculated as the average position of the pixels of the light (r i , where i = x, y, z), weighted by the velocity information (v i ) as follows: where i represents the pixels within the vessel's outline, and j represents the spatial orientation of the pixels relative to the location of the center of the vessel. This method is described in more detail in the literature [39,40]. The flow displacement was calculated for the maximum systolic phase only, as shown in Figure 24. The maximum systolic phase was determined for each transmitter plane and determined using the time slot at the top of the 2D graph. The maximum systolic phase was indicated as a marker of phase slider. A description of the measurements is presented in Table 6.  Antegrade flow is used to quantify the blood pumped in the positive direction measured on the plane in a cardiac cycle, while retrograde flow is used to measure the amount of blood pumped in the positive direction on the plane in a cardiac cycle. The pump function of the valve was measured using forward flow to the back. The regurgitation fraction is the ratio of the backward flow (mL) to the forward flow (mL) and the area under the curve of the negative (backward flow) portion of the cardiac cycle. The gap between the vessel's axis and the center of the eccentric flow was normalized to the total size of the vessel to enhance the heat transfer displacement [39]. It is possible to draw from this section the interest of studying the fifth dimension of blood flow to quantify and estimate both the regurgitation rate as well as a simulation of the minimum negative sign velocity of −81.4 (cm/s) to predict stenosis at the level of the aortic valve.

Discussion
A medical decision support method for cardiac imaging in MRI focused on 5D modeling (3D anatomical structure, temporal dimension, and blood flow dimension) for the study of a promising case of aortic coarctation with extreme valve narrowing was described. Our object of interest was a segment of the descending aorta for prediction, identifying the presence of aortic stenosis. The results show that in the first stage, the 3D modeling provides a very interesting index for the experts, which makes it possible to estimate the occlusion rate of 80.5% compared to what was manifested in the clinical assessment with a rate of 82%. The measurements extracted for the thickness and the local curvature with respect to the geometry of the aorta mark the zone of stenosis and the degree of deflection of the narrowing portion. In the second stage, the presentation of the 5D approach was performed through a combination of the 3D model and the size of the circulating blood flow as a function of time. The experiments that were performed for the fifth dimension provided high accuracy for the location of the stenosis zone of 15 mm in length compared to the clinical prognosis, indicating that aortic stenosis is extended over 10 mm, which is deduced from the solver of fluid (NS). During the first 10 iterations, a significant decrease in the flux mass was reported with −0.0050 (kg/s), as well as high blood turbulence in vortex field lines and low geometry Reynolds cells, which is based on the understanding and observation of a negative pressure value of −3.735 ×10 5 [Pa]. For negative velocity recognition (−81.4 cm/s), the fifth dimension was managed separately to assess velocity at the aortic valve with shone syndrome.

Conclusions
An approach for diagnosing medical decisions using cardiac imaging with MRI was described previously. This method necessitates the creation of a 5D model, which is composed of five dimensions: the anatomical structure of the heart in 3D, temporal dimensions, and a functional dimension of blood flow for the diagnosis of valve stenosis. A comparative study of practical technologies that leverage the fifth dimension of flow for the derivation of medical inference in clinical routine was also included. This contribution was studied for an aortic stenosis, and it comprised creating a 3D model and solving the Navier-Stokes equations for laminar and viscous blood fluid to arrive at the proposed 5D model (3D + time + flow). The region and degree of stenosis can be classified by extracting measures (vortex field, flow masses, static pressure, and Reynolds number) based on the fifth dimension. With the increasing need for high-resolution simulations, it is critical to investigate the cost and reaction time of digital solvers that could benefit from recent architectures, including multicore processors, in the future.