虚拟机上openfoam10和11没有问题

@李东岳 在 非定常CFL数相关问题 中说：

就是简单的一个传输的标量，附加一个s的传输方程

哦哦，谢谢老师！！

@bestucan 谢谢大佬

只有一个图。那就修网格吧。

最近编写了一个拉格朗日DPMFoam下面的颗粒长大的求解器，在OpenFOAM/OpenFOAM-3.0.0/src/lagrangian/intermediate/lnInclude对源文件进行了修改。
源文件，求解器均编译成功，使用新求解器时候不能运行stochasticDispersionRAS的随机游走模型，有如下报错

Starting time loop Courant Number mean: 0 max: 0 Time = 0.01 Evolving kinematicCloud Solving 3-D cloud kinematicCloud #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigSegv::sigHandler(int) at ??:? #2 ? in "/lib64/libc.so.6" #3 std::string::rfind(char, unsigned long) const in "/lib64/libstdc++.so.6" #4 Foam::IOobject::group() const at ??:? #5 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::kModel() const at ??:? #6 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::cacheFields(bool) at ??:? #7 ? at ??:? #8 ? at ??:? #9 ? at ??:? #10 __libc_start_main in "/lib64/libc.so.6" #11 ? at ??:? Segmentation fault (core dumped）

修改了一下算法，算简单一点也有报错

Courant Number mean: 0 max: 0 Time = 0.001 Evolving kinematicCloud Solving 3-D cloud kinematicCloud [2] #[3] #0 0 Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)[1] #0 Foam::error::printStack(Foam::Ostream&)[0] #0 Foam::error::printStack(Foam::Ostream&) at ??:? at ??:? at ??:? at ??:? [3] #1 [1] #[2] #1 Foam::sigSegv::sigHandler(int)1 Foam::sigSegv::sigHandler(int)Foam::sigSegv::sigHandler(int)[0] #1 Foam::sigSegv::sigHandler(int) at ??:? [3] #2 ? at ??:? [2] #2 at ??:? [1] #2 ?? at ??:? [0] #2 ? in "/lib64/libc.so.6" [1] #3 std::string::rfind(char, unsigned long) const in "/lib64/libc.so.6" [3] #3 std::string::rfind(char, unsigned long) const in "/lib64/libc.so.6" [2] #3 std::string::rfind(char, unsigned long) const in "/lib64/libc.so.6" [0] #3 std::string::rfind(char, unsigned long) const in "/lib64/libstdc++.so.6" in "/lib64/libstdc++.so.6" [1] #4 Foam::IOobject::group() const[3] #4 Foam::IOobject::group() const in "/lib64/libstdc++.so.6" [2] #4 Foam::IOobject::group() const in "/lib64/libstdc++.so.6" [0] #4 Foam::IOobject::group() const at ??:? [1] #5 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::kModel() const at ??:? [3] #5 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::kModel() const at ??:? [2] #5 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::kModel() const at ??:? [0] #5 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::kModel() const at ??:? [1] #6 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::cacheFields(bool) at ??:? [3] #6 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::cacheFields(bool) at ??:? [2] #6 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::cacheFields(bool) at ??:? [0] #6 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::cacheFields(bool) at ??:? [1] #7 at ??:? [2] #7 at ??:? [3] #7 at ??:? [0] #7 ???? at ??:? [1] #8 at ??:? [2] #8 at ??:? [3] #8 at ??:? [0] #8 ???? at ??:? [1] #9 at ??:? [2] #9 at ??:? [3] #9 at ??:? [0] #9 ???? at ??:? [1] #10 __libc_start_main at ??:? [3] #10 __libc_start_main at ??:? [2] #10 __libc_start_main at ??:? [0] #10 __libc_start_main in "/lib64/libc.so.6" [1] #11 in "/lib64/libc.so.6" [3] #11 in "/lib64/libc.so.6" [2] #11 in "/lib64/libc.so.6" [0] #11 ???? at ??:? at ??:? at ??:? at ??:?

看起来是和Collisionparcel文件有关系？？？？？为什么我改一下颗粒直接变化还和碰撞模型的文件产生了影响？？我修改完了之后，不用随机游走模型倒是能跑，但是之后我是算例是需要运动随机游走模型的，所以我想搞想清楚原因，是不是我还要对CollisonParcel文件进行修改

@coolhhh 谢谢老师！

下面是我需要计算的方程：

fvScalarMatrix TbEqn ( - fvm::div(-phi, Tb) + fvm::laplacian(DT, Tb) ); TbEqn -= Q; TbEqn.solve();

其中，DT为热传导系数，Q为数值1的标量场，phi通过SIMPLE求解。
设置出口边界条件为：

outlet { type groovyBC; variables "Tk=DT;h=U&normal();Tinf=0;f=1/(1+Tk/(h*mag(delta())));"; valueExpression "Tinf"; gradientExpression "0"; fractionExpression "f"; value uniform 0; }

设置的divSchemes为：

div(-phi,Tb) Gauss linearUpwind grad(Tb);

使用的solvers为：

"(Tb)" { solver PBiCGStab; preconditioner DILU; tolerance 1e-6; relTol 0; }

上述方程是用于拓扑优化求解的，但是优化过程中，经常出现浮点数错误“Signal: Floating point exception (8)”，如何有效解决这个问题呢？
尝试过在h*mag(delta()中添加一个小数，但是也同样出现浮点数报错。

@jlx 请问您在fe40中安装好wave2foam了没

@youhaoyu 我是使用compressibleVoF中的throttle算例更改了模型和边界条件

@sxz0823 您好，可以加个联系方式请教您一下turbinefoam吗，我想向您请教一些问题，我的qq是：935437658 :140:

大佬们，想请教下为什么OpenFOAM自定义并行计算时，初始化（这边初始化是指下述的代码）有时候很快2-3s，有时候需要很久1-2min（测试的时候没有其他程序运行）。
下述是我自定义的并行计算初始化代码：

scalarField MyField(5*N,0); List<scalarField> AllMyField_List(Pstream::nProcs()); AllMyField_List[Pstream::myProcNo()] = MyField; Pstream::gatherList(AllMyField_List); Pstream::scatterList(AllMyField_List); scalarField AllMyField(ListListOps::combine<scalarField>(AllMyField_List, accessOp<scalarField>()));

其中，N为所有网格单元的数量，N=86400。程序是有个变量，我将其分配到各个线程进行计算，在总线程中进行收集。
此外，程序中不止这一个变量，有五六个需要同样操作的变量。之前测试过，如果N很小的时候，初始化过程很快，但是当N增大时，初始化时间就完全不一样了。想请教下如何解决这个bug。

@chenboyao 天哪，膜拜大佬，成功解决了TAT

换课题对于学生来说是非常不友好的。有时候被换了也没招..

接触角不是在三相(固/液/气)的接触点才有吗? 两相流体的界面处应该没有接触角吧?

楼上说的对。版本差的太大，移植相对麻烦。

@SSSSK 这两个图横纵坐标分别是什么？看起来确实是从20s开始统计的话收敛很快，150s开始统计收敛就好慢了，而且150-250峰值也没到2.5以上。有参考的DNS数据吗？和DNS比较一下。
你是用的哪一种驱动槽道的方式？是meanVelocityForce的话，就把pressureGradient随时间变化的图像贴一下。

@bestucan :146: :146: :146:

@李东岳 已经没招了，一直溢出:136:

在改变并行计算的数量之后，问题可解决

@foolfishzju1990 :146: :146: :146: