CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   pimpleFoam suddenly get exit (https://www.cfd-online.com/Forums/openfoam-solving/218914-pimplefoam-suddenly-get-exit.html)

Milan Shrestha July 9, 2019 07:51

pimpleFoam suddenly get exit
 
Hi everyone,

I have been trying to run the pimpleFoam into my propeller case but there has been some error on executing the pimpleFoam after some minutes of execution.

Code:

PIMPLE: Iteration 1
AMI: Creating addressing and weights between 37933 source faces and 37933 target faces
AMI: Patch source sum(weights) min/max/average = 0.802181, 2.00213, 1.00903
AMI: Patch target sum(weights) min/max/average = 0.000569305, 2.01386, 1.0088
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.00652399, No Iterations 3
time step continuity errors : sum local = 3.64691e-10, global = 1.50977e-11, cumulative = -2.62063e-09
smoothSolver:  Solving for Ux, Initial residual = 0.000880174, Final residual = 7.935e-07, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.000120996, Final residual = 1.37146e-07, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.000929456, Final residual = 6.0459e-07, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0059716, Final residual = 3.85954e-05, No Iterations 4
time step continuity errors : sum local = 5.12401e-11, global = 1.72928e-11, cumulative = -2.60334e-09
GAMG:  Solving for p, Initial residual = 0.00229053, Final residual = 4.39153e-07, No Iterations 6
time step continuity errors : sum local = 5.94828e-13, global = 1.61902e-13, cumulative = -2.60317e-09
PIMPLE: Iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.000274223, Final residual = 2.82053e-07, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 3.76589e-05, Final residual = 4.82803e-08, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.000326422, Final residual = 2.12023e-07, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.00601082, Final residual = 3.15255e-05, No Iterations 4
time step continuity errors : sum local = 4.2683e-11, global = -1.4213e-11, cumulative = -2.61739e-09
GAMG:  Solving for p, Initial residual = 0.00213069, Final residual = 5.31767e-07, No Iterations 6
time step continuity errors : sum local = 7.03291e-13, global = -1.76353e-13, cumulative = -2.61756e-09
smoothSolver:  Solving for epsilon, Initial residual = 1.84801e-05, Final residual = 6.1009e-08, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 4.95208e-05, Final residual = 1.68352e-07, No Iterations 1
ExecutionTime = 803.14 s  ClockTime = 803 s

forces forces write:
    sum of forces:
        pressure : (-2.94193 36.0195 -30.8047)
        viscous  : (-0.0142583 -0.0344043 -0.0187688)
        porous  : (0 0 0)
    sum of moments:
        pressure : (-0.341646 -29.0562 -3.71919)
        viscous  : (4.3322e-05 -0.0234912 0.00595161)
        porous  : (0 0 0)

Courant Number mean: 0.000539523 max: 0.998274
deltaT = 3.63716e-06
Time = 0.000137992

PIMPLE: Iteration 1
AMI: Creating addressing and weights between 37933 source faces and 37933 target faces
AMI: Patch source sum(weights) min/max/average = 0.796835, 2.02346, 1.00927
AMI: Patch target sum(weights) min/max/average = 0, 2.043, 1.00903
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.00642757, No Iterations 3
time step continuity errors : sum local = 3.65445e-10, global = 1.45006e-11, cumulative = -2.60306e-09
smoothSolver:  Solving for Ux, Initial residual = 0.000874066, Final residual = 7.71885e-07, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.00011986, Final residual = 1.32681e-07, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.000931321, Final residual = 5.87728e-07, 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  Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5  ? at ??:?
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  ? at ??:?
Floating point exception (core dumped)

Also, on executing the checkMesh, I got
***Max skewness = 7.68338, 60 highly skew faces

I have been trying to run it in both serial and parallel but non of this work. Same problem reoccur.


what could be the possible error to this? I don't know whether the problem arise because of skewness or something else.


I am working in a project and there is no one to help me in OpenFOAM. This is only one option where I could ask for some help.


Thank you.


Best regards,

Krao July 9, 2019 10:50

Hi Milan,

Go through this link https://www.cfd-online.com/Forums/op...-problems.html. Here they have mentioned what are the effects of having skewness and also how to visualize and overcome it. Try to always produce good and quality mesh.

Also, go through this link you will learn a lot to tackle some basic issues in OpenFOAM http://www.wolfdynamics.com/wiki/OFtipsandtricks.pdf. Hope it helps.

Krao

Milan Shrestha July 11, 2019 01:07

Quote:

Originally Posted by Krao (Post 738441)
Hi Milan,

Go through this link https://www.cfd-online.com/Forums/op...-problems.html. Here they have mentioned what are the effects of having skewness and also how to visualize and overcome it. Try to always produce good and quality mesh.

Also, go through this link you will learn a lot to tackle some basic issues in OpenFOAM http://www.wolfdynamics.com/wiki/OFtipsandtricks.pdf. Hope it helps.

Krao




Thank you for the support and reminding me to go through wolfDyanmics files. :)


All times are GMT -4. The time now is 07:40.