FLUENT软件在圆柱绕流模拟中的应用

使用计算流体力学软件FLUENT，模拟均匀来流绕固定圆柱的流动，模拟雷诺数为20，40，100时的绕流流动，得到流场的流函数等值线图和速度矢量图。计算结果表明：当雷诺数增加时，流动表现出一系列不同的构造。在雷诺数约为40前后流场有明显变化。小于这个数时，存在一对位置固定的旋涡。大于40时，流场开始变得不稳定，旋涡扩大、脱落、又生成，逐渐发展成两排周期性摆动和交错的旋涡。并与实验及数值模拟结果比较，确认FLUENT能够很好地预测流动结构。     

并列旋转双圆柱流动特性的数值模拟

该文基于k - ε湍流模型,采用Galerkin有限元法对并列旋转双圆柱的绕流特性进行了数值模拟,计算的雷诺数为 1550.为了考查两圆柱旋转和间距的相互作用,文中采用三种间距比分别是T/D = 1.2,1.6和3.0 (T 为两圆柱中心之间的距离,D为圆柱直径) 和一系列不同的旋转速度比 (|α| ≤ 2).计算显示,当 |α| = 0,即圆柱不转动时,对应三种间距有三种典型的流型,单钝体流型对应小间距、偏流对应中等间距和对称流对应大间距;当 |α| 达到临界值时,涡脱落得到了有效的抑制,流动趋于稳定,升力系数和阻力系数的脉动值趋于零;平均升力系数和阻力系数随着 |α| 的增大分别增大和减小.     

底部加热的两板间水平流动进口段流动特性分析

采用二维流体力学基本方程组对普朗特数Pr=0.72的具有水平流动的底部加热两板之间的流动特性进行数值模拟。结果表明,对于给定相对瑞利数r,随着雷诺数Re的增加,底部加热的两板之间依次出现定常对流滚动,均匀行波对流滚动,具有水平流动进口段的局部行波对流滚动及水平流动。对于相对瑞利数r=3,在雷诺数9.0＜Re＜21.0范围内,系统是局部行波对流滚动,存在水平流动进口段。随着雷诺数增加,水平流动进口段长度增加。对于雷诺数Re=25.0,在相对瑞利数4＜r＜12范围内,两板之间是局部行波对流滚动,出现水平流动进口段。随着相对瑞利数增加,水平流动进口段长度减小。最后给出了进口段长度随雷诺数或相对瑞利数变化的经验式。     

流动浆液中宾汉剪切力变化规律^＊

河道中洪水高含沙水流和工业管道输送的高浓度浆体,多数属于宾汉型流体.宾汉剪切力大小对于高含沙水流或浆体的流动特性有很大影响.本文通过水槽试验,研究了流动浆液中宾汉剪切力随着水流的有效雷诺数的增加而减少的变化规律.在层流向紊流过渡时(有效雷诺数按2100计)宾汉剪切力约为静止条件下的60％.当有效雷诺数达到3.5×10~4时,宾汉剪切力接近为零,浆体可近似地按牛顿流体考虑.     

多圆柱绕流旋涡脱落和流场形态概论

当流体绕多个圆柱流动时,圆柱的个数、相对位置和大小都明显影响着绕流流场结构。基于已有研究成果,归纳总结了流场中放置双圆柱、三圆柱、四圆柱和流场中主圆柱附近放置1个或2个附属小圆柱时,随圆柱相对位置变化,流场结构和旋涡脱落状态的演变规律。双圆柱相对位置分串列、并列和错置3种,其绕流流态按间距比和水流攻角对旋涡脱落形态的影响划分,流态在临界间距比上下变化很大;三圆柱呈“品”字布置,当间距比小于临界值时,出现明显的偏流现象,大于临界值后不再出现偏流现象;三圆柱呈倒“品”字布置,在任何间距比下,流动状态相对于来流几乎都是对称分布的;四圆柱绕正方形排列,实验研究发现上游两圆柱后没有旋涡形成和脱落,下游两圆柱旋涡脱落反相同步,而数值模拟得到上下游柱后都有明显的涡脱落;主圆柱尾流中存在附属小圆柱时,调整小圆柱的位置能够有效抑制尾流中的旋涡脱落。     

矩形管道中平板障碍物后旋涡发射特性研究

通过理论推导及流动显示实验，研究了矩形管道中平板障碍物后旋涡发射频率及尾流特征。结果得到了旋涡发射频率与雷诺数、平板阻塞率之间的关系，特别是发现了在高阻塞率条件下，当雷诺数增至某个临界值时旋涡发射被抑制，以及尾流存在周期性旋涡发射和无明显周期性的回流区两种不同的流动状态。本文还分析了上述现象的流动机理，即尾流特性主要取决于两边的分离剪切流之间的相互作用。     

电加热圆柱尾迹涡的流动显示

本文在低湍流度风洞及水槽中，就低雷诺数范围内对加势圆柱尾迹涡进行流动显示研究，实验表明在水介质中均匀加热圆柱时可促进其尾涡的脱落、而沿展向局部加热时，其局部涡脱落频率增大，导致原准二维卡门涡街三维化；相反空气中加热圆柱对旋涡脱落有扼制作用，本文首次观察到加热诱发的圆柱尾迹不稳定性的展向波，表现为Ｒｅ＾＊〈Ｒｅ（临界涡脱落雷诺数）时双子涡展向断裂成间距为展向波长λ的涡段。圆柱加热增强了尾迹近区的不稳     

串列梯形柱绕流水动力特性数值分析

该文采用数值模拟分析了雷诺数为150和间距比为3的串列双梯形柱绕流水动力特性。当梯形柱背流面与迎流面长度比为d/D=0、0.3和0.5时,双柱的上侧旋涡融合,在尾部形成了不对称的P+S泄涡模式,导致下游柱体承受的平均升力偏离零值。d/D=0.7的双梯形柱泄涡转变为2S模式,表现为上游柱体旋涡与下游柱体旋涡融合后交替脱落。双方柱(d/D=1)间隙以及下游方柱尾部均形成了准静态涡,剪切层有小幅摆动,但未卷曲形成旋涡。d/D=0和0.3时,上游柱体后方由于边界层再附着形成了次生涡,增大了局部压力。随着d/D的增大,上游柱体旋涡形成长度增长,双柱间隙流动趋于稳定,下游柱体前缘的横向速度脉动逐渐减弱,其后方的尾迹也趋于对称。上游柱体的遮蔽效应使其平均阻力大于下游柱体,而间隙流的不稳定性与边界层的交替再附着使得下游柱体升力偏大。双方柱绕流场处于拟稳定状态,其水动力系数远小于三棱柱和梯形柱。相较于双圆柱,三棱柱和梯形柱的边界层分离点靠前,增大了绕流的升阻力。     

通气空泡重力效应研究

该文基于VOF方法和有限体积法,求解气水多相流动的RANS方程及标准湍流模式,研究了通气空泡的重力效应.结果表明模拟得到的绕圆盘通气空泡的重力效应与Buyvol的分类准则吻合,通过数值模拟复演了通气空泡的双涡管泄气模式;证实由于与尾空泡的相互作用,重力作用下绕航行体通气空泡在尾部附近出现下沉,航行体上的空泡截面呈现出"月牙"形;阻力随空化数而减小,并在空泡完全包裹住航行体前部锥段时出现突降.     

绕轴对称体三维非定常空化流动的数值与实验研究

采用数值和实验相结合的方法研究了绕轴对称体三维非定常空化流动现象。实验中,采用高速摄像的技术和动态测力系统对绕半球型和平头型轴对称体的非定常空化流场及其动力特性进行了分析,观察了在不同空化数下,绕半球型和平头型轴对称体的空泡形态;测量了轴对称体受到的阻力,并对阻力信号进行了时频分析,得到了在非定常空化阶段,轴对称体动力特征频率。基于实验现象,发展了一种基于密度修正的分域空化流动的计算方法,并与实验结果进行了对比,结果表明:该方法可以较好的模拟绕轴对称体三维非定常空化流动现象。绕不同头型轴对称体的空化流场结构均存在明显的三维非定常流动特性,并且,空化流场形态与动力特征频率存在高度的相关性,不同头型轴对称体的非定常空化流场结构存在较大的差异,半球型轴对称体空泡流动的脉动主要是空泡尾部的高频小脱落引起的,而平头型轴对称体的空泡流脉动成分主要是大尺度的漩涡空泡团的周期性脱落。     

不同空泡脱落模式下空化多尺度特性研究

云状空化是一种以相变为核心,包含多种尺度空泡生长演变行为的强瞬态多相流动现象。该文采用双向耦合欧拉-拉格朗日多尺度空化模型,模拟了空化数σ=1.2时,绕8°攻角三维NACA66水翼空化流动,通过VOF(volume of fluid)界面捕捉法求解欧拉体系下大尺度空穴演化,利用拉格朗日体系下的离散空泡模型(DBM)追踪亚网格尺度离散空泡运动及生长溃灭过程。结果表明：绕水翼云状空化流动,包含附着型片状空泡生长与发展、回射流产生与推进、云状空泡脱落和宏观片空泡再次生长以及压力波推进4个阶段。离散空泡数量及尺度随大尺度空穴周期性演变而变化,主要分布在强湍流脉动区。在各空化发展阶段,离散空泡集中分布在水翼中后部,微空泡直径概率密度函数均符合单峰Gamma分布,最高空泡数密度泡径为50μm。在湍动程度较小的附着型空穴生长阶段,微空泡数密度相对于空泡直径呈-5/3幂率特征;在湍动程度较大的回射流发展阶段、云状空穴脱落溃灭及压力波推进阶段,较小尺度的微空泡数密度相对于空泡直径呈-2/3幂率特征,较大尺度的微空泡数密度相对于空泡直径呈-6幂率特征。     

二维叶栅通道直接数值模拟涡系分析

该文基于有限体积法对NACA0065叶栅内二维不可压缩流动进行直接数值模拟,得到不同雷诺数下叶栅通道内的瞬时涡系结构和时均物理量分布,进而表明了涡系结构对叶片气动系数的影响。在低雷诺数情况下,尾迹旋涡的周期性脱落频率导致升阻系数的周期性变化。利用傅立叶变换,气动系数谱中的主频与涡脱落频率精准对应,进一步分析表明次频是由脱落周期内相邻区域的不同涡核低压区交替出现而导致。在高雷诺数下,气动力系数频谱中出现了更多的高频谐波。此外,该文利用Liutex理论中的R-S分解和涡量输运方程解释了NACA0065叶片表面剪切层产生旋涡的机理,并研究了强剪切区和弱剪切区处的涡核特性。该研究结果可使工程界对NACA0065叶栅中的相应旋涡结构及其相互作用有了更深入的物理理解。     

方柱绕流数值模拟

利用FLUENT中的RNG k-ε算法对长宽比H/B=1.5方柱绕流进行数值模拟,分别模拟了在雷诺数Re=800的亚临界条件下,4种不同情况下的流动。分别为:流体流过单个方柱,流体流过间距为1H,2H,3H的两并列方柱的流动。根据数值模拟结果,分别得出了4种流动状态下流场,壁面的升阻力系数。并对升力系数进行傅里叶变换,分析及归纳了方柱绕流旋涡脱落的规律。     

不同边壁条件下矩形明渠水力特性实验研究

为了研究明渠水流水力特性变化规律,在室内实验的基础上,采用矩形明渠变坡水槽,通过改变边壁条件、坡度和流量进行实验,得到了较好的实验数据。实验结果表明:边壁条件和坡度相同的情况下,随着流量的增加,水深呈线性增加,平均流速和雷诺数呈现幂函数增大;当Fr1时,弗劳德数随着流量的增加而增加;当1Fr1.526时,弗劳德数随着流量的增加先增加而后减小;当Fr1.526时,弗劳德数随着流量的增加而减小;曼宁糙率系数随着流量的增加呈对数减小趋势。边壁条件和流量相同的情况下,随着坡度的增大,水深呈减小趋势,平均流速、雷诺数、弗劳德数和曼宁糙率系数均呈现增加的趋势。并由实验数据得出了雷诺数与坡度和流量关系的经验公式。     

混合对流绕圆柱涡旋脱落及对传热影响的数值研究

利用本文作者改进的关于流函数一涡量方向的数值求解方法，模拟了混合对流绕圆柱涡旋脱落及对传热的影响，得到了雷诺数Ｒｅ与格拉肖夫数Ｇｒ双重作用下圆柱尾迹涡型转变的规律，并数值研究耻存在小扰动的自由来流、圆柱旋转（γ）及浮力与来流成不同角度（δ）等因素对旋涡脱落，对传热的影响。除Ｒｅ数处，δ，Ｇｒ数对尾涡流型及传热产生较大影响，可利用这些参数来控制绕流流动及强化传热作用。     

C形钝体绕流的尾涡结构与水动力特性分析

该文以不同攻角(0°<α<180°)的C形钝体在低雷诺数(40相似文献   

动网格技术在非稳态空化流计算中的应用

该文基于对钝锥头航行体的空化流动的实验观测,利用商业软件FLUENT[1]考核了分块网格模型对降低动网格方法误差影响、提高空化流计算精度的可行性和有效性。通过与实验的对比表明,采用分块网格模型能够使基于动网格方法的非稳态空化流计算很好的捕捉到空泡生长,回射,翻卷,脱落,溃灭这一过程,并获得较为准确的航行体的运动轨迹和所受的水动力。进而分析了流场涡结构与空泡演化的机理。     

绕二维对称水翼的通气空泡流数值研究

本文通过有限体积离散方法求解Reynolds平均Navler—Stokes方程，采用VOF法追踪气一液两相交界面。研究了在高压低速来流，即本质上不发生自然空化的条件下，绕二维对称水翼的非定常通气空泡流动。分析了物理参数——通气方向角和攻角，对通气空泡外形、稳定性及水动力参数的影响。计算结果表明：在通气空泡的形成过程中，通气方向角存在着一个临界值，使得空泡形态由片状空泡过渡到超空泡；特定的参数组合将会诱发空泡表面出现波动现象。而且，当攻角做周期性变化时，水动力参数的变化也随之满足周期性规律，且出现阶段性振荡。     

水平进水口前临界淹没水深试验研究

旋涡是水工建筑物进水口前常见的水力现象,会给水工建筑物和水电设备带来严重的危害。消除旋涡最广泛采用的方法是保证进水口前的水深超过临界淹没水深。从实际工程应用计算考虑,利用水力模型试验,重点分析了弗劳德数Fr、雷诺数Re对临界淹没水深的影响,并由试验结果得到旋涡相对临界淹没水深H/D随弗劳德数Fr和雷诺数Re的变化规律。结果表明:临界淹没水深随着弗劳德数和雷诺数的增大而增大,并且与两者均有幂函数关系。综合弗劳德数Fr和雷诺数Re对H/D的影响规律,得到了相对临界淹没水深与弗劳德数Fr和雷诺数Re的定量关系式,关系式拟合度为0.863,可以满足工程实际计算要求。     

喷气推进航行体空泡流的尺度效应研究

基于Reynolds平均N-S方程和Mixture模型,耦合自然空化模型,运用气、汽、液多相空泡流动的数值模拟方法,模拟了水下航行体自然空泡和尾部燃气射流发展的非定常过程。对于六种模型尺度,计算分析了自然空泡和燃气射流的尺度效应。计算结果表明:随着模型尺度减小,相应的雷诺数也缩小,水流紊动减弱,空泡的脱体点后移,航行体头部低压区范围减小,导致压差阻力系数增大。同时,雷诺数越小意味着黏性力影响越显著,相应的黏性阻力系数增大。模型缩尺对于尾部燃气射流的发展影响较小,不同尺度模型的马赫数分布与阻力系数分布基本相同。