@李东岳 老师，MRF区域应该没啥问题，沿Y轴旋转，650rmp。8e9a8d80-2f54-4d1b-b6f3-37f0d071d10d-image.png
然后速度场看上去也是在转的：20cb61fd-47da-449b-bc33-cbcfbfbd63b0-image.png

那个就是发散了。

把log贴出来看一下

@zhangK 客气，祝顺利！

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

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

:136:

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

我尝试了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

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

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

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

速度方程不需要solve

http://dyfluid.com/icoFoam.html

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

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

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

需要挂参数 -region fluid类似

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

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

@allanZhong 材料非常好，我看看

你要求的是标量phi对吧，它收敛了么？

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

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

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

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

算垂直平板绕流势流，入流速度为1、压力接近0（出口压力为0）。前驻点压力应该为1吧？咋是1.427？
微信图片_20240415110333.jpg
手机拍的图，搞不会paraview显示点数据时咋截图 →_→

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

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

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

@huiyin 谢谢~

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

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

打开试了试，步长极小，残差也一直保持在一个值不下降，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.2）的情况，在设置重叠部分也即船舶部分网格时，blockmesh需要画的比较小，底部与船舶底部很接近，这就造成了snappyhexmesh后船舶下部的网格质量很差，请问有没有什么办法或者思路来改善网格质量呢，如有相似问题的朋友欢迎讨论，谢谢!261e3ad06976a58e400a1fc15243dc1.png 78823dc1464db8ab58735e54dcd17d2.png efa88c4c50ed4254badb8dbc8db3089.png

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

我也在算类似的问题，是不是给定进口水位和流量的设置方法不对

用interFoam计算带自由液面的倒虹吸，参考其他人边界设置，往往是入口给定流量（流速），出口给定水位（alpha.air)。这种设置方式如果通过手动设置（setFields)往往因初始流场不好极易发散，有没有办法能跟Fluent一样通过混合初始化得到一个初始场保证计算能够进行下去。（openFoam好像可以通过potentialFoam初始化场，但是这个对于interFoam也可以使用吗）79538e5a-b907-420a-8777-c7a72dd6011f-image.png

@cccrrryyy谢谢前辈

你好，请问你编写的NSRFG，在使用时计算域是多大，网格尺度是多大，能否在整个计算域上保持低耗散甚至无耗散么

参考之前大神发的帖子就可以解决了
https://www.cfd-china.com/topic/6750/mapfield报错

@jxzs_zy 在开发，好像还没正式发布。Added dynamicPolyRefinementFvMesh

问题找到了，对于这个算例，在滤波尺度的选择中，不能够使用剪切层调整亚格子长度尺度修正，但是别的文献中有人使用过，目前没有找到原因。

你好，请问这个问题解决了吗？

您好，您解决了吗

您好，您解决了吗

@BznW 您好。打扰一下，您解决了吗

前一段时间搭建了流化床求解器（燃烧、密相颗粒流动求解器，下称newSolver），目前需要对这个求解器进行冷态颗粒流动的验证（常温，关闭化学反应、燃烧和传热）。通过验证发现颗粒轴向速度和实验吻合较好，但是气相压降差距很大，而使用MPPICFoam计算相同的case可以得到准确的压降（稳定时newSolver比MPPICFoam压降高两倍多）。

为了找到问题所在，于是关闭颗粒计算纯气体流动的速度和压降，发现newSolver和MPPICFoam计算得到的气相速度相同但是压降有一定差距。因此猜测是气相压力方程计算的问题。

由于newSolver是基于coalChemistryFoam流体方程加入alpha的修改的，纯气体流动情况下alpha=1。为了排除我修改方程所导致的错误，因此用coalChemistryFoam（关闭化学反应、燃烧和传热）计算纯气体流动看能否得到准确的结果。结果发现压降仍然有差距。

难道coalChemisryFoam变密度的求解器无法计算常密度的情况吗？coalChemistryFoam是反应所导致的密度变化，关闭化学反应按理说密度应该是不变的，为何计算结果与实际差距大？常温纯气体流动压降都不准，加上化学反应燃烧后压降会准吗？

MPPICFoam压力方程

fvScalarMatrix pEqn ( fvm::laplacian(alphacf*rAUcf, p) == fvc::ddt(alphac) + fvc::div(alphacf*phiHbyA) );

