OpenFOAM——高空腔内的湍流自然对流

本算例来自《ANSYS Fluid Dynamics Verification Manual》中的VMFL052: Turbulent Natural Convection Inside a Tall Cavity

顶部和底部壁面为绝热边界，左侧壁面（COLD_WALL）为288.25K，右侧壁面为（HOT_WALL）为 307.85K

流体的物性参数为：

密度：1.17kg/m³

粘度系数：1.81×10^-5N·s/m²

运动粘度：1.54×10^-5 m²/s

定压比热容：1005J/(kg·K)

导热系数：0.0261W/(m·K)

普朗特数：0.705

针对自然对流的方法一般有三种：

①布辛尼斯克(Boussinesq)近似

②不可压理想气体

③密度线性差分方法

本算例采用布辛尼斯克(Boussinesq)近似，布辛尼斯克假设中密度作为常数考虑，浮力作为源项添加，但仅适用于小温差,密度变化小于20%的情况,一般壁面与流体的温差超过200K，布辛尼斯克近似就不适用了。

首先进行建模操作，任何建模软件均可，本算例采用ICEM直接建模，生成网格，然后利用OpenFOAM下转化网格，划分完成的网格局部如下：

网格比较密集，高度方向设置了2000个节点，宽度方向上设置了100个节点，在导出网格的同时，我们对网格进行了缩放。建议在导出网格的时候就缩放网格，避免在后面的操作中，由于忘记缩放网格，导致计算出错

接下来转入OpenFOAM的操作：

首先新建一个文件夹，名字任取，用来作为算例文件夹，本算例中我将该文件夹命名为：Tall_Cavity

然后进入OpenFOAM的安装目录，将安装目录下的hotRoom算例（我的目录为/opt/openfoam5/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom）下的0文件夹、constant文件夹和system文件夹拷贝到Tall_Cavity文件夹下，然后删除system目录下的blockMeshDict文件，因为我们利用OpenFOAM的命令转化.msh文件为OpenFOAM能接受的网格文件。

然后我们将刚才我们生成的.msh网格拷贝到Tall_Cavity文件夹下。在算例文件夹下打开终端，输入fluentMeshToFoam命令：

我们打开constant文件夹下的transportproperties文件，内容修改如下：

说明一下：

nu表示层流粘度，这里根据上面的物性填写1.81e-05即可，注意下量纲，这里是运动粘度，需要用动力粘度除以密度

beta表示热膨胀系数，一般取为三百分之一（0.003）

Tref表示参考温度，这里我们取为两壁面温度的算术平均值，这里填写298.05

Pr表示普朗特数，这里我们根据上面的物性参数填写0.705

最后的湍流普朗特数根据实验得到，一般不用改动

接下来，修改turbulenceProperties文件的内容如下：

本算例我们采用k-ωSST湍流模型来进行计算（后面需要在0文件夹下添加omega文件）

g文件无需修改，默认设置超Y轴方向重力，重力加速度为9.81m/s²

然后对初始边界条件进行设置，下面转入0文件夹下进行操作：

删除不需要的文件

alpha文件当中的内容如下：

p文件当中的内容如下：

p_rgh文件当中内容如下：

U文件当中的内容如下：

将T.org的文件名修改为T,当中的内容如下：

k文件当中的内容如下：

说明一下：

在壁面处，k应该为0。这里我们使用了壁面函数。

k表示湍动能，计算公式为：

其中：u_avg为平均速度

I为湍流强度

湍流强度的计算公式为：

式中：Re表示雷诺数

为何一般初始湍流强度设置为5%，可参考Launder B E, Spalding D B. The numerical computation of turbulent flow Computer methods in applied mechanics and engineering[J]. Надійшла до редакції 04.12. 2009 р, 1974.

MLA(文章链接：http://www.sciencedirect.com/science/article/pii/0045782574900292)

omega文件当中的内容如下：

说明一下：

omega表示比耗散率，计算公式如下：

式中：C_μ为0.09

k为湍动能

l为特征长度

nut文件内容如下：

说明一下：

nut表示湍流粘度，其计算公式为：

接着我们设置controlDict文件

fvSchemes文件修改如下：

fvSolution文件修改如下：

为了加快计算速度，我们进行并行计算，回到算例文件夹下的system目录，新建一个decomposeParDict字典，其中的内容：

打开终端，输入decomposePar进行分块

由于我安装了PyFoam来实时输出残差，所以在终端中输入pyFoamPlotRunner.py --clear mpirun -np 4 buoyantBoussinesqSimpleFoam -parallel开始计算：

等到计算结束

将计算结果导入paraview里面进行处理

计算结果与实验数据对比：

转载于:https://www.cnblogs.com/liusuanyatong/p/11259648.html

OpenFOAM——高空腔内的湍流自然对流相关推荐

python反射和高阶内置方法
1.isinstance:判断对象和类的关系 #判断结果返回bool类型 class A:pass class B(A):pass a = A() print(isinstance(a,A)) #Tr ...
Python学习（十）——高阶内置函数
内置函数内置函数高阶内置函数 abs().sum().round() map() filter() sorted() max().min() 内置函数内置函数:python提供的,可以直接使用 ...
积分球高漫反射内表面涂层
积分球-光谱仪系统一般由光谱仪.球体.光纤.光色度探头.软件等组成.积分球是一个内壁涂有高漫反射材料的空腔球体,又称光度球,光通球等.球壁上开一个或几个窗孔,用作进光孔和放置光接收器件的接收孔.球内壁 ...
OpenFOAM——过渡管中的湍流
本算例来自<ANSYS Fluid Dynamics Verification Manual>中的VMFL016:Turbulent Flow in a Transition Duct 一 ...
国外长高书籍内有详细锻炼方法
为什么80%的码农都做不了架构师?>>> 这个章节中介绍的运动本身目的就是释放这些压力,从而你的身体可以继续长高并且达到潜在的最大高度.有一个观点是我们普遍接受的,那就是物理运 ...
我们可观测的宇宙或许仅是更高维度内的一个泡沫
新理论提出,我们的宇宙也许位于一个正在向更高维度扩张的泡沫之中. 新浪科技讯北京时间2月22日消息,据国外媒体报道,就像海浪最高点的泡沫一样,我们可观测的宇宙或许仅仅是一个不断向更高维度扩张的泡沫内 ...
高云片内时钟OSC的调用
一.高云FPGA器件内部内嵌了一个硬件晶振OSC,可以为用户提供时钟源. 原本时钟频率:根据不同器件OSC晶振频率取值不同,GW1N-4,GW1NR-4,GW1N-4B,GW1NR-4B,GW1NRF ...
[html] 可以给内联元素设置宽和高吗？为什么？
[html] 可以给内联元素设置宽和高吗?为什么? 1.正常情况下不能设置宽高(可替换内容元素除外例如img)2.元素被浮动,绝对定位,固定定位后,可以设置宽高3.display:block,disp ...
【基线,内容区, 行高/行间距,行距, 行内框,行框的区分说明】
基线,内容区, 行高,行距, 行内框,行框的区分说明 ⑴ 基线 ⑵ 内容区 ⑶ 行高/行间距 ⑷ 行距 ⑸ 行内框 ⑹ 行框 ♣ 结束语和友情链接 ⑴ 基线基线沿着文本最底部的一行线, ...

OpenFOAM——高空腔内的湍流自然对流

OpenFOAM——高空腔内的湍流自然对流相关推荐

最新文章

热门文章