<?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">OJMSi</journal-id><journal-title-group><journal-title>Open Journal of Modelling and Simulation</journal-title></journal-title-group><issn pub-type="epub">2327-4018</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojmsi.2017.53013</article-id><article-id pub-id-type="publisher-id">OJMSi-77788</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Simulation Analysis of Flow Characteristics in the Vortex Diode with Double Tangential Tube
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Lijun</surname><given-names>Zhang</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>Yuelai</surname><given-names>Dai</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>Yulei</surname><given-names>Li</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>Qun</surname><given-names>Wang</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>Xiaoyu</surname><given-names>Jin</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>Lan</surname><given-names>Ma</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>Wentao</surname><given-names>Huang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Ningbo Sub-Academy, National Weapons Science Research Academy, Ningbo, China</addr-line></aff><pub-date pub-type="epub"><day>23</day><month>06</month><year>2017</year></pub-date><volume>05</volume><issue>03</issue><fpage>174</fpage><lpage>181</lpage><history><date date-type="received"><day>April</day>	<month>29,</month>	<year>2017</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>July</month>	<year>18,</year>	</date><date date-type="accepted"><day>July</day>	<month>21,</month>	<year>2017</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-NonCommercial International License (CC BY-NC).http://creativecommons.org/licenses/by-nc/4.0/</license-p></license></permissions><abstract><p>
 
 
  The paper analyzed flow characteristics in vortex diode with double tangential tube by numerical simulation. Compared to traditional vortex diode with one tangential tube, the structure of double tangential tube can reduce the forward resistance and improve overall performance. The symmetrical design adopted in vortex diode with double tangential tube; the internal flow field in chamber showed symmetrical distribution in the reverse flow, which can improve the flow stability and accelerate the convergence speed in simulation.
 
</p></abstract><kwd-group><kwd>Vortex Diode</kwd><kwd> Double Tangential Tube</kwd><kwd> Performance Parameters</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>As a new maintenance-free device, vortex diode pump device begins to be used in nuclear industry for spent fuel transportation, and is extended to the chemical and pharmaceutical industries. The basic working principle of vortex diode is that its working fluid through the pros and cons vortex diode has different drag coefficients, which is similar to the role of “diode” in the circuit. The structure of Vortex diode includes three parts: tangential tube, axial tube and chamber. In the case of forward flow, the working fluid flows into the chamber from the axial tube, and outflows from the tangential tube. Its internal flow field is relatively simple. In the case of reverse flow, the working fluid flows into the chamber from tangential tube and outflows from the axial tube. The internal field is rendered as complex three-dimensional turbulent flow field, in which a strong vortex is formed near the center of chamber. The most important parameter to measure its overall performance is impedance ratio which is defined as the ratio of reverse resistance coefficient to the forward resistance coefficient.</p><p>Many scholars have done a lot of work in studying the structure of vortex diode. Tapered and divergent structure was used in the axial tube. Taper angle of tangential tube was changed and the chamber in configuration was designed. All above is in order to improve its overall performance or one-way performance. In addition, some calculation models were considered in reverse flow. Malhotra A used standard k-ε model to calculate swirl problem, which had a degree of error [<xref ref-type="bibr" rid="scirp.77788-ref1">1</xref>] . Realizable k-ε model and RNG k-ε model were used by Speziale C G to calculate turbulence model for separated flows; the results were better than standard k-ε model, but not ideal [<xref ref-type="bibr" rid="scirp.77788-ref2">2</xref>] . Reynolds stress model was very close to the experimental results, verified by Launder B E [<xref ref-type="bibr" rid="scirp.77788-ref3">3</xref>] . Large eddy simulation is the best model in swirl flow, verified by Pierce C D, but quite time-consuming [<xref ref-type="bibr" rid="scirp.77788-ref4">4</xref>] . In overall consideration, Reynolds stress model is suitable to calculate swirl problem. However, less research is done on vortex diode with double tangential tube.</p><p>This paper analyzed flow characteristics in vortex diode with double tangential tube by numerical simulation. Vortex diode with double tangential tube shows better performance in the forward flow and overall performance, compared with traditional vortex diode. What’s more, the convergence performance is improved, which can reduce computing time.</p></sec><sec id="s2"><title>2. Working Mechanism and Characteristic Parameters</title><p>The traditional vortex diode structure is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>, whose working mechanism is as follows. When the fluid flows from the tangential tube into the axial tube, the working fluid jets into the chamber due to the contraction of the pipe, so that a vortex is formed by the action of the wall. This leads to a great flow resistance, which is called reverse flow. When the working fluid flows from the axial tube into the tangential tube, the flow, known as the forward flow, is relatively simple. In this condition, no swirl flow is generated so that the flow resistance is relatively small. The main performance parameters of vortex diode include the resistance coefficient ξ (Euler number Eu) and the impedance ratio ε:</p><disp-formula id="scirp.77788-formula100"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2860120x2.png"  xlink:type="simple"/></disp-formula><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> The 3D model of vortex diode</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2860120x3.png"/></fig><disp-formula id="scirp.77788-formula101"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2860120x4.png"  xlink:type="simple"/></disp-formula><p>where ∆P is the difference between input pressure and output pressure and ρ is the density of working fluid. u is the velocity of working fluid. Eu<sub>F</sub> is the forward flow resistance coefficient and Eu<sub>R</sub> is the reverse flow resistance coefficient.</p></sec><sec id="s3"><title>3. The Model</title><sec id="s3_1"><title>3.1. Geometry Model</title><p>The vortex diode structure of double tangential tube is shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>. The structure of the two tangential tubes set symmetrically are the same. Cavity diameter of chamber is 150 mm, and the tube diameter of import and export are both 25 mm. Pre-treatment of finite model is carried out by ICEM software. Hybrid meshes are generated using structured grids and unstructured grids, and tetrahedral meshes are applied to joints between chamber and tangential tube. The rest part is meshed by hexahedral meshes. Mesh quality is good, and the total number of grids is more than 80 million.</p></sec><sec id="s3_2"><title>3.2. Calculation Model</title><p>In Forward flow, for it is in isotropic turbulence state, so the calculation model uses the standard k-epsilon model. In reverse flow, the fluid inside the chamber presents vortex motion, which shows the performance of anisotropy and instability. For the swirl problem, standard k-ε model, realizable k-ε model, RNG k-ε model, Reynolds stress model and large eddy simulation model can all be considered. Although realizable k-ε model and RNG k-ε model are more reliable than standard k-ε model, they are still not ideal. Reynolds stress model directly models individual items, which is better than the two models above in simulation results. Large eddy simulation has the best performance, but it is time-consuming during computation. Given simulation accuracy and economy, Reynolds stress model is appropriate to be used to calculate in reverse flow.</p></sec><sec id="s3_3"><title>3.3. Boundary Conditions and Solution Control</title><p>The import condition of forward and reverse flow are set as speed of working</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> The 3D model of vortex diode with double tangential tube</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2860120x5.png"/></fig><p>fluid. Outlet pressure is set as the constant 0, while a given hydraulic diameter and turbulence intensity are taken as import turbulent conditions. In reverse</p><p>flow simulation, import Reynolds stress values meet<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2860120x6.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2860120x7.png" xlink:type="simple"/></inline-formula>, and export Reynolds stress values satisfy<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2860120x8.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2860120x9.png" xlink:type="simple"/></inline-formula> [<xref ref-type="bibr" rid="scirp.77788-ref5">5</xref>] . Due to a long pe-</p><p>riod of swirl at the outlet, the outlet pressure produces reflux phenomenon, so the outlet boundary condition is set as outlet-vent and reflux is eliminated by increasing the resistance loss coefficient at the outlet. Equations are solved in Fluent, and the reference pressure (Pressure-Based) algorithm is adopted. In addition, optimization polyhedron tetrahedral mesh grid technology is applied, and PRESTO is used to discrete items of pressure as well. The SIMPLE algorithm is applied to format, pressure and velocity coupling, and second order upwind scheme is applied to Momentum and turbulence kinetic energy dissipation rate too [<xref ref-type="bibr" rid="scirp.77788-ref6">6</xref>] . Fluid field are considered as steady flow, and temperature influence is ignored. Working media is water for simulation. Convergence criteria is that importers mass flow flux difference ratio is less than 10<sup>−</sup><sup>5</sup>.</p></sec></sec><sec id="s4"><title>4. Numerical Results</title><sec id="s4_1"><title>4.1. Contrast of Performance Parameters</title><p><xref ref-type="fig" rid="fig3">Figure 3</xref> and <xref ref-type="fig" rid="fig4">Figure 4</xref> show that the forward and reverse resistance coefficient of vortex diode with double tangential tube are both reduced compared to the traditional vortex diode, so is its impedance ratio. In other words, its overall performance is improved.</p></sec><sec id="s4_2"><title>4.2. Internal Flow Field Analysis</title><p>The pressure and speed contours of the vortex diode with double tangential tube are shown in <xref ref-type="fig" rid="fig5">Figure 5</xref> and <xref ref-type="fig" rid="fig6">Figure 6</xref>. When import velocity of two tangential tubes is 5 m/s, both tangential tubes show symmetric distribution in pressure and speed, so the internal flow field is stable. The import and export pressure are obviously different, and the maximum velocity is near the center of the chamber because of the rotation. <xref ref-type="fig" rid="fig7">Figure 7</xref> shows that a vortex structure of velocity vector is formed in the central position.</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Forward and reverse resistance coefficient of the two different structures</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2860120x10.png"/></fig><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Impedance ratio of the two different structures</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2860120x11.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Contours of static pressure</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2860120x12.png"/></fig><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Contour of velocity</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2860120x13.png"/></fig><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Vector of velocity in the chamber</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2860120x14.png"/></fig></sec><sec id="s4_3"><title>4.3. Contrast of Reverse Flow Simulation Time</title><p>Due to the stability of the internal flow field, the convergence speed has improved significantly in reverse flow. <xref ref-type="fig" rid="fig8">Figure 8</xref> shows the computation time of two</p><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> Computation time of two different structures in the same computer</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2860120x15.png"/></fig><p>different structures in the same computer.</p></sec></sec><sec id="s5"><title>5. Conclusion</title><p>At the same Reynolds number, impedance ratio of the vortex diode with double tangential tube is significantly higher than ordinary vortex diode. With one more exit tube set than ordinary vortex diode, forward flow resistance decreases significantly, and the reverse flow resistance reduces at the same time. Thus, the comprehensive performance is improved. Due to the symmetry of the structure, the internal flow field of reverse flow exhibits symmetrically, which improves the flow stability and speeds up the computing rate significantly as well. Multi-tube vortex diode structure can provide reference for improving the performance of vortex diode.</p></sec><sec id="s6"><title>Acknowledgements</title><p>Sincere thanks to the members of OJMSi for their professional performance and high quality work.</p></sec><sec id="s7"><title>Cite this paper</title><p>Zhang, L.J., Dai, Y.L., Li, Y.L., Wang, Q., Jin, X.Y., Ma, L. and Huang, W.T. (2017) Simulation Analysis of Flow Characteristics in the Vortex Diode with Double Tangential Tube. Open Journal of Modelling and Simulation, 5, 174-181. https://doi.org/10.4236/ojmsi.2017.53013</p></sec></body><back><ref-list><title>References</title><ref id="scirp.77788-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Malhotra, A., Branion, R.M.R. and Hauptmann, E.G. (1994) Modelling the Flow in a Hydrocyclone. The Canadian Journal of Chemical Engineering, 72, 953-960.  
https://doi.org/10.1002/cjce.5450720603</mixed-citation></ref><ref id="scirp.77788-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Speziale, C.G. and Thangam, S. (1992) Analysis of an RNG Based Turbulence Model for Separated Flows. International Journal Engineering Science, 30, 1379-1388.  
https://doi.org/10.1016/0020-7225(92)90148-A</mixed-citation></ref><ref id="scirp.77788-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Launder, B., Reece, G.J. and Rodi, W. (1975) Progress in the Development of a Reynolds Stress Turbulence Closure. Journal of Fluid Mechanics, 68, 537-566.  
https://doi.org/10.1017/S0022112075001814</mixed-citation></ref><ref id="scirp.77788-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Pierce, C.D. and Moin, P. (1998) Large Eddy Simulation of a Confined Coaxial Jet with Swirl and Heat Release. Proceedings of the 29th AIAA Fluid Dynamics Conference, Albuquerque, NM, 15-18 June 1998, AIAA 98-2892.  
https://doi.org/10.2514/6.1998-2892</mixed-citation></ref><ref id="scirp.77788-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Jordan, K., Zahrai, S., Macchion, O. and Vomhoff, H. (2006) Numerical Modeling of Highly Swirling Flows in a Through-flow Cylindrical Hydrocyclone. AIChE Journal, 52, 3334-3344. https://doi.org/10.1002/aic.10955</mixed-citation></ref><ref id="scirp.77788-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Kulkarni, A.A., Ranade, V.V., Rajeev, R. and Koganti, S.B. (2008) CFD Simulation of Flow in Vortex Diodes. AIChE Journal, 54, 1139-1152.  
https://doi.org/10.1002/aic.11439</mixed-citation></ref></ref-list></back></article>