请问：固体颗粒的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:

你好，我也想在rhocenteralFoam中添加燃烧，做超音速和跨音速燃烧，不止能够请假该如何将燃烧模块添加到rhocenteralFoam中？期待您的答复。

@Samuel-Tu 自己写一个也不难，OF的数学库都是针对GeometricField的，估计是没有

@Samuel-Tu 哦哦 有点启发了 多谢了

大家好，我最近在使用interfoam求解器通过openmpi算一个算例，但是每次算一会儿就会出现以下错误提示：
e86f3c41-75e6-4103-a136-858021b7b253-image.png
有人能解答一下么，万分感谢。

在v1812里，使用了 cyclic boundary

preservePatches (LEFT RIGHT) (TOP DOWN);

decompose时报错：

--> FOAM Warning : Using 'preservePatches' constraint specification. This constraint keyword is deemed to be 54 months old. --> FOAM FATAL IO ERROR: Entry 'preservePatches' has 4 excess tokens in stream 8(( LEFT RIGHT ) ( TOP DOWN ))

怎么回事:chouchou:

请问先用粗网格跑至充分发展湍流再用mapfield映射到细网格上，那粗网格可以有多粗？
我目前DNS网格是7500W，那我可否就用五六百万的网格先跑至充分发展？
另外既然是用于DNS的初始场，那粗网格也是设成laminar跑是吗？谢谢。

@Samuel-Tu 谢谢

Kornev N, Shchukin E, Taranov E, et al. Development and implementation of inflow generator
for LES and DNS applications in OpenFOAM. Proceedings of the Open Source CFD International
Conference; 2009 Nov 12–13; Barcelona, Spain; 2009

@人来人往 relaxationFactors实在PIMPLE下面的，稍稍调整一下，另外删除一些没必要的东西，GAMG不需要preconditioner。
具体上网找找每个求解器的定义吧，都需要什么，什么是default的，到底哪个参数需要修改

GAMG generalised geometric-algebraic multi-grid. { solver GAMG; tolerance 1e-07; relTol 0; cacheAgglomeration true; switch specifying caching of the agglomeration strategy (default true). nCellsInCoarsestLevel 2; approximate mesh size at the most coarse level in terms of the number of cells (default 10). directSolveCoarset false use a direct solver at the coarsest level (default false). agglomerator faceAreaPair; mergeLevels 2; keyword controls the speed at which coarsening or refinement is performed. simple mesh use 2 preSweepsLevelMultiplier multiplier for the the number of sweeps between each coarsening level (default 1). maxPreSweeps maximum number of sweeps as the algorithm is coarsening (default 4). nPostSweeps number of sweeps as the algorithm is refining (default 2). postSweepsLevelMultiplier multiplier for the the number of sweeps between each refinement level (default 1). maxPostSweeps maximum number of sweeps as the algorithm is refining (default 4). nFinestSweeps number of sweeps at finest level (default 2). smoother DILUGaussSeidel; nPreSweeps 0; number of sweeps as the algorithm is coarsening (default 0). nPostSweeps 2; minIter 1; }

@yfclark 在 undefined symbol: _ZN4Foam23ControllerStateProvider8positionE错误求助 中说：

ic或sta

我是用的静态数据变量static 然后也是public成员