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

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