Message Board

Dear readers, authors and reviewers,you can add a message on this page. We will reply to you as soon as possible!

2021 Volume 43 Issue 8
Article Contents

REN Pu, LONG Wei. Numerical Analysis of Flow Field Characteristics and Hydrodynamic Cavitation Mechanism of Tesla Valves[J]. Journal of Southwest University Natural Science Edition, 2021, 43(8): 161-166. doi: 10.13718/j.cnki.xdzk.2021.08.020
Citation: REN Pu, LONG Wei. Numerical Analysis of Flow Field Characteristics and Hydrodynamic Cavitation Mechanism of Tesla Valves[J]. Journal of Southwest University Natural Science Edition, 2021, 43(8): 161-166. doi: 10.13718/j.cnki.xdzk.2021.08.020

Numerical Analysis of Flow Field Characteristics and Hydrodynamic Cavitation Mechanism of Tesla Valves

More Information
  • Corresponding author: LONG Wei
  • Received Date: 09/04/2021
    Available Online: 20/08/2021
  • MSC: O359+.1

  • In this study, based on the Fluent simulation platform, from the perspective of fluid dynamics, the flow field at the crossing exit of the Tesla valve reverse guide is divided into five regions, and the flow field characteristics are analyzed. The cause of the cavitation phenomenon and the influence of the cavitation model on flow field characteristics are analyzed. The results show that there are two low-speed fluid regions, one double-vortex region, one cavitation region and one high-speed fluid region in the channel when the reverse direction of the Tesla valve is running. At the same time, without the cavitation model, the minimum flow field pressure is -3.66×105 Pa, and after adding the cavitation model, the minimum pressure is -9.92×104 Pa. With the addition of the cavitation model, the vortex at the cross exit is smaller and closer to the exit region. At the same time, the reduction of eddy current and its left shift can greatly reduce the intensity of turbulent kinetic energy at the left side of the cross exit, and the peak value of turbulent kinetic energy is 5.6×10 m2/s2.
  • 加载中
  • [1] KUBAR A A, CHENG J, KUMAR S, et al. Strengthening Mass Transfer with the Tesla-Valve Baffles to Increase the Biomass Yield of Arthrospira Platensis in a Column Photobioreactor[J]. Bioresource Technology, 2021, 320: 124337. doi: 10.1016/j.biortech.2020.124337

    CrossRef Google Scholar

    [2] ARUNACHALA U C, RAJAT A C, SHAH D, et al. Stability Improvement in Natural Circulation Loop Using Tesla Valve-An Experimental Investigation[J]. International Journal of Mechanical and Production Engineering Research and Development, 2019, 9(6): 13-24. doi: 10.24247/ijmperddec20192

    CrossRef Google Scholar

    [3] QIAN J Y, WU J Y, GAO Z X, et al. Hydrogen Decompression Analysis by Multi-Stage Tesla Valves for Hydrogen Fuel cell[J]. International Journal of Hydrogen Energy, 2019, 44(26): 13666-13674. doi: 10.1016/j.ijhydene.2019.03.235

    CrossRef Google Scholar

    [4] QIAN J Y, CHEN M R, GAO Z X, et al. Mach Number and Energy Loss Analysis Inside Multi-Stage Tesla Valves for Hydrogen Decompression[J]. Energy, 2019, 179: 647-654. doi: 10.1016/j.energy.2019.05.064

    CrossRef Google Scholar

    [5] WANG C T, CHEN Y M, HONG P A, et al. Tesla Valvesin Micromixers[J]. International Journal of Chemical Reactor Engineering, 2014, 12(1): 397-403. doi: 10.1515/ijcre-2013-0106

    CrossRef Google Scholar

    [6] 冯鲁文, 程泽堃. 特斯拉阀门的优化与改进仿真[J]. 技术与市场, 2019, 26(10): 33-34, 37.

    Google Scholar

    [7] PORWAL, PIYUSH. Thermal and Fluidic Characterization of Tesla Valve Via Computational Fluid Dynamics[D]. Mississippi State University, 2016.

    Google Scholar

    [8] 刘海洋, 马佳, 王森, 等. 特斯拉阀单向导通性研究[J]. 物理与工程, 2020, 30(1): 120-124.

    Google Scholar

    [9] 周润中, 乔宇杰, 张钰翔, 等. 特斯拉阀性能的仿真研究[J]. 物理实验, 2020, 40(9): 44-50.

    Google Scholar

    [10] 张之禾, 陈真, 陈婷. 基于有限元仿真与3D打印技术的特斯拉阀结构研究[J]. 大学物理, 2020, 39(4): 73-77.

    Google Scholar

    [11] THOMPSON S M, PAUDEL B J, JAMAL T, et al. Numerical Investigation of Multistaged Tesla Valves[J]. Journal of Fluids Engineering, Transactions of the ASME, 2014, 136(8): 081102. doi: 10.1115/1.4026620

    CrossRef Google Scholar

    [12] PORWAL P R, THOMPSON S M, WALTERS D K, et al. Heat Transfer and Fluid Flow Characteristics in Multistaged Tesla Valves[J]. Numerical Heat Transfer, Part A: Applications. An International Journal of Computation and Methodology, 2018, 73(6): 347-365.

    Google Scholar

    [13] QIAN J Y, CHEN M R, LIU X L, et al. A Numerical Investigation of the Flow of Nanofluids Through a Micro Tesla Valve[J]. Journal of Zhejiang University: Science A, 2019, 20(1): 50-60. doi: 10.1631/jzus.A1800431

    CrossRef Google Scholar

    [14] 钱锦远, 陈珉芮, 李晓娟, 等. 特斯拉阀型微通道热沉换热强化的数值研究[J]. 流体机械, 2020, 48(12): 36-42. doi: 10.3969/j.issn.1005-0329.2020.12.006

    CrossRef Google Scholar

    [15] 王涛, 王昊文, 林煜琛, 等. 特斯拉阀水力空化的数值研究[J]. 高校化学工程学报, 2020, 34(4): 884-889. doi: 10.3969/j.issn.1003-9015.2020.04.005

    CrossRef Google Scholar

  • 加载中
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(9)  /  Tables(1)

Article Metrics

Article views(2124) PDF downloads(1042) Cited by(0)

Access History

Other Articles By Authors

Numerical Analysis of Flow Field Characteristics and Hydrodynamic Cavitation Mechanism of Tesla Valves

    Corresponding author: LONG Wei

Abstract: In this study, based on the Fluent simulation platform, from the perspective of fluid dynamics, the flow field at the crossing exit of the Tesla valve reverse guide is divided into five regions, and the flow field characteristics are analyzed. The cause of the cavitation phenomenon and the influence of the cavitation model on flow field characteristics are analyzed. The results show that there are two low-speed fluid regions, one double-vortex region, one cavitation region and one high-speed fluid region in the channel when the reverse direction of the Tesla valve is running. At the same time, without the cavitation model, the minimum flow field pressure is -3.66×105 Pa, and after adding the cavitation model, the minimum pressure is -9.92×104 Pa. With the addition of the cavitation model, the vortex at the cross exit is smaller and closer to the exit region. At the same time, the reduction of eddy current and its left shift can greatly reduce the intensity of turbulent kinetic energy at the left side of the cross exit, and the peak value of turbulent kinetic energy is 5.6×10 m2/s2.

  • 特斯拉阀因其特殊的结构和控制流动特性在微流体控制、微机电系统和生物化学等领域得到广泛应用[1-5].

    近年来,国际上许多学者对特斯拉阀的结构和导通性能开展了大量研究[6-10]. Thompson等[11]利用三维CFD技术,对多级特斯拉阀(multi-stage Tesla valve,MSTV)的有效性进行了量化,并研究了低雷诺数和不同阀阀距离下多级特斯拉阀的有效性. Porwal等[12]通过数值分析研究了进口雷诺数(25~200)对单个特期拉阀和MSTV(10级)整流和热增强能力的影响. Qian等[13]通过数值分析对Al2O3-水纳米流体微尺度T45-R型特斯拉阀的流动特性进行了研究. 钱锦远等[14]采用数值模拟的方法,对不同参数和级数的特斯拉阀型微通道热沉的强化换热进行了研究.

    但目前对于特斯拉阀的研究多倾向于结构和正向导通时的研究,对于特斯拉阀逆向导通时的空化现象[15]及其形成机理的研究还不多见. 因此,本文从特斯拉阀逆向导通时流场特性入手,对特斯拉阀交叉出口处流场特性进行详细分区分析,研究逆向导通时特斯拉阀发生空化的原因,同时对比加入空化模型前后,空化效应对流场特性的影响.

