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成员

@Samuel-Tu 嗯嗯

@linhan-ge 请问您考虑了什么作用，现在颗粒润湿是不是还没人实现

想要用OpenFoam计算螺旋桨或者风扇的声传播问题。所以想问一下各位在OpenFoam中有没有求解FW——H方程的求解器？

网上的有些命令可能不对，我照着安装了，没法运行

pyFoamPlotWatch.py cavity.log，而正确的命令应该为：

pyFoamPlotWatcher.py cavity.log

@wwzhao
真的非常感谢。
抱歉， 因offline没及时发现您的回复。

@东岳 老师好，如何在指定时间范围内对新变量求平均，我尝试在controldict添加fieldaverage，但是计算时没有Mean和prime2mean文件输出，应该如何解决？

TIM图片20191111100118.png

@东岳 对，我后来是用List<type>解决的。最后用List来初始化UList

List<label> IBMarks(50,(label)1); labelUList affectedList(IBMarks); Info << affectedList<<endl; Info << affectedList[25]<<endl;

结果如下

50{1} 1

v_感觉是类似数组头的指针，运算符[]可能被重载了，可以直接取值，就像我上面的affectedList[25]，输出为1。
另外我原来

大家好，我最近在用openfoam中的interfoam进行并行计算时出现了以下的错误。

068f8746-20e6-4aab-8566-00ba4f1862d0-image.png

我看相关论坛之后说是初值设置有问题或者求解器有问题，但修改之后还是出错误，希望大家能够解答一下，万分感谢。

看同济硕士开题有人开通量重构算法的课题。。这个大家有了解吗。说是高阶的算法，实现简单。想了解一下。另外OpenFOAM可能实现吗。。

不好意思大家，图片未能上传成功，浪费大家时间了，非常抱歉。

@yfclark 方便私聊么？我也想咨询一下你关于平面 测parcel的直径问题，
我邮箱在资料

@win 如果涉及到parcel到边界的话，注意cell size 和parcel size的关系~

@ZYj

interpolationSchemes//添加离散格式 { rho.air cell; U.air cellPoint; mu.air cell; DUcDt cellPoint; } pressureGradient//添加虚拟质量力 { U U.air; } virtualMass { Cvm 0.5; }

