自定义Solver出现“浮点数例外，核心已转储”

[VamosZhi]

经过试算发现，应该是求解Pp的那段程序导致的发散，但是为什么我还不清楚，有什么别的写法使它稳定吗？

[李东岳]

fvScalarMatrix PpEqn
       (
         (rho/c21)*fvm::ddt(Pp)
        +fvc::div(F)
        );

这是求解公式8么？

[VamosZhi]

@李东岳 嗯嗯是的老师，求解公式8

[wwzhao]

文献中说的是compressible，为什么要用 icoFoam 呢？

[VamosZhi]

@wwzhao 您好，声场之前的流场计算是不可压缩的。还有，可能是我不太了解，我这段程序发散：

  fvScalarMatrix PpEqn
       (
         (rho/c21)*fvm::ddt(Pp)
        +fvc::div(F)
        );
Info<<"point4.1"<<endl;
       solve(PpEqn==-(rho/c21)*fvc::ddt(p));

有哪些可能呢？

[youmengtian]

gdb跑一下，这一个错误提示根本看不出来个啥

[李东岳]

嗯，最好也贴一下发散前3步后的log

[wwzhao]

@vamoszhi 我的意思是，文献的方法是用于可压缩流体，你应用到不可压缩流不会有问题么？

[VamosZhi]

@wwzhao 哦哦好的，我再想想看。不过这种方法的流场部分本来是按照不可压缩流算得，这个应该是没有问题的，我的工作是把这部分方程加到计算流场的大循环中。

[VamosZhi]

各位老哥，是我蠢了，我昨天把时间步长降到0.00001,居然算过去了。。。谢谢各位大佬帮助，我现在新问题是，如何在场内某一点上加上随时间三角函数变化的波动项。我自己查查，各位了解也可以提示我用哪种函数

[xuebao1989]

@vamoszhi 您的声波扰动方程添加成功了吗？我也准备将声学扰动方程加入OF，请问是否方便交流下