Research Article  Open Access
An Ultrasound Simulation Model for the Pulsatile Blood Flow Modulated by the Motion of Stenosed Vessel Wall
Abstract
This paper presents an ultrasound simulation model for pulsatile blood flow, modulated by the motion of a stenosed vessel wall. It aims at generating more realistic ultrasonic signals to provide an environment for evaluating ultrasound signal processing and imaging and a framework for investigating the behaviors of blood flow field modulated by wall motion. This model takes into account fluidstructure interaction, blood pulsatility, stenosis of the vessel, and arterial wall movement caused by surrounding tissue’s motion. The axial and radial velocity distributions of blood and the displacement of vessel wall are calculated by solving coupled NavierStokes and wall equations. With these obtained values, we made several different phantoms by treating blood and the vessel wall as a group of point scatterers. Then, ultrasound echoed signals from oscillating wall and blood in the axisymmetric stenoticcarotid arteries were computed by ultrasound simulation software, Field II. The results show better consistency with corresponding theoretical values and clinical data and reflect the influence of wall movement on the flow field. It can serve as an effective tool not only for investigating the behavior of blood flow field modulated by wall motion but also for quantitative or qualitative evaluation of new ultrasound imaging technology and estimation method of blood velocity.
1. Introduction
Atherosclerosis is one of the most common types of cardiovascular disease and is usually related to the existence of stenosis in an artery. The stenosis is caused by a plaque formed by the accumulation of lipid substances, cholesterol, cellular waste products, calcium, and fibrin in the inner lining of an artery. The plaque may partially or totally make arteries clog up, which finally may result in a heart attack and stroke [1–3]. In addition, it has been well accepted that once a mild stenosis is developed, the resulting flow disorders further influence the development of the disease and arterial deformability and change the regional blood rheology [4]. Thus, efficient imaging of the blood flowing through a stenosed vessel is necessary for the fundamental understanding, early diagnosis, and treatment of many cardiovascular diseases.
As a consequence of technological advances in ultrasound imaging, the ultrasonic echography has become one of the most important noninvasive diagnostic methods for detection and monitoring of cardiovascular disease. Compared with other available imaging methods, it is quick, safe, and relatively inexpensive when visualizing arterial wall and flow in vivo [5]. However, these ultrasound images, such as Bmode, Mmode Doppler, and color flow imaging (CFI) ultrasound images, cannot always accurately describe the real flow behaviors. The weak signals backscattered from the flowing blood, especially the signals echoed from the area near the vessel wall, will be contaminated by the strong backscattered clutter signals originating from stationary tissues, slowly moving vessels, and noise. However, it is believed that the echoed ultrasonic signals from slowly moving blood flow near the wall may provide useful information to improve the diagnostic efficiency.
Thus, many researchers have made remarkable efforts to develop new velocity estimation method, clutter filter, signal processing, and visualization techniques to obtain more accurate signals of blood flow. Therefore, the evaluation of these new techniques is important. It has been generally accepted that the validation using in vivo data is not desirable, because the signal components from the surrounding tissue and blood are not easily distinguishable, and the exact behaviors of the vessel wall are unknown [6]. Conversely, the computer simulation is a useful validation method because all parameters can be well defined.
Hence, some researchers constructed simulation models to imitate the behaviors of pulsatile blood flow based on the assumption that the wall of the artery is rigid [7–10]. Recently, taking into account the interaction between blood flow and the vessel wall, several investigators set up more realistic and complex models with the help of computational fluid dynamics (CFD) software [11–17]. They studied the effects of mechanical parameters of the elastic vessel such as stenosis severity, stenosis geometry, and length on the flow field. In these studies, models were built up considering wall moving only along with blood pulsation. However, it was reported that the motion of the vessel wall is influenced by not only blood pulsation, but also other physiological actions such as heartbeat, breathing, and body posture [6, 18]. These actions will change the movement of the vessel wall, which will interact with blood flow. The interaction is not yet clearly understood and will deteriorate the blood velocity estimation due to the strong echoes from the vessel wall.
More recently, a coupled method integrating the computational fluid dynamics (CFD) and ultrasound simulation was introduced in some studies [19–23]. Investigators modeled blood and vessel as a collection of point scatterers whose positions are updated using the interpolated velocity generated by CFD software. This coupled method provides a better understanding of the relationship between ultrasound images and the actual flow dynamics, as well as a better investigation of blood flow behaviors. However, for the method based on CFD software, setting and changing the boundary and initial conditions are inconvenient, and computing is timeconsuming, which restricts its scope of application.
The purpose of this study is to present a simulation model which can generate more realistic ultrasound echoed signals, considering the interaction between the vessel wall and blood flow. It simulates the behaviors of pulsatile blood flow, modulated by wall motion, passing through a stenosed common carotid artery (CCA). By combining numerical RungeKutta method and ultrasound simulation, it offers considerable flexibility in the ability to change mechanical and ultrasonic parameters, such as the stenosis shape and the characteristics of the ultrasonic probe. Furthermore, compared with methods based on CFD software, this model requires fewer computer resources and shorter computing time. Through this model, the successive Bmode, Mmode, and Doppler US images can be obtained efficiently. It would be a useful tool not only for investigating the behaviors of flow field modulated by the wall motion and the changes of ultrasonic images caused by the wall motion but also for quantitative or qualitative evaluation of new ultrasound imaging technology and the estimation method of blood velocity.
This paper is organized as follows. In Section 2, we firstly describe the general workflow of our study. And then, we present the geometric and mathematical model of pulsatile blood flowing through a stenosed vessel, considering fluidstructure interactions (FSI). The rest of this section elaborates on the processes of numerical RungeKutta simulation, coupling, and ultrasound simulation. In Section 3, this simulation model was applied to several imaging applications. The Bmode, Mmode, and Doppler US images were obtained. Then, to demonstrate the capability of this model, we made two types of comparisons: here one is among the results with different wall movements, for illustrating the model which can reflect the influence of wall movement on the ultrasound images; and the other is between the numerical and ultrasound simulation results, for demonstrating the similarity between the reference values and results of simulation. In Section 4, the discussion and conclusions of this work are presented.
2. Material and Methods
Figure 1 demonstrates the general workflow of our study. Firstly, we defined the geometric and physiological parameters of stenotic arteries in MATLAB, such as the stenosis length, the radius of the artery, the fundamental frequency of blood flow, and Young’s modulus. Then, based on these parameters, the fluidstructure interaction model, described by several secondorder ordinary differential equations (ODEs), was established. By solving the ODEs with RungeKutta method, we obtained numerical solutions of the biomechanical properties of blood flow passing through a constricted vessel, which will be used as input for succeeding ultrasound simulation. Then, the positions and amplitudes of the scatterers were updated, according to the displacements and locations calculated based on obtained blood velocity field and vessel radius changes. Through this combined approach, the acoustic analysis was performed by Field II (an ultrasonic field simulation software). Finally, the Bmode, Mmode, and CFI images were generated and compared with numerical simulation results. Compared with CFD coupled method mentioned, this method requires fewer computing resources and less computing time and is much more convenient for changing model parameters.
2.1. Geometrical Model
In most cases, the forms of stenosis can be approximately simplified as axisymmetric shapes, although the development of arterial structure does not follow any clearly defined geometrical pattern, and the geometry of the growth is changeable. As shown in Figure 2, the shape of the plaque, in our study, is assumed as the most frequently used cosine curves along the longitudinal direction, which is initially proposed by Young and Tsai [24]. The geometrical parameters of the CCA are created with a radius of 4 mm and an arterial wall thickness of 0.3 mm, which are close to physical values according to recent researches [25, 26]. The length of stenosis is set to be twice of the artery radius, which is corresponding to the average length of the carotid plaques [14]:
The geometry of the stenosis, whose mathematical description is given by (1) with the cylindrical polar coordinate system, is assumed to be axisymmetric. The maximum height of stenosis represents the severity of arterial obstruction. In our study, the motion of the arterial wall, assumed as a sinusoidal vibration, is considered. Accordingly, product term in this formula gives the mathematical expression of the motion of the stenotic arterial wall conducting sinusoidal oscillation at a certain frequency , while the amplitude of the oscillating arterial stenotic part is determined by a coefficient . To demonstrate the influence of wall movement, four different vibration frequencies and three different stenosis degrees (0%, 15%, and 20%) have been investigated, where is the fundamental frequency of blood flow. Hence, indicates that the wall oscillating frequency equates to the product of the fundamental frequency of blood flow multiplied by zero, and so on.
2.2. Mathematical Modeling
2.2.1. Fluid and Solid Models
In order to simplify the calculation, several relatively unimportant features such as nonNewtonian viscosity, slurry particles in the fluid, and temperature of biological flows may be neglected [27]. Thus, in this study, blood is assumed to be laminar, Newtonian, fully developed, viscous, and incompressible pulsatile flow [28]. For an arterial system, the wave speed of pulsatile blood flow is greater than flow speed, and the radius of the vascular tube is much smaller than the wavelength of the pulse [29]. Consequently, NavierStokes Equations can be linearized as follows [30]:
2.2.2. Wall Modeling
The arterial wall is assumed as isotropic and elastic material obeying Hooke’s law. Accordingly, the behavior of periodically pulsatile blood flow and the motion of the arterial wall are interactional and interdependent:
Noslip and nopenetration boundary conditions are supposed between the blood flow and the arterial wall. Therefore, the displacements of the fluid and solid part must be consistent. These coupling conditions are given in
If the highorder influence of the inertia items of wall and the viscous term of blood are ignored, the equation of wall motion at certain frequency can be approximated as follows [31]:
2.3. Solutions for the Equations
Because pulsatile blood flow is physically periodic, the parameters of blood flow, such as pressure, pressure gradient, and velocity, can be expressed by a sum of Fourier components in the following:where , , and are stationary items and , , and are th harmonic items of pulsatile blood flow. is the total number of Fourier decomposition items, which is set to fourteen here for a better approximation. The boundary conditions equations (4) are represented in the form of the following Fourier series:
After substituting (8)(9) into (2) with boundary conditions equations (9), the following set of equations are obtained [32]:where is the kinematic viscosity of blood, , , and and are the first kind Bessel functions with zeroorder and firstorder, respectively. Substituting (5) into (14), the secondorder variable coefficient differential equation can be obtained:
2.4. Numerical Method
In our study, the function ode45 in MATLAB is utilized to solve (15), which is the builtin solver for ODEs and implements a RungeKutta method with a variable time step for efficient computation. Firstly, (15) is reduced to a series of firstorder equations. After the defined integral span, precision and initial conditions are substituted into these firstorder ODEs, the numerical results of pressure harmonic components are obtained. Finally, by substituting into (14) and (15), the parameters, such as pressure, pressure gradient, and velocity of blood flow, can be calculated.
Figure 3(a) shows the centerline velocity waveform during one cardiac cycle in the stenotic arterial segment upstream at the position (), which is measured as the ensemble averaged mean velocity waveform from 50 realizations of a normal CCA (the carotid artery without stenosis) by using a portable Doppler ultrasound system (KJ2V2U, Nanjing KeJin Industrial Limited, Nanjing, China) and a 4MHz transducer [33]. The arterial pressure wave during one cardiac cycle, shown in Figure 3(b), is deduced from the distance between the heart and the CCA. Here, a dimensionless parameter is introduced to represent the horizontal direction, where is the specific time point and is the period of the cardiac cycle. As shown in both figures, the red line segments represent the systolic acceleration phase ( = 0.22–0.26). The parameters used in our computation are listed in Table 1.

