OpenFOAM

4.5k 主题 26.4k 帖子
  • openFoam与Fluent模拟简单的Laval蒸汽喷管

    21 帖子
    3k 浏览

    @李东岳 东岳老师,这个用rhoCentralFoam的case可以发上来参考一下吗

  • 8 帖子
    261 浏览

    @wsy11 好的,稍等,我提交了

  • 燃烧求解在AMR之后报错

    5 帖子
    219 浏览

    破案了。程序没问题,是我自定义的用于AMR的coded量忘记clear了

  • OpenFOAM outflow边界条件

    6 帖子
    1k 浏览

    @shrine 是只在出口下降了吗,你把出口压强改成零梯度试一下,可以看一下案例里面那个open channel的设置

  • 边界条件tableFile读取错误问题

    8 帖子
    3k 浏览

    @Hope 楼主有没有ESDU的document呀,风剖速度和湍流强度以及风谱,没订阅,也搜索不到

  • toposet 得到的新边界怎么加边界层

    1 帖子
    77 浏览

    各位前辈好,请问我使用toposet创建的新的边界或者面之后,我应该怎么在这个边界或者说面上加上边界层呢,是不是相加边界层只能用snappyHexMesh

  • 3 帖子
    163 浏览

    @学流体的小明 恩呢,谢谢您的回答。
    /--------------------------------- C++ -----------------------------------
    | ========= | |
    | \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
    | \ / O peration | Version: v2212 |
    | \ / A nd | Website: www.openfoam.com |
    | \/ M anipulation | |
    *---------------------------------------------------------------------------*/
    FoamFile
    {
    version 2.0;
    format ascii;
    class dictionary;
    location "constant";
    object dynamicMeshDict;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

    dynamicFvMesh dynamicRefineFvMesh;

    // How often to refine
    refineInterval 1;

    // Field to be refinement on
    field alpha.water;

    // Refine field inbetween lower..upper
    lowerRefineLevel 0.001;
    upperRefineLevel 0.999;

    // If value < unrefineLevel unrefine
    unrefineLevel 1;

    // Have slower than 2:1 refinement
    nBufferLayers 6;

    // Refine cells only up to maxRefinement levels
    maxRefinement 1;

    // Stop refinement if maxCells reached
    maxCells 15000000;

    // Flux field and corresponding velocity field. Fluxes on changed
    // faces get recalculated by interpolating the velocity. Use 'none'
    // on surfaceScalarFields that do not need to be reinterpolated.
    correctFluxes
    (
    (phi none)
    (nHatf none)
    (rhoPhi none)
    (alphaPhi_ none)
    (ghf none)
    (phi0 none)
    (dVf_ none)
    (alphaPhi0.water none)
    (alphaPhiUn none)
    );

    // Write the refinement level as a volScalarField
    dumpLevel true;

    // ************************************************************************* //
    这是我的dynamicMeshDict文件。
    我按您说的方法去试一下。

  • 6 帖子
    333 浏览

    @Prometheus10
    我看了一下那个帖子,可能因为实现方法的差异,两种边界会有细微的差别,但应该不会这么大。

    你的其他边界是怎么设置的?symmetry不用特别设置,每一个变量都是sym就行。但是,slip你是具体怎么设置的?速度肯定是slip,其他边界是什么?如果搭配不合适的话,和sym不一致是很正常的事情。

    再一个,是不是因为雷诺数特别小,以致实现方法带来的微小区别被放大了?或者说,在低雷诺数下,实现方法带来的差异特别大?建议可以加点流速试试,看一下高流速下是不是还差这么多?

  • openFoam postProcess报错

    2 帖子
    127 浏览

    看起来好像你得OpenFOAM没编译完整

  • wsl2中安装ubuntu18.04后只有64个CPU分配

    1 帖子
    87 浏览

    wsl2中安装ubuntu18.04后只有64个CPU分配,按照教程新建并修改./wslconfigure文件后,关闭wsl并重启,仍然还是64个核,请问有人遇到到这种情况吗,谢谢!!!图片.png

  • 7 帖子
    292 浏览

    @coolhhh
    大佬您好,我可能忘了介绍我的研究时段了,我研究的是较长时间的模拟,理想模拟时间应该是2-3天;模拟区域也比较大,大概为长宽高为6km6km600m。
    1、因为研究的时间较长,所以每个边界的风速肯定是变化的(风速输入用的是每个边界所设置的非均匀分布风塔的观测数据,每个面都有数个观测点,观测数据风速利用timeVaryingMappedFixedValue每隔十分钟输入)因为之前学习了windaroundbuildings这个算例(也是建筑扰流),所以设置求解器为simpleFoam,这个算例风速确实是时不变的,这一点没有考虑到,等换成瞬时求解器试试;
    2、目前还没有试过五个入流面,因为我也认为可能不太符合通量守恒,所以先设置了3个速度入口和2个出口,我先做一下大佬在4中提供的建议;
    3、LES模拟也考虑过,正在学习,等做完2中试验后尝试一下大佬提供的条件,就不用计算k和epsilon文件的设置了。

  • 1 帖子
    77 浏览

    在ParaView 中用 Temporal Particles To Pathlines 追踪粒子轨迹,会出现多余的线条是怎么回事,求指导。

    abc25b9f-320d-4e32-8d77-4c45ac5ea775-image.png

  • 14 帖子
    559 浏览

    @coolhhh 我明白了,您不妨试试uniformFixedValue。另外dsmcFoam+的dsmcCloud.C中没有correctBoundaryConditions,我估计是包装在别的函数里了。

  • 1 帖子
    71 浏览

    在进行模拟时,运行136s出现上述报错:Primary job terminated normally, but 1 process returned
    a non-zero exit code. Per user-direction, the job has been aborted.
    2024-03-15_21-17.png

    附上我的算例边界条件,是一个基于真实地形,使用codefixedVal函数实现非均匀边界条件的模拟。1.zip
    请哪位大神给点建议,谢谢

  • 1 帖子
    81 浏览
  • chemkin中的therm.dat文件理解

    7 帖子
    291 浏览

    还有一个问题就是我在chemkin反应机理中添加了这个反应:
    1.png
    对于hv这个光子应该在哪里定义呀?我在chem.inp以及therm.date文件中没有定义光子,也不清楚如何定义,但是采用chemkinToFoam时显示了如下错误:
    屏幕截图 2024-03-15 190833.png
    hv改成大写以后就成了这样:
    1.png
    恳请各位老师,大佬解惑,感谢各位了!

  • divShemes离散格式的顺序问题

    6 帖子
    494 浏览

    :135: :135: :135:

  • 56 帖子
    10k 浏览

    @Prometheus10 对 如果这样理解的话 为什么symmetry计算结果和理论解更相近呢

  • interfaceHeight 运行报错

    8 帖子
    456 浏览

    @tidedrinker @李东岳 找到原因了,是因为我输入命令的问题,我之前输入的是postProcess,应该用 interFoam -postProcess , 。。。。新手不太会

  • MPPICFoam/denseParticleFoam 对于重力的考虑

    32 帖子
    2k 浏览

    @李东岳 李老师,最近在做实验,昨天又试验了一下,p场稳定了一些但是还在波动,但是U.air场变得混乱了,李老师有空可帮我看下吗,还有一个问题719e060c-5757-402d-b40a-0d978c2be0a7-tp_of10.zip
    b39e6d46-adc8-4f48-b766-281d3efcc0b7-image.png

    DPMFoam中的UcEqn.H:

    fvVectorMatrix UcEqn
    (
    fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc)

    fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc) continuousPhaseTurbulence->divDevRhoReff(Uc)
    ==
    (1.0/rhoc)*cloudSU
    );

    UcEqn.relax();

    volScalarField rAUc(1.0/UcEqn.A());
    surfaceScalarField rAUcf("Dp", fvc::interpolate(rAUc));

    surfaceScalarField phicForces
    (
    fvc::flux(rAUccloudVolSUSu/rhoc) + rAUcf(g & mesh.Sf())
    );

    if (pimple.momentumPredictor())
    {
    solve
    (
    UcEqn
    ==
    fvc::reconstruct
    (
    phicForces/rAUcf - fvc::snGrad(p)*mesh.magSf()
    )
    );
    }

    ///////还有就是这一项代表了什么意思呢////////
    surfaceScalarField phicForces
    (
    fvc::flux(rAUccloudVolSUSu/rhoc) + rAUcf(g & mesh.Sf())
    );

  • reactingFoam求解器修改化学反应机理

    3 帖子
    181 浏览

    @LiuWenchao 好嘞,我去看看,谢谢大佬!

  • 关于nut和mut的问题

    5 帖子
    282 浏览

    谢谢李老师,那我是不是直接用turbulence->nut()就可以获得混合湍流运动粘度了?我是看到createFields.H定义了这样的指针,对应找到他的 incompressible::turbulenceModel 类,就发现mut和nut是一样的,对于这种多相流的,它实际用的是哪个类,在哪里体现了密度的影响?对openFOAM不太了解,谢谢老师指点
    image.png

  • DPMFoam压力计算失稳

    8 帖子
    1k 浏览

    您好,我在修改了入射条件后依然没有解决这个问题,不知道您有什么解决办法吗

  • OpenFOAM外挂libtorch

    2 帖子
    260 浏览
  • LES定义入口速度的问题(DSRFG方法)

    82 帖子
    66k 浏览

    @xzdbjx 大佬,可以提供一下编译进Fluent里面的udf吗? CFD小白,自己写的总是报错:zoule:

  • 关于 rhoCentralFoam 中 fvSchemes 设置的疑问。

    3 帖子
    175 浏览

    @李东岳 多谢李老师,解析写的很清晰,基本弄明白了

  • 【分享+搬运】自定义非均匀inlet U

    11 帖子
    3k 浏览

    没事,找到错误就好。这个方法其实比较笨,因为每次你修改网格,就要重新处理一遍数据。
    有的时候会麻烦一些。

    祝收敛

  • 编译求解器的时候出现undedined reference错误

    4 帖子
    214 浏览

    倒数第二个图,这种提示都看不出来具体咋回事
    倒数第一个图,我挂libtorch的时候用过那个,你那个挂libtorch了么

  • 求教边界层生成问题

    17 帖子
    3k 浏览

    @上级 大佬,想问下cfMesh上手难吗?我在做相似的研究,想转成cfMesh试一下

  • openfoam如何中途调整更改网格?

    5 帖子
    287 浏览

    @学流体的小明 我也觉得难搞,我重组了某个时间步的网格文件,发现只有一个point文件,可能要一个点一个点的加😫,而且还有各种场文件貌似也要改😂

  • 鼓泡床参数设置的一个问题

    19 帖子
    1k 浏览

    好的,谢谢李老师!! 这个问题我纠结了好几天

  • 6 帖子
    1k 浏览

    @李东岳 明白了,谢谢!

  • 对时间和网格无关性的研究

    16 帖子
    794 浏览

    @crispy 那你看其他库朗数的结果怎么样,1以下的,是时间步长越大,cd越好么

  • 5 帖子
    291 浏览

    @李东岳 李老师,我发现应该是dev2(T(……))这一小段出的问题,删掉这个转置操作就能编译成功。因为我是在传热求解器的基础上进行修改的,T在createFields.H里声明了是温度场,所以是不是程序将T识别成了温度而进行不了转置操作了呢?

  • 关于openfoam边界条件的选择

    1 帖子
    124 浏览

    请教各位大佬一个边界条件的选择问题
    我的研究区域是一个三维区域,长宽高分别为(6000m*6000m*600m),如下图(俯视图)所示:
    08ae4b29-0f40-4ce5-b8b2-db687c3fc3fb-image.png

    图中的红色点是在研究区域内的观测站点,每一个观测点位可以得到从地面到600m高度处的风速、气压、温度等数据,
    请问应该选择什么边界条件来把这些观测点位的数据输入到openfoam中?

  • fvSolution中PIMPLE算法的参数设置问题

    5 帖子
    2k 浏览

    @WYingfvSolution中PIMPLE算法的参数设置问题 中说:

    @李东岳 谢谢李老师的回复!抱歉我的上个残差图绘制有问题,正确的应该是这样

    Residuals_Ur5.png

    从上到下依次是速度残差,压力残差和湍流统计量残差。

    上面那个code是我随意找的,我目前使用的code是:

    PIMPLE { correctPhi yes; nOuterCorrectors 2; nCorrectors 2; nNonOrthogonalCorrectors 1; turbOnFinalIterOnly false; } relaxationFactors { fields { p 0.3; } equations { "(U|k|omega)" 0.7; "(U|k|omega)Final" 1.0; } }

    我主要做了以下修改:
    (1)调整 nOuterCorrectors=2,结果几乎不变
    (2)turbOnFinalIterOnly=false,即每次PIMPLE loop都计算湍流,结果几乎不变
    (3)去掉全部的松弛因子,结果几乎不变。

    在模拟中,CFL保持在0.5-1.5之间,但是和实验结果对的不好。李老师觉得还有可能是哪里出了问题呢?

    请问同学,你的问题解决了吗

  • 关于InterFOAM曲率计算的问题

    5 帖子
    241 浏览

    @李东岳 嗯对,感谢李老师!

  • 服务器卡死关机重启后不能并行

    5 帖子
    265 浏览

    2、经验教训
    不知道说的对不对,前辈们可以建议一下,在新编译软件的时候需要注意,尤其是新加PATH的时候,编译后可能并不会直接出现问题,但是重启电脑将原有.bashrc重新source会引起错误

  • 6 帖子
    366 浏览

    !感谢大神的回答!

  • cfdem中cfd流场重组reconstructPar

    7 帖子
    392 浏览

    @Shihang-Chen 哇!感谢,确实解决了!

  • OpenFOAM python PINN tensorflow gpu (小白,莫笑)

    6 帖子
    1k 浏览

    有数据的PINN效果还行,但是完全的无数据的不太行

    我不太懂这个。确认一下:你说的这个意思,是说PINN需要CFD算出多套参数的结果,然后这些多套参数和结果就是数据。比如先用CFD算100套进口条件,得出100套结果,把这100套进口参数+结果给到PINN,这样的话效果还行?

  • 3阶以上矩阵定义与求逆

    9 帖子
    475 浏览

    绝对到位!嘎嘎一嘎子!

  • VOF添加斥力模型

    24 帖子
    2k 浏览

    另一种方法:仅添加到UEqn当中,压力方程不做处理,斥力的贡献依靠UEqn.H参与到后续的计算当中。结果就没有压力振荡了。

    :146: :146: :146: :146:

    Interesting!

  • 圆柱计算中无关性验证的问题

    34 帖子
    4k 浏览

    好的老师,我算完给您回复。

  • 3 帖子
    222 浏览

    谢谢岳哥,成了,牛牛牛

  • 关于cfdem计算发散

    2 帖子
    139 浏览

    只要有浮力就会发散,为什么呢

  • MPPICInterFOAM中粒子的interface力

    8 帖子
    456 浏览

    好的,谢谢李老师,我试一下

  • 14 帖子
    1k 浏览

    0时间步的U,是用湍流生成方法直接生成的,可以计算得到div(U)。然后就可以计算div(phi): U -> phi = fvc::flux(U) -> div(phi)。测试结果是,这种方式计算的div(U)和div(phi)结果基本一样的,区别是小数点后五六位。此时div(U)=div(phi)=0

    那这个算法很屌。不错。

  • CodedFixedValue边界设置问题

    4 帖子
    280 浏览

    @李东岳CodedFixedValue边界设置问题 中说:

    http://dyfluid.com/code.html 来这查一下

    好的,谢谢东岳哥

  • 用snappyHexMeshDict生成周期性边界问题

    14 帖子
    4k 浏览

    @Calf-Z-DNS用snappyHexMeshDict生成周期性边界问题 中说:

    @东岳 谢谢东岳老师
    我在cfd-online上面找到一种方法:
    下面是原文:
    “I encountered the same problem and I could solve it by defining rotational periodicity and then by declaring the corresponding vertices in ICEM as periodic. Then Right Click on Faces in the Blocking tree -> Periodic Faces -> ensure that really every periodic face is colored accordingly. createPatch went fine after this procedure.” ---by Fabian E.
    然后我勾选了ICEM里面blocking里的faces的Periodic,再导入openfoam,发现就不会出现上述错误。
    但我不太清楚为何这样会消除错误?

    请问同学可以告知是怎么设置的吗,icem里面有face pericdic吗?我这里只设置了点对称