CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   porousInterFoam (http://www.cfd-online.com/Forums/openfoam/125764-porousinterfoam.html)

 hesamgh October 30, 2013 12:14

porousInterFoam

hi,

i am trying to compile the porousInterFoam solver, but whem i run the example, this error occurs ,

MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321967 Max(alpha1) = 1.62518
MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321915 Max(alpha1) = 1.62518
DICPCG: Solving for p_rgh, Initial residual = 0.924755, Final residual = 0.0281868, No Iterations 257
time step continuity errors : sum local = 2.13828e-05, global = 1.74228e-07, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.645343, Final residual = 0.030643, No Iterations 259
time step continuity errors : sum local = 0.000109586, global = 3.36113e-07, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.748674, Final residual = 9.91322e-08, No Iterations 446
time step continuity errors : sum local = 2.27742e-09, global = 3.95198e-11, cumulative = -1.00629
ExecutionTime = 3.25 s ClockTime = 3 s

Courant Number mean: 0.146919 max: 224.824
Interface Courant Number mean: 0.00106774 max: 6.39375
deltaT = 8.82449e-15
--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 1029
Increased the timePrecision from 10 to 11 to distinguish between timeNames at time 0.00459495
Time = 0.0045949473993

--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 1029
Increased the timePrecision from 11 to 12 to distinguish between timeNames at time 0.00459495
MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321899 Max(alpha1) = 1.62518
--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 1029
Increased the timePrecision from 12 to 13 to distinguish between timeNames at time 0.00459495
MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321884 Max(alpha1) = 1.62518
DICPCG: Solving for p_rgh, Initial residual = 0.997524, Final residual = 0.0482074, No Iterations 119
time step continuity errors : sum local = 1.98589e-05, global = -2.09955e-06, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.109669, Final residual = 0.00416443, No Iterations 291
time step continuity errors : sum local = 3.66002e-06, global = -1.02672e-08, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.0447736, Final residual = 9.61283e-08, No Iterations 531
time step continuity errors : sum local = 8.9133e-11, global = 3.26094e-14, cumulative = -1.00629
ExecutionTime = 3.91 s ClockTime = 4 s

Courant Number mean: 0.000801206 max: 1.18367
Interface Courant Number mean: 1.45436e-05 max: 0.066757
deltaT = 3.72759e-15
Time = 0.004594947399295

MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321872 Max(alpha1) = 1.62518
MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.32186 Max(alpha1) = 1.62518
DICPCG: Solving for p_rgh, Initial residual = 0.751844, Final residual = 0.0271748, No Iterations 206
time step continuity errors : sum local = 1.02263e-05, global = -4.9812e-08, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.138542, Final residual = 0.00650515, No Iterations 255
time step continuity errors : sum local = 4.51487e-06, global = -1.06425e-08, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.0540249, Final residual = 9.98206e-08, No Iterations 410
time step continuity errors : sum local = 7.93589e-11, global = 3.34275e-14, cumulative = -1.00629
ExecutionTime = 4.54 s ClockTime = 5 s

Courant Number mean: 0.00291016 max: 2.35784
Interface Courant Number mean: 7.37114e-05 max: 0.327241
deltaT = 7.90465e-16
Time = 0.004594947399296

MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321613 Max(alpha1) = 1.62518
--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 1029
Increased the timePrecision from 13 to 14 to distinguish between timeNames at time 0.00459495
MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321366 Max(alpha1) = 1.62518
DICPCG: Solving for p_rgh, Initial residual = 0.510742, Final residual = 0.010461, No Iterations 272
time step continuity errors : sum local = 1.17433e-05, global = -5.60236e-09, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.379362, Final residual = 0.0155277, No Iterations 196
time step continuity errors : sum local = 3.53471e-05, global = -1.73532e-06, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.270411, Final residual = 5.46214e-08, No Iterations 468
time step continuity errors : sum local = 1.62082e-10, global = 2.19395e-13, cumulative = -1.00629
ExecutionTime = 5.2 s ClockTime = 5 s

Courant Number mean: 0.00519716 max: 4.91789
Interface Courant Number mean: 0.000167649 max: 0.556441
deltaT = 8.03662e-17
Time = 0.0045949473992963

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 void Foam::MULES::limit<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double, int, bool) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/interFoam"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/interFoam"
Floating point exception (core dumped)

:):):):):)

 ngj October 30, 2013 14:22

Hallo Hesam,

Your time steps are in the order of machine precision, so something went wrong long time before. I suspect that your boundary conditions or the definition of your porous zone(s) are not correct.

Kind regards

Niels

 hesamgh November 15, 2013 04:05

porous

hi dear friend
excuse me for delay

do you know how can i identify the porous medium in openfoam when import the shape from gambit to openfoam?
i make the porous media in gambit from this address:
operation,zones,specify continuum types