OpenFOAM

4.5k Topics 26.2k Posts
  • 3 Posts
    279 Views

    @李东岳MULES::correct() 和 MULES::explicitSolve() 有什么区别? 中说:

    MULES要看基金会的版本,这是基金会那面提出的。我记得CMULS已经被删了。我在无痛苦ns方程笔记里面说的很齐全。你看看参考一下

    好的,谢谢李老师

  • case外包,求助大佬们

    3 Posts
    333 Views

    导师经费很多,可详谈

  • 如何增加0时刻的采样参数输出

    8 Posts
    730 Views

    按照V8的release说明,除了某些特别的需要依赖时间进行平均的函数,其余的比如sample采样等函数会在0时刻(如果startTime是0)先执行得到初始时刻的值。这样在最后的dat中就会有0时刻的数据了。
    commit原文在这。

  • 求助各位老师,谢谢

    1 Posts
    142 Views

    我想在openfoam里面引入液体tait状态方程,该怎样实现啊?make文件夹怎么编写?

  • 《OpenFOAM用户指南》勘误

    Moved
    31 Posts
    14k Views

    @kuangyonglin 谢谢!下一版本更新

  • 1 Posts
    111 Views

    ee6ad707-d3e6-48f0-b4cf-ab27fb7423f0-image.png
    计算虚拟机和本地计算机共享的文件夹里的算例时候卡顿,尤其是并行计算后进行文件合并,打开终端需要等半分钟,有谁有解决办法吗?

  • OvermultphaseInterDyMFoam

    7 Posts
    3k Views

    @wzlr2023 我看你这个算例里只有两相流,正常用 interFoam足够了。这个相关资料还挺多的。这个多相流的求解器正如我前面所说我没有去验证其是否可以正常求解。就你这个算例来讲,里面涉及到了六自由度运动求解。造成两相交界面求解发散的问题也可能是六自由度求解引起的流场发散。建议可以先设置一个简单的动网格算例,如动边界造波来验证多相流求解是否稳定。如发现是多相流求解出问题了,欢迎进行补充修正。

  • 水平管中大雷诺数颗粒沉降问题

    1 Posts
    151 Views

    我使用CFDEM中的解析CFD-DEM模型+RANS+k epsilon计算水平圆管中颗粒沉降的时候(管径115mm,入口流速3m/s,粘度0.01Pa·s,密度1000kg/m3;颗粒密度1100kg/m3,直径15mm;雷诺数几万),发现颗粒越靠近圆管底部,下沉速度越慢,甚至在距离圆管中心0.4mm的时候,颗粒没有下沉,反而向上移动了。想请教一下各位老师,这种情况是不是不对呢?
    印象中层流里面,管壁附近重颗粒才会往管道中间移动,不知道我这个认识对不对。

  • 5 Posts
    760 Views

    @luca请教在自开发AMR求解器中应用cyclic边界条件报错 中说:

    单核是可以运行的,但是一旦并行运算就会报错。

    我在使用并行计算时也遇见了相同的问题,题主解决了吗?

  • 20 Posts
    2k Views

    学流体的小明 在 interFoam计算气泡槽道流时的压力问题 中说:

    应该是 并行 + 自适应网格 的问题?

    更正一下,是 并行+自适应网格+压力参考点 三个要素的共同作用,去掉任何一个都可以算出来好的结果。

    一个新的发现是壁面上会有奇异点,如下图,着色是p_rgh,这些奇异点都处于并行分区的各个界面上。更奇怪的是,上壁面并没有这样的问题。Note:参考点放置在流场中间。
    01292e97-19c8-45be-889b-189c0ff6fe54-image.png

    现在暂时放弃自适应加密了,直接全场加密算,压力就是正确的。
    😂

  • 自定义梯度边界条件

    10 Posts
    904 Views

    楼上说的对,嗯,这是面法向梯度

  • 2 Posts
    244 Views

    可不可以结合codestream实现旋转呢?

  • outfile精度好像不够

    4 Posts
    954 Views

    @Samuel-Tu 请问你解决了这个问题了吗?我是在controlDict中添加code进行计算得到的一些自定义的变量,但是结果在写入到文件中时,除了时间步数之外,其他的参数就一直是6位数或更少,小数点后边的就更少了。也尝试改了李老师说的这个,但是好像没用。

  • Coupled level set-VOF方法

    4 Posts
    2k Views

    @Zhujh 您好,也采用了与您基本相同的CLSVOF方法进行计算,然而对于毛细张力主导的问题,这个方法表现出了更强的虚假流动(寄生流动)的问题,如无重力状态下的水滴。请问您遇到过类似的问题么?同时我也参考fluent里面的两种抑制虚假流动方法(密度和H函数),但是收效甚微,想问下您有什么建议吗

  • 6 Posts
    597 Views

    @zhoushan 简单的几何可以用blockMesh直接生成网格,复杂的就要绘制stl文件,用blockMesh搭配snappyHexMesh生成

  • 求助求助 foam extend 4.0安装fsi失败

    1 Posts
    98 Views

    fluid@zhoushan:~/$WM_PROJECT_USER_DIR/FluidSolidInteraction/src$ ./Allwmake

    wmake libso fluidSolidInteraction
    ./Allwmake: 4: wmake: not found wmake solvers/fsiFoam
    ./Allwmake: 6: wmake: not found wmake solvers/ampFsiFoam
    ./Allwmake: 7: wmake: not found wmake solvers/weakFsiFoam
    ./Allwmake: 8: wmake: not found wmake solvers/fluidFoam
    ./Allwmake: 9: wmake: not found wmake solvers/solidFoam
    ./Allwmake: 10: wmake: not found wmake solvers/thermalSolidFoam
    ./Allwmake: 11: wmake: not found wmake libso utilities/functionObjects/pointHistory
    ./Allwmake: 13: wmake: not found wmake libso utilities/functionObjects/patchAvgTractionHistory
    ./Allwmake: 14: wmake: not found wmake libso utilities/functionObjects/centrifugalBodyForce
    ./Allwmake: 15: wmake: not found
    请问上述的not found怎么解决呀 在网上参考了各种方法 都无法成功
  • 1 Posts
    132 Views

    最近在做多相流的仿真验证,但是一直找不到LS方法的多相流案例,想请教以下各位大佬openfoam中有没有关于LS方法的多相流案例

  • 13 Posts
    3k Views

    Stop refinement if maxCells reached maxCells 200000 这个参数可以调大吗? 单核计算到一定时间后也报错了,这个报错是什么意思@李东岳

    image.png

  • 1 Posts
    124 Views

    如上所述

  • 5 Posts
    438 Views

    @TM 好的,感谢感谢👍👍

  • openfoam添加湍流入口方法?

    45 Posts
    14k Views

    @dyb 谢谢你的回答,我还有个问题你这边L一般就是取特征长度的10%来做是么,这个是按边界网格为基准计算么

  • open foam空化流出现负压

    2 Posts
    245 Views

    这个问题我也遇到过,比如对于同一个流场,同样的空化数,来流速度10m/s就不会产生负压,来流速度20m/s就有负压产生。增加SS空化模型的系数n的值或许能解决问题。

  • 1 Posts
    156 Views

    大家好,请问大家在使用sixDoFRigidBodyMotion进行流固耦合的时候,在仿真完毕之后如果想要续算是怎么设置的?

    我使用的是OpenFOAM V2306,使用的求解以下是我的dynamicMeshDict文件:

    /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v1912 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object dynamicMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dynamicFvMesh dynamicOversetFvMesh; motionSolverLibs (sixDoFRigidBodyMotion); motionSolver sixDoFRigidBodyMotion; dynamicFvMesh dynamicOversetFvMesh; sixDoFRigidBodyMotionCoeffs { patches (cylinder); innerDistance 100; outerDistance 101; mass 10; centreOfMass (0.0 0.0 0.0); momentOfInertia (1 1 1.5625); g (0 0 0); rho rhoInf; rhoInf 1; report on; solver { type Newmark; } constraints { yLine { sixDoFRigidBodyMotionConstraint line; direction ( 0 1 0); } rotation { sixDoFRigidBodyMotionConstraint axis; axis (0 0 1); } } restraints { spring1 { sixDoFRigidBodyMotionRestraint linearSpring; anchor (0.0 0.0 0.0); refAttachmentPt (0.0 0.0 0.0); stiffness 10.96622711232151; damping 0.6283185307179585; restLength 0; } spring2 { sixDoFRigidBodyMotionRestraint linearAxialAngularSpring; axis (0 0 1); stiffness 1.7134729863002358; damping 0.09817477042468103; } } } // ************************************************************************* //

    以下是正常仿真300s之后的log.overPimpleDyMFoam中最后一步和sixDoFRigidMotion和overset有关的输出:

    Restraint spring1: attachmentPt - anchor (0 -0.367945 0) spring length 0.367945 force (-0 4.0584 -0) Restraint spring2: angle 0.516284 moment (0 0 -0.883591) 6-DoF rigid body motion Centre of rotation: (0 -0.368716 0) Centre of mass: (0 -0.368716 0) Orientation: (0.869769 -0.493459 0 0.493459 0.869769 0 0 0 1) Linear velocity: (0 -0.372264 0) Angular velocity: (0 0 -0.107739) inverseDistance : detected 2 mesh regions zone:0 nCells:50721 voxels:(40 40 1) bb:(16.9781 3.77358 -0.500006) (20.9906 8.00001 0.500006) zone:1 nCells:2280 voxels:(40 40 1) bb:(-72 -40 -2.50004) (-40 -24 -1.49996) Overset analysis : nCells : 53001 calculated : 52701 interpolated : 189 (from local:0 mixed local/remote:0 remote:189) hole : 111

    如果我想续算,直接修改controlDict中的endTime显然是会出现问题的,它会使用原来的dynamicMeshDict中的数据作为300s的初值,导致计算无法正常进行。

    而如果我参照最后一步的输出修改dynamicMeshDict文件:

    /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v1912 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object dynamicMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dynamicFvMesh dynamicOversetFvMesh; motionSolverLibs (sixDoFRigidBodyMotion); motionSolver sixDoFRigidBodyMotion; sixDoFRigidBodyMotionCoeffs { patches (cylinder); innerDistance 100; outerDistance 101; mass 10; // centreOfMass (0.0 0.0 0.0); centreOfMass (0 -0.368716 0); momentOfInertia (1 1 1.5625); orientation (0.869769 -0.493459 0 0.493459 0.869769 0 0 0 1); velocity (0 -0.372264 0); angularMomentum (0 0 -0.1683421875); // angularVelocity (0 0 -0.107739); g (0 0 0); rho rhoInf; rhoInf 1; report on; solver { type Newmark; } constraints { yLine { sixDoFRigidBodyMotionConstraint line; direction ( 0 1 0); } rotation { sixDoFRigidBodyMotionConstraint axis; axis (0 0 1); } } restraints { spring1 { sixDoFRigidBodyMotionRestraint linearSpring; anchor (0.0 0.0 0.0); // anchor (0 -0.368716 0); // refAttachmentPt (0.0 0.0 0.0); refAttachmentPt (0 -0.368716 0); stiffness 10.96622711232151; damping 0.06283185307179585; restLength 0; } spring2 { sixDoFRigidBodyMotionRestraint linearAxialAngularSpring; axis (0 0 1); stiffness 1.7134729863002358; damping 0.009817477042468103; // referenceOrientation 0; } } } // ************************************************************************* //

    首先根据最后一步的数据修改dynamicMeshDict文件我不知道会不会带来新的误差,然后这种方法也没有明显效果,以下是修改后续算的第一步输出:

    Restraint spring1: attachmentPt - anchor (0 -0.368716 0) spring length 0.368716 force (-0 4.06681 -0) Restraint spring2: angle 0.516062 moment (0 0 -0.883201) 6-DoF rigid body motion Centre of rotation: (0 -0.369639 0) Centre of mass: (0 -0.369639 0) Orientation: (0.869901 -0.493226 0 0.493226 0.869901 0 0 0 1) Linear velocity: (0 -0.371839 0) Angular velocity: (0 0 -0.10837) inverseDistance : detected 2 mesh regions zone:0 nCells:50721 voxels:(230 230 1) bb:(-8.00004 -8.00004 -0.500036) (24 8.00004 0.500036) zone:1 nCells:2280 voxels:(230 230 1) bb:(-0.999904 -1.00093 -0.500003) (1.0001 0.99908 0.500003) Overset analysis : nCells : 53001 calculated : 52683 interpolated : 202 (from local:202 mixed local/remote:0 remote:0) hole : 109

    6-DoF相关的输出,通过修改dynamicMeshDict已经和前面比较接近了,但是zone 0, zone 1相关的输出有些差异,同时力系数突然变得很大:

    forceCoeffs forceCoeffs write: Coefficient Total Pressure Viscous Internal Cd: -734104 -733933 -171.164 0 Cd(f): -367052 -366966 -85.5821 0 Cd(r): -367052 -366966 -85.5821 0 Cl: -458471 -457934 -537.556 0 Cl(f): -295374 -294999 -375.408 0 Cl(r): -163097 -162935 -162.148 0 CmPitch: -66138.8 -66032.2 -106.63 0 CmRoll: -7.89329e-15 -7.89337e-15 7.8177e-20 0 CmYaw: 8.40485e-15 8.40299e-15 1.86337e-18 0 Cs: 2.00416e-14 2.00378e-14 3.83127e-18 0 Cs(f): 1.84257e-14 1.84219e-14 3.779e-18 0 Cs(r): 1.61596e-15 1.61591e-15 5.22695e-20 0 writing force and moment coefficient files.

    导致后续计算无法正常进行(时间步长变得很短)。

    所以我想请教一下大家在使用sixDoFRigidMotion进行流固耦合仿真的时候,如果想要续算有什么好方法吗?因为似乎openfoam没有保存sixDoFRigidMotion的相关参数(如质心位置、速度、加速度等)。

  • 1 Posts
    265 Views

    最近在植入一种DES混合模型。期望我植入的DES模型的时间平均速度,与RANS一致。对比后发现不一样。debug几个小时没发现问题。索性采用OpenFOAM原生的Smagorinsky算了一下。发现原生的Smagorinsky模型计算的时间平均的速度跟RANS也不一样。

    简述:OpenFOAM的pitzDaily算例,用2DLES计算与2DRANS计算时间结果不一致,我怀疑是不是2D LES的问题?

    因为我这跑不动3D LES,所以没法做3D LES与2D RANS的结果对比,大半夜的我感觉啪的一下就起来了,很快啊!问问各位大佬有何高见。

    屏幕截图 2023-12-15 080858.png

    下图是2D RANS的结果:图片出处
    替代文字

    下面是我用原生的2D Smagorinsky算的时间平均速度:

    屏幕截图 2023-12-14 171725.png

    下图是我自己植入的2D DES算的时间平均速度:

    屏幕截图 2023-12-14 175456.png

    目前没有3D LES的结果。纵然我有万核CFD服务器在闲着,但debug起来太慢了主要。。可以理解为跑(debug)不动,没时间持续研究这个事情。

    :136:

  • Optuna + OpenFOAM 寻优

    1 Posts
    180 Views

    今天偶尔看到一个文章。好像意思是通过Optuna + OpenFOAM,可以实现最小的计算时间。感觉类似通过Optuna跑好几次,然后寻找出来最合适的松弛因子。不知道有没有哪个大佬感兴趣玩一玩。

    https://ss1.xrea.com/penguinitis.g1.xrea.com/study/OpenFOAM/relax_factor_optimization.html

  • paraview 不显示坐标轴数值

    5 Posts
    1k Views

    @wangfei9088paraview 不显示坐标轴数值 中说:

    @liujm 假设你用的版本是OpenFOAM-10,其他版本也一样。
    路径:OpenFOAM-10/etc/config.sh/
    打开paraview文件,将69行取消注释,注释掉72行。
    ./makeParaView重新编译就是5.0.1版本了。

    #export ParaView_VERSION=5.0.1 //69行 #export ParaView_VERSION=5.4.0 #export ParaView_VERSION=5.5.0 export ParaView_VERSION=5.6.3 //72行

    感谢!

  • 19 Posts
    3k Views

    @LGY 请问一下你这边能够在wave2Foam中生成非定常风了吗

  • 8 Posts
    729 Views

    @gooseEastrhoCentralFoam发散,出现Maximum number of iterations exceeded 中说:

    在thermoI.H文件第46行标量f是什么

    OpenFOAM代码中,这个f是Cp,定压比热容,单位是[J/(kg K)]。这个T函数是通过定压比热容的值迭代求出温度。limit,F,dFdt分别对应的是hConstThermoI.H文件里limit,Cp,dCpdt三个函数。

  • 流体热物理属性是温度的函数

    1 Posts
    139 Views

    我在使用buoyantSimpleFoam时需要将流体密度、黏度、热导率、定压比热容等设置成温度的插值函数,通过如下的thermophysicalProperties字典文件能够实现依照温度插值的效果吗?

    /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2306 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object thermophysicalProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // thermoType { type heRhoThermo; mixture pureMixture; transport tabulated; thermo hTabulated; equationOfState icoTabulated; specie specie; energy sensibleEnthalpy; } mixture { specie { molWeight 1; } equationOfState { rho ( (300 1) (400 1) ... (1000 1) ); } thermodynamics { Hf 0; Sf 0; Cp ( (300 1) (400 1) ... (1000 1) ); } transport { mu ( (300 1) (400 1) ... (1000 1) ); kappa ( (300 1) (400 1) ... (1000 1) ); } } // ************************************************************************* //
  • 3 Posts
    240 Views

    @李东岳 老师,如果不使用refineMesh的话,结果会好很多,而且和论文中的Cd, Cl指数可以基本对上,但是在overset的边缘处还是会有一点问题。
    648bb924-7f1b-4e8c-a39d-48cc7c40906e-image.png

    我的oversetInterpolation的method使用的是inverseDistance.

  • 关于waves2foam造波衰减问题

    6 Posts
    2k Views

    首选请明确是否调用了湍流模型!

    Larsen证明了两方程模型会导致准势流区域的湍动能会指数增长,波高沿程衰减。

    建议调用Larsen改进的stablized版本的额湍流模型。

  • paraview中矢量的表示出现了问题

    2 Posts
    183 Views

    求大佬帮我看一下,跪求

  • 6 Posts
    429 Views

    @李东岳
    当边界不做判断直接叠加之后是这样的
    fc3c9e32-1d39-4a41-b44a-85c2594da99b-image.png

    bottom
    {
    type codedFixedValue;
    value uniform 101325;
    name face1P;
    code #{
    scalarField& field = *this;
    const label id = this->patch().boundaryMesh().findPatchID("bottom");//在网格中找到inlet
    const volScalarField& pp = this->db().objectRegistry::lookupObject<volScalarField>("p");
    const scalarField& ppp=pp.boundaryField()[id].patchInternalField();
    forAll(field, faceI)
    {

    field[faceI]=ppp[faceI]+100; } #}; }
  • HPC 上编译 waves2foam 的问题

    5 Posts
    632 Views

    @peng23 https://cfd-china.com/topic/4499/openfoam第三方造波包waves2foam安装 你看看这个。waves2foam要注意版本一致性的问题。

    另外,GLIBCXX_3.4.29' not found应该是你的ubuntu系统太高。

  • 关于晃荡惯性力的植入

    11 Posts
    824 Views

    怎么样 后来有进展么

  • 3 Posts
    294 Views

    @李东岳 多谢李老师指导,我去试一下

  • 1 Posts
    143 Views

    我想把图片中的红色框的一部分输出,然后在paraview中做后处理,为了方便,我建立了一个同样的参数(黄色框),想用类似湍流模型里边的方法试一下,但是写不出来b910a92ec6b2db806f3cf7ed6ac073a.png
    我还尝试了另一种方法:在createfield.h里面定义了一个参数,这样会写出一个文件,但是写出的文件不会更新。55cdd8962f572210b4507aa805be817.png
    其中m是自定义的一个参数453c8ef0-3b5a-481d-90df-3fb52ad4320e-image.png

  • 3 Posts
    936 Views

    感觉跟setExprFieldsDict的语法一样

  • 边界积分

  • 推板造波-孤立波

    2 Posts
    161 Views

    请查看上海交通大学,宣瑞韬、张金牛、荣一毅等人的硕士论文。

  • 1 Posts
    117 Views

    在solids4foam中运行验证算例(孤立波下斜坡上的直墙变形和受力),但是当波浪到固体交界面时,就会出现浮点报错问题,修改了网格和固体交界面的问题,还是没有改善这个问题。请问有大佬们知道怎么解决这个问题吗?
    报错时间步alpha场:(每次到这个地方就出现报错)
    2a67acf0-d296-4f56-82eb-f8a159853cca-image.png

    报错部分log文件:

    Time = 8.933, iteration: 9 Current fsi under-relaxation factor (Aitken) of interface: 3.24861e-06 Maximal accumulated displacement of interface 0: 2.02709e-06 GAMG: Solving for cellMotionUx, Initial residual = 0.27053, Final residual = 4.88471e-05, No Iterations 2 GAMG: Solving for cellMotionUz, Initial residual = 0.276424, Final residual = 7.60589e-05, No Iterations 2 GAMG: Solving for cellMotionUx, Initial residual = 5.19396e-05, Final residual = 4.62925e-07, No Iterations 4 GAMG: Solving for cellMotionUz, Initial residual = 8.68504e-05, Final residual = 8.18649e-07, No Iterations 3 Evolving fluid model: wave2solids PIMPLE: iteration 1 MULES: Solving for alpha1 alpha1 BC on patch inlet alpha1 BC on patch inlet Liquid phase volume fraction = 0.259489 Min(alpha1) = -2191.92 Max(alpha1) = 2360.17 MULES: Solving for alpha1 alpha1 BC on patch inlet alpha1 BC on patch inlet Liquid phase volume fraction = 0.259489 Min(alpha1) = -158062 Max(alpha1) = 163488 Velocity BC on patch inlet 3D_2D Absorption BC on patch outlet "Correction Levels" 1( 0 ) DICPCG: Solving for pd, Initial residual = 0.10778, Final residual = 0.00369688, No Iterations 8 DICPCG: Solving for pd, Initial residual = 6.32049e-06, Final residual = 2.94799e-07, No Iterations 644 DICPCG: Solving for pd, Initial residual = 4.65247e-05, Final residual = 2.1618e-06, No Iterations 637 DICPCG: Solving for pd, Initial residual = 2.0431e-07, Final residual = 9.74494e-08, No Iterations 110 Velocity BC on patch inlet 3D_2D Absorption BC on patch outlet "Correction Levels" 1( 0 ) DICPCG: Solving for pd, Initial residual = 4.64318e-05, Final residual = 2.21746e-06, No Iterations 638 DICPCG: Solving for pd, Initial residual = 2.16117e-07, Final residual = 8.27667e-08, No Iterations 132 Velocity BC on patch inlet 3D_2D Absorption BC on patch outlet "Correction Levels" 1( 0 ) Setting traction on solid interfaces Mapping face values using directMap Total force on fluid interface 0: (4.3321e+10 -2.82466e-19 -0.673933) Total force on solid interface 0: (-4.3321e+10 2.82466e-19 0.673933) Evolving solid solver Solving the updated Lagrangian form of the momentum equation for DD -------------------------------------------------------------------------- mpirun noticed that process rank 6 with PID 168600 on node node03 exited on signal 8 (Floating point exception).
  • piso和pimple

    3 Posts
    3k Views

    @李东岳piso和pimple 中说:

    数值稳定性

    李老师您好,我在用buoyantPimple求解器的PSIO模式计算流场压力、密度和速度出现明显的数值振荡,就是随时间步长周期性振荡发展,但是振幅确实是一直在减小的,可否理解为这就是数值振荡问题,通过调整PIMPLE求解器外部动量预测修正次数转化为PIMPLE增强稳定性?
    计算时序图如下:
    8055ef71-35b7-40e6-b692-ac6aa27a4f2c-image.png
    58fdb9e2-331b-4f6b-9a4a-39cb7b94f1bf-image.png
    064a77b7-1f57-4077-9480-e6395c761979-image.png

  • 1 Posts
    138 Views

    paraview以列表形式导出流场所有数据时导出太慢(每个时间步的数据比较小时无影响)。如果用自动写每个时间步的数据时,则每个时间步的数据都不全。请问各位有遇到过吗?有无解决办法。

  • 9 Posts
    2k Views

    学流体的小明 在 压力基求解器在OpenFOAM中的植入问题 中说:

    CaseB中,同样的方法,压力p基本都是10000,或者有是9999.99,这和我限制了writePrecision 6有关系。

    所以输出流场的时候,最好以binary二进制的格式输出。
    后处理的文件,由于要使用其它软件读取,就用ASCII。

  • 怎么用ICEM画楔形体网格

    21 Posts
    3k Views

    最近也在折腾做旋转对称圆柱型计算域,找到个这:https://www.fangzhenxiu.com/post/1826059?code=031TTpFa1IW6xG0skbHa1OsJy23TTpFt&state=afterWeixin
    根据这个操作画出来的网格虽然是楔形块,但是楔面并不参照哪一坐标平面对称,依然不能用在OpenFOAM上。
    又找到一个这个:https://www.fangzhenxiu.com/post/1829359
    作者好像通过两次旋转实现了对称的楔面,但我没有复现出来。
    最后还是打算用blockMesh来画了。

  • 13 Posts
    1k Views

    @ShaneHEEE 是的,和出口设置一样,就可以了

  • DSMC时间平均

    1 Posts
    121 Views

    dsmc求解中,由于是使用统计学相关的算法进行求解,计算结果具有较大的统计噪音,使用时间平均法是降噪的一种有效方法.

    作者在学习相关软件时,发现tecplot具有pytceplot拓展模块,能够实现此功能.

    首先需要在系统中装python环境,并且安装pytecplot模块

    直接使用pip install pytecplot命令即可安装

    在官方github仓库中找到tputils.py, tpmath.py, TimeAverage.py3个文件,将其放到OpenFOAM的结果文件夹下

    使用tecplot打开contralDict文件,并且打开软件的PyTecplot Connections功能

    在终端使用python环境运行TimeAverage.py脚本,待程序运行完成,在zone style中将会出现时间平均后的计算结果

    以某微纳通道为例,进行时间平均,并做对比:

    时间平均前的结果:
    c03304cf-34df-40a0-b708-a63631857ecf-image.png

    时间平均后的结果:
    3247da89-b756-4a9c-a8ed-8bde44fd4f0c-image.png

    显然,时间平均起到了降噪的作用

  • 2 Posts
    227 Views

    各位老师,我主要是想消除这种异常情况,然后刚刚测试了一下发现对于压力设置totalPressure等压力控制入口,速度与zeroGradient相关设置的边界都可能出现这个问题。想请问一下,这样的情况可能与什么有关呢?

  • 10 Posts
    7k Views
  • 18 Posts
    5k Views

    @HITSC30 你好,请问你这个问题解决了嘛?我也遇到了和你一样的报错。