@东岳 老师好，如何在指定时间范围内对新变量求平均，我尝试在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,

@东岳 多谢李老师了，是时间步长的原因，我开始按库伦数取得 听了您的建议往小了试，没有出现发散

@东岳 感谢感谢

那位OpenFOAM经验丰富的前辈能讲讲 OpenFOAM的算例之间如何耦合？

还有很神奇的事情，就是440核的话我用20个节点，每个节点22核测是可以正常跑，但是换成22个节点每个节点20核跑也出现上述问题，不知道具体问题是出在内存不够呢还是说编译器并行方式不同呢？一个节点满核28核，上面的两个我都是节点上满核跑的，也就是504核用了18个节点，644核用了23个节点。

我没仔细看，你的代码思路上没啥问题，addparticle传入一个particletype类型的指针，就是你前面用new创建的指针，但是好像参数不太对吧，new sprayparcel（一些参数），应该不是传入this指针，你可以参考一下injection model里面入射粒子部分的写法

@星星星星晴 傻了 知道了

怎么没法写公式了 = = 是不是方法用错了

67069363-6745-40a9-8a27-6aa33e64e681-image.png
质量守恒很简单，但是根据code，child的速度和数量就没办法设置。。。

我再看看哪有方法解决吧

哈啊，喝 :xiabanle: