圆柱绕流的数值模拟
低雷诺数下圆柱绕流减阻的数值模拟研究
硕士学位论文低雷诺数下圆柱绕流减阻的数值模拟研究A NUMERICAL STUDY OF FLUID FORCES REDUCTION ON A CIRCULAR CYLINDER AT LOW REYNOLDS NUMBER黄玮哈尔滨工业大学2005 年 3 月国内图书分类号：O357.1国际图书分类号：621工学硕士学位论文低雷诺数下圆柱绕流减阻的数值模拟研究硕士研究生：黄玮导师：周超英副教授申请学位：工学硕士学科、专业：机械电子工程所在单位：深圳研究生院答辩日期： 2005 年 3 月授予学位单位：哈尔滨工业大学Classified Index: O357.1U.D.C.: 621Dissertation for the Master Degree of EngineeringA NUMERICAL STUDY OF FLUID FORCES REDUCTION ON A CIRCULAR CYLINDER AT LOW REYNOLDSNUMBERCandidate:Supervisor:Academic Degree Applied for: Specialty:Affiliation:Date of Defence: Huang WeiAssociate Prof. Zhou Chaoying Master of Engineering Mechatronics Engineering Shenzhen Graduate School March, 2005Degree-Conferring-Institution: Harbin Institute of Technology哈尔滨工业大学工学硕士学位论文摘要研究圆柱绕流的减阻控制对于减小流动诱发振动、提高结构体抗疲劳强度有着非常重要的作用。
数值模拟能够克服实验的不足，高效的给出相关流场的具体信息。
正是因为数值模拟具有较多的优点，所以计算机数值模拟已逐渐成为流体力学强有力的研究工具。
不同雷诺数下倾斜圆柱绕流三维数值模拟研究
不同雷诺数下倾斜圆柱绕流三维数值模拟研究近年来，研究倾斜圆柱绕流特性引起了学界的广泛关注。
圆柱绕流可分为水平和垂直两类，其中倾斜圆柱绕流为一种特殊的二维绕流状态，它在一定雷诺数范围内具有更复杂的流场结构特性，并且受水文物理过程的影响更为显著，研究其特性更为重要。
本研究使用时间和空间设置，以带边界流作为边界条件，利用基于六边形网格的数值模拟方法研究不同雷诺数下的倾斜圆柱绕流特性。
实验参数包括：倾斜角度α=20°，Re=1000 ~ 10000，向心轴比例范围为0.5 ~ 2.0，圆柱入口处外提升速度Um=0.3 ~ 0.8，及空气密度ρ=1。
有鉴于此，本研究根据不同雷诺数和向心轴比例，计算出倾斜圆柱绕流特性。
首先，主要考察不同雷诺数Re下倾斜圆柱绕流的流态特性，包括在不同位置的压力梯度，流场动量，温度梯度，流态结构以及涡度等信息。
其次，重点考察不同向心轴比例和轴向外提升速度下倾斜圆柱绕流的流态特性，包括压降，动量，温度梯度，以及不同方向的涡度分布。
结果表明，不同的雷诺数和向心轴比例会对倾斜圆柱绕流的流动特性产生明显不同的影响。
随着雷诺数的增大，压力梯度增大，动量梯度减小，温度梯度增大，涡度明显减少，圆柱内部的流场会变得更加复杂，气泡变小，而且其会从一种混合流场演变为一种逆流的流场结构。
另外，随着向心轴比例的增加，轴向外提升速度的变化会出现显著影响，但随着向心轴比例的增加，压力梯度会逐渐减小，动量梯度增大，温度梯度变化不大，涡度分布也会有较大变化。
研究结果表明，在不同雷诺数和向心轴比例范围内，倾斜圆柱绕流的流动特性会发生明显的变化。
本研究对于进一步理解流动特性和确定流动行为有重要的理论意义，同时也为实际工程的设计提供了参考。
总的来说，本研究通过应用数值模拟方法研究不同雷诺数下倾斜圆柱绕流特性，得出上述结论。
未来可以将此模拟实验方法应用于建立更复杂物理系统的研究，以更深入地理解绕流特性和其流动性质。
基于H-O型网格的圆柱绕流的数值模拟
基于H-O型网格的圆柱绕流的数值模拟摘要：本文用GAMBIT针对相同的模型划分了不同结构的网格，用FLUENT进行了圆柱绕流的数值模拟，并将模拟结果跟文献的实验数值做了对比，数值计算的结果与相关实验结果具有合理的一致性。
同时也证明了基于H-O型网格的圆柱绕流的模拟精度要高于普通的四边形网格。
关键词：H-O型网格圆柱绕流数值模拟Numerical Simulation of Fluid Flow around Circular-cylinder Based on theH-O GridAbstract: This paper supplies different grids for the same model based on the GAMBIT. Numerical simulation of the fluid flow around circular-cylinder is executed by the software of FLUENT. Compared the result with the tested data, the result of the numerical simulation is coincidental with that of experimental. The numerical simulation results also approve that the result based on the H-O grid is more exacted than the traditionally quadrangular grid.Key words: H-O Grid；Fluid Flow around Circular-cylinder；Numerical Simulation1 引言静止地圆柱绕流早期的研究由于受到理论和实验条件的限制,圆柱绕流主要侧重于对圆柱体受力的实验测定,随着计算流体力学和现代实验技术的迅速发展，数值模拟及实验测定方法被广泛地应用于这些方面。
流体力学Fluent报告——圆柱绕流
亚临界雷诺数下串列单圆柱与圆柱绕流的数值模拟之阳早格格创做目要：原文使用Fluent硬件中的RNG k-ε模型对付亚临界雷诺数下二维串列圆柱战圆柱绕流问题举止了数值钻研，通过截止对付比，分解了雷诺数、柱体形状对付柱体绕流阻力、降力以及涡脱频次的效率.普遍而止，Re数越大，圆柱的阻力越大，圆柱体则可则；而Re越大，二种柱体的降力均越大.相对付于圆柱，共种条件下，圆柱受到的阻力要大；好异天，圆柱涡脱降频次要小.Re越大，串列柱体的Sr数越靠近于单圆柱体的Sr数.闭键字：圆柱绕流、降力系数、阻力系数、斯特劳哈我数正在工程试验中，如航空、航天、航海、体育疏通、风工程及大天接通等广大的本量范畴中，绕流钻研正在工程本量中具备要害的意思.当流体流过圆柱时, 由于漩涡脱降，正在圆柱体上爆收接变效率力.那种效率力引起柱体的振荡及资料的疲倦，益坏结构，成果宽沉.果此，近些年去，稠稀博家战教者对付于圆柱绕流问题举止过细致的钻研，特天是圆柱所受阻力、降力战涡脱降以及涡致振荡问题.沈坐龙等[1]鉴于RNG k⁃ε模型，采与有限体积法钻研了亚临界雷诺数下二维圆柱战圆柱绕流数值模拟，得到了圆柱战圆柱绕流阻力系数Cd与Strouhal 数随雷诺数的变更顺序.姚熊明等[2]采与估计流体硬件CFX中LES模型估计了二维不可压缩匀称流中孤坐圆柱及串列单圆柱的火能源个性.使用非结构化网格六里体单元战有限体积法对付二维N- S圆程举止供解.他们着沉钻研了下雷诺数时串列单圆柱正在分歧间距比时的压力分集、阻力、降力及Sr数随Re数的变更趋势.费宝玲等[3]用FLUENT硬件对付串列圆柱绕流举止了二维模拟，他们采用间距比L/D(L为二圆柱核心间的距离，D为圆柱直径)2、3、4共3个间距举止了数值分解.估计均正在Re = 200 的非定常条件下举止.估计了圆柱的降阻力系数、尾涡脱降频次等形貌绕流问题的主要参量，分解了分歧间距对付圆柱间相互效率战尾流个性的效率.圆柱绕流的一个要害个性是震动形态与决于雷诺数.Lienhard[4]归纳了洪量的真验钻研截止并给出了圆柱体尾流形态随雷诺数变更的顺序.当Re<5时，圆柱上下游的流线呈对付称分集,流体本去不摆脱圆柱体，不旋涡爆收.此时与理念流体相似，若改变流背，上下游流形仍相共.当5<Re<40时,鸿沟层爆收分散,分散剪切层正在圆柱体里前产死一对付宁静的“附着涡”.当40<Re<150时，震动脆持层流状态而且流体旋涡接替天从圆柱后部做周期性的脱降并正在尾流中产死二列接叉排列的涡，即卡门涡街.从150<Re<300启初，旋涡里里启初由层流背湍流转捩，直至减少至3x105安排，此时圆柱体表面附近的鸿沟层仍为层流，所有涡街渐渐转化成湍流，及e<3xl05称为亚临界天区.当3xl05<Re<3.5x106时,鸿沟层的震动也渐渐趋于湍流状态，尾流中不明隐的涡街结构，称为临界状态.[5]圆柱绕流的另一个隐著个性是斯特劳哈我数是雷诺数的函数.早正在1878年，捷克科教家Strouhal[6]便对付风吹过金属丝时收出鸣喊声做过钻研，创造金属丝的风鸣音调与风速成正比,共时与弦线之细细成反比，并提出估计涡脱降频次f的体味公式:式中即斯特劳哈我数Sr由Re所唯一决定.原文使用Fluent硬件中的RNG k-ε模型对付亚临界雷诺数下二维串列圆柱战圆柱绕流问题举止了数值钻研，通过截止对付比，分解了雷诺数、柱体形状对付柱体绕流阻力、降力以及涡脱频次的效率.1.数教模型1.1统造圆程对付于停止圆柱绕流，原文钻研对付象为二维不可压缩震动.正在直角坐标系下，其疏通顺序可用N-S圆程去形貌，连绝性圆程战动量圆程分别为:其中ui为速度分量；p为压力；ρ为流体的稀度；ν为流体的能源黏性系数.对付于湍流情况，原文采与RNG k⁃ε模型，RNG k⁃ε模型是k⁃ε模型的矫正规划.通过正在大尺度疏通战建正后的粘度项体现小尺度的效率，而使那些小尺度疏通有系统天从统造圆程中去除.所得到的k圆程战ε圆程，与尺度k⁃ε模型非常相似，其表白式如下：其中Gk为由于仄衡速度梯度引起的湍动能的爆收项，，，体味常数=0.084 5，==1.39，=1.68.相对付于尺度k⁃ε模型，RNG k⁃ε模型通过建正湍动粘度，思量了仄衡震动中的转动及转动震动情况，RNG k⁃ε模型不妨更佳的处理下应变率及流线蜿蜒程度较大的震动.1.2相闭参数圆柱绕流的相闭参数主要有雷诺数Re、斯特劳哈我数Sr、降力系数Cl战阻力系数Cd，底下给出各个参数的估计公式战物理意思.雷诺数Re与圆柱绕流的状态战雷诺数有很大闭系，雷诺数代表惯性力战粘性力之比：其中U为去流速度；L为个性少度，原文与圆柱直径或者圆柱边少；为流体稀度；、分别为流体介量能源粘度战疏通粘度.斯特劳哈我数Sr是Strouhal 指出圆柱绕流后正在圆柱后里不妨出现接替脱降的旋涡，旋涡脱降频次、风速、圆柱直径之间存留一个闭系：式中：Sr为斯托罗哈数，与决于结构的形状断里；f 为旋涡脱降频次；L为结构的个性尺寸； U 为去流速度.阻力系数战降力系数是表征柱体阻力、降力的无量目参数.定义为：，式中ρ为流体稀度；V为去流速度；A为迎流截里里积；战.由于涡脱降的闭系，阻力系数将爆收振荡，原文采用仄衡脉动降力去钻研，即与圆均根值去钻研.2.数值估计2.1物理模型二维数值模拟单圆柱流场估计天区的采用如图1所示，圆柱绕流以圆柱体直径为个性尺度D，采用圆柱半径为1.5 mm，估计天区为9D×32D的矩形天区.柱1距上游少度图 1 串列圆柱战圆柱的估计天区5D，下游少度27D，脆持二柱间距 L/D= 2. 5D稳定 (L是二圆柱核心连线少度)，二柱到上下鸿沟距离相等.对付于圆柱绕流，采用圆柱边少为个性少度，D=30mm.2.2网格区分估计天区采与分块结构化网格，柱体表面网格干加稀处理，鸿沟区网格相对付稠稀.简直网格区分情况睹图2.其中串列圆柱网格31116个节面，30615个四边形里单元；串图 2 圆柱绕流与圆柱绕流估计域的网格区分列圆柱46446个节面，46550个四边形里单元.2.3鸿沟条件管讲壁里战柱体表面均采与无滑移的停止壁里条件.而出心采用速度出心，出心采用自由出流.去溜速度大小根据Re去树坐，雷诺数分300、3000、12000、30000四个等第，速度大小依次为0.1m/s、1m/s、4m/s、10m/s.2.4估计模型原文湍流模型采与尺度壁里函数的RNG k-ε模型.采与有限容积法供解二维不可压缩粘性流体非定常震动统造圆程，即把估计天区分成很圆柱近壁里网格多小的统造体，对付每个统造体的各个变量举止积分.统造圆程的对付流项采与二阶迎风圆法失集，速度战压力采与SIMPLE算法耦合供解，将所有天区瞅成一个完全举止耦合估计.动量、湍动能战湍动耗集率均采与二阶迎风圆法.先定常估计流场，再用定常估计的截止动做非定常迭代的初初值举止估计.根据初略估计的涡脱频次，牢固树坐时间步少为0. 002s, 正在每个时间步内树坐迭代次数为20.流体介量为液态火.3.估计截止3.1网格模型考证为考证网格独力性，原文估计了网格节面数为8346，里单元为8932的细网格、节面数为31116，里单元为30615的稀网格、节面数为63432，里单元为67434的细稀网格下Re=200、L/D=2的串列网格的Sr数，截止隐现三套网格的估计截止分别为0.143、0.133、0.133.故稀网格可用.而圆柱绕流则采与共级别网格.[7]的估计数据相比较，比较图像如图3所示，最大缺面为2.2%.图3串列圆柱分歧间距的Sr数估计对付比3.2流线与涡量图图 6 Re=3000圆柱绕流流线图图 7 Re=3000圆柱绕流涡量等值线图图 4 Re=3000圆柱绕流流线图图 5 Re=3000圆柱绕流涡量等值线图原文给出了估计历程中雷诺数Re=3000，t=1s时的流线图战涡量图.3.3阻力系数图 9 Re=3000圆柱绕流脉动阻力系数图 8 Re=3000圆柱绕流脉动阻力系数原文给出了Re=3000时，圆柱绕流战圆柱绕流的脉动阻力系数图如下.由图9战错误!未找到引用源。
圆柱绕流数值模拟资料
圆柱绕流的数值模拟研究摘要：选取直径为D=10mm的圆柱及6D×3D的计算区域，利用GAMBIT进行模型的创建模型，对计算区域采用分块网格划分与结构化网格划分相结合的技术进行网格划分。
对0.03m/s~1.0m/s的低流速情况下的圆柱绕流进行模拟研究，结果发现在速度达到0.1m/s前圆柱后侧没有出现明显的漩涡，在速度大于0.1m/s后漩涡开始出现，当速度达到0.5m/s时漩涡的范围最大。
最后利用FLUENT的网格自适应技术对入口速度为0.5m/s的情况进行了网格加密，发现网格自动加密可以改进网格分布情况，但对计算结果的影响程度有限。
关键词：网格划分；圆柱绕流；涡量；网格自适应钝体绕流中尤其以圆柱体的绕流问题最为经典和引起人们的注意。
圆柱绕流属于非定常分离流动问题，在工业工程中的应用非常广泛。
圆柱绕流同时也是一个经典的流体力学问题，流体绕圆柱体流动时，过流断面收缩，流速沿程增加，压强沿程减小，由于黏性力的存在，就会在柱体周围形成附面层的分离，形成圆柱绕流。
而由于圆柱的存在，会在圆柱迎水面产生壅水现象，同时也增加了圆柱的受力，使得圆柱绕流问题变得十分复杂。
研究圆柱绕流问题在工程实际中也具有很重要的意义。
如在水流对桥梁、海洋钻井平台支柱、海底输运管线、桩基码头等的作用中，风对塔建筑、化工塔设备、高空电缆等的作用中，都有重要的工程应用背景。
因此，对圆柱绕流进行深入研究，了解其流动机理和水动力学规律，不仅具有理论意义，还具有明显的社会经济效益。
1数学模型与计算方法1.1几何模型结合本文研究目标，取圆柱直径D=10mm，计算区域为6D×3D的矩形区域，如图1所示。
上游尺寸1.5D，下游尺寸4.5D。
使用GAMBIT建模软件按照图1所示的计算域建立了二维的计算模型。
图1计算区域1.2网格划分及边界条件设置为提高模拟精度，计算区域采用分块网格划分与结构化网格划分相结合的技术。
计算区域共分两块，尺寸见图1所示。
圆柱绕流的三维数值模拟
式中Ｒ．ｆ０分别为作用于圆柱上总的升力和阻力，沿柱长方向各断面的升力和阻力系数分
别定义为ｃ￡＝２凡／州２Ｄ，ＣＤ＝２％／ｐｕ ２Ｄ，式中兄，ｆｊ分别为作用于圆柱断面单位
宽度周线上的升力和阻力。为了与三维数值结果比较，本文对圆柱绕流进行了二维数值模拟。 图５（ａ）为两种数值模拟的升力系数和阻力系数曲线，图５（ｂ）为沿柱长平均的升力系数与各断 面的升力系数曲线，图５（ｃ）为沿柱长平均的阻力系数与各断面的阻力系数曲线．给出１／２柱 长三个断面的图形。从图中可咀看出．雷诺数在亚临界区范围内，圆柱升力系数和阻力系数
圆柱周围的流场利用ＣＦＸ－４计算流体力学软件求得，计算采用有限体积法，为解决动量 方程中速度分量和压力耦台问题．采用ＳＩＭＰＬＥ算法。数值计算中，由于交错网格受到限制， 一方面是有限体积法中各变量占用不同的控制体，每个控制体可能会有不同的几何形状，使 问题复杂化：另一方面，网格信息的存储量特别太，因此，计算中采用非交错网格，并将压 力和速度定义在网格的中心上，这样能显著减少对结点变量的存储。流场的数值计算采用贴 体网格，它把具有复杂边界的物理区域变换到规则的计算区域，从而在数值计算中能够比较 准确地满足边界条件，提高解的精度。数值求解Ｎ．ｓ方程时，为克服数值模拟在高雷诺数时 的数值不稳定性，计算中采用了ＱＵＩＣＫ迎风格式（Ｌｅｏｎａｒｄ １９８１［７１，Ｆｅｆｉｔａｓ ｅ￡ａ１．１９８５）其对 流项为三阶精度，其余项如扩散项等为二阶精度。
ｇｏｌｖｎＩｎ：Ｐｍ刚ｉｎｐ ｏｆＴｈｅＴｅｎｔｈｌｎｌｅｒｒＩａｔｉｏｎａｌＯｆｆＢｈｏｒｅＡｎｄＰｏｌａｒＥｎｇ［ｎ＝ｅｔｉｎ８Ｃｏｎｋｒ∞ｃｅ．Ｓｅａｔｔｌｅ：ＩＳＯＰＥ Ｐｒｅｓｓ，２０００（ＩＨ）＝４７１４７７
三维波浪作用下大直径圆柱绕流的数值模拟
三维波浪作用下大直径圆柱绕流的数值模拟祝兵;宋随弟;谭长建【摘要】In order to research interaction between three-dimensional waves and structures, a numerical model for the interaction of three-dimensional waves and a structure was established based on the concept of two-phase flow, the volume of fluid (VOF) method and incompressible viscous fluid equations in large eddy simulation ( LES). Interaction between three-dimensional waves and a large-diameter circular cylinder was simulated numerically. The virtual boundary force method (VBFM) and two-step boundary-pinpoint method were adopted to determine the boundaries of the structure. The research show that the numerical simulation result of the circumfluence factor for a large-diameter circular cylinder has a good agreement with the analytical solution. The established numerical model can successfully simulate the interaction between three-dimensional waves and structures.%为探讨三维波浪与结构物的相互作用,以两相流概念、大涡模拟的不可压缩粘性流体运动方程和自由水面追踪分段线性近似的流体体积(VOF)法为基础,建立了三维波浪与结构物相互作用的数学模型；对三维波浪作用下大直径圆柱绕流进行了数值模拟,用两步边界定位法和虚拟边界力法确定波浪与结构物接触面.结果表明:大直径圆柱绕流系数的数值计算结果与理论解吻合,所建立的数学模型能够很好地模拟三维波浪与结构物的相互作用.【期刊名称】《西南交通大学学报》【年(卷),期】2012(047)002【总页数】6页(P224-229)【关键词】三维波浪与结构物相互作用;流体体积(VOF)法;大涡模拟;两步边界定位法;虚拟边界力法【作者】祝兵;宋随弟;谭长建【作者单位】西南交通大学土木工程学院,四川成都610031;西南交通大学土木工程学院,四川成都610031;上海建筑科学研究院上海市工程结构新技术重点实验室,上海200032【正文语种】中文【中图分类】TV139.2随着海洋开发和交通事业的发展，修建了大量海洋建筑物，如海洋平台、跨海大桥等．因此，波浪与结构物相互作用的研究越来越受到各国学者的重视．研究波浪与结构物相互作用不但可以了解结构物周围流场的特性，而且通过仿真计算可以获得波浪荷载，以便研究结构物在波浪作用下的振动响应．国内外学者对波浪与结构物的相互作用进行了大量研究：Lin和Li采用大涡模拟技术，在垂向进行σ坐标变换，用有限差分法离散方程，研究了波浪或波浪-流与方柱的相互作用［1-3］;Palma、Zhang等用虚拟边界力法求解流体与结构物的相互作用［4-5］;李玉成等利用波浪弥散关系迭代计算求波向与流向的夹角，并用有限元法求解含流的缓坡方程，得到在缓变地形和定常流场共同影响下的大尺度圆柱周围波流场的耦合解，并对波流共同作用下大尺度圆柱墩群周围的波流场进行了数值研究［6-7］．本文基于两相流概念，采用大涡模拟的纳维-斯托克斯方程作为波浪运动的基本控制方程，采用两步边界定位法和虚拟边界力法确定波浪与结构物的接触面，自由水面追踪采用分段线性近似的流体体积法，建立了三维波浪与结构物相互作用数学模型，对三维波浪作用下大直径圆柱的绕流现象进行了数值研究．1 波浪运动的基本控制方程计算中，自由表面采用两相流模型，即水和空气同时参与计算，并且假定水和空气都是不可压缩的．采用纳维-斯托克斯方程组描述不可压缩流体运动：式中：i，j=1，2，3，代表三维流体的 3 个方向;ui为i方向的流速;p为压力;gi为i方向的重力加速度;ρ为流体密度(空气中，ρ=ρa;水中，ρ=ρw);t为时间;τij=2ργsij为分子粘性应力张量的分量，其中γ为动力粘度(空气中，γ=γa;水中，γ=γw)，sij为应变率张量的分量，采用大涡模拟方法［8］，首先对方程(1)和(2)进行空间滤波，经过Box滤波器［9］滤波，滤波后的方程去掉滤波符号为：式中，γt为湍流动力粘性系数．进而采用 Smagorinsky 假设［10］，有式中：为局部应变率;C s为常数;Δ为滤波宽度，通常取Δ =(ΔxΔyΔz)1/3，其中Δx、Δy 和Δz分别为在x、y和z方向上的网格宽度．连续方程仍为式(1)．2 结构物界面处理采用两步边界定位法和虚拟边界力法对结构物界面进行处理．在虚拟边界力法中，结构物被移出计算区域，在结构物与流体交界处，采用虚拟边界力f VBF作用到流体上，如图1所示．图1 虚拟边界力法中用虚拟边界力代替任意形状结构物表面的边界条件Fig．1 Visual boundary force replacing the boundary condition on body surface in the virtual boundary force method对于形状简单的物体，可以通过简单识别定位边界线，而对于形状复杂的物体，则较难确定．为此，本文提出了两步边界定位法：第1步：采用结构边界函数初步定义结构物边界．第2步：根据第1步定义的结构物边界，开发了计算识别程序，判断各物理量所在位置是否在结构物边界上——若在边界上，则设虚拟边界力标记点;否则，通过插值计算各物理量．从而实现了流体与复杂结构物边界的精确处理，简单阐述如下．采用一系列二次函数来定义结构物边界面，这些二次函数定义的截面可以相交或叠加，以便能够模拟任意结构形状．二次函数定义如下：式(7)中，除自变量x、y和z外，其余均为构造函数所用的参数，根据结构的实际形状进行调整．从式(7)可知，要与结构物边界吻合，应有f(x，y，z)=0．当 f(x，y，z) ＜0 时，表明计算点在结构物内部，为非流体区域;当 f(x，y，z)＞0时，则表明计算点在结构物外部，为流体区域．为初步定义流体与结构区域，引入函数ψ，ψ在x、y、z方向分别定义．以在x方向定义为例：ψx以结点(x，y，z)所在网格位置y和z方向的网格中心计算得到．用同样的方法，可以求得ψy和ψz．根据ψ值，便可以初步判定计算结点是否为流体区域内的结点．但一些形状复杂的结构物或计算域中存在多个结构时，有较多边界网格，边界网格里既有水，也有结构物，而计算参数如速度、压力等不是在网格中心，就是在网格面的中心．因此，第2步需要开发计算识别程序进行判断，以便精确确定界面的位置，求解结构物边界处的各物理量，详见文献［11］．3 数值实现3．1 控制方程的差分离散及虚拟边界力求解通过以上计算，准确定位了结构物的界面，现在关键的问题是如何求解虚拟边界力．在方程(4)中，加入虚拟边界力f VBF i，即有计算区域的空间离散采用交错网格，压强p、湍流动力粘性系数γt定义在控制体的中心，速度定义在控制体面的中心，采用两步映射法［12-13］求解控制方程．不考虑压强梯度对速度场的贡献，即在动量方程中略去压力梯度项，引入中间速度˜ui，且时间项采用向前差分，得式中：Δt为时间步长;上标(n)、(n+1)表示时间迭代步;i，j=1，2，3，此处x1=x，x2=y，x3=z，u1=u，u2=v，u3=w，分别代表三维流体的3个方向和3个方向的速度．一般来讲，中间速度˜ui不再满足连续方程．将中间速度˜ui映射到速度散度为0的平面上，从而可以得到该时间步最终的速度u(n+1)i ：对式(11)取散度，并应用式(12)，可得修正的压力泊松方程：方程(13)中，若f VBF i已知，则可以进行求解．由于f VBF i只在单元与结构物边界交界处不为0，因此，可以利用式(11)定义f VBF i：式中为满足非滑移边界条件的速度．式(14)中，为了在与结构物边界交界的单元上施加非滑移的速度条件，用代替了式(11)中的，将式(14)代入式(13)离散求解，即可求得f VBF i的值．由于f VBF i是结构物对流体的反作用力，因此，作用在结构物上总的流体作用力F t i可以通过绕结构物表面进行体积积分求得：式中：Ω为结构物表面体积(数值模型中，表示与结构物表面相交的单元的体积)．3．2 流体体积(VOF)法追踪自由表面流体密度的输运方程为采用流体体积(volume of fluid，VOF)法追踪自由面变化．对于不可压缩的流体，VOF函数F(x，t)(简写为F)可以看作是量纲为1的流体体积比，即将式(17)代入式(16)，得VOF方法追踪自由表面，采用精度较高的Youngs方法［14］，它采用分段线性逼近自由面，在每一个控制体内用一个斜线段表示自由面．3．3 边界条件的处理入流边界自由面的位移和速度根据波浪理论的解析解赋值，为消除二次反射波，采用弱反射边界条件．压力边界采用第二类边界条件：式中，g w x为左侧入流处x方向的重力加速度．出流边界为辐射边界条件，压力边界仍采用第二类边界条件(式(19))，其中重力加速度改为右侧出流处x方向的重力加速度即可．平行波浪传播的两侧边界的速度采用滑移边界条件，压力边界仍采用第二类边界条件(式(19))，其中重力加速度改为两侧边界处y方向的重力加速度．计算域上部设置速度边界时，根据VOF方法的计算值采用滑移边界．压力边界采用第一类边界条件，即p=0．底部速度边界同样采用滑移边界，压力边界采用第二类边界条件(式(19))，其中重力加速度改为底部边界处z方向的重力加速度．4 大直径圆柱波浪绕射4．1 理论分析波浪传播过程中与大尺度结构物相遇时，由于结构物的影响，会产生绕射现象，即波浪的一部分能量被结构物阻挡，另一部分能量从结构物侧面绕过，传入结构物后面的水域中．对于大尺度结构物对波浪场的影响以及波浪力的计算等，都需要研究波浪绕射问题．MacCamy和Fuchs给出了大直径圆柱波浪绕射速度势的理论解［15］：式中：φ为大直径圆柱波浪绕射速度势;s为自由水面高度;h为静水深度;δm为待定系数(m=0时，δm=1;m≥1时，δm=2);H为汉克尔函数;H(1)m 为第一类m阶汉克尔函数;Jm为第一类m阶贝塞尔函数;k为波数;ω为波的频率;r和θ为极坐标;a为大直径圆柱的半径．根据式(20)，可以得到波浪绕圆柱的绕流系数式中，H0为0阶汉克尔函数．4．2 参数设置图2 大直径圆柱波浪绕流的平面布置及网格划分Fig．2 Local mesh deployment and plane layout of wave diffraction around a large-diameter circular cylinder大直径圆柱波浪绕射参数：圆柱直径为80 m，计算水域水平面取1 080 m×1 000 m，计算水深20 m，线性波波高0．3 m，周期为8 s．圆柱位于计算域中心(图2)，圆柱足够高，计算中柱顶不淹没．计算网格划分时，在水平方向覆盖圆柱的区域附近，采用2 m×2 m的均匀网格，网格数为272×92;竖向37个网格，自由水面附近采用均匀网格．空间步长为0．03 m;时间步长可自动调整，计算总时间为200 s．波浪为线性波，采用边界造波法造波，从左侧传入．出流边界采用辐射边界条件，竖向取22 m．圆柱距离左、右边界的距离均为500 m，大于3倍波长．两侧边界离圆柱足够远，从而使边界造波和幅射边界对波浪和圆柱的绕射影响较小．计算域剖分时，空间3个方向均采用非均匀网格．4．3 结果分析对边界造波法和自由面追踪方法进行验证，数值计算结果与线性波理论解的比较见图3．图4为三维波面．从图3(a)可见，在边界处，数值计算结果与理论解非常吻合．图3(b)给出了整个计算域波浪自由面的数值计算结果与理论解的比较，可见，除刚起波时稍有差别外，整个计算域波面位置的数值计算结果与理论解吻合很好．图3 线性波数值计算结果与理论解的比较Fig．3 Comparison between numerical result and theoretical solution of linear wave图4 线性波的三维波面(t=180．012 9 s)Fig．4 Three-dimension free surface for linear wave边界造波法和自由面追踪方法验证后，研究大直径圆柱三维波浪的绕射现象．图5给出了t=70．003 1和180．000 6 s时的波面图．从图 5(a)可见，波浪刚刚传播到圆柱时，圆柱的存在影响了波浪的波形，在圆柱正对来波的方向波浪产生爬高现象;从图5(b)可见，圆柱周围有很明显的反射波纹，以圆柱为中心向两边扩散，圆柱正对来波方向波浪的波幅明显增大，靠近圆柱处波浪产生爬高现象，圆柱背浪侧波浪的波幅减小，圆柱后侧形成了清晰的绕射波纹．图5 三维波面Fig．5 Three-dimension free surface为了与理论上的绕流系数比较，在圆柱周围设置了一系列波高测点，计算波浪充分发展后绕流系数随θ的变化(θ见图2)．用式(21)计算理论绕流系数，并与数值计算得到的绕流系数进行比较．从图6可见，绕流系数的数值计算结果与理论值非常接近，表明本文建立的波浪数学模型不仅能够反映波浪的反射特性，还可以反映绕射波高的变化，可以成功地模拟三维波浪与结构物的相互作用．图6 绕流系数理论值与数值计算结果的比较Fig．6 Numerical and analytical results of circumfluence factor for circular cylinder5 结论以两相流概念、大涡模拟的不可压缩粘性流体运动方程和自由水面追踪分段线性近似的流体体积法为基础，建立了三维波浪与结构物相互作用的数学模型，并对该模型进行了验证，得到以下结论：(1)对边界造波法和流体体积法追踪自由面进行数值计算验证，所得波浪的数值计算结果与理论解非常吻合，实现了波浪自由面的高精度跟踪．(2)结合虚拟边界力法，提出了一种新的流体与结构物接触面的数值模拟方法——两步边界定位法，实现了流体与复杂结构物边界的精确模拟．(3)大直径圆柱波浪绕流系数的数值计算结果与理论值非常接近，所建立的数学模型不仅能够反映波浪的反射特性，还可以反映绕射波高的变化，成功地模拟三维波浪与结构物的相互作用．【相关文献】［1］ LIN P．Numerical modeling of water waves［M］．New York：Taylor&Francis，2008：324-384．［2］ LIN P，LI CW．Wave-current interaction with a vertical square cylinder［J］．Ocean Engineering，2003，30(7)：855-876．［3］ LI Chiwai，LIN Pengzhi．Numerical study of threedimensional wave interaction with a square cylinder［J］．Ocean Engineering，2001，28(12)：1545-1555．［4］ PALMA P D，TULLIO M D．PASCAZIO G，et al．An immersed-boundary methodfor compressible viscous flows［J］．Computers&Fluids，2006，35(7)：693-702．［5］ ZHANG N，ZHENG Z C．An improved direct-forcing immersed-boundary method for finite difference applications［J］．Journal of Computational Physics，2007，221(1)：250-268．［6］李玉成，刘德良，陈兵，等．大尺度圆柱周围的波流场的耦合计算［J］．力学研究与进展，2004，19(2)：168-173．LI Yucheng，LIU Deliang，CHEN Bing，et al．Wave diffraction in a current around a large-scale cylinder［J］．Journal of Hydrodynamics，2004，19(2)：168-173．［7］李玉成，刘德良，陈兵，等．大尺度圆柱墩群周围的波流场的数值模拟［J］．海洋通报，2005，24(2)：1-12．LI Yucheng， LIU Deliang， CHEN Bing， et al．Interaction of waves with vertical large-scale cylinders in a current［J］．Marine Science Bulletin，2005，24(2)：1-12．［8］ SAGAUT P．Large eddy simulation for incompressible flows［M］．Berlin：Springer，2001：5-18．［9］ DEARDORFF J W． A numerical study of threedimension turbulent channel flow at large Reynolds numbers［J］． Journal of Non-Newtonian Fluid Mechanics，1970，41(2)：453-480．［10］ SMAGORINSKY J S．General circulation experiments with the primitive equations：Part I［J］．The Basic Experiment：Monthly Weather Review，1963，91：99-163．［11］谭长建．波浪效应下大跨度斜拉桥索桥耦合振动研究［D］．成都：西南交通大学土木工程学院，2009．［12］ CHORIN A J．Numerical solution of the Navier-Stokes equations ［J］．Mathematics of Computation，1968(22)：745-762．［13］ CHORIN A J． On the convergence of discrete approximations of the Navier-Stokes equations［J］．Mathematics of Computation，1969(23)：341-353．［14］ YOUNGS D L．Time-dependent multi-material flow with large fluid distortion［C］∥Numerical Methods for Fluid Dynamics．New York：Academic，1982：273-285．［15］ MACCAMY R C，FUCHSR A．Wave forces on piles：a diffraction theory［J］． Beach Erosion Board Technical Memorandum，1954，69：1-17．。
圆柱振荡绕流LBM数值模拟的开题报告
圆柱振荡绕流LBM数值模拟的开题报告一、研究背景及意义圆柱绕流问题是流体力学中经典的研究问题之一，在空气动力学中有着广泛的应用，如飞机、桥梁、烟囱、塔楼等建筑结构受风振动问题。
目前对于圆柱绕流问题的研究主要采用实验和数值模拟相结合的方法。
在数值模拟方面，近年来LBM（Lattice Boltzmann Method）方法因其简单低耗，适用于各种欧拉和纳维-斯托克斯方程问题而受到广泛研究。
LBM方法通过在非均匀格子上演化局域平衡分布函数来求解流动问题，实现了传统CFD方法的高效计算，并且在高雷诺数下仍能保持很好的数值稳定性。
本研究旨在采用LBM方法对圆柱绕流问题进行数值模拟，并研究圆柱粘性系数、柱直径以及雷诺数对流动特性的影响，为实际应用提供理论依据。
二、研究方法1.建立数学模型：通过欧拉方程、斯托克斯方程及纳维-斯托克斯方程对圆柱绕流问题进行数学建模。
2.确定求解方法：采用LBM方法对圆柱绕流问题进行数值模拟，并通过MATLAB编程实现模型求解。
3.数据分析：采用图像分析方法对模拟结果进行可视化分析，并将模拟结果与理论计算结果进行比较分析。
三、研究内容及进度安排研究内容主要包括：1.建立圆柱绕流LBM模型，并进行数值模拟；2.分析圆柱粘性系数、柱直径以及雷诺数对流动特性的影响；3.比较与理论计算的结果，并对模拟结果进行可视化分析；研究进度安排如下：1.文献调研及分析（1个月）；2.建立圆柱绕流LBM模型（2个月）；3.进行数值模拟并分析结果（3个月）；4.编写论文及进行答辩（1个月）。
四、研究预期成果1.建立了圆柱绕流LBM模型，通过数值模拟研究了圆柱粘性系数、柱直径以及雷诺数对流动特性的影响。
2.编写论文，将研究结果公开发表，并与同行分享。
3.为圆柱绕流问题的理论研究及实际应用提供理论依据。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
圆柱绕流的数值模拟 张玉静 20070360204 过控（2）班 化工与能源学院 摘 要:使用计算流体力学软件FLUENT,模拟均匀来流绕固定圆柱的流动,模拟雷
诺数为5，20,40,100时的绕流流动, 得到流场的流函数等值线图和速度矢量图。计算结果表明:当雷诺数增加时,流动表现出一系列不同的构造。当Re=5时，流动不发生分离,其后未形成旋涡，当Re=20,40,100时，流体发生分离，其后形成旋涡，且旋涡随着Re的增大而增大。利用计算流体力学软件FLUENT可以成功地模拟圆柱绕流问题,反映出流动特性。 关键词:圆柱绕流；FLUENT；雷诺数 Abstract:Uniform flow around a mounting cylinder is simulated with the application of FLUENT software while Reynolds number is 5,20,40,100. Stream function and velocity vector distributions are indicated. The results show that a series of construction appears as Reynolds number increases. When Re is 5, Flow separation does not occur, and it does not form vortex . When Re is 20,40,100, Flow separation occurs, and it forms vortex. Vortex increases with the increase of Re. Using computational fluid dynamics software FLUENT can successfully simulate flow around cylindrical, reflect the flow characteristic. Key words：Flow around a circular cylinder；FLUENT；Reynolds number 1 圆柱绕流理论分析研究的状况 一个世纪以来,圆柱绕流一直是众多理论分析、实验研究及数值模拟对象。但迄今对该流动现象物理本质的理解仍是不完整的。圆柱绕流中,起决定作用的是雷诺数,但还受到许多因素,如阻塞比,来流湍流度,下游边界条件等的影响。随着雷诺数的增加,粘性不可压缩流体绕圆柱的流动会呈现各种不同的流动状态,在小雷诺数时,流动是定常的,随着雷诺数的增加,圆柱后会出现一对尾涡。当雷诺数较大时,尾流首先失稳,出现周期性的振荡。而后附着涡交替脱落,泻入尾流形成Karman涡街,随着雷诺数的增加,流动变得越来越复杂,最后发展为湍流。White认为圆柱涡流具有经典性的重要意义。 一般认为圆柱绕流有2种定常的流动图案:雷诺数为较小时,圆柱后无尾涡；当雷诺数为较大时,圆柱后有一对对称的尾涡。关于定常流失稳以及出现湍流的临界雷诺数主要是通过应用流场显示技术观察流动形态得到的,所以不是准确值。对于分界点雷诺数就有不同的见解,Kovasznay,Roshko等认为定常流动失稳的临界雷诺数大约为40。而从周期性尾流到湍流的详细的转变过程的实验研究似乎还是空白。对均匀来流绕固定圆柱的二维平面流动,国内外许多学者进行过大量的研究。决定圆柱绕流流态的是雷诺数(Re)的值,Re<5时,流动不发生分离, 5一个旋涡开始脱落,接着另一个也脱落,在圆柱体后面又生成新的旋涡,这样逐渐发展成两排周期性摆动和交错的旋涡,即Karman涡街。Re<150,涡街是层流,150柱体表面上的边界层为层流,而柱体后面的涡街已完全转变为湍流,并按一定的频率发放旋涡。3×105目前,对圆柱绕流的三维数值研究则不多,然而试验表明,即使均匀来流垂直流过等截面的圆柱体,当Re足够大时,也会呈现三维的流动状态。圆柱绕流属于非定常分离流动问题,在工业工程中的应用非常广泛,数值模拟是研究这类问题的有力工具。但是,控制方程的高度非线性以及边界条件的多样性严重阻碍了这类问题的解决。到目前为止,有限差分、有限元、大涡模拟(LES)、直接数值模拟(DNS)等方法都先后被应用。但是,这些方法都需要增加网格生成的附加工作,尤其是在遇到复杂几何区域时更需要巨大的网格数量,大大影响计算效率。 本文利用N-S方程,对固定圆柱绕流进行了三维数值模拟,利用计算流体力学软件FLUENT,选取Re=20,40,100,模拟了圆柱周围的流场,得到流场的涡量等值线和速度矢量图。 2 数学模型与数值方法 对不可压缩粘性流体,在直角坐标系下,其运动规律可用N-S方程来描述,连续性方程和动量方程分别为： （1）连续性方程
0zwyvx
u
（1）
式中，u、v、w表示速度在x、y、z方向上的分量。 （2）动量方程 jiiiijjixxuxpfxuu
21)( （4）
式中，xi（i=1,2,3）为坐标系坐标，xj(j=1,2,3)为坐标系坐标，ui
(i=1,2,3)为
沿i方向的速度分量，fi
为沿i方向的质量力，P是静压，ρ为空气密度。
3 模型建立与网格划分
一个无穷长直径为2.0m的圆柱体,放置在无穷远来流速度为1.0m/s不受干扰的均匀横流中,计算区域如图1所示。网格划分采用四边形网格，网格数为16000。Gambit中网格划分结果如图2所示。
图1 计算区域 图2 网格划分 4边界条件与求解设置
进口为Velocity-inlet，出口为Outflow，其他边界默认为Wall，计算区域设为fluid，在gambit中设置界面如图3所示。
图3边界条件置界面 图4指定流体区域
求解设置如下： （1）控制方程采用连续方程和动量方程，而不考虑能量方程。 （2）求解器选择COUPLED IMPLICIT SOLVER，选择Steady状态。 （3）为了提高计算精度，差分格式采用二阶迎风格式。 （4）速度场和压力场的求解基于Simple算法，松弛因子默认。 （5）模型选择层流流动。 （6）残差设为e-6。 在fluent中的操作具体如下： （1）读入从gambit中导出的msh文件，检查网格质量，并确定模型单位，如图5所示。 （2）设置求解器，如图6所示。 （3）选择层流模型，如图7所示。 图5 确定单位 图6 选择求解器 图7 选择层流模型 （4）改变流体介质，如图8所示。 （5）设置入口速度，如图9所示。 （6）设置残差，如图10所示。 （7）初如化流场，进行迭代求解。
图8 改变流体介质 图9 设置入口速度
图10 设置残差 5数值计算与结果讨论 雷诺数由圆柱体直径和自由来流速度确定。通过改变粘度,得到不同的雷诺值(5，20,40,100),Re =ρUD/μ。 边界条件:进口和侧面边界始终为自由来流条件U=1.0m/s,流动出口为零法向梯度出口边界。 流体特性假定为常数,结果见表1。 表1 流体特性 Re 5 20 40 100
密度ρ/(kg/m3) 1 1 1 1
动力粘度μ/(Pa.s) 0．4 0.1 0.05 0.02 Re=20时的残差曲线如图11所示。
图11 Re=20时计算的残差曲线 (a) Re=5时的速度分布 （b）Re=20时的速度分布 (c) Re=40时的速度分布 (d) Re=100时的速度分布 图12不同Re下的速度分布 （a）Re=5时的流线分布 (b)Re=20时的速度分布 (c) Re=40时的流线分布 (d)Re=100时的流线分布 图13不同Re下的流线分布 图12为不同Re下的速度分布，图13为不同Re下的流线分布，由图12、
图13可以看出流体流过圆柱体时，受圆柱体影响，圆柱后面会形成一低速区。当Re=5时，流动不发生分离,其后未形成旋涡，当Re=20,40,100时，流体发生分离，其后形成旋涡，且旋涡随着Re的增大而增大。 6 结论 本文使用计算流体力学软件FLUENT,模拟均匀来流绕固定圆柱的流动,模拟雷诺数为5，20,40,100时的绕流流动,得到流场的速度分布云图和流线图。计算结果表明: （1）当雷诺数增加时,流动表现出一系列不同的构造。当Re=5时，流动不发生分离,其后未形成旋涡，当Re=20,40,100时，流体发生分离，其后形成旋涡，且旋涡随着Re的增大而增大。 （2）利用计算流体力学软件FLUENT可以成功地模拟圆柱绕流问题,反映出流动特性。
参考文献: [1]酆庆增.圆柱绕流的非线性动力学[J].力学进展,1994, 24(4):525~544. [2]苏铭德,康钦军.亚临界雷诺数下圆柱绕流的大涡模拟[J].力学学报,1999,31(1):100~105. [3]叶春明,吴文权.数值模拟圆柱绕流旋涡生成、分离及演化[J].华东工业大学学报,1995,17(4):25~30. [4]叶春明,吴文权.数值模拟圆柱绕流旋涡运动及尾流不稳定性分析[J].工程热物理学报,1997,18(2):169~172. [5]陈斌,郭烈锦,杨晓刚.圆柱绕流的离散涡数值模拟[J].自然科学进展,2002,12(9):964~969. [6]万德成.用水深平均雷诺方程模拟有限长直立圆柱绕流 [J].上海大学学报,1995,1(3):259~268.