OpenFOAM

• Z

  在颗粒流中的动量转移大概有多少？
  • zhangxc0223  

  1
  1
  Posts
  71
  Views

  Z

  抱歉问的问题可能太泛了，具体举个例子来说。比如说在直径D=0.01m的空气管流中，速度Uf=15m/s，雷诺数Re=10000。颗粒的斯托克斯数St=10，直径是大概dp=50×10^(-6)m，密度为ρp=1000kg/m^3，颗粒-空气的质量比是Φ=0.1，双向耦合，只考虑颗粒曳力。 在这种情况下，我查看OpenFOAM计算得到的kinamaticCloude:UTrans文件中，动量转移都是在10^(-14)量级（单位是kg×m/s），请问这合理吗，为什么只有这么小。。？ 据我所知，这个UTrans文件表征的就是颗粒→流体相的动量转移，也就是颗粒相通过源项Sp，施加在流体相的曳力作用。不知道这个理解对不对？
• S

  修改文件内容的shell
  • Samuel-Tu  

  8
  8
  Posts
  154
  Views

  S

  最近多次使用foamDictinary感觉真是太用了，自己在文件里面加的keys也能识别出来并修改
• W

  拉格朗日颗粒喷射
  • wangtao121  

  1
  1
  Posts
  83
  Views

  W

  各位大佬问一下有关拉格朗日颗粒喷射的问题，我如何往喷射模型中（比如manualInjection模型）动态的传递喷射点坐标，坐标每一个时间步都会发生变化，怎么实现内部传递呢。
• 

  请问，MPPIC这个方法的particle distribution function的运作
  • zhe  

  11
  11
  Posts
  253
  Views

  

  @zhe 在 请问，MPPIC这个方法的particle distribution function的运作 中说： 固相是由particle distribution function表示的，并且固相的相分数也是通过它进行求解的。那么这个方程在OpenFOAM里有体现吗？ 没有，你说的这个方程是一个非常基本的控制方程，其实就是玻尔兹曼方程，很多模型都是从这个模型上展开对他进行封闭的。比如 双流体模型，就是粒子分布函数传输方程取2次矩的模型。 矩方法就是取高阶矩的模型。 MPPIC在这个方程上对碰撞进行BGK封闭。可以更简单的处理碰撞。看这个文献的方程2 A model for collisional exchangeing as /liquid/solidfluidizedbeds 《无痛苦N-S方程笔记》4,8节也有讨论。这个方程没人去解。自由度太多了。 可以彼此融会贯通 左右逢源一下
• S

  of输出continuity error残差
  • Samuel-Tu  

  3
  3
  Posts
  105
  Views

  S

  @马乔 有没有什么办法输出这个误差呢，一般大家是怎么判断pisoFoam或者pimpleFoam是否收敛呢? 我发现误差降到１e-5以下就很困难了。。。所以我是设置的监测点，发现关注的变量变化不大了就认为收敛。。
• 

  关于程序调用问题
  • win  

  4
  4
  Posts
  141
  Views

  

  你这是什么版本啊，为嘛那一通运算没有量纲检查。forAll这个只会遍历内场。
• 

  如何根据残差判断是否收敛？
  • Izumi  

  12
  12
  Posts
  4431
  Views

  S

  看论文发现有作者说fluent稳态计算残差降到1e-7则认为收敛停止计算，但是我用of的simpleFoam的时候发现残差根本降不到那么低，连降到１e-６都已经是极限了。我都不敢在论文里写残差判断标准了，一般你们怎么说收敛标准呢
• N

  请教一下，各位前辈如何在paraview中多个窗口中输出不同时间的状态
  • nanxuan  

  5
  5
  Posts
  149
  Views

  S

  @bestucan 好滴，我试一下
• S

  shm的stl文件要求
  • Samuel-Tu  

  3
  3
  Posts
  90
  Views

  S

  @xpqiu 好滴
• 

  cfMesh
  • 同学博  

  5
  5
  Posts
  159
  Views

  

  @Samuel-Tu 用了，确实快，但是有时候分网不成功，而且网格质量很难保证
• X

  浅水方程模型
  • xiaolin  

  2
  2
  Posts
  66
  Views

  

  本来想说试试加深水槽。 静水流深嘛，一看是浅水方程…… 如果边界条件设定为带反射的，那就应该不停波动。如果是自由边界，不停波动不正常
• S

  snGrad非正交修正corrected和fvSolution里的非正交迭代
  • Samuel-Tu  

  2
  2
  Posts
  66
  Views

  S

  在cfd online 里逛了一下，有些收货。 似乎non-orthogonal corrector 是在组装矩阵Ax=b的时候用的，它会把扩散项中的正交部分放入系数矩阵Ａ中，而把非正交部分放入方程右边的源项b中，使得系数矩阵更稀疏且主对角线控制，使得计算更稳定？？cfd online
• 胖

  关于RAS Model里k-Epsilon中k和Epsilon的取值问题~
  • 胖小帅  

  5
  5
  Posts
  139
  Views

  胖

  @cccrrryyy 受教受教，授我以渔，万分感谢~ 帖子链接里开篇就阐明了我的疑惑，也难怪k-epsilon可以在工程领域运用的如此广泛~ 我就先对照着openfoam tutorials和CFDEM案例里的k-epsilon取值作参考，计算后再做调整，再次感谢指点~
• M

  添加源项fvOptions编译出错
  • merlinizlf  

  10
  10
  Posts
  249
  Views

  M

  @东岳 非常感谢老师 ！ 就是这样迂回的话好像有些难，我还要再思考思考，再尝试尝试
• S

  bluecfd
  • Samuel-Tu  

  10
  10
  Posts
  208
  Views

  S

  @东岳 好滴
• K

  source term only in boundary cell.
  • kimy  

  2
  2
  Posts
  55
  Views

  K

  Already solved.
• K

  New scalar in phaseModel.c
  • kimy  

  30
  30
  Posts
  288
  Views

  K

  @bestucan Thanks for your attention.
• V

  twoPhaseEulerFoam并行出错
  • veen  

  7
  7
  Posts
  148
  Views

  V

  @东岳 好的，谢谢李老师！
• 

  关于createField.H的问题
  • win  

  3
  3
  Posts
  82
  Views

  

  谢谢老师
• U

  openFoam5.x中，MPPICFoam后处理时无法读入粒子文件
  • upc_ngh  

  5
  5
  Posts
  141
  Views

  U

  @东岳 这个当时手一抖就点上了，发现没法删掉，真是尴尬
• 搬

  如何显示涡量图
  • 搬运工不好当  

  6
  6
  Posts
  3058
  Views

  小

  我用fluent做的也可以做出这种图来吗，可以的话用什么做后处理呢 万分感谢
• 

  OF7 更新的buoyantPimpleFoam 使用Boussinesq不正确
  • sibo  

  5
  5
  Posts
  121
  Views

  

  @东岳 还没毕业。。都是泪 结冰也在搞，也和别人合作搞搞别的。希望早日上岸。
• 莫

  OpenFOAM platForms 问题
  • 莫晓柒  

  2
  2
  Posts
  60
  Views

  X

  @莫晓柒 编译之后得到的库和求解器程序都在 platforms 下面，所以只能重新安装了。
• S

  of指定CPU计算
  • Samuel-Tu  

  5
  5
  Posts
  151
  Views

  S

  @东岳 我不算运营人员，只是参与了号主的一些活动，有时候帮帮忙。。。
• 

  particleforce类内ErgunWenYuDrag的表达式感觉有错误
  • 上级  

  22
  22
  Posts
  1983
  Views

  

  @upc_ngh 1.是的，$\mathbf{U}$代表的是KinematicParcel中的$\mathbf{U_c}$，也就是连续相的在颗粒所处位置的速度插值，具体的插值方式OF里面有cell，cellpoint那些，$\mathbf{U}_\rm{d}$表示的就是颗粒的速度。 2.是的，系数乘以速度差。
• 

  paraview显示问题
  • 李东岳  

  5
  5
  Posts
  207
  Views

  M

  新知识点
• 沐

  代发：贴底圆柱绕流升力和阻力系数不正确
  • 沐风  

  5
  5
  Posts
  157
  Views

  S

  @沐风 因为是百度云，所以没有下载 不过我有几个小建议： 计算域是否够大，圆柱和上下前后边界的距离是否够远，去除边界的影响； snGrad建议用corrected； 我印象中上下边界应该用symmetry，用noslip是有什么原因或者依据吗。
• S

  Gauss linear upwind是二阶迎风格式吗
  • Samuel-Tu  

  3
  3
  Posts
  128
  Views

  S

  @xpqiu 感谢
• J

  Paraview 编译
  • Jacob Z  

  14
  14
  Posts
  305
  Views

  J

  @chien 是的谢谢🙏，放弃编译了😹
• 小

  怎样限制流场中最大速度
  • 小考拉  

  2
  2
  Posts
  59
  Views

  

  @小考拉 maxU会突然增大只是发散的迹象，不是发散的原因，所以这种问题很难debug
• Y

  paraview设置使用用户设置
  parav • • yfclark  

  5
  5
  Posts
  130
  Views

  闻

  我也是每次都重新设置，很麻烦，前来学习
• 桎

  interIsoFoam求解问题
  • 桎梏  

  2
  2
  Posts
  52
  Views

  

  min(alpha.water) 都出现负值了。 时间步自己调整的特别小是为了让库朗数降低，结果时间步那么小，库朗数还那么大，可能网格有问题
• B

  有没有人用颗粒求解器计算过带AMI面的问题的？
  • bubbleEL  

  7
  7
  Posts
  449
  Views

  液

  @bubbleEL 谢谢你
• 波

  请教：不同方向的旋涡相遇能保持继续保持自己的旋转方向？
  • 波流力  

  2
  2
  Posts
  54
  Views

  

  尝试解释一下 就是0.3图里的涡在0.4不见了，然后1里又出来了？ 0.3里刚过柱体，柱体后面压力小，然后左右流体往中间挤，挤不动了，要么往上要么往下，往下是墙，只能往上冒，就是0.3里的涡 0.4里不见了，因为这个柱体下面有个空隙，流体在过空隙的时候很挤，过了空隙就释放压力，下边是墙，是能往上释放，然后空隙里出来的流体就把0.3里的涡往上托，托到液面托没了就成0.4了 然后0.6，为了释放压力往上托涡，但又因为惯性，释放过头了，涡把流体都往上带，下面就出现低压区了，然后左右的流体都往中间流，在这里碰头，可以从0.6看见他们碰头后很纠结，不知道往上走还是往下走 然后是1，往上的那部分变成了0.3的涡尾，往下的那部分，又卷出来个涡。 要是水足够深，应该能同时看见0.3和1的涡共存。 再形象点，就是0.3里的涡是柱体后的低压区造成。1里的涡就像柱体下的狭缝放的“屁”，一个屁把0.3里的涡挤走了……
• 液

  液固两相流动+滑移网格颗粒运动求解的问题
  • 液固两相流  

  1
  1
  Posts
  70
  Views

  液

  计算区域为旋转区域与非旋转区域，颗粒由非旋转区域入射，但是后期处理颗粒轨迹时候发现，颗粒无法通过动静交界面（cyclicAMI）. 在颗粒求解文件中设置交界面为none 请问大家有碰到过这问题吗？有什么好的解决方法吗？
• S

  请教在dam break里如何计算获得 wavefront的速度
  • Sloan  

  4
  4
  Posts
  116
  Views

  S

  @桎梏 好的，我试一下，谢谢建议！
• 

  使用interDyMFoam时，correctPhi应该在什么情况下打开？
  • tidedrinker  

  5
  5
  Posts
  118
  Views

  

  @bestucan 头文件里说的是如果非保守通量（ non-conservative fluxes，应该是翻译成非保守通量吧？）对算法的预测部分有影响的话是需要correctPhi的，但是什么情况下会影响，我完全没有概念。至于后半说多相流的部分，4.1多相流里不开correctPhi的还是有几个的。所以我完全没提取出什么有效信息。
• 

  关于二维模型的计算，这两种snappyHexMesh方法有什么不同，应该采用哪一个？
  • DY大世界  

  16
  16
  Posts
  348
  Views

  

  @xpqiu 非常感谢大佬的回答，那不用extrudeMesh命令所得到的确实是不正确的二维网格。而且我对比了使用extrudeMesh与不使用extrudeMesh的计算案例，得到的结果也是不完全吻合。
• C

  自适应网格技术生成液滴
  • conghongchuan  

  3
  3
  Posts
  144
  Views

  桎

  可以用paraview直接算出自适应网格得出的球体体积，然后跟理论值对比一下，求出两者的相对误差就可以了
• C

  后处理如何获得破碎液滴的粒径和速度分布？
  • conghongchuan  

  3
  3
  Posts
  77
  Views

  C

  @bestucan 好的，非常感谢老师的回复
• S

  请教divSheme中设置对溃坝自由液面形态的影响, RNG k-e湍流模型
  • Sloan  

  36
  36
  Posts
  455
  Views

  S

  @波流力不好意思后来赶ddl耽搁了。后来我测试了下，您的算例运行的不错，液面形态很好。然后参考您的设置更改了我的计算，液面有提升，但是在撞到右侧壁，上升时候还是出现了裂隙。然后在您的算例里改成我的case的几何参数，裂隙现象就没了。。。很奇怪，我再研究看看
• X

  movingConeTopoFvMesh使用的一些疑惑
  • xiaolin  

  6
  6
  Posts
  97
  Views

  史

  @xiaolin leftEdge没有意义，leftObstacleEdge（rightObstacleEdge）是指运动的左（右）侧边界的位置 这个类其实是一个不太好用的类，你自己要改一下。这个建议你先去看一下源码，了解里面的操作流程，然后修改一下就可以自己用了，亲测可以
• J

  Not enough slots available
  • Jacob Z  

  1
  1
  Posts
  53
  Views

  J

  大家好！运行一个CFDEM案例报错”There are not enough slots available in the system to satisfy the 20 slots that were requested by the application “，意思大概就是说分配的slot达不到这个算例运行的要求……。有老师遇到过这个情况吗？如何解决的呢？求指教🙏!
• 

  paraview处理进口段波动
  • chengan.wang  

  4
  4
  Posts
  69
  Views

  

  估计是一个吹出来的波，一个荡出来的波，叠一块了，想把荡出来的波滤掉，只看吹出来的波。 滤波器一般按频率滤波，不是按看见的特征滤波。 你花出来的曲线已经是叠加波了。如果不加吹的算出来的才是你要滤掉的部分。这部分应该很规律，能写成三角函数，那么用calculator 把物理量减去三角函数就行了，三角函数的自变量是x的坐标。 如果就要减去这个曲线图，得拟合出来这个函数，在 paraview 里处理才方便。 其他麻烦的方法倒是一大堆，paraview 里有 python。或者导出，用matlab之类的处理。
• M

  请教一个周期性边界的设置问题
  • Mania  

  11
  11
  Posts
  237
  Views

  F

  @东岳 我在速度文件力设置好内部速度为定值（0.6），压力文件里设置内部压强（1e5）通过一段时间计算后，流速越来越大。像我需要指定初始速速，只需要设置速度文件，压力文件中设置为0就可以，或者将指定的初始速度转化为相应的压力添加在压力文件中，速度文件设置为0呢!
• 林

  paraview能否显示一定速度区间的网格部分
  • 林之流风  

  9
  9
  Posts
  157
  Views

  林

  @xpqiu 感谢老师详细指点，试了下这个方法确实简洁
• M

  求助：付费咨询，微观多孔介质两相流代码问题
  • mutilphase flow  

  14
  14
  Posts
  543
  Views

  S

  这是没有编译成功吧，我在foam-extend-4.0+2.4.0上面编译成功过这套程序
• 

  v2f模型求解射流时，有时计算发散，有时计算正常，但射流集中在壁面附近，求问
  • strive_DUT  

  2
  2
  Posts
  41
  Views

  

  计算发散报错如下：
• 

  被后处理工具取代？
  • 刀尔東  

  5
  5
  Posts
  1372
  Views

  Y

  @刀尔東 你好，我想问一下，在后处理中想对一个量做加减运算，这个在后处理里应该是什么指令呢？比如对速度U的模取反，postProcess -func "-mag (U)"，这样是报错的，不知道应该怎么处理？感谢
• Y

  壁面边界条件零梯度的物理意义
  • Yu_Tian  

  5
  5
  Posts
  166
  Views

  Y

  @cccrrryyy 好的，感谢解答，我再思考一下