Skip to content

OpenFOAM

OpenFOAM交流区

5.3k 主题 31.4k 帖子
  • sensibleEnthalpy与 sensibleInternalEnergy

    1
    1 帖子
    2k 浏览
    羽之下

    thermophysicalProperties中energy可以选择sensibleEnthalpy, sensibleInternalEnergy 和 absoluteEnthalpy。前两者与后者的区别,后者使用绝对能量表示包含热源的情况。而sensibleEnthalpy与sensibleInternalEnergy的区别,前者是利用焓值求解温度,后者是利用内能求解温度。热⼒学中分析闭式系统⽤内能,分析开式系统⽤焓值。那么在实际使用中,该如何选择这两种模式?例如高压离心风机。tutorials中的案例似乎找不出规律。
    eb324645-1dd5-499a-bfc0-6d77667f7fa5-image.png

  • paraview如何获得三维流线图

    1
    1 帖子
    2k 浏览
    D

    大家好,请问paraview如何获得三维流线图?如下图这样的,用户指南上关于流线图的说明仅仅是寥寥几笔,请各位不吝赐教,谢谢!
    feaadcd8-3566-4e3f-831c-5e9a11ca59e0-图片.png

  • 动网格求解中correctPhi的作用是什么

    5
    5 帖子
    5k 浏览
    李东岳

    是的,一般情况下你1)可以忽略correctPhiH,2)pcorr也可以不求解

  • 2 帖子
    2k 浏览
    李东岳

    在代码里加入

    volScalarField rho ( IOobject ( "rho", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh );

    然后直接用就行

  • 4 帖子
    4k 浏览
    李东岳

    @zp气固两相流 只有考虑面的时候才有一个own和一个nei,考虑cell的时候一个own,其他都是nei,所以考虑cell的时候,一个own,很多nei很正常啊

  • 瞬态问题残差一直保持水平的疑问

    1
    1 帖子
    1k 浏览
    C

    我正在做加热管道空气流动例子的DNS,有新旧两个算例,其中旧算例的残差情况是:

    初始.png

    然后我把中间结果文件copy到新算例里面,接着继续跑,然后现在残差一直都比较水平:

    中间.png

    但是量级还是有点高,看了平均壁温,也没有达到均匀和期望效果。

    我目前用的是基于buoyantPimpleFoam改编的求解器,设置成piso模式来跑,内迭代两次。U跟h的离散格式是QUICK.

    不知道出现这种情况暗示着什么?有什么办法可以改变这种残差一直水平的情况呢?谢谢!

  • 求OpenFOAM中MULES算法的解析

    4
    4 帖子
    5k 浏览
    李东岳

    @史浩 他那个应该是个笔误。写反了。

    我这周刚把原始的FCT算法植入进去了,OpenFOAM的算法跟原始的不太一样,主要在于

    最大值和最小值的确定,这对应不同的方法,参考Zalesak1979 342页末

    反扩散系数的确定(你贴的按个公式)

    这也是他之所以叫做MULES的原因吧,不过思路都是一样的,Weller为了加快计算速度也优化的非常好。原始FCT已经完全可以对变量进行有界了,我对比原始的FCT和MULES差异非常小,原始FCT也不需要对$\lambda$进行迭代。我要用原始FCT验证几个算例分析一下。

    在测试可压缩传输的时候,我发现可压缩FCT和不可压缩FCT的密度处理需要特殊处理。也就是说 http://www.cfd-china.com/topic/1029/mules算法中的const-rhotype-rho/3 帖子中的可压缩不能直接这么处理,要把rho和alpha分开。不能看做对rhoAlpha的传输。具体的我还在研究

    下图FCT1 和FCT2对应原始FCT和MULES,差别微小。抽取的damBreak的alpha的最小值

    捕获.JPG

  • MULES::explicitSolve求解相方程

    2
    2 帖子
    3k 浏览
    李东岳

    phi是用来调用判断速度方向的: surfaceScalarField phiBD(upwind<scalar>(psi.mesh(), phi).flux(psi));

  • 并行运算出现错误

    4
    4 帖子
    2k 浏览

    一般什么原因比较多

  • 5 帖子
    3k 浏览
    B

    @东岳 老师人太好了!再次感谢!

  • 关于自定义边界条件

    3
    3 帖子
    3k 浏览

    @adcpk 在 关于自定义边界条件 中说:

    codestream

    多谢

  • 基类中用克隆代替复制

    1
    1 帖子
    1k 浏览
    马乔

    在构建类层次结构是,我们通常都是想要获得多态行为。我们需要对象一旦创建就能保持其类型和标识。但是在想复制一个对象时会发现C++中的对象复制语义与这一目标是冲突的,因为构造函数都是与类名相同的,也就是构造函数都是每个类特有的,复制构造函数不能设为虚拟的。
    为了用基类指针(或引用)正确的将对象进行复制,更为通用的惯用法解决方案是,将基类的复制构造函数设为protected甚至是public,并改而依赖于一个虚拟的clone函数。例如particle类中

    //- Construct as a copy particle(const particle& p); //- Construct a clone virtual autoPtr<particle> clone() const { return autoPtr<particle>::New(*this); } //Constructor Foam::particle::particle(const particle& p) : mesh_(p.mesh_), coordinates_(p.coordinates_), celli_(p.celli_), tetFacei_(p.tetFacei_), tetPti_(p.tetPti_), facei_(p.facei_), stepFraction_(p.stepFraction_), origProc_(p.origProc_), origId_(p.origId_) {} 在派生类KinematicParcel中 //- Construct as a copy KinematicParcel(const KinematicParcel& p); //- Construct and return a (basic particle) clone virtual autoPtr<particle> clone() const { return autoPtr<particle>(new KinematicParcel(*this)); }

    通过虚拟函数clone就能正确的实现复制功能啦!

  • 对于数字一定要谨慎!!!!

    3
    3 帖子
    3k 浏览
    Y

    不同的语言在这块不太一样,估计是写惯了其它语言,转到c++就惯性思维了:duang:

  • 非稳态计算向压力方程中添加源项

    12
    12 帖子
    10k 浏览
    C

    @Tong 好的,谢谢啊!

  • 请问Info 命令如何设置特定精度?

    5
    5 帖子
    4k 浏览

    @东岳 我在controldict 试过一次 好像不管用

  • 请问open foam中解方程的方法

    4
    4 帖子
    4k 浏览
    李东岳

    嗯,牛顿迭代法,Ridder迭代法之类

  • 如何设置水体密度随深度线性变化?

    2
    2 帖子
    2k 浏览
    李东岳

    interFoam加不了,compressibleInterFoam你在物性里面修改

    thermoType { type heRhoThermo; mixture pureMixture; properties liquid;//here energy sensibleInternalEnergy; }

    里面的properties改成不同的状态方程试试,

  • 结果无法输出

    2
    2 帖子
    2k 浏览

    请问你的
    adjustTimeStep yes/no?
    maxCo?

  • openfoam中压力梯度计算不均

    3
    3 帖子
    3k 浏览
    Z

    非常感谢李老师,cellLimited Gauss linear 1 格式比Gauss linear好很多,运行程序之前运行renumberMesh -overwrite改善网格质量也有作用

  • OpenFOAM 非对称矩阵解法器

    2
    2 帖子
    2k 浏览
    李东岳

    法器算得总是不太对, 误差有点大.

    从何说起?

    OpenFOAM这里面不会有bug,这个太重要了,所有求解器都要用,如果有bug那结果就全都乱糟糟啦