coalChemistryFoam压力方程

fvScalarMatrix pEqn ( fvm::ddt(psi, p) + fvc::div(phiHbyA) - fvm::laplacian(rhorAUf, p) == coalParcels.Srho() + fvOptions(psi, p, rho.name()) );

纯气体流动验证的几何模型和网格如下图所示
图片1.png

请各位大佬指教！

这段用来统计网格单元内粒子个数的代码，编译的路径是什么，求东岳老师指导

这是我用flunet做的多相流模拟，把一个水喷向液态铅铋，如果想用opemFOAM,请问用什么分割两相间的分界面，使用cyclic边界条件吗。5baa602a-c4ca-473c-925a-631e18c31623-image.png

@sxz0823 如果您方便的话，可以加个联系方式嘛。v:13920535789/qq:2660317708。

基于interFoam求解器，添加了标量场S，在模拟了一段时间后，得到了两相流体分布，再使用swak4foam中的funkySetFields，把标量赋值到alpha1流体中，另一相不考虑，设置如下：
b8c93503-efa8-4a6d-9b18-7852dcb2d072-image.png
运行后报错:
15a01aae-36ee-42e4-8223-4fc9eba68dd9-image.png
请问各位大神该如何解决呢？

@liuyf 我在求解时候也遇到了这个问题，导致鼓泡塔的液面比实验高很多，因为上边部分存在液体速度导致求解一开始液体会一下子冲上去，最后液面高度就会变得非物理，您解决这个问题了吗

g++ -std=c++11 -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -O3 -DNoRepository -ftemplate-depth-100 -I/usr/local/include -I./windProfile/lnInclude -I./extrudedPatch/lnInclude -I/home/windlab_lzx/openfoam8/OpenFOAM-8/src/finiteVolume/lnInclude -I/home/windlab_lzx/openfoam8/OpenFOAM-8/src/meshTools/lnInclude -I/home/windlab_lzx/openfoam8/OpenFOAM-8/src/sampling/lnInclude -I/home/windlab_lzx/openfoam8/OpenFOAM-8/src/surfMesh/lnInclude -I/opt/intel/oneapi/mkl/2024.1/include -IlnInclude -I. -I/home/windlab_lzx/openfoam8/OpenFOAM-8/src/OpenFOAM/lnInclude -I/home/windlab_lzx/openfoam8/OpenFOAM-8/src/OSspecific/POSIX/lnInclude -fPIC -fuse-ld=bfd -Xlinker --add-needed -Xlinker --no-as-needed /home/windlab_lzx/openfoam8/OpenFOAM-8/platforms/linux64GccDPInt32Opt/applications/utilities/preProcessing/DFSRTurb/DFSRTurb.o -L/home/windlab_lzx/openfoam8/OpenFOAM-8/platforms/linux64GccDPInt32Opt/lib \ -L/usr/local/include -lwindProfile -lextrudedPatch -lmeshTools -lsampling -lsurfMesh -lfiniteVolume -lfftw3 -L/opt/intel/oneapi/mkl/2024.1/lib/intel64 -lmkl_intel_lp64 -L/opt/intel/oneapi/mkl/2024.1/lib/intel64 -lmkl_core -L/opt/intel/oneapi/mkl/2024.1/lib/intel64 -lmkl_sequential -lm -ldl -lOpenFOAM -ldl \ -lm -o /home/windlab_lzx/OpenFOAM/windlab_lzx-8/platforms/linux64GccDPInt32Opt/bin/DFSRTurb /usr/bin/ld.bfd: 找不到 -lwindProfile /usr/bin/ld.bfd: 找不到 -lextrudedPatch collect2: error: ld returned 1 exit status /home/windlab_lzx/openfoam8/OpenFOAM-8/wmake/makefiles/general:140: recipe for target '/home/windlab_lzx/OpenFOAM/windlab_lzx-8/platforms/linux64GccDPInt32Opt/bin/DFSRTurb' failed make: *** [/home/windlab_lzx/OpenFOAM/windlab_lzx-8/platforms/lDFSRinux64GccDPInt32Opt/bin/DFSRTurb] Error 1在编译DFSR方法最后一步wmake时候报错，为什么呀，找了很久没解决

@shirouchen 谢谢，我暂时是用matlab排列的