手机版
您的当前位置: 好论文网 > 科技论文 > 最佳东方_最佳鳃片倾斜角下分离鳃内部流场的数值模拟

最佳东方_最佳鳃片倾斜角下分离鳃内部流场的数值模拟

来源:科技论文 时间:2018-09-13 点击: 推荐访问:流场模拟软件

  摘要:采用均匀试验设计理论对鳃片倾斜角α与β进行数值模拟方案的设计,通过采用Fluent软件中的欧拉模型,以鳃片3上表面的泥沙速度作为考核指标,对不同鳃片倾斜角下的分离鳃速度流场进行模拟。运用投影寻踪回归分析法(PPR)对数值计算结果进行了分析,得出鳃片β倾斜角对分离鳃水沙分离效果影响是最大的,其次是鳃片α倾斜角的影响;当鳃片倾斜角的组合为α=60°与β=45°时,分离鳃中的速度流场有利于水沙的快速分离且水沙分离效果最佳。
  关键词:数值模拟;鳃片倾斜角;欧拉模型;均匀试验;分离鳃;泥沙速度;仿真分析
  中图分类号:TV14 文献标志码:A 文章编号:1672-1683(2015)06-1105-05
  Abstract:The scheme design of numerical simulation on tilt angle α and tilt angle β of gill-piece applied the uniform experimental design theory.With the sediment velocity on the upper surface of gill-piece 3 as the evaluation index,velocity flow field was simulated under different tilt angles of gill-piece through the adoption of Eulerian model provided by the Fluent software.The projection pursuit regression (PPR) was used to analyze the numerical simulation results.The results showed that the influence of angle β on the effect of water-sediment separation in the gill-piece separation device is the highest,followed by the tilt angle α.The velocity field of gill-piece separation device is beneficial to increase water- sediment separation and the effect of water-sediment separation is the best under the combination of α= 60 ° and β=45 °.
  Key words:numerical simulation;tilt angle of gill-piece;Eulerian model;uniform test;gill-piece separation device;sediment velocity;simulation analysis
  我国是世界上河流泥沙问题最严重的国家之一。据观测资料统计,我国西北和华北地区的大部分河流都具有含沙量高且泥沙颗粒细的特点,水沙分离难度大。目前,采用沉沙池与过滤器来处理浑水,但沉沙池水沙分离周期长且投资大,而经过过滤器过滤后的水往往不符合规范的要求[1-2]。故研究一种新型水沙分离装置势在必行。分离鳃具有结构简单、水沙分离快、工程造价低及运行管理方便等特点,该装置对解决分散居住的农牧民生活用水问题具有积极意义,同时分离鳃作为灌区泥沙前期处理设施,可为灌区滴、喷灌提供优质水源,促进节水灌溉技术在多沙河流地区的推广应用[3]。目前,分离鳃已通过物理模型试验及数值模拟的研究方法得到了以下研究成果:(1)分离鳃中泥沙通道淤堵的临界宽度是1 mm[4];(2)分离鳃在处理含有黏性泥沙的浑水时,能够体现出优势[5];(3)鳃片间距为d=5 cm时,水沙分离效果最佳[6];(4)含沙量为10~80 kg/m 时,分离鳃的水沙分离效果更佳[7];(5)Fluent软件中的欧拉模型适合模拟分离鳃中的水沙两相流流场[8]。鳃片倾斜角对分离鳃水沙分离效果的影响不容忽视,而采用数值模拟的方法模拟静水中不同鳃片倾斜角下的分离鳃内部流场是一种非常有效的方法,为减少计算次数,数值模拟方案设计采用均匀试验设计理论,并运用投影寻踪回归(PPR)的方法对数值计算结果进行分析与讨论,来确定最佳鳃片倾斜角。
  1 分离鳃结构及静水沉降工作原理
  从分离鳃结构图(图1)可看出该装置由鳃片与鳃管构成,从分离鳃最上端至最下端为鳃片1至鳃片5。静水沉降时,分离鳃的上端为开口状,下端呈封闭状;而在动水沉降时,把多个分离鳃布置在沉淀池中,此时分离鳃上下端均为开口状,以便清水的上升和泥沙的排出,并在分离鳃长度方向的壁面上布置小孔,以便分离鳃内外的浑水交换。采用的分离鳃特征尺寸如下:长度a=10 cm、宽度b=5 cm、高度c=40 cm、鳃片间距d=5 cm、清水流上升通道宽度e=1 cm及泥沙流下降通道宽度f=1 cm,而鳃片与鳃管长度侧壁构成的α倾斜角及鳃片与鳃管宽度侧壁构成的β倾斜角的数值见表1。
  图2表示静水条件时分离鳃的工作原理。将浑水注入到分离鳃中,泥沙颗粒在重力作用下先自由沉降,当沉降至鳃片上表面后便汇集成泥沙流,其沿着鳃片上表面的高端流动至低端,而清水流沿着鳃片的下表面由低端流动至高端,从而鳃片上表面的泥沙流运动与鳃片下表面的清水流运动在鳃片间形成一个逆时针方向的横向异重流,见图2中带箭头的虚线;当各鳃片低端的泥沙流滑落至泥沙通道时便一起垂直向下运动,而各鳃片高端的清水流运动至清水通道时便一起垂直向上运动,从而形成沿分离鳃边壁的顺时针垂向异重流,见图2中带箭头的实线,从而达到水沙分离的目的。   2 数值模拟方案设计
  2.1 数值模拟的均匀试验设计
  采用Fluent软件中的欧拉模型模拟不同倾斜角下的分离鳃中的水沙两相流流场,欧拉模型的各方程见文献[9]。分离鳃的边界条件包含2个,其设置见文献[8]。初始条件如下:将水定义为主相,密度ρw=998.2 kg/m3;泥沙颗粒定义为第2相,假定泥沙颗粒为球形,平均粒径为0.015 mm(物理模型试验时采用的中值粒径),密度ρs=2 650 kg/m3,考虑重力作用。初始化时,假定分离鳃内水沙混合均匀,含沙量为S=60 kg/m 。
  均匀设计可以用较少的试验次数,安排多因素、多水平的析因试验,是在均匀性的度量下最好的析因试验设计方法。它可以使试验点在试验范围内充分地均匀分散,不仅可大大减少试验次数,而且仍能得到反映试验体系主要特征的试验结果。均匀设计也是电脑仿真试验设计的重要方法之一,从而可用较少的试验次数获得比较满意的结果[10]。鳃片倾斜角α与鳃片倾斜角β对分离鳃的水沙分离效果影响很大,根据这2个影响因素的特征,每个因素对应取6个水平,详见表1。2.2 考核指标
  鳃片上表面的泥沙速度大小代表泥沙从浑水中分离出来的快慢,速度越大,说明泥沙从浑水中脱离的越快,水沙分离效果越好,反之,速度越小水沙分离效果就越差。各鳃片上表面的泥沙速度相差不大[6],可取鳃片3上的泥沙速度作为分离鳃水沙分离效果的考核指标。进行数值计算前,在Fluent软件中设置监测鳃片3上表面的泥沙速度与时间的关系,以获得不同时刻的泥沙速度。
  2.3 数值计算结果
  因选取的是2因素6水平,故数值模拟的均匀试验可选用U*6(64)均匀设计表。每个均匀设计表都附有一个使用表,从U*6(64)的使用表(表2)可知,2个因素时,应选用1,3两列来安排试验,则均匀度的偏差D=0.187 5,该值表示均匀分散性较好。均匀试验的数值计算方案及数值计算结果见表3,其中V表示时间t=110 s时鳃片3上表面的泥沙速度值。
  3 PPR分析
  PPR是用于处理和分析非正态非线性总体高维数据的一种方法,其不需要对试验数据作任何人为假定、分割或变量处理,是统计学、应用数学及计算机技术三者交叉的前沿领域,已成功运用到众多学科和相关部门,正日益受到人们的重视[11-14]。该方法具有以下3个特点:(1)能有效克服“维数灾难”;(2)客观充分地挖掘数据的各种非正态非线性规律,用数值函数加以记忆;(3)模型具有稳定性。
  设y是因变量,x是p维自变量,PPR 模型如式(1)
  =E(y|x)=y+∑MUi=1δifi(ζTix)(1)
  式中:MU为数值函数最优个数;δi为第i个数值函数的相应权重;fi为数值函数; ζTi为i方向的投影值,其中‖ζTi‖=1,i=1,2,…,MU;ζTi=(ζi1,ζi2,…,ζip)。
  运用PPR方法将α倾斜角与β倾斜角这两个影响因素综合在一起考虑,对表2中的6组数值计算结果进行PPR分析,PPR的算法及操作程序见文献[15]。反映投影灵敏度指标的光滑系数取0.90,投影数值函数个数初始值M=3,最终投影数值函数个数取MU=2。模型参数为:N(要求数据大于9),P,Q,M,MU=6×2,2,1,3,2。PPR建模结果为:δ1=1.0013,ζT1=(0.2773,0.9608);δ2=0.0290,ζT2=(-0.9742,-0.2254)。不同倾斜角下鳃片3上表面的拟合值、计算值以及两者的绝对误差见表4。由表4可知数值模拟的计算值与拟合值两者吻合较好,泥沙速度的绝对误差≤±0.52 cm/s,相对误差在3%左右;但鳃片倾斜角α=15°,β=35°时,相对误差为24.30%,鳃片倾斜角α=55°,β=15°时,相对误差为36.67%,分析原因是网格数量及质量造成的,但一定程度上PPR模型可以反映α倾斜角与β倾斜角对分离鳃水沙分离效果的影响。
  4 单因子分析
  利用鳃片3上表面泥沙速度的PPR模型可以进行各种类型仿真。把各因子按其变幅分成11级水平,每次考虑变动一个因子,其它因子保持在中水平,即将α倾斜角与β倾斜角都分成如下的11级水平:15°,20°,25°,30°,35°,40°,45°,50°,55°,60°,65°,每次变动α倾斜角(或β倾斜角),而β=40°(或α=40°),可知单因子α倾斜角(或β倾斜角)对鳃片上表面泥沙速度的影响,用PPR仿真后的数据绘制各因子对泥沙最大速度的影响,见图3。从图3中可直观地看出两个因子对泥沙速度影响大小的排序:β倾斜角大于α倾斜角(因β倾斜角的曲线斜率大于α倾斜角的斜率),这同PPR分析的因子贡献权重是一致的,即β倾斜角的贡献权重为1.000 00,α倾斜角的贡献权重为0.260 43。从图3中还可看出α倾斜角为55°~60°时,α倾斜角对泥沙速度影响最大,且在α=60°时泥沙速度达到峰值,即该倾斜角下分离鳃的水沙分离效果最好。β倾斜角为15°~25°时,泥沙最大速度随β倾斜角的增大而增幅很少,V<1 cm/s;β倾斜角为25°~60°时,随着β倾斜角的增大泥沙速度快速地增大;而β角在60°~65° 时可发现泥沙速度相差不大。β倾斜角对分离鳃的水沙分离效果影响包括两个方面。一方面,倾斜角过小,则泥沙的下滑力不能克服水流的阻力与摩擦力,不利于泥沙的下滑,且泥沙速度较小;另一方面,当倾斜角过大,虽然可以缩短泥沙颗粒在鳃片上的聚结时间,加快泥沙沉降,但同时随着β倾斜角增大,鳃片的长度也会增加,则泥沙颗粒与清水滑动的路程相应地也增加,此时有效鳃片数也会减少,相应的沉淀面积也会减少,对水沙分离过程起着负面作用,因此,存在一个合理的β倾斜角。所以,单因子分析也为进一步进行的因子全局优化仿真提供了参考信息。
  5 多因子优化仿真   综合考虑α倾斜角与β倾斜角同泥沙速度的关系,将PPR仿真的数据,利用Surfer8.0软件绘制出泥沙速度等值线图。由图4可知当β在45°~65°,α在30°~65°取值时,泥沙速度在4.22 cm/s以上,据文献[6]可知,鳃片间距对分离鳃的水沙分离效果影响最大,故为了能在分离鳃中多布置鳃片,减少鳃片间距,以加大沉淀面积,最佳倾斜角的组合为α=60°,β=45°,该倾斜角下对应的泥沙速度为5.69 cm/s。
  6 最佳鳃片倾斜角下的速度场
  6.1 鳃片上、下表面的速度场
  图5表示计算时间为110 s时,倾斜角α=60°,β=45°下的鳃片1、3及5的上、下表面速度矢量分布图。从图5(a)可知,鳃片1、2及3上表面的流动主要为泥沙流运动,泥沙流先沿鳃片斜面的最大倾斜角63.43°向下滑动,当滑至鳃片长边时同其它泥沙流一起以β=45°向下运动,但鳃片2与3上表面的泥沙流不仅受到α与β倾斜角的影响,还受到来自泥沙通道中下沉的泥沙流影响,使得这两个鳃片上表面右上角的速度呈向内放射状,而鳃片1靠近分离鳃上端,其上表面的泥沙流仅受α与β倾斜角的影响,则无此特征,故鳃片1上表面的速度分布规律与其它两个鳃片上表面不同。从图5(b)可知,鳃片1、2及3下表面的流动主要为清水流运动,清水先沿鳃片斜面的最大倾斜角63.43°向上滑动,当流至鳃片长边时同其它清水流以β=45°向上流动,但鳃片1与2下表面的清水流不仅受到α与β倾斜角的影响,还受到来自清水通道中上升的清水流影响,使得鳃片1与2下表面左下角的速度呈向内放射状,而鳃片3 靠近分离鳃下端,其仅受α与β倾斜角的影响,则无此特征,故鳃片3下表面的速度分布规律与其它两个鳃片下表面不同。
  6.2 宽方向的速度场
  图6表示计算时间为110 s时,倾斜角α=60°,β=45°下的宽方向(Z断面)的速度矢量分布图。从图6(a)中可知,5个鳃片下表面的水流从鳃片的低端以β=45°沿着鳃片长边滑至鳃片高端,后汇入清水通道向上运动。从图6(c)可知,5个鳃片上表面的水流从鳃片的高端以β=45°沿着鳃片长边滑至鳃片低端,后汇入泥沙通道向下运动。图6(a)与图6(c)的速度分布规律在图6(b)中得到了体现,且从图6(b)中还可看出鳃片上表面的水流、鳃片下表面的水流、泥沙通道中的水流及清水通道中的水流相互干扰较少,各自按着相应的运动轨迹流动;鳃片间的局部涡流对水流和泥沙具有卷吸作用,加速了泥沙的沉降和清水的上升,从而达到水沙分离的目的。
  7 结论
  (1) PPR模型能够较好地反映鳃片的α与β倾斜角对分离鳃水沙分离效果的影响。
  (2)影响分离鳃水沙分离效果的各因素主次排列为鳃片的β倾斜角,鳃片的α倾斜角。
  (3)当α=60°,β=45°时,分离鳃的水沙分离效果最佳。
  (4)鳃片倾斜角为α=60°,β=45°时的速度流场有利于分离鳃中的水沙快速分离。
  参考文献(References):
  [1]韩丙芳,田军仓.微灌用高含沙水处理技术研究综述[J].宁夏农学院学报,2001,22(2):63-69.(HAN Bing-fang,TIAN Jun-cang.General introduction on the technology of purifying water for micro-irrigation with high-silt content water[J].Journal of Ningxia Agricultural College,2001,22(2):63-69.(in Chinese))
  [2]李琳,邱秀云,聂境,等.基于滴灌浑水水力分离装置的水沙分离试验研究[J].水资源与水工程学报,2010,21(6):37-40.(LI Lin,QIU Xiu-yun,NIE Jin,et al.Experiment on the water-sediment separation by turbid water hydraulic separation device Based on drip irrigation[J].Journal of Water Resources & Water Engineering,2010,21(6):37-40.(in Chinese))
  [3]邱秀云,龚守远,严跃成,等.一种新型水沙分离装置的研究[J].新疆农业大学学报,2007,30(1):68-70.(QIU Xiu-yun,GONG Shou-yuan,YAN Yao-cheng,et al.The research on a new separator device for water-sediment[J].Journal of Xinjiang Agricultural University,2007,30(1):68-70.(in Chinese))
  [4]严跃成,邱秀云,张翔,等.两相流分离鳃泥沙下沉通道宽度对水沙分离影响的试验研究[J].新疆农业大学学报,2011,34(6):526-528.(YAN Yao-cheng,QIU Xiu-yun,ZHANG Xiang,et al.Experimental study on influence of sediment settling access witdh in two-phase flow separation device on water-sediment separation.Journal of Xinjiang Agricultural University,2011,34(6):526-528.(in Chinese))
  [5]朱超,邱秀云.垂向异重流式分离鳃适用泥沙的试验研究[J].人民长江,2009,40(5):60-61.(ZHU Chao,QIU Xiu-yu.Experimental study of sediment on vertical density flow gill-piece separation device[J].Yangtze River,2009,40(5):60-61.(in Chinese))   [6]陶洪飞,邱秀云,李巧,等.不同鳃片间距下的分离鳃内部流场三维数值模拟[J].农业机械学报,2014,45(6):183-189.(TAO Hong-fei,QIU Xiu-yun,LI Qiao,et al.3-D numerical simulation of internal flow field in gill-piece separation device under different gill-piece spacing[J].Transactions of the Chinese Society for Agricultural Machinery,2014(6):183-189.(in Chinese))
  [7]陶洪飞,邱秀云,何照青,等.含沙量对分离鳃的水沙分离影响试验[J].水电能源科学,2013,31(7):93-95.(TAO Hong-fei,QIU Xiu-yun,HE Zhao-qing,et al.Sediment concentration on water-sediment separation effect test of separation device[J].Water Resources and Power,2013,31(7):93-95.(in Chinese))
  [8]陶洪飞,邱秀云,赵丽娜,等.水沙分离鳃内部流场数值模拟与验证[J].农业工程学报,2013,29(17):38-46.(TAO Hong-fei,QIU Xiu-yun,ZHAO Li-na,et al.Numerical simulation and PIV display of internal flow field in gill-piece separation device[J].Transactions of the Chinese Society of Agricultural Engineering,2013,29(17):38-46.(in Chinese))
  [9]江帆,黄鹏.FLUENT高级应用与实例分析[M].北京:清华大学出版社,2010.(JIANG Fan,HUANG Peng.FLUENT Advanced Applications and Case Analysis[M].Beijing:Tsinghua Press,2010.(in Chinese))
  [10]任露泉.试验优化设计与分析[M].第二版.北京:高等教育出版社,2003:193-207.(REN Lu-quan.Optimum design and analysis of experiments[M].Second Edition.Beijing:Higher Education Press,2003:193-207 (in Chinese))
  [11]侯杰,赵涛,牧振伟,等.悬栅消能率的PPR因子贡献率分析及优化仿真[J].水力发电,2005,31(2):38-40.(HOU Jie,ZHAO Tao,MU Zhen-wei,et al.Analysis and optimization simulation of PPR factors contribution ratio of suspension girder dumping ratio[J].Water Power,2005,31(2):38-40.(in Chinese))
  [12]冯勇,何建新,刘亮,等.冻融循环作用下细粒土抗剪强度特性试验研究[J].冰川冻土,2008,30(6):1013-1017.(FENG Yong,HE Jian-xin,LIU Liang,et al.Experimental Study of the Shear Strength Characteristics of Fine-Grained Soil under Freezing and Thawing Cycles[J].Journal of glaciology and geocryology,2008,30(6):1013-1017.(in Chinese))
  [13]何建新,郭鹏飞,刘录录,等.阳离子乳化沥青混凝土配合比设计的优选方法研究[J].水利与建筑工程学报, 2013,11(3):96-98.(HE Jian-xin,GUO Peng-fei,LIU Lu-lu,et al.Study on Optimal Method for Mix Design of Cation Emulsified Asphalt Concrete[J].Journal of Water Resources and Architectural Engineering,2013,11(3):96-98.(in Chinese))
  [14]吴素芬,韩萍,李燕,等.塔里木河源流水资源变化趋势预测[J].冰川冻土,2003,25(6):708-711.(WU Shu-fen,HAN Ping,LI Yan,et al.Predicted Variation Tendency of the Water Resources in the Headwaters of the Tarim River[J].Journal of glaciology and geocryology,2003,25(6):708-711.(in Chinese))
  [15]周新民,李嘉,赵文谦,等.确定固液两相流扩散系数的投影寻踪模型[J].水动力学研究与进展:A辑,2001,16(2):175-180.(ZHOU Xin-min,LI Jia,ZHAO Wen-qian,et al.Project pursuit regression model for diffusion coefficient of suspended particles in liquid-solid two-phase flow[J].Journal of Hydrodynamics,Ser A,2001,16(2):175-180.(in Chinese))

好论文网 www.ho59.com

Copyright © 2002-2018 . 好论文网 版权所有 京ICP备10015900号

Top