Skip to content

OpenFOAM

OpenFOAM交流区

5.4k 主题 31.8k 帖子
  • OpenFOAM 结冰计算

    4
    4 帖子
    4k 浏览

    @sibo @东岳 多谢!

  • OF4.1这行代码释义,动网格问题

    4
    4 帖子
    5k 浏览
    V

    您好,您现在还有做SIXDOF得问题吗?有的话我们可以交流一下,谢谢啦

  • LPT-VOF coupling in OF 7

    3
    3 帖子
    3k 浏览
    Q

    @东岳 但是那个貌似是在v1906这种版本,我现在用的是OpenFoam 7, 貌似没有MPPICInterFoam

  • 关于interfoam出口压力设置的问题

    3
    3 帖子
    3k 浏览
    F

    @东岳 老师您好!是由于实际模型太大,为了减少计算量,截取了关键的部分作为计算域,但是由于是水下注汽,在实际模型中,该计算域最上方还有40mm的水深,我想设置一下该压力。

  • 进口通量等于出口通量的边界条件

    10
    10 帖子
    8k 浏览

    @东岳 嗯嗯嗯,我去在这个方向上面试一试,这些李老师

  • DNS求解VOF

    9
    9 帖子
    9k 浏览
    linhan.geL

    @东岳 感谢!

  • 请问OF计算的结果可以用其他后处理软件么

    4
    4 帖子
    4k 浏览
    N

    我转换ensight完成之后发现,只能读出几何模型,里面并没有数据,这是怎么回事

  • 11 帖子
    12k 浏览
    Jack QQJ

    您好!我遇到类似问题,准备采用上述的第二种方法即topoSet + createPatch,但不太清楚具体应该如何操作,希望能指点一下,谢谢!@赵一铭 @小龙

  • 怎么把mapfields的code写进solver里?

    2
    2 帖子
    2k 浏览
    L

    您好,这里两种网格求解的方程形式一样吗?您是研究多尺度模拟的吗?我也想尝试OpenFoam里植入一种多尺度算法,是否方便交流一下

  • reactingtwophaseeulerfoam?

    23
    23 帖子
    18k 浏览
    李东岳

    alpha was solved in twoPhaseSystem.C

  • createPatch创建cyclic边界条件的问题

    7
    7 帖子
    7k 浏览
    cfd_xfxC

    @一颗橙子 您好,想问一下您这个问题解决了嘛?后续是如何操作的。。我似乎遇到了和您完全一样的问题

  • Rhie-Chow interplation

    1
    1 帖子
    1k 浏览
    C

    Original RC interplation
    No time dependency, No relaxation of the velocity field:

    // Rhie-Chow interplation phi = (U_avg_f & mesh.Sf()) - ( (DUf*( gradp_f - gradp_avg_f)) & mesh.Sf() );

    Time dependency, relaxation of the velocity field:

    // Rhie-Chow interplation phi = (U_avg_f & mesh.Sf()) - ( (DUf*( gradp_f - gradp_avg_f)) & mesh.Sf() ) + (scalar(1) - URFU)*(phi.prevIter() - (U_avg_prevIter_f & mesh.Sf())) + DTf*( phi_old - (U_old_f& mesh.Sf()));

    OpenFOAM
    No time dependency, with or without relaxation of the velocity field (check simpleFoam solver):

    surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(HbyA) & mesh.Sf());

    Questions: Where is the second term of the original RC? Where is the third term of the original RC if taking the relaxation of the velocity field into consideration?

    Time dependency, relaxation of the velocity field (check pisoFoam solver):

    surfaceScalarField phiHbyA ( "phiHbyA", (fvc::interpolate(HbyA) & mesh.Sf()) + fvc::ddtPhiCorr(rAU, U, phi) );

    Questions: ddtPhiCorr takes the function of (the third term+the fourth term of the original RC)? I guess NO. Let's recall the definition of ddtPhiCorr:

    tmp<fluxFieldType> ddtPhiCorr ( new fluxFieldType ( ddtIOobject, rDeltaT * this->fvcDdtPhiCoeff(U.oldTime(), phiAbs.oldTime()) * ( fvc::interpolate(rA*rho.oldTime())*phiAbs.oldTime() - ( fvc::interpolate(rA*rho.oldTime()*U.oldTime()) & mesh().Sf() ) ) ) );

    Look like ddtPhiCorr returns the following codes:

    fvcDdtPhiCoeff*DTf*( phi_old - (U_old_f& mesh.Sf()));

    Last question: How to take the effect of relaxation of the velocity field into ddtPhiCorr?
    谢谢

  • 请教OpenFOAM不同网格之间场映射的原理

    2
    2 帖子
    2k 浏览
    李东岳

    meshToMeshInterpolate.C 可以看看这个内容,对,就是插值

  • interFoam droplet

    4
    4 帖子
    4k 浏览
    同学博

    @东岳 :xinlei:

  • 请问OpenFOAM中如何在入口添加边界层厚度

    15
    15 帖子
    12k 浏览
    疾风GAVIN

    @东岳
    老师这是我目前算出来的一些结果:
    链接: https://pan.baidu.com/s/163fY1uXF_sCCcN9nt1MiCQ
    提取码: w3gp
    网格量很少的,我现在只是想通过调试把路给走通,方便以后的使用

  • 《无痛苦N-S方程笔记》求指导

    3
    3 帖子
    3k 浏览
    Z

    @东岳 谢谢老师,这个问题当时思维定式了,把U_p提出来,立马就懂了。

  • TopoSet关于action的疑问

    2
    2 帖子
    3k 浏览
    B

    @D-Benjamin 求交集的意思

  • 湍流速度脉动

    7
    7 帖子
    7k 浏览
    Y

    @小考拉 你好这个问题我和我师兄讨论过,你不可能输出瞬时的脉动速度场,因为湍流模型根本就没有求解脉动速度,只能输出二阶矩这种统计值,建模的话可以假设三个方向的速度大小相同,给出脉动速度的大小,但是不能给出速度的方向

  • 关于negSumDiag()的一个问题

    11
    11 帖子
    11k 浏览
    G

    @wwzhao

    是的,OpenFOAM中面的方向定义永远是从小编号指向大编号,Gauss定理中面的方向指向外。也许把这些面分为大面和小面区分一下更容易理解。

    对于编号大于$P$的单元$U$,他们之间面上的量为:

    $$\phi_ f = \varpi \phi_{P} + (1 - \varpi) \phi_{U}$$

    $$\dot { m } _ { f } = ( \rho \mathbf { u } ) _ { f } \cdot \mathbf { S } _ { f }$$

    对于编号小于$P$的单元$L$,他们之间面上的量为:

    $$\phi_ f = (1 - \varpi) \phi_{P} + \varpi\phi_{L}$$

    $$\dot { m } _ { f } =- ( \rho \mathbf { u } ) _ { f } \cdot \mathbf { S } _ { f }$$

    $$
    \sum _ { f \sim n b ( P ) } ( \rho \mathbf { u } \phi ) _ { f } \cdot \mathbf { S } _ { f }
    = \sum _ { N \in L ( P ) } \left( - \dot { m } _ { f } \phi _ { f } \right) + \sum _ { N \in U ( P ) } \dot { m } _ { f } \phi _ { f }
    $$
    $$
    = \sum _ { N \in L ( P ) } - \dot { m } _ { f } \left[ ( 1 - \varpi ) \phi _ { P } + \varpi \phi _ { N } \right] + \sum _ { N \in U ( P ) } \dot { m } _ { f } \left[ \varpi \phi _ { P } + ( 1 - \varpi ) \phi _ { N } \right]
    $$

    $$
    = \left( \sum _ { N \in L ( P) } - \dot { m } _ { f } ( 1 - \varpi ) + \sum _ { N \in U ( P ) } \dot { m } _ { f } \varpi \right) \phi _ { P } + \sum _ { N \in L ( P ) } - \dot { m } _ { f } \varpi \phi _ { N } + \sum _ { N \in U ( P ) } \dot { m } _ { f } ( 1 - \varpi ) \phi _ { N }
    $$

    这样可以看出lower存的是 $- \varpi _ { f } \dot { m } _ { f } $, upper存的是$\dot { m } _ { f } ( 1 - \varpi ) _ { f } $.并且diag存的是
    $$a _ { P } = -[\sum _ { N \in L ( P ) } \dot { m } _ { f } ( 1 - \varpi ) + \sum _ { N\in U ( P ) } (-\dot { m } _ { f } \varpi )]$$

  • 有人用过of里的fft吗

    2
    2 帖子
    2k 浏览
    W

    写个类调用下fftw的函数接口就好了