本算例采用的是气液两相流，想要给空气相加一个速度边界条件，而水相不加这个边界条件，下面给出的这个边界条件是空气相的，是依时的边界条件，按照网站上相似的做法编了以下，但是始终出现如下错误。请各位老师同学有时间帮忙解答一下，非常感谢！
下面是边界条件：

afInlet { type codedFixedValue; value $internalField; redirectType ramp; name sinewave; code #{ #include <math.h> scalar t = this->db().time().value(); fvPatchField<vector> u ( patch().lookupPatchField<volVectorField,vector>("U.air") ); u[t] = vector(0,0,2*3.1416*10/1000/0.77*1*sin(2*3.1416*1*t)); (*this) == u; #};

下面是报错的内容：

--> FOAM FATAL ERROR: request for volVectorField U from objectRegistry region0 failed available objects of type volVectorField are 8 ( HbyA.air U.air U.water U.air_0 KdUByA.air U.water_0 KdUByA.water HbyA.water ) From function const Type& Foam::objectRegistry::lookupObject(const Foam::word&) const [with Type = Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>] in file /home/dyfluid/OpenFOAM/OpenFOAM-8/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 211. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const& Foam::objectRegistry::lookupObject<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> >(Foam::word const&) const in "/home/dyfluid/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam" #3 Foam::totalPressureFvPatchScalarField::updateCoeffs() at ??:? #4 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/home/dyfluid/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam" #5 ? in "/home/dyfluid/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam" #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #7 ? in "/home/dyfluid/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam" Aborted (core dumped)

@number44 在 foam-extend的fsiFoam并行效率问题 中说：

foam-extend的icoFsiFoam

您好，请问下，下面两者有什么区别呢？
foam-extend中自带的icoFsiFoam
跟
单独的fsiFoam

@samuel-tu 这个可以有。如果可以的话用结构网格是最方便的了，直接在该区域定义一个block，单独进行网格划分

@hongjiewang 请问你是在这步Clapeyron方程用局部压力计算饱和温度吗？还有这里的饱和压力pSat是个变量吗？

我有的时候画了一个非常复杂的网格，用各种block切分，然后导入报错。然后就是重新画，每个步骤导入一次看看有没有错误。然后搞得特别细。就没错误了。目前也不知道哪里有问题。

@cfdngu 楼主您好，我现在也有你之前的问题，算了一个和你类似的例子，波浪作用在矩形防波堤下三个自由度的响应，在沿波浪传播方向的那个自由度，无论是加密重叠区域网格，还是改进时间步和库郎数的关系，都没有好转，有了解贴体网格，但还没有采用做验证。以及您说的测点如何设置能在重叠网格上随物体一起运动。期待您的回答，祝您一切都好。

请教各位使用过，浸入边界法IBM，解决流固耦合fsi问题的大佬。

最近在用foam-extend中的浸入边界法，算圆柱绕流（刚性问题） 如果我想考虑圆柱的流固耦合fsi，请问设置上的区别在哪呢？
即：浸入边界法，做流固耦合问题，和刚性问题，设置上区别在哪呢？
万分感谢！

7ffe581e-acb0-42ea-abb7-64507e4652c1-image.png

@oitocfd
您好，请问foam-extend，浸入边界法，做流固耦合问题，和刚性问题，
设置上区别在哪呢?

@micro 自挖自填。
大家可以参考东岳老师关于icoFoam解析的文章，说的很正确的，一个时间步骤动量方程的系数是不变的，因为动量方程每个时间步中仅被组建了一次。但是每个时间步中，U可是被更新了好几次（在压力泊松循环中），两者并不矛盾！另外，我们同样可知，一个时间步中，我们反复求解压力泊松方程（一般是3次），压力泊松方程的系数同样是不变的，变化的是右端项（因为右端项是用速度组建的）。

在设置探针的时候，有个fixedLocation的条目可以设置为false，这样的话对于动网格探针应该是会随着网格拓扑更新移动的。

@xpqiu 谢谢您的回复！关于2，我设置的phi的tolerance是10^-6，solver是GAMG。把nNonOrthogonalCorrectors设为10之后，最后几次计算的No. Iteration 就随计算次数减小了，最后为0，我想应该是残差已经降到最低。谢谢！

@cresendo 那就是纯对流方程了

@xpqiu 收到，还是非常感谢！！解释的已经很清晰了！

@五好青年 tutorials里面的wingMotion

目前也不确定是openMPI的问题还是OpenFOAM太挑剔，所以各位使用集群计算的老师同学们，能不能看看你们集群的hostname是不是有下划线的，能让我定位一下问题的位置。:xiezuoye:

@李东岳 了解了，谢谢东岳老师

关于OpenFOAM中CIP方法（Constrained Interpolation Profile）紧致差分插值曲线方法

请问OpenFOAM的哪个版本中，有CIP方法的代码，想学习下。

目前参考了一些文献（胡长洪、刘成、赵西增等老师），但是不知道怎么植入。请问谁做过这方面的植入（或者有现成的开源代码）么

感谢同在研究CIP的朋友回复，感谢！

实在没空弄了 我看看有没有人给你弄 https://www.cfd-china.com/topic/4861/领取任务免费赠送cfd中文网注册资格

@wying 要非常高的精度，比如普方法、有限差分那种7 8 9阶，有限体积法大体都不行

如题，根据probes的字母意思解释，是Set of locations to sample（设置位置集），我看了下probes.H和probes.C的源程序，没有看懂，probes.C的部分代码如下：

const label celli = mesh.findCell(location); elementList_[probei] = celli; if (celli != -1) { const labelList& cellFaces = mesh.cells()[celli]; const vector& cellCentre = mesh.cellCentres()[celli]; scalar minDistance = great; label minFaceID = -1; forAll(cellFaces, i) { label facei = cellFaces[i]; vector dist = mesh.faceCentres()[facei] - cellCentre;//是指面心和体心的插值么？ if (mag(dist) < minDistance) { minDistance = mag(dist); minFaceID = facei; } } faceList_[probei] = minFaceID; }

@李东岳 是这样啊，我还以为是软件系统文件缺少呢，感谢感谢

在OpenFOAM文件家下打开终端输入

find src applications -name "*.L" -type f | xargs sed -i -e 's=\(YY\_FLEX\_SUBMINOR\_VERSION\)=YY_FLEX_MINOR_VERSION < 6 \&\& \1='

@cfdngu 好的，非常感谢您快速的解答，我再尝试尝试。

@samuel-tu 主要这文章看看也很有帮助，它的角度偏代码整体架构而不是力学，对理解openfoam的设计还挺好的。

@小狗狗 谢谢您，我学习一下

@Yongbo :high: 谢谢

@bestucan 后来我也是这么解决的。。相当于直接安了新的Paraview

可以参考邱小平大佬的博客
链接文本

@starDust boundaryFieldRef()是非const的可以赋值

@bestucan 找原作者拿到了 ` 谢谢回复

@五好青年 emmm……这就不知道了，我没研究过。这只是一个线索。现成的应该是没有的。但是照着pdf里的线索查比从头做要快多了。

@Micro 关注过，但是没研究过。好像是8刚支持外挂矩阵求解器。还没机会学习:shangxue:

@shiyu 一个有点冗余的办法，定义一个vector场，每个网格都等于你那个gradP。这样就可以在别处lookup到了。

@李东岳 李老师，请问发散后继续计算又不发散了，比如如下报错：
a0c00201-ea8c-4d8d-bdad-e4e033c1ad9e-image.png
或者报错与solver有关（gauss-seidel），这证明植入的模型还是不容易收敛是吧？

pisoFoam不可以考虑重力

interFoam可以考虑中立，你这么操作就是对的，没有别的操作了。

检查下网格边界条件吧。感觉有点问题。。

最后就想请教下，各位大佬

在进行OpenFOAM并行计算时，不管是虚拟机（还是真Ubuntu系统），CPU的核心可以100%满载吗？

比如电脑16核，设置16核。有人尝试过么？实际会乱跳么？

我这并行计算，实际只占用5%，这也太坑了叭！:136:

@Samuel-Tu: 好， 但仍谢谢Tu兄的回复，Tu兄，辛苦了

@浪迹天大 谢谢，我研究一下，这个file 我就是读了一下文件名字，然后放到IOobject里，原来在of4就完全没问题

positionsFile_(this->coeffDict().lookup("positionsFile")), positions_ ( IOobject ( positionsFile_, owner.db().time().constant(), owner.mesh(), IOobject::MUST_READ, IOobject::NO_WRITE ) ), PDFFileName_(this->coeffDict().lookup("PDFFile")), PDFsamples_ ( IOobject ( PDFFileName_, owner.db().time().constant(), owner.db(), IOobject::MUST_READ, IOobject::NO_WRITE ) ),

各位简直就是大佬中的大佬 厉害厉害 :chitang:

东岳老师，我最近在做一个类似的模型，颗粒粘度与其体积分数有关的，用的of8的multiphaseEulerFoam，现在完全没有思路，不知道要将其放在哪里:136: 望您指点一二

查到了一个类似的问题，是fluent用tecplot处理的时候也发现壁面速度不为0：
出现这样问题的原因为：fluent的计算数据是存储在每个单元的中心位置，而tecplot的云图显示是按照节点来显示的。如果是直接导入fluent的.cas文件和.dat文件的话，那么在tecplot中，节点上的值会通过单元中心点的值插值得到，从而造成壁面速度不为0的情况。
不知道OF这面有无特殊处理方式

是我没有解释清楚么？这里的genotype[i]就是输入的字符，因此minus 1.0 2.0 对应三个字符串，从i=0开始，对应minus, i++变为1， 返回第一个递归函数evalGeno对应的i=1,也就是进入else, i++ 为2，numValue为1.0，继续减号，然后到第二个递归evalGeno，此时i为2，所以按程序流程就是1-2，但是程序测试结果是2-1的结果，为什么会这样呢？请各位老师指点一下

问题已解决~因为变量不能从字典中正确读取，修改读取语句后可以正常运行

@bestucan 好的哈！谢谢版主！我仔细看看您发的网址

@bestucan 好的哈，谢谢教授回复！不过您说“但是经过相邻面的速度也反了”，速度的方向是（绝对）唯一的吧，就是说该是多少就是多少。感觉还是从高斯定理出发，法向量永远取外侧为正

@OItoCFD 在 新写的程序单核计算无问题，但无法并行计算？ 中说：

mpi并行那些需要考虑吗？

一般不需要考虑

除非你的surface场进行一些特殊的插值

第31页（页码13）最上面那段是不是你要的

https://www.jianguoyun.com/p/DYNwnmQQsLv2BRiS3fED

@Yongbo 你说的两种方式是对的，我验证均通过编译了。东岳老师应该是笔误了

@cresendo 好的好的，我去试试，谢谢！！