Skip to content
  • buoyantBoussinesqSimpleFoam求解器

    OpenFOAM
    3
    0 赞同
    3 帖子
    2k 浏览
    李东岳

    https://cfd-china.com/topic/8018 可以上传看看

  • 0 赞同
    5 帖子
    8k 浏览
    F

    可以设计一个量纲的类,里面包含一个整形数组来记录量纲的分量,如[L, M, T]还要定义量纲的运算。然后设计一个标量的类,里面除了浮点数据之外,再加一个量纲,定义标量的动作时也要注意其中的量纲计算。同理还可以定义向量和张量。比如距离的量纲是[1, 0, 0],速度的量纲是[1, 0, -1]。

  • 0 赞同
    2 帖子
    3k 浏览
    S

    您好,请问您发散的问题解决了吗,我现在也遇到了相似的问题。就是PBM模型里面的相间作用力我只能考虑曳力,添加其他作用力就会发散,卡了好久了。

  • checkMesh failed,求助

    Meshy
    6
    0 赞同
    6 帖子
    9k 浏览
    Y

    @东岳 老师,我重新画网格后,decomposePar经常会报错,但是不影响后面计算,没有查到原因0_1537492983872_0403f9aa-ee80-4134-a228-ac6289ba58de-image.png

  • 0 赞同
    1 帖子
    1k 浏览
    P

    我在引射器模拟中,RP_Get_Real("operating-pressure") + C_P(c, t)为什么小于 fluent云图中的absolute pressure的值
    0bc4d2ecfd45dabedb282d4c8ac5477.png 528e7412a62c8d3c22a9ed7c2706a23.png

  • 0 赞同
    5 帖子
    8k 浏览
    李东岳

    这个其实是个好话题,很多大佬在sci里面也明确说了,LES之类的模拟,如果对于这种类似的圆管,施加cyclic边界条件,会对湍流产生抑制。但是我也没见过大佬们进行对比。

    另一方面,感觉边界条件,直觉感觉,单单的一个patch引起的差异,不至于那么大。LES模型参数的影响更大。如果谁有类似的文章,可以发来看看。

  • 0 赞同
    4 帖子
    6k 浏览
    C

    @allanzhong 看这个意思,应该和你说的差不多,理论上肯定为,但实际上从理论转换为数值算法的时候总会有些偏差,算法在具体实现的时候又会有一些偏差。如果是这样,这些偏差我认为是可以接受的。

  • 0 赞同
    3 帖子
    5k 浏览

    这个问题已经解决,麻烦大家了

  • Omega涡识别算法

    Algorithm
    3
    0 赞同
    3 帖子
    5k 浏览

    @李东岳 里面的(u)和(y)我想选用特征速度和特征长度,就是物体的运动速度和物体的特征长度,然后du/dy想在流场中提取

  • VOF问题,自由液面无理由偏移

    Fluent
    1
    0 赞同
    1 帖子
    1k 浏览
    F

    FFF--001--0728(1).jpg
    一个容器加一个U型管,上面两个出口为自由流出口,下部为一个液相速度进口,初始液面为U型管出口位置所在的水平面(图中红线位置),随着计算进行,液面开始倾斜,后期计算就无限制的倾斜发散,重力设置方向为Z方向-9.81。
    请各路大神指教,在线等

  • 纳米尺度下该如何模拟两相流动?

    CFD彩虹条
    4
    0 赞同
    4 帖子
    6k 浏览
    李东岳

    理论上了解一些求解方法,软件没有用过 :chitang:

  • 0 赞同
    1 帖子
    3k 浏览
    袁宝强

    请问各位大佬,对于可压缩流大涡模拟的总湍动能如何定义呢?

    对于不可压缩,大涡模型采用的亚格子模型采用的雷诺时均的形式,因此可以用雷诺时均计算大涡湍动能和亚格子湍动能,两者的和就是总湍动能。
    $K=0.5*(\overline{uu}-\bar{u}\bar{u}+\overline{vv}-\bar{v}\bar{v}+\overline{ww}-\bar{w}\bar{w})+\bar{k}_{sgs}$ 但是对于可压缩流的大涡模拟,亚格子模型采用的为法夫尔时均形式,此时如何定义总湍动能呢?我的理解是
    $\widetilde{K}=0.5*(\overline{\rho uu}-\bar{\rho}\widetilde{u}\widetilde{u}+\overline{\rho vv}-\bar{\rho}\widetilde{v}\widetilde{v}+\overline{\rho ww}-\bar{\rho}\widetilde{w}\widetilde{w})+\widetilde{k}_{sgs}$
    请问各位大佬,我的理解对么
  • 0 赞同
    6 帖子
    9k 浏览
    S

    @liuyk 神医啊,关闭超线程马上就解决了

  • interFoam 模拟水滴运动

    OpenFOAM
    10
    0 赞同
    10 帖子
    12k 浏览
    1

    学到了

  • 加班文化

    CFD彩虹条
    5
    0 赞同
    5 帖子
    8k 浏览
    李东岳

    今天青年节,下午放假 :chigua:

  • 0 赞同
    3 帖子
    3k 浏览
    李东岳

    coalChemisryFoam变密度的求解器无法计算常密度的情况吗?

    我确定这个可以。http://dyfluid.com/class.html#id4 我在这个算例里面做的就是一个基本常密度的,附加密度的求解器,压力方程跟coalChemistryFoam压力方程类似,。

    具体coalChemisryFoam这个我不太清楚了。这求解器没玩过。

  • 向CFD之父致敬,Brain Spalding

    CFD彩虹条
    2
    0 赞同
    2 帖子
    6k 浏览
    L

    缅怀。。作为CFD萌新,已经无数次地看过这位老教授的名字

  • 判定发散

    CFD彩虹条
    2
    0 赞同
    2 帖子
    4k 浏览
    李东岳

    并行计算时候可以,单核运算CPU占用太低,看不清楚。

  • 0 赞同
    2 帖子
    3k 浏览
    Prometheus10P

    博主您好,我用v1812版本得wedge网格旋转一周生成360°的3D网格似乎很顺利,没有遇到博主您说的问题,可能差别在于axis这里:

    nLayers 36; expansionRatio 1.0; sectorCoeffs { axisPt (0 0 0.0); axis (0 0 -1); angle 360; // For nLayers=1 assume symmetry so angle/2 on each side }

    但是仿真的时候就有点问题了。请教一下博主,我后续的计算是用wedge网格生成的3D圆管状网格进行两相流动计算,在对称轴位置处,界面波动非常剧烈,(举例来说就像包子的褶皱),请问博主您是否有解决类似问题的经验?

    2c3fb71d-127a-479f-b43f-1a41360f7552-image.png
    9d8b3854-1345-4cb7-bf4c-5159ab3e4ea2-image.png

    感谢。

  • 0 赞同
    32 帖子
    54k 浏览
    李东岳
    # output particle properties dump dmp all custom/vtk 5000 post/liggghts_init_*.vtk id type x y z vx vy vz fx fy fz omegax omegay omegaz radius ERROR: Invalid dump style (/home/dyfluid/CFDEM/LIGGGHTS/src/output.cpp:589)

    这个看起来是后处理输出什么变量导致的
    # output particle properties