[uckmhnds]

Hello dear all,


I try to simulate tip leakage flow over a low pressure turbine blade for my thesis so that i set up my case in OpenFoam in order to proceed it. I completed my pre-processing steps (creating domain, mesh, initial files on OpenFoam etc.) and i used Pointwise to create my hybrid mesh. However, when i begin to run simulation i get the error, that is "floating point exception (core dumped)" all the time. I changed my mesh, reset the boundary conditions but still the same error for WEEKS. Could you please help me? I ve read many posts about the problem and they say there might be something which is divided by zero but i don't know what that is.

Courant Number mean: 7.53016e-08 max: 0.000123026
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 3.47215e-05, No Iterations 1000
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 3.51367e-05, No Iterations 1000
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 4.29613e-05, No Iterations 1000
#0 Foam::error; printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const at ??:?
#6 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#10 ? at ??:?
#11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12 ? at ??:?
Floating point exception (core dumped)

[uckmhnds]

Ok. I have another problem now. I ve noticed that i must specify "patch" for the periodic boundaries instead of "cyclic" on Pointwise. Then, i must create a "createPatchDict" file in the "system" folder and run "createPatch" on terminal to define those periodic boundries properly. However, i get another error while doing that.
Here is my createPatchDict

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object createPatchDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// Do a synchronisation of coupled points after creation of any patches.
// Note: this does not work with points that are on multiple coupled patches
// with transformations (i.e. cyclics).
pointSync false;

// Patches to create.
patches
(
{
//- Master side patch
name periodic1;
patchInfo
{
type cyclic;
matchTolerance 0.001;
neighbourPatch periodic2;
transform translational;
separationVector (0 0.158 0);
}
constructFrom patches;
patches (up);
}

{
//- Slave side patch
name periodic2;
patchInfo
{
type cyclic;
matchTolerance 0.001;
neighbourPatch periodic1;
transform translational;
separationVector (0 -0.158 0);
}
constructFrom patches;
patches (down);
}

)

Here is the boundary file in "polyMesh";

FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}

7
(
blade
{
type wall;
nFaces 5278;
startFace 1215527;
}
bottom
{
type wall;
nFaces 14653;
startFace 1220805;
}
down
{
type patch;
nFaces 6003;
startFace 1235458;
}
front
{
type symmetryPlane;
nFaces 14653;
startFace 1241461;
}
inlet
{
type patch;
nFaces 1421;
startFace 1256114;
}
outlet
{
type patch;
nFaces 2581;
startFace 1257535;
}
up
{
type patch;
nFaces 6003;
startFace 1260116;
}
)

Here is the error i'm getting

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 3.0.1-d8a290b55d28
Exec : createPatch
Date : Aug 13 2018
Time : 21:42:02
Host : "yavuzabd"
PID : 6964
Case : /home/yavuzabd/Desktop/coarse
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Reading createPatchDict



--> FOAM FATAL IO ERROR:
"ill defined primitiveEntry starting at keyword 'patches' on line 24 and ending at line 59"

file: /home/yavuzabd/Desktop/coarse/system/createPatchDict at line 59.

From function primitiveEntry::readEntry(const dictionary&, Istream&)
in file lnInclude/IOerror.C at line 132.

FOAM exiting

I barely comprehended how createPatch command works but still don't know why i get that error. Actually i ve changed my mesh many times so that i'm really bored of doing same things again again. At the end, i got just stupid different error which i don't know what they are about. OpenFOAM's error handling is terrible, just few lines but no clear explanations

[uckmhnds]

I solved the problem but do not know how to vanish that post. Or, i could help if someone has the same problem

[divyathokala01@gmail.com]

Originally Posted by uckmhnds

I solved the problem but do not know how to vanish that post. Or, i could help if someone has the same problem

Hi, Can you please help me? How could you resolve the problem? I am getting almost a similar error as follows.

Courant Number mean: 0.000152296 max: 0.0566087
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 4.18724e-06, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 5.88284e-06, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 3.85448e-06, No Iterations 3
Pressure gradient source: uncorrected Ubar = 4.39516e-05, pressure gradient = 24.0577
GAMG: Solving for p, Initial residual = 1, Final residual = 2.11191e+50, No Iterations 1000
time step continuity errors : sum local = 3.67866e+46, global = 4.86877e+39, cumulative = 4.86877e+39
Pressure gradient source: uncorrected Ubar = -1.03328e+44, pressure gradient = 1.24319e+46
GAMG: Solving for p, Initial residual = 0.233875, Final residual = 3.60805e+45, No Iterations 1000
time step continuity errors : sum local = 6.99856e+92, global = 5.74806e+87, cumulative = 5.74806e+87
Pressure gradient source: uncorrected Ubar = -9.96494e+92, pressure gradient = 1.19893e+95
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4 Foam::GaussSeidelSmoother::smooth(Foam::Field<doub le>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/pisoFoam"
#8 Foam::fvMatrix<double>::solve() in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/pisoFoam"
#9 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:?
#10 Foam::RASModels::realizableKE<Foam::Incompressible TurbulenceModel<Foam::transportModel> >::correct() at ??:?
#11 ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/pisoFoam"
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/pisoFoam"
Floating point exception (core dumped)