<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">ENG</journal-id><journal-title-group><journal-title>Engineering</journal-title></journal-title-group><issn pub-type="epub">1947-3931</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/eng.2021.138031</article-id><article-id pub-id-type="publisher-id">ENG-111370</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Design Optimization Analysis of Venturi Tube for Medium Conveying in Strengthen Grinding Process
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jinrui</surname><given-names>Xiao</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Zhongwei</surname><given-names>Liang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Xiaochu</surname><given-names>Liu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Zhuan</surname><given-names>Zhao</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Xincheng</surname><given-names>Xie</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib></contrib-group><aff id="aff3"><addr-line>School of Electromechanical Engineering, Guangdong University of Technology, Guangzhou, China</addr-line></aff><aff id="aff2"><addr-line>Guangzhou Key Laboratory of High-Performance Metal Grinding Processing, Guangzhou, China</addr-line></aff><aff id="aff1"><addr-line>School of Mechanical and Electrical Engineering, Guangzhou University, Guangzhou, China</addr-line></aff><pub-date pub-type="epub"><day>18</day><month>08</month><year>2021</year></pub-date><volume>13</volume><issue>08</issue><fpage>431</fpage><lpage>447</lpage><history><date date-type="received"><day>12,</day>	<month>July</month>	<year>2021</year></date><date date-type="rev-recd"><day>16,</day>	<month>August</month>	<year>2021</year>	</date><date date-type="accepted"><day>19,</day>	<month>August</month>	<year>2021</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Venturi tubes are widely used in manufacture and chemical industry and at
  t
  ract broad attention. To improve the transmission efficiency, the optimization of Venturi tube has been carried out, the transport substances involved in either liquid, solid, gas or two of them. The Venturi tubes that used for conveying of three-phase flow, which is consisted of solid, liquid and gas, are poorly investigated. In this paper, a Venturi tube employed in one surface treatment equipment was proposed. The velocity and pressure distribution in the flow field of Venturi tube with different key geometric parameters including diameter ratio, convergence angle and diffuser angle have been studied. One optimal solution has been selected. The negative pressure of suction port can reach a value as high as -
  1550 Pa to -
  1600 Pa. Such negative pressure offers the opportunities for transporting mixed abrasive.
 