(a) Blood velocity
(b) Pressure
2.5. Simulating Ultrasound Coupling Numerical Simulation
In this study, the ultrasound simulation is conducted with the Field II software which is created by Jensen [34, 35]. Based on linear systems theory, this software simulates ultrasound field using the spatial impulse response method proposed by Tupholme and Stepanishen [36, 37]. Since any transducer can be simulated by splitting the aperture into some kind of small shaped subapertures, and any transducer excitation and apodization can be included in the calculation, it is feasible to obtain realistic simulated ultrasound images [38].
In Field II program, the tissue is treated as a collection of randomly distributed point scatterers. Each ultrasound echoed beam, collected by the receiving transducer, can be calculated by summing the responses from these scatterers, and the scattering strength is determined by the density and speed of sound perturbations in the tissue [39]. As shown in Figure 1, the ultrasound simulation is performed based on the numerical outputs such as the velocity vector of blood flow, vessel radius, and other geometrical and physiological parameters. Through this combined approach, these scatterers in an ultrasound simulation move with velocities generated by the aforementioned numerical simulation. In this study, a full coupling method is used to capture the dynamic arterial flow and the estimated velocity field. Thus, the position of the scatterers can be updated at any moment over the cardiac cycle, and the dynamic ultrasound Bmode, Mmode, and CFI images can be achieved.
2.5.1. Setup Phantom
In this study, the phantoms are angled 45° with respect to the ultrasound beam which is in the axial direction. Considering the length of ultrasound simulation time, the scatterer density is set to ten scatterers per resolution cell, although it is related to the resolution of the images. According to its specific location (blood, arterial wall, and surrounding tissues) and simulation type (Bmode, Mmode, and CFI), these point scatterers are given different amplitude values whose distribution follows a Gaussian distribution [40]. The amplitude values used in the simulations are listed in Table 2.

