Skip to content

OpenFOAM

OpenFOAM交流区

5.2k 主题 30.9k 帖子
  • 求助:结构入水问题

    8
    8 帖子
    8k 浏览
    Mrc.YoungM

    @enthusiasm qq168167704 加下好友交流一下

  • 关于paraview的问题(Ubuntu17.04)

    8
    8 帖子
    12k 浏览
    bestucanB

    应该是正在写“4.9”,心里打着草稿后面要写安装5,然后5就溜出来。。。一边聊天一边写字很容易把说的听的话写出来:confused:

  • 虚拟机这个边框有办法去掉么?

    3
    3 帖子
    3k 浏览
    李东岳

    我刚换virtualBox,这几天Vmware总是有问题,我这个是4K显示器,感觉水土不服。我试试virtualBox怎么样。

  • Ubuntu 关掉特效

    1
    1 帖子
    1k 浏览
    李东岳
    sudo apt-get install compizconfig-settings-manager compiz-plugins-extra

    然后运行。

  • OpenFOAM中为何要这样定义雷诺应力模型?

    8
    8 帖子
    9k 浏览
    zym604Z

    @wwzhao 谢谢。这个问题我现在已经搞清楚了,是为了增加矩阵的“对角占优”特性,提高收敛性。经过实验,假如修改成不加nu_t的公式,对于一些Co数较大的case就有可能计算发散。

  • 动态铺层网格遇到的问题

    1
    1 帖子
    2k 浏览
    昝磊

    大家好,我在用动态铺层的网格做模拟,模拟的物理过程为:
    在一个用来冷却的容器中,每单位时间内有固定质量且固定温度(1640)的液体进入该容器中(认为不带入动量)并进行冷却。简单的模型如下:
    0_1517960921550_图片1.png

    为了实现这一过程,我使用了动态铺层的网格模拟液体的增加,并且当有新增网格出现时,将新增网格的internalFieldValue进行人为的定义:
    0_1517961307920_1517960362(1).png
    从而避免新增网格的值是周围网格值差值计算得来,实现每单位时间内有固定质量且固定温度液体进入的物理过程。
    但是,我计算出来的结果显示:流场内的最高温度超过了1640度,并且随着网格数目的增加,最高温度越来越高,我做的是冷却过程啊!
    哪位大神知道我这是出什么问题了吗?

  • internalField() 和 Ref()

    3
    3 帖子
    4k 浏览
    W

    这个改动主要是为了区别 const 和 non-const 访问。

  • 中间变量的输出

    3
    3 帖子
    3k 浏览
    C

    @李东岳 谢谢 看来得一步步挖到底层了

  • of中的多相流求解器

    34
    34 帖子
    41k 浏览
    李东岳

    如果单纯的mules求解每一相是无法保证alpha1+alpha2=1的

    现在我就想确定,alpha2通过对应方程推导,用MULES求解不可行的吗?

    MULES的出发点是保证变量的有界。所以如果你用MULES求解alpha1,那么理论上alpha1是有界的。alpha2的求解没有调用MULES,他是通过1-alpha1算出来的,如果alpha1有界,alpha2也有界。

    如果你用算法同时求解alpha1和alpha2,如何处理耦合?这俩个变量是耦合在一起的,就像速度和压力。你不能单独的去分离求解,如果分离求解就就需要迭代。迭代就导致计算速度变慢。因此现存大厂据我所知都是只求解alpha1,然后alpha2=1-alpha1。

    我在去你年底验证了MULES,理论上是可以保证有界,但是真实计算的时候,还是会越界。尤其是物理模型比较复杂的时候。个人觉得这方面内容搞出来,绝对是个好文。

  • icoFoam的公式推导是不是有点问题

    4
    4 帖子
    4k 浏览
    李东岳

    邱老师现在玩的大项目,都是上亿网格的工业级别 :cheeky:

  • 2 帖子
    3k 浏览
    李东岳

    OpenFOAM对于这种切应力有没有什么计算的方法。

    剪切应力的定义是固定的,用什么软件的方程都是一样的。你要从流场中获得$\tau$?印象中OpenFOAM里面没有计算tau的现成的工具。如果你不确定,你可以和Fluent比较其他量看看。

    如果你要计算tau,你可以写个后处理小工具。比如这个

    bool Foam::functionObjects::Lambda2::calc() { if (foundObject<volVectorField>(fieldName_)) { const volVectorField& U = lookupObject<volVectorField>(fieldName_); const tmp<volTensorField> tgradU(fvc::grad(U)); const volTensorField& gradU = tgradU(); const volTensorField SSplusWW ( (symm(gradU) & symm(gradU)) + (skew(gradU) & skew(gradU)) ); return store ( resultName_, -eigenValues(SSplusWW)().component(vector::Y) ); } else { return false; } }

    https://coding.net/u/dyfluid/p/OpenFOAM-5.x/git/blob/master/src/functionObjects/field/Lambda2/Lambda2.C

  • 表面张力与温度的函数

    4
    4 帖子
    4k 浏览
    李东岳

    - 这个操作符重载后你进行的操作不匹配-操作符的参数,你可以把你的代码贴在这里,我看看。

  • 国内coding上的OpenFOAM

    1
    1 帖子
    2k 浏览
    李东岳

    https://coding.net/u/dyfluid/project

    没啥特殊的,就是把github上面的OpenFOAM拽过来了。国内访问比较快。有时候紧急看一下代码,打开网站随便看一下就行了。不用开linux

    注册用户可以创建2个私有项目,我用了半个多月,还挺满意。速度非常不错。主要是github私有项目要收费。感兴趣的可以试试这个平台。不了解的可以百度一下版本控制

  • OF中如何计算得到壁面摩擦应力呢?

    9
    9 帖子
    9k 浏览
    A

    @aeronastro 您好,使用wallgradu求出的结果是个张量,要提取法向的乘以miu吗?还有一个问题请教您,wallgradu的源代码在哪里?

  • 一个不得不承认的事实

    10
    10 帖子
    10k 浏览

    @李东岳 对,我做喷雾燃烧,喷雾用拉格朗日,所以计算量不大,主要的计算量就在化学反应矩阵的求解上。很多都在在化学加速的算法,比如最近一年OF更新的TDAC,就是一个加速的方法,用于几千个反应的化学动力学机理。

  • Ubuntu 16 把工作栏放在左边

    7
    7 帖子
    6k 浏览
    X

    @bestucan 多谢指点,已经按这个方法弄好了!
    @李东岳 也感谢东岳

  • 圆管内固液两相流动

    5
    5 帖子
    5k 浏览
    J

    找到了浮力项了

  • interFoam怎么嵌入到icoFoam中使用

    4
    4 帖子
    4k 浏览
    李东岳

    哦,这个,OpenFOAM编程,我计划加在OKS课程中。那个pdf可能不会写了,当时读博士时间很多,现在时间很紧 :big_mouth:

  • paraview中怎么画虚线

    1
    1 帖子
    2k 浏览
    A

    paraview中contour中能画虚线么,求大神指导

  • 请问thermoType中的Hf是什么含义?

    5
    5 帖子
    5k 浏览
    Y

    谢谢啦。