[jackdaniels83]

Hi,
i have modified the turbfoam solver in a way that i can calculate forces in selected patches.
the problem is that the solver aborts and gives the following message:

Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model kEpsilon
Reading the patch force dictionary
Reading force solver data

Starting time loop

Time = 0.01

Courant Number mean: 0 max: 0
DILUPBiCG: Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 2.34397e-13, No Iterations 1
AMG: Solving for p, Initial residual = 1, Final residual = 0.0483882, No Iterations 1
time step continuity errors : sum local = 0.000553499, global = -8.75052e-07, cumulative = -8.75052e-07
AMG: Solving for p, Initial residual = 0.314146, Final residual = 0.0218164, No Iterations 1
time step continuity errors : sum local = 0.000303637, global = 4.96421e-05, cumulative = 4.8767e-05

Calculating Pressure Forces....
Area of patch cock = 820.124
Pressure force on patch cock = -1.02047e+08 N
Total Pressure Force in cock moving direction = -1.02047e+08 N

Calculating Wall Shear Forces....
Shear Force on patch cock = 0 N
Total Shear Force in cock moving direction = 0N
Total Force in cock moving direction = -1.02047e+08N
#0 Foam::error::printStack(Foam:stream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 Uninterpreted: [0xffffe420]
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&)
#4 void Foam::divide<foam::fvpatchfield,>(Foam::GeometricF ield<double,>&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&)
#5 Foam::tmp<foam::geometricfield<double,> > Foam::operator/<foam::fvpatchfield,>(Foam::tmp<foam::geometricfie ld<double,> > const&, Foam::GeometricField<double,> const&)
#6 Foam::turbulenceModels::kEpsilon::correct()
#7 main
#8 __libc_start_main
#9 Foam::regIOobject::readIfModified()

Can someone tell me what i have done wrong, because the problem is the same when i use the normal tubfoam solver, too.

Daniel