OpenFOAM

4.5k 主题 26.4k 帖子
  • 合集:OpenFOAM安装准备工作贴

    已固定
    7 帖子
    15k 浏览

    @sdugem 已解决

  • OpenFOAM教程 初级入门建议(2020年更新版)

    已固定
    47 帖子
    97k 浏览

    @黑山雁 处理了 :-)

  • 11 帖子
    54 浏览

    supersonic_injet-stead.tar.xz

    更新在此: http://dyfluid.com/openfoam.html

    对了,我没加湍流,湍流你自己加一下

  • 振荡流入口压力边界条件相关问题

    7 帖子
    26 浏览

    @李东岳 嗯嗯,好的,谢谢老师

  • 5 帖子
    19 浏览

    @lwjetmann 感谢老师,组里同学大家都不太会OpenFOAM,我马上也看一下。再次感谢。

  • icoFoam中压力泊松方程的求解

    9 帖子
    1k 浏览

    更新到了5.9.1节

  • potentialFoam计算的驻点压力系数为啥不是1?

    8 帖子
    97 浏览

    @李东岳 我用的v2006,Phi是速度势,用-grad(Phi)计算的速度和输出的U有些差异,尤其是在远场边界。

  • 60 帖子
    16k 浏览

    @lizhisongsjtu 这个图太酷炫了。 有doi不

    有些大佬对liutex法的态度偏消极(审稿意见)。

    是,我也看到了。为什么呢

  • OpenFOAM 浮体铰接可以计算吗

    7 帖子
    655 浏览

    @YuPeng 我也遇到了相同的问题,请问一下您现在解决了吗

  • fvSchemes中interpolationSchemes设置

    3 帖子
    28 浏览

    感谢李老师!经过修改,已经可以成功运行了。:146:

  • openfoam添加湍流入口方法?

    46 帖子
    14k 浏览

    请问在使用turbulentDFSEMInlet入口类型计算过程中,随着计算的进行,后半部分的风场慢慢的变得不是对数廓线了,请问是什么原因,有没有大佬知道。
    我的初始边界信息:
    计算域:xyz为1011(m)
    U:

    inlet { type turbulentDFSEMInlet; delta 1; U { type mappedFile; mapMethod nearest; } R uniform (0.8 0 0 0.4 0 0.4); L uniform 0.004; d 1; nCellPerEddy 1; scale 1; value $internalField; } outlet { type advective; value $internalField; } "(topWall|front|back)" { type zeroGradient; } "(bottomWall|Model)" { type noSlip; }

    输入文件:
    point:

    ( ( 0 0 0.01 ) ( 0 0 0.05 ) ( 0 0 0.1 ) ( 0 0 0.15 ) ( 0 0 0.2 ) ( 0 0 0.25 ) ( 0 0 0.3 ) ( 0 0 0.35 ) ( 0 0 0.4 ) ( 0 0 0.45 ) ( 0 0 0.5 ) ( 0 0 0.55 ) ( 0 0 0.6 ) ( 0 0 0.65 ) ( 0 0 0.7 ) ( 0 0 0.75 ) ( 0 0 0.8 ) ( 0 0 0.85 ) ( 0 0 0.9 ) ( 0 0 0.95 ) ( 0 0 0.99 ) )

    U:

    ( ( 0 0 0 ) ( 1.048097372 0 0 ) ( 1.499488341 0 0 ) ( 1.763535131 0 0 ) ( 1.95087931 0 0 ) ( 2.096194744 0 0 ) ( 2.2149261 0 0 ) ( 2.31531203 0 0 ) ( 2.402270279 0 0 ) ( 2.47897289 0 0 ) ( 2.547585713 0 0 ) ( 2.609653562 0 0 ) ( 2.666317069 0 0 ) ( 2.718442442 0 0 ) ( 2.766702999 0 0 ) ( 2.811632503 0 0 ) ( 2.853661247 0 0 ) ( 2.893141184 0 0 ) ( 2.930363858 0 0 ) ( 2.965573488 0 0 ) ( 2.992431707 0 0 ) )

    计算3600s后,平均风速廓线,U1U2U3U4U5U6U7U8U9的位置为0.1,2,4,5,6,7,8,9,9.9
    e595f49f-347a-4180-82e2-b6295cd45a15-图片.png

  • 升阻力系数计算问题

    3 帖子
    40 浏览

    lref设置为直径D=1,Aref=D×径向长度,改成2.68328就ok了,你使用什么软件导入的网格,fluent.mesh导入的时候会告诉你径向拉伸了多少,你再尝试一下吧https://blog.csdn.net/weixin_39124457/article/details/120187874
    你可以看看这个网站设置的东西

  • 5 帖子
    116 浏览

    @李东岳 就是单核跑也不行,现在解决了,应该是openfoam10里面的pitzDaily哪里要调一下,换成其他算例就能跑了,谢谢东哥

  • 后向台阶流动的入口速度给定问题

    1 帖子
    44 浏览

    我目前想要计算一个后向台阶流动问题,实验中台阶上方为充分发展的湍流边界层,因此入口速度需要额外的处理。目前认为最准确的方法是单独计算一个槽道流的充分发展湍流并提取速度平面作为入口条件。然而这两个计算需要满足对应关系。Edmond Shehadi博士的博士“Large Eddy Simulation of Turbulent Flow over a Backward-Facing Step”中给出了通过后向台阶流动动量耗散雷诺数和台阶高度雷诺数得到对应槽道流参数的流程,然而根据他们的流程,我计算他们的台阶流参数得到的槽道流摩擦雷诺数为730.7073,他们最终的结果为388。我检查过公式没有问题,这说明有额外的变化在这里并没有交代,想请问各位老师和同学能否指点迷津,非常感谢。下方是论文中给出的换算流程。槽道流的半高他们说应该取为台阶流中边界层厚度δ。
    屏幕截图 2024-04-20 212245.png
    屏幕截图 2024-04-20 212331.png

  • Openfoam使用SST k-omega湍流模型相关问题。

    1 帖子
    50 浏览

    我在做射流模拟的时候,使用了SST k-omega湍流模型,控制了y+~1,但模拟出来的边界层和自由剪切层与实验相比厚了两倍,没有使用壁面函数,这种情况下该从哪个方面着手去解决呢?湍流边界条件omega是根据文献中建议的值去给定的。

  • 4 帖子
    101 浏览

    root文件夹是不是没空间了?有一些系统root单独空间划片了。我看你跑在了root下。操作不是很常规

  • openfoam fvDOM中的太阳辐射模型

    2 帖子
    128 浏览

    这个问题解决了,主要是因为OpenFOAM在solarCalculator.C中有出纰漏
    修改如下并编译就好了,可以模拟在useExternalBeam下给定太阳辐射值(随时间变化)的效果

    void Foam::solarCalculator::correctSunDirection() { if (sunDirectionModel_ == mSunDirTracking) { calculateBetaTheta(); calculateSunDirection(); //directSolarRad_ = A_/exp(B_/sin(max(beta_, ROOTVSMALL))); if (sunLoadModel_ == mSunLoadTimeDependent) { directSolarRad_ = directSolarRads_->value(mesh_.time().value()); } else if (sunLoadModel_ == mSunLoadFairWeatherConditions) { directSolarRad_ = A_/exp(B_/sin(max(beta_, ROOTVSMALL))); } } } void Foam::solarCalculator::correctDirectSolarRad() { if (sunLoadModel_ == mSunLoadTimeDependent) { directSolarRad_ = directSolarRads_->value(mesh_.time().value()); } else if (sunLoadModel_ == mSunLoadFairWeatherConditions) { directSolarRad_ = A_/exp(B_/sin(max(beta_, ROOTVSMALL))); } }
  • 5 帖子
    111 浏览

    @李东岳 好的,李老师,那我自己再修改修改定义步长的代码,感谢李老师,有什么参考的资料吗qaq

  • openfoam如何设置圆形相域?

    9 帖子
    3k 浏览

    半圆柱没有

  • 3 帖子
    99 浏览

    @李东岳 谢谢李老师的回答,我再找找原因!

  • Tobias书中关于pimple的final松弛因子的设定

    3 帖子
    512 浏览

    @沐风 你上传的作者联系我删除你上传的pdf文件,我将在明天将其进行删除

  • 搅拌桨叶附近压力异常

    18 帖子
    296 浏览

    @知行 ok,我的网格都是ICEM画的,我加密试试,谢谢

  • 多相流interface

    2 帖子
    154 浏览

    两相界面的处理是用contour,然后切alpha=0.5的分界面

  • interFoam的问题

    16 帖子
    9k 浏览

    那个就是发散了。

  • 请教一下关于OpenFOAM计算速度的问题

    2 帖子
    105 浏览

    把log贴出来看一下

  • interFoam 模拟水滴运动

    9 帖子
    991 浏览

    @zhangK 客气,祝顺利!

  • snappyHexMesh画网格边界层

    3 帖子
    166 浏览

    @李东岳 好滴,谢谢李老师~

    发现这个参数匹配好麻烦啊,层数、第一层高度等等之类的,感觉按照计算得到的画一下基本不出来

    :136:

  • OpenFoam 5.x 计算颗粒是不是有bug?

    9 帖子
    3k 浏览

    CFDEM3.8及OF5.X是用这个reconstructPar -noLagrangian,或者直接在solver中把postpro改成true

  • olaFlow 造波为什么差别这么大

    1 帖子
    173 浏览

    我尝试了olaFlow 来进行造波,水池深0.5米,波高0.024,周期1.26s,水池长12米,以下是各个位置的波高数据:
    x= 0.png
    x= 1.png
    x= 2.png
    x= 3.png
    x= 4.png
    x= 5.png
    x= 6.png
    x= 7.png
    x= 8.png
    x= 9.png

  • 200万网格并行算力测试(OpenFOAM版本)

    199 帖子
    213k 浏览
    CPU型号: Intel(R) Xeon(R) Platinum 8475B 系统:ubuntu 22.04 版本:openfoam10 内存:128G 16G 3200MHz*8 48 51 32 60 16 91 8 117 4 226 2 410 1 848
  • 有偿!江湖救急

    6 帖子
    1k 浏览

    @liujunCFD 最后用matlab怎么解决的呢,我目前也遇到这个问题了,没人带,只能自己摸索,所以大佬可以告知一下方法吗?

  • 5 帖子
    148 浏览

    @李东岳 老师,我有一个问题,压力梯度是指沿着流体流动的方向压力的变化,如果我的流体入口与水平面有45°的夹角,但我的水是沿着x方向进入的,即只有x方向的速度,我的压力梯度应该是把x方向的压力梯度转换到面法向方向上去么?

  • TwoPhaseEularFoam

    2 帖子
    150 浏览

    速度方程不需要solve

    http://dyfluid.com/icoFoam.html

  • 偏导数的问题

    2 帖子
    143 浏览

    你那个第一\partial (rho f0)/\partial (rho) 是不是可以继续展开?

  • 自然对流温度场振荡(异常?)的问题

    2 帖子
    226 浏览

    纯粹的震荡了。看不出其他问题。

  • 2 帖子
    90 浏览

    找了很多手册和视频,貌似都没有解决方法。也有想到用localfoamtotecplot360转换,但也没有找到有v2312适用的版本

  • paraFoam读取不到场信息

    2 帖子
    105 浏览

    需要挂参数 -region fluid类似

  • OpenFOAM第三方造波包waves2Foam安装

    7 帖子
    4k 浏览

    https://www.cfd-china.com/topic/6835/waves2foam安装-2024年更新版本/6

    这个链接里面的wavesfoam里面包含了oceanwave3D

  • CFD与机器学习基本概念入门

    14 帖子
    1k 浏览

    @allanZhong 材料非常好,我看看

  • 2 帖子
    108 浏览

    你要求的是标量phi对吧,它收敛了么?

    我是建议把源项去掉,只保留对流、扩散和时间项,看看标量的输运是否合理。

    印象中Fluent不会对uds做什么限制,它也不好限制,因为它不清楚要加的标量是什么量级。

  • 请问如何得到颗粒迹线?

    11 帖子
    2k 浏览

    @calrel 您好,您现在解决这个问题了吗

  • openfoam中的重叠网格浮体位移的提取

    2 帖子
    226 浏览

    你好,我也在做这一项工作,请问您现在寻找到解决该问题的办法了吗?

  • FPVFoam

    1 帖子
    71 浏览

    请问各位大佬,OpenFOAM v1912中怎么植入FPVFoam(新手)

  • blockmesh生成网格总报错没有相邻面

    9 帖子
    983 浏览

    @liurui 面的编号要用右手定则指向计算域外侧,所以就比如第一个hex应该是(0 7 6 1 8 9 14 15)

  • snappyHexMesh网格贴合问题

    2 帖子
    287 浏览

    您好,请问解决了嘛?我现在也碰到了同样的问

  • snappyHexMesh以及paraview显示问题

    12 帖子
    3k 浏览

    @huiyin 谢谢~

  • interFoam吐核

    3 帖子
    195 浏览

    @学流体的小明 好的,谢谢

  • 3 帖子
    232 浏览

    @李东岳 感谢老师回复!动量预测是关掉了的。

    打开试了试,步长极小,残差也一直保持在一个值不下降,cumulative则很大,如下:

    Courant Number mean: 7.953414508e-05 max: 0.4993036972 deltaT = 4.979675437e-35 Time = 6.00001e-05 PIMPLE: iteration 1 DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 0.04941709914, No Iterations 29 DICPCG: Solving for pcorr, Initial residual = 0.03122316034, Final residual = 0.001556737496, No Iterations 121 DICPCG: Solving for pcorr, Initial residual = 0.006804280077, Final residual = 0.0003314734696, No Iterations 359 DICPCG: Solving for pcorr, Initial residual = 0.002308369425, Final residual = 9.709852492e-11, No Iterations 911 time step continuity errors : sum local = 1.793718362e-20, global = 1.178518414e-23, cumulative = 1.037419916e+19 smoothSolver: Solving for Ux, Initial residual = 3.985498772e-05, Final residual = 2.227155701e-11, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 6.302253167e-05, Final residual = 2.385619185e-11, No Iterations 1 GAMG: Solving for p, Initial residual = 0.004038713667, Final residual = 0.0001722133227, No Iterations 18 GAMG: Solving for p, Initial residual = 0.0005583071565, Final residual = 2.673668049e-05, No Iterations 34 GAMG: Solving for p, Initial residual = 9.485975434e-05, Final residual = 4.53317788e-06, No Iterations 33 GAMG: Solving for p, Initial residual = 2.596708738e-05, Final residual = 1.231466526e-06, No Iterations 38 time step continuity errors : sum local = 7.929103562e-15, global = -8.302368146e-16, cumulative = 1.037419916e+19 GAMG: Solving for p, Initial residual = 0.003717218882, Final residual = 0.000164557501, No Iterations 17 GAMG: Solving for p, Initial residual = 0.0005378364039, Final residual = 2.53035408e-05, No Iterations 34 GAMG: Solving for p, Initial residual = 9.265726415e-05, Final residual = 4.535601234e-06, No Iterations 36 GAMG: Solving for p, Initial residual = 2.556961887e-05, Final residual = 9.918214454e-11, No Iterations 548 time step continuity errors : sum local = 6.482264028e-19, global = 6.603865839e-20, cumulative = 1.037419916e+19 PIMPLE: iteration 2 ... PIMPLE: iteration 50 smoothSolver: Solving for Ux, Initial residual = 1.226945138e-10, Final residual = 6.6691269e-16, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 1.651845516e-10, Final residual = 3.177048458e-16, No Iterations 1 GAMG: Solving for p, Initial residual = 0.000949721577, Final residual = 4.516299707e-05, No Iterations 41 GAMG: Solving for p, Initial residual = 0.0004457930337, Final residual = 2.206483819e-05, No Iterations 23 GAMG: Solving for p, Initial residual = 7.977007558e-05, Final residual = 3.747481321e-06, No Iterations 36 GAMG: Solving for p, Initial residual = 2.418915618e-05, Final residual = 1.162060008e-06, No Iterations 48 time step continuity errors : sum local = 7.482194647e-15, global = -7.221446192e-16, cumulative = 1.037419916e+19 GAMG: Solving for p, Initial residual = 0.003705959955, Final residual = 0.0001646823634, No Iterations 17 GAMG: Solving for p, Initial residual = 0.0005377574684, Final residual = 2.535476181e-05, No Iterations 34 GAMG: Solving for p, Initial residual = 9.241219827e-05, Final residual = 4.525895839e-06, No Iterations 36 GAMG: Solving for p, Initial residual = 2.542969083e-05, Final residual = 9.774400044e-11, No Iterations 548 time step continuity errors : sum local = 6.388632178e-19, global = 6.502735317e-20, cumulative = 1.037419916e+19 PIMPLE: not converged within 50 iterations ExecutionTime = 21923.2 s ClockTime = 21963 s Courant Number mean: 7.964991508e-05 max: 0.4989139005 deltaT = 4.99051583e-35 Time = 6.00001e-05

    看到站内其它的帖子,试了试有量纲/无量纲计算(此前圆柱直径0.05,nu=50;改为圆柱直径0.0005,nu=0.5),也无法收敛。

    走投无路试了一下StarCCM,多低雷诺数都能算,和OpenFOAM那些能收敛的算例比得还贼好

  • 限制区域内提高重叠网格质量

    1 帖子
    71 浏览

    各位好,我现在设想一个工况是在水深非常浅,同时水深与船舶吃水非常小(大约1.2)的情况,在设置重叠部分也即船舶部分网格时,blockmesh需要画的比较小,底部与船舶底部很接近,这就造成了snappyhexmesh后船舶下部的网格质量很差,请问有没有什么办法或者思路来改善网格质量呢,如有相似问题的朋友欢迎讨论,谢谢!261e3ad06976a58e400a1fc15243dc1.png 78823dc1464db8ab58735e54dcd17d2.png efa88c4c50ed4254badb8dbc8db3089.png

  • overinterdymfoam的一些疑问

    6 帖子
    284 浏览

    @学流体的小明 好的,非常感谢解答