[vianne]

Hi to all,

it would be really really really great if someone had an idea why this case isn't working: I'm simulating a micromixer with two inlets and one outlet, whereby the mixing task contains a reaction. I attached all the files I'm using for this case. blockMesh works, but I get the following error using reactingFoam and I'm currently running out of ideas how to solve it...

Create time
Create mesh for time = 0
// using new solver syntax:
rho
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-09;
relTol 0;
}
// using new solver syntax:
Yi
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
h
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
Reading chemistry properties
Reading g
Reading thermophysicalProperties
Selecting psiChemistryModel ODEChemistryModel<gasThermoPhysics>
Selecting thermodynamics package hPsiMixtureThermo<reactingMixture<gasThermoPhysics >>
Selecting chemistryReader chemkinReader
Selecting chemistrySolver ode
Selecting ODE solver SIBS
ODEChemistryModel: Number of species = 5 and reactions = 1
Reading field U
Reading/calculating face flux field phi
Creating turbulence model.
Selecting turbulence model type laminar
Creating field DpDt
Courant Number mean: 10.0785 max: 160
Starting time loop
Courant Number mean: 0.00629904 max: 0.0999994
deltaT = 6.24996e-07
Time = 6.24996e-07
Solving chemistry
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/steffi/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/steffi/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/steffi/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/steffi/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/reactingFoam"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/steffi/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/reactingFoam"
#6 main in "/home/steffi/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/reactingFoam"
#7 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#8 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/i386/elf/start.S:122
Floating point exception