最后还是用了手动赋值
a9878440-a6cc-4560-a822-8956f3fbc06d-image.png

@东岳 谢谢老师建议! 回头再遇到类似的问题, 我多看看报错, 这样子更有针对性

@bestucan 谢谢回复。我看了一下，OpenFOAM7中是指定一个点的，通过location；OpenFOAM1912中有一个option可选的是locations，可指定多个点

@bestucan

intermediateNew$ wmake libso wmake libso . 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 -O3 -DNoRepository -ftemplate-depth-100 -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/lagrangian/basic/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/lagrangian/distributionModels/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/transportModels/compressible/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/specie/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/basic/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/thermophysicalProperties/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/reactionThermo/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/SLGThermo/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/radiation/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/transportModels -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/regionModels/regionModel/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/regionModels/surfaceFilmModels/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/dynamicFvMesh/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/sampling/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/finiteVolume/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/meshTools/lnInclude -IlnInclude -I. -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/OpenFOAM/lnInclude -I/home/hoyu/OpenFOAM/OpenFOAM-5.x/src/OSspecific/POSIX/lnInclude -fPIC -c parcels/derived/basicKinematicParcel/defineBasicKinematicParcel.C -o Make/linux64GccDPInt32Opt/parcels/derived/basicKinematicParcel/defineBasicKinematicParcel.o In file included from lnInclude/KinematicParcel.H:692:0, from parcels/derived/basicKinematicParcel/basicKinematicParcel.H:40, from parcels/derived/basicKinematicParcel/defineBasicKinematicParcel.C:26: lnInclude/KinematicParcelI.H:127:8: error: prototype for ‘Foam::KinematicParcel<ParcelType>::KinematicParcel(const Foam::polyMesh&, const barycentric&, Foam::label, Foam::label, Foam::label, Foam::label, Foam::scalar, Foam::scalar, Foam::scalar, const vector&, const vector&, const Foam::KinematicParcel<ParcelType>::constantProperties&)’ does not match any in class ‘Foam::KinematicParcel<ParcelType>’ inline Foam::KinematicParcel<ParcelType>::KinematicParcel

sorry我看眼花了，多加了个点这下输出真正的编译错误了

@东岳 谢谢东岳老师！:xiexie:

想在cavity算例模型内部加入一单位大小的矩形悬浮块，想观察它随着周围流态变化的运动过程，请问动网格这块应该设置什么类型，或者of有什么算例能参照下，正在初学动网格，好些不懂。

@bestucan 好的，我试一下，谢谢前辈

@bestucan 好的谢谢。感觉你说的挺有道理的。之前尾部网格是很扁的矩形（因为存在边界层），结果微小的幅度旋转都会使网格拉扯出现负体积，所以我才逐步减小尾部网格的横纵比，让网格看起来不会那么扁，成功实现了8°范围的俯仰，然后为了做更大的俯仰角度，继续减小网格横纵比，试图可以让他多转一些角度，当网格的横纵比减小到现在的情况时，并没有出现网格畸形，也出现了发散，就不知道是为什么了。

@星星星星晴 可以留个联系方式

@cccrrryyy 非常感谢~

@东岳 好的，谢谢李老师。今天早上破案了，OpenFOAM7成功编译了。

@东岳 感谢李老师！！是我边界条件设的不合理加上网格密度取值的问题。已经更正。谢谢李老师耐心解答！！不可压的alpha就是groMetricOneField，我想太简单了，以为稳流模型把空气自动过滤掉了。

@东岳 好的，谢谢老师～～那主要还是从模拟结果上分析最为有效吧～～

@东岳 好的东岳老师，谢谢东岳老师，实在不好意思，上次断网没发出去刚发现。

@东岳 对的，我也正在琢磨在cfd-online和stackoverflow上推一下的，此插件支持双语版本，根据用户的vscode界面语言选择使用中文还是英文显示

@东岳 谢谢东岳老师！:ok3:

@东岳 顶部#include已加

我找着我第一次发的方程改得，就是第二个block的代码，算出来的结果和MPPICInterfoam的基本一样。当然有可能我的粒子直径比较小，对计算结果影响不大？

我没明白您下面写的这个是什么意思。如果不看mrf的话， 你这个和interfoam相比就是减上了一项

fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)

这一项就是65293cb6-a527-49eb-bdc4-2355ec82c507-image.png
我不知道怎么推出的这一项？

https://github.com/OpenFOAM/OpenFOAM-dev/blob/master/src/OpenFOAM/fields/GeometricFields/GeometricField/GeometricBoundaryField.C

这是边界场

@星星星星晴 谢谢！！