</p></abstract><kwd-group><kwd>Venturi Tube</kwd><kwd> Strengthen Grinding</kwd><kwd> Medium Conveying</kwd><kwd> Parameter Optimization</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Over the past decades, Venturi tube has grown of importance in many applications such as cavitation bubbles [<xref ref-type="bibr" rid="scirp.111370-ref1">1</xref>], liquid flow metering [<xref ref-type="bibr" rid="scirp.111370-ref2">2</xref>], powder and particle transportation [<xref ref-type="bibr" rid="scirp.111370-ref3">3</xref>]. It’s worth pointing out that Venturi tube plays an extremely important role in the transportation of medium (i.e. powder, particles and liquid) in manufacturing industry, such as shot peening, water jet, sand blasting [<xref ref-type="bibr" rid="scirp.111370-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref8">8</xref>]. Strengthen grinding technology [<xref ref-type="bibr" rid="scirp.111370-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref10">10</xref>], which employs a three-phase jet flow that consists of abrasive powders, steel balls and grinding fluids to shoot workpiece, emerges as a new solution to improve the fatigue life of metals (See <xref ref-type="fig" rid="fig1">Figure 1</xref>). Such technology can reduce surface roughness and form a harden layer, which is believed to enhance the fatigue life of the metals since the micro-cutting effect and shock processing of abrasive powders and steel balls. Moreover, the grinding fluid helps clear dust and rust. Therefore, strengthen grinding technology opened a new window for processing metal parts such as bearing, gear, and mould.</p><p>Venturi tubes are among the most important components in strengthen grinding equipment, but they usually suffered from poor conveying and non-continuity feeding of mixed abrasive during working condition. The conveying efficiency is mainly impacted by the flow resistance. The flow resistance, which may weaken the conveying efficiency of the Venturi tube, is highly impacted by geometry parameters. It mainly reflects in velocity and pressure distribution of the flow field. To address this issue, many studies have been carried out [<xref ref-type="bibr" rid="scirp.111370-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref12">12</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref13">13</xref>]. An example of Venturi tube has been optimized to evaluate its effect on the conveying properties and the pressure reduction by flowing gas-solid mixtures [<xref ref-type="bibr" rid="scirp.111370-ref14">14</xref>]. Li [<xref ref-type="bibr" rid="scirp.111370-ref15">15</xref>] demonstrates that the outlet tangle and the diameter ratio of Venturi tubes are the two key parameters in controlling the throat pressure and power consumption. However, these works are mainly focused on either solid-gas or fluid transportation, the suction pressure inside is relatively low. Venturi tube used for transporting of three-jet flow has a higher requirement on suction pressure and the investigation is still absent.</p><p>In this paper, a Venturi tube that used for transporting mixed abrasive has been proposed and optimized. Finite element method (FEM) was taken into consideration and the pressure and velocity distribution in the flow field have been studied. The optimal parameters including the ratio between throat diameter and inlet diameter (k<sub>1</sub>), the ratio between suction diameter and throat diameter (k<sub>2</sub>), and the convergent angle (α) as well as the diffuser angle (β) were selected. Through further verification, the optimal solution performs well and has great potential to be used in transporting mixed medium.</p></sec><sec id="s2"><title>2. Modeling and Principle of the Tube</title><p>The proposed simplified Venturi tube is shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>. It includes two parts: main tube and suction tube, the suction tube located in the throat of main tube. <xref ref-type="table" rid="table1">Table 1</xref> reported the definitions of the geometric parameters.</p><p>When medium flow through main tube, the gas inside is carried away and the static pressure drops as the flow velocity increases. A negative pressure effect is formed once the flow velocity reaches a critical value, which help to absorb material from the suction port. The velocity increases rapidly in the process from inlet to throat since the diameter convergence. The maximum velocity appears in the throat tube and thus causes a maximum negative pressure. Therefore, the solid particles such as steel balls and abrasive powders can be sucked into the main tube by the negative pressure effect. Finally, the transportation of the solid material or mixed abrasive can be realized.</p><p>The pipeline pressure versus flow velocity is described by the following Bernoulli equation</p><p>P = C − 1 2 ρ v U 2 − ρ g h (1)</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Definition of geometric parameters</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Parameter</th><th align="center" valign="middle" >Definition</th></tr></thead><tr><td align="center" valign="middle" >D<sub>1</sub></td><td align="center" valign="middle" >The entrance diameter, which is equal to 20 mm.</td></tr><tr><td align="center" valign="middle" >D<sub>2</sub></td><td align="center" valign="middle" >The throat diameter</td></tr><tr><td align="center" valign="middle" >D<sub>3</sub></td><td align="center" valign="middle" >The outlet diameter, which is equal to 20 mm.</td></tr><tr><td align="center" valign="middle" >d</td><td align="center" valign="middle" >Diameter of the suction port</td></tr><tr><td align="center" valign="middle" >d<sub>i</sub></td><td align="center" valign="middle" >The length of the intersection of the suction pipe and the throat on the symmetry plane.</td></tr><tr><td align="center" valign="middle" >α</td><td align="center" valign="middle" >Convergent angle</td></tr><tr><td align="center" valign="middle" >β</td><td align="center" valign="middle" >Diffuser angle</td></tr><tr><td align="center" valign="middle" >AB</td><td align="center" valign="middle" >The intersecting line of suction pipe and throat pipe on the symmetry plane.</td></tr><tr><td align="center" valign="middle" >CD</td><td align="center" valign="middle" >A line that is perpendicular to the midpoint of line AB.</td></tr></tbody></table></table-wrap><p>where P is the pressure, U is the flow velocity, ρ is the fluid density, g is the gravitational acceleration, h is the height of the flow and C is a constant.</p><p>The suction force of the Venturi tube can be expressed as [<xref ref-type="bibr" rid="scirp.111370-ref16">16</xref>]</p><p>F d = 1 8 π C d ρ U 2 D 2 (2)</p><p>where C<sub>d</sub> is the drag coefficient of the flow around the medium, ρ is the density of the ambient flow, D is the cross-section diameter.</p><p>The drag coefficient of the flow around the medium can be expressed as [<xref ref-type="bibr" rid="scirp.111370-ref17">17</xref>]</p><p>C d = 24 R e ( 1 + b 1 R e b 2 ) + b 3 R e b 4 + R e (3)</p><p>where R<sub>e</sub> is the Reynolds number, b<sub>x</sub> is the coefficient that depends on the shapes of solid and the location of parts.</p><p>The mixed abrasive used for strengthen grinding is a mixture of steel balls and abrasive powder, and its shape is irregular. Here, a variate τ is introduced to describe the surface area of the mixed abrasive. And it’s defined as</p><p>τ = S 1 S 2 (4)</p><p>where S<sub>1</sub> is the surface area of the equivalent sphere, and S<sub>2</sub> is the surface area of the mixed abrasive. Then, the coefficient of b<sub>x</sub> can be expressed as [<xref ref-type="bibr" rid="scirp.111370-ref17">17</xref>]</p><p>b x = { exp ( 2.33 − 6.46 τ + 2.54 τ 2 )                                                   ( x = 1 ) 0.09 + 0.56 τ                                                                                         ( x = 2 ) exp ( 4.91 − 13.84 τ + 18.42 τ 2 − 10.26 τ 3 )                 ( x = 3 ) exp ( 1.47 + 12.26 τ − 20.73 τ 2 + 15.89 τ 3 )                 ( x = 4 ) (5)</p><p>The airflow density is related to the pressure and the velocity according to the Bernoulli equation, shown in Equation (6).</p><p>ρ = 2 C − 2 P U 2 (6)</p><p>Substituting Equation (3) and Equation (6) into Equation (2), the suction force can be expressed as</p><p>F d = [ 3 R e ( 1 + b 1 R e b 2 ) + b 3 R e 8 b 4 + 8 R e ] ( 2 C − 2 P ) D 2 (7)</p><p>with joint Equation (5) and Equation (7), the suction force of the Venturi tube can be got.</p><p>It can be seen that the suction force of Venturi tube is highly depend on the flow velocity, which was impacted by the throat diameter (D<sub>2</sub>), the suction port diameter (d), the convergent angle (α), the diffuser angle (β) and the inlet velocity. Therefore, the geometry parameters were taken into consideration in this study. To simplify the analysis, k<sub>1</sub> and k<sub>2</sub> were introduced, which were defined as k<sub>1</sub> = D<sub>2</sub>/D<sub>1</sub> and k<sub>2</sub> = d/D<sub>2</sub>, respectively.</p></sec><sec id="s3"><title>3. Methods and Simulation</title><p>This work was carried out in the fluent module of ANSYS Workbench software using FEM method. It was regarded as an effective tool for the design and optimization of Venturi tubes [<xref ref-type="bibr" rid="scirp.111370-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref19">19</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref20">20</xref>]. As mentioned before, a three-jet flow needs a higher suction pressure that caused by negative pressure effect, while negative pressure effect is highly depends on flow velocity. Hence, pressure and velocity were simulated. To improve simulation efficiency, a single medium of gas was selected and a 3-D model was considered.</p><p>The geometric parameters were selected from <xref ref-type="fig" rid="fig2">Figure 2</xref>, the boundary condition is shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>. The initial pressure in both the suction port and outlet port was set as zero. The flow velocity in the inlet was set as 20 m/s. The flow field is in a stable state and work at room temperature. The adaptive method was adopted for the meshing of model (See <xref ref-type="fig" rid="fig3">Figure 3</xref>).</p><p>The governed equations in the simulation are shown in the following:</p><p>∂ ρ ∂ t + ∇ ⋅ ( ρ U ) = 0 (8)</p><p>v = ∑ ​ 1 2 α 1 ρ 1 U 1 ρ (9)</p><p>ρ = ∑ ​ 1 2     α 1 ρ 1 (10)</p><p>where α 1 , ρ 1 and U 1 are the fraction, density and velocity, respectively.</p><p>The velocity dependence of pressure can be obtained by discrete momentum equation, which is shown in Equation (11).</p><p>ρ = ∑ ​ 1 2     α 1 ρ 1 (10)</p><p>The velocity dependence of pressure can be obtained by discrete momentum equation, which is shown in Equation (11).</p><p>a P U P r → + ∑ N a N U N r → = S → − ∇ P 0 (11)</p><p>where a P and a N are constants, S is a microunit on the surface of the spatial control body, U P r → is the predicted center velocity of the control body, U N r → is the predicted velocity around the control body and P 0 is the initial pressure.</p><p>The Semi-implicit method for pressure linked equations consistent (SIMPLEC) algorithm, which is based on continuity equation of Equation (12) and Equation (13) [<xref ref-type="bibr" rid="scirp.111370-ref21">21</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref22">22</xref>], was employed for the calculation. And the solving process of velocity and pressure in the Venturi tube was illustrated in <xref ref-type="fig" rid="fig4">Figure 4</xref>.</p><p>∇ ⋅ ( 1 a P + ∑ N a N ∇ P n ) = ∇ ⋅ [ H b y → A r + ( 1 a P + ∑ N a N − 1 a P ) ∇ P r ] (12)</p><p>U p n → = H b y → A r + ( 1 a P + ∑ N a N − 1 a P ) ∇ P r − 1 a P + ∑ N a N ∇ P n (13)</p><p>The variation parameters include k<sub>1</sub>, k<sub>2</sub>, α, β. A single variable method was used. k<sub>1</sub> and k<sub>2</sub> ranging from 0.4 to 0.9 with a step of 0.1, α ranging from 5˚ and 30˚ with a step of 5˚, β ranging from 4˚ to 14˚ with a step of 2˚. Details see in <xref ref-type="table" rid="table2">Table 2</xref>. After finish the first part of our simulation, a verification mode is used to assess its performance.</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> The parameters and simulation arrangement</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >No.</th><th align="center" valign="middle" >k<sub>1</sub></th><th align="center" valign="middle" >k<sub>2</sub></th><th align="center" valign="middle" >α (˚)</th><th align="center" valign="middle" >β (˚)</th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.6</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.7</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.8</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.9</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >8</td><td align="center" valign="middle" >0.6</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >9</td><td align="center" valign="middle" >0.7</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle" >0.8</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >11</td><td align="center" valign="middle" >0.9</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >13</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >10</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >14</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >15</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >16</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >25</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >17</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >30</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >18</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >6</td></tr><tr><td align="center" valign="middle" >19</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >8</td></tr><tr><td align="center" valign="middle" >20</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >21</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >12</td></tr><tr><td align="center" valign="middle" >22</td><td align="center" valign="middle" >0.4</td><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >14</td></tr></tbody></table></table-wrap></sec><sec id="s4"><title>4. Results and Discussion</title><p><xref ref-type="fig" rid="fig5">Figure 5</xref> reported the pressure distribution along the flow direction. It can be seen that static pressure initially falls as medium flow in the tubes, when arrives at the end of throat, the static pressure reaches its valley −1090 Pa. It is interesting to observe that the minimum static pressure achieves after medium pass through the throat. This could be explained by the absorption of extra air in the suction port, which caused oscillation. <xref ref-type="fig" rid="fig5">Figure 5</xref>(b) reported the static pressure distribution in the suction port projected to the axis of main tube. The pressure is overall in a stable state, ranging from −840 to −949 Pa. The pressure in the intersection part of suction pipe and throat pipe was significantly lower than the other parts of suction port. Hence, absorbing materials is become possible.</p><p>As mentioned above, the airflow velocity was closely related to the ratio of throat diameter to inlet diameter (k<sub>1</sub>), which controlled the cross-sectional area, under steady flow conditions. And it was confirmed again in <xref ref-type="fig" rid="fig6">Figure 6</xref>(a), where k<sub>1</sub> increased from 0.4 to 0.9 and the maximum airflow velocity of Venturi tube</p><p>decreased from 145 m/s to 27.2 m/s. But other variables have little effect on maximum velocity compare with k<sub>1</sub>. Moreover, the maximum velocity increases as k<sub>2</sub> increases (See <xref ref-type="fig" rid="fig6">Figure 6</xref>(b)), which is opposite to the trend in <xref ref-type="fig" rid="fig6">Figure 6</xref>(a). When k<sub>2</sub> increases, the diameter of suction pipe will also increase. Then the fluid drag force that courses by the suction pipe will decrease and leading to the increasing of maximum velocity in Venturi tube. An interesting fact that the velocity distribution of diffuser section is symmetrical when the convergent angle is 20˚ (See <xref ref-type="fig" rid="fig6">Figure 6</xref>(c)), whereas the high velocity of the others spread over the position that is opposite to the side of suction port when parameters change. No obvious velocity change can be found with the change of diffuser angle from 4˚ to 14˚ (See <xref ref-type="fig" rid="fig6">Figure 6</xref>(d)). Hence, a conclusion can be drawn that the value of diffuser angle has the least effect on the airflow velocity of Venturi tube among the four factors in this study.</p><sec id="s4_1"><title>4.1. Diameter Ratio Comparison</title><p>Axial airflow velocity in the intersection of suction pipe and throat pipe versus k<sub>1</sub> was reported in <xref ref-type="fig" rid="fig7">Figure 7</xref>. k<sub>1</sub> ranges from 0.4 to 0.9 with a step of 0.1. It can be seen that the airflow velocity is increasing as much as k<sub>1</sub> decreasing. This can be explained that when the inlet velocity and flow rate kept constant, the increase of k<sub>1</sub> will lead to an increase in the throat diameter, which in turn causes a decrease on airflow velocity. It is noteworthy that the airflow velocity is lower than 30 m/s and decreases with increasing distance from A when k<sub>1</sub> = 0.7, which is an idealess outcome. When k<sub>1</sub> ranges from 0.7 to 0.9, the airflow velocity is at a low value without any changing. But when k<sub>1</sub> = 0.4, the airflow velocity is significantly improved, which can reach a highest value of 84.5 m/s. Hence, k<sub>1</sub> = 0.4 is the best choice for the Venturi tube that used for mixed abrasive transportation in current study.</p><p>As indicated by <xref ref-type="fig" rid="fig8">Figure 8</xref>, axial airflow velocity in the intersection of suction pipe and throat pipe as function of d<sub>i</sub> has been studied. The flow velocity is increasing as much as the distance from point A increasing. Moreover, we can observe that airflow velocity is almost a constant and always at a low value when the distance from A more than 2.5 mm and k<sub>2</sub> of 0.4 and 0.9 have been selected. However, when the other distance from A less than 2.5 mm, the airflow velocity increases very quickly. A k<sub>2</sub> of 0.5 or 0.6 has a significant impact on the airflow velocity, which results in a quickly increasing of airflow velocity. Moreover, a highest airflow velocity of 76.3 m/s appears when a k<sub>2</sub> equal to 0.5. Based on the above analysis, k<sub>2</sub> = 0.5 is the best choice, which can create largest negative pressure and improve the efficiency of the proposed device.</p><p>Except for airflow velocity, the pressure distribution is also an important index to evaluate the performance of such device. The minimum pressure on CD line varies with k<sub>1</sub> and k<sub>2</sub>, which ranges from 0.4 to 0.9, is illustrated in <xref ref-type="fig" rid="fig9">Figure 9</xref>. When x = 1, then k<sub>1</sub> is 0.4 to 0.9 with a step of 0.1 and k<sub>2</sub> is 0.4. While x = 2, then k<sub>1</sub> is 0.4, and k<sub>2</sub> is 0.4 to 0.9 with a step of 0.1. The minimum pressure is gradually</p><p>increasing as k<sub>1</sub> increasing. It can be explained that the increase in k<sub>1</sub> will increase the diameter of the throat when k<sub>2</sub> remains unchanged, which in turn causes the flow velocity to decrease and the pressure to rise [<xref ref-type="bibr" rid="scirp.111370-ref23">23</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref24">24</xref>]. But with the increasement of k<sub>2</sub>, the minimum pressure decreases slightly and then increases. The change trend of minimum pressure is in agreement with the airflow velocity distribution trend of <xref ref-type="fig" rid="fig6">Figure 6</xref>. At the k<sub>2</sub> of 0.5, the minimum pressure reaches its lowest value ~ −4414.09 Pa. From here, we can observe that k<sub>2</sub> equal to 0.7 is a compromise. Furthermore, 0.4 and 0.5 are found to be the optimal selection of k<sub>1</sub> and k<sub>2</sub>, respectively.</p></sec><sec id="s4_2"><title>4.2. Convergent and Diffuser Angle Comparison</title><p>Furthermore, convergent angle (α) and diffuser angle (β) also influence the performance of the proposed Venturi tube (See in <xref ref-type="fig" rid="fig6">Figure 6</xref>(a) and <xref ref-type="fig" rid="fig6">Figure 6</xref>(b)), which also need to take into consideration. Then the airflow velocity versus different α and β have also been investigated (See <xref ref-type="fig" rid="fig1">Figure 1</xref>0). The α is range from 5˚ to 30˚ with a step of 5˚ and the β is range from 4˚ to 14˚ with a step of 2˚. Compared with other parameters, α and β have a relatively small impact on the airflow velocity, which is similar to the results reported by Haifeng Lu et al. [<xref ref-type="bibr" rid="scirp.111370-ref25">25</xref>].</p><p>We can observe that the airflow velocities increase while the distances from point A increase (See <xref ref-type="fig" rid="fig1">Figure 1</xref>0), which is similar to <xref ref-type="fig" rid="fig7">Figure 7</xref> and <xref ref-type="fig" rid="fig8">Figure 8</xref>. Moreover, when α ranges from 5˚ to 20˚, the airflow velocity is increasing as much as β increasing, while the opposite trend is found with α ranges from 20˚ to 30˚. The higher velocity is more benefit for particle transportation in the venturi tube [<xref ref-type="bibr" rid="scirp.111370-ref26">26</xref>]. Hence, α = 20˚ has been selected since it possesses the highest airflow velocity of 90.5 m/s (See <xref ref-type="fig" rid="fig1">Figure 1</xref>0(a)). The airflow velocity is significantly higher than the others when β equal to 4˚ (See <xref ref-type="fig" rid="fig1">Figure 1</xref>0(b)). This is due to that the larger diffuser angle, the greater the energy loss, which causes by the boundary layer separation and eddy effects [<xref ref-type="bibr" rid="scirp.111370-ref27">27</xref>]. Even so, the larger diffuser angles also induce velocity sharply fluctuation that may reduce the kinetic energy. Since the higher velocity is more beneficial for mixed material conveying, the selection of diffuser angle should also take higher airflow velocity into account firstly. Finally, a convergent angle of 20˚ and a diffuser angle of 4˚ were selected.</p></sec><sec id="s4_3"><title>4.3. Verification of the Optimal Parameters</title><p>In addition, comparative analysis between the optimal parameters and the others were then conducted for verifying the validity of the selection. As mentioned previously, k<sub>1</sub> = 0.4, k<sub>2</sub> = 0.5, α = 20˚ and β = 4˚ are the optimal selections. To confirm the effectiveness of the optimization results, the selected optimal parameters were compared with the others in the previous analysis according to the airflow velocity and static pressure distribution (See in <xref ref-type="fig" rid="fig1">Figure 1</xref>1 and <xref ref-type="fig" rid="fig1">Figure 1</xref>2). We</p><p>can find that the optimal one obtained the maximum airflow velocity of 154 m/s, which is only 1 m/s lower that can be ignored comparing with the highest velocity. Moreover, the airflow velocity distribution curve on AB line of the optimal one is higher than that of the others (See in <xref ref-type="fig" rid="fig1">Figure 1</xref>2(a)). As a result, reasonable negative pressure appears on the CD line of the optimal one. In addition, the optimal one possesses the lowest negative pressure, which is in the range of −1550 Pa and −1600 Pa, in the suction port (See in <xref ref-type="fig" rid="fig1">Figure 1</xref>3). It is consistent with the results shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>2. Moreover, to further verify the reliability of the optimal parameters, the acceleration coefficient (ε) of the proposed Venturi tube is introduced, which is defined as shown in Equation (14).</p><p>ε = U s U o &#215; D 2 D 1 (14)</p><p>where U<sup>s</sup> is the maximum velocity on AB.</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Comparisons of acceleration coefficient</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Proposed structure</th><th align="center" valign="middle" >ε</th></tr></thead><tr><td align="center" valign="middle" >The optimal one in this study</td><td align="center" valign="middle" >1.856</td></tr><tr><td align="center" valign="middle" >Reference [<xref ref-type="bibr" rid="scirp.111370-ref16">16</xref>]</td><td align="center" valign="middle" >0.542</td></tr><tr><td align="center" valign="middle" >Reference [<xref ref-type="bibr" rid="scirp.111370-ref28">28</xref>]</td><td align="center" valign="middle" >1.067</td></tr><tr><td align="center" valign="middle" >Reference [<xref ref-type="bibr" rid="scirp.111370-ref29">29</xref>]</td><td align="center" valign="middle" >1.057</td></tr></tbody></table></table-wrap><p>It can be seen from Equation (14) that the acceleration coefficient (ε) represents the ratio of the product of the maximum velocity of the suction port and the diameter of the throat to the product of the inlet velocity and diameter. When the acceleration coefficient is larger, it means that the Venturi tube has a more significant speed-up effect, which is more conducive to the suction and transportation of mixed medium. Then the comparation of the optimal one and the others that proposed by references [<xref ref-type="bibr" rid="scirp.111370-ref16">16</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref28">28</xref>] [<xref ref-type="bibr" rid="scirp.111370-ref29">29</xref>] were carried out, as shown in <xref ref-type="table" rid="table3">Table 3</xref>.</p><p>Evidence can be found that the acceleration coefficient of the optimal one is 1.856, which is the greatest among the four. Furthermore, it’s 73.4% higher than the next largest, which proposed in reference [<xref ref-type="bibr" rid="scirp.111370-ref28">28</xref>]. Obviously, the optimal one meets the requirements of higher speed. Hence, the selected parameters are approved to be credible for the Venturi tube of strengthen grinding equipment.</p></sec></sec><sec id="s5"><title>5. Conclusions</title><p>The geometric optimization of Venturi jet tubes for medium conveying in strengthen grinding equipment was performed by FEM. The following conclusions can be drawn based on the analysis.</p><p>When the diameter ratios of k<sub>1</sub> and k<sub>2</sub> are equal to 0.4 and 0.5, a better acceleration effect and negative pressure effect achieve, which is more conducive to the transportation of medium during the strengthen grinding process. A convergence angle of 20˚ and a diffuser angle of 4˚ had selected for the optimal angle of the proposed Venturi tube in this study. The suction force mathematical model for medium conveying of different sizes particle by Venturi tube had been established, which could help the parameter design of Venturi tube.</p><p>A single-phase flow is used to simulate the flow field of the Venturi, and the rationality of the structure is evaluated by the flow velocity and pressure, which constitutes the limitations of this study. And the multiphase flow simulation and experiment need to be carried out in the future.</p></sec><sec id="s6"><title>Acknowledgements</title><p>This research was supported by the Science and Technology Planning Project of Guangdong China (2019B020404) and the Science and Technology Planning Project of Guangzhou China (202102080225).</p></sec><sec id="s7"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s8"><title>Cite this paper</title><p>Xiao, J.R., Liang, Z.W., Liu, X.C., Zhao, Z. and Xie, X.C. (2021) Design Optimization Analysis of Venturi Tube for Medium Conveying in Strengthen Grinding Process. Engineering, 13, 431-447. https://doi.org/10.4236/eng.2021.138031</p></sec></body><back><ref-list><title>References</title><ref id="scirp.111370-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Long, X.P., Zhang, J.Q., Wang, J., Xu, M.S., Lyu, Q. and Ji, B. (2017) Experimental Investigation of the Global Cavitation Dynamic Behavior in a Venturi Tube with Special Emphasis on the Cavity Length Variation. International Journal of Multiphase Flow, 89, 290-298. https://doi.org/10.1016/j.ijmultiphaseflow.2016.11.004</mixed-citation></ref><ref id="scirp.111370-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Brinkhorst, S., Lavante, E.V. and Wendt, G. (2015) Numerical Investigation of Cavitating Herschel Venturi-Tubes Applied to Liquid Flow Metering. Flow Measurement and Instrumentation, 43, 23-33. https://doi.org/10.1016/j.flowmeasinst.2015.03.004</mixed-citation></ref><ref id="scirp.111370-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Bertoldi, D., Dallalb, C.C. and Barbosa, J.R. (2015) Experimental Investigation of Two-Phase Flashing Flows of a Binary Mixture of Infinite Relative Volatility in a Venturi Tube. Experimental Thermal and Fluid Science, 64, 152-163. https://doi.org/10.1016/j.expthermflusci.2015.02.011</mixed-citation></ref><ref id="scirp.111370-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Kumar, R.K., Kumaran, P.S., Seetharamu, S., Kumar, S.A., Pramod, T. and Naveen, G.J. (2019) Investigation of Shot Peening Effect on Titanium Alloy Affecting Surface Residual Stress and Roughness for Aerospace Applications. Procedia Structural Integrity, 14, 134-141. https://doi.org/10.1016/j.prostr.2019.05.018</mixed-citation></ref><ref id="scirp.111370-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Zhang, Z., Lin, M., Seng, D.H.L., Teo, S.L., Wei, F.X., Tan, H.R., Cheong, A.K.H., Lim, S.H., Wang, S.J. and Pan, J.S. (2020) Fatigue Life Enhancement in Alpha/Beta Ti-6al-4v after Shot Peening: An EBSD and Tem Crystallographic Orientation Mapping Study of Surface Layer. Materialia, 12, Article ID: 100813. https://doi.org/10.1016/j.mtla.2020.100813</mixed-citation></ref><ref id="scirp.111370-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Urazmetov, O., Cadet, M., Teutsch, R. and Antonyuk, S. (2021) Investigation of the Flow Phenomena in High-Pressure Water Jet Nozzles. Chemical Engineering Research and Design, 165, 320-332. https://doi.org/10.1016/j.cherd.2020.10.030</mixed-citation></ref><ref id="scirp.111370-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Kanesan, D., Mohyaldinn, M.E., Ismail, N.I., Chandran, D. and Liang, C.J. (2019) An Experimental Study on the Erosion of Stainless Steel Wire Mesh Sand Screen Using Sand Blasting Technique. Journal of Natural Gas Science and Engineering, 65, 267-274. https://doi.org/10.1016/j.jngse.2019.03.017</mixed-citation></ref><ref id="scirp.111370-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Liu, X.C., Huang, J.F., Liang, Z.W., Huang, W.F., Zhu, R., Gao, W.L. and Xiao, J.R. (2020) Preparation and Properties of the Composite Enhancement Layer of Bearing Ring. World Journal of Mechanics, 10, 139-153. https://doi.org/10.4236/wjm.2020.1010010</mixed-citation></ref><ref id="scirp.111370-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Liang, Z.W., Xie, B.H., Liao, S.P. and Zhou, J.H. (2015) Concentration Degree Prediction of AWJ Grinding Effectiveness Based on Turbulence Characteristics and the Improved ANFIS. The International Journal of Advanced Manufacturing Technology, 80, 887-905. https://doi.org/10.1007/s00170-015-7027-0</mixed-citation></ref><ref id="scirp.111370-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Liu, X.C., Liang, Z.W., Wen, G.L. and Yuan, X.F. (2019) Waterjet Machining and Research Developments: A Review. The International of Advanced Manufacturing Technology, 102, 1257-1335. https://doi.org/10.1007/s00170-018-3094-3</mixed-citation></ref><ref id="scirp.111370-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Wu, H.T., Xu, Y., Wang, J.H., Zhang, T. and Wang, H.X. (2019) Study on the Similarity of Wet Gas Pressure Drop in Long-Throat Venturi. Flow Measurement and Instrumentation, 68, Article ID: 101580. https://doi.org/10.1016/j.flowmeasinst.2019.101580</mixed-citation></ref><ref id="scirp.111370-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Kurimoto, R., Nakazawa, K., Minagawa, H. and Yasuda, T. (2017) Prediction Models of Void Fraction and Pressure Drop for Gas-Liquid Slug Flow in Microchannels. Experimental Thermal and Fluid Science, 88, 124-133. https://doi.org/10.1016/j.expthermflusci.2017.05.014</mixed-citation></ref><ref id="scirp.111370-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Wu, H.T., Xu, Y., Xiong, X.Q., Mamat, E., Wang, J.H. and Zhang, T. (2020) Prediction of Pressure Drop in Venturi Based on Drift-Flux Model and Boundary Layer Theory. Flow Measurement and Instrumentation, 71, Article ID: 101673. https://doi.org/10.1016/j.flowmeasinst.2019.101673</mixed-citation></ref><ref id="scirp.111370-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Liu, K., Guo, X.L., Peng, G. and Xin, H.F. (2017) Design Optimization of a Venturi Tube Geometry in Dense-Phase Pneumatic Conveying of Pulverized Coal for Entrained-Flow Gasification. Chemical Engineering Research &amp; Design: Transactions of the Institution of Chemical Engineers, 120, 208-217. https://doi.org/10.1016/j.cherd.2017.02.020</mixed-citation></ref><ref id="scirp.111370-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Li, M., Bussonni`ere, A., Bronson, M., Xu, Z. and Liu, Q.X. (2019) Study of Venturi Tube Geometry on the Hydrodynamic Cavitation for the Generation of Microbubbles. Minerals Engineering, 132, 268-274. https://doi.org/10.1016/j.mineng.2018.11.001</mixed-citation></ref><ref id="scirp.111370-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Ji, D., Zhang, M.L., Xu, T.G., Wang, K.Y., Li, P.C. and Ju, F. (2015) Experimental and Numerical Studies of the Jet Tube Based on Venturi Effect. Vacuum, 111, 25-31. https://doi.org/10.1016/j.vacuum.2014.09.010</mixed-citation></ref><ref id="scirp.111370-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Haider, A. and Levenspiel, O. (1989) Drag Coefficient and Terminal Velocity of Spherical and Nonspherical Particles. Powder Technology, 58, 63-70. https://doi.org/10.1016/0032-5910(89)80008-7</mixed-citation></ref><ref id="scirp.111370-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Soyama, H. (2021) Luminescence Intensity of Vortex Cavitation in a Venturi Tube Changing with Cavitation Number. Ultrasonics Sonochemistry, 71, Article ID: 105389. https://doi.org/10.1016/j.ultsonch.2020.105389</mixed-citation></ref><ref id="scirp.111370-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Song, Y.C., Shen, Y., Qian, Y.L., Yin, J.L. and Wang, D.Z. (2021) Experiment and Modeling of Liquid-Phase Flow in a Venturi Tube Using Stereoscopic PIV. Nuclear Engineering and Technology, 53, 79-92. https://doi.org/10.1016/j.net.2020.06.027</mixed-citation></ref><ref id="scirp.111370-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Wu, H.T., Xu, Y., Wang, J.H., Yang, Y.G., Li, T., Zhang, T., Li, J.L. and Hao, C.J. (2021) Gas-Liquid Two-Phase Flowrate Measurement in Pseudo-Slug Flow with Venturi. Flow Measurement and Instrumentation, 78, Article ID: 101887. https://doi.org/10.1016/j.flowmeasinst.2021.101887</mixed-citation></ref><ref id="scirp.111370-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Serra, N. and Semiao, V. (2021) ESIMPLE, a New Pressure-Velocity Coupling Algorithm for Built-Environment CFD Simulations. Building Environment, 204, Article ID: 108170. https://doi.org/10.1016/j.buildenv.2021.108170</mixed-citation></ref><ref id="scirp.111370-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Aguerre, H., Venier, C.M., Pairetti, C.I., Damián, S.M. and Nigro, N.M. (2020) A SIMPLE-Based Algorithm with Enhanced Velocity Corrections: The COMPLEX Method. Computers &amp; Fluids, 198, Article ID: 104396. https://doi.org/10.1016/j.compfluid.2019.104396</mixed-citation></ref><ref id="scirp.111370-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Mishra, C. and Peles, Y. (2006) An Experimental Investigation of Hydrodynamic Cavitation in Micro-Venturis. Physics of Fluids, 18, Article ID: 103603. https://doi.org/10.1063/1.2360996</mixed-citation></ref><ref id="scirp.111370-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Sharifi, A. and Mohebbi, A. (2014) A Combined CFD Modeling with Population Balance Equation to Predict Pressure Drop in Venturi Scrubbers. Research on Chemical Intermediates, 40, 1021-1042. https://doi.org/10.1007/s11164-013-1018-2</mixed-citation></ref><ref id="scirp.111370-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Lu, H.F., Guo, X.L., Li, P., Liu, K. and Gong, X. (2017) Design Optimization of a Venturi Tube Geometry in Dense-Phase Pneumatic Conveying of Pulverized Coal for Entrained-Flow Gasification. Chemical Engineering Research and Design, 120, 208-217. https://doi.org/10.1016/j.cherd.2017.02.020</mixed-citation></ref><ref id="scirp.111370-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Novo, P.G., Kyozuka, Y. and Villamayor, M.J.G. (2019) Evaluation of Turbulence-Related High-Frequency Tidal Current Velocity Fluctuation. Renewable Energy, 139, 313-325. https://doi.org/10.1016/j.renene.2019.02.035</mixed-citation></ref><ref id="scirp.111370-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Liu, K., Lu, H.F., Guo, X.L., Sun, X.L., Tao, S.L. and Gong, X. (2014) Experimental Study on Flow Characteristics and Pressure Drop of Gas-Coal Mixture through Venturi. Powder Technology, 268, 401-411. https://doi.org/10.1016/j.powtec.2014.07.019</mixed-citation></ref><ref id="scirp.111370-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Fang, L., Li, W., Li, Q. and Wang, Z.B. (2020) Numerical Investigation of the Cavity Shedding Mechanism in a Venturi Reactor. International Journal of Heat and Mass Transfer, 156, Article ID: 119835. https://doi.org/10.1016/j.ijheatmasstransfer.2020.119835</mixed-citation></ref><ref id="scirp.111370-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Ahmed, S., Mohsin, H., Qureshi, K., Shah, A., Siddique, W., Waheed, K., Irfan, N., Ahmad, M. and Farooq, A. (2018) Investigation of Dust Particle Removal Efficiency of Self-Priming Venturi Scrubber Using Computational Fluid Dynamics. Nuclear Engineering and Technology, 50, 665-672. https://doi.org/10.1016/j.net.2018.01.016</mixed-citation></ref></ref-list></back></article>