To combine the numerical and ultrasound simulations, several works must be done. First of all, the temporal and spatial interpolations of outputs derived from the RungeKutta simulation are indispensable because of the difference of resolution. Secondly, the coordinate of every randomly distributed scatterer must be mapped onto the spatially interpolated coordinate of the obtained velocity vector. This operation will inevitably introduce approximate error. In this way, the radial and axial velocities of each scatterer in ultrasound simulation are obtained. Then, the new position of any scatterer can be acquired aswhere , and , are the radial and axial coordinates of the scatterer at current and previous time step, respectively, , are the axial and radial velocity of the mapped scatterer, and is the emitting interval of successive ultrasound beams.
For tissue motion simulation, these scatterers imitate wall and tissue movement only with radial velocity. The displacement within can be obtained from the interpolated output of the RungeKutta simulation. Consequently, the updated radius of scatterer in the wall or tissue is described aswhere is the radial scan area of the simulation.
Owing to the timevarying behavior of blood flow, files including the position and the amplitude of each scatterer point will be generated after each update. Thereby, a large number of files are involved here, and the exact number depends on the pulse repetition frequency.
2.5.2. Ultrasound Simulation
With the help of Field II software, the generated phantom files are then used to calculate RF lines individually for each different imaging direction and for each time step. Further, by representing these calculated amplitudes of echo beams with brightness and assembling these sequential signals, Bmode and Mmode images are obtained. To simulate velocity of blood flow or tissue motion, we applied an autocorrelation method, which was introduced by Kasai et al. [41] for ultrasound application. This approach computes mean axial velocity based on the phase shift of the center frequency of the RF signal. It is mathematically described aswhere is the wave velocity, is the central frequency, is the pulse repetition period, and are the real and imaginary parts of echo signal, respectively, and is the number of samples, which is related to accuracy. Thus, an estimated velocity in one direction is obtained. Doing this for all directions in an image yields a mapping of the velocity, displayed in color flow mapping (CFM) systems. The acoustic parameters of ultrasound simulation in Field II software are summarized in Table 3.

