reactingFoam修改化学反应机理后遇到了问题

[yuky]

目前研一刚入门阶段，针对reactingfoam进行一些改动，用的是自己的构型，边界条件也相应做了一些调整，原本的参考案例是用CH4，我自己改成了C7H16，也把化学机理加了进来，现在遇到了这个问题不知道怎么办（报错如下所附），请求各位大神指点一下

--> FOAM FATAL ERROR: 
Sum of mass fractions is zero for species 
544
(
各种组分
）
    From function void Foam::basicSpecieMixture::correctMassFractions()
    in file mixtures/basicSpecieMixture/basicSpecieMixture.C at line 205.

FOAM exiting

[yuky]

这个问题解决了，但是出现了新的问题

No MRF models present

No fvModels present
No fvConstraints present
#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::nutkWallFunctionFvPatchScalarField::nut() const at ??:?
#4  Foam::nutWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
#5  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
#6  Foam::RASModels::kEpsilon<Foam::CompressibleMomentumTransportModel<Foam::dynamicTransportModel> >::correctNut() at ??:?
#7  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
#8  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
Floating point exception (core dumped)

[李东岳]

epsilon不能给0的初始场。你检查一下

[yuky]

@李东岳 好的，谢谢东岳老师，我去看看

[yuky]

又又又遇到了新问题，一步都没算下去，请问这种问题一般出在哪里
Create time

Create mesh for time = 0

PIMPLE: No convergence criteria found

PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode

Using LTS
Reading thermophysical properties

Selecting thermodynamics package
{
type heRhoThermo;
mixture multiComponentMixture;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState perfectGas;
specie specie;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Prt 0.85;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

Creating thermophysical transport model

Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity
Creating reaction model

Selecting combustion model EDC
Selecting chemistry solver
{
solver ode;
method standard;
}

standardChemistryModel: Number of species = 544 and reactions = 2446
Selecting ODE solver seulex
Creating field dpdt

Creating field kinetic energy K

No MRF models present

No fvModels present
No fvConstraints present

Starting time loop

Time scales min/max:
Flow = 9.728569846e-08, 0.0001
Temperature = 4.494232837e+307, 4.494232837e+307
Composition = 4.494232837e+307, 4.494232837e+307
Overall = 9.728569846e-08, 0.0001
Time = 1e-07

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 0.0001092756986, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.006862564884, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 1, Final residual = 0.0001167898495, No Iterations 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 double Foam::standardChemistryModel<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGasFoam::specie >, Foam::sensibleEnthalpy> > >::solve<Foam::Field<double> >(Foam::Field<double> const&) at ??:?
#4 Foam::combustionModels::EDC::correct() at ??:?
#5 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
#6 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
Floating point exception (core dumped)