Skip to content

OpenFOAM

OpenFOAM交流区

5.3k 主题 31.2k 帖子
  • 超音速下液膜演化问题--可压缩多相流求解器

    7
    7 帖子
    1k 浏览
    E

    明白了,谢谢东岳老师,我接下来试试3D模型。老师那您觉得是用层流还是湍流比较好呢

  • snappyhexmesh网格不是正方形

    5
    5 帖子
    1k 浏览
    J

    {5531D08E-7EEB-4a66-8334-0712B6413B99}.png
    对slice 和clip 看网格用crinkle 看 会显示真实网格 不用的话 切面在网格内布 就会有杂线

  • 4 帖子
    544 浏览

    @李东岳 李老师,计算过程中没发散啊,第一幅图我是在 paraview 中用了 ISO volume 过滤器出来的液相图,也搞不懂是啥原因提取出来的东西并不连续

  • 城市工厂流场CFD流场数据集

    3
    3 帖子
    922 浏览
    李东岳

    我这个8个图,最笨的方法纯自己手搞的。你用宏更好,我没时间研究。有时候研究了半天就过去了。直接开8个算例完事了。

  • 实时神经网络代理模拟技术的补充

    5
    5 帖子
    1k 浏览

    可能因为您用的激活函数是ReLU,正好可以很好地拟合甚至泛化y=x,但如果学习一个高次多项式,损失就不一定很小了。

  • 大步长下变形网格发散

    13
    13 帖子
    3k 浏览
    L

    十分感谢您的分享!

    @Mina-Lee 在 大步长下变形网格发散 中说:

    换句话说,是不是overset边界,根本不影响变形。fixed边界,才会影响变形。

    确实是啊!displacementLaplacianFvMotionSolver类中组建cellDisplacement场时,只对pointDisplacement为fixedValue类型的patch进行了处理(改为cellMotion类型,并对face各点求平均),其它类型的没做任何处理;最后cellDisplacement再插值回pointDisplacement时,点位移场只有原本的fixedValue边界没被修改,其它的都被求解结果修改了。

    我看代码时没深想到这一层,只在用的时候发现有时远场使用fixedValue和slip效果不同,我在这个贴子中的回复属于是瞎猫碰上死老鼠了,看来远场只要不是fixedValue,给任意都行。

    最后我的问题还是用的单一区域网格以及求扩散方程变形成功实现的,为了验证还用StarCCM算了比过,吻合得很好。倒是在StarCCM中试过重叠网格,overset交界面网格难匹配,常常挖洞出错,算出来的结果锯齿也很厉害。至于rbf,有些懒得折腾了,以后遇到非用不可的情况再说。

    这句话应该是您的“泣血之言”,指给后来的人看看:146: :146:

    @Mina-Lee 在 大步长下变形网格发散 中说:

    如果你要用重叠网格,记住,不要用逆距离,OF的这个方法存在BUG!

  • 关于计算船舶动网格相关参数问题

    6
    6 帖子
    1k 浏览
    L

    @李东岳 好的谢谢东岳老师

  • openfoam动网格加密后 计算很慢

    5
    5 帖子
    1k 浏览
    李东岳

    那是不太合理,可能你67万的要发散了

  • 如何基于blockMesh生成的圆柱网格添加边界层

    7
    7 帖子
    2k 浏览
    李东岳

    你分成了5块,并且4块都是边界的block,你要是改的话,那4个应该一致

    要不你就用商软做网格吧,3分钟出来了

  • 可压缩两相流问题(超声速下液膜演化)

    9
    9 帖子
    2k 浏览
    E

    2dd00f31c9f8005c2ff9c532bf41087.png
    主流速度较低的时候还是有液膜破碎分离的,上图是300m/s的工况,但是速度提高之后破碎现象就渐渐消失了,变得十分平滑

  • 3 帖子
    560 浏览
    R

    @李东岳 感谢李老师!原来没有特意写出来的都是同一个:zoule:

  • CFDEM粒子间黏附力怎么表示

    1
    1 帖子
    307 浏览
    C

    在CFDEM中,颗粒间相互吸引的作用力是使用skjr模型来实现的吗,还是说通过力来实现。sjkr模型中内聚能密度又是什么呢?

  • 3 帖子
    666 浏览
    I

    @lwjetmann 多谢,我自己增加了一个solver,确实是可行的。

  • bound 函数的实现方法

    7
    7 帖子
    4k 浏览
    C

    @李东岳

    帮李老师补充一下OpenFOAM里面average的实现,对于vol场要先插值成surface场

    template<class Type> tmp<GeometricField<Type, fvPatchField, volMesh>> average ( const GeometricField<Type, fvPatchField, volMesh>& vtf ) { return fvc::average(linearInterpolate(vtf)); }

    所以|3|0|2|到average函数里面以后要先插值成

    | 3 | 0 | 2 | 1.5 1

    再基于

    $v_{celli}=\frac{\sum (|\bfS_f| v_f)}{\sum |\bfS_f|}$
    得到中间网格的值为1.25

  • 5 帖子
    1k 浏览
    G

    @尚善若水 谢谢大佬,笔芯

  • 8 帖子
    5k 浏览
    S

    @李东岳 谢谢李老师,我按照您的方法已经解决了

  • 4 帖子
    741 浏览
    李东岳

    自己写的求解器别人很难debug
    有可能是边界条件没处理好

  • CFDEM运行时中断signal 9 (Killed)

    5
    5 帖子
    1k 浏览

    @李东岳 多谢老师解惑!:xinxin:

  • 关于dynamicMeshDict中的参数

    3
    3 帖子
    622 浏览
    Z

    @guanguan0922 非常感谢!

  • 势流算法问题若干

    1
    1 帖子
    341 浏览
    Y

    各位大佬,我在看东子老师的势流算法文章,有些地方思考了很久还是不得解。问题①:出口边界条件为速度梯度=0,方程退化成方程(5),求解困难;假设出口足够远,又可以通过方程(5)来求解了,这俩表述看的我一脸懵62cfaf3b-86df-414e-974f-0351bf0a52b5-image.png 问题②:“方程(5)和方程(6)可以实现相同的效果”,是指啥效果呀?我手动算了一下(见图),怎么保证两个c一样呢?∅_outlet的值如何确定呢?50cc50b2-4934-4b37-a443-2d507b865049-image.png 问题③:方程(10)我自己算了一遍,是个二次函数,和之前(5)、(6)都是一次函数相差太大了吧,不清楚(10)是如何能等于(6)8405b555-e186-4a2f-8ef1-9e5cc41fb6ec-image.png