您好，您解决了吗

@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

请各位大佬指教！

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

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

大家好！想请教各位在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的信息。

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

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

已解决。

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

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

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

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]; } }

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

已经搞清楚了，谢谢！