Skip to content

OpenFOAM

OpenFOAM交流区

5.4k Topics 31.8k Posts
  • 4 Posts
    4k Views
    李东岳

    替代文字

    你说上面这个图。这个算力对应的不是这个图。所以不一样。不过如果你增加近气速度的话,应该也是可以的。你加大10倍试试。

  • 3 Posts
    3k Views
    Z

    您好,请问这个问题您是否有好的解决办法了吗?
    强制设为0可能不太对,反映出来的可能不是压实效果。。

  • 关于量纲对齐

    5
    5 Posts
    4k Views

    @李东岳 东岳老师,这个主要是因为我的是粘弹性流体,在动量方程里边加入了弹性相,不只有雷诺数,还有魏森贝格数,可能是这两个无量纲数导致的:135:

  • owner neighbour 网格寻址

    3
    3 Posts
    4k Views
    李东岳

    4年过去了,当时发完这个贴就出去玩了。:jingya:

  • 4 Posts
    3k Views
    李东岳

    @ruibin_tongji 你看你的log是不是有些迭代次数太多了

  • 6 Posts
    4k Views
    Z

    @是嘿嘿牛啊 密度是被动地,是根据流体种类、温度和压力变化的,你可以设置上下两层温差或不同的流体种类

  • 湍流模型增加一laplacian项问题请教

    3
    3 Posts
    3k Views
    W

    @李东岳 噢噢,谢谢东岳老师。

  • 从中间某一时间步开始计算发散

    4
    4 Posts
    3k Views
    李东岳

    @shaneheee 那就不用管了。

    有时候加上一行啥用没有的代码,算例就能用。那就不用管了,放在这。

  • openfoam 流向做空间平均

    3
    3 Posts
    3k Views

    @xpqiu 太感谢您了 您教的好细, 非常感谢 ,对我帮助非常大 ,谢谢您:140:

  • 1 Posts
    1k Views

    各位老师好!
    报错信息:
    Screenshot from 2021-11-20 14-31-59.png

    solidParticleCloud.C的代码:

    solidParticle::trackingData td(*this, rhoInterp, UInterp, nuInterp, g.value()); Cloud<solidParticle>::move(*this, td, mesh_.time().deltaTValue()); if(mesh_.time().value()> td.cloud().tInjStart_ && // mesh_.time().value()< td.cloud().tInjEnd_) {this->inject(td);} void Foam::solidParticleCloud::inject(solidParticle::trackingData &td) //. { //. label cellI = -1; //. label tetFaceI = -1; label tetPtI = -1; mesh_.findCellFacePt(td.cloud().posB_, cellI, tetFaceI, tetPtI); if(cellI > 0){ solidParticle* ptr1 = new solidParticle(mesh_, td.cloud().posB_, cellI, tetFaceI, tetPtI,td.cloud().dB_, td.cloud().UB_); Cloud<solidParticle>::addParticle(ptr1); } }
  • 关于并行计算问题

    5
    5 Posts
    4k Views
    F

    @bestucan 嗯嗯,谢谢老师,我试试:146:

  • awk提取某一行的信息

    5
    5 Posts
    4k Views
    李东岳

    @bestucan 你们都溜的很啊!

  • topoSet的faceSet问题

    3
    3 Posts
    3k Views

    @是嘿嘿牛啊 按照我的理解,这样分的size不应该都是10吗?为什么size会变来变去呢?求大佬指点一下,拜托啦!

  • 三角形面积计算公式

    6
    6 Posts
    5k Views
    李东岳

    好像我问的有问题,$A=0.5(\overline{p_1}-\overline{p_0})\times(\overline{p_2}-\overline{p_0})$这个东西看起来是个矢量。但是面积不是矢量啊。我回去对对OpenFOAM的代码。

    我看了下,OpenFOAM里面mesh.Sf()对应的是这个:$0.5(\overline{p_1}-\overline{p_0})\times(\overline{p_2}-\overline{p_0})$,mesh.magSf()对应的是这个:$\left|0.5(\overline{p_1}-\overline{p_0})\times(\overline{p_2}-\overline{p_0})\right|$,后面这个应该就是面积了。

    @evensun $\left|0.5(\overline{p_1}-\overline{p_0})\times(\overline{p_2}-\overline{p_0})\right|$ 这个面积可以理解不

  • 1 Posts
    1k Views
    M

    大家好,我想问下在rhoCentralFoam中想植入一种新的重构格式,文献中以标量守恒律作为控制方程,对其的未知量进行了界面左右值的重构,这个应该在rhoCentralFoam中哪一块进行植入呢,我不太清楚,望解答一下 :yes: ,谢谢!

  • icoFoam 计算不出涡

    6
    6 Posts
    4k Views
    Z

    @冠竹 有可能是入口速度太小了

  • openfoam网格求教

    3
    3 Posts
    2k Views
    H

    @tens 好的!非常感谢您

  • 3 Posts
    3k Views

    @李东岳 可以把日常的流体都看成稀薄流体呀。从各种benchmark来看的话,精度还不错,处于同一个水平。
    现在主要的问题是,显式LBM并行效率特别的高,然后我就寻思着有限体积也可以显式呀,谁怕谁。结果性能卡在了压力泊松方程的求解上:136:

  • 调用实际运算步长

    3
    3 Posts
    2k Views
    疏影横斜水清浅

    @tens 好的,谢谢老师

  • 2 Posts
    865 Views
    L

    已解决,postProcess -func "components“ 应该是双引号:mihu: