@东岳 东岳老师高抬了- -我做OpenFOAM才做了不到两年，你才是大佬啊

很古老的OpenFOAM里面存在pd，我已经记不清什么概念了，和p_rgh不是一个东西

好的，我认真看看，非常感谢您的回答，谢谢

谢谢东岳老师。
我的研究目的是看使用毛巾往复擦桌子过程中，由于毛巾的往复运动，引起的周围气流变化，毛巾速度是0.2m/s，如下图，先向左前进1s，再向右运动1s，回来。d49e6d5e-44a9-45c2-add6-7afe01a3bee3-image.png
这种匀速的往复匀速，速度边界条件应该如何给定？

谢谢您。

@zhe 直接设置水体区域，没有另外处理

dyfluid@dyfluid:~/twoWayHyQMOMFoam/tutorials/PTC_2D/testCase2$ oneWayHyQMOMFoam oneWayHyQMOMFoam: symbol lookup error: oneWayHyQMOMFoam: undefined symbol: _ZN4Foam6HyQMOMC1ERKNS_7PtrListINS1_INS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEEEEEE

原因：当时调整了src下的的构造函数，紧跟着solver那面的构造函数也加了一个参数进去。然后重新编译了src，忘记了编译solver，出错。
方案：重新编译solver后解决。

#0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigSegv::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"

原因：在class里面增加了一个私有成员，编译库之后，求解器可以对类进行初始化，随即报错。报错只有3行，没有任何其他信息
方案：编译class之后，尤其这种增加私有成员的时候，记得重新编译solver一下就可以解决

@东岳 ok:huahua:

嗯嗯，好的，谢谢李老师指点:chitang:

边界条件、网格、格式

不可压缩LES模型算法非常成熟了，一般都是设置不正确带来的发散

按照twoPhaseEulerFoam中的EEqns.H写出能量方程如下（算例文件中能量形式是sensibleInternalEnergy）大神们先看看这么写对吗？

fvScalarMatrix E1Eqn ( fvm::ddt(alpha1, rho1, he1) + fvm::div(alphaRhoPhi1, he1) - fvm::Sp(contErr1, he1) + fvc::ddt(alpha1, rho1, K1) + fvc::div(alphaRhoPhi1, K1) - contErr1*K1 + ( he1.name() == thermo1.phasePropertyName("e") ? fvc::div(fvc::absolute(alphaPhi1, alpha1, U1), p) + p*fvc::ddt(alpha1) : -alpha1*dpdt ) - fvm::laplacian ( fvc::interpolate(alpha1) *fvc::interpolate(thermo1.alphaEff(phase1.turbulence().mut())), he1 ) ); E1Eqn.relax(); E1Eqn -= ( Kh*(thermo2.T() - thermo1.T()) + Kh*he1/Cpv1 - fvm::Sp(Kh/Cpv1, he1) + alpha1*rho1*(U1&g) + fvOptions(alpha1, rho1, he1) );

a8753d25-8ea0-4cf2-8fae-10a71aa5cdfd-image.png
那如果不考虑传热的话，只把导热系数设为0和温度设为一样，这样做对吗？
我为了不考虑传热，是直接把twoPhaseEulerFoam.C中的EEqns.H给注释掉了，我这种做法对吗？

请问：固体颗粒的turbulent viscosity 怎么编程让其等于流体的？

请教一下大家，我现在取一小的truncated domain作为LES计算域，初始RANS结果是在20倍的弦长上计算的，截断区域是法向距离为0.5倍弦长，尾迹区域为1倍的弦长，如图所示。
truncated domain.jpg
入口设置为RANS结果插值速度和压力，通常OF入口压力都设为ZG，我不清楚这种插值速度设置是否可行，大家如果有知道的是否可以提供一些建议参考？
出口速度设置为advective，请问此时压力该如何设置？我设置为ZG, 计算结果直接从出口开始发散并向内部区域传递。

我在使用solidDisplacementFoam求解器进行热应力计算时，需要对位移向量场重新设置边界条件：其中的x,y方向为zeroGradient边界，z方向为fixedValue边界。请问OpenFOAM里有类似的边界条件或者知道该怎么设置的吗？？

李老师，我之前在of6和of4.1中的 controlDict里面设置为 边计算边输出空穴体积变化，现在2.3中我目前没有找到 边计算边监控体积变化的方法。。。

老师，您可否指点一个实现上述过程的思路，我去试试:xinxin:

貌似需要两边都懂才能讲出区别... 就像你明白你那面，我明白openfoam这面，但好像联系不起来，等待大佬分享..

建立了一个长宽高（x-y-z）为3015的立方体，目的是建立水槽。
1.jpg
使用setFieldDict字典设置水域范围，目标水域高度3米

FoamFile { version 2.0; format ascii; class dictionary; location "system"; object setFieldsDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // defaultFieldValues ( volScalarFieldValue alpha.water 0 ); regions ( boxToCell { box (-1.0 2.0 -1.0) (30.0 -2.0 3.0); fieldValues ( volScalarFieldValue alpha.water 1 ); } );

只计算到第5个时间步电脑就不运行了，有的时候可以一直运行，但是结果是全流域全部都是空气没有set的water field，而且我的建的流域模型X方向上的长度变成了15m，被砍掉了一半。
21.jpg
请老师们指教！感谢！

先对路过的大佬表示下感谢！
我看到不同算例中使用约束有constraints和restraints两种不同关键字，它们有区别吗？其中fixedPoint和fixedLine能理解，但fixedAxis代表什么呢？另外，我模拟的是浮体在水面上自由漂浮，不受约束，设置时直接将这个关键字删除了，但是运行过程发现浮体不动了，这是怎么回事呢？

5b9ac92e-ef00-485c-92bc-983ef6807c88-image.png
我想在求解浓度标量传输方程的时候将浓度的边界条件设置成图片的形式，其中速度是根据我0文件下的速度进行提取的，这个边界条件openfoam里可以直接设定吗，还是要我自己改变边界条件？

用tecplot得到圆管中切片x=15D的速度流场
small.JPG
希望能够得到如下图的速度曲线
案例.png
目前导出所有数据，摸索用matlab取一些半径上的平均速度得到的折线太丑了:xinlei: ；(https://www.cnblogs.com/code-saturne/p/10491817.html)的方法没怎么搞懂。请问用tecplot或paraview该怎样进行处理，恳求各位老师指点一下:xiexie:

Hello foamers:
如题
比如我想画水体的流场矢量图
那么应该怎样只把液相（体积分数>0.5）的速度场数据提取出来呢？
谢谢！

@范准
兄弟
这个问题解决了吗:chouchou:

@东岳 谢谢东岳老师，就是下面A beta Ta三个参数求化学反应速率常数，化学反应速率中反应物质浓度次幂直接在方程里面表现出来了哇，没有是默认为0次吗？

这个问题解决了
phaseProperties文件中的type设置错误

东岳老师好，大家好：
RT
背景网格全是正方形
然后用一条斜线来切割
用SHM画完之后总是不能贴合得很光滑
像阶梯一样（见下图）
请问怎样才能完全贴合成一条光滑的斜线呢
请指教，谢谢！
1.png

感谢李老师！:xinxin:

@一颗鸭蛋 谢谢你的回答，我用的重叠网格，贴体网格和背景网格的长宽都是1:1且一样大，然后就不知道问题出在哪里了:zoule:

我用的Codelite，查找头文件和某个函数都很方便。

补充：

根据网上说的，加入对应代码之后，编译了ThirdParty-$versionNum/Allwmake，已经成功

然后编译src/Pstream文件下面的 Allmake文件，会报错如下：

'/home/zywang2/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/dummy/libPstream.so' is up to date. + case "$WM_MPLIB" in + set +x Note: ignore spurious warnings about missing mpicxx.h headers wmake libso mpi make: mpicc: Command not found SOURCE=UOPwrite.C ; g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-100 -DOMPI_SKIP_MPICXX -IlnInclude -I. -I/home/zywang2/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude -I/home/zywang2/OpenFOAM/OpenFOAM-2.3.0/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOptSYSTEMOPENMPI/UOPwrite.o UOPwrite.C:29:17: fatal error: mpi.h: No such file or directory #include "mpi.h" ^ compilation terminated. make: *** [Make/linux64GccDPOptSYSTEMOPENMPI/UOPwrite.o] Error 1

错误是 没有找到对应的mpi.h文件，，，在以前的论坛上也没有找到有效的解决办法，，，

这个问题应该怎么解决呢？

恳请哪位热心朋友解答一下，谢谢

谢谢教授的夸奖，我的是弄的比较简单:mihu:

@东岳 是的，李老师，感谢李老师！两个版本的定义不一样！OpenFOAM-4中GeometricField字典中是这样定义的：TIM截图20191122141846.png
而foam-extend-4.0中是这样定义的：TIM截图20191122141934.png

@东岳 感谢李老师解答，我再试试

可不可以给速度场给两个入口边界类型呢

@Dingcy 哈哈哈哈，我就瞎看看，我主要做CFDEM方面

谢谢教授

@bestucan 在 openfoam5 安装问题 中说：

sudo apt install paraview

好的谢谢~

ESI版本1906，有个bug，在compressibleInterFoam下，将turbulence models换成twoPhaseTransport时，除了laminar外都不能算，直接报错。bug提交在8个月之前，但是并没有人回应，不知其他人有没有研究。

blockmesh需要在case里运行，这是非常基础的操作。找找东岳流体翻译的openfoam手册看看。

B336D52677F66CD2291AFAA97B39D47C.png
如图所示，不太了解这个错误是由于什么原因引起的，求解器是overOlaDyMFlow，恳请路过的大佬帮忙解答一下！

@东岳 谢谢您，是边界条件的问题。

为何不固定物体让水面上移？

@Jhhu 你好，我最后没有成功，由于时间很紧张，还是使用了商用软件。

。。。。。。自问自答 ，汗

phiHbyABf不是边界流量。。。。

所以以下代码应该是符合动量方程，猜测可能是进入压力PISO循环之前修改下边界条件的系数

( phiHbyABf[patchi] - rho.boundaryField()[patchi] *MRF.relative(SfBf[patchi] & UBf[patchi], patchi) ) /(magSfBf[patchi]*rhorAUBf[patchi])

找到了一个方法
用动态对象数组

IB **myIB=new IB *[3]; for(int i=0;i<3;i++) { myIB[i]=new IB(i+1); }

注意最后delete。

@nbyjn 界面重构是为了消除界面的数值耗散，vof都需要重构，只不过代数重构不需要重构准确的界面位置

谢谢两位大佬@史浩 @东岳

各位大佬，of里面除了时均值有面均值吗？就是假如有个1m长的水管，网格dx=0.001m，能不能在某个时间步下输出1000个值代表着一共1000个剖面上的某个物理量的面平均值？如果要自己编程解决，那下面这25284个网格怎么和实际坐标联系起来？谢谢

dimensions [0 0 0 0 0 0 0];//整个计算域所有网格的空气相分数 internalField nonuniform List<scalar> 25284 ( 0.878566 0.854193 0.876398 0.860644 0.879329 0.859877 0.878446 0.860776

@东岳 谢谢，改完您建议的就没有那个警告了

把原来流场的小眼睛点开

#include "fixedValueFvPatchFields.H" #include "zeroGradientFvPatchFields.H" #include "wallPolyPatch.H" if isA<wallPolyPatch> ( mom[0][0].boundaryField()[P] ) else if ( isA<fixedValueFvPatchScalarField> ( mom[0][0].boundaryField()[P] ) )

@东岳 嗯嗯 好的 谢谢啦:shangxue: