群组信息 私有

OpenFOAM讲师

被批准进入本群组,则表示群组成员较为理解OpenFOAM的算法和编程。OpenFOAM教授、副教授、讲师岗位授予比例不超过OpenFOAM总徽章中的10%,30%,60%。

 

成员列表

  • 有人尝试过在OF中实现MMIM方法么?

    如题,MMIM即modified momentum interpolation method

    原始文献:DISCUSSION ON MOMENTUM INTERPOLATION METHOD FOR COLLOCATED GRIDS OF INCOMPRESSIBLE FLOW

    OF中的动量插值方法为OMIM,其中还加入了一项ddtCorr,该项可以消除时间步的影响。

    但是在MMIM文献的分析中,表示ddtCorr并不能消除时间步长的影响。而且这个ddtCorr还会带来额外的耗散,如果做LES模拟会降低精度。

    最近我尝试实现MMIM,按着文章中所说,发现实现不了。问题主要是更新phi和U的时候二者相差甚远,会发散

    有没有有兴趣的同学一起来搞一下?

    发布在 OpenFOAM
  • RE: 关于bubble flow中的通量计算问题

    @东岳 如此说来,除了加密网格,没有好的解决方案啦:zoule: :zoule:

    发布在 Algorithm
  • RE: 关于bubble flow中的通量计算问题

    @东岳 如果面上的alpha插值不当,就没法同时保证体积守恒和质量守恒了,我看这对带相变问题计算精度影响很大:zoule: :zoule:

    发布在 Algorithm
  • 关于bubble flow中的通量计算问题

    在网格分辨率不足的情况下研究气泡流,一般会用Mixture or TFM model。计算相方程时,将体积分数插值到面上,选择满足CBC条件的scheme。
    这么插出来的通量是不对的,实际情况下,面上的体积分数不一定非要介于相邻两个cell,这天然不满足CBC。有没有处理这种问题的特殊算法呢?

    发布在 Algorithm
  • RE: share with you

    为社区做贡献😄

    发布在 OpenFOAM
  • share with you 发布在 OpenFOAM
  • RE: FLuent使用两方程模型如何统计流场的雷诺应力分布

    只要得到了速度场随时间变化的数据,你自己就可以求出雷诺应力。

    雷诺分解就是把瞬时速度重写成两项:时均速度与速度的波动(瞬时速度减去时均速度)。用公式写出来就是:
    \begin{equation}
    u = \bar{u} + u^{\prime}, \
    v = \bar{v} + v^{\prime}, \
    w = \bar{w} + w^{\prime}.
    \end{equation}

    把雷诺应力的定义及其展开形式:
    $$ \tau_{ij}^{\prime} \equiv \rho
    \begin{bmatrix}\overline{u_x^{\prime}u_x^{\prime}}&\overline{u_x^{\prime}u_y^{\prime}} &\overline{u_x^{\prime}u_z^{\prime}} \\ \overline{u_y^{\prime}u_x^{\prime}}&\overline{u_y^{\prime}u_y^{\prime}} &\overline{u_y^{\prime}u_z^{\prime}} \\ \overline{u_z^{\prime}u_x^{\prime}}&\overline{u_z^{\prime}u_y^{\prime}} &\overline{u_z^{\prime}u_z^{\prime}}
    \end{bmatrix},
    $$

    其中 (举个例子),$ \overline{u_x^{\prime}u_y^{\prime}} $的意思是 $u_x^{\prime}u_y^{\prime}$ 的时均值。

    具体的算法 ($ \overline{u_x^{\prime}u_y^{\prime}} $) 就是:

    1. 计算 $\overline{u_x}$
    2. 计算 $u_x^{\prime}$
    3. 计算 $\overline{u_y}$
    4. 计算 $u_y^{\prime}$
    5. 计算每个时间点的 $ {u_x^{\prime}u_y^{\prime}} $
    6. 计算 $ \overline{u_x^{\prime}u_y^{\prime}} $

    将最后的结果乘以流体的密度,就得到了雷诺应力。 如果有错,还请大家指正。

    0_1537558530579_ReynoldStress.png

    发布在 ANSYS Fluent
  • RE: 懂网络的专家进来看看

    @东岳 意思是要中转N多次?为什么不能用远程桌面?Linux上没有X?

    发布在 C斯达克
  • RE: 懂网络的专家进来看看

    采用NAT方式网段一般不同,需要在虚拟机上做端口转发。此外,Windows通过路由器上网,那么路由器上也需要做端口转发。

    这么做比较麻烦,直接用TeamViewer吧。

    发布在 C斯达克