Skip to content

OpenFOAM

OpenFOAM交流区

5.3k 主题 31.4k 帖子
  • openFOAM自适应网格求助!!

    5
    5 帖子
    1k 浏览

    @李东岳 好的好的,谢谢老师回复。

  • 边界条件设定

    1
    1 帖子
    431 浏览
    A

    各位大佬,我要模拟一个圆柱体,上表面逐步加载,上面表面没有流入和流出。水从下表面流出,算出来的结果很奇怪,要求计算10S,前3S都算正常,水向下流,但算到5S就停止了,最后一步水反向流出了。以下是流速的计算结果,帮忙看看,边界条件设置是否有问题,为什么结果怪怪的,谢谢。

    circle.jpg

    circle2.jpg

    压力边界条件
    inlet
    {
    type fixedValue;
    value uniform 0.01;
    }

    pipe
    {
    type zeroGradient;
    }
    outlet
    {
    type zeroGradient;

    }

    流速边界条件
    inlet
    {
    type zeroGradient;

    }

    pipe
    {
    type noSlip;
    }
    outlet
    {
    type inletOutlet;
    inletValue uniform (0 0 0.1);
    value uniform (0 0 0.1);
    }

  • 如何建一个圆柱体

    14
    14 帖子
    6k 浏览
    A

    我有几个问题啊,inlet上表面压力为什么是固定值啊? 我要是跑着跑着里面的压力变得大于入口压力了,那岂不是反流了?
    此外,速度边界条件的Outlet为什么给uniform (0 0 0.1);这个值啊?

  • DNS求取yplus

    2
    2 帖子
    1k 浏览
    X

    wallShearStress可以计算tau_w。然后用来算y+就好啦。

  • Openfoam安装报错失败

    8
    8 帖子
    2k 浏览
    T

    感谢🙏,虚拟机总共给的4个,CSDN也查到说是并行太多,单J无后缀就是拉满,等会儿用J2重新编

  • fvc::grad(CenPot)为什么算出来不对??

    3
    3 帖子
    831 浏览
    C

    @东方白杨 不清楚什么原因,你看在程序最后再重新计算一遍保存,看看结果有无区别

    GradCenPot = - fvc::grad(CenPot); GradCenPot.write();
  • postProcess和并行冲突

    4
    4 帖子
    921 浏览
    李东岳

    这个是自己写的求解器,所以有可能差别确实很大

  • Openfoam1912编译waves2Foam失败,求教

    11
    11 帖子
    2k 浏览
    L

    https://www.cfd-china.com/topic/7301/waves2foam的一键安装解决方案?_=1738765131135

  • 小白入openFOAM

    3
    3 帖子
    740 浏览
    L

    应该是偏算法吧,线下课应该也是以算法为主吧?那我还是学C++吧,到时候看赶不上这个三月底的线下课了

  • 添加第三个区域到chtMultiRegionFoam

    3
    3 帖子
    593 浏览
    W

    @李东岳 感谢大佬的回复!!
    您的意思是我需要创建cellZones来划分这些区域吗,那么关于这些区域的控制方程以及其他属性我需要如何进行分配,之前是在求解器内部进行布置的。

  • blockMeshDict报错相邻面数量不一致怎么解决?

    10
    10 帖子
    2k 浏览
    E

    @李东岳 一个普通的打工研究生罢了。。。

  • limitedLinear格式

    7
    7 帖子
    3k 浏览
    李东岳

    div(phi,U) Gauss limitedLinear 0;

    我还第一次见有人这样用。我下周找个小网格测试一下。

    div(phi,U) Gauss linear;

    另外你这个结果看起来不太对,有震荡

  • GAMG 计算发散

    8
    8 帖子
    2k 浏览
    A

    你网格变密了,库朗数也得跟着改吧,是自适应时间步长吗? 可以考虑把库朗数限制在0.5以下?

  • 以温度表示的能量方程到底是k/Cp还是k/Cv?

    3
    3 帖子
    797 浏览
    东方白杨

    @Amadeus 是的,一个是Cp一个是Cv,但是目前看的各种版本推导的温度方程,有的用Cp有的用Cv不太一致。。。。

  • OpenFOAM SRF求解器有大问题呀?

    6
    6 帖子
    2k 浏览
    东方白杨

    @Amadeus @aurora16 感谢各位回复,最后发现是由于对流-扩散方程中,离心力项太大了(高达几百量级),而其它项还在10的量级,离散求解方程时,量级差异太大了,导致发散。计算还是要保证方程中各项量级差异不能太大

  • 请教boxTurb生成初始场问题

    4
    4 帖子
    1k 浏览
    李东岳

    不需要potentialFoam

  • rhoSimpleFoam 求解可压缩流动时一个奇怪的现象

    54
    54 帖子
    79k 浏览
    B

    我之前算过高亚声速叶栅,用rhoSimpleFoam求解器收敛不了,换成LUSGS就很好收敛

  • 1 帖子
    335 浏览
    C

    如题,现在已经计算完成了,在paraview中也可以查看云图,怎么绘制最大值随时间变化曲线

  • 如何设置运行期自动变化写数据的时间间隔?

    16
    16 帖子
    15k 浏览
    Prometheus10P

    @浪迹天大 老师您好,按照这个方法,我需要更改controlDict中的maxDeltaT。但是会发现文件更新不起作用。

    /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2012 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application interFlowVb1; startFrom latestTime; startTime 0; stopAt endTime; endTime 0.3; deltaT 0.00005; writeControl adjustableRunTime; writeInterval 0.05; purgeWrite 0; writeFormat binary; writePrecision 8; writeCompression off; timeFormat general; timePrecision 8; runTimeModifiable yes; adjustTimeStep yes; maxCo 0.5; maxAlphaCo 0.5; maxCapillaryNum 10; // stability criteria is 1 maxDeltaT 2e-5;//controlDict.a; // maxDeltaT 1e-5;//controlDict.b; // maxDeltaT 5e-6;//controlDict.c; functions { fileUpdate1 { type timeActivatedFileUpdate; libs ("libutilityFunctionObjects.so"); writeControl timeStep; writeInterval 1; fileToUpdate "$FOAM_CASE/system/controlDict"; timeVsFile ( (-1 "$FOAM_CASE/system/controlDict.a") (1e-4 "$FOAM_CASE/system/controlDict.b") (2e-4 "$FOAM_CASE/system/controlDict.c") ); } }

    三个文件只有maxDeltaT有差别,请老师捞捞@李东岳 @浪迹天大

    运行后结果提示:

    timeActivatedFileUpdate: copying file "~/orifice3Da/system/controlDict.c" to: "~/orifice3Da/system/controlDict" Courant Number mean: 2.0171174e-08 max: 0.3369988 Interface Courant Number mean: 2.978204e-10 max: 0.22052861 Capillary Number: 2.4526254 deltaT = 2e-05 Time = 0.00022

    文件有更新,并且maxDeltaT也更新为了1e-5; 但是它没有发挥作用,可以看到文件更新后的deltaT 依然是 2e-05,如果是我自己运行过程中手动改小maxDeltaT,它会马上从deltaT = 2e-05,便为deltaT = 1e-05。

  • 在interfoam两相流中,如何计算相对渗透率

    6
    6 帖子
    1k 浏览
    1

    @JC 在 在interfoam两相流中,如何计算相对渗透率 中说:

    @JC 在 在interfoam两相流中,如何计算相对渗透率 中说:

    有一个标准里面写有计算方法,记录每个时刻流出的两相体积,以及压力梯度,应该就可以算出来了。对了还需要这个多孔介质的初始渗透率

    《SY/T 5345-1999》一个很老的标准,可以参考下

    谢谢,后面我也是直接用出口流量来算每一相渗透率再计算相对渗透率,部分文献里也存在根据饱和度计算相对渗透率的模型,不过不好确定误差。