Skip to content

OpenFOAM

OpenFOAM交流区

4.8k 主题 28.5k 帖子
  • rhoReactingFoam求解化学反应

    1
    1 帖子
    2k 浏览
    C

    请问有谁用rhoReactingFoam算过变热物性的例子,能否给提供的例子作为参考?
    热物理模型的形式是
    0_1493025160357_upload-bf30d0c9-ae7c-4d4f-847b-ce5438bb68c2

  • 如何得到processor 的编号?

    7
    7 帖子
    5k 浏览
    qjh888Q

    @小葱小虫
    是第三方的,要使用foam-extend才能安装。
    transonicMRFDyMFoam可以计算稳态

  • 1 帖子
    2k 浏览

    如题,貌似OF的BC都是施加在原始非守恒变量上的,怎么整成特征变量呢?

  • 管道流动

    4
    4 帖子
    4k 浏览

    @zhanghan 你的压力场进出口不建议设置为全部为zeroGradient。可以把出口设置为固定值。

  • 1 帖子
    2k 浏览
    C

    上次提了问题,就是发现利用不可压二相动网格求解器计算波浪冲击直墙时,压力出现震荡,地址如下:
    http://www.cfd-china.com/topic/784/求教关于计算开始时压力震荡的问题/5

    多谢@李东岳 和 @wwzhao 两位版主的回复,我回去又研究了一下,发现以前有不少算例都出现震荡情况,只不过后来的压力变化相对于之前的震荡值量级比较大,因此掩盖了这个问题,但是问题依旧存在,所以我加了一个直墙,设置的边界条件与右边边界相同:p_rgh中为 fixedFluxPressure ,U中为fixedValue,大小都为0.

    再交代一下,我采用动网格造波,即通过边界的运动模拟实现类似于推板造波机的功能。之前造波效果都不错,直到最近发现计算压力存在震荡现象,这次我又算了几个例子,首先我的计算区域是这样的:
    算例长21m,设置一个直墙,位于x=20m。右端设置为直墙,左端动边界,距离动边界越远,网格运动的幅度越小。
    0_1492567693235_0s-chushi.jpg
    初始速度场:
    0_1492567614566_u-0s.jpg
    计算开始后,我发现速度场在之直墙两端有突变:
    t=0.5s,整体速度场:
    0_1492567862778_u-0.5s.jpg
    直墙附近压力震荡:
    0_1492567891321_u-0.5s-jubu.jpg
    t=1s,整体速度场:
    0_1492567920059_u-1s.jpg
    直墙附近压力震荡:
    0_1492567927302_u-1s-jubu.jpg

    基本上计算开始几个时间步后,就会在直墙的水、气交界面开始出现这种问题,最后就会在得到震荡的压力。
    我还特意试了一下wave2foam速度入口造波,发现没有这个问题。虽然我知道不可压缩流体的速度和压力瞬间更新全场,但是这个地方怎么这么诡异?还是动网格导致的呢?如何能够优化或者避免这种情况?

    请大牛们不吝赐教,多谢了:)

  • mapFields 太慢,想在里面加入cputime这样的输出

    4
    4 帖子
    4k 浏览
    H

    @李东岳
    所以每一个map时间步后我把所有的0文件重命名成相应的时间步再进行下一步。
    下面是相关的附件:

    我用的脚本:
    0_1492507241431_Allrun1.txt
    mapFieldsDict:
    0_1492507250371_mapFieldsDict.txt

    Source case 转子部分
    0_1492509534847_SourceMesh_with_blades_rotwall_blade.jpg
    还有cyclicAMI类型的rotinterface_rmax _zmax _zmin没有显示出来

    Target case 转子部分,cuttingPatches的两个patch
    0_1492508448117_TargetMesh_no_blades_cuttingPatches.jpg

  • 关于reactingFoam计算化学反应问题

    1
    1 帖子
    2k 浏览
    C

    利用reactingFoam计算正癸烷在水平圆管中的裂解反应(18组分生成物)、热物理模型采用
    0_1492507996639_upload-7ce04876-1271-4e49-b10f-e26cf46d3989
    计算结果中出口速度偏小、正癸烷在圆管中的停留时间过长,导致出口处的正癸烷质量分数偏小,考虑到与其密度有关;因为正癸烷裂解反应是在超临界条件下进行的,算例中采用的热物理模型不能满足计算要求,故考虑采用rhoReactingFoam求解器,热物理模型采用
    0_1492508023620_upload-b90bd2f7-b8c7-4863-8558-4089770b0ad8
    各个热物性参数通过与温度之间的多项式关系进行计算。例:
    0_1492507617167_upload-2786160a-cf82-4436-b6a6-f40fb8d70e2f
    但计算一步就出现错误
    0_1492507781160_upload-d1e18a85-91ef-448b-a634-cacd5cda8fc0
    目前不知道什么原因导致的,请问哪位大神能给指点一下?

  • 采用reactingFoam求解化学反应

    1
    1 帖子
    2k 浏览
    C

    利用reactingFoam计算正癸烷在水平圆管中的裂解反应(18组分生成物)、热物理模型采用0_1492507194722_upload-7d5ad516-4737-4153-b8ec-01e682a983c3 ,计算结果中出口速度偏小、正癸烷在圆管中的停留时间过长,导致出口处的正癸烷质量分数偏小,考虑到与其密度有关;因为正癸烷裂解反应是在超临界条件下进行的,算例中采用的热物理模型不能满足计算要求,故考虑采用rhoReactingFoam求解器,热物理模型采用0_1492507500590_upload-d139e8d7-727d-4acd-8d0b-71cc7303f484
    各个热物性参数通过与温度之间的多项式关系进行计算。例:
    0_1492507617167_upload-2786160a-cf82-4436-b6a6-f40fb8d70e2f
    但计算一步就出现错误
    0_1492507781160_upload-d1e18a85-91ef-448b-a634-cacd5cda8fc0
    目前不知道什么原因导致的,请问哪位大神能给指点一下?

  • 后处理:tecplot360 画矢量图问题

    1
    1 帖子
    2k 浏览
    H

    在用tecplot360画速度矢量图时,箭头分布比较密,想把它调稀疏点,调point只能改变I方向上的点数,J,K方向的无法调,请问有遇到类似的吗?

  • OpenFOAM + Cantera

    7
    7 帖子
    8k 浏览
    Z

    楼主将cantera移植到最新OpenFOAM的努力成功了吗?还望分享

  • Trailing edge 绕流产生checkerboarding

    4
    4 帖子
    4k 浏览
    qjh888Q

    或者说是,在使用backwardDual physDeltaT CoDeltaT这个ddtScheme的时候,如何去确定我想要的deltaT 和求解的subCycle 的数量呢?

  • 求教关于计算开始时压力震荡的问题

    5
    5 帖子
    5k 浏览
    C

    @wwzhao 这个是我自己写的边界,没写别的造波方式。。。。

  • 用simpleFoam求解器时出现的错误,求教如何解决

    13
    13 帖子
    13k 浏览
    C

    @温宝霖 恭喜。
    你模拟的是什么问题?能否发个网格截图出来看看?听起来挺有意思。

  • 如何模拟移动的内热源?

    4
    4 帖子
    4k 浏览
    C

    @piteqiu 使用laplacianFoam,添加#include "fvOptions.H",方程中添加源项(参考scalarTransportFoam),在fvOptions中定义codedSource。

    src/fvOptions/sources/general/codedSource/CodedSource.H

    你需要依据时间确定你的区域添加源项就可以了。可能需要一些网格的东西:mesh.C()返回网格位置,mesh().time()返回当前时间。

  • 关于SpalartAllmaras模型

    4
    4 帖子
    4k 浏览
    H

    OF有标准版,什么都不需要改,如果仔细看程序的话。
    Wall shear在有fv3的情况下不是特别好而已。

  • paview怎么计算管道截面的平均温度

    3
    3 帖子
    4k 浏览
    G

    在https://www.cfd-online.com/Forums/openfoam-paraview/122425-error-calculators-equation.html#post447832
    中有提到用Paraview计算平均温度的方法,

    0_1491533074551_upload-3aaf8874-b0d3-40b7-822c-db8ddb6b5aa6
    不知道这种面积平均得到的温度是不是就是平均温度。

  • 关于边界条件

    6
    6 帖子
    6k 浏览
    L

    @xpqiu 边界条件刚上来的那个图是什么意思呢?

  • 3 帖子
    3k 浏览
    Y

    @xpqiu :big_mouth:

  • 2 帖子
    2k 浏览

    发散的问题很常见,并不一定要换求解器。跨音速用密度基确实合理。提供更多信息看看是如何发三的。

  • openfoam初场插值的问题

    2
    2 帖子
    2k 浏览

    可以,你试试mapFields,具体如何使用你看OpenFOAM的用户指南里面有使用方法,icoFoam下面有个算例可以专门练习使用mapFields。你先试试,不成功再联系!