cfdonline2mohsen |
October 21, 2013 09:28 |
and the final step is running pimpleFoam:
5th Step: Running pimpleFoam for My 1st Case (TwinLundWithCyclic2):
Code:
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic2$ pimpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : pimpleFoam
Date : Oct 21 2013
Time : 11:22:37
Host : "moh3n"
PID : 3554
Case : /home/moh3n/Desktop/920729/TwinLundWithCyclic2
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 turbulence model type LESModel
Selecting LES turbulence model dynamicSmagorinsky
dynamicSmagorinskyCoeffs
{
filter simple;
ce 1.048;
}
No field sources present
PIMPLE: Operating solver in PISO mode
Starting time loop
Courant Number mean: 0.00064134496 max: 0.00066433568
Time = 1.5e-05
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 5.4453587e-07, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.99448474, Final residual = 2.7934192e-08, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.98434205, Final residual = 5.1759288e-08, No Iterations 2
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.041080613, No Iterations 184
time step continuity errors : sum local = 1.9183968e-07, global = -3.7234788e-10, cumulative = -3.7234788e-10
DICPCG: Solving for p, Initial residual = 0.0036086778, Final residual = 9.0427381e-07, No Iterations 255
time step continuity errors : sum local = 3.3297187e-10, global = -2.0421325e-12, cumulative = -3.7439002e-10
bounding k, min: -3.4336254e-06 max: 0.93980124 average: 0.001494962
ExecutionTime = 23.49 s ClockTime = 24 s
Courant Number mean: 0.00095867724 max: 0.0032413083
Time = 3e-05
DILUPBiCG: Solving for Ux, Initial residual = 0.81314678, Final residual = 5.0166127e-06, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.31087511, Final residual = 2.1087205e-08, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.32819123, Final residual = 7.1596872e-07, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.014583725, Final residual = 0.00072210423, No Iterations 193
time step continuity errors : sum local = 2.9425828e-07, global = -5.9706547e-09, cumulative = -6.3450447e-09
DICPCG: Solving for p, Initial residual = 0.0022751552, Final residual = 9.142224e-07, No Iterations 240
time step continuity errors : sum local = 3.0511796e-10, global = -8.1926347e-13, cumulative = -6.345864e-09
bounding k, min: -0.00020715488 max: 0.88381996 average: 0.0024604238
ExecutionTime = 43.1 s ClockTime = 43 s
Courant Number mean: 0.0012224251 max: 0.003201858
Time = 4.5e-05
DILUPBiCG: Solving for Ux, Initial residual = 0.86187261, Final residual = 6.7931013e-06, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.45940243, Final residual = 2.8525143e-08, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.38062296, Final residual = 6.9428457e-07, No Iterations 1
And it works just fine for my 1st case both in Parallel & Serial.
5th Step: Running pimpleFoam for My 2nd Case (TwinLundWithCyclic3):
Code:
moh3n@moh3n:~/Desktop/920729/TwinLundWithCyclic3$ pimpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : pimpleFoam
Date : Oct 21 2013
Time : 11:22:54
Host : "moh3n"
PID : 3555
Case : /home/moh3n/Desktop/920729/TwinLundWithCyclic3
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 turbulence model type LESModel
Selecting LES turbulence model dynamicSmagorinsky
dynamicSmagorinskyCoeffs
{
filter simple;
ce 1.048;
}
No field sources present
PIMPLE: Operating solver in PISO mode
Starting time loop
Courant Number mean: 0.00064134496 max: 0.00066433568
Time = 1.5e-05
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::scaledMappedVelocityFixedValueFvPatchField::updateCoeffs() in "/home/moh3n/OpenFOAM/moh3n-2.1.1/platforms/linux64GccDPOpt/lib/libscaledMappedVelocityFixedValueKataoka.so"
#4 Foam::fvMatrix<Foam::Vector<double> >::fvMatrix(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pimpleFoam"
#5 Foam::fv::gaussLaplacianScheme<Foam::Vector<double>, double>::fvmLaplacianUncorrected(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6 Foam::fv::gaussLaplacianScheme<Foam::Vector<double>, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fv::laplacianScheme<Foam::Vector<double>, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so"
#9 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so"
#10 Foam::incompressible::LESModels::GenEddyVisc::divDevBeff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so"
#11 Foam::incompressible::LESModel::divDevReff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so"
#12
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pimpleFoam"
#13 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/pimpleFoam"
Floating point exception
And it gave Floating point exception Error!!!:mad::mad::mad:
The 0, system and constant(except polyMesh) folders are the same for both cases. and the geometries are the same for both cases except little differences in checkMesh ( using e-17 or e-16 by OpenFOAM instead of 0).
My 1st case runs just fine and my 2nd case gives Floating point exception error!!
I've also uploaded both cases:
http://www.mediafire.com/?m6gd1m06ro88bmj
You can do all the above steps yourself to see that I'm right!
It drove me crazy!
Some body please tell me what's the reason???
|