1.   模型的建立
  • 本文以图 1(a)所示结构的特斯拉阀为研究对象,结构参数如图 1(b)所示. 通道截面形状参数为10 mm×10 mm的矩形.

  • 当流体中局部压力低于流体的空气分压力时,流体中融入的气体会析出,发生空化现象. 特斯拉阀反向导通时,阀体的交叉出口处会产生局部低压,当阀体进出口压差较大时,交叉出口处就发生空化. 采用mixture模型中Zwart-Gerber-Belamri(ZGB)空化模型对流场进行表达计算,液相-汽相质量传递由蒸汽传输方程表示如式(1)所示:

    式中v表示气相,α为蒸汽体积分数,ρv为蒸汽密度,vv为气相速度,ReRc分别为与蒸汽泡生长、破裂相关的传质源项. 在ZGB模型中,气相和液相之间的总相间传质率R

    式中RB为气泡半径,取值10-6m; PB为气泡压力; P为流体当前压力; ρl为液相密度. 建立空化模型如下式所示:

    PPv,则

    PPv,则

    式中Fvap为蒸发系数,取值50;Fcond为冷凝系数,取值0.01.

  • 以六面体网格为主对几何模型进行网格划分,网格尺寸限制为0.4 mm,模型总网格数为57万,如图 2所示. 边界条件设置如表 1所示.

