python做圆柱绕流_基于snappyHexMesh生成网格的圆柱绕流算例

本文档详细介绍了如何在OpenFOAM中使用snappyHexMesh工具生成圆柱绕流的网格。首先,通过3Dbuilder或3DCAD创建STL几何文件,接着创建背景网格,使用surfaceAutoPatch划分多个表面。然后复制并修改snappyHexMeshDict等配置文件,最后执行表面特征提取、blockMesh、snappyHexMesh等步骤生成网格。对比snappyHexMesh与blockMesh生成的网格效果,发现snappyHexMesh在壁面附近的网格更加精细。
摘要由CSDN通过智能技术生成

基于snappyHexMesh生成网格的圆柱绕流算例

1.创建几何体STL文件

在了解了原理之后,要生成网格第一步就是要创建一个几何体的STL文件,可以使用3Dbulider,3DCAD等软件生成几何体然后以STL格式输出,然后将文件拷贝到虚拟机中。可以使用paraview打开几何体进行查看,确认没有问题。

2.创建背景网格

修改blockMeshDict字典文件来创建一个背景网格,在这个算例中背景网格是一个非常简单的block,背景网格的几何大小与第二章的圆柱绕流算例的尺寸一致。

需要注意的是,在blockMeshDict中所有面的type都为patch。 生成的背景网格如下:

将几何体放置到网格中,发现没有问题,继续之后的步骤。

3.为几何体创建多个表面

默认情况下,几何体STL文件由一个表面组成,如果要创建多个表面,可以使用surfaceAutoPatch。在几何体STL文件所在的目录下终端键入:

surfaceAutoPatch geo.stl cylinder.stl 130程序surfaceAutoPatch将读取原始STL文件geo.stl(一个面组成),并将使用130的角度标准将几何体划分为多个面。

这个角度划分标准类似于表面特征细化角,如果几何体相邻面角度大于resolveFeatureAngle,则标记相邻STL面进行划分(细化),反之不进行划分(细化)。

cylinder.stl最终要拷贝到constant/triSurface中。

这时候打开cylinder.stl,会发现里面定义了三个表面:patch0(通过部分代码上的点位置可以确定patch0是圆柱体的侧表面)

patch1(可以判断为后端面)

patch2(同理可以判断为前端面)

在paraview中可以直观的看到三个面分别是不同的颜色:

4.复制snappyHexMeshDict,meshQualityDict和surfaceFeatureExtractDict文件

在完成上面的操作后,就可以进行网格细化了,网格细化需要snappyHexMeshDict,meshQualityDict和surfaceFeatureExtractDict这三个文件ÿ

这是一个关于使用Python圆柱绕流的Fluent学习笔记。卡门涡街是在定常不可压缩势流圆柱绕流中观察到的一种稳定的涡街现象。在这个问题中,我们可以使用Fluent软件模拟流体流动,并使用Python脚本来自动化模拟和数据分析。 以下是一个Python脚本示例,用于设置和运行Fluent模拟,并将结果保存到文件中: ```python import os # 设置Fluent环境变量 os.environ['FLUENT_HOSTNAME'] = 'localhost' os.environ['FLUENT_ARCH'] = 'lnamd64' os.environ['FLUENT_INC'] = '/usr/local/Fluent.Inc/fluent/fluent17.2.0' os.environ['FLUENT_LIB'] = '/usr/local/Fluent.Inc/fluent/fluent17.2.0/lib' # 导入Fluent模块 from fluent import * # 创建Fluent进程 fluent = Fluent() # 打开Mesh文件 fluent.meshRead('cylinder.msh') # 设置模拟参数 fluent.timeScheme('steady-state') fluent.physics('viscous') fluent.model('inviscid') fluent.material('air') fluent.viscosity(1.789e-5) fluent.density(1.225) fluent.referenceValues(1.0,1.0,1.0,1.0) # 定义边界条件 fluent.boundaryCondition('inlet', 'velocity-inlet', u=0.2, v=0.0, w=0.0) fluent.boundaryCondition('outlet', 'pressure-outlet', p=0.0) fluent.boundaryCondition('cylinder', 'wall') # 运行模拟 fluent.solve() # 保存结果 fluent.writeCase('cylinder.cas') fluent.writeData('cylinder.dat') ``` 在这个示例中,我们首先设置了Fluent环境变量,并导入了Fluent模块。然后,我们创建了一个Fluent进程,并打开了一个Mesh文件。接下来,我们设置了模拟参数、定义了边界条件,并运行了模拟。最后,我们保存了结果到文件中。 使用Python进行数据分析也非常方便。我们可以使用pandas库来读取和处理数据文件,并使用matplotlib库来绘制图形。以下是一个简单的Python脚本示例,用于读取和绘制圆柱绕流结果数据: ```python import pandas as pd import matplotlib.pyplot as plt # 读取数据文件 data = pd.read_csv('cylinder.dat', delim_whitespace=True, header=None, names=['x', 'y', 'z', 'u', 'v', 'w', 'p', 'c']) # 绘制速度云图 plt.contourf(data['x'].values.reshape(-1, 201), data['y'].values.reshape(-1, 201), data['u'].values.reshape(-1, 201)) plt.colorbar() plt.title('Velocity Contour') plt.xlabel('x') plt.ylabel('y') plt.show() # 绘制压力分布图 plt.plot(data['y'][::201], data['p'][::201]) plt.title('Pressure Distribution') plt.xlabel('y') plt.ylabel('p') plt.show() ``` 在这个示例中,我们使用pandas库读取了数据文件,并使用matplotlib库绘制了速度云图和压力分布图。这些图形可以帮助我们更好地理解流体流动的特性。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值