[adarsh tiwari]

Hi all,

While simulating for the problem, i got the following error,

eatin@EAT-Standalone:~/ADARSH/cavity$ rhoSimpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.2 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.2-9240f8b967db
Exec : rhoSimpleFoam
Date : Feb 03 2014
Time : 16:21:26
Host : "EAT-Standalone"
PID : 7620
Case : /home/eatin/ADARSH/cavity
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


SIMPLE: convergence criteria
field p tolerance 1e-06
field U tolerance 1e-06
field h tolerance 1e-06
field k tolerance 1e-06
field omega tolerance 1e-06

Reading thermophysical properties

Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
alphaK1 0.85034;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.85616;
Prt 1;
gamma1 0.5532;
gamma2 0.4403;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
}

No finite volume options present


Starting time loop

Time = 0.005

GAMG: Solving for Ux, Initial residual = 1, Final residual = 9.76443e-07, No Iterations 115
GAMG: Solving for Uy, Initial residual = 1, Final residual = 9.71495e-07, No Iterations 61
GAMG: Solving for Uz, Initial residual = 1, Final residual = 9.75421e-07, No Iterations 116
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 8.10833e-07, No Iterations 69
GAMG: Solving for p, Initial residual = 0.999137, Final residual = 0.00126243, No Iterations 1000
time step continuity errors : sum local = 6.93746e-12, global = 3.97247e-14, cumulative = 3.97247e-14
rho max/min : 1.5 1.36693
GAMG: Solving for omega, Initial residual = 0.213657, Final residual = 2.02712e-08, No Iterations 2
GAMG: Solving for k, Initial residual = 1, Final residual = 1.7107e-07, No Iterations 2
ExecutionTime = 72.46 s ClockTime = 72 s

Time = 0.01

GAMG: Solving for Ux, Initial residual = 0.852602, Final residual = 9.3075e-07, No Iterations 59
GAMG: Solving for Uy, Initial residual = 0.866457, Final residual = 8.52692e-07, No Iterations 57
GAMG: Solving for Uz, Initial residual = 0.875269, Final residual = 9.25048e-07, No Iterations 63
DILUPBiCG: Solving for h, Initial residual = 0.116407, Final residual = 8.73549e-07, No Iterations 61
GAMG: Solving for p, Initial residual = 0.660759, Final residual = 6.79783e-07, No Iterations 1000
time step continuity errors : sum local = 6.85625e-12, global = -9.17279e-15, cumulative = 3.05519e-14
rho max/min : 1.5 1.37284
GAMG: Solving for omega, Initial residual = 0.173706, Final residual = 4.94879e-10, No Iterations 2
GAMG: Solving for k, Initial residual = 0.979579, Final residual = 3.70192e-10, No Iterations 2
ExecutionTime = 139.7 s ClockTime = 140 s

Time = 0.015

GAMG: Solving for Ux, Initial residual = 0.885481, Final residual = 9.79136e-07, No Iterations 68
GAMG: Solving for Uy, Initial residual = 0.924134, Final residual = 9.63321e-07, No Iterations 65
GAMG: Solving for Uz, Initial residual = 0.90357, Final residual = 8.69081e-07, No Iterations 69
DILUPBiCG: Solving for h, Initial residual = 0.0601542, Final residual = 4.73406e-07, No Iterations 64
GAMG: Solving for p, Initial residual = 0.814731, Final residual = 9.26881e-09, No Iterations 143
time step continuity errors : sum local = 3.72959e-11, global = 1.38386e-13, cumulative = 1.68938e-13
rho max/min : 1.5 1.37284
GAMG: Solving for omega, Initial residual = 0.208416, Final residual = 1.47666e-08, No Iterations 2
GAMG: Solving for k, Initial residual = 0.970247, Final residual = 9.02517e-08, No Iterations 2
ExecutionTime = 160.3 s ClockTime = 160 s

Time = 0.02

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#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) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::GaussSeidelSmoother::smooth(Foam::Field<doub le>&, Foam::Field<double> const&, unsigned char, int) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 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::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#7
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/rhoSimpleFoam"
#8
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/rhoSimpleFoam"
#9
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/rhoSimpleFoam"
#10
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/rhoSimpleFoam"
#11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/rhoSimpleFoam"
Floating point exception (core dumped)




Initially i was solving with smoothSolver but later I have chaged the solver from smoothSolver to GAMG, also reduced the relaxation factor, but still getting the same O/P

[alexeym]

Hi,

Lots of reasons may lead to this error (bad IC/BC, bad settings for GAMG solver etc).

Try using PBiCG instead of GAMG (or maybe first use DICGaussSeidel as a smoother for GAMG, as the error was during smooth method call).

[adarsh tiwari]

greetings Alexeym,

First of all thanks for your quick response.
I have tried all the ways you mentioned, and again faced the same problem. later I changed the pressure solver to PCG, with preconditioner DIC, but again the same issue arises. I have pointed out few interesting things in the solution:

Create time

Create mesh for time = 0


SIMPLE: convergence criteria
field p tolerance 1e-06
field U tolerance 1e-06
field h tolerance 1e-06
field k tolerance 1e-06
field omega tolerance 1e-06

Reading thermophysical properties

Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
alphaK1 0.85034;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.85616;
Prt 1;
gamma1 0.5532;
gamma2 0.4403;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
}

No finite volume options present


Starting time loop

Time = 0.005

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 9.4729e-07, No Iterations 49
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 8.8146e-07, No Iterations 56
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 8.02978e-07, No Iterations 48
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 8.10833e-07, No Iterations 69
DICPCG: Solving for p, Initial residual = 0.999137, Final residual = 9.13645e-07, No Iterations 211
DICPCG: Solving for p, Initial residual = 0.00897964, Final residual = 8.85922e-07, No Iterations 166
DICPCG: Solving for p, Initial residual = 0.00164253, Final residual = 9.66848e-07, No Iterations 157
DICPCG: Solving for p, Initial residual = 0.00077649, Final residual = 9.73811e-07, No Iterations 145
DICPCG: Solving for p, Initial residual = 0.000255178, Final residual = 9.7144e-07, No Iterations 144
time step continuity errors : sum local = 4.81219e-11, global = 1.99058e-13, cumulative = 1.99058e-13
rho max/min : 1.5 1.36693
smoothSolver: Solving for omega, Initial residual = 0.213657, Final residual = 4.49362e-13, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 3.25225e-08, No Iterations 2
ExecutionTime = 30.38 s ClockTime = 31 s

Time = 0.01

DILUPBiCG: Solving for Ux, Initial residual = 0.853572, Final residual = 9.17802e-07, No Iterations 73
DILUPBiCG: Solving for Uy, Initial residual = 0.867465, Final residual = 9.01878e-07, No Iterations 76
DILUPBiCG: Solving for Uz, Initial residual = 0.878958, Final residual = 9.35019e-07, No Iterations 74
DILUPBiCG: Solving for h, Initial residual = 0.116405, Final residual = 8.7348e-07, No Iterations 61
DICPCG: Solving for p, Initial residual = 0.635366, Final residual = 9.3124e-07, No Iterations 217
DICPCG: Solving for p, Initial residual = 0.0104651, Final residual = 9.48305e-07, No Iterations 168
DICPCG: Solving for p, Initial residual = 0.00204681, Final residual = 9.85206e-07, No Iterations 158
DICPCG: Solving for p, Initial residual = 0.000909866, Final residual = 9.84431e-07, No Iterations 149
DICPCG: Solving for p, Initial residual = 0.000311013, Final residual = 9.3295e-07, No Iterations 144
time step continuity errors : sum local = 7.74558e-10, global = -4.21626e-11, cumulative = -4.19635e-11
rho max/min : 1.5 1.37284
smoothSolver: Solving for omega, Initial residual = 0.173706, Final residual = 2.81439e-14, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.979579, Final residual = 2.70355e-09, No Iterations 2
ExecutionTime = 60.54 s ClockTime = 61 s

Time = 0.015

DILUPBiCG: Solving for Ux, Initial residual = 0.882737, Final residual = 7.38937e-07, No Iterations 104
DILUPBiCG: Solving for Uy, Initial residual = 0.923639, Final residual = 9.97728e-07, No Iterations 107
DILUPBiCG: Solving for Uz, Initial residual = 0.910762, Final residual = 6.60213e-07, No Iterations 102
DILUPBiCG: Solving for h, Initial residual = 0.0603949, Final residual = 9.20886e-07, No Iterations 62
DICPCG: Solving for p, Initial residual = 0.791538, Final residual = 9.74006e-07, No Iterations 219
DICPCG: Solving for p, Initial residual = 0.0109548, Final residual = 9.2436e-07, No Iterations 170
DICPCG: Solving for p, Initial residual = 0.00235946, Final residual = 9.63775e-07, No Iterations 159
DICPCG: Solving for p, Initial residual = 0.00100188, Final residual = 9.7107e-07, No Iterations 150
DICPCG: Solving for p, Initial residual = 0.000361423, Final residual = 9.80694e-07, No Iterations 144
time step continuity errors : sum local = 2.57303e-07, global = 8.75624e-09, cumulative = 8.71428e-09
rho max/min : 1.5 1.37284
smoothSolver: Solving for omega, Initial residual = 0.208342, Final residual = 1.51846e-10, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.970282, Final residual = 9.42625e-08, No Iterations 2
ExecutionTime = 93.65 s ClockTime = 94 s

Time = 0.02

DILUPBiCG: Solving for Ux, Initial residual = 0.751383, Final residual = 34.257, No Iterations 1001
DILUPBiCG: Solving for Uy, Initial residual = 0.822324, Final residual = 857.331, No Iterations 1001
DILUPBiCG: Solving for Uz, Initial residual = 0.777107, Final residual = 71.3373, No Iterations 1001
DILUPBiCG: Solving for h, Initial residual = 0.999996, Final residual = 8.43945e-07, No Iterations 173
DICPCG: Solving for p, Initial residual = 0.994674, Final residual = 1.2802, No Iterations 1001
DICPCG: Solving for p, Initial residual = 0.44538, Final residual = 1.03715, No Iterations 1001
DICPCG: Solving for p, Initial residual = 0.622246, Final residual = 4.07268, No Iterations 1001
DICPCG: Solving for p, Initial residual = 0.815726, Final residual = 3.39158, No Iterations 1001
DICPCG: Solving for p, Initial residual = 0.849077, Final residual = 1.77735, No Iterations 1001
time step continuity errors : sum local = 15810.8, global = 0.0563862, cumulative = 0.0563862
rho max/min : 1.5 0.5
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 at smoothSolver.C:0
#4 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/rhoSimpleFoam"
#7 Foam::fvMatrix<double>::solve() in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#8 Foam::compressible::RASModels::kOmegaSST::correct( ) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#9
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/rhoSimpleFoam"
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/rhoSimpleFoam"
Floating point exception (core dumped)



here we can see that at time t=0.02 time step continuity error increases dramatically and hence everything messed up, also the number of iterations are also increases drastically.

I have also tried to solve by keeping turbulence model off, followed by reducing the relaxation factor and increasing the orthogonalCorrectors


and very sad to say that nothing works, it is facing the same problem

[alexeym]

Hi,

Surely I can keep guessing, you can keep checking if my guess is wrong or not. But simpler solution is:

1. Post your case files
2. If 1 is for some reason unacceptable, post (in CODE blocks or as attached files):
- short case description
- checkMesh output
- boundary conditions
- initial conditions
- fvSchemes & fvSolution

[adarsh tiwari]

hi

actually i had imported the mesh file from fluent, later i modified the case files of the 'cavity' to make it for my case. i have uploaded the new case files here. after looking at the simulation i can say that the problem may be with pressure entry only.

i have changed the solvers from GAMG to smoothSolver and what not .
i have also switched off the turbulence for each solvers and later again turned on. but i don't know what is preventing it to get the solutions.

the checkMesh output is written here.


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

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 341856
faces: 1000065
internal faces: 974907
cells: 329162
faces per cell: 6
boundary patches: 4
point zones: 0
face zones: 0
cell zones: 0

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

Checking topology...
Boundary definition OK.
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
HOUSING 13658 13832 ok (non-closed singly connected)
VALVE 5076 5106 ok (non-closed singly connected)
INLET 5060 5178 ok (non-closed singly connected)
OUTLET 1364 1394 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-12.15 -31 0.108183) (12.15 11 32.25)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-2.43187e-16 1.67233e-17 -4.02264e-16) OK.
Max cell openness = 1.92396e-14 OK.
Max aspect ratio = 624.06 OK.
Minimum face area = 7.0046e-05. Maximum face area = 0.767025. Face area magnitudes OK.
Min volume = 1.51164e-05. Max volume = 0.291075. Total volume = 16674.7. Cell volumes OK.
Mesh non-orthogonality Max: 65.3276 average: 15.3553
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 2.55632 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End

[alexeym]

Hi,

I don't know why you've decided that it will be very comfortable to read files from screenshots. It was actually rather annoying.

You've got mesh with max non-orthogonality = 65, so
1. Increase the number of nNonOrthogonalCorrectors
2. Use cellMDLimited and faceMDLimited schemes for gradient schemes
3. Switch from corrected to limited in laplacian and snGrad schemes
4. Reduce time step

[adarsh tiwari]

Hi Alexeym,

First of all sorry for the inconvinience caused.
as per the instructions i have changed the required entries and the files are mentioned here

[controlDict]

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

application rhoSimpleFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 1;

deltaT 1e-4;

writeControl timeStep;

writeInterval 100;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

runTimeModifiable yes;


// ************************************************** *********************** //]

[fvSolutions]

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

solvers
{
p
{
solver GAMG;
tolerance 1e-08;
relTol 0;
smoother GaussSeidel;
cacheAgglomeartion true;
nCellsInCoarsestLevel 100;
agglomerator faceAreaPair;
mergeLevels 1;
}

U
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-06;
relTol 0;
}

h
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}

k
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-08;
relTol 0;
}

omega
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-08;
relTol 0;
}

}

SIMPLE
{
nNonOrthogonalCorrectors 10; //0;
nCorrectors 10;
rhoMin rhoMin [1 -3 0 0 0] 0.5;
rhoMax rhoMax [1 -3 0 0 0] 1.5;

residualControl
{
p 1e-6;
U 1e-6;
h 1e-6;
k 1e-6;
omega 1e-6;
}
}


relaxationFactor
{
p 0.02;
rho 0.02;
U 0.02;
k 0.02;
omega 0.02;
h 0.02;
}
// ************************************************** *********************** //]

[fvSchemes]

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

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default cellLimited; //Gauss linear;
grad(U) cellMDLimited; //Gauss linear;
grad(p) faceMDLimited; //Gauss linear;
}

divSchemes
{
div(phi,U) bounded Gauss;
div((muEff*dev2(T(grad(U))))) Gauss linear;
div(phi,h) Gauss upwind;
div(phi,omega) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,K) Gauss upwind;
div(phid,p) Gauss upwind;
div(U,p) Gauss upwind;

//div(phi,U) bounded Gauss upwind;
//div((muEff*dev2(T(grad(U))))) Gauss linear;
//div(phi,e) bounded Gauss upwind;
//div(phi,epsilon) bounded Gauss upwind;
//div(phi,k) bounded Gauss upwind;
//div(phi,Ekp) bounded Gauss upwind;
}

laplacianSchemes
{
laplacian(muEff,U) Gauss linear limited; //corrected;
laplacian(alphaEff,h) Gauss linear limited; //corrected
laplacian((rho|A(U)),p) Gauss linear limited; //corrected
laplacian((rho*rAU),p) Gauss linear limited; //corrected
laplacian(DomegaEff,omega) Gauss linear limited; //corrected
laplacian(DkEff,k) Gauss linear limited; //corrected
laplacian(1,p) Gauss linear limited; //corrected
laplacian((rho*(1|A(U))),p) Gauss linear limited; //corrected

//laplacian(alphaEff,e) Gauss linear corrected;
//laplacian(muEff,U) Gauss linear corrected;
//laplacian(alphaEff,e) Gauss linear corrected;
//laplacian((rho*(1|A(U))),p) Gauss linear corrected;
//laplacian(DepsilonEff,epsilon) Gauss linear corrected;
//laplacian(DkEff,k) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
div(U,p) upwind phi;
}

snGradSchemes
{
default limited;
}

fluxRequired
{
default no;
p ;
}

// ************************************************** *********************** //]

while simulating I got the following error

[o/p in terminal][/eatin@EAT-Standalone:~$ cd ADARSH/cavity/
eatin@EAT-Standalone:~/ADARSH/cavity$ rhoSimpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.2 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.2-9240f8b967db
Exec : rhoSimpleFoam
Date : Feb 05 2014
Time : 16:43:18
Host : "EAT-Standalone"
PID : 14673
Case : /home/eatin/ADARSH/cavity
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


SIMPLE: convergence criteria
field p tolerance 1e-06
field U tolerance 1e-06
field h tolerance 1e-06
field k tolerance 1e-06
field omega tolerance 1e-06

Reading thermophysical properties

Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kOmegaSST


--> FOAM FATAL IO ERROR:
Grad scheme not specified

Valid grad schemes are :

9
(
Gauss
cellLimited
cellMDLimited
edgeCellsLeastSquares
faceLimited
faceMDLimited
fourth
leastSquares
pointCellsLeastSquares
)


file: /home/eatin/ADARSH/cavity/system/fvSchemes.gradSchemes.grad(U) at line 26.

From function gradScheme<Type>::New(const fvMesh& mesh, Istream& schemeData)
in file /home/opencfd/OpenFOAM/OpenFOAM-2.2.2/src/finiteVolume/lnInclude/gradScheme.C at line 54.

FOAM exiting

]

I have changed several times the same entries, also sometimes the schemes mentioned but the error got worse <sad >

thanks in advance,
Adarsh

[alexeym]

Hi,

If you read a little bit of documentation (for example here http://openfoam.org/docs/user/fvSchemes.php), you'd know that you should put in fvSchemes not just

cellLimited

but

cellLimited Gauss linear 1;

Also you can just attach your files instead of posting them in the body of a message.

[adarsh tiwari]

hi Alexeym,

initially it seems to be working but at time-step no. 50 it stopped and gave a message of 'core dumped'.
I have tried to solve the issue simply by editing the tutorial files and you would be happy to know that its workig fine.

one most importat thing to mention, THANK YOU FOR YOUR QUICK AND VALUABLE SUGGESTIONS.


with best regards,
Adarsh

[tevid21]

Hi, i'm trying to simulate a wing and i have this problem...
I do checkMesh and is ok, so i don't know where is the problem. I try de solutions that are in this post but i can't simulate. Trank you.






/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 5.x-197d9d3bf20a
Exec : simpleFoam
Date : May 17 2018
Time : 12:53:20
Host : "David-HP"
PID : 25925
I/O : uncollated
Case : /home/david/Escritorio/airFoil2D
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


SIMPLE: convergence criteria
field p tolerance 1e-05
field U tolerance 1e-05
field nuTilda tolerance 1e-05

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model SpalartAllmaras
Selecting patchDistMethod meshWave
RAS
{
RASModel SpalartAllmaras;
turbulence off;
printCoeffs on;
sigmaNut 0.66666;
kappa 0.41;
Cb1 0.1355;
Cb2 0.622;
Cw2 0.3;
Cw3 2;
Cv1 7.1;
Cs 0.3;
}

No MRF models present

No finite volume options present


Starting time loop

Time = 1

#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 ? at ??:?
#7 ? at ??:?
#8 ? at ??:?
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10 ? at ??:?
Excepción de coma flotante (`core' generado)


For ControlDict:


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

application simpleFoam;

startFrom latestTime;

startTime 0;

stopAt endTime;

endTime 500;

deltaT 1;

writeControl timeStep;

writeInterval 2;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable true;


// ************************************************** *********************** //


fvschemes;


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

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default cellLimited Gauss linear;
grad(U) cellMDLimited Gauss linear 1;
grad(p) faceMDLimited Gauss linear 1;
}

divSchemes
{
default none;
div(phi,U) Gauss linear;
div(phi,nuTilda) Gauss linear;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
div(U,p) upwind phi;
}

snGradSchemes
{
default limited;
}

wallDist
{
method meshWave;
}


// ************************************************** *********************** //


fvSolutions:


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

solvers
{
p
{
solver GAMG;
tolerance 1e-08;
relTol 0.1;
smoother GaussSeidel;
cacheAgglomeartion true;
nCellsInCoarsestLevel 100;
agglomerator faceAreaPair;
mergeLevels 1;
}

U
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-06;
relTol 0.1;
}

nuTilda
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-08;
relTol 0.1;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 10;
nCorrectors 10;
pRefCell 0;
pRefValue 0;

residualControl
{
p 1e-5;
U 1e-5;
nuTilda 1e-5;
}
}

relaxationFactors
{
fields
{
p 0.3;
}
equations
{
U 0.7;
nuTilda 0.7;
}
}


// ************************************************** *********************** //




I do checkMesh and is ok, so i don't know where is the problem. I try de solutions that are in this post but i can't simulate. Trank you.