@yfclark 真是太不好意思，又来打扰您，今天我又尝试了bluecfd的方法，按照他给的步骤，wmake完成后把库也链接到算例下的controlDict了，按照他的描述应该在每个时间步下/lagrangian/kinematicCloud`内出现一个position.orig但是无论我怎么尝试，也都出现不了，终端显示也是读取了这个命令了

writeCloudOldStyle writeCloudOldStyle1: applying to clouds: kinematicCloud

这中间哪里出了问题呢？？？

@东岳 解决了，确实是内存的问题，把并行编译核数减少就可以了

求解器编译成功了，，，
将Options中动态链接库中的：-lincompressibleRASModel 去掉就编译通过了。

这个看不出什么问题啊

主要是D的先进性起了决定性作用

你好，我也遇到了你一样斗问题，请问你最后是怎么设置成功的呢，是给定压力进出口边界条件，压力固定值进口，固定值出口，速度全部零法向梯度吗

请问入口处有压力阀控制气压，已知压力值，入口边界条件该怎么设置?压力边界可以用fixValue,速度边界用哪种啊?
出口跟外界连通，出口压力和速度边界又该如何设置呢?

用icem合并生成网格，之后在openfoam里定义ami边界条件（水槽中间有个封闭的矩形边界，里面是动网格）后，用waveDyMfoam求解器求解时，原本设置好的探针波高检测器（cellpointface）无法产生实际效果，运行后没有生成对应的surfaceElevation文件

前辈，您好，在这个设置区域的浓度方面这个fvoptions与setfields的用法上，fvoptions能实现随时间的变化，其他还要什么呢。

麻烦各位大佬看一下，小白问的可能比较蠢的问题:mihu:
Q1:最近在用DPMFoam做颗粒流，想在paraview中显示出粒子的轨迹有没有相关的例子可以参考下；
Q2:在计算管道流动过程中，流体还未流出，为什么在outLet处就可以看到相应的速度，这个地方的速度不应该是０吗

@东岳 按上诉方法已经成功删除了，还成功重装了of1806。删除主要是想腾出空间。

@秋山结月
1.对的，最好另外再建立新的几何体，用这个新的几何体包裹要refine的局部区域。
2.这个问题倒是没有遇到过，是不是stl文件的问题？
3.snap结束后不应该直接没掉啊，即便是没有需要施加边界层的区域，也会有“No layers to generate ...”，这样的信息。看看你后面有没有这样的提示信息？

感谢楼主分享

@星星星星晴
谢谢你，我已经把这个错误解决了，是因为我在files文档的第一行kOmega231/kOmega231.C应该写成kOmega231.C

new一次，delete一次，完美

cyclic和cyclicAMI的value的值没有意义，这是类的成员，必须赋值不管有没有用，随便飞翔吧老铁！

现在明白了，谢谢东岳老师讲解

@东岳 李老师，关掉考虑自转之后出现倒流现象更严重了

@星星星星晴 好的谢谢，我目前网格有7500W，质量0.72以上，结构是简单的圆管，在算DNS，ncorrector设置的2，您觉得可以吗？

大家好，我参考openfoam自带的求解器，在求解器中添加自适应网格的功能。我的算例网格总数是50000，运行算例的时候，遇到下面的错误
1.png
以下是reactMinMax.H 17-30行

17 const label tautableDim = 3; 18 19 forAll(p,i) 20 { 21 List<scalar> lookupCoordinates(tautableDim,0.0); 22 lookupCoordinates[0]=fH[i]; 23 lookupCoordinates[1]=pHI[i]; 24 lookupCoordinates[2]=THI[i]; 25 26 const List<scalar> &LOG_t_ign_temp = tignTable.lookUpBounded(lookupCoordinates,false); 27 tIgnHI[i] = min(max(1e-14,pow(10,LOG_t_ign_temp[0])),1e9); 28 29 } 30 tIgnHI.correctBoundaryConditions();

请教大家问题在哪呢

请问PBiCGStab在更稳定，收敛速度更快的情况下，会不会牺牲精度呢？还是两者都会给出正确结果？我目前正在算DNS，对于计算准确性比较关心。谢谢。

大家好，问一个很弱智的问题。。。用了openfoam的sampling surface，输出格式为raw。请问用什么软件处理呢？

此贴基本完结，按着源代码琢磨了一下午，基本弄清楚了是如何计算的。nearestPoint()主要使用了实时碰撞检测算法，可以计算出网格中心点到一个面的最短距离。这个最短距离是通过寻找面上离网格中心点的最近点来计算的。等我过段时间总结好了，到时候再补充到这个贴子下面。

@东岳 谢谢李老师！我去试试

@anubis 哦那我记错了，talphaPhi1Un是高阶的，alphaPhi10是低阶的。

在头文件找到了描述：The low-Re correction is activated by setting the entry lowReCorrection to 'on'; in this mode the model switches between laminar and turbulent functions based on the laminar-to-turbulent y+ value derived from kappa and E. When the lowReCorrection is inactive, the wall function operates in high-Re mode.可能是在constant文件里使用这epsilonWallFunction的时候能够选择打不打开LowRe

之前用interFoam，VOF做液滴在气体中的碰撞（湍流LES模型），现在想做水滴在油里的碰撞，除了要改油的物性参数（transportproPerties）以外还需要改哪些文件？
另外，我想对比层流中液滴碰撞和湍流中的液滴碰撞，但是没有看出区别是什么原因

https://github.com/OpenFOAM/OpenFOAM-2.3.x

有区别，最好下载这个官方版

if (Pstream::parRun()) { if (Pstream::master()) { Info << "balabala..." << endl; } } else { Info << "balabala..." << endl; }

这样串行并行就都OK了，可能需要#include "Pstream.H"。希望能帮到你。

@麦迪文 那您觉得应该从哪方面入手，小白我还没有改过边界

equation.png
大气边界层风力机的控制方程为上图（比较前沿的文章基本都采用上式）。我觉得OpenFOAM中的boussinesq浮力项与文章的并不一致，因此我并不敢直接用，有知道其中原理的朋友可以解释一下吗？假如要自己写一个浮力项，<potential temperature>这个带了“<>”的势温按照文献中的意思为"horizontal average"，在求解器应如何实现？是同一水平面的势温取平均还是邻近网格取平均？

请问大家 我用如下的并行命令
mpirun -np 20 mapFields sourceDirectory -parallelTarget
但是在执行该命令时出现下面的情况
Screenshot from 2019-10-22 12:49:34.png
但是map过程没有结束, 请问这种情况如何解决?还是说并行命令有问题?

@Hungryandfool 在 有没有哥们做破碎波浪的，可以交流下么？我用REEF3D做，可以和OpenFOAM的哥们讨论 中说：

2106936

OK,