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

VamosZhi

@李东岳嗯嗯，那我现在能改吗？公式带‘的项都是波动项，对应于我代码中加了p的变量，如u'对应于Up，p'对应于Pp,也是我要求解的东西，公式中的g相当于一个中间变量，没有什么实际意义，对应于代码中F。
我通过计算出的那几步调用paraview发现，波动项在几步计算后便达到了e23次量级，应该是这个原因出现“浮点数除外”
我把带有F的那几项换成FVM隐式后，求解器就编译不成功了。

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，请问是否方便交流下

CFD中文网

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