首页/文章/ 详情

【经典教材翻译】9-边界层和流动分离

8天前浏览405

1.外流中的粘性效应


     

   

实际上,忽略了粘性流动纯粹是理论上的。在圆柱体周围的流动中,粘性会完全改变圆柱体后部的流动。自然界中任何真实的流动都包含粘性。正是粘性产生了许多流动的有趣物理特性。

靠近表面的流动层受粘性支配。这个层叫做边界层,将是这里的重点。


2.边界层流动


     

   

对于平板上的流动,正好在表面上发生了“无滑移”条件。这意味着流体粘附在表面上,相对于表面不移动。这是粘性的典型效应。

图 1:边界层的形成

在图 1 中,考虑了平板前均匀流动,速度为 U∞。当流动接触到平板时,发生了无滑移条件。结果,物体表面的速率变为零。由于粘性的作用是抵抗流体运动,速度将被表面的零速率流动所减慢。这种作用随着远离平板的距离而减少,直到速度再次等于主流值 U∞。流体中建立了一个垂直于流动方向的速度梯度。因此,靠近墙壁的一层建立了自己。这被称为边界层

边界层并不是一种静态现象。它是动态的。边界层的厚度(从固体表面到流动再次达到主流速度 99% 的高度)沿着平板下游不断增加。固体表面上发展出剪切应力。正是这种剪切应力导致平板上的阻力。

边界层对于任何沉浸在流体中并移动的物体都有显著影响。飞机或船只上的阻力以及管道中的摩擦是边界层的一些常见表现。可以理解,边界层分析已成为流体动力学研究的一个重要分支。


3.层流和湍流边界层


     

   

边界层可以是层流或湍流。层流边界层是流动顺畅且分层的,即每层滑过相邻层。这与下图图 2 所示的湍流边界层形成对比,其中有层之间的混合。

在层流边界层中,任何质量或动量的交换只能在微观尺度上在相邻层之间进行,肉眼看不见。因此,分子粘性 μ 能够预测产生的剪切应力。层流边界层只有在雷诺数很小的时候才会出现。

图 2:层流和湍流边界层的典型速度剖面

另一方面,湍流边界层的特点是流体在多个层中随机混合。混合是在宏观尺度上进行的。可以看到流体包在层中上下移动。因此,与层流边界层相比,有更大尺度上的质量、动量和能量的交换。在较大的雷诺数下形成湍流边界层。混合的尺度不能单独由分子粘性处理。湍流的计算依赖于一个称为湍流粘性Turbulence Viscosity或涡粘性Eddy Viscosity的对流参数。这个没有确切的表达式,所有情况都需要进行数值建模。为了简化这种统计变化的流动,已经开发了几种模型。

图 3:层流和湍流边界层的典型速度剖面

湍流层可以使用平均速度剖面来描述。由于强烈的混合,湍流边界层在壁面处有较大的速度梯度,因此剪切应力较大。此外,热传递率也很高。图 3 显示了典型的层流和湍流边界层剖面。基于这些近似速度剖面,可以计算边界层的生长速率。生长速率表示边界层厚度 δ 沿板的长度增加的速率。层流边界层的生长速率较小。对于平板,它由下式给出:

其中 Re x 是基于平板前缘距离的雷诺数。

对于湍流,它由下式给出:

壁面剪切应力是边界层的另一个重要参数。它通常表示为壁面摩擦系数Skin Friction Coefficient,并定义为:

其中 τw 是壁面剪切应力,由下式给出:

U∞ 是主流速度。

层流和湍流的壁面摩擦分别为:


4.流动分离


     

   

压力梯度是影响边界层生长的另一个因素。由粘性引起的剪切应力对流动有阻碍作用。然而,如果施加负压力梯度,这种效果可以被克服。负压力梯度被称为有利压力梯度Favourable pressure gradient,代表边界层外流动加速的情况。

正压力梯度的效果相反,被称为不利压力梯度Adverse Pressure Gradient。在这种情况下,外部流动正在减速。流体发现很难通过不利的压力梯度。不利压力梯度显著增加了边界层的生长速率。

图 4:曲表面上流动的分离

不利压力梯度的一个严重影响是使流动从物体表面分离。如图 4 所示的曲表面上的流动,不利压力梯度是由表面形状引起的。表面几何形状是这样的,最初有一个有利的压力梯度,直到点 P。负压力梯度将抵消边界层中剪切应力的阻碍作用。然后,在 P 下游有一个不利的压力梯度。

不利压力梯度开始阻碍流动。这种效果在靠近壁面的区域感受到更强,那里的动量比主流区域的动量低。如图中所示,壁面附近的速率降低,边界层变厚。持续的流动减速使壁面上点 S 处的壁面剪切应力降至零。在这里,壁面附近的流动完全停止了,不仅仅是表面的无滑移流动。从这一点开始,剪切应力变为负值,流动反向,形成了一个回流区域。流动不再沿着物体的轮廓移动。流动已经分离了。剪切应力为零的点 S 称为分离点Point of Separation

根据流动条件,回流可能终止,然后流动可能重新附着在物体上。形成了一个分离泡。有多种因素可能影响这种重新附着。由于物体几何形状的进一步变化,压力梯度可能变得有利。另一个因素是最初是层流的流动可能在泡内过渡并变得湍流。湍流流动在近壁区域产生更多的能量和动量。这可以消除回流区域,流动可能重新附着。

图 5:翼型上的分离泡

在低角度的翼型上,分离可能发生在前缘附近,并通过过渡到湍流产生短泡。在高角度下,分离的流动可能不会重新附着,而是简单地与尾迹区域合并。这将导致翼型的失速(升力损失,阻力增加)。


5.阻力


     

   

阻力是阻碍运动的力。飞行中的飞机必须克服它所受的阻力,飞行中的球,高速行驶的帆船和汽车等是其他一些例子。显然,粘性是导致阻力的介质。它在固体表面上引起边界层。边界层中的剪切应力对物体施加了阻碍力。这在图 6 中的翼型表面上显示,并称为壁面摩擦阻力Skin friction Drag

图 6:物体上的剪切应力

另一种可能导致阻力的机制是流动上的压力差。这可能是由于几何效应和分离造成的。这被称为压力阻力Pressure Drag形状阻力Form Drag,因为它是由于物体的几何形状。

压力阻力和皮肤摩擦阻力的总和构成了物体周围的阻力或轮廓阻力Profile Drag

图 7:物体厚度对阻力的影响

物体的形状决定了阻力分量的相对大小。如图 7 所示,薄物体(小 t/l 比率)显然会引起较小的压力阻力。几乎所有阻力都来自皮肤摩擦。厚物体(大 t/l 比率)容易分离并产生相当大的压力阻力。通过流线化物体以避免分离将大大降低压力阻力。像圆柱体或球体或垂直于流动放置的平板这样的钝物体会导致分离并导致压力阻力远远大于皮肤摩擦阻力。在平板平行于流动的情况下(图 8),是皮肤摩擦阻力占主导地位。另一方面,如果平板垂直于流动放置,是压力阻力占主导地位。在后一种情况下,平板的行为像一个钝物体,在角落处分离,留下后方的低压区域,有助于阻力。

图 8:平板周围的阻力


6.阻力系数


     

   

阻力力被无量纲化为

其中 C D 被定义为阻力系数Drag Coefficient。U∞ 是主流速度,ρ∞ 是主流密度,A 是面积。在此计算中使用的面积取决于应用。在圆柱体的情况下,它是垂直于流动的投影面积。对于薄平板的流动,它将是暴露于流动的板面面积。

两种阻力的相对重要性在圆柱体或球体的流动中非常明显。流动强烈依赖于雷诺数,如图 9 所示。当雷诺数很小(1 或以下)时,流动表现为势流。没有分离。阻力完全由皮肤摩擦引起。随着雷诺数的增加,这种阻力会减小。在雷诺数约为 2 - 30 时,边界层会分离,但尾迹长度有限。形成的涡旋固定在圆柱体后面。对于雷诺数约为 40 - 70,尾迹会有周期性摆动。对于更高的雷诺数,涡旋会从圆柱体上脱落。随着雷诺数的增加,涡旋会不断从圆柱体上脱落并被冲走。形成了两排涡旋,称为涡流街Vortex Street。现在压力阻力几乎贡献了总阻力的 90%。CD 的值在雷诺数约为 2000 时达到最低值,约为 0.9