3. Results and Discussions
This section presents several applications of this model, which includes two parts. The first part illustrates the applications related to the movement of the vessel wall. In this part, Bmode and Mmode figures were generated and compared. Furthermore, we obtained the wall displacement with an image segmentation technique and compared it with reference values. The second part deals with applications related to the blood velocity estimation. In this part, CFI figures and velocity profiles for all cases were compared with each other and to reference numerical results. The simulations for all cases were performed in MATLAB 8.3 and on a workstation with two Intel Xeon 4 cores 64bit processors (1.8 GHz) and 16 G RAM.
3.1. Wall Movement Application
3.1.1. BMode Simulation
Figure 4 shows the Bmode images which are samples extracted from the simulated frames with different physical factor when the inlet velocity almost reaches the maximum (). Each row illustrates arterial stenosis with different wall oscillating frequencies and same narrow degree, while each column illustrates that with different stenosis degrees and the same frequency.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
It can be clearly observed in these subfigures that the arterial wall and stenosis parts are hyperechoic, demonstrated with a brighter color, while blood flow is hypoechogenic, illustrated with black. Thus, the narrow degree of the blood vessel can be detected from the hyperechoic area of these subfigures. Further, the differences can be easily found by comparing these subfigures with each other. As shown, we can learn that the position of vessel wall changes with the wall oscillating frequencies. However, the changes are imperceptible when the contracted degree is not serious, in spite of the different wall oscillating frequencies. Nevertheless, it gets more noticeable with the increase of narrow severity. For example, as shown in Figure 4(f) (with 25% stenosis, wall oscillating frequency), the hyperechoic area of the vessel is larger than that with other frequencies. Hence, these Bmode images can qualitatively demonstrate the movement of the vessel wall. Moreover, this model can generate successive Bmode images, which will dynamically show the wall movement.
3.1.2. MMode Simulation
As Figure 5 shows, the simulated Mmode images for different cases illustrate the movement of the vessel wall during one whole cardiac cycle. We just chose the echoed RF data at the axial central position of stenosis () to show the wall movement over time. These Mmode images clearly demonstrate that the wall vibrates almost in a sinusoidal pattern with a certain frequency, which coincides with the mathematical formula defined by (1). Moreover, it can be learned from these figures that the vibration amplitudes increase with the narrow degree. However, the vessel dilatation and constriction caused by pulsatile blood flow are hardly discerned in these figures.
A segmentation algorithm, speckle reducing anisotropic diffusion (SRAD) [42], was used to demonstrate the application of this model. Here, we only chose Figure 5(f) to show the result of SRAD process. As shown in Figure 6, the speckle noise is reduced, and wall boundaries are sharper. Then, the wall displacement can be easily acquired from these processed Mmode images, based on a simple thresholding method.
Figure 7 demonstrates the wall displacement from the RungeKutta numerical and the ultrasound simulation. The blue curves in these subfigures illustrate that the wall displacements got from the ultrasound simulation, which are the average values of 30 simulation results. All these curves were obtained from SRAD processed output with the same threshold. As shown, the ultrasound simulation results are in agreement with red curves, which depict the data got from the RungeKutta method. However, there is a small underestimation in some subfigures, such as in Figures 7(c) and 7(d). For further quantitative analysis, we calculate the mean and the maximum relative errors between the two types of wave forms and list them in Table 4. As shown, the mean error in all cases is minor, while the maximum relative error is 7.48% (stenosis degree = 25%, ).

3.2. Flow Velocity Simulation
As mentioned before, we can estimate the blood flow velocity from autocorrelation. The parameters of ultrasound simulation in Field II software are listed in Table 3. Figure 8 demonstrates color flow images (CFI) of the systolic acceleration phase (the red line segment shown in Figure 3, = 0.22–0.26) of blood flowing through the vessel with different stenosis degree and wall oscillating frequencies. Several features can be obtained by comparing these figures with each other. First, the maximal velocity of blood flow in the vessel with oscillating wall is greater than that in tubes without wall oscillating (the first column, ). Second, the maximum estimated velocity of blood decreases along with the increase of wall oscillating frequency. In particular, when the wall oscillating frequencies equate to the fundamental frequency of blood flow, these estimated maximum velocities (as shown in the second column of Figure 8) are higher than the maximum velocity with the same stenosis. Moreover, by comparing the rows, we can see that the increased value of the estimated velocity with a higher degree of stenosis is much more striking than that with lower degrees.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
As Figure 9 shows, the color flow image with data from the RungeKutta numerical simulation (a) is generated and compared with that obtained from ultrasound simulation (b). In spite of the fact that the simulated ultrasonic CFM image is more irregular to the velocity estimator statistics, it can be observed from these color flow images that the velocity field acquired from ultrasound simulation agrees well with the RungeKutta simulation results. For further quantitative analysis, the value of each point in the ultrasound simulated velocity field is compared with that of the corresponding reference velocity.
(a)
(b)
The velocity profiles obtained from the lateral center are illustrated in Figure 10, in which the blue curves depict the average velocity of 30 autocorrelation results while the red curves represent the reference velocity, and the light blue shaded areas in these figures represent the variance of obtained velocity. As can be observed, the ultrasound simulated velocity profile shows the almost same shape as the reference velocity. However, the difference between the two curves increases with the growth of the maximum velocity because of Nyquist’s limit. In general, the simulated velocity is underestimated, especially near the central area. Furthermore, because these ultrasonic CFM figures (Figure 8) are illustrated without filtering out lowfrequency and highamplitude clutters echoed by slowly moving tissues, the estimated blood flow velocities, especially higher velocities, are contaminated. For example, as Figures 10(d) and 10(e) show, the estimated velocities in the area near the wall significantly differ from the reference velocities. The presence of these discrepancies reflects the features of real signals in clinic, which is the reason why the clutter filter is indispensable. Meanwhile, it suggests that our model can be used to assess the clutter filter.
To further demonstrate capability of this model, we compared the simulated CFM figures with clinical data, as shown in Figure 11. The clinical CFM figures were obtained from the right CCA of a patient with a mild plaque. The angle between the probe and artery is about 120°, so the blood in Figure 11 is displayed in blue. Figures 11(a) and 11(c) illustrate CFM images in systole, and the middle parts in the vessels of both images are light blue which means the blood with higher velocity (the maximum blood velocity in Figure 11(c) is about 1.4 m/s.), while the blood velocity of outer parts is lower, represented in deeper blue. Figures 11(b) and 11(d) represent images in diastole. The blood in both figures is deep blue and there is a little change in color, which means the blood is stable and the velocity is lower. By comparison, the results illustrate that the simulated images show good consistency with clinical ones.
(a)
(b)
(c)
(d)
3.3. Discussions
With this model, as shown in Figures 4, 5, and 8, we can easily obtain Bmode, Mmode, and CFM ultrasound images. From these generated images, by comparing with different cases, it is clear that this model can qualitatively reflect the impact of wall motion on the blood field. Moreover, it can generate Bmode figures at any time, or successive Bmode and CFM figures during any period, as well as Mmode figures at any location. Accordingly, these figures can dynamically show the motions of blood flow and vessel wall, which is of more help in understanding the behaviors of blood flow modulated by wall motion. Meanwhile, the results shown in Figure 7, Table 4, and Figures 9 and 10 demonstrate that the ultrasound simulated results provide a good match with corresponding theoretical values. Furthermore, by comparison with clinical data, Figure 11 demonstrates that this model can simulate realistic ultrasonic figures. Therefore, we suggest that this model can be used as an analysis tool for quantitative evaluation of new ultrasound imaging technology and the estimation method of blood velocity, such as image segmentation or clutter filter.
In view of the fact that the object of this study is the stenosed CCA, the shape of the vessel in this study is 3D tubular instead of bifurcate or curved. Further, as mentioned, the main aim of this model is to generate ultrasound echoed signals, which can be used for investigating the changes of ultrasound images caused by wall motion and for evaluating ultrasound imaging and the estimation method of blood velocity. Thus, for this model, the most important requirement is the ability to accurately reflect the behaviors of blood flow modulated by wall motion. This is demonstrated by those abovementioned figures and comparison results. Moreover, with the tubular geometry, this model will not need to deal with complicated flows and wall motions caused by complex shapes, which will increase the complexity of ultrasound imaging and estimation method of blood velocity. Furthermore, the pipelike shape will simplify the simulation model and decrease the calculation time.
Based on the velocity field obtained from the RungeKutta numerical simulation, Field II software was used to generate Bmode, Mmode, and CFI images to visualize the motions of the blood flow and the vessel wall. As mentioned, for this fully coupled simulation method, the amplitude and position of every scatterer need to be updated after each RF calculating, which accordingly make the ultrasound simulation results show good agreement with the reference model. However, a great amount of computation time is undoubtedly involved, which is related to the scatterer density, the pulse repetition frequency, and scan line number. For example, the simulation time for any CFM case was almost four hours when fprf = 16 k and scan line = 20, while it was reduced to about half an hour when parallel computing technology was used on a multicore workstation.
To a certain extent, the implementation in our study is similar to the ones proposed by Balocco et al. and Swillens et al., in which the biomechanics and ultrasonic analysis are integrated based on US simulation software (Field II). The computational fluid dynamics (CFD) method was used in their studies to generate a new amplitude and position for each scatterer, while our model uses the RungeKutta method. Compared with CFDUS coupled method, this model needs fewer computer resources and less computing time and is much more convenient for changing model parameters. Furthermore, by taking the movement of tissue into account, this model can be used not only as an evaluation tool, but also as a simulation environment for investigating the interaction between the blood flow and arterial wall. For example, in Section 3, the comparison between results of all cases illustrates that the movement of the stenosed vessel wall will significantly affect the behavior of blood flow. It can be listed as follows. The vibrating motion of the arterial wall increases contraction and further increases the maximal velocity of blood flow. The closer to blood fundamental frequency the wall vibrating frequency is, the higher peak velocity becomes. The increased value of blood velocity with a higher stenosis degree is much striking than that with lower degrees.
In this study, the arterial wall motion is treated as a simple sinusoidal motion, while a more realistic motion should be taken into consideration in future research. For example, the wall motion can be treated as a complex movement influenced by breathing and the heart beating. Moreover, due to the assumption of taking arterial wall as isotropic and elastic material obeying Hooke’s law, the model cannot accurately depict the behaviors of all biological tissues. Hence, considering the anisotropy of vessel membrane and the viscoelasticity of the vessel wall may improve this model. However, although some simplifications have to be done in this study, this model is explicit on quantitatively assessing the ultrasound image changes induced by the changes of the vessel wall movement, such as the wall vibrating frequency, amplitude, and stenosis degree.
4. Conclusion
An ultrasound simulation model for the pulsatile blood flow, modulated by the motion of a stenosed vessel wall, is presented in this paper. Considering fluidstructure interaction, blood pulsatility, stenosis of the vessel, and arterial wall movement caused by surrounding tissue’s motion, it can generate more realistic ultrasonic signals. This proposed model couples ultrasound and numerical RungeKutta method and is characterized by its considerable flexibility in the ability to change mechanical and acoustic parameters. The simulated results show greater consistence with corresponding theoretical values and reflect the influence of wall movement on the flow field. It can serve as an effective tool not only for investigating the behavior of blood flow field modulated by wall motion but also for quantitative or qualitative evaluation of new ultrasound imaging technology and estimation method of blood velocity.
Nomenclature
:  The radius of stenosed artery 
:  The healthy artery radius 
:  The thickness of the healthy arterial wall 
:  The radial coordinate 
:  The axial coordinate 
:  The angular coordinate 
:  The maximum height of stenosis 
:  The oscillating amplitude of wall 
:  The oscillation frequency of stenosis 
:  The halflength of the stenosis 
:  The axial component of velocity 
:  The radial component of velocity 
:  The pressure of blood 
:  The blood density 
:  The blood viscosity 
:  The time 
:  The radial vascular displacements 
:  The axial vascular displacements 
:  The wall density 
:  Young’s modulus 
:  Passion’s ratio. 
Competing Interests
There is no conflict of interests.
Acknowledgments
This work was supported by the Grant (61261007 and 61561049) from the National Natural Science Foundation of China and the Grant (2013FA008) from the Key Project of Yunnan Natural Science Foundation.
References
 M. Texon, “A hemodynamic concept of atherosclerosis, with particular reference to coronary occlusion,” AMA Archives of Internal Medicine, vol. 99, no. 3, pp. 418–427, 1957. View at: Publisher Site  Google Scholar
 C. G. Caro, J. M. FitzGerald, and R. C. Schroter, “Atheroma and arterial wall shear observation, correlation and proposal of a shear dependent mass transfer mechanism for atherogenesis,” Proceedings of the Royal Society of London, Series B: Biological sciences, vol. 177, no. 46, pp. 109–159, 1971. View at: Publisher Site  Google Scholar
 D. L. Fry, “Acute vascular endothelial changes associated with increased blood velocity gradients,” Circulation Research, vol. 22, no. 2, pp. 165–197, 1968. View at: Publisher Site  Google Scholar
 D. Liepsch, “An introduction to biofluid mechanics—basic models and applications,” Journal of Biomechanics, vol. 35, no. 4, pp. 415–435, 2002. View at: Publisher Site  Google Scholar
 G. Von Reutern and H. Von Büdingen, Ultrasound Diagnosis of Cerebrovascular Disease, vol. 46, Angiología, Barcelona, Spain, 1994.
 M. Schlaikjer, S. TorpPedersen, and J. A. Jensen, “Simulation of RF data with tissue motion for optimizing stationary echo canceling filters,” Ultrasonics, vol. 41, no. 6, pp. 415–419, 2003. View at: Publisher Site  Google Scholar
 H. W. Hoogstraten, J. G. Kootstra, B. Hillen, J. K. B. Krijger, and P. J. W. Wensing, “Numerical simulation of blood flow in an artery with two successive bends,” Journal of Biomechanics, vol. 29, no. 8, pp. 1075–1083, 1996. View at: Publisher Site  Google Scholar
 H. Oung and F. Forsberg, “Doppler ultrasound simulation model for pulsatile flow with nonaxial components,” Ultrasonic Imaging, vol. 18, no. 3, pp. 157–172, 1996. View at: Publisher Site  Google Scholar
 R. Botnar, G. Rappitsch, M. Beat Scheidegger, D. Liepsch, K. Perktold, and P. Boesiger, “Hemodynamics in the carotid artery bifurcation: a comparison between numerical simulations and in vitro MRI measurements,” Journal of Biomechanics, vol. 33, no. 2, pp. 137–144, 2000. View at: Publisher Site  Google Scholar
 J. P. Ortiz, K. L. Bessa, D. F. Legendre, and R. H. A. Prado, “Physiological pulsatile waveform through axisymmetric stenosed arteries: numerical simulation,” in ABCM Symposium Series in Bioengineering, vol. 1, ABCM, Rio de Janeiro, Brazil, 2006. View at: Google Scholar
 B. Buriev, T. Kim, and T. Seo, “Fluidstructure interactions of physiological flow in stenosed artery,” Korea Australia Rheology Journal, vol. 21, no. 1, pp. 39–46, 2009. View at: Google Scholar
 D. Tang, C. Yang, H. Walker, S. Kobayashi, and D. N. Ku, “Simulating cyclic artery compression using a 3D unsteady model with fluidstructure interactions,” Computers and Structures, vol. 80, no. 2021, pp. 1651–1665, 2002. View at: Publisher Site  Google Scholar
 D. Tang, C. Yang, S. Kobayashi, J. Zheng, and R. P. Vito, “Effect of stenosis asymmetry on blood flow and artery compression: a threedimensional fluidstructure interaction model,” Annals of Biomedical Engineering, vol. 31, no. 10, pp. 1182–1193, 2003. View at: Publisher Site  Google Scholar
 T. Belzacq, S. Avril, E. Leriche, and A. Delache, “A numerical parametric study of the mechanical action of pulsatile blood flow onto axisymmetric stenosed arteries,” Medical Engineering and Physics, vol. 34, no. 10, pp. 1483–1495, 2012. View at: Publisher Site  Google Scholar
 Z. Li and C. Kleinstreuer, “Blood flow and structure interactions in a stented abdominal aortic aneurysm model,” Medical Engineering and Physics, vol. 27, no. 5, pp. 369–382, 2005. View at: Publisher Site  Google Scholar
 A. Valencia and F. Baeza, “Numerical simulation of fluid–structure interaction in stenotic arteries considering two layer nonlinear anisotropic structural model,” International Communications in Heat and Mass Transfer, vol. 36, no. 2, pp. 137–142, 2009. View at: Publisher Site  Google Scholar
 A. Valencia and M. Villanueva, “Unsteady flow and mass transfer in models of stenotic arteries considering fluidstructure interaction,” International Communications in Heat and Mass Transfer, vol. 33, no. 8, pp. 966–975, 2006. View at: Publisher Site  Google Scholar
 M. Schlaikjer, S. T. Petersen, J. A. Jensen, and P. F. Stetson, “Tissue motion in blood velocity estimation and its simulation,” in Proceedings of the International Ultrasonics Symposium, pp. 1495–1500, October 1998. View at: Google Scholar
 A. Swillens, L. Løvstakken, J. Kips, H. Torp, and P. Segers, “Ultrasound simulation of complex flow velocity fields based on computational fluid dynamics,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 56, no. 3, pp. 546–556, 2009. View at: Publisher Site  Google Scholar
 A. Swillens, J. Degroote, J. Vierendeels, L. Lovstakken, and P. Segers, “A simulation environment for validating ultrasonic blood flow and vessel wall imaging based on fluidstructure interaction simulations: ultrasonic assessment of arterial distension and wall shear rate,” Medical Physics, vol. 37, no. 8, pp. 4318–4330, 2010. View at: Publisher Site  Google Scholar
 A. Swillens, T. De Schryver, L. Løvstakken, H. Torp, and P. Segers, “Assessment of numerical simulation strategies for ultrasonic color blood flow imaging, based on a computer and experimental model of the carotid artery,” Annals of Biomedical Engineering, vol. 37, no. 11, pp. 2188–2199, 2009. View at: Publisher Site  Google Scholar
 S. Balocco, O. Basset, J. Azencot, P. Tortoli, and C. Cachard, “3D dynamic model of healthy and pathologic arteries for ultrasound technique evaluation,” Medical Physics, vol. 35, no. 12, pp. 5440–5450, 2008. View at: Publisher Site  Google Scholar
 A. Swillens, P. Segers, and L. Lovstakken, “Twodimensional flow imaging in the carotid bifurcation using a combined speckle tracking and phaseshift estimator: a study based on ultrasound simulations and in vivo analysis,” Ultrasound in Medicine and Biology, vol. 36, no. 10, pp. 1722–1735, 2010. View at: Publisher Site  Google Scholar
 D. F. Young and F. Y. Tsai, “Flow characteristics in models of arterial stenoses—I. Steady flow,” Journal of Biomechanics, vol. 6, no. 4, pp. 395–410, 1973. View at: Publisher Site  Google Scholar
 F. J. H. Gijsen, F. N. Van de Vosse, and J. D. Janssen, “The influence of the nonNewtonian properties of blood on the flow in large arteries: steady flow in a carotid bifurcation model,” Journal of Biomechanics, vol. 32, no. 6, pp. 601–608, 1999. View at: Publisher Site  Google Scholar
 K. Perktold and G. Rappitsch, “Computer simulation of local blood flow and vessel mechanics in a compliant carotid artery bifurcation model,” Journal of Biomechanics, vol. 28, no. 7, pp. 845–856, 1995. View at: Publisher Site  Google Scholar
 D. N. Ku, “Blood flow in arteries,” Annual Review of Fluid Mechanics, vol. 29, pp. 399–434, 1997. View at: Publisher Site  Google Scholar
 D. A. McDonald, Blood Flow in Arteries, Williams & Wilkins, Baltimore, Md, USA, 1974.
 D.K. Yao, Y.B. Yan, and Z.R. Liu, “The effects of blood viscoelasticity on the pulse wave in arteries,” Applied Mathematics and Mechanics, vol. 21, no. 9, pp. 1058–1065, 2000. View at: Publisher Site  Google Scholar
 L. ZhaoRong, “Pulsatile flow in arteries with external restraint of viscoelastic tissue,” Applied Mathematics and Mechanics, vol. 5, no. 4, pp. 1485–1491, 1984. View at: Publisher Site  Google Scholar
 Z. Liu, Cardiovascular Hydrodynamics, Fudan University Presses, Shanghai, China, 1986.
 X. Shixiong, J. Lin, and Y. Chew, “Hemodynamics on pulsatile flow through a vessel with mild stenosis amplitude oscillating radically,” Chinese Journal of Applied Biomechanics, vol. 14, pp. 202–207, 1999. View at: Google Scholar
 L. Gao, Y. Zhang, K. Zhang, G. Cai, J. Zhang, and X. Shi, “A computer simulation model for Doppler ultrasound signals from pulsatile blood flow in stenosed vessels,” Computers in Biology and Medicine, vol. 42, no. 9, pp. 906–914, 2012. View at: Publisher Site  Google Scholar
 J. A. Jensen, “Field: a program for simulating ultrasound systems,” in Proceedings of the 10th NordicBaltic Conference on Biomedical Imaging Published in Medical & Biological Engineering & Computing, vol. 4, supplement 1, part 1, pp. 351–353, 1996. View at: Google Scholar
 J. A. Jensen and N. B. Svendsen, “Calculation of pressure fields from arbitrarily shaped, apodized, and excited ultrasound transducers,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 39, no. 2, pp. 262–267, 1992. View at: Publisher Site  Google Scholar
 G. E. Tupholme, “Generation of acoustic pulses by baffled plane pistons,” Mathematika, vol. 16, no. 2, pp. 209–224, 1969. View at: Publisher Site  Google Scholar
 P. R. Stepanishen, “The timedependent force and radiation impedance on a piston in a rigid infinite planar baffle,” The Journal of the Acoustical Society of America, vol. 49, no. 3, pp. 841–849, 1971. View at: Publisher Site  Google Scholar
 J. A. Jensen and P. Munk, “Computer phantoms for simulating ultrasound BMode and CFM Images,” in Acoustical Imaging, vol. 23, pp. 75–80, Springer, 1997. View at: Publisher Site  Google Scholar
 J. A. Jensen, “Linear description of ultrasound imaging systems,” Notes for the International Summer School on Advanced Ultrasound Imaging, Technical University of Denmark, 1999. View at: Google Scholar
 J. A. Jensen, Users' Guide for the Field II Program, vol. 2800, Technical University of Denmark, 2001.
 C. Kasai, K. Namekawa, A. Koyano, and R. Omoto, “Realtime twodimensional blood flow imaging using an autocorrelation technique,” IEEE Transactions on Sonics and Ultrasonics, vol. 32, pp. 458–464, 1985. View at: Google Scholar
 Y. Yu and S. T. Acton, “Speckle reducing anisotropic diffusion,” IEEE Transactions on Image Processing, vol. 11, no. 11, pp. 1260–1270, 2002. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Qinghui Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.