CFD Online Discussion Forums - porousSimpleFoam

Welcome Foamers ;)

I spent all day on searching solutions for this error. I use porousSimpleFoam and I took mesh from Fluent. look's like a simple problem - the one input one output plus the porous zone. after 20 iterations the simulation breaks down..

I would ask for any suggestions what might be wrong where I made a mistake?

[plgsebaa@n1-4-4 wydob_1]$ porousSimpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.7.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.7.1-03e7e056c215
Exec : porousSimpleFoam
Date : Jan 04 2012
Time : 19:23:06
Host : n1-4-4.local
PID : 7370
Case : /mnt/auto/people/plgsebaa/test_case/wydob_1
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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 RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
sigmaEps 1.3;
}

Creating porous zone: POROUS_ZONE

Starting time loop

Time = 0.01

smoothSolver: Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.068853, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0364543, No Iterations 9
time step continuity errors : sum local = 3.87347e-12, global = 3.30007e-15, cumulative = 3.30007e-15
smoothSolver: Solving for epsilon, Initial residual = 0.999623, Final residual = 1.86974e-14, No Iterations 2
smoothSolver: Solving for k, Initial residual = 1, Final residual = 4.7752e-11, No Iterations 2
ExecutionTime = 0.74 s ClockTime = 1 s

Time = 0.02

smoothSolver: Solving for Ux, Initial residual = 0.342943, Final residual = 0.0107729, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.638351, Final residual = 0.0246273, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.592757, Final residual = 0.0367471, No Iterations 2
GAMG: Solving for p, Initial residual = 0.121501, Final residual = 0.00582584, No Iterations 3
time step continuity errors : sum local = 2.19959e-09, global = 3.49719e-12, cumulative = 3.50049e-12
smoothSolver: Solving for epsilon, Initial residual = 0.99988, Final residual = 1.16313e-15, No Iterations 2
smoothSolver: Solving for k, Initial residual = 0.999905, Final residual = 2.64026e-09, No Iterations 2
ExecutionTime = 1.2 s ClockTime = 1 s

...

Time = 0.21

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 9.29518e-50, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 1.21793e-49, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 3.46985e-50, No Iterations 2
GAMG: Solving for p, Initial residual = 0.00823832, Final residual = 0.000275484, No Iterations 1
time step continuity errors : sum local = 5.75209e+56, global = -7.33984e+51, cumulative = -7.33984e+51
smoothSolver: Solving for epsilon, Initial residual = 0.00156525, Final residual = 87816.5, No Iterations 1000
bounding epsilon, min: 8.27141e-24 max: 1.32026e+116 average: 1.68431e+112
smoothSolver: Solving for k, Initial residual = 3.36682e-29, Final residual = 3.36682e-29, No Iterations 0
ExecutionTime = 18.27 s ClockTime = 19 s

Time = 0.22

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 7.10764e-85, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 6.8068e-85, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 2.53065e-85, No Iterations 2
#0 Foam::error::printStack(Foam::Ostream&) in "/software/local/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/software/local/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/software/local/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/software/local/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/software/local/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/software/local/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/software/local/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#8 main in "/software/local/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/porousSimpleFoam"
#9 __libc_start_main in "/lib64/libc.so.6"
#10 _start in "/software/local/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/porousSimpleFoam"
Floating point exception

and files with msh

Code:

http://speedy.sh/pbF98/porousSimpleFoam.tar.gz

thanks for reply,
Sebaj