| 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? |
