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排列的

大家好！想请教各位在tecplot中对三维模型进行Mirror操作后，使用slices提取出一个面，可以观察到提取出来的面云图都还是对称的,但画出来的流线为什么不是对称的了？原计算案例对称轴所在的面本身使用的也是symmetryPlane型边界条件。
QQ图片20240407133422.png
$DR1%SP2~_Y1_A{0ZO8ZND8.png
TUG[A$3E@2}8%R85`77T]3K.png

那我不知道咋办了，要不你用apt删了重装一下，或者编译源代码安装。
或者你电脑提前改变了apt的安装路径，安装到了其它地方？
下面这个里面安装的openfoam9就是在/opt目录下面。
https://blog.csdn.net/sagjhdj/article/details/123435344

@学流体的小明
1、谢谢大佬！就是这个原因，mesh.Cf()改成mesh.C()就没问题了！
2、当时注释这个是因为我发现有问题后用Info输出坐标和mesh的size想看看问题出在哪，但是只输出了processer0的网格信息，所以有了这个疑问。后来我查到了是Info只会输出master processor的信息，我改成Pout就可以输出所有processor的信息。

最近在做一个模型，有一个面既按照一定的角速度旋转，同时又有流体射出，但是不知道该怎么设置这个面的边界条件，有大佬知道吗

@尚善若水 问题解决了，多谢多谢

已解决。

这么做尖角网格是不是看起来不太对呀，这类尖角网格有没有什么好一点的做法呢？
虽然看着网格质量还行，但总感觉不太对
Snipaste_2024-04-03_21-48-45.png

在轴对称模型中，可以只计算两个方向（轴向和径向）的量吗，在写方程时，如何不计算第三个方向的偏导数，可以像volTensorField一样定义一个Tensor2D的体积场吗

@coolhhh 太细致了！！！

@huangyuhui723 我明白你的意思了。

@huangyuhui723 在 关于在求解器中使用热物理库中的函数 中说：

通过其它计算更新了T之后想对he赋值

如果是求解TEqn，可以是得到温度T后更新he。OpenFOAM里的求解器一般都是求解EEqn，得到he后通过thermo.correct()更新T。

如果已知边界温度T后要对he赋值，应该在求解EEqn之前对he边界赋值，然后求解EEqn。大致代码是这样：

volScalarField& he = thermo.he(); //****************************************//加这一段 forAll(p.boundaryField(), patchi) { forAll(p.boundaryField()[patchi], facei) { he.boundaryFieldRef()[patchi][facei] = function(p.boundaryField()[patchi][facei], T.boundaryField()[patchi][facei]); } } //****************************************// fvScalarMatrix EEqn ( ...... );

这个function函数就是T与he的关系，类似上面，混合物的he通过各组分的显焓的质量分数加权得到，显焓用janaf里的关系确定。

供参考。

@学流体的小明 感谢老师，我找时间试试看

我在做单根不锈钢丝切割液滴的模拟，发现液滴触碰到固体后，会快速铺展覆盖固体，这明显与实际情况不符：
dc6202db-0a4b-4bee-8dfd-d8db7a8ba44e-image.png
9f370870-b0be-4f65-9075-57dd718ce5cd-image.png
求解器用的v2212的MPPICInterFOAM,固体壁面的边界条件为：U noSlip；P fixedFluxPressure；Alpha constantAlphaContactAngle，接触角为155疏水界面

文献和实际可视化拍摄中，液体不会有如此快速的铺展，还请各位指点一二
c833896f-62e9-4161-ae81-1b226e382477-image.png

@尚善若水 感谢，我去试试

1711709658497.png阿里云的VFC链接可以安装图形界面，按照这个教程来就可以，非常简单https://help.aliyun.com/zh/ecs/user-guide/install-a-gui-desktop-for-ubuntu-18

可以从边界上的face返回这个face属于的cell

const Foam::fvBoundaryMesh &procBoundary = mesh.boundary(); forAll(procBoundary, patch) { forAll(procBoundary[patch], facei) { const label bcell = boundaryMesh[patch].faceCells()[facei]; } }

你可以得到所有的边界单元编号，然后做你其它向做的事就应该比较方便了

@pengdi 没有收缩段，RANS结果就是很差的，LES即使入口只加平均风，也比RANS结果精确很多

已经搞清楚了，谢谢！

openfoam自带有interfaceHeight的function