基于Fluent二自由度圆柱涡激振动数值模拟与应用研究
- 作者优秀
- 优秀教师/意见领袖/博士学历/特邀专家/独家讲师
- 平台推荐
- 内容稀缺
导读:大家好,我是仿真秀专栏作者-姜工,流体开发工程师,研究生阶段主要研究海洋管柱涡激振动抑制装置的优化问题,参与多项项目,申请发明专利数项,发表涡激振动方面学术论文SCI两篇,中文核心一篇;工作后从事化工流体机械方面的的设计开发与仿真优化方面的工作,对旋转机械、多相流方面有丰富的仿真经验。近日,我计划在仿真秀官网和APP原创精品课《基于Fluent的二自由度圆柱涡激振动数值模拟与应用研究8讲》,欢迎专家学者批评指正。以下是正文:
一、案例描述
本案例源自美国康奈尔大学经典圆柱水槽试验,在实际的工程应用中,管柱被认为是一个有弹性支撑的刚性结构,结合切片理论可将计算模型简化至二维,可以看作是一个质量弹簧阻尼系统,如图中所示二、采用Designmodeler进行几何模型建立
流体计算域模型如图所示,模型使用的计算域为30D×20D矩形区域。入口边界距离圆柱中心10D处,出口边界距离圆柱中心20D处,上部边界和下部位于距圆柱体中心10D处,并被定义为对称边界。其中圆柱位于入口下游10D处,位于横流方向的中心线(y=0)。阻塞率为圆柱直径与流场域宽度的比值,本文阻塞率为0.05,不超过0.05,故满足要求。采用“刚性运动区域+动网格区域+静止网格区域”多区域划分的方法来建立圆截面涡激振动流场区域,整体流场分为外流域部分和内流域部分。注:为保证流场在内外流域交界地方一致性,选中内外流域,鼠标右击from a part,可保证后续画网格时具有网格一致性。
三、采用Meshing进行网格划分
Workbench中导入Fluid Flow
- 双击Meshing,详细网格划分步骤在这里不介绍了
整体流场在进行网格划分时,采用“刚性运动区域+动网格区域+静止网格区域”的 方法来建立圆截面涡激振动网格划分区域,如图中所示,具体划分区域如下:(1) 动网格区域采用非结构化三角形网格进行网格划分,Fluent动态网格模型算法采用Smoothing+Remeshing的方式更新动网格,此区域为变形区,为圆柱的振动提供了一个运动的空间,其网格会因圆柱的运动而进行重新划分,因此网格的初始质量好但随着重新划分可能有所下降,但可通过参数设置使重新划分后的网格质量尽可能好。此区域的大小由圆柱的振幅大小决定, 考虑到文献实验中的最大振幅为1.5D左右,因此动网格区域尺寸取为6D×6D,留有充足的空间允许圆柱进行大范围的的振动,能够很好地满足数值计算的需要(避免振幅超过容许空间从而挤压外边界形成负网格);(2)刚性运动区域主要包括圆柱流固耦合边界与边界层区域,圆柱近壁面边界层网格采用膨胀层网格处理,使在主要流速下预估的Y+=1,边界层第一层高度可按经验公式计算,近壁面划分15层四边形结构化边界层网格,网格增加率设置为1.2,向外过渡至与动网格的交界处。边界层区域网格尺寸较小、网格划分质量较好并具有与圆截面相同的运动速度,因此将随着圆截面的振动而振动,由此来保证圆截面附近的网格一直处于初始的高质量状态从而确保各时间步都具有足够的计算精度。 (3) 静止网格区域采用四边形结构化网格,整体网格尺寸网格尺寸保持不变,用以为模拟尾流中的涡街提供更好的计算精度。四、Fluent 基本设置
点击左侧模型树节点 General ,出现 General 面板,如图所示。
点击scale…
将单位修改为与几何建模时尺寸单位保持一致
点击 Check 按钮,屏幕输出网格信息如图所示,确保网格最小体积为正值由于涡激振动属于周期性振动,Time选项设置为transient,进行瞬态研究选择模型树节点 Models ,在右侧 Models 列表项中双击,弹出湍流模型设置对话框。设置参数:选择模型树Materials,分为流体材料fluid和固体材料solid,双击材料空气air进行设置,弹出材料设置窗口
点击Fluent Datebase…
下拉菜单至最后,找到材料水(water)
点击Copy
选择模型树Cell Zone Conditions,双击弹出区域设置窗口,如图所示- 分别选择part-dongquyu和solid-part-jignzhiyu,将属性设置为流体
- 将material name 选择为 water-liquid
其余保持默认,点击close关闭
步骤六:Boundary Conditions
点击模型树节点 Boundary Conditions ,参数面板如图所示。
1、Inlet 边界设置
鼠标双击 inlet 列表项。设置参数:
- 设置 Velocity Magnitude 为 0.09m/s
- 设置 Specification Method 为 Intensity and Hydraulic Diameter
- 设置 Turbulent Intensity 为5%
- 设置Hydraulic Diameter 为 0.0381
选择OK按钮关闭对话框
2 、 Outlet 设置
鼠标选中 outlet 项,设置 Type 为 Pressure-outlet,Gauge pressure设置为0,其他参数保持默认。
步骤七 :Dynamic Mesh
点击模型树节点 Dynamic Mesh ,参数面板如图所示。
导入外部UDF文件
- 点击Functions,选择第二个Compiled
点击Build,点击加载Load
边界层网格分离
- 设置number of cells 为相应边界层网格的层数
点击close关闭
- 点击Domain 中的 separate,选择cell…
- 选择boundary_0和part-dongquyu
点击close关闭
边界运动或变形的指定
- 点击动网格部分的dynamic mesh zones的create
- 将分离出的边界层网格part-dongquyu:013进行相同的操作
其余默认,点击close关闭
步骤八 Reference Values
选择模型树节点 Reference Values,如图所示
步骤九 Solution Methods
此模型保持默认参数
Step 11 :Solution Controls
保持默认参数
步骤十 Monitors
选择模型树节点 Monitors,如图所示,选择Report Files,鼠标右键,点击edit‘’
- 勾选report file和report plot
点击OK,关闭窗口即可
步骤十一 Initialization
点击模型树中的 initialization,出现如图所示窗口
- 选择 Standard Initialization
点击 initialize
步骤十二 Calculation activities
双击 solution animations,如图所示
- 云图Contours of 选择 为velocity
- 鼠标滚轮调整到合适窗口,选择contour2,点击use active
点击OK
步骤十三 Run calculation
点击左侧模型树节点Run calculation ,出现 Run calculation 面板,如图所示。
- 将number of time steps 设置为需要时间步数
- 每一个时间步时长 time step size设置为0.01
点击calculate,进行计算
步骤十四 查看计算结果
点击模型树Results,可查看残差曲线、各种云图(速度、压力等),速度矢量图、动画等本案例结束,更多内容欢迎大家订阅我的精品课。
五、我的精品课
《Fluent涡激振动数值模拟与应用研究8讲》全过程仿真复现经典管柱涡激振动试验。本课程通过对国际顶尖高等学府康奈尔大学研究人员的经典二自由度管柱涡激振动试验过程的深入了解及复现,快速掌握ANSYS workbench的仿真方法及建模技巧,熟练开展包括几何模型建模、网格划分、fluent参数设置、后处理等全过程仿真。以下是课程安排
此外,笔者为所有订阅用户提供知识圈答疑服务和VIP用户交流群。并附赠课程相关资料等(平台支持自行开具电子发票)。
3、读者福利(扫码立即领取)
声明:本文首发仿真秀App,部分图片和内容转自网络,如有不当请联系我们,欢迎分享,禁止私自转载,转载请联系我们。
获赞 10000粉丝 21452文章 3502课程 218