进一步增加雷诺数会导致大角度速度和涡旋退化成湍流。

图 9:不同雷诺数下圆柱体周围的流动

图 10:不同雷诺数下圆柱体周围的流动,续。

在雷诺数范围 104105 中 在圆柱体的垂直中心线上存在层流边界层(图 10)。流动在这一点 S 分离,与圆柱体中心大约成 80角。下游可以看到宽阔的尾迹。分离区域的压力几乎恒定。图 31 中的潜在流动部分显示了观察到的 C P 分布。净压力差 P A - P B 有助于压力阻力。

当雷诺数在 2x10 5 左右时发生了戏剧性的变化 当边界层在分离前变成湍流。现在分离被推迟了 因为湍流边界层比层流能够持续更长的时间。分离点 S 现在在 130 0 处找到,如图 10 所示。尾迹现在已经变窄了。C P 分布表明尾迹中的压力现在比层流情况更高。结果是现在 C D 减少到大约 0.3

在高尔夫球中利用了绕圆柱体阻力的减少。提供高尔夫球上的凹痕的目的是使湍流发生,以减少阻力。板球投手,特别是那些投“摆动”球的人,希望球的一侧比其他一侧更光滑。想法是保持球的一侧流动层流,另一侧湍流。球将从层流向湍流侧摆动。

通过延迟或停止流动分离,可以减少压力阻力。开发的一种策略是使物体流线化。翼型表面是一个很好的例子。在自然界中,鸟翼和鱼是减少阻力物体的例子。

图 11 给出了一些熟悉的二维形状的 C D 数值。很明显,C D 取决于物体对流动的方向。图 12 给出了一些三维物体的 CD 值。

图 11:熟悉的二维物体的 CD 值。

图 12:熟悉的三维物体的 CD 值。





来源:CFD饭圈
FluentCFX燃烧湍流汽车Polyflow理论ParaViewParticleWorks管道
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-09-08
最近编辑:8天前
CFD饭圈
硕士 分享CFD文章,感谢关注
获赞 9粉丝 4文章 345课程 0
点赞
收藏
作者推荐

Fluent仿真实例-载人舱以超高音速再入大气层的仿真1网格剖分

本教程的目的是计算高超音速速度下绕过再入舱的流动。模拟的速度、轨迹和环境条件代表了这样一种飞行器在大约50 [km]高度通过地球大气层时的情况。本教程演示了如何执行以下操作:使用无缝隙几何工作流创建网格。使用高速数值、双温度模型进行能量建模,以及适用于空气属性的适当模型来模拟高超音速流动。设置外部空气动力学的边界条件。使用k-ω SST湍流模型。使用基于密度的耦合求解器计算解决方案。1. 问题描述 问题考虑了在攻角α=-25°和自由流马赫数17.0下绕过再入舱的流动。几何形状如下图所示,该图还显示了给定情况下的升力和阻力方向。对于本教程,可以假定绕过再入舱的流动是对称的。2.启动Fluent 在文章末尾下载reentry_capsule.zip文件(腾讯微云)。将reentry_capsule.zip解压到您的工作目录。可以在文件夹中找到SpaceClaim CAD文件CapsuleFlow.scdoc。使用Fluent启动器启动ANSYS Fluent。在左上角的选择列表中选择Meshing以启动Fluent的网格模式。在选项下启用双精度。在并行(本地机器)下将网格过程和求解器过程设置为4。3.启动网格工作流 a. 在工作流选项卡中,选择水密几何工作流。b.工作流的面板如下。4.导入CAD几何 a. 选择导入几何任务。b. 对于单位,选择米。c. 对于文件名,输入要导入的CAD几何的路径和文件名(CapsuleFlow.scdoc)。d. 点击导入几何。这将更新任务,在图形窗口中显示几何,并允许您继续进行工作流中的下一个任务。5.添加局部尺寸 使用面尺寸在墙表面添加局部网格尺寸控制,并在再入舱周围的流动体积中使用影响体(BOIs),在流动感兴趣的区域。请注意,所有网格尺寸,局部和全局,都比工业用途典型尺寸更粗糙,以确保网格对于教程目的不会太大。a. 在添加局部尺寸任务中,通过选择是来向面几何添加局部尺寸控制:i. 在添加局部尺寸提示时选择是。ii. 为尺寸控制的名称输入capsule。iii. 指定增长速率为1.1。iv. 保留面尺寸作为尺寸控制类型。v. 指定目标网格尺寸为0.1。vi. 选择区域作为选择方式。vii. 选择origin-capsule。b. 点击添加Add Local Sizing。i. 选择Body Of Influence作为Size Control Type。ii. 保留尺寸控制名称的默认值 "boi_1"。iii. 指定增长速率为1.1。iv. 指定目标网格尺寸为0.2。v. 选择“按标签选择”作为选择方式。vi. 选择 "capsuleflow-boi1"。c.点击Add Local Sizingi. 保留尺寸控制名称的默认值 "boi_2"。ii. 指定增长速率为1.1。iii. 指定目标网格尺寸为0.05。iv. 选择 "capsuleflow-boi2"。v. 点击Add Local Sizing以完成此任务,并继续进行工作流中的下一个任务。6.生成表面网格 按照上述设置的局部尺寸,全局表面网格尺寸仅定义了再入舱远离其他表面、流入、流出和对称表面上的最大元素。a. 在Generate the Surface Mesh任务中,您可以为面几何设置各种表面网格属性。b. 指定最小尺寸为0.1。c. 指定最大尺寸为1。d. 指定增长速率为1.1。e. 点击Generate the Surface Mesh以完成此任务并继续进行工作流中的下一个任务。7.描述几何 当您选择Describe Geometry任务时,系统会提示与导入几何的性质相关的问题。a. 在Geometry Type选择“几何仅由流体区域组成,没有空隙The geometry consists of only fluid regions with no voids”选项,因为此模型仅包含流体区域。b. 保留此任务的其余默认设置。c. 点击Describe Geometry以完成此任务并继续进行工作流中的下一个任务。8.确认并更新边界 a. 选择Update Boundaries任务,您可以在此检查网格边界,并确认并相应更改任何指定的边界。ANSYS Fluent尝试自动确定正确的边界排列。b. 为inflow边界选择pressure-far-field。c. 为outflow边界选择pressure-outlet。d. 为sym边界选择symmetry。e. 所有建议的边界都是正确的,点击Update Boundaries并继续进行下一个任务。9.更新区域 a. 选择Update Regions任务,您可以在此查看并更改从导入的几何生成的各种区域的名称和类型,并根据需要进行更改。我们可以看到,唯一定义的区域是流体区域。b. 建议的区域类型是正确的,因此点击Update Regions以更新您的设置。10.添加边界层 a. 选择Add Boundary Layers任务,您可以在此设置边界层网格的属性。b. 选择last-ratio作为Offset Method Type类型。c. 保留尺寸控制的默认名称last_ratio_1。d. 指定层数为30。为了模拟接近壁面的流动,需要许多边界层。e. 指定过渡比率为1。f. 指定第一高度为0.001。g. 点击Add Boundary Layers。11.生成体积网格 按照上述设置的局部尺寸,包括影响体,全局体积网格尺寸仅定义了流场中的最大元素。在这种情况下,最大值设置为与指定的全局表面网格尺寸一致。a. 选择Generate the Volume Mesh任务,以设置体积网格的属性。b. 选择多面体polyhedra以填充。c. 指定最大单元长度为1。d. 保留启用并行网格的默认选择。e. 点击Generate the Volume Mesh。ANSYS Fluent将应用您的设置,并继续为机翼几何生成体积网格。网格显示在图形窗口中,自动插入一个剪切平面,绘制一层单元,以便您可以快速查看体积网格的细节。12.检查和保存网格 Mesh → Check。保存网格文件(CapsuleFlow.msh.gz)。File → Write → Mesh...13.切换到求解Solution模式 现在,使用ANSYS Fluent在网格模式下生成了网格,您现在可以切换到求解器模式以完成仿真的设置。具体的仿真设置请看下一篇教程文章。来源:CFD饭圈

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