1、孑沟槽过嫌哉廊Mf丿丿刃J图2CFD模型表1三种沟槽表面上网格点数列表形状顶角sh沟槽表面网格总数沟槽表面积网格数三角形900.15x 10-2141810234.53 X 10-6扇贝形-5X 10-2161978214.95 X 10-6刀刃形-25X 1021240005-66.00 X 10在ANSYS中建立的几何模型如图2( a)所示,计算域展向长度为0. 8mm,流 向长度为4mm,垂向高度为4mm图2( b) 2( f ) 为离散后的CFD模型,扇贝形和三角形离散的六面体网格图相似,在此只列出三角形的CFD模型,沟槽表 面所划分的网格点数列于表1中,刀刃形沟槽表面形状比较特殊,很
2、难划分六 面体网格,只能划分四面体网格。2.2边界条件流场分析时认为流场稳定,各参数不随时间变化;流体为不可压缩流,密 度和黏性等物理性质不随时间变化;不涉及传热问题;采用各向异性假设。用有 限体积法求解N-S方程,近壁区采用B-L两层湍流模型,远离壁面区采用雷诺应 力湍流模型(RSM),边界条件为光滑表面与沟槽表面施加固壁无滑移条件,即 上下壁面X , Y, Z 三个方向速度为零,入口处给定质量流率m = 0. 2 kg/ s, 流 向和垂直于流向的方向施加周期边界条件。 计算温度t = 20C,流体介质为水,动 力黏度卩=0. 001N x s/ m2 ,密度 p = 998 kg / m
3、, 运动黏度 u = 1X 10-6 m / s, 雷诺数Re= 6.0 x 105 ,雷诺数计算时取过流断面水力直径为特征长度。3三种仿生非光滑沟槽表面流场分析3.1剪应力分析计算结果表明三种沟槽表面剪应力均小于光滑表面,具有减阻效果,但由 于沟槽形状的不同减阻效果不同,三角形沟槽减阻效果最差,刀刃形沟槽减阻 效果最佳。I I.4&+04 |P35e+O4 1. 27e+CM 1 無祕B L:1.06e+049.7, 7,館网3 6. 37c+O3 5.67c+0314.脸心I -L烁卜昭 .:kI m詣 2. .I2e+O3M . 42c+0:J o. OfeKOO-MII-36e+O41
4、,舐卜041.L ff7e+04I.Ole+OI 9.45e-K13 8.皈心 19e+O37. 56e*(KJ 一 fi, 93cO3 16. 30f+035- 67403 5P 04e*034. 41 M31 76e+033. ifie+OS 込 5203L26e*03 6.0. IWe+OO图3沟槽表面应力云图光滑表面单元平均剪应力为1.2 X 104 N/ m2。由沟槽表面剪应力云图3可以看出三角形沟槽表面图4 半个沟槽表面节点剪应力随沟槽高度变化曲线图最大剪应力为1.07 X 104 N/ m 2 , 比光滑表面减小10.8%;扇贝形沟槽表面最大剪应力为5.1 X 103 N/m2
5、,比光滑表面减小57.5% ;刀刃形沟槽表面最大剪应 力为4. 5 X 103 N/ m2 ,比光滑表面减小62.5%。由图3 4可见三种沟槽表面 高应力区出现位置及面积不同,三角形沟槽高应力区起始位置最低并且面积最 大,壁面总阻力最大,减阻效果最不明显;刀刃形沟槽高应力区出现在刀刃尖 端非常窄的带形区域内,起始位置最高并且面积最小,壁面阻力最小,减阻效 果最好。a二侍跻rrna h隔贝刑商憎衷面图5三种沟槽表面及光滑表面展向及垂向速度矢量图3. 2速度场分析图5给出了三角形、扇贝形和刀刃形沟槽表面和光滑表面展向和垂向速度矢量图,通过对比可见沟槽表面对湍流结构产生影响。 三种沟槽表面流向涡的展
6、向 运动都产生分离并生成二次涡,但是二次涡的位置、近壁区流体上冲和下扫运动 的强度不同。Choi曾指出沟槽减阻与否与沟槽面漩涡有很大关系12。三角形沟 槽表面二次涡位置最低,壁面受高强度流体上冲和下扫运动强度最大,动量传 输最强,故壁面受高应力影响区域的面积最大,由于沟槽谷底间距及壁面曲线 斜率最小,谷底附近流体受壁面阻碍并且安静流区域面积最小 ,所以层流底层厚度最小,壁面总阻力最大,减阻效果最差。刀刃形沟槽表面二次涡停留在沟槽 尖顶附近,位置最高,与壁面相接触的流体运动强度最弱,壁面几乎不受高强 度流体上冲和下扫运动影响,壁面流体的动量传输较低,由于沟槽谷底间距及 壁面曲线斜率最大,谷底附近
7、安静流体的运动比较自由并且区域最大 ,所以层流底层厚度最大,整个非光滑表面剪应力最低,总阻力最小,减阻效果最好。-8i IB!V V 00 6 4 2 0-27 6 瞎工点齐-Jr岸-r/rtnb刀刃形图6雷诺剪应力曲线图&63. 3雷诺应力分析由图6雷诺剪应力曲线可知三种沟槽表面的雷诺应力均小于光滑表面。 光滑表面、三角形、扇贝形和刀刃形沟槽雷诺应力最大值分别为 8 m2 /s 2、7. 8 m 2 /s2、7 m2 / s 2和6. 2 m 2 / s 2。雷诺应力是横向脉动在剪切场中的体现,横向脉 动体现了平均剪切率产生的流向脉动动能的传递。 故三角形沟槽表面附近横向脉动最大,对流向脉动
8、动能的传递功能最强,流体对壁面的作用力最大,壁面摩擦 阻力最大,减阻效果最差;刀刃形沟槽表面横向脉动最小并且流向脉动动能的 传递功能最弱,流体对壁面的作用力最小,壁面总阻力最小,减阻效果最好。3. 4湍流统计分析图7为光滑表面与三种沟槽涡量曲线图。由计算结果可知光滑表面流向涡 量、展向涡量和垂向涡量平均值分别为1. 5 X 1051/ s、3. 8 X 1061/ s和5 X1051/ s。三角形沟槽表面流向涡量、展向涡量和垂向涡量的最大值分别为 1. 25X 105 1/ s 、2. 5 X 1061/ s 和2. 0 X 106 1/ s; 扇贝形沟槽三个方向涡量的最大值分别为 1.05
9、X 1051/ s、2. 0 X 106 1/ s 和 1.5 X 106 1/ s, 小于三角形沟槽表面上相应涡量值;刀刃形沟槽三个方向涡量的最大值分别为1.0 X 1051/ s、1.8 X 106 1/ s 和2. 4 X 106 1/ s, 在三种情况下最小。由此可见三种沟槽表面流向涡量和展向涡量的最大值小于光滑表面附近相应涡量 的最大值,垂向涡量最大值大于光滑表面,即沟槽抑制流向涡的展向运动。由于 三角形沟槽流向涡量和展向涡量的最大值在三种沟槽中最大 ,与光滑表面相应 涡量值的差距最小,流向和展向涡量被三角形沟槽表面减小的程度最小 ,故其对流向涡的展向运动的抑制作用最不明显,对湍流拟
10、序结构的影响程度最差,减 阻率最低13。刀刃形沟槽流向涡量和展向涡量的最大值在三种沟槽中最小 ,与光滑表面相应涡量值差距最大,流向和展向涡量值被刀刃形沟槽表面减小的程 度最大,故其对流向涡的展向运动的抑制作用及对湍流拟序结构的影响程度最大减阻率最高。-非光滑衣面 光滑农面|5.0*102.010L0xl0 l非光滑表面 光滑哀面t 5.0x101/ 、5.0101常-1心1叶2. OxlO-f-2后 i(r :2.0x10-,l.SxltfJjLOxlO非光滑表面-m 5.0x10 負 0-5.0x107mif-L 55;-2. 0x10+p. 2 3 4F/mm一一1.10-5,OxlO3-
11、光滑农面2 K/mnd扇贝形展向涡量團-喟咚怎L 2xi012八6.0x10 4.0x10s2. 0x10非光滑衣面“2.0皿-tOxlCT-8. 0x10- 0.OMltf- .DxlO5*- -3,0x10- -九(MW/光滑表丙-5. ()X10 LQxlO -1,510-2. Bl护 r Q6. OxlQq5 0xl0$42 0刈叫3.04非光滑表面2. (M叫/1. OxlQJ1理炖非光滑i0-r I I I rT1 =5 O* I.22.0皿l.OxlO55.0x10秋滑表面光滑表面123K/nuig闻贝形垂向涡絃图-1y/nnf三角形展向涡量圏L 26x10 L 0010T 7.60x10 _ a. ooxio 耳 2,50xltf WBE-2p 50* 10* 怎临0010b -7. 50畑 -I. 0010 -L2&1071 2 3y/mo刀刃形垂向涡量图尹滑表血光滑表面、5.01054.0xltf声兀滑表面3. OxlO12 Oxl(fI.OIO1U-t QmI-2. Oxlff玄 Oxg11 * * 、 | | i *0 12 3 4Y /ituui三角形展向涡量图图7沟槽三个方向涡量图
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1