各位老师好！
报错信息：
Screenshot from 2021-11-20 14-31-59.png

solidParticleCloud.C的代码：

solidParticle::trackingData td(*this, rhoInterp, UInterp, nuInterp, g.value()); Cloud<solidParticle>::move(*this, td, mesh_.time().deltaTValue()); if(mesh_.time().value()> td.cloud().tInjStart_ && // mesh_.time().value()< td.cloud().tInjEnd_) {this->inject(td);} void Foam::solidParticleCloud::inject(solidParticle::trackingData &td) //. { //. label cellI = -1; //. label tetFaceI = -1; label tetPtI = -1; mesh_.findCellFacePt(td.cloud().posB_, cellI, tetFaceI, tetPtI); if(cellI > 0){ solidParticle* ptr1 = new solidParticle(mesh_, td.cloud().posB_, cellI, tetFaceI, tetPtI,td.cloud().dB_, td.cloud().UB_); Cloud<solidParticle>::addParticle(ptr1); } }

@bestucan 嗯嗯，谢谢老师，我试试:146:

@bestucan 你们都溜的很啊！

@是嘿嘿牛啊 按照我的理解，这样分的size不应该都是10吗？为什么size会变来变去呢？求大佬指点一下，拜托啦！

好像我问的有问题，$A=0.5(\overline{p_1}-\overline{p_0})\times(\overline{p_2}-\overline{p_0})$这个东西看起来是个矢量。但是面积不是矢量啊。我回去对对OpenFOAM的代码。

我看了下，OpenFOAM里面mesh.Sf()对应的是这个：$0.5(\overline{p_1}-\overline{p_0})\times(\overline{p_2}-\overline{p_0})$，mesh.magSf()对应的是这个：$\left|0.5(\overline{p_1}-\overline{p_0})\times(\overline{p_2}-\overline{p_0})\right|$，后面这个应该就是面积了。

@evensun $\left|0.5(\overline{p_1}-\overline{p_0})\times(\overline{p_2}-\overline{p_0})\right|$ 这个面积可以理解不

大家好，我想问下在rhoCentralFoam中想植入一种新的重构格式，文献中以标量守恒律作为控制方程，对其的未知量进行了界面左右值的重构，这个应该在rhoCentralFoam中哪一块进行植入呢，我不太清楚，望解答一下 :yes: ，谢谢！

@冠竹 有可能是入口速度太小了

@tens 好的！非常感谢您

@李东岳 可以把日常的流体都看成稀薄流体呀。从各种benchmark来看的话，精度还不错，处于同一个水平。
现在主要的问题是，显式LBM并行效率特别的高，然后我就寻思着有限体积也可以显式呀，谁怕谁。结果性能卡在了压力泊松方程的求解上:136:

@tens 好的，谢谢老师

已解决，postProcess -func "components“ 应该是双引号:mihu:

https://www.cfd-online.com/Forums/openfoam-solving/120990-small-great-rootvsmall-what.html
里面写了small常数还有这些以及大小

static const doubleScalar doubleScalarGREAT = 1.0e+15; static const doubleScalar doubleScalarVGREAT = 1.0e+300; static const doubleScalar doubleScalarROOTVGREAT = 1.0e+150; static const doubleScalar doubleScalarSMALL = 1.0e-15; static const doubleScalar doubleScalarVSMALL = 1.0e-300; static const doubleScalar doubleScalarROOTVSMALL = 1.0e-150;

以及在OpenFOAM-7数值改变了(https://github.com/OpenFOAM/OpenFOAM-7/commit/c902c7a396c43ff733b33294f79acab03d7b92d8#diff-118c1e97b96036fb1f9d7f1ee5912713)
d8bec1f9-3829-4c38-873f-62207d498415-image.png
测试了一下OpenFOAM-7的VSMALL结果：
430638be-f3c0-4d6c-a853-fd82f9b03400-image.png

@carmelosun 如果流体是理想气体之类，热物性是根据温度和压力计算得到，入口温度低，初始场温度高。其实入口速度还是0.4，但是由于入射的流体密度比初始场内密度大，就会有个加速效果，看起来像是入口速度大

@evensun 是boundary的patch上的值，不是边界第一层网格的值

parcel。parcel守恒，然后改变particle数

@李子橙 明白了 谢谢

@李东岳 现在的瓶颈是单个节点内存不足的问题，但优势是节点数够多。我这边一直想通过parallel mesh generation的方法生成上亿网格。步骤跟之前的帖子您提到的差不多，blockMesh->decomposePar->snappyHexMesh...原来以为这种方法，能绕开内存不足的问题，结果发现processor0通讯开销太大导致内存不足。单节点64G内存的话，千万网格大概占20%内存。因此，上亿网格需要100G以上的节点支持。但这个设计是非常浪费的，因为千万网格->亿网格加密过程中，其他processor只占内存0.4%，processor0直接超过100%。不知道东岳老大有没有其他的方法绕开内存不足的问题，多次refine的话，我认为最终refine那一次所占的内存，跟直接blockMesh所需要的内存差别不大

@dzw05 我是这么使用这个通量的，surfaceScalarField phiVp = fvc::flux(Vp); 计算时是fvm::div(phiVp, C)，然后使用求解器时要求我设置div(interpolate(Vp),C)的格式。如果使用surfaceScalarField phiVp = fvc::interpolate(Vp) & mesh.Sf();要求我设置div((interpolate(Vp)&S),C)的格式，所以我想问他们的区别。
谢谢老师

我以前也遇到过，我觉得可能是之前setfields错误的问题，后来我好像是重新删一下相文件，就没了。

@氕氘-川
时间可以接受，支持并行转换。