CFD中文网

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

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()中添加一个小数，但是也同样出现浮点数报错。

大佬们，想请教下为什么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。

@李东岳 李老师 动图中iter是指迭代步数。目前设置的是计算速度场计算200步完成后，再继续伴随速度场，这边为了减少gif的大小，我只生成了前50步的结果。
在10-15步的时候，此时计算的Ua是正常的，但是随着迭代步增加，计算就不收敛了。一直没搞明白为什么。

这是其他密度场下计算得到的Ua场：
2da7895d-d452-4997-bef2-ce5c1f02afb9-image.png

现在的问题是，有些密度场下可以计算，有些就异常了。（密度场是通过0时刻读入的，作为一个初始值）

volScalarField xh ( IOobject ( "xh", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), mesh, scalar(voluse), zeroGradientFvPatchScalarField::typeName ); xh.read();

@李东岳 在 不使用湍流模型，而是直接求解器中实现湍流计算 中说：

有一些代码不知道有啥作用，那就留着，留着总比删了更保险。

好像是的，对结果不影响的都会跟着作者留下来了。:papa:

@李东岳 在 不使用湍流模型，而是直接求解器中实现湍流计算 中说：

这个在extend那面很常见。

extend我只是听过，但没有用过。现在我常使用的版本一个是Org 6.0的，另外一个是ESI v2312的。

@李东岳 好勒好勒，谢谢李老师回答。感觉又接触到很多没学过的知识了:xiezuoye:

想请假一个问题，我在使用OpenFOAM进行拓扑优化，其中设置了时间步为每次的迭代优化步。每次迭代中比如NS方程采用Simple求解200次（for循环），但我想在下一步求解中重新计算这个NS方程，而不是采用上一步求解的迭代信息继续计算（想要保证每次迭代优化时开始求解的initial residual都和第一次迭代优化求解的initial residual一致），这个OpenFOAM可以实现吗？