CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Same SimpleFOAM Case converges with openFOAM 2.1 but diverges with openFOAM 2.0.1 (http://www.cfd-online.com/Forums/openfoam-solving/108372-same-simplefoam-case-converges-openfoam-2-1-but-diverges-openfoam-2-0-1-a.html)

alsdia October 21, 2012 23:07

Same SimpleFOAM Case converges with openFOAM 2.1 but diverges with openFOAM 2.0.1
 
I have two machines
1) OpenSUSE 12.1 with OpenFOAM 2.1.x
[http://susestudio.com/a/2qtLK2/geekocfd?version=3.1.0]
2) CAELINUX2011 launched inside VirtualBOX on a XP 64 machine with openFOAM v2.0.1
http://www.caelinux.com/CMS/index.ph...d=69&Itemid=58


I launch the following case of a sphere invested by a flow of air t 14.6 m/s
https://www.dropbox.com/sh/gg1ypj0c1...here.case.test
https://www.dropbox.com/sh/gg1ypj0c1...lome.ball.jpeg

1) in OpenSUSE with with OpenFOAM 2.1.x, I get convergence and no error

Code:

/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.1.x                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 2.1.x-63f6b7e20176
Exec  : simpleFoam
Date  : Oct 22 2012
Time  : 03:00:21
Host  : "linux-baef"
PID    : 18519
Case  : /home/serena/sphere.medium.40000vertex.withCoeffs.USB.test
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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 kOmegaSST
kOmegaSSTCoeffs
{
    alphaK1        0.85034;
    alphaK2        1;
    alphaOmega1    0.5;
    alphaOmega2    0.85616;
    gamma1          0.5532;
    gamma2          0.4403;
    beta1          0.075;
    beta2          0.0828;
    betaStar        0.09;
    a1              0.31;
    c1              10;
}

No field sources present


SIMPLE: no convergence criteria found. Calculations will run for 105 steps.


Starting time loop

Reading surface description:
    yNormal

Time = 1

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.0576103, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.051907, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.0754043, No Iterations 3
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0981085, No Iterations 5
time step continuity errors : sum local = 0.100654, global = -0.00926128, cumulative = -0.00926128
smoothSolver:  Solving for omega, Initial residual = 0.0394307, Final residual = 0.00164586, No Iterations 3
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 0.0804227, No Iterations 3
ExecutionTime = 1 s  ClockTime = 1 s

forceCoeffs output:
    Cd = 0.945971
    Cl = 0.00211079
    Cm = -0.426151
    Cl(f) = 0.427207
    Cl(r) = -0.425096

....
Time = 104

smoothSolver:  Solving for Ux, Initial residual = 0.00150213, Final residual = 7.9976e-05, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.00335438, Final residual = 0.000166074, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 0.00404789, Final residual = 0.000185933, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0189557, Final residual = 0.00182748, No Iterations 1
time step continuity errors : sum local = 0.00550264, global = -3.908e-05, cumulative = -0.0017971
smoothSolver:  Solving for omega, Initial residual = 0.000176399, Final residual = 1.75706e-05, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.00183807, Final residual = 6.93889e-05, No Iterations 3
ExecutionTime = 35.54 s  ClockTime = 36 s

forceCoeffs output:
    Cd = 0.00226424
    Cl = -0.00123669
    Cm = -0.000745832
    Cl(f) = 0.000127488
    Cl(r) = -0.00136418

Time = 105

smoothSolver:  Solving for Ux, Initial residual = 0.00146514, Final residual = 7.74717e-05, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.00327494, Final residual = 0.000160581, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 0.00394729, Final residual = 0.000181077, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0187608, Final residual = 0.00179919, No Iterations 1
time step continuity errors : sum local = 0.00541478, global = -0.000123, cumulative = -0.0019201
smoothSolver:  Solving for omega, Initial residual = 0.000174777, Final residual = 1.71551e-05, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.0018046, Final residual = 6.76978e-05, No Iterations 3
ExecutionTime = 35.83 s  ClockTime = 36 s

forceCoeffs output:
    Cd = 0.00225507
    Cl = -0.00125446
    Cm = -0.000737782
    Cl(f) = 0.000110551
    Cl(r) = -0.00136501

End

2) Now same files launched in CAELINUX con OF 2.0.1 ( same files nothing changed!!) I get divergence...
Code:

/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.0.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.com                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 2.0.1-51f1de99a4bc
Exec  : simpleFoam
Date  : Oct 22 2012
Time  : 09:43:12
Host  : orion
PID    : 2050
Case  : /media/sf_My_Documents/linux.Ferries.Carrier/sphere.medium.40000vertex.withCoeffs.USB.test
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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 kOmegaSST
kOmegaSSTCoeffs
{
    alphaK1        0.85034;
    alphaK2        1;
    alphaOmega1    0.5;
    alphaOmega2    0.85616;
    gamma1          0.5532;
    gamma2          0.4403;
    beta1          0.075;
    beta2          0.0828;
    betaStar        0.09;
    a1              0.31;
    c1              10;
}


SIMPLE: no convergence criteria found. Calculations will run for 105 steps.


Starting time loop

Reading surface description:
    yNormal

Time = 1

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.0997099, No Iterations 223
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0996995, No Iterations 23
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.09999, No Iterations 336
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0867693, No Iterations 4
time step continuity errors : sum local = 0.37954, global = -0.0213585, cumulative = -0.0213585
smoothSolver:  Solving for omega, Initial residual = 0.370863, Final residual = 0.0369602, No Iterations 151
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 0.0978377, No Iterations 63
ExecutionTime = 10.48 s  ClockTime = 14 s

forceCoeffs output:
    Cd = 1.60633
    Cl = 0.00632664
    Cm = -0.724241
....since the beginning the Cd is different

Time = 8

smoothSolver:  Solving for Ux, Initial residual = 0.999943, Final residual = 0.0719959, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.999943, Final residual = 0.0524752, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 0.999946, Final residual = 0.0524753, No Iterations 3
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  in "/lib/libc.so.6"
#3  Foam::GAMGSolver::scalingFactor(Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/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 "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::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 "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8  Foam::fvMatrix<double>::solve() in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/simpleFoam"
#9 
 in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/simpleFoam"
#10  __libc_start_main in "/lib/libc.so.6"
#11 
 in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/simpleFoam"
Floating point exception

I did not change nothing in the input files. But the results on the two machines are different. in openSuse converge, in Caelinux diverges.
Why? is it because in Opensuse I have Build : 2.1.x-63f6b7e20176 and in CAELINUX Build : 2.0.1-51f1de99a4bc? Or is something else?

alberto October 22, 2012 02:24

Hi Serena,

first of all, thank you for using GeekoCFD. Does it work well for your needs?

I am not sure why you do not obtain a converged solution with openfoam 2.0.1, but it is possible that some bug was fixed in 2.1, which might explan the difference you notice. That's why it is often convenient to use the latest version of OF, or the git version.

Best,

alsdia October 22, 2012 02:58

Alberto

thank you for your quick reply.
Yes your GeekoCFD works beautifully, it is a very convenient tool and I love also the choice of OpenSuse instead of Ubuntu, more sharp and well focused environment.
By the way, I would like to ask you, when you submit a new Geeko release: could you include in the scientific software also :
1) SALOME (a nice alternative to gmsh)
http://www.salome-platform.org/
2) SCILAB (a nice alternative to R)
http://www.scilab.org/

Bye

Serena

alberto October 22, 2012 11:25

Hi Serena,

thank you for your feedback.

Quote:

Originally Posted by alsdia (Post 387825)
1) SALOME (a nice alternative to gmsh)
http://www.salome-platform.org/

When I started building GeekoCFD I thought of including SALOME. However, building it from source is a major effort. SALOME also provides a very large binary package that I could add "as is" to GeekoCFD, since it works out of the box. I thought it would be better to leave that choice to users, if they need SALOME, since adding it increases the image size significantly. However, if there is enough demand, I will re-consider :-)

Quote:

...
2) SCILAB (a nice alternative to R)
http://www.scilab.org/
I did not include SCILAB because of its license. However I have just read they changed the license to a GPL 2 compatible one. I will consider its inclusion.

Thank you for your suggestions


All times are GMT -4. The time now is 03:25.