Skip to content
  • triSurfaceMesh对象切割cell单元

    OpenFOAM
    2
    0 赞同
    2 帖子
    122 浏览
    A

    是想用基于cut CELL的浸没边界法? 可以参考FOAM Extend 4.1 之后的做法, 找ibPatch

  • 求助centos安装openfoam

    OpenFOAM
    2
    0 赞同
    2 帖子
    79 浏览
    A

    CentOS 7 上是能无root权限安装OpenFOAM的, 就是需要自己编译某些依赖, 比如说flex, bison之类的. 主要看想要安装OpenFOAM的哪些版本吧, 毕竟G++可能需要不同的版本

  • 关于wave2foam造波问题

    OpenFOAM
    8
    0 赞同
    8 帖子
    1k 浏览
    C

    @liujm 对的,通常我们计算整个风机都是采用komegaSST湍流模型,毕竟这个湍流模型基本上是做风机的都会用的,很少用其他湍流模型。

  • 0 赞同
    1 帖子
    31 浏览
    R

    大家晚上好
    1136009f-4795-4f3b-8f6f-9fdd3b5900fa-image.png
    图中红线对应的是thermicity peak,类似的我想实时获得每个y对应的热释率最大值所在的位置,想问一下有没有大佬有经验呀!

  • 离散伴随和连续伴随的优劣

    Algorithm
    8
    0 赞同
    8 帖子
    8k 浏览
    N

    @木木枭 我在做离散伴随,大佬有机会的话交流哇

  • 请教大佬们关于后处理中的阻力计算

    OpenFOAM
    9
    0 赞同
    9 帖子
    232 浏览
    N

    @xpqiu 感谢感谢

  • 网格y+求助

    OpenFOAM
    2
    0 赞同
    2 帖子
    86 浏览
    C

    网格只是全都正交的空计算域网格,还是有放几何的网格。记得之前用ICEM的 o-block 画的方柱网格,看起来质量还行,其实网格质量很差,用fluent的PISO求解器都发散

  • ubuntu 24 使用OF4

    OpenFOAM
    4
    0 赞同
    4 帖子
    366 浏览

    @Noob 这个要看你超算的singularity设置了。这个超算不是我们搞得,所以不太清楚,至少目前一个node是没问题的。
    不过我在运行的时候提示过inifiniteband 相关的library缺失。

    我不确定你的情况,具体你得问超算admin怎么设置singularity的。理论上应该是没问题的

  • 0 赞同
    1 帖子
    68 浏览

    近期尝试使用compressibleInterDyMFoam计算近壁面空化泡的溃灭,但在计算一定时间后出现了一些与实际不符的现象。
    算例设置如下:
    初始气泡半径R0=0.1mm,距离壁面距离λ=2,原计算域大小为5mm5mm1mm,实际计算域取其1/4。草图如下:
    1265c60d-3e45-457e-84df-081f40874abd-image.png
    初始网格密度为Δx=Δy=Δz=0.025mm。后续用自适应网格进行加密,dynamicMeshDict:

    FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object dynamicMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dynamicFvMesh dynamicRefineFvMesh; //dynamicFvMesh staticFvMesh; dynamicRefineFvMeshCoeffs { // How often to refine refineInterval 1; // Field to be refinement on field alpha.water; // Refine field inbetween lower..upper lowerRefineLevel 0; upperRefineLevel 0.999; // If value < unrefineLevel unrefine unrefineLevel 10; // Have slower than 2:1 refinement nBufferLayers 3; // Refine cells only up to maxRefinement levels maxRefinement 5; // Stop refinement if maxCells reached maxCells 6000000; // 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) (alphaPhiUn none) (alphaPhi10 none) (alphaPhi1Corr0 none) (ghf none) ); // Write the refinement level as a volScalarField dumpLevel true; }

    计算到10.2μs后,气泡突然从溃灭状重新开始增大,而并没有完全溃灭,如图所示。
    10.2μs:
    a833b3e2-ebc9-43b4-8363-d153ec8db0e0-10.2um.png
    10.4μs:
    e248f2ea-78d2-4100-a3a9-551f88c499e3-10.4um.png

    想问问大家这是什么原因所导致的,可能会是初始场设置不合理或是边界条件设置不合理?
    初始场由于初始网格较疏以及后续加密的缘故,几个时间步后会变成这样:
    170ce1d8-861d-4cfb-857b-2e2da4ed06eb-初始场.png

    边界条件如下所示:
    alpha.water:

    boundaryField { wall { type zeroGradient; } outlet { type zeroGradient; } symmetry { type symmetry; } }

    p:

    boundaryField { wall { type calculated; value uniform 100000; } outlet { type calculated; value uniform 100000; } symmetry { type symmetry; } }

    p_rgh:

    boundaryField { wall { type fixedFluxPressure; gradient uniform 0; value uniform 100000; } outlet { type fixedFluxPressure; gradient uniform 0; value uniform 100000; } symmetry { type symmetry; } }

    T:

    boundaryField { wall { type zeroGradient; } outlet { type zeroGradient; } symmetry { type symmetry; } }

    T.air:

    boundaryField { wall { type zeroGradient; } outlet { type zeroGradient; } symmetry { type symmetry; } }

    T.water:

    boundaryField { wall { type zeroGradient; } outlet { type zeroGradient; } symmetry { type symmetry; } }

    U:

    boundaryField { wall { type noSlip; } outlet { type zeroGradient; value uniform (0 0 0); } symmetry { type symmetry; } }
  • 0 赞同
    2 帖子
    121 浏览

    我现在想明白了,我要自定义的热流条件是基于流固界面两侧物理量的,也就是说在不coupled的条件下,也体现了高温风流对低温相变板的加热作用,此时流体一侧的shadow面绝热是没有问题的,湿空气与相变板之间传热的相互作用体现在两者温差上。

  • 0 赞同
    4 帖子
    3k 浏览
    WeakFormW

    我参考HZDR的代码弄了个setbubble小工具,让气泡初始化就是圆的:5be8a963-f837-4652-9c92-a44ffd016afd-image.png
    但是我的问题是进行计算之后貌似界面的细化很差,很多界面细化网格级别没达到最低,不知道您怎么解决的这个问题:
    2964e223-615d-4d60-81f0-7a70f7db825d-image.png

  • 0 赞同
    142 帖子
    581k 浏览
    C

    @random_ran 老师您好 请问您论文题目叫啥呢 我打开链接找不到您的论文

  • 0 赞同
    1 帖子
    86 浏览

    非均匀温度边界条件可以通过codedFixedValue轻松实现,
    带壁厚的温度边界条件可以通过externalWallHeatFluxTemperture实现。
    如代码所示:

    \\非均匀温度边界条件 wall { type codedFixedValue; name dummy; code #{ const vectorField& Cf = patch().Cf(); const scalar a0 = 0.1; forAll(Cf, faceI) { const scalar y = Cf[faceI].y(); (*this)[faceI] = a0*y; } #}; } \\薄壁温度条件 wall { type externalWallHeatFluxTemperature; mode coefficient; h uniform 1e10; Ta 1620; thicknessLayers (0.015); kappaLayers (4); kappaMethod fluidThermo; value $internalField; }

    但是,当我想像下面这样编程使Ta为非均匀值时,会报错,即使我使用了codeInclude codeOptions codeLibs之后,仍会报错 ‘patch’ was not declared in this scope,似乎在externalWallHeatFluxTemperature中,无法访问到网格信息。请问各位大佬有解决办法吗?或者有没有别的办法实现带薄壁的非均匀温度边界呢?

    ‘vectorField’ does not name a type ‘Cf’ was not declared in this scope Ta coded; code #{ const vectorField& Cf = patch().Cf(); const scalar a0 = 0.1; forAll(Cf, faceI) { const scalar y = Cf[faceI].y(); (*this)[faceI] = a0*y; } #};
  • 0 赞同
    2 帖子
    2k 浏览
    余正东

    能给出这篇文章的名字吗

  • fluidsolver和isothermalfluid的继承关系问题

    OpenFOAM
    1
    0 赞同
    1 帖子
    82 浏览
    M

    各位老师好,之前我都在of10中修改求解器,但最近换到了of13,我希望去更改buossinesq假设时的求解器代码,但是对于fluidsolver和isothermalfluid的继承感到困惑,在官网看到fluidsolver继承isothermalfluid,但在fluidsolver中我并未找到相关代码,因此想请教各位在模块化代码下要怎么去修改求解器,谢谢大家!

  • 两相流自由液面处速度场异常

    OpenFOAM
    1
    0 赞同
    1 帖子
    140 浏览
    liujmL

    各位老师好,

    最近在使用K-omega SST湍流模型处理两相流的时候发现自由液面处的速度异常,但是不影响波面,如图所示。请问,这种情况是什么引起的,会影响计算结果吗?要如何解决?

    万分感谢!

    662a54c8-2f33-4bff-ac92-c8c5059e396d-image.png

  • 0 赞同
    13 帖子
    6k 浏览
    L

    @李东岳 很奇怪,又画了个很稀的6万网格量的网格,先前会发散的极低雷诺数,现在却不发散了。看来雷诺数极低,网格就不能太密?

    通过网盘分享的文件:2D_21cylinders_Re0.0001.zip
    链接: https://pan.baidu.com/s/1NdfFZrV2k_T8mAY7fA7ktg?pwd=8888 提取码: 8888

  • 使用cavitatingfoam模拟气泡溃灭

    OpenFOAM
    7
    0 赞同
    7 帖子
    3k 浏览

    @李东岳 老师,请问我想算空化泡溃灭,但似乎compressibleInterFoam并不包含传质项?是要自己改成compressibleInterPhaseChangeFoam嘛?

  • 如何使用moveDynamicMesh初始化变形网格

    OpenFOAM
    4
    0 赞同
    4 帖子
    237 浏览
    C

    @liujm 。。。。。。。。。,除了你画图的时候浮体的重心位置变化了,其他设置都没啥变化的。

  • 如何验证网格的二阶收敛精度

    OpenFOAM
    1
    0 赞同
    1 帖子
    131 浏览
    L

    simpleFoam求解的定常问题,用的空间离散格式都是二阶的,审稿人让验证网格的二阶收敛精度。不知有没有人用过paraview将稀网格数据插值到密网格上并算过残差E和收敛阶次p?另外对于很复杂的三维网格,网格的非均匀性对收敛阶次公式的计算又该如何考虑呢?
    35fb0dd0-05d4-4fcf-b850-3b3e8d317042-image.png

    be5eddf7-3dd1-41c7-b320-e9fa9531fffc-image.png

京ICP备15017992号-2