VOF方法的原理是什么？如何设置浪高仪？

NO.1

数据处理-各种平均值的含义

四种平均：facet，area，mass，vertex应该用哪个？

在Report definitions定义时，可以看到surface report中有四种平均，都是什么意思呢？又该选择哪个？

1，Facet 是算数平均，计算的是面上的值求和然后除以面的个数

2，Area 是面积加权平均，就是每一项×面积相加后再除以总面积

3，Mass 是质量加权平均，

4，Vertex 是节点平均，顶点平均的解释指定字段变量在曲面上的顶点平均，是通过所选变量的顶点值之和除以顶点总数来计算的，计算的是面上的点的值求和然后除以点的个数，我理解的就是计算每个网格中心点的值相加以后再除以点的总数，浪高仪的设置就是用的这个平均。

每个平均得到的结果也有可能不同，看自己研究对象选择哪个吧。

1

NO.2

VOF方法的原理

在 ANSYS Fluent 中采用了两相流模型 体积分数法 （VOF）来追踪空气与水的自由界面，该方法对两相流体，各自定义了在每一个计算单元中所占的体积分数函数。

什么意思呢，一开始我也不太懂，现在我觉得是这个意思

比如我们看见的现实的波浪，它的波动都是在水面上，也就是水和空气的交界面上，越往水下，其实波浪就不是太明显了。所以把波浪看作表面波的话，那么波面的vof就等于0.5，即一半水，一半空气的比例，后面设置 波面检测 ，就是这个原理。

VOF=Volume fraction ，即体积分数（我理解的）

2

NO.3

如何设置浪高仪

前面铺垫的两个，就是为了理解浪高仪为什么要这样设置。网上应该有教程，但我一开始不理解为什么要这样。

假设坐标系如上图，下面来介绍一下如何设置，以及为何要这样设置浪高仪。因为查看的是水面上某一点随着时间流动，它的波面的变化。所以最终要建立一个点。而建立点的顺序为面-线-点

1、设置等值面ISO-surface

为什么不能设置一个平面或者直接设置一个点呢，因为plane和point，是不会随着波浪移动的，它就是固定的一个点，一个面，那就没办法监测它的波高。（我理解的）

点击surface->create->ISO-surce，弹出以下界面，起个名字，然后选择phase(相)，相选择空气和水都没有影响，将ISO-values值设置为0.5，即水气交界面，建立的就是一个水和空气体积分数始终都为0.5的等值面。基于模型建立，->create

2,建立线line--基于mesh

建立vof面后右边就会出现vof面,在这个面上找一条线,基于mesh,基于z方向,就会建立一条线,可以compute以下,会出现z的最大值和最小值坐标,我一般选择面上中心的那条线,就是z=0.4处了.

3,建立点-基于mesh

有了线,基于这条线,只要指定x的坐标,就能创建沿着x向不同的点了.

为什么要基于mesh ,一开始我不理解,mesh 就是网格的意思,网格又是怎么存储的呢,以坐标的形式,我觉得,每个网格中点的坐标计算机其实是已知的,就可以根据xyz坐标,建立各种数据监测点.

建立万可以display一下,看一下是不是正确的,这个图上就明显的有面-线-点了.

下面就是定义数据监测了Report Definitions->New->Surface Report->Vertex Average,打开下个界面

基于mesh ,Y-coodinate,因为监测的波高,波高方向是y,所以要看要研究点的波高,即Y方向的变化的值.->ok就可以了,后面计算的时候,窗口会显示这些数据的变化的.

这些数据都是.out文件,可以将后缀改成.txt文件,打开excel,打开->选择该.txt文件,选择分隔符号,下一步,再选择Tab键和空格,下一步,下一步,完成.数据就导入进去了,删除不用的文字就可以处理了.

3

**NO.4 **

参考压力,总压,表压

--这个可以看B站:计算传热学大叔,讲的特别好,还通俗易懂

表压 ，也叫相对压力，fluent中存储的是表压，

参考压力 ，就是fluent中的操作压力。

绝对压力 ，如果是不可压缩流动，不用管这个参数，当计算理想流体时，则必须要用到绝对压力，它等于相对压力＋参考压力。

对不可压流动，绝对压力没有意义，压差才有意义。

对于不涉及任何压力边界条件的不可压缩流动,ANSYS FLUENT在每次迭代后要调整表压值。这个过程通过使用参考压力位置处(或该位置附近)节点的压力完成。因此,参考压力位置处的表压应一直为0。如果使用了压力边界条件,则不会使用到上述关系,因此参考压力位置不被使用。

参考压力位置默认为等于或接近(0,0,0)的节点中心位置。实际计算中可能需要设置参考压力位置到绝对静压已知的位置处。在Operating Conditions对话框中的Reference Pressure Location选项组中设置新的参考压力位置的x,y,z的坐标即可。

如果要考虑某一方向的加速度,如重力,可以勾选Gravity复选框。

对于VOF计算,应当选择SpecifiedOperating Density,并且在OperatingDensity 下为最轻相设置密度。这样做排除了水力静压的积累,提高了round-off精度为动量平衡。同样需要打开 Implicit Body Force,部分平衡压力梯度和动量方程中体积力,提高解的收敛性。

ReferencePressure Location(参考压强位置)应是位于流体永远是100%的某一相(空气)的区域,光滑和快速收敛是其基本条件。

参考压力是为了减少舍入误差计算误差设置的 。可以把流体域内某一点设置成参考压力点，该点的压力为0，这个参考点的压力在每一步迭代后都会变成0，如果没有这个锚点，每次计算迭代一步会飘。每一次迭代，其他点的压力都要以该点作为参考。如果没有压力边界，就要设置压力参考点。当然也可以设置在流体域外，fluent会自动将离该点最近的流体域的点设为参考点。

1，压力出口边界有两个参数，表压默认为0，默认压力出口处的静压为0

静压和总压的区别：对于有粘性的流体来说，根据伯努利定律，总压驱动流动，总压=静压＋动压，如果左右两侧的总压，左边大于右边，那流体一定是从左向右流动的。但静压不一定是左边大于右边，因为还有动压。但对于无粘性流体，没有粘性损失，左右的总压一定是一样的。

2，对于回流出口边界，总压就是环境压力。

4

**NO.5 **

二维mesh网格划分,三角形和四边形同时存在

我发现在画二维网格时，在边界添加尺寸后，画出来的网格只在边界上加密，而内部区域非常差，也不是均匀分布的四边形网格，我发现这是因为没有添加面网格尺寸

Constrain Boundary 选择yes就可以画出均匀网格了。

自己瞎试试,或许也能找到解决办法.

5

NO.6

什么是源项

对固体传热什么的，只用考虑能量方程就可以了，因此固体只有一个方程,能量方程.但流体就不一样了.

流体力学三大方程：

质量方程（连续方程） kg/m3.s

**动量方程 **N/m3（N.s/m3.s）

**能量方程 **W/m3（J/m3.s）

动量方程： 每立方米空间内加载的体积力

能量方程： 每立方米空间内每秒的发/吸热量。造波不涉及

质量方程： 每立方米空间内每秒产生的物质质量，一般用于传质过程。连续方程。

对于不可压缩流体，密度为常数，所以就有

源项造波包括质量源项和动量源项造波 ，原理就是在动量方程和质量方程右边添加一个源项。

源项的方程要根据自己的研究来定，我不太理解什么意思。我以为的可能就是比如质量源造波，连续方程的右端本来是0，也就是他是处于一个平衡状态的，这时候如果在右边添加一个不为零的项，无异于给了它一个“推力”。看文献说，质量源造波类似于水下爆炸，这个爆炸的能量使得水体不再静止，开始流动，从而形成波。海绵层消波的思想是在动量方程的右边添加一个源项，由于要消波，所以他的思想是希望产生和传播过来的波相对抗的能量，也就是“1-1=0”的思想。明白了这些，UDF就可以编写了，但这也是个大工程。(我理解的,可能不对)

不管是什么源项，都是针对cell来说的，所以他的定义实在cell-zone-conditions那里，而不是在boundary-conditions那里， 不同的区域添加的源项可能不同，这个可以在UDF时指定，也可以在DM画图时就已经把区域分好，但还是建议在UDF编译，这样好改，如果模型已经划分好了固定了，后面改参数就比较麻烦了。

6

NO.7

扩散项，非稳态项，源项的理解——NS方程

N_S方程根据牛顿第二定律推导而来，实际上就是动量的守恒，微原体左边的动量等于右边，N-S描述的流体运动的普遍规律，就是动量方程。只是有不同的形式。

即其对流项均采用散度形式表示的形式，这种控制方程的形式称为控制方程的守恒形式，这种方程称为守恒型的控制方程，即通用形式，分别为非稳态项，对流项，扩散项和源项。从微元体的角度考虑，守恒型控制方程等价于非守恒型控制方程

对流项的意义就是“速度运输速度本身”，或者速度运输动量。扩散项是由梯度驱动的，和扩散定律及热传导一致，作用效果在于空间上重新分配动量，从而减少速度梯度，使得速度场更光滑。

这可能就类似于,一个商场本来关着门,营业的时候,大家都冲进去了,那我们对这个商场来说就是外来物,就是"源项",人与人之间的流动是对流项,如果这时候商场有个孕妇生了个孩子,那这个孩子就是扩散项,假设这时候商场达到一个平衡,但又有人进来了(源项),商场又引起了波动,就产生了"波".(只是我的理解,可能不对).

我觉得ZhiHu上的这个回答更好理解.

7

**NO.8 **

DM中添加材料和添加冻结的区别

在DM中，没有像ABAQUS一样的装配体功能，模型在默认条件下是激活状态（解冻）

所有的接触在一起的解冻（又称激活）状态下的几何体会自动合并为一个零件，要创建多个零件，需要用到冻结功能.

在DM****中建模时

Add Material****表示添加材料，建立的几何体将合并到与它相接触的零件中；这时候即使你建立了两个零件,如果它俩接触了,那也会变成一个part.

Add Frozen****表示添加冻结体，建立的几何体将成为独立的零件，不管有没有零件与它相互接触。

8

NO.9

文件名后缀的意思

.h5文件：

目前新版的fluent默认保存为.cas.h5文件，这个文件的优势是在于可以减小文件本身所占用的大小，缺点也十分明显，不能读入旧版tecplot中，就十分的麻烦。目前一共有三种解决方法；使用新版本

（2）输入tui命令就可以很好的修改保存格式。输入代码为

f c n wcd 3

就可以非常快速的得到命名为123.cas和123.dat的文件。

动画存储时，若不在Workbench环境中，通过桌面快捷方式单独启动CFD-Post，则需要读取.dat格式结果，不能读取dat.h5格式结果

（3）点击File -> Preference... -> Default Format for I/O -> legacy" 把CFF改成legacy即可 。以后保存就默认.cas和.dat了。

CFD-POST里面导入的是.dat文件，所以这一步还是很有必要的，要不然没办法后处理。我当时算了结果却导不进POST里后处理,真的服了.

.cxa文件：点击results--animations--animation play back-read ,读取cxa文件，点击播放即可观看动画，也可用看图软件打开hmf文件

也可以在动画后处理的编辑对话框导出MPEG格式的一系列图片，然后可以用PS做成GIF动图，然后这些文件就没用了。

.HSF文件： 动画文件，开头的.cax文件相当于一个目录，计算的视频顺序都在里面。视频文件。

9

NO.10

Fluent单位制

默认为国际单位

m s kg

所以对于数据单位制、网格尺度等信息需要进行检查，保证单位制的统一。对于非国际单 位制的输入，FLUENT 表达式在运行过程中会自动且强制地转换为国际单位制。

首先，在flunet材料设置中，这里的粘度指的是流体的“动力粘度”也就是说这里的kg/m-s是流体的动力粘度，单位Pa.s等价于kg/m-s：