OpenFOAM

4.4k 主题 25.9k 帖子
  • 13 帖子
    2k 浏览

    时至今日,我是刚刚接触这方面小白,问下风力机致动线模型是否有相关算例可以学习,如有打扰,非常抱歉

  • 3 帖子
    164 浏览

    @李东岳非结构化网格如何加显式滤波 中说:

    OpenFOAM有一些LES模型需要调用显示滤波。比如动态k模型。OpenFOAM已经处理好了。你可以看看源代码。在我那个笔记里面也写的非常详细,不管是理论,还是方程方面。

    多谢老师

  • 4 帖子
    265 浏览

    @李东岳 谢谢李老师,也谢谢楼上几位。目前找到一些问题所在了。在fluent中进出口的温度不是周期性的而流场是,每个周期都会重置入口温度(之前用fluent的人不清楚),所以这一点是我们自己的问题所在。另外还有一个问题就是在计算之前fluent进行初始化为300K,但是计算一段时间后会发现有小于300K的,这个问题就不知道是什么原因了。目前计划是对温度场T+、y+了。和实验的验证结果如下图所示,个人感觉能凑合用吧。
    9a7bebfa-8bf1-450e-8768-5e386dd453a2-微信截图_20240115155859.png

  • OpenFOAM循环壁面生成

    2 帖子
    130 浏览

    感谢分享!!:146: :146: :146:

  • parafoam打开加密网格

    3 帖子
    166 浏览

    @cresendo 但是我每次运行前都会allclean。是在什么时候删除时间步呢,toposet之前还是refinemesh之后呢?

  • 不可压缩多相流计算出现负压

    5 帖子
    282 浏览

    谢谢大佬,我的Courant数设置的是0.5,时间步长使用的自适应时间步长。我的想法是可不可能是因为alpha.water场某个单元计算出来比较小,导致得到该单元的U特别大,然后求解压力方程后在该单元附近出现负压?

  • 1 帖子
    65 浏览

    使用OpenFOAM-v2006版本的setExprBoundaryFieldsDict在设置某个坐标方向上边界值符合分段函数分布规律的初值时,比如z方向上再0<z<11000m 时符合T = 288.15 - 0.0065*z,在z>11000m时T = 216.65,分段函数如下:
    7332c5a5-eb8a-48b7-b4d0-b2effae3bbdc-image.png
    这样的边界条件设置我一开始模仿了setExprFields的语法,使用condition语法限定分段区域,但是查看0/T文件发现后面的函数值会覆盖前面的赋值,也就是说condition语句在这里被忽略了,tutorials仅给出了2个示例,都比较简单。请问各位道友怎么实现这种功能??

  • levelset初始条件

    3 帖子
    190 浏览

    @fireztw 好的谢谢!

  • tecplot加载OpenFOAM文件报错且速度特别慢

    10 帖子
    1k 浏览

    @myheart 好的,谢谢您

  • Floating point exception

    8 帖子
    4k 浏览

    @Alvin 请问前辈最后有没有找到解决方法??

  • 液体与固体表面的接触角怎么定义和操作?

    4 帖子
    2k 浏览

    @wwzhao液体与固体表面的接触角怎么定义和操作? 中说:

    接触角是指气-液界面的切线穿过液体与固-液交界线之间的夹角θ。θ小于90°代表固体表面亲水,大于90°代表固体表面疏水。

    可以通过对wall类型的边界设置alpha边界条件来设置接触角。

    参考constantAlphaContactAngle、dynamicAlphaContactAngle以及timeVaryingAlphaContactAngle。

    想请教一下,如果像定义个类似于timeVaryingAlphaContactAngle的边界条件,只是把时间换成坐标x,y或者z,该如何创建新的边界条件呢?

  • 1 帖子
    70 浏览

    Openfoam中的多相流问题,如果考虑壁面润湿性,需要将alpha.water的wall设置为接触角边界条件,最常用的就是constantAlphaContactAngle和dynamicAlphaContactAngle。特别地,如果接触角随时间变化,有一个对应的timeVaryingAlphaContactAngle;现在想做一个接触角对空间位置变化的,如θ=θ1+(θ2-θ1)*x。也就是将timeVaryingAlphaContactAngle中的time改为x,该如何创建这一新的边界条件呢?求赐教

  • 高压下费托合成反应发散

    3 帖子
    155 浏览

    @bestucan 感谢您的指导,我去试着改一下

  • matllab数据导入paraview

    1 帖子
    85 浏览

    如何在matlab中将二维流场写入vtk文件,找了一圈都只有三维的,有大佬知道吗

  • 3 帖子
    185 浏览

    @李东岳OpenFOAM-dev 支持二维AMR吗?测试计算一直报错。 中说:

    不支持2维AMR,但是可以把z方向处理为slip

    原来如此,多谢李老师

  • Openfoam 远程并行计算问题

    4 帖子
    259 浏览

    可以试试screen,tmux,byobu

  • MPPICInterFoam支持自适应网格技术吗?

    1 帖子
    74 浏览

    各位大佬好!在MPPICInterFoam的源代码中没看到关于动网格的源文件,所以MPPICInterFoam求解器支持自适应网格技术吗?

  • 关于fvSchemes字典中wallDist的疑问

    4 帖子
    2k 浏览

    @tidedrinker关于fvSchemes字典中wallDist的疑问 中说:

    哪些模型最简单的方法是空置这一项,或者故意填个错的,然后报错的时候会提示你有那些可选
    确实k-omega SST需要,没有就会报错

  • 浮点数例外 (核心已转储)问题

    3 帖子
    201 浏览

    加密前网格由blockMesh生成blockMeshDict

    FoamFile { version 2.0; format ascii; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // scale 1; vertices ( (-0.15242 0 -0.005) ( 0 0 -0.005) ( 0.3048 0.081670913853 -0.005) (-0.15242 0.1524 -0.005) ( 0 0.1524 -0.005) ( 0.3048 0.1524 -0.005) (-0.15242 0 0.005) ( 0 0 0.005) ( 0.3048 0.081670913853 0.005) (-0.15242 0.1524 0.005) ( 0 0.1524 0.005) ( 0.3048 0.1524 0.005) ); blocks ( hex (0 1 4 3 6 7 10 9 ) (80 80 1) simpleGrading (1 1 1) hex (1 2 5 4 7 8 11 10) (160 80 1) simpleGrading (1 1 1) ); edges ( ); boundary ( inlet { type patch; faces ( (0 6 9 3) ); } outlet { type patch; faces ( (2 5 11 8) ); } bottom { type symmetryPlane; faces ( (0 1 7 6) ); } top { type symmetryPlane; faces ( (3 9 10 4) (4 10 11 5) ); } obstacle { type patch; faces ( (1 2 8 7) ); } ); mergePatchPairs ( ); // ************************************************************************* //

    加密后是由fluent.msh转换的fluent.msh
    加密前:controlDict
    加密后:controlDict
    其他设置基本是一样的
    物性参数:thermophysicalProperties
    湍流模型laminarturbulenceProperties
    离散格式:fvSchemes
    求解器控制:fvSolution
    来流条件:p T U

  • OpenFOAM中的sphereDrag相关参数修改

    3 帖子
    193 浏览

    @李东岳 谢谢李老师

  • paraview处理试验数据

    2 帖子
    160 浏览

    能否根据一些散点的坐标xyz,即这些散点的风压值,生成paraview可以读取的格式让其能插值显示曲面风压?

  • 关于setFields工具设置的一些疑问

    17 帖子
    3k 浏览

    @wangwang 谢谢前辈,是的发现问题了

  • OpenFOAM温度越界

    2 帖子
    164 浏览

    温度越界这个问题在可压缩那面出现的太多了。没有一个普适性的答案。这就是发散了。

  • openfoam波浪与结构相互作用,结构受力发散

    1 帖子
    110 浏览

    使用waves2foam造波,interfoam求解器。运行到后面时间步变得无限小,结构受力无限大。我这里有十多个构件单独输出受力,每个构件的受力全都变得很大。我觉得可能是网格质量不好的原因(网格用snappy生成),但是调整了很多参数,依旧会发散。以下是log文件、snappy文件(由于杆件太多,只展示部分杆件,所有杆件设置一致)以及网格局部图片,请问大家有没有调整的建议或者除了网格导致发散的其他原因。86a8b480-7d79-4850-86d4-6140dfede200-image.png 3db69d5e-0327-426e-a358-54a61e4a5a3a-image.png 25c1ee00-0922-4a62-bf33-b7c3741cc376-image.png f112809c-4b63-43d3-b173-d40aca264e39-image.png 22f99966-bf52-4050-8b3e-271e6bcfcb07-image.png db211ec2-f12a-4897-85af-b21f259fd754-image.png 6601262c-29d1-4f28-aaa7-ffd719e3f306-image.png 1287beb0-7bab-4333-b0af-a96fe9381bb0-image.png b6ae8373-35bf-4f95-b812-caa392da8516-image.png 65fcf0a7-e5cf-4748-8487-0d9cdfbeb344-image.png

  • OpenFOAM的两个燃烧求解器

    11 帖子
    4k 浏览

    @cccrrryyy 大佬,请教个问题,最近这两个求解器也困惑了我许久,拜读这个帖子后仍有疑问。
    问题1:按照李老师所说,声速较大会导致rhoReactingFoam中采用的压力方程无法求解,那是不是意味着reactingFoam适用于声速可能无穷大也就是低马赫数时密度不随压力改变的情况?
    问题2:同时由于rhoReactingFoam采用rhoPimpleFoam的压力方程,其为压力基求解器,捕捉激波的能力没有rhoCentralFoam强,所以其实际上也不适用于高马赫数流动?
    问题3:既然这俩都不能计算高速可压流,那基于psi和基于rho的热物性带来的差异主要体现在哪?

  • OF 中 颗粒的 sizedistribution models 问题?

    7 帖子
    4k 浏览

    @chpjz0391 请问这个multiNormal的文件在哪个tutorials里面?

  • openfoam set fields

    12 帖子
    9k 浏览

    关于应用setFields和codestream的组合定义初始场,可以参考大神的一个贴,https://www.topcfd.cn/12245/

  • rhoPimpleFoam在大网格下的发散问题

    8 帖子
    436 浏览

    绘图1.png
    李老师,上面算的算例因为网格问题要重新计算了,我参照您发我的算例在ICEM中画网格。如图所示A的严格对称点应该是B点,您画在了C点,所以要在openfoam里面用cyclicAMI来周期上下两个面。出现了面不匹配的老问题,我尝试在ICEM中设定周期,但是好像这种叶珊上下面应该没法设置节点的周期把?随后我仔细观察了您的算例,发现在boundary文件中的cyclicAMI的容差设置为1,随后我在createPatchDict采用了相同的容差。这一步算是可以通过了,但是在decomposePar中采用scotch方法进行分块时,还是在cyclicAMI这里报错了,麻烦李老师给看看我上述的操作是漏了什么嘛?
    另附一下 网格处理命令:
    fluentMeshToFoam ICEM.mesh
    autoPatch 45 -overwrite
    createPatch -overwrite

    屏幕截图 2023-12-30 235244.png

  • OpenFOAM第三方造波包waves2Foam安装

    5 帖子
    4k 浏览
  • 3 帖子
    214 浏览

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

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

    好的,谢谢李老师

  • divShemes离散格式的顺序问题

    3 帖子
    202 浏览

    这个就很奇怪。你要是能准备一个小算例,我可以看看,如果是bug,可以提bug。但是这个连bug都不像。像是你还是改了什么,然后忘记了

  • case外包,求助大佬们

    3 帖子
    254 浏览

    导师经费很多,可详谈

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

    8 帖子
    549 浏览

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

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

    32 帖子
    2k 浏览

    @李东岳 老师,请问现在关于时间步长无关性的研究有结论了吗?

  • 求助各位老师,谢谢

    1 帖子
    101 浏览

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

  • 《OpenFOAM用户指南》勘误

    已移动
    31 帖子
    13k 浏览

    @kuangyonglin 谢谢!下一版本更新

  • 1 帖子
    78 浏览

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

  • OvermultphaseInterDyMFoam

    7 帖子
    3k 浏览

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

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

    1 帖子
    104 浏览

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

  • 5 帖子
    640 浏览

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

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

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

  • interFoam计算气泡槽道流时的压力问题

    20 帖子
    1k 浏览

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

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

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

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

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

  • 自定义梯度边界条件

    10 帖子
    651 浏览

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

  • openfoam如何结合主动控制和被动控制?

    2 帖子
    179 浏览

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

  • outfile精度好像不够

    4 帖子
    852 浏览

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

  • Coupled level set-VOF方法

    4 帖子
    1k 浏览

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

  • 6 帖子
    399 浏览

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

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

    1 帖子
    63 浏览

    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怎么解决呀 在网上参考了各种方法 都无法成功
  • openfoam中有没有LS方法模拟多相流的案例

    1 帖子
    90 浏览

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

  • OpenFOAM自适应网格和并行计算冲突的问题

    13 帖子
    3k 浏览

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

    image.png

  • 1 帖子
    96 浏览

    如上所述