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

2023-07-11
关注

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:

  • 高仪
您觉得本篇内容如何
评分

相关产品

成都迈铄智感 MSD-300 高精度数字波高仪|浪高仪 波高仪

MSD-300 数字波高仪是一款高精度和高可靠性的波高测量仪器,它与TIA/EIA-485数字总线接口相结合,消除了模拟数据采集系统中常见的许多问题,例如长电缆上的噪声和接地问题。这使得MSD-300非常适合用于高质量的波高测量和主动波吸收系统的集成。

GMC-Instruments 高美测仪 M248A 频率计

"The METRAHIT ULTRA M248A is a precision high resolution TRMS system digital multimeter with IR remote control capability, temperature measurement and data logging for demanding, universal applications in the laboratory or out in the field on service calls.Instruments included in the Gossen high-resolution series are extremely rugged, reliable digital multimeters with housings made of impact resistant ABS. With a resolution of 310,000 digits and roughly 30 different measuring functions, they’ve been developed for professional use.Triple display with a resolution of ±310,000 digits,Intrinsic error: ± (0.02% rdg. + 20 digits),Measuring functions: 1 µV to 600 VDC and 10 µV to 600 VAC\/AC+DC TRMS (bandwidth: 15 Hz to 100 kHz), 1 nA to 10 ADC and 10 nA to 10 AAC\/AC+DC TRMS (bandwidth: 15 Hz to 10 kHz),Resistance, frequency, capacitance and temperature measurement: 1 mΩ to 30 MΩ, 0.001 Hz to 1 MHz, 1 pF to 3 mF, °C\/°F Pt100, Pt1000, type K thermocouple: -270 to +1372 ºC,Automatic range selection and battery cutoff,Min-Max measured value storage,Integrated 2 MByte measurement data memory for approx. 300,000 measured values, sampling rate: 100 ms to 9 hrs.,IR data interface (this instrument can be controlled via the IR interface with optional accessory Z216C - RS Stock No 606-3776),CAT IV @ 300 V or CAT III @ 600 V per IEC 61010-1,DAkkS calibration certificate as standard featureMultimeter Type = HandheldAbsolute Maximum AC Current Measurement = 10A acAbsolute Maximum DC Current Measurement = 10A dcAbsolute Maximum AC Voltage Measurement = 600V acAbsolute Maximum DC Voltage Measurement = 600V dcAbsolute Maximum Resistance Measurement = 30MOhmsAbsolute Maximum Capacitance Measurement = 3mFAbsolute Maximum Temperature Measurement = +100 (PT100) °C, +1372 (Type K) °C, +80 (Internal) °C, +850 (PT1000) °CFunctions Measured = AC Current, AC Voltage, Capacitance, DC Current, DC Voltage, Diode, Frequency, Resistance, TemperatureTrue RMS = Yes"

鹏城半导体 高真空磁控溅射仪 控制器及系统

鹏城半导体 高真空磁控溅射仪

Taiken 泰肯光电 TC-AMBox系列高稳定电光强度调制仪 其它

TC-AMBox系列电光强度调制仪是TC拥有自主知识产权的高度集成化产品,该仪器将电光强度调制器、微波放大器及其驱动电路等必要集成于一体,这不但方便了用户的使用,而且大大提高了MZ强度调制器的可靠性,并可根据用户要求提供定制化服务

Cubic 四方光电 Gasboard-3000UV 紫外NOx排放分析仪

Gasboard-3000UV紫外NOx排放分析仪是基于紫外吸收光谱气体分析技术,自主研发的新一代氮氧化物分析仪。采用独特算法,高精度直通短气室,抗干扰能力强,测量精度高,满足GB 3847和GB 18285标准00级测量精度要求。

深圳深国安 SGA-606-C4H8O2 气体仪器

高灵敏便携式丁酸检测仪-5秒响应结果

汉威科技 WP160 燃气调压箱柜压力在线监测

WP160型无线智能压力监测仪为我公司生产的具备无线上传功能、用于监测表压压力的一款先进仪器。本监测仪采用高性…

Laoying Haina 崂应海纳 LY-NDUV-NO 气体仪器

本产品专为NDUV氮氧化物分析仪的NO气体分析开发,可用于烟气分析、汽车尾气分析、船舶尾气分析等,体积小巧、测量精度高,方便集成于各种NO气体分析仪。

评论

您需要登录才可以回复|注册

提交评论

感知论坛

这家伙很懒,什么描述也没留下

关注

点击进入下一篇

新基建新动力,光储充一体化电站怎么赚钱?

提取码
复制提取码
点击跳转至百度网盘