OpenFOAM

4.5k Topics 26.2k Posts
  • 1 Posts
    1k Views

    RT,我在做管流的模拟,现在想在一个横截面上做出雷诺剪切力(u'v' VS r/D)的曲线图象。现在用openfoam的sample功能提取了大量的速度-位置数据(100个时间步长, 每个时间步长在该横截面上有360个径向直线,在一共36000组,每组直线取了200个点的速度(u,v,w)的数据,所以现在在MATLAB上得到了一个200×36000的矩阵。现在想根据这个矩阵得到雷诺剪切力的图像,请问怎么操作呀?

  • 做wedge网格做了一个下午

    2 Posts
    1k Views

    这个是我画了一下午的wedge网格,先给所需,在这个基础上改吧........

  • 关于开发buoyantPimpleDyMFoam

    3 Posts
    2k Views

    @东岳 谢谢大神回答!

  • Ubuntu 18.04 LTS有人尝鲜么

    1 Posts
    935 Views
  • 求解器中如何使用私有变量

    5 Posts
    2k Views

    @浪迹天大 非常感谢您的帮助,我已经完美解决这个问题。

  • 2 Posts
    1k Views

    http://www.cfd-china.com/topic/1548/

    好像跟这个帖子有关。气动噪声算法上我还不太理解。目前只知道是计算压力小波的传递?

  • 观察到一个很有意思的局部奇点

    3 Posts
    2k Views

    你用的是非结构网格?很有意思 哈哈

  • 静止水中气泡上升流线图输出

    4 Posts
    3k Views

    在你画出来流线图之后,在你原始的.foam基础上,再做contour plot,然后同时显示流线和contour。

  • 小白求助OF安装问题

    4 Posts
    2k Views

    你试试运行这个命令:

    sudo chmod u+x /home/lnq/OpenFOAM/
  • non reflective可用于壁面吗

    1 Posts
    1k Views

    有时候相间的物质转换会造成局部高压,然后碰到wall就会反弹。。。
    一般non reflective boundary都是用于出口入口的,那么它
    可以用在wall上吗。

  • OpenFOAM如何截图

    4 Posts
    2k Views

    他是不是想截这个图 :mihu:

    替代文字

  • 3 Posts
    2k Views

    而且用的是压力基 两相流

  • 3 Posts
    2k Views

    谢谢!经过测试,完美地解决了我的问题!

  • 如何获取相邻节点的标量值

    3 Posts
    2k Views

    谢谢大神,我看到你原来的代码里有相关的。再次感谢。
    http://www.cfd-china.com/topic/1228/owner-neighbour-网格寻址/2

  • 虚拟机扩容

    6 Posts
    3k Views

    在外挂硬盘上计算非常慢,不写log可以,写log的话,在机械硬盘上写log严重耽误了计算的时间。

    另外,log文件有时候变的特别大,我这个算例时间步长比较小,写文件间隔比较大,导致结果没有多大,但是log上G。可以通过 > /dev/null来禁止写log,当然如果你确定不需要看log的话。

  • OpenFOAM安装大小

    1 Posts
    935 Views

    0_1525840011230_2018-05-09 12-25-51屏幕截图.png

    约2G左右

  • 并行计算后如何reconstruct结果

    3 Posts
    2k Views

    https://www.sciencedirect.com/science/article/pii/S0010465517303375

    See this dsmcFoam+ published on Computer Physics Communications.

  • 这种网格在openfoam适用吗?

    4 Posts
    3k Views

    5.5.1 fluentMeshToFoam
    ......
    There is currently no support for embedded interfaces and refinement trees.
    貌似对这个是不支持的。我也用过snappyHexMesh, 这种refinement的过渡处是会生成非结构网格的,也就是没有interface。现在snappyhexmesh是做了一些改进吗?

  • 9 Posts
    9k Views

    我的理解是fvc::div(phiHbyA)已经将原控制方程中的散度表示为高斯积分的求和形式了

    是的。正确。

  • 2D的DNS有意义吗

    2 Posts
    2k Views

    @yuan_neu 4如果从数学意义上来说没有错误。考虑到物理机理,如果模拟的是湍流,那么所谓的2D是错误的,因为湍流本身是三维的。

  • momentumPredictor 开关对比

    1 Posts
    1k Views

    momentumPredictor为动量预测步骤,对应链接文本中的方程14。动量预测步骤可有可无。从经验来看,在粘度比较大的时候,方程比较刚性,若附加可能会引起发散。在多相流以及动网格算法中,通常此项关闭。具体原因,目前并没有见过相关文献讨论动量预测对收敛的严格证明。

    目前尚不清楚动量预测的普适性方案,后续还有待系统研究。下面的算例是一个表示开启动量预测较好的算例。但不具有普适性意义。算例模拟一个静止的流体。纯测试算法,理论上模拟结果速度应该为0。模拟中考虑了重力。采用p方程计算。

    首先,俩种结果的压力场(右图)均正确。均体现了由于重力导致的静水压力的渐进变化。但速度场都不是很正确。这是因为在附加重力场的时候,用的是p方程而不是p_rgh方程。类似讨论: http://www.cfd-china.com/topic/394/如果用dpmfoam求解稀相流会怎么样-误差大么/37

    然后,momentumPredictor开启的时候,假拟速度更好一些,如下图。

    momentumPredictor off

    0_1525236812110_pic2.png

    momentumPredictor On

    0_1525236812110_pic2.png

  • 4 Posts
    3k Views

    谢谢东岳老师考古!讲解非常清楚!

  • 谢谢李老师

    3 Posts
    2k Views

    @东岳 谢谢您的回复。我的边界条件是进口速度fixvalue,进口压力zeroGradient。出口速度inletoutlet,出口压力0。因为我之前AMI几何形状不一致,后来才改的。是不是除以0就是因为AMI有一部分没有贴合上。

  • OF下怎么设置弹性边界?

    3 Posts
    2k Views

    @东岳 谢谢东岳老师,这个据我所看文献现在of并没有植入,植入很麻烦。现在也不是我研究重点,等有时间在弄这个。

  • 2 Posts
    2k Views

    虚拟机肯定是没问题,跟Ubuntu真系统感觉没什么区别。
    但是BSL下面的OpenFOAM能感觉出来有点慢,之前对比过,主要是来自于磁盘读写。如果你有兴趣的话,可以试试对比下,谢谢。

  • OpenFOAM3.0.0的一个小bug

    9 Posts
    4k Views

    @浪迹天大 问题我已经解决了,就是说明一种现象

  • 关于AMI的疑问

    2 Posts
    2k Views

    是的,cyclicAMI是CFD中的滑移网格,试试添加关键词lowWeightCorrection 0.2;。可以上传你的图片看看

  • 4 Posts
    4k Views

    参考interDyMFoam下面floatingObject算例,用topoSet生成中间的小方块

  • 11 Posts
    6k Views

    可以用输入边界法造波,比如IHFoam,waves2Foam等,现在的of也自带边界输入法造波了,只是功能还比较单一。

  • 12 Posts
    7k Views

    @东岳 在 DPMFoam求解器颗粒计算有颗粒跑出计算域外! 中说:

    粒速度持续

    假如设定反弹,颗粒轨迹与壁面有交叉,但还是跑出了计算域,正如之前碰到的情况,还是算法本身不够完善,openfoam5.0以经把颗粒追踪那一部分的算法完全更新了,颗粒追踪不会失效了

  • 3 Posts
    2k Views

    @东岳 对,后来查看了一下userguide 中的solver部分,里面PCG的全称是preconditioned conjugate gradient,就一下明白了

  • 8 Posts
    4k Views

    因为你肯定没回头去把手册或者视频看一遍:xinlei:

    从哪算取决于starttime的设定和有没有那个时间节点的文件夹。

    看是不是只能从0.8算,那就mapfields过来后,从0看看能不能算;算了之后0.8有没有更新。

    至于看0用没用上,把0里的边界条件瞎胡改看看抱不报错不就知道了

    高中学习和大学考试还是那种“我一恼”“我一生气”就做到了的模式。往后就是“临渊羡鱼不如退而结网”。重点是“退”,死磕的话。。。。不是有个成语叫“磕(客)死他乡”:wocao:

  • simpleFoam求解二维翼型Cd不变

    19 Posts
    8k Views

    @dzw05 湍流模型中有一个运动黏度nu,正常情况下在后面赋值即可。我当时在网上查询如何提取升力系数和阻力系数时有篇帖子提示要在nu后面先再输入一个nu字符再赋值,由于之后的计算我都用的是这个错误的湍流模型文件,导致阻力系数求解一直出现问题。后来一个一个文件排查的时候发现是这里的问题。

  • 2 Posts
    1k Views

    paraivew据我所知不能。但在C++里一切皆有可能,但是可能你需要自己写点代码。

  • 怎样输出interFoam中的项分数?

    3 Posts
    2k Views

    @dzw05 明白了,也就是说,water并不是一个真正原始的量,这个量的值实际上是由alpha1来指定的,那么我想输出water,直接就输出alpha1就行了。我刚刚试了一下,确实是这样的。多谢啦~~~~

  • OpenFOAM的2.0.x版本是什么意思?

    7 Posts
    4k Views

    @bestucan 那也是学长,厉害了~

  • 4 Posts
    3k Views

    @东岳 这个模型是根据实际的物体建的,所以模型重心就直接取了实际物体的重心坐标,没有数学公式。

  • 有关p_rgh收敛的问题

    17 Posts
    7k Views

    @dzw05 谢谢您的解答,讲解的很明白,我还要看一看,如果有新想法,会回复您的。

  • overInterDyMFoam求解器报错

    7 Posts
    4k Views

    明白了,谢谢!

  • 使用mapfields时出现了这个报错

    7 Posts
    3k Views

    @东岳 谢谢东岳的耐心:baobao: :baobao: 感动!

  • 6 Posts
    2k Views

    @肖恩曼 看了下paraview手册,实现的逻辑一样,就是用python实现。画图可以直接调用paraview,比用其他绘图软件方便。

  • 6 Posts
    3k Views

    但是我设置数组变量的时候也是您这样设置的,然后这样怎么能够分别设置边界条件呢?

    我知道了,你要在创建数组的时候,自动给outlet,inlet之类的创建边界条件?这没法实现。只能设定为MUST_READ,然后手动在0文件夹里面设定。或者,就是在求解器内部强制赋值某边界条件。

    而且我要的变量场是因变量,就是需要在0文件中设置初始时刻的边界条件。

    这和上文不是相符的么?

  • OpenFOAM 中的refineMesh

    5 Posts
    3k Views

    @东岳 好的,李老师。我仔细看看ICEM。

  • OpenFoam 如何设置稳态算例?

    3 Posts
    2k Views

    嗯嗯,好的,谢谢岳哥!

  • 6 Posts
    3k Views

    @bestucan 谢谢你~我去仔细研读一下~

  • 段错误 (核心已转移)

    2 Posts
    1k Views

    只能看出来是发散了,无法更深的debug,对结果影响很严重。

  • simpleFoam求解NACA0012翼型

    1 Posts
    1k Views

    学生在用simpleFoam 求解NACA0012翼型的升力系数和阻力系数时,所得结果Cd100多,Cl 20多。学生采用速度1.6m/s,攻角8度,弦长1m,SA模型。U文件中入口和出口类型均为freestream,翼型为fixedValue;P文件中入口和出口类型均为 freestreamPressure;nut和nuTlida文件中参数为0.14。controlDict文件密度输入1.225,lref1,aref0.08。学生是刚学习openfoam,很多东西还不是很清楚。请各位老师帮忙看一下问题出在哪里,不胜感谢!

  • OpenFOAM对标量求散度?

    3 Posts
    3k Views

    @东岳 大致明白了,这里相当于是用新的算法代替了传统的散度操作,从而在能达到原有求散度目的的情况下避免了直接求散度造成的麻烦。多谢前辈!!!

  • 关于无粘流

    8 Posts
    3k Views

    @东岳
    谢谢,我也做了类似的仿真,相比sonicfoam, rhocentralfoam似乎结果得到的时间更短。

    对于简单的问题,例如无黏流流经wedge一类,两个求解器的结果其实是相似的,而且较为满意。

    但是,不知道您有没有发现,在这个喷管的模拟中,在打洞的地方,喷管的收缩段似乎出现了速度为0的区域,很像边界层的感觉。

    然而我做的是inviscid flow我设置mu=0。我看了rhocentralfoam的代码,感觉rhocentralfoam是能做无粘的,但相似的现象仍然有。不知道是哪里的问题。