CFD Online Discussion Forums - studying a valve case

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)

- - studying a valve case (https://www.cfd-online.com/Forums/openfoam/121522-studying-valve-case.html)

Quote:

Originally Posted by mina.basta (Post 447638)

Dear Martin,

here is the test case

Could you add to the test case the "system" folder. I would like to see settings of the numerics also.

Quote:

Originally Posted by novakm (Post 447645)

Could you add to the test case the "system" folder. I would like to see settings of the numerics also.

Here is the system folder of the simple case (simulation inside cylinder)

Quote:

Originally Posted by mina.basta (Post 447652)

Here is the system folder of the simple case (simulation inside cylinder)

the boundary conditions does not correspond to the geometry.

Prepare the test case properly and post it here. After that I may take a look on it.

BR

Martin

Quote:

Originally Posted by novakm (Post 447655)

the boundary conditions does not correspond to the geometry.

Prepare the test case properly and post it here. After that I may take a look on it.

BR

Martin

Dear Martin,

Which case do you mean because at the begining i posted my valve studying case problem . After that i posted just a simple case to generalize my problem which is just a cylinder see post #18.
So i thought that you told me about the simple case which is just a cyliner and i had the same problem at the interface between the inlet and the wall with elements with very high pressure.

Quote:

Originally Posted by mina.basta (Post 447665)

Just the simplest case that shows the problem you have.

Dear Martin,

Excuse me , I attached a wrong case. Here are the write one.

Best regards,
Mina

Hello all,

Did anybody have a solution for my problem?

Best regards,

Mina

Hi Mina
I changed divScheme by adding bounded to the terms with phi and it seems run without problem.

Code:

divSchemes

{

    default                          none;

    div(phi,U)                       bounded Gauss GammaV 0.5; //was 0.5 1.0

    div(phi,k)                       bounded Gauss limitedLinear 1;

    div(phi,epsilon)                 bounded Gauss limitedLinear 1;

    //div((nuEff*dev(grad(U).T())))    Gauss linear;

    div((nuEff*dev(T(grad(U)))))       Gauss linear;

    //div((nuEff*dev(grad(U).T())))    Gauss linear faceMDLimited Gauss linear 0.85;

    //div(phi,U)                       Gauss linearUpwindV faceMDLimited linear 0.68; //was 0.5 1.0

    //div(phi,k)                       Gauss linearUpwind faceMDLimited Gauss linear 0.2;

    //div(phi,epsilon)                 Gauss linearUpwind faceMDLimited Gauss linear 0.91;

    //div((nuEff*dev(grad(U).T())))    Gauss linearUpwind faceMDLimited Gauss linear 0.85;

}

Quote:

Originally Posted by immortality (Post 448560)

Hi Mina
I changed divScheme by adding bounded to the terms with phi and it seems run without problem.

Code:

divSchemes

{

    default                          none;

    div(phi,U)                       bounded Gauss GammaV 0.5; //was 0.5 1.0

    div(phi,k)                       bounded Gauss limitedLinear 1;

    div(phi,epsilon)                 bounded Gauss limitedLinear 1;

    //div((nuEff*dev(grad(U).T())))    Gauss linear;

    div((nuEff*dev(T(grad(U)))))       Gauss linear;

    //div((nuEff*dev(grad(U).T())))    Gauss linear faceMDLimited Gauss linear 0.85;

    //div(phi,U)                       Gauss linearUpwindV faceMDLimited linear 0.68; //was 0.5 1.0

    //div(phi,k)                       Gauss linearUpwind faceMDLimited Gauss linear 0.2;

    //div(phi,epsilon)                 Gauss linearUpwind faceMDLimited Gauss linear 0.91;

    //div((nuEff*dev(grad(U).T())))    Gauss linearUpwind faceMDLimited Gauss linear 0.85;

}

Dear Ehsan,

I got this error after modifying divScheme:

Code:

--> FOAM FATAL IO ERROR: 

Unknown convection scheme bounded



Valid convection schemes are :



1

(

Gauss

)





file: /home/fayez/Bureau/ehsan_test/system/fvSchemes::divSchemes::div(phi,U) at line 30.



    From function convectionScheme<Type>::New(const fvMesh&, const surfaceScalarField&, Istream&)

    in file /home/opencfd/OpenFOAM/OpenFOAM-2.1.1/src/finiteVolume/lnInclude/convectionScheme.C at line 97.



FOAM exiting

what do you think ?

it works fine for me!

Code:

Time = 293



smoothSolver:  Solving for Ux, Initial residual = 0.00302422, Final residual = 1.63714e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00286406, Final residual = 1.40895e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00163624, Final residual = 9.37947e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00816459, Final residual = 4.26288e-05, No Iterations 2

time step continuity errors : sum local = 4.67106, global = -1.04718, cumulative = -1602.9

smoothSolver:  Solving for epsilon, Initial residual = 0.000209883, Final residual = 1.19376e-06, No Iterations 2

bounding epsilon, min: -138.041 max: 2.9235e+08 average: 9.09496e+06

smoothSolver:  Solving for k, Initial residual = 0.000215752, Final residual = 8.36107e-07, No Iterations 3

ExecutionTime = 1.22 s  ClockTime = 1 s



Time = 294



smoothSolver:  Solving for Ux, Initial residual = 0.00299998, Final residual = 1.35978e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00283594, Final residual = 1.3551e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00161325, Final residual = 8.1245e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00804053, Final residual = 1.87376e-05, No Iterations 2

time step continuity errors : sum local = 2.05042, global = -1.73826, cumulative = -1604.64

smoothSolver:  Solving for epsilon, Initial residual = 0.000211165, Final residual = 1.19876e-06, No Iterations 2

smoothSolver:  Solving for k, Initial residual = 0.000214596, Final residual = 8.11414e-07, No Iterations 3

ExecutionTime = 1.22 s  ClockTime = 1 s



Time = 295



smoothSolver:  Solving for Ux, Initial residual = 0.00269395, Final residual = 1.35094e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00268798, Final residual = 1.34244e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00138372, Final residual = 8.11865e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00734547, Final residual = 2.72228e-05, No Iterations 2

time step continuity errors : sum local = 2.98289, global = -1.95487, cumulative = -1606.59

smoothSolver:  Solving for epsilon, Initial residual = 0.000209892, Final residual = 1.16917e-06, No Iterations 2

smoothSolver:  Solving for k, Initial residual = 0.000211634, Final residual = 8.11507e-07, No Iterations 3

ExecutionTime = 1.23 s  ClockTime = 1 s



Time = 296



smoothSolver:  Solving for Ux, Initial residual = 0.00260183, Final residual = 1.34072e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00259837, Final residual = 1.34191e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00141887, Final residual = 8.04426e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00717931, Final residual = 2.6446e-05, No Iterations 2

time step continuity errors : sum local = 2.89427, global = -1.88796, cumulative = -1608.48

smoothSolver:  Solving for epsilon, Initial residual = 0.000205624, Final residual = 1.19057e-06, No Iterations 2

smoothSolver:  Solving for k, Initial residual = 0.000209961, Final residual = 7.95193e-07, No Iterations 3

ExecutionTime = 1.23 s  ClockTime = 1 s



Time = 297



smoothSolver:  Solving for Ux, Initial residual = 0.00246914, Final residual = 1.34119e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00247703, Final residual = 1.31678e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00133655, Final residual = 8.0358e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00733201, Final residual = 2.25581e-05, No Iterations 2

time step continuity errors : sum local = 2.47292, global = -1.16714, cumulative = -1609.65

smoothSolver:  Solving for epsilon, Initial residual = 0.000209202, Final residual = 1.17054e-06, No Iterations 2

bounding epsilon, min: -15.779 max: 2.92343e+08 average: 9.09302e+06

smoothSolver:  Solving for k, Initial residual = 0.000209555, Final residual = 8.24086e-07, No Iterations 3

ExecutionTime = 1.24 s  ClockTime = 1 s



Time = 298



smoothSolver:  Solving for Ux, Initial residual = 0.00262038, Final residual = 1.38132e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00253905, Final residual = 1.35259e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00141747, Final residual = 8.14806e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00717604, Final residual = 2.40071e-05, No Iterations 2

time step continuity errors : sum local = 2.62906, global = -0.296388, cumulative = -1609.95

smoothSolver:  Solving for epsilon, Initial residual = 0.00020553, Final residual = 1.17951e-06, No Iterations 2

smoothSolver:  Solving for k, Initial residual = 0.000212137, Final residual = 7.99386e-07, No Iterations 3

ExecutionTime = 1.24 s  ClockTime = 1 s



Time = 299



smoothSolver:  Solving for Ux, Initial residual = 0.00262036, Final residual = 1.34373e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00243577, Final residual = 1.31502e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.0013605, Final residual = 8.10666e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00723721, Final residual = 1.94311e-05, No Iterations 2

time step continuity errors : sum local = 2.13106, global = 0.691298, cumulative = -1609.26

smoothSolver:  Solving for epsilon, Initial residual = 0.000209559, Final residual = 1.17209e-06, No Iterations 2

bounding epsilon, min: -24.6726 max: 2.92349e+08 average: 9.09526e+06

smoothSolver:  Solving for k, Initial residual = 0.000210269, Final residual = 8.30651e-07, No Iterations 3

ExecutionTime = 1.24 s  ClockTime = 1 s



Time = 300



smoothSolver:  Solving for Ux, Initial residual = 0.00257644, Final residual = 1.3703e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.0025284, Final residual = 1.34594e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00143646, Final residual = 8.13869e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00720673, Final residual = 2.29444e-05, No Iterations 2

time step continuity errors : sum local = 2.51267, global = 1.1237, cumulative = -1608.13

smoothSolver:  Solving for epsilon, Initial residual = 0.000206882, Final residual = 1.18584e-06, No Iterations 2

smoothSolver:  Solving for k, Initial residual = 0.000211105, Final residual = 7.97838e-07, No Iterations 3

ExecutionTime = 1.25 s  ClockTime = 1 s



End

I tried also with your folders (0 , constant, system) but I got the same error :( don't know why

Quote:

Originally Posted by immortality (Post 448685)

it works fine for me!

Code:

Time = 293



smoothSolver:  Solving for Ux, Initial residual = 0.00302422, Final residual = 1.63714e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00286406, Final residual = 1.40895e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00163624, Final residual = 9.37947e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00816459, Final residual = 4.26288e-05, No Iterations 2

time step continuity errors : sum local = 4.67106, global = -1.04718, cumulative = -1602.9

smoothSolver:  Solving for epsilon, Initial residual = 0.000209883, Final residual = 1.19376e-06, No Iterations 2

bounding epsilon, min: -138.041 max: 2.9235e+08 average: 9.09496e+06

smoothSolver:  Solving for k, Initial residual = 0.000215752, Final residual = 8.36107e-07, No Iterations 3

ExecutionTime = 1.22 s  ClockTime = 1 s



Time = 294



smoothSolver:  Solving for Ux, Initial residual = 0.00299998, Final residual = 1.35978e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00283594, Final residual = 1.3551e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00161325, Final residual = 8.1245e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00804053, Final residual = 1.87376e-05, No Iterations 2

time step continuity errors : sum local = 2.05042, global = -1.73826, cumulative = -1604.64

smoothSolver:  Solving for epsilon, Initial residual = 0.000211165, Final residual = 1.19876e-06, No Iterations 2

smoothSolver:  Solving for k, Initial residual = 0.000214596, Final residual = 8.11414e-07, No Iterations 3

ExecutionTime = 1.22 s  ClockTime = 1 s



Time = 295



smoothSolver:  Solving for Ux, Initial residual = 0.00269395, Final residual = 1.35094e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00268798, Final residual = 1.34244e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00138372, Final residual = 8.11865e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00734547, Final residual = 2.72228e-05, No Iterations 2

time step continuity errors : sum local = 2.98289, global = -1.95487, cumulative = -1606.59

smoothSolver:  Solving for epsilon, Initial residual = 0.000209892, Final residual = 1.16917e-06, No Iterations 2

smoothSolver:  Solving for k, Initial residual = 0.000211634, Final residual = 8.11507e-07, No Iterations 3

ExecutionTime = 1.23 s  ClockTime = 1 s



Time = 296



smoothSolver:  Solving for Ux, Initial residual = 0.00260183, Final residual = 1.34072e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00259837, Final residual = 1.34191e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00141887, Final residual = 8.04426e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00717931, Final residual = 2.6446e-05, No Iterations 2

time step continuity errors : sum local = 2.89427, global = -1.88796, cumulative = -1608.48

smoothSolver:  Solving for epsilon, Initial residual = 0.000205624, Final residual = 1.19057e-06, No Iterations 2

smoothSolver:  Solving for k, Initial residual = 0.000209961, Final residual = 7.95193e-07, No Iterations 3

ExecutionTime = 1.23 s  ClockTime = 1 s



Time = 297



smoothSolver:  Solving for Ux, Initial residual = 0.00246914, Final residual = 1.34119e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00247703, Final residual = 1.31678e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00133655, Final residual = 8.0358e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00733201, Final residual = 2.25581e-05, No Iterations 2

time step continuity errors : sum local = 2.47292, global = -1.16714, cumulative = -1609.65

smoothSolver:  Solving for epsilon, Initial residual = 0.000209202, Final residual = 1.17054e-06, No Iterations 2

bounding epsilon, min: -15.779 max: 2.92343e+08 average: 9.09302e+06

smoothSolver:  Solving for k, Initial residual = 0.000209555, Final residual = 8.24086e-07, No Iterations 3

ExecutionTime = 1.24 s  ClockTime = 1 s



Time = 298



smoothSolver:  Solving for Ux, Initial residual = 0.00262038, Final residual = 1.38132e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00253905, Final residual = 1.35259e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00141747, Final residual = 8.14806e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00717604, Final residual = 2.40071e-05, No Iterations 2

time step continuity errors : sum local = 2.62906, global = -0.296388, cumulative = -1609.95

smoothSolver:  Solving for epsilon, Initial residual = 0.00020553, Final residual = 1.17951e-06, No Iterations 2

smoothSolver:  Solving for k, Initial residual = 0.000212137, Final residual = 7.99386e-07, No Iterations 3

ExecutionTime = 1.24 s  ClockTime = 1 s



Time = 299



smoothSolver:  Solving for Ux, Initial residual = 0.00262036, Final residual = 1.34373e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.00243577, Final residual = 1.31502e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.0013605, Final residual = 8.10666e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00723721, Final residual = 1.94311e-05, No Iterations 2

time step continuity errors : sum local = 2.13106, global = 0.691298, cumulative = -1609.26

smoothSolver:  Solving for epsilon, Initial residual = 0.000209559, Final residual = 1.17209e-06, No Iterations 2

bounding epsilon, min: -24.6726 max: 2.92349e+08 average: 9.09526e+06

smoothSolver:  Solving for k, Initial residual = 0.000210269, Final residual = 8.30651e-07, No Iterations 3

ExecutionTime = 1.24 s  ClockTime = 1 s



Time = 300



smoothSolver:  Solving for Ux, Initial residual = 0.00257644, Final residual = 1.3703e-05, No Iterations 3

smoothSolver:  Solving for Uy, Initial residual = 0.0025284, Final residual = 1.34594e-05, No Iterations 3

smoothSolver:  Solving for Uz, Initial residual = 0.00143646, Final residual = 8.13869e-06, No Iterations 3

GAMG:  Solving for p, Initial residual = 0.00720673, Final residual = 2.29444e-05, No Iterations 2

time step continuity errors : sum local = 2.51267, global = 1.1237, cumulative = -1608.13

smoothSolver:  Solving for epsilon, Initial residual = 0.000206882, Final residual = 1.18584e-06, No Iterations 2

smoothSolver:  Solving for k, Initial residual = 0.000211105, Final residual = 7.97838e-07, No Iterations 3

ExecutionTime = 1.25 s  ClockTime = 1 s



End

you use 2.1.1 vesion? download 2.2.0 or 2.2.1 version so that can run.

Quote:

Originally Posted by immortality (Post 448691)

you use 2.1.1 vesion? download 2.2.0 or 2.2.1 version so that can run.

Dear Ehsan,

I used OpenFoam 2.2 and finally, I didn't have this error. But after launching Paraview, I realised that my problem didn't been solved. Here are the attachement that i got with elements at the interface between the inlet and wall with high pressure !

why youe mesh is so coarse? use a finer mesh and see if changes.