openfoam11中多组分求解器混合规则调用失败
-
各位老师好,
我最近在做一个气固耦合模拟,在气相中,使用了多组分模拟,下面是我在多组分混合时选用的数据模型
type heRhoThermo;
mixture multicomponentMixture;
transport polynomial;
thermo hPolynomial;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;最初测试,transport使用的是sutherland,而 thermo使用的是janaf模型,是没有问题的,但因为缺少数据,所以,我改为了 polynomial 和 hPolynomial,程序可以跑1步,然后开始报错
报错信息如下
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 Foam::GAMGSolver::Vcycle(Foam::PtrListFoam::lduMatrix::smoother const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8 Foam::fvMatrix<double>::solve() at ??:?
#9 Foam::solvers::isothermalFluid::correctBuoyantPressure() at ??:?
#10 Foam::solvers::isothermalFluid::pressureCorrector() at ??:?
#11 ? in "/opt/openfoam11/platforms/linux64GccDPInt32Opt/bin/foamMultiRun"
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? in "/opt/openfoam11/platforms/linux64GccDPInt32Opt/bin/foamMultiRun"
./Allrun:行 22: 157260 浮点数例外 (核心已转储) foamMultiRun请问应该如何处理?谢谢