OpenFOAM
文章平均质量分 60
ambu1230
这个作者很懒,什么都没留下…
展开
-
解决OpenFOAM颗粒计算输出文件Paraview无法打开问题(二)
第二个方案的源是在CFD中文网上看到的一篇帖子,其具体链接忘了。这个帖子给了一个github的链接,就是将OpenFOAM输出的颗粒位置信息转变为真实的位置信息的脚本。。原创 2023-08-23 11:37:29 · 720 阅读 · 0 评论 -
集群安装OpenFOAM
在个人电脑上安装比较简单,在服务器上安装多出几个步骤,需要首先安装boost库以及openmpi库下面步骤主要参考,这里我自己重复一遍。原创 2023-08-11 14:17:15 · 339 阅读 · 2 评论 -
OpenFOAM的fvOptions
对此,我在进行槽道流仿真时,进行了验证。对比了中间截面的平均速度和出口截面的平均速度,发现它们几乎一致。这个fvOptions影响的是方程中的源项,类似重力项,通过调整体积力的大小来设定速度的大小。采用OpenFoam中的fvOptions。进行如下设置,可以大致设定流动的平均速度。原创 2023-08-19 19:24:36 · 261 阅读 · 0 评论 -
DPMFoam给粒子添加周期性边界条件
DPMFoam是OpenFOAM中的一个颗粒两相流求解器,在流体的运算中,我们经常采用周期性边界条件,在OpenFOAM中即为cyclic,那么对于颗粒,能否在边界也设置成周期性边界呢,在CFDEM中是可以的,因为liggghts的设置项中即有颗粒的周期性边界条件设置,p代表周期性边界。其中rebund为反弹边界,stick为吸附边界,escape为逃逸边界,这几个都有比较明确的意义,但是没有出现周期性边界。首先我们需要设置流体区域的上下表面为cyclic边界条件,然后设置粒子的相应边界为none。原创 2023-08-10 16:51:07 · 413 阅读 · 2 评论 -
OpenFOAM常用命令(自用)
【代码】OpenFOAM常用命令(自用)原创 2023-06-20 15:24:26 · 428 阅读 · 0 评论 -
MATLAB读取OpenFOAM的二进制文件
ny*nz是数据的大小,nx,ny,nz是每个方向的网格数,由于速度有三个分量,所以这里乘以3。忘了,不重要,知道是double就行)。前面的括弧就是只有一个字节。最后附上一个我个人写的读取数据的函数,该函数有四个输入参数,路径(fpatch),以及各个方向的网格数(gridNumY,gridNumX,gridNumZ)上面是OpenFOAM二进制文件的格式,我们可以看出,前面21行都是无关的说明文件,22开始时除了一个括号之外,其它的都是数据。读取数据的思路非常简单,忽略不需要的,读取需要的。原创 2023-06-17 22:42:47 · 1566 阅读 · 0 评论 -
如何计算努塞尔数
其实这里还有一点不是很理解,就是为什么温度梯度那一项经过平均后除以分母会等于1,如果是pure conduction的话,全场内任意处的温度梯度应该是一致的,因为温度会线性分布,但是现在是对流,全场温度梯度经过平均后会等于pure conduction的情况吗?在传热传质基础一书中,努塞尔数的定义如下,对流(convection)传热与纯热传导(pure conduction)传热之比。这里的下标V,t代表空间平均和时间平均。但是似乎在热对流研究领域,努赛尔数有另一个定义,即总传热与热传导传热之比。原创 2023-02-13 20:39:38 · 3743 阅读 · 0 评论 -
关于瑞利数等于0的讨论
一般来说,对于浮力驱动的热对流流动,其控制参数有Pr和Ra,由相似理论,我们只要控制这两个参数一致,那么经过无量纲化的流动参数也一定是一致的,也就是Ra和Pr唯一确定了流动状态。是一个比较可靠的做法,因为这样不会导致无量纲温度的值奇异,当然还有一个重要的原因是对于理想槽道流流动来说,是没有重力的,这样的话。,如果上下壁面没有温差,并且我们设定初始场为上下壁面温度,那么根据能量方程,这个场会一直的保持均一下去。可以看出,Ri为0的话只能是Ra为0。,但是其最总的结果确是不同的。,那么无量纲的温度则是奇异的。原创 2023-02-10 23:14:00 · 683 阅读 · 0 评论 -
OpenFOAM里fixedJump命令注意事项
在使用这个命令时,需要注意一点,即并行计算时需要将Jump的两个面划分在同一个分块里,不然会出现压力云图的异常情况。之前运行算例的时候,没有注意到这一点,过了好多天之后,今天终于发现了这个解决办法。以三维槽道流为例,我们选用最简单的scotch分类方法,decomposeParDict文件如下最总运行结果如下如下图所示:上面的现象,我分别采用pimpleFoam,icoFoam求解器分别运行过,都会得到上图所示的压力云图异常,后来暂时没找到解决办法,修改了求解器,采用体积力来驱动流体流动,不在使用fi原创 2023-01-26 16:54:00 · 144 阅读 · 0 评论 -
OpenFOAM里fixedJump命令
在这里,周期性边界条件的意义是进口和出口面上物理量完全一致,简单的来说,即出口流出的继续从入口流入,采用这种边界条件的意义主要是用来模拟无限长的通道。这样的话会产生一个问题,对于压力p来说,进口和出口完全一致的话,就没有压力差了,那我们如何用压力差进行驱动流体流动呢?事实上,fixedJump命令就是解决这个问题的,我们可以设置一个Jump量,使得入口处的压力等于出口处的压力加上一个常数,这样就可以产生压力差。可以看出,达到了我们想要的效果,进口壁出口处压力大一个定值,且正好等于我们设置的值。原创 2023-01-14 23:50:38 · 573 阅读 · 0 评论 -
OpenFOAM里fieldAverage工具
resetOnRestart的值决定当solver继续运行时,是否要读取最近一个时间步的meanField的值来计算接下来时刻的时均值;3、由2中,我们可以得出,如果设置的deltaT为1,timeStart也为1,那么1时刻生成的U文件和UMean文件是相同的。1、如果不设置window命令,base后的time和iteration效果是一样的,都是从timeStart开始进行平均。但是不重要,我们知道不对其进行设置,计算出的结果是对的。这两个命令不对其进行设置。链接中对如下两个命令的解释不是很明白。原创 2023-01-12 21:36:21 · 764 阅读 · 2 评论