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

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:

@李东岳 好嘞，谢谢李老师

@User0221 df37a77c-beac-4b74-ab72-89cd53ef4d8f-1728959891654.png 我在计算到这一步的时候直接停住不继续计算了，也不会报错退出，cpu占用也是满的，请问你也是这种问题吗

@学流体的小明 我仔细看看感谢老师分享

@Sk 我的问题是与phiIB计算时的局部参考点有关。由于默认设置为0号网格，对每个核来说，都会在角点位置出现一个异常的结果。我采用的方式是对phiIB这个量的计算避免使用非迪利克雷边界，进而避免引入参考点场值。

感觉应该考虑时间步长/库朗数/边界条件设置的问题

topoSet不能删除。只能把这些网格处理到某个形状里面。subset可以删除。