[imani]

hi
who can help me
what is this error...


Time = 1

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.00995599, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0104898, No Iterations 1
DILUPBiCG: Solving for T, Initial residual = 1, Final residual = 0.0591834, No Iterations 1
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 7.1898, No Iterations 1001
time step continuity errors : sum local = 84.0209, global = 0.0235848, cumulative = 0.0235848
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 at kEpsilon.C:0
#6 Foam::incompressible::RASModels::kEpsilon::correct () in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#7
in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/buoyantBoussinesqSimpleFoam"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/buoyantBoussinesqSimpleFoam"
Floating point exception (core dumped)



i think my problem is that my mesh is very big....because my room is very large ,when i fine my meshs it can not make blockMesh successfully
what should i do
thanks alot

[Lieven]

Hi imani,

What do you mean exactly when you say "my mesh is very big". Do you mean a large number of cells or very big dimensions?
Can you post the output of checkMesh?

Regarding the error you get, there is a division by zero somewhere in the run. Since you are using the kEpsilon model, my first guess will be that epsilon is set to 0 in your initial conditions. Set it to any other positive value (but best is to make it small, 1e-3 or so).

Which solver are you using? Just asking since your pressure equation is messed up...

Cheers,

Lieven

[imani]

hi dear
my dimension is a room 8*6*2.7 meter
this room has windows with gap
these gap are 2 mm and around the windows
when i fine my meshes i have problem
my solver is buoyantBossinesqueSimpleFoam

[Lieven]

Hi Imani,

Did you try to modify the epsilon-initial conditions? Could you please give the output of checkMesh? At the moment, you really didn't give enough to be able to help you properly...

[imani]

hi dear lieven
thank you

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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 : checkMesh
Date : Jun 21 2014
Time : 08:24:10
Host : "ubuntu"
PID : 2679
Case : /home/dell/Desktop/zohreh
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 polyMesh for time = 0

Time = 0

Mesh stats
points: 139468
faces: 398460
internal faces: 379140
cells: 129600
boundary patches: 7
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 129600
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0

Checking topology...
Boundary definition OK.
***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface topology
fixedWall-1 1620 1708 ok (non-closed singly connected)
fixedWall-2 2160 2268 ok (non-closed singly connected)
fixedWall-3 1620 1736 ok (non-closed singly connected)
fixedWall-4 2160 2268 ok (non-closed singly connected)
floor 4800 4981 ok (non-closed singly connected)
ceiling 4800 4981 ok (non-closed singly connected)
defaultFaces 2160 2268 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (0 0 0) (8 2.7 6)
Mesh (non-empty, non-wedge) directions (1 1 0)
Mesh (non-empty) directions (1 1 0)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (2.80429e-16 7.81092e-16 5.17203e-16) OK.
Max cell openness = 1.66463e-16 OK.
Max aspect ratio = 1.44416 OK.
Minumum face area = 0.00479307. Maximum face area = 0.0191723. Face area magnitudes OK.
Min volume = 0.000479307. Max volume = 0.00191723. Total volume = 129.6. Cell volumes OK.
Mesh non-orthogonality Max: 0 average: 0
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 5.2182e-12 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End

dell@ubuntu:~/Desktop/zohreh$ buoyantBoussinesqSimpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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 : buoyantBoussinesqSimpleFoam
Date : Jun 21 2014
Time : 08:24:20
Host : "ubuntu"
PID : 2680
Case : /home/dell/Desktop/zohreh
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 g
Reading thermophysical properties

Reading field T

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Creating turbulence model

Selecting RAS turbulence model kEpsilon
--> FOAM Warning :
From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(c onst fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
in file groovyBCFvPatchField.C at line 119
No value defined for k on fixedWall-4 therefore using 2160{0}
bounding k, min: 0 max: 0.1 average: 0.1
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
sigmaEps 1.3;
}

Reading field kappat

Calculating field g.h


SIMPLE: convergence criteria
field p_rgh tolerance 0.01
field U tolerance 0.0001
field T tolerance 0.01
field "(k|epsilon|omega)" tolerance 0.001


Starting time loop

Time = 1

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.00995599, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0104898, No Iterations 1
DILUPBiCG: Solving for T, Initial residual = 1, Final residual = 0.0591834, No Iterations 1
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 7.1898, No Iterations 1001
time step continuity errors : sum local = 84.0209, global = 0.0235848, cumulative = 0.0235848
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 at kEpsilon.C:0
#6 Foam::incompressible::RASModels::kEpsilon::correct () in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#7
in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/buoyantBoussinesqSimpleFoam"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/buoyantBoussinesqSimpleFoam"
Floating point exception (core dumped)
dell@ubuntu:~/Desktop/zohreh$


have good day
imani