CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Same SimpleFOAM Case converges with openFOAM 2.1 but diverges with openFOAM 2.0.1

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   October 21, 2012, 23:07
Default Same SimpleFOAM Case converges with openFOAM 2.1 but diverges with openFOAM 2.0.1
  #1
New Member
 
edy
Join Date: Jul 2012
Location: Tokyo
Posts: 6
Rep Power: 0
alsdia is on a distinguished road
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


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?
alsdia is offline   Reply With Quote

Old   October 22, 2012, 02:24
Default
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
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,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   October 22, 2012, 02:58
Default
  #3
New Member
 
edy
Join Date: Jul 2012
Location: Tokyo
Posts: 6
Rep Power: 0
alsdia is on a distinguished road
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
alsdia is offline   Reply With Quote

Old   October 22, 2012, 11:25
Default
  #4
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Hi Serena,

thank you for your feedback.

Quote:
Originally Posted by alsdia View Post
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
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
SimpleFoam case with SpalartAllmaras turbulence model implemented nedved OpenFOAM Running, Solving & CFD 2 November 30, 2014 23:43
64bitrhel5 OF installation instructions mirko OpenFOAM Installation 2 August 12, 2008 18:07
Is Playstation 3 cluster suitable for CFD work hsieh OpenFOAM 8 March 12, 2008 05:58
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 12:24


All times are GMT -4. The time now is 08:41.