2.   流场特性分析
  • 利用Fluent平台,以图 1所示的模型为研究对象,反向导通条件下对阀内流场进行分析. 采用多相流欧拉空化(mixture)和标准k-ε湍流方程相结合的模型,设定流体密度为998.2 kg/m3,动力黏度为0.001 003 Pa·s,空化气相(vapor)密度为0.554 2 kg/m3,动力黏度为1.34×10-5 Pa·s,空化模型选用Zwart-Gerber-Belamri模型,具有较高的计算精度以及较好的收敛性,空化压力为4 133 Pa,压力-速度耦合器选择Coupled算法,计算结果收敛残差小于10-6.

  • 首先对特斯拉阀反向导通时流场特性进行分析,明确流场内部的高速、低速和高压、低压区域的形成原因.

    图 3(a)速度云图可知,特斯拉阀逆向导通时,直通道中共有1、2两处低速流体区域.

    区域1(射流末端区域)的流体流速较低是由于弯管出口处的高速射流与交汇处直管中流出的流体发生冲击,产生的速度梯度较大,速度迅速降低,且由于两个方向流体的冲击作用,弯管出口处的高速流体产生沿垂直方向的转向.

    区域2流体速度较低,一方面是由于分叉入口进入区域2的流体较少,流速较低,另一方面是由于弯管中的流体高速射流流出弯管后与直管中水平流动的流体发生冲击,阻碍了区域2中流体的流动,因此区域2中流体流速较低,压力较大.

    区域3的流体流速较高是由于此区域流体冲击后发生转向,流速较高,且由于射流冲击和流体流向的转变,在区域3中形成如图 4所示的2个涡流.

    受区域1和区域3共同影响,弯管中的高速射流绝大部分在1、3区域左侧发生垂直转向流入左侧水平通道中,即图 3(b)中的区域4. 区域4的流体流速最大,一方面是弯管中的高速射流转向后流入,另一方面是区域4右侧均为高压区域,左侧联通出口为标准大气压,形成较大压差.

    图 5压力云图可知,弯管-直管交叉处右侧压力较大,左侧压力较低,且区域4、5最低压力低于液态水的空气分离压4 133 Pa. 其中区域4压力较低是由于流体流量和流速较大.

  • 由水力空化效应可知,当通道内流体压强降低至空气分离压时,流体内部气核析出形成气泡. 由于区域5最低压力低于液态水的空气分离压,因此在计算模型中加入多相流的欧拉空化模型(mixture)对计算模型进行修正. 数值分析发现:如果忽略空化效应,会在流场内部形成过低负压与实际情况不符.

    图 6所示为未考虑空化时的流场压力,由于计算模型中没有空化模型,压力会随流体流速增大不断下降,压力最低值为-3.66×105 Pa. 而图 5中,考虑空化时,由于空化效应的存在,流体发生空化,压力最小值为-9.92×104 Pa. 因此,在对特斯拉阀逆向导通过程进行数值分析时,考虑空化效应的影响是必要的.

    图 7可知,空化效应集中发生在低压区域5. 由上述分析可知,弯管-直管交叉处右侧压力较大,左侧压力较低,形成巨大压差; 由图 3(c)速度矢量图可知,直管中的少量流体和弯管中绝大部分的高速流体转向之后集中流向区域4. 当入口压力为0.5Mp时,区域4的最高速度可达30.9 m/s,由于流体本身流量、流速很大,加之较大的压差,使汇聚后的流体高速流过区域4,因此会在区域5形成低压区域,当压力降低至流体空气分力压时,区域5流体发生空化,且由于特斯拉阀出口处压力为标准大气压,因此空化区域并未延伸至出口处.

    对比空化模型对交叉出口之后的流场特性的影响,图 8所示分别为不考虑空化和考虑空化的后半段流场特性.

    图 8(a)可知,无空化模型时,流场内部不会发生空化,区域5因区域4的高速转向流体所产生的牛顿内摩擦力形成涡流,涡流强度取决于转向流体的速度大小. 加入空化模型后,由上述分析可知会将区域5流体中气体析出,产生不同程度的空化现象,因此区域5几乎不会产生大的涡流,而在区域5之后由于流体流速减小流体垂直扩散而产生较小的涡流,涡流更接近出口区域.

  • 由流场特性和空化效应分析可知,考虑空化效应时,交叉出口左侧区域流场会由于空化效应导致涡流更小更接近出口区域,对特斯拉阀出口处流体湍动能和空化区域温度将产生影响.

    图 9湍动能的对比,结合图 8中流场迹线涡流位置及大小可知,空化模型的加入对特斯拉阀弯管-直管交叉处左侧区域湍动能产生极大的影响,无空化模型中,交叉出口左侧为涡流产生区域,此处湍动能极大,峰值为8.4×104 m2/s2; 由于空化效应的存在,交叉出口左侧空化区域发生空化,涡流更小且更接近出口区域,因此湍动能极小且发生湍流的主要区域也更接近出口区域,湍动能峰值为5.6×10 m2/s2.

3.   结论
  • 1) 基于Fluent仿真平台,采用多相流欧拉空化和标准k-ε湍流方程相结合的模型,从流体动力学角度出发,将特斯拉阀反向导通时交叉出口处的流场分为5个区域,通道内存在2个低速流体区域,1个双涡流区域,1个空化区域和1个高速流体区域.

    2) 对比有无空化模型的影响,无空化模型流场压力最低值为-3.66×105 Pa,通道内部压力峰值为4.6×105 Pa; 加入空化模型后,压力最低值为-9.92×104 Pa,压力峰值为4.6×105 Pa,空化模型对压力峰值无明显影响,但会显著提升最低压力值; 加入空化模型,交叉出口左侧的涡流更小且更接近出口区域; 同时涡流的减小与左移,可以极大减小交叉出口左侧湍流动能强度,湍动能峰值为5.6×10 m2/s2.

Figure (9)  Table (1) Reference (15)

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return