Skip to content
  • 0 赞同
    14 帖子
    18k 浏览
    李东岳

    其他你的理解大体正确,·volScalarField和scalarField的区别仅仅在于volScalarField附加了边界的信息,scalarField`只是网格点的信息不包含边界。

  • 0 赞同
    1 帖子
    2k 浏览
    A

    OpenFOAM-3.0.x/src/lagrangian/basic/particle文件夹下particleTemplates.C中line 202——line 530描述了trackToFace函数,即对颗粒在欧拉网格下的追踪过程。对lambdaMin小于0或者大于1的情况出现的原因不是很清楚,这种情况下记录的颗粒所在的网格id号和颗粒真实的位置不一致,对应于源程序中501——506行The particle can be 'outside' the tet. This will yield a lambda larger than 1, or smaller than 0. For values < 0 ...............,程序下面给出了解决方案
    但是究竟是什么原因导致的The particle can be 'outside' the tet?希望大神能够给些建议。

    目前个人对程序理解到的程度是:

    1)在一个拉格朗日时间步长内,需要对颗粒运动路径通过的每一个欧拉网格的id号以及颗粒在那个网格中的停留时间进行识别和计算;

    2)颗粒下一时步所在的网格id号是通过时间逐步推进来获取相邻网格id号的办法实现,这里最重要的是判断颗粒穿透了网格的哪一个面以及那个面对应的相邻网格;

    3)通过颗粒当前所在网格(实际应用的时候是将颗粒所在网格划分为符合一定质量要求的四面体网格,细化颗粒追踪过程)面信息,颗粒当前位置,四面体网格中心点位置,颗粒一个时步内终点位置来计算lambdaMin的值(并记录了此值对应的网格面编号),个人理解程序通过lambdaMin的值是否在0与1之间来判断出颗粒是否“真正”穿过了那个网格面,对于此范围之外的lambdaMin的值,对应于颗粒实际位置在它所记录的网格之外,可以认为是颗粒跟踪失败,程序给出了“补救措施”

    简而言之,lambdaMin小于0或者大于1的情况从程序的角度是如何产生的?如何通过输出有关的变量来进行验证(同时也为了解决这样一个疑惑:颗粒仿真时有时会有颗粒穿透边界而落到计算域的外部,个人感觉和这个过程很相关)

  • 0 赞同
    2 帖子
    4k 浏览
    I

    没有遇到过,您探索一下,然后总结一下岂不是更好!
    功德无量!

  • 0 赞同
    12 帖子
    16k 浏览
    J

    您好,我也在学习用OpenFOAM实现源造波。
    请问质量源造波是只能通过在interFoam的pEqu.H中的压力方程的等号右边加入源项吗?
    我了解到为了方便源项有fvOptions,就是在case里增加一个fvOptions文件,在这个文件里说明源项等相关设置,请问楼主在方面是否有尝试?
    我目前在做源项造jonswap不规则波,但是在引入的源项方面有些不确定,原理是Lin和Liu的1999年的文献中的这一段:
    5c5d8285-847a-4bed-bfed-54cb967a6b2f-image.png
    所以想了解楼主加入源项的源项是怎么表达的?加入源项后是通过在求解器里wmake编译就可以在case里使用源造波功能了吗?
    希望能得到老师和各位前辈们的指导,感激~

  • 基金会那个造波的怎么样?

    OpenFOAM
    6
    0 赞同
    6 帖子
    10k 浏览
    C

    @maoyanjun_dut 我写了一个推板造波的solver,也写了主动消波,阻尼区什么的,目前正在做质量源造波,这样就可以算浮体了

  • CFL3D开源了

    OpenFOAM
    8
    0 赞同
    8 帖子
    15k 浏览

    现在这个V6.7版本有在Windows下编译好了的吗?

  • LES和RANS中,离散的对象是相同的吗?

    Algorithm
    31
    0 赞同
    31 帖子
    59k 浏览
    H

    如果使用standard KE,计算二维流场,在某个高雷诺数下把Cd, Cl, St都算对,那么说明碰对了,意思就是好运。这并不能说明standard KE可以算对漩涡脱离。

  • Foam4.0编译出错请教

    OpenFOAM
    5
    0 赞同
    5 帖子
    7k 浏览
    Z

    @队长别开枪 就是单独编译后,paraview也正常了。不提示有错是不是就算可以了?

  • 关于2D axisymmetric边界条件的设定

    OpenFOAM
    4
    0 赞同
    4 帖子
    6k 浏览
    D

    对称轴可以不设任何边界条件,因为轴上没有网格面。

  • 0 赞同
    9 帖子
    14k 浏览
    B

    @xpqiu 好的,谢谢老师

  • 一阶迎风格式数值扩散的原因

    Algorithm
    6
    0 赞同
    6 帖子
    11k 浏览
    H

    0_1501229627502_figure1.PNG

  • ANSYS Fluent理论指南开源文档

    Fluent
    6
    0 赞同
    6 帖子
    12k 浏览
    A

    不知道现在是否发布了理论指南的开源文档?

  • 周期网格角部处理问题

    Meshy
    6
    0 赞同
    6 帖子
    11k 浏览
    piteqiuP

    我这个是一个圆柱,分成1/26,成为周期对称模型。中间是空气,不能省略。

  • swak4Foam 和 OpenFOAM4.1

    OpenFOAM
    3
    0 赞同
    3 帖子
    6k 浏览
    X

    多谢指点!

  • LES介绍的文章的一个公式

    OpenFOAM
    14
    0 赞同
    14 帖子
    21k 浏览
    G

    @一二 在 LES介绍的文章的一个公式 中说:

    嗯,下面是我推的
    \begin{equation}
    \overline{\mathbf{S}}:\mathbf{\tau}+C_{e}\frac{k_{sgs}^{1.5}}{\Delta}=0 \
    left = \overline{\mathbf{S}}:[\frac{2}{3}k_{sgs}\mathbf{I}-2\mu_{sgs}\mathbf{dev}(\overline{\mathbf{S}})]+k_{sgs}^{1.5}\frac{C_{\epsilon}}{\Delta}
    \end{equation}
    因为
    \begin{equation}
    \mu_{sgs}=C_{k}\Delta\sqrt{k_{sgs}} \
    = \overline{\mathbf{S}}:[\frac{2}{3}k_{sgs}\mathbf{I}-2C_{k}\Delta\sqrt{k_{sgs}}]+k_{sgs}^{1.5}\frac{C_{\epsilon}}{\Delta} \
    = \sqrt{k_{sgs}}(\frac{C_{\epsilon}}{\Delta}k_{sgs}+\frac{2}{3}\mathbf{tr}(\overline{\mathbf{S}})\sqrt{k_{sgs}}-2C_{k}(\mathbf{dev}(\overline{\mathbf{S}}):\overline{\mathbf{S}})) \
    = ak_{sgs}+b\sqrt{k_{sgs}}-c\
    = right = 0
    \end{equation}
    其中
    \begin{equation}
    a=\frac{C_{\epsilon}}{\Delta} \
    b=\frac{2}{3}tr(\overline{\mathbf{S}}) \
    c=2 \Delta C_{k}(\mathbf{dev}(\overline{\mathbf{S}}):\overline{\mathbf{S}} \
    \sqrt{k_{sgs}}=\frac{-b+\sqrt{b^2+4ac}}{2a}
    \end{equation}
    当为不可压缩流体时$tr{\overline{\mathbf{S}}}=0$,那么$b=0$、$c=2 \Delta C_{k} (\mathbf{dev}(\overline{\mathbf{S}}):\overline{\mathbf{S}})$,就可以得到$k_{sgs}=\frac{c}{a}=\frac{2C_{k}||\overline{\mathbf{S}}||^{2}{}\Delta^2}{C_{\epsilon}}$

  • 有关压力速度耦合算法的文章

    Algorithm
    2
    0 赞同
    2 帖子
    4k 浏览
    W

    @李东岳 coupled solver比segregated solver能更快收敛,但在处理复杂几何形状物体时不太容易收敛[1]。

    [1] http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2012/KlasJareteg/KlasJareteg-OF2012-Project-3-Version-2.0.pdf

  • snappyhexmesh生成的网格质量不好

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

    @搬运工不好当 shm生成的网格有的时候非常的tricky,不错如果你采用shm成功的生成了mesh okay的网格,openfoam大部分求解器都能算收敛。of里面单相求解器已经非常稳健了。

    如果对于现有几何,一定要生成一个非常好的网格,我建议你

    简化几何 换专业的网格生成程序 慢慢的,充分的调节shm里面的网格控制参数,主要是网格密度,迭代次数,甚至删除边界层,这个步骤需要很长时间。

    话说回来,shm对于复杂几何生成的是混合网格。虽然单相求解器对混合网格非常稳健,但某些多相求解器还算不了。

  • FLUENT ERROR

    Fluent
    3
    0 赞同
    3 帖子
    5k 浏览
    书生kao

    并行计算的问题,一般是计算过程出错就会出现这些error文件。
    删掉没啥影响,重新打开并行计算就可以了。

  • of-v1706中的IsoAdvector

    OpenFOAM
    25
    0 赞同
    25 帖子
    37k 浏览
    李东岳

    我也忘了..

  • 按照随机顺序对计算结果进行采样

    OpenFOAM
    7
    0 赞同
    7 帖子
    11k 浏览
    WayneW

    @程迪 厉害。确实有一些计算机方面的方法可以搞,不过用vtk解决了,多谢!

京ICP备15017992号-2