Research Article  Open Access
Numerical Analysis on Flow Behavior of Molten Iron and Slag in Main Trough of Blast Furnace during Tapping Process
Abstract
The threedimensional model was developed according to number 4 of the main trough of blast furnace at China Steel Co. (CSC BF4). The equations and volume of fluid (VOF) were used for describing the turbulent flow at the impinging zone of trough, indicating fluids of liquid iron, molten slag, and air in the governing equation, respectively, in this paper. The pressure field and velocity profile were then obtained by the finite volume method (FVM) and the pressure implicit with splitting of operators (PISO), respectively, followed by calculating the wall shear stress through Newton’s law of viscosity for validation. Then, the operation conditions and the main trough geometry were numerically examined for the separation efficiency of iron from slag stream. As shown in the results, the molten iron losses associated with the slag can be reduced by increasing the height difference between the slag and iron ports, reducing the tapping rate, and increasing the height of the opening under the skimmer.
1. Introduction
Liquid phase separation plays a central role in a variety of production technologies [1]. It is important to reduce metal losses associated with the slag in the main trough of the blast furnace [2, 3]. Insufficient separation gives rise to technological problems in the subsequent processing of hot metal and slag in addition to the losses themselves. The trough, or main iron runner, is made of refractory material and designed to be long enough with suitable skimmer arrangements to obtain good separation of the two liquid phases of iron and slag. Depending on operating conditions, the mixture of molten iron and slag [4], which are continually tapped by the internal pressure of the blast furnace from the blast furnace hearth through a taphole, is separated by the density difference of fluids [5] in the downstream of main trough. As sketched in Figure 1(a), the separation of iron from slag proceeds as the liquids flow down the main trough; consequently, nearly complete separation is expected in the skimmer region. After separation through the skimmer arrangements, the slag is directed to subsequent slag processing facilities, and hot metal is sent to a desulphurization process prior to the steel making process [6].
(a)
(b)
A few investigations were carried out to study the fluid flow characteristics to enhance the separation efficiency. In the laminar flow regime, where Reynolds number [7] is approximately less than 0.1, the settling velocity can be expressed by Stoke’s law [8] in liquid Al during grain refinement that can be applied for the estimation of settling rate of metal droplets in the liquid slag [9]. For the case of static flow condition, the maximum velocity of droplets is a function of interfacial tension between two liquid phases, density of two liquid phases, and viscosity of continuous phase [10]. In the agitated mixing vessel, the mean diameter of droplets is a function of Weber number [11], viscosity ratio between two phases, and the geometrical factors [12]. An increase in the slag volume leads to a reduction in metal losses because of the greater pressure of the slag column on the pig iron [13]. However, there may exist a certain lower limit of the level when the slag penetration into molten iron [14] is taken into account. For separating two phases concurrent flow in a rectangular channel, the flow behavior was studied by Zhao et al. [15]. They reported that fractional height of upper phase stream (more viscous) is a function of the viscosity ratio and volume fraction of each liquid. Basically, this mechanism behavior will be applicable for the blast furnace trough system; accordingly, the interface level will be determined by the competition between gravitational force balance [16] in the trough and viscous force balance between two phases [17].
In the mold of a continuous caster [18], which is similar to blast furnace trough system, the slag entrainment mechanism has been studied. Feldbauer et al. [19] reported that there exists a critical velocity of molten steel necessary to cause slag entrainment and the velocity is a function of physical properties of fluids such as density, viscosity, and interfacial tension. Kim et al. [20] carried out a series of experiments, including changing the dimensions of the trough and various process variables such as total discharge rate from the furnace, volume fraction of each liquid, and viscosity of oil. Based on the results, the separation mechanism and the effect of each variable on the separation efficiency were identified. Furthermore, Luomala et al. [21] modeled fluid flows in the main trough of blast furnace to analyze the effects of the height difference between the slag and iron ports, the height of the bath level, and the height of the opening under the skimmer on the amount of entrained slag. As resulting from the calculation, the bigger the height of the iron layer above the opening under the skimmer, the smaller the amount of entrained slag.
As viewed from the above, the active modulation and control of the liquid phase separation are still challenge for hightemperature mixtures of iron and slag. In the present study, we aim to numerically analyze the separation efficiency of melt iron and slag in the main trough of blast furnace through the computer fluid dynamics (CFD) with the different operating conditions and main trough geometry. From the methodology point of view, an effective route to manipulate the separation of iron and slag in main trough of blast furnace depends on the geometry of trough, tapping rate, and the ratio of iron to slag in the fluids. At the present, according to our best knowledge, the computational fluid dynamics are firstly to systematically the parameters of main trough for studying the separation efficiency of molten iron from slag of blast furnace during tapping process.
2. Problem Formulation
2.1. Physical Model
The main trough of blast furnace number 4 at China Steel Co (CSC BF4) with 7200 tons/day of production capacity was modeled for the threedimensional computational fluid dynamics (CFD) to investigate the separation efficiency of molten iron and slag in this work. As expressed in the Cartesian (, , ) system, Figure 1 shows the side view and cross section of trough with the dimension in the unit of meter. The molten iron and slag were tapped into the main trough from the taphole tilted by the angle of 10 degrees and drained out at the slag port and iron dam, respectively. In the upstream portion of the trough, a considerable amount of turbulence of the bulk molten iron and slag arises from the kinetic energy of the falling stream. It is clear that the turbulent energy is variable and depends on the discharge rate of liquid, increasing with enlarging the taphole size. The slag separated from iron while traveling along the trough equipped with the skimmer due to the density difference and flowed out via slag port. Additionally, the thermal heat transfer was considered for forced convection induced by the difference temperature of tapping stream and the vertical wall of trough as well as natural convection due to the gravity force in the normal direction. In addition, according to literature [22], Table 1 lists the physical properties of liquid in this study.

2.2. Conservation Equation
For conservation equations, the following assumptions were made: (1) the chemical reaction [23], thermal radiation [24], and Marangoni effect [25] are ignorable; (2) air, iron, and slag are incompressible fluid [26]; (3) the interface of ironslag in the taphole was kept as a horizontal surface, and slag was floated on iron because of lighter density; (4) the determination of the size of liquid droplets in another immiscible liquid was not done; and (5) the main trough of blast furnace was unattacked. Accordingly, the equations of continuity and motion were expressed in the rectangular coordinates [27] asIn (1), the mean density iswith , where indicates the fraction of volume and is subscript noticing phase number, is the time, indicates the gravity force, and is the density of individual phase, which should be associated with the volume fraction equation [28],where is vector differential operator, also called nabla operator [29], and expresses the velocity vector including the velocities components of , , and , responding to , , and direction, respectively, and the effective viscosity is , where mean viscosity was defined bywhere is the viscosities of air, slag, and iron, respectively, denoted by phase number, , and the turbulent viscosity is , where the turbulent kinetic energy, , and energy dissipation rate, , were obtained by solving the equations [30],where the turbulence kinetic energy due to velocity gradient [31] is and due to the mean velocity fluctuation is , and the turbulence kinetic energy for thermal buoyancy [32] is , with thermal expansion coefficient, , subject to the momentum equation and temperature, .
In addition, conservation of energy is written aswhere the mean specific heat, , and mean heat conductivity, , were expressed in terms of volume of fluid aswhere and are the specific heat and thermal conductivity of individual fluid, respectively.
2.3. The Initial Condition
Under the conditions of large capacity of blast furnace equipped with multiple tapholes, interchanging tapping is done at time intervals for batch production of main trough. In the final discharging stage at first batch, the molten iron and slag cannot flow because the liquid levels are not over the iron dam and slag port, respectively, stagnating in the main trough. Figure 2 schemes the initial level of melt slag and iron in the main trough before next tapping, as . Since the pressure on both sides of skimmer, as denoted by the red start, is virtually the same, the layer thicknesses of melt slag and iron are m and m at the initial stage, which are calculated by (8) and (9), respectively [21].The nomenclatures used in (8) and (9) are, respectively, illustrated in the symbol list.
2.4. The Boundary Conditions
To simulate that the tapping stream is injected into the main trough from the taphole of blast furnace hearth and solve (1) to (6), the boundary conditions were specified as (1) the pressure which is contact surface under the atmosphere of 1 atm.; (2) the noslip conditions were imposed, and the wall temperature is set as 1773 K; (3) the rate ratio of molten iron to slag in the taphole was 3 : 1; (4) the velocity of iron at the inlet was expressed by , with, , , and being the unit vector at the , , and directions, respectively, whose turbulence intensity is 3.07%. (5) The velocity of slag at the inlet was expressed by , whose turbulence intensity is 6.63%, and the liquid temperature at the inlet is 1823 K. The turbulence intensities were calculated by the formula from the documentation [33]. The velocities of molten iron and slag at the inlet were calculated from the actual tapping rate, which was assumed to be 5 tons/min of molten iron on average.
3. Numerical Solution
The Eulerian model for a multiphase flow as implemented in the Fluent 6.2 [33] was used to solve conservation equations (1)~(6) that were separated using the finite volume technique in the computational domain ranging from the inlet at the taphole to the outlet at slap port and iron dam, respectively, as sketched in Figure 1, associated with the initial and boundary conditions. PressureImplicit SplitOperator (PISO) [34] was applied in pressurevelocity coupling in the segregated solver to adjust the velocity fields by correcting the pressure field. The second upwind scheme [35] was chosen to separate the convection term of governing equation. Figures 3 and 4 demonstrate the computational domain and grid of blast furnace main trough, respectively. In Figure 3, the virtual space was added to aid for simulating the actual tapping stream from the taphole to the main trough of blast furnace.
(a)
(b)
The resulting numerical solution was used to calculate the wall shear stress () and separation efficiency () of molten iron from slag at slag port by (10) and (11), respectively, at steady state.where is the normal vector.where the mass flow rate of melted iron at the slag port, , was obtained byand the mass flow rate of iron at the taphole, , was calculated byIn (10) and (11), and express the cross area at slag port and inlet, respectively. During computational period, the steady state was determined by the criteria, , which is defined aswhere and are the total mass flow rates at inlet and outlet, respectively, including the mass flow rates of molten iron and slag.
4. Results and Discussion
4.1. Validation
For the validation of computational grid, the wall shear stress profile in the direction at the position of m in the main trough affected by the mesh number was presented in Figure 5. As shown in the figure, the maximum wall shear stress at the distance from the taphole being about 3 m would be nearly a constant value when the mesh number is larger than 276705. Therefore, for the accuracy of the solution and saving computational time, the case C of mesh number was considered in this study.
Figure 6 shows the effect of natural convection on the velocity vector in the cross section at the buffer region of main trough. It was clearly seen that swirl flows of melt slag and iron occur in close to wall of main trough due to buoyancy dominated by the density of fluid as function of temperature, inducing the wall shear stress on the vertical wall in the downstream of main trough, as displayed in Figure 7.
In addition, it was found that the loss rate of molten iron from slag port varying with the flow rate of tapping stream obtained by numerical results is similar to the experimental result, whose water and oil were used to simulate the melted iron and slag, respectively, from literature [20], as shown in Figure 8. Therefore, the numerical model of this study could be applied for the sensitivity study on the separation efficiency of molten iron and slag by the different operation conditions and the geometry of blast furnace trough.
4.2. Effect of Discharge Rate
The discharge rates of the molten iron and slag would be significantly changed as the cast proceeds. This is mainly because the size of taphole was enlarged during tapping due to erosion of the refractory material. Figure 9 demonstrates the path line of molten iron and slag in terms of the density phase to be affected by the velocity of tapping stream near the slag port. As seen from the figure, the amount of iron trapped in the slag stream was increased with the velocity of tapping stream. When the metal droplets are under the influence of the bulk flow motion, their settling will be hindered and retention time would be shorted in the downstream of main trough. Therefore, the flow velocity of the liquids is an important factor to determine the separation efficiency of molten iron and slag, meaning that the enlarged tapping stream rate would increase the losses of molten iron at slag port.
(a)
(b)
(c)
4.3. Effect of Main Trough Geometry
The effect of main trough geometry on separation efficiency was studied by changing the high difference between iron dam and slag port and the depth of skimmer. The volume rate in main trough is considered for the geometry effect. In the main trough, the interface level of molten iron and slag was influenced by the flow velocity in each stream and the fluid dynamic balance between two phases of fluids, respectively. The molten iron lost from the slag port will be increased when the high difference between the slag port and iron dam is reduced, as shown in Figure 10, resulting from elevating the liquid level of molten iron in the downstream of main trough. Figure 11 displays that the separation efficiency of molten iron from the slag could be 100% when the height difference between iron dam and slag port is 0.2 m in this work. In contrast, when the height of iron dam is increased, the interface level will be elevated in the trough, resulting in the molten iron flowing into the slag port easily.
(a)
(b)
(c)
In addition, as sketched in Figure 1(a), the depth of skimmer is deeply leading to the interface level being risen and enhances the molted iron flowing into slag port. Figure 12 shows that the depth of skimmer is close to the bottom of the trough, where the turbulent kinetic energy is larger, strongly induced by the slag circulation flow in the upstream to mix the molten iron. This means that the more the skimmer is deeper to the trough bottom, the worst the separation efficiency of molten iron from slag at the slag port is. Figure 13 presents that the separation efficiency can be improved by 5.7%, while the depth of skimmer to the main trough bottom of blast furnace was changed from 0.1 to 0.3 m in the this study
(a)
(b)
(c)
5. Conclusion
In this study, the full threedimensional computational fluid dynamics composed from finite volume method and volume fraction equation in the present of natural thermal convection have been successfully developed to analyze the separation efficiency of molten iron and slag in the main trough of blast furnace. As resulting from the sensitivity analysis, particularly, it was evident that the lost rate of molten iron from the slag port is linearly increasing with increasing the velocity of tapping stream due to insufficient retention time of molten iron in the buffer region of main trough, and when the high difference between slag port and iron dam is less than 0.2 m, the separation efficiency of molten iron from slag could be more and more worse, resulting from elevating the liquid level of molten iron in the downstream of the trough. This suggests that the precise control of the high difference between iron dam and slag port is the priority of importance for maintaining the effectiveness of the slag and metal separation in the main trough of blast furnace during tapping process.
Symbols
:  The cross area at the inlet [m^{2}] 
:  The cross area at slag port [m^{2}] 
:  The mean specific heat [Jkg^{−1} K^{−1}] 
:  The specific heat of individual fluid [Jkg^{−1} K^{−1}] 
:  An empirical constant specified in the turbulence model (0.09) [—] 
:  The gravity acceleration [ms^{−2}] 
:  The turbulence kinetic energy for thermal buoyancy [kgm^{2}s^{−2}] 
:  The turbulence kinetic energy due to velocity gradient [kgm^{2}s^{−2}] 
:  The thicknesses of molten iron layers at the initial stage [m] 
:  The distance from the skimmer bottom to the top of iron dam [m] 
:  The thicknesses of molten slag at the initial stage [m] 
:  Unit vector in the direction [—] 
:  Unit vector in the direction [—] 
:  The turbulent kinetic energy [kgm^{2}s^{−2}] 
:  Unit vector in the direction [—] 
:  The distance from iron dam to the top of main trough [m] 
:  The depth of slag port [m] 
:  The depth of skimmer [m] 
:  Mass flow rate of molten iron at the taphole [kgs^{−1}] 
:  Total mass flow rate of molten iron at the taphole [kgs^{−1}] 
:  Mass flow rate of molten iron at the slag port [kgs^{−1}] 
:  Total mass flow rate of molten iron at the iron and slag ports [kgs^{−1}] 
:  Normal vector [m] 
:  Pressure [Pa] 
:  Index of fluid phase [—] 
:  Time [s] 
:  Temperature [K] 
:  Velocity in the direction [ms^{−1}] 
:  Velocity vector [ms^{−1}] 
:  Velocity of molten iron at the inlet [ms^{−1}] 
:  Velocity of molten slag at the inlet [ms^{−1}] 
:  Mean velocity fluctuation [ms^{−1}] 
:  Velocity in the direction [ms^{−1}] 
:  Velocity in the direction [ms^{−1}] 
, , :  Coordinates [m] 
:  The criteria defined by equation (14) [—] 
:  Vector differential operator [m^{−1}]. 
:  The fraction of volume [—] 
:  Thermal expansion coefficient [K^{−1}] 
:  Heat conductivity of individual fluid [Wm^{−1} K^{−1}] 
:  Mean heat conductivity [Wm^{−1} K^{−1}] 
:  Separation efficiency [—] 
:  The density of molten iron [kgm^{−3}] 
:  The mean density [kgm^{−3}] 
:  Density of individual fluid [kgm^{−3}] 
:  The effective viscosity [kgm^{−1}s^{−1}] 
:  The mean viscosity [kgm^{−1}s^{−1}] 
:  Viscosity of individual fluid [kgm^{−1}s^{−1}] 
:  The turbulent viscosity [kgm^{−1}s^{−1}] 
:  The wall shear stress [Pa]. 
:  Phase number indicating air, slag, and iron, respectively [—]. 
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
The authors would like to express sincere appreciation to Ministry of Science Technology (MOST: 1032221E005089) and New Materials Research & Development Department, China Steel Corporation for funding of this study.
References
 T. Letcher, J. Scot, and D. A. Patterson, Chemical Processes for a Sustainable Future, Royal Society of Chemistry, 2014.
 M. J. Luomala, T. T. Paanen, M. J. Köäykkä, T. M. J. Fabritus, H. Navla, and J. Härrhhi, “Modelling of fluid in the blast furnace trough,” Steel Research, vol. 4, pp. 130–135, 2001. View at: Google Scholar
 A. Kumar, S. Ali Khan, S. Biswas, and A. Pal, “Strategic steps towards longer and reliable blast furnace trough campaign—Tata Steel experience,” Ironmaking and Steelmaking, vol. 37, no. 1, pp. 15–20, 2010. View at: Publisher Site  Google Scholar
 A. E. A. Nogueira, M. B. Mourão, C. Takano, and D. M. dos Santos, “Effect of slag composition on iron nuggets formation from carbon composite pellets,” Materials Research, vol. 13, no. 2, pp. 191–195, 2010. View at: Publisher Site  Google Scholar
 N. I. Fedorova, S. A. Semenova, and Z. R. Ismagilov, “Material composition of fractions with different densities separated from slurry coal,” Solid Fuel Chemistry, vol. 47, no. 4, pp. 243–248, 2013. View at: Publisher Site  Google Scholar
 J. Lasota, “Refractory linings of pig iron transfer ladles,” Ceramic Materials, vol. 63, no. 3, pp. 692–695, 2011. View at: Google Scholar
 A. S. Semykina, V. I. Shatokha, S. Seetharaman, and O. Y. Gorobets, “Physical simulation of magnetite particle motion in the molten steelmaking slag,” Metallurgical and Mining Industry, vol. 3, no. 3, pp. 82–86, 2011. View at: Google Scholar
 G. S. V. Kumar, B. S. Murty, and M. Chakraborty, “Settling behaviour of TiAl_{3}, TiB_{2}, TiC and AlB_{2} particles in liquid Al during grain refinement,” International Journal of Cast Metals Research, vol. 23, no. 4, pp. 193–204, 2010. View at: Publisher Site  Google Scholar
 R. H. Eric, “Slag properties and design issues pertinent to matte smelting electric furnaces,” The Journal of the South African Institute of Mining and Metallurgy, vol. 104, no. 9, pp. 499–510, 2004. View at: Google Scholar
 P. Scheller and R. Hagemann, “Model investigations on slag entrainment in continuous casting,” Archives of Metallurgy and Materials, vol. 57, no. 1, pp. 283–289, 2012. View at: Publisher Site  Google Scholar
 S. SaintJean, A. Testa, L. V. Madden, and L. Huber, “Relationship between pathogen splash dispersal gradient and Weber number of impacting drops,” Agricultural and Forest Meteorology, vol. 141, no. 2–4, pp. 257–262, 2006. View at: Publisher Site  Google Scholar
 R. K. Thakur, C. Vial, G. Djelveh, and M. Labbafi, “Mixing of complex fluids with flatbladed impellers: effect of impeller geometry and highly shearthinning behavior,” Chemical Engineering and Processing: Process Intensification, vol. 43, no. 10, pp. 1211–1222, 2004. View at: Publisher Site  Google Scholar
 A. A. Antonov, M. F. Mar'yasov, and V. V. Khoroshikov, “Design changes to main troughs to reduce loss of pig iron with slag,” Metallurgist, vol. 28, no. 3, pp. 85–86, 1984. View at: Publisher Site  Google Scholar
 N. Fukami, R. Wakamatsu, N. Shinozaki, and K. Wasai, “Wettability between porous MgAl_{2}O_{4} substrates and molten iron,” Materials Transactions, vol. 50, no. 11, pp. 2552–2556, 2009. View at: Publisher Site  Google Scholar
 Y. Zhao, G. Chen, and Q. Yuan, “Liquidliquid twophase flow patterns in a rectangular MicroChannel,” AIChE Journal, vol. 52, no. 12, pp. 4052–4060, 2006. View at: Publisher Site  Google Scholar
 M. J. Golding, J. A. Neufeld, M. A. Hesse, and H. E. Huppert, “Twophase gravity currents in porous media,” Journal of Fluid Mechanics, vol. 678, pp. 248–270, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 O. Srámek, Y. Ricard, and D. Bercovici, “Simultaneous melting and compaction in deformable twophase media,” Geophysical Journal International, vol. 168, no. 3, pp. 964–972, 2007. View at: Publisher Site  Google Scholar
 Q. Yuan, B. G. Thomas, and S. P. Vanka, “Study of transient flow and particle transport in continuous steel caster molds: part I. Fluid flow,” Metallurgical and Materials Transactions B: Process Metallurgy and Materials Processing Science, vol. 35, no. 4, pp. 685–702, 2004. View at: Publisher Site  Google Scholar
 S. Feldbauer, J. Harmanand, and A. W. Cramb, “Slag entrainment in continuous casting,” CISR Reports, Carnegie Mellon University, 1994. View at: Google Scholar
 H. Kim, B. Ozturk, and R. J. Fruehan, “Slagmetal separation in the blast furnace trough,” ISIJ International, vol. 38, no. 5, pp. 430–439, 1998. View at: Publisher Site  Google Scholar
 M. J. Luomala, T. T. Paananen, M. J. Köykkä, T. M. Juhani Fabritius, H. Nevala, and J. J. Härkki, “Modelling of fluid flows in the blast furnace trough,” Steel Research, vol. 72, no. 4, pp. 130–135, 2001. View at: Publisher Site  Google Scholar
 I. Jimbom and A. W. Cramb, “The density of liquid ironcarbon alloys,” Ironmaking & Steelmaking, vol. 15, pp. 55–63, 1993. View at: Google Scholar
 H. Efendy, M. Safarudin, and H. Sihombing, “Molten metalslagrefractory reactions during converting process,” International Journal of Engineering & Technology, vol. 10, pp. 40–44, 2010. View at: Google Scholar
 A. Kumar, S. N. Bansal, and R. Chandraker, “Computational modeling of blast furnace cooling stave based on heat transfer analysis,” Materials Physics & Mechanics, vol. 15, pp. 46–56, 2012. View at: Google Scholar
 S. I. Karakashev and R. Tsekov, “Electromarangoni effect in thin liquid films,” Langmuir, vol. 27, no. 6, pp. 2265–2270, 2011. View at: Publisher Site  Google Scholar
 B. G. Thomas, “Modeling of continuouscasting defects related to mold fluid flow,” in Proceedings of the Congress on Science & Technology of Steel making, pp. 847–851, 2005. View at: Google Scholar
 R. B. Bird, W. E. Stewart, and E. N. Lightfoot, Transport Phenomena, John Wiley & Sons, New York, NY, USA, 2007.
 M. Owkes and O. Desjardins, “A computational framework for conservative, threedimensional, unsplit, geometric transport with application to the volumeoffluid (VOF) method,” Journal of Computational Physics, vol. 270, pp. 587–612, 2014. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. J. Pahl and R. Damrath, Mathematical Foundations of Computational Engineering, Springer, Berlin, Germany, 2001. View at: Publisher Site  MathSciNet
 A. K. A. Shati, S. G. Blakey, and S. B. M. Beck, “A dimensionless solution to radiation and turbulent natural convection in square and rectangular enclosures,” Journal of Engineering Science and Technology, vol. 7, no. 2, pp. 257–279, 2012. View at: Google Scholar
 R. Absi, “Analytical solutions for the modeled k equation,” Journal of Applied Mechanics, vol. 75, no. 4, 4 pages, 2008. View at: Publisher Site  Google Scholar
 R. Kuma and A. Dewan, “Assessment of buoyancycorrected turbulence models for thermal plumes,” Engineering Applications of Computational Fluid Mechanics, vol. 7, no. 2, pp. 239–249, 2013. View at: Publisher Site  Google Scholar
 S. Edwards, Fluent Aphasia, Cambridge Studies in Linguistics, Cambridge University Press, Cambridge, UK, 2005. View at: Publisher Site
 M. A. Sohail and R. Ullah, “CFD of oscillating airfoil pitch cycle by using PISO algorithm,” Engineering and Technology, vol. 60, pp. 1929–1933, 2011. View at: Google Scholar
 J. W. Banks and W. D. Henshaw, “Upwind schemes for the wave equation in secondorder form,” Journal of Computational Physics, vol. 231, no. 17, pp. 5854–5859, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Li Wang 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.