圆柱绕流理论研究和数值模拟

圆柱绕流的数值模拟张玉静 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. V ortex 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 number1 圆柱绕流理论分析研究的状况一个世纪以来,圆柱绕流一直是众多理论分析、实验研究及数值模拟对象。

不同雷诺数下的圆柱绕流数值模拟研究引言：圆柱绕流是流体力学领域中一个经典的、被广泛研究的问题。

在众多的工业应用中，圆柱绕流的研究对于风力发电机组的设计优化、管道内部液体运动的控制等方面具有重要实际意义。

雷诺数是描述流体流动的一个无量纲参数，它与流体的流速、流体的粘性有关。

本文将对不同雷诺数下的圆柱绕流进行数值模拟研究。

方法：数值模拟是一种有效的研究流体力学问题的方法，它能够通过计算机模拟得到流体的速度场、压力场等关键参数，从而进一步分析流体的特性。

在本文中，我们将使用计算流体力学方法进行圆柱绕流的数值模拟研究。

结果与讨论：我们选取了不同雷诺数的圆柱绕流作为研究对象，分别为200、400、600、800和1000，通过数值模拟得到了不同雷诺数下的圆柱绕流的速度场和压力场等关键参数。

首先，我们分析了速度场的分布。

通过数值模拟可以得到圆柱绕流过程中流体速度的分布情况。

随着雷诺数的增加，流体速度场呈现出不同的特征。

在雷诺数较低的情况下，流体绕圆柱流动的速度场分布较为简单，流速主要集中在圆柱前部和尾部。

随着雷诺数的增加，流体速度场呈现出更复杂的结构，流速分布更加均匀。

其次，我们研究了压力场的分布。

通过数值模拟可以得到圆柱绕流过程中流体压力的分布情况。

在不同雷诺数下，圆柱周围存在不同的压力区域。

当雷诺数较低时，圆柱前后表面存在较大的压差，压力分布较为不均匀。

而当雷诺数增加时，压力分布更加均匀，圆柱表面的压力变化较小。

最后，我们研究了绕流过程中的阻力情况。

通过数值模拟得到了不同雷诺数下圆柱绕流过程中的阻力系数。

我们发现，随着雷诺数的增加，阻力系数逐渐增大。

这是因为当雷诺数较低时，流体绕圆柱流动的速度较低，阻力较小；而当雷诺数增加时，流体流动速度较高，阻力也逐渐增大。

结论：本文通过数值模拟的方式研究了不同雷诺数下的圆柱绕流问题。

通过分析速度场、压力场和阻力系数等关键参数，我们得出了以下结论：随着雷诺数的增加，流体速度场更加复杂，流速分布更加均匀；压力场分布更加均匀，圆柱表面的压力变化较小；阻力系数随着雷诺数的增大而增加。

流体力学Fluent报告——圆柱绕流亚临界雷诺数下串列双圆柱与方柱绕流的数值模拟摘要：本文运用Fluent软件中的RNG k-ε模型对亚临界雷诺数下二维串列圆柱和方柱绕流问题进行了数值研究，通过结果对比，分析了雷诺数、柱体形状对柱体绕流阻力、升力以及涡脱频率的影响。

一般而言，Re数越大，方柱的阻力越大，圆柱体则不然；而Re越大，两种柱体的升力均越大。

相对于圆柱，同种条件下，方柱受到的阻力要大；相反地，方柱涡脱落频率要小。

Re越大，串列柱体的Sr数越接近于单圆柱体的Sr数。

关键字：圆柱绕流、升力系数、阻力系数、斯特劳哈尔数在工程实践中，如航空、航天、航海、体育运动、风工程及地面交通等广泛的实际领域中，绕流研究在工程实际中具有重大的意义。

当流体流过圆柱时 , 由于漩涡脱落，在圆柱体上产生交变作用力。

这种作用力引起柱体的振动及材料的疲劳，损坏结构，后果严重。

因此，近些年来，众多专家和学者对于圆柱绕流问题进行过细致的研究，特别是圆柱所受阻力、升力和涡脱落以及涡致振动问题。

沈立龙等[1]基于RNG k⁃ε模型，采用有限体积法研究了亚临界雷诺数下二维圆柱和方柱绕流数值模拟，得到了圆柱和方柱绕流阻力系数C与 Strouhal 数d随雷诺数的变化规律。

姚熊亮等[2]采用计算流体软件CFX中LES模型计算了二维不可压缩均匀流中孤立圆柱及串列双圆柱的水动力特性。

使用非结构化网格六面体单元和有限体积法对二维N- S方程进行求解。

他们着重研究了高雷诺数时串列双圆柱在不同间距比时的压力分布、阻力、升力及Sr数随Re数的变化趋势。

费宝玲等[3]用FLUENT软件对串列圆柱绕流进行了二维模拟，他们选取间距比L/D(L为两圆柱中心间的距离，D为圆柱直径)2、3、4共3个间距进行了数值分析。

计算均在 Re = 200 的非定常条件下进行。

计算了圆柱的升阻力系数、尾涡脱落频率等描述绕流问题的主要参量，分析了不同间距对圆柱间相互作用和尾流特征的影响。

圆柱绕流的数值模拟一、问题简介我们考虑一个固定的无限长圆柱体，其直径为10mm，空气以均匀的速度由远处而来绕过圆柱，气流会在圆柱后发展为复杂的流动。

这是一个经典的流体力学问题，随雷诺数的增加，柱体后的流动形态会由对称向不对称转变，并产生卡门涡街。

我对不同雷诺数下的流动进行了数值模拟，并对计算所得流场进行了比较和分析。

二、文献综述圆柱绕流作为最为常见的钝体绕流现象，演绎出了大量的流体控制工程技术和理论研究课题。

这类问题常见的有风掠过建筑物，气流对电线的作用，海流冲击海底电缆，河水对桥墩的冲击，气流经过冷凝器中的排管、空中加油机的油管以及飞行器上的柱体等等，具有很高的工程实践意义。

同时圆柱绕流又是流体力学的经典问题，其蕴含了丰富的流动现象和深刻的物理机理，长久以来一直是众多理论分析、实验研究及数值模拟的研究对象。

流体绕圆柱体流动时，过流断面收缩，流速沿程增加，压强沿程减小，由于黏性力的存在，就会在柱体周围形成附面层的分离，形成圆柱绕流。

在圆柱绕流问题中，流体边界层的分离与脱落、剪切层的流动和变化、尾迹区域的分布和变动，以及它们三者之间的相互作用等因素，使得该问题成为了一项复杂的研究课题。

圆柱绕流的流动状态主要由雷诺数(Re)决定,根据不同的Re范围，流动会经历多种流动状态，在我们流体力学的教材上，就可以查到不同雷诺数下圆柱绕流的形态变化，而下表更加完整详细。

表一在使用CFD方法对圆柱绕流进行求解时，早期使用求解二维定常N-S方程的方法来模拟绕流流场。

然而，由于圆柱尾部涡脱落的存在，绕流流场随时间在不断改变，具有非定常特性，因此就需要求解非定常N-S方程。

目前，在低雷诺数层流条件下，多以求解二维非定常N-S方程来研究圆柱绕流。

但随着雷诺数的增加，绕流流场中沿展向的三维特性越来越显著，如果还使用二维计算模型求解流场，必然不能正确的解析流场结构，获得正确的流场参数。

所以在大雷诺数条件下就需要求解三维的N-S方程。

小雷诺数下圆柱绕流数值模拟
小雷诺数下圆柱绕流数值模拟是一种常用的研究绕流流动特性的方法。

它是一种基于小雷诺数下流体运动的数值模拟，用于研究圆柱体绕流的结构及其流动特征。

小雷诺数下圆柱绕流数值模拟的基本思想是，利用小雷诺数的概念，通过模拟圆柱体的绕流流动特性，求得其相应的流速、压力、温度等各种流动参数。

小雷诺数的定义是在恒定的空气压力下，流体的绕流流动速度和其局部粘性系数之比，其值与空气压力以及温度有关。

圆柱绕流数值模拟主要分为三个步骤：（1）建立基本流场模型，确定模型空间的尺寸和流体的流动参数；（2）使用计算流体力学（CFD）技术计算圆柱体绕流流动的模拟结果；（3）分析模拟结果，确定圆柱体绕流流动的特性以及绕流流动的特性参数。

小雷诺数下圆柱绕流数值模拟是一种有效的研究圆柱体绕流流动特性的方法。

它可以用来模拟圆柱体绕流流动的特性，并及时获取流动参数，从而更好地分析绕流流动特性，为设计提供参考依据。

亚临界雷诺数下串列单圆柱与圆柱绕流的数值模拟之阳早格格创做目要：原文使用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战错误!未找到引用源。

二维圆柱绕流摘要:采用有限体积法求解二维N -S 方程，对雷诺数1,10,100的二维圆柱非定常流场进行了数值模拟，对比各雷诺数下其流动情况发现，在Re=1时，圆柱上下游的流线前后对称，此Re 数范围的绕流称为斯托克斯区；随着Re 的增大，圆柱上下游的流线逐渐失去对称性。

当Re=10时，沿圆柱表面流动的流体在到达圆柱顶点（90度）附近就离开了壁面，分离后的流体在圆柱下游形成一对固定不动的对称漩涡（附着涡），涡内流体自成封闭回路而成为“死水区”；随着Re 的增大，死水区逐渐拉长圆柱前后流场的非对称性逐渐明显，此Re 数范围称为对称尾流区。

圆柱下游流场不再是定常的，圆柱后缘上下两侧有涡周期性地轮流脱落，形成规则排列的涡阵，这种涡阵称为卡门涡街。

1. 圆柱绕流研究圆柱绕流是一个经典的流体力学问题，流体绕圆柱体流动时，过流断面收缩，流速沿程增加，压强沿程减小，由于粘性力的存在，，就会在柱体周围发生边界层的分离，形成圆柱绕流。

圆柱绕流现象比较复杂，因此，对圆柱绕流研究具有重要的基础理论意义。

研究圆柱绕流问题在工程实际中也具有非常重要的意义。

如水流对桥梁、海上钻井平台支柱、海上输运管线、桩基码头等的作用中，风对塔设备、化工塔设备，高空电缆等的作用中，都有着重要的工程应用背景。

因此，对圆柱绕流进行深入研究，对其流动机理进行分析，不仅具有理论意义，还有明显的社会经济效益。

2. 数值方法因为本文主要求解雷诺数Re=1,10,100时的圆柱绕流情况，需要求解二维非定常不可压的N—S 方程组：本文采用有限体积法对上述微分方程进行离散，然后用SIMPLE 算法对离散方程进行求解，计算中时间推进采用二阶隐式格式，空间离散采用三阶精度的QUICK 格式。

控制方程如下：0jju x ∂=∂ (1) 1()()ji j i j j j ju u p u u t x u x x νρ∂∂∂∂∂+=-+∂∂∂∂∂ (2) 3. 网格划分及模拟工况3.1计算网格计算的区域大小为上游边界距圆柱圆心为2.5D ，下游边界距圆柱圆心10D ，顶部和底部边界距圆柱圆心2.5D ，如图1所示。

圆柱绕流的数值模拟研究摘要：选取直径为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所示。

二维圆柱绕流摘要:采用有限体积法求解二维N-S 方程，对雷诺数1,10,100的二维圆柱非定常流场进行了数值模拟，对比各雷诺数下其流动情况发现，在Re=1时，圆柱上下游的流线前后对称，此Re 数范围的绕流称为斯托克斯区；随着Re 的增大，圆柱上下游的流线逐渐失去对称性。

当Re=10时，沿圆柱表面流动的流体在到达圆柱顶点（90度）附近就离开了壁面，分离后的流体在圆柱下游形成一对固定不动的对称漩涡（附着涡），涡内流体自成封闭回路而成为“死水区”；随着Re 的增大，死水区逐渐拉长圆柱前后流场的非对称性逐渐明显，此Re 数范围称为对称尾流区。

圆柱下游流场不再是定常的，圆柱后缘上下两侧有涡周期性地轮流脱落，形成规则排列的涡阵，这种涡阵称为卡门涡街。

1. 圆柱绕流研究圆柱绕流是一个经典的流体力学问题，流体绕圆柱体流动时，过流断面收缩，流速沿程增加，压强沿程减小，由于粘性力的存在，，就会在柱体周围发生边界层的分离，形成圆柱绕流。

圆柱绕流现象比较复杂，因此，对圆柱绕流研究具有重要的基础理论意义。

研究圆柱绕流问题在工程实际中也具有非常重要的意义。

如水流对桥梁、海上钻井平台支柱、海上输运管线、桩基码头等的作用中，风对塔设备、化工塔设备，高空电缆等的作用中，都有着重要的工程应用背景。

因此，对圆柱绕流进行深入研究，对其流动机理进行分析，不仅具有理论意义，还有明显的社会经济效益。

2. 数值方法因为本文主要求解雷诺数Re=1,10,100时的圆柱绕流情况，需要求解二维非定常不可压的N —S 方程组：本文采用有限体积法对上述微分方程进行离散，然后用SIMPLE 算法对离散方程进行求解，计算中时间推进采用二阶隐式格式，空间离散采用三阶精度的QUICK 格式。

控制方程如下：0j ju x ∂=∂ (1)1()()j i j i j j j ju u p u u t x u x x νρ∂∂∂∂∂+=-+∂∂∂∂∂ (2)3. 网格划分及模拟工况3.1计算网格计算的区域大小为上游边界距圆柱圆心为2.5D ，下游边界距圆柱圆心10D ，顶部和底部边界距圆柱圆心2.5D ，如图1所示。

二维圆柱绕流数值仿真模拟报告
二维圆柱绕流数值仿真模拟是流体力学领域的一个重要研究课题。

在这个问题中，我们将从多个角度来讨论这个课题，包括数值
模拟的背景和意义、数值方法、模拟结果及其分析等方面。

首先，让我们来谈谈数值模拟的背景和意义。

在工程和科学研
究中，流体力学的数值模拟在许多领域都发挥着重要作用，例如飞
行器设计、汽车空气动力学、海洋工程等。

而二维圆柱绕流作为流
体力学中经典的基本问题，对于理解流体运动、探索流体特性具有
重要意义。

通过数值仿真模拟，我们可以更好地理解圆柱绕流的流
场特性，为相关工程和科学问题提供重要参考。

其次，数值方法是进行二维圆柱绕流数值仿真模拟的关键。

常
用的数值方法包括有限元法、有限体积法和边界元法等。

这些方法
可以通过数值计算求解流体力学方程，如Navier-Stokes方程，从
而得到圆柱绕流的流场分布、压力分布等信息。

在选择数值方法时，需要考虑计算精度、计算效率以及模拟的适用范围等因素。

接下来，我们来讨论模拟结果及其分析。

通过数值仿真模拟，
我们可以得到圆柱绕流的流速场、压力分布、升力和阻力等重要参
数。

通过分析这些参数，我们可以深入理解圆柱绕流的特性，比如卡门涡的形成与演变、升力和阻力的变化规律等。

这些分析结果对于优化工程设计、改进流体力学理论具有重要意义。

总的来说，二维圆柱绕流数值仿真模拟是一个复杂而有意义的课题，通过深入研究和全面分析，可以为相关领域的工程和科学问题提供重要的参考和指导。

希望以上讨论能够对你有所帮助。