CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPO

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By alexeym

Reply
 
LinkBack Thread Tools Display Modes
Old   February 3, 2014, 07:07
Default Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPO
  #1
Member
 
adarsh tiwari's Avatar
 
adarsh tiwari
Join Date: Feb 2014
Location: Bangalore
Posts: 40
Rep Power: 3
adarsh tiwari is on a distinguished road
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
adarsh tiwari is offline   Reply With Quote

Old   February 3, 2014, 08:45
Default
  #2
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,133
Rep Power: 20
alexeym will become famous soon enoughalexeym will become famous soon enough
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).
alexeym is offline   Reply With Quote

Old   February 3, 2014, 12:58
Unhappy
  #3
Member
 
adarsh tiwari's Avatar
 
adarsh tiwari
Join Date: Feb 2014
Location: Bangalore
Posts: 40
Rep Power: 3
adarsh tiwari is on a distinguished road
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
adarsh tiwari is offline   Reply With Quote

Old   February 3, 2014, 16:39
Default
  #4
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,133
Rep Power: 20
alexeym will become famous soon enoughalexeym will become famous soon enough
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
alexeym is offline   Reply With Quote

Old   February 4, 2014, 00:14
Default
  #5
Member
 
adarsh tiwari's Avatar
 
adarsh tiwari
Join Date: Feb 2014
Location: Bangalore
Posts: 40
Rep Power: 3
adarsh tiwari is on a distinguished road
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
Attached Images
File Type: jpg controlDict.jpg (17.8 KB, 16 views)
File Type: jpg fvSchemes.jpg (21.3 KB, 11 views)
File Type: jpg mut.jpg (18.1 KB, 7 views)
File Type: jpg p.jpg (18.7 KB, 9 views)
Attached Files
File Type: txt fvSolution.txt (1.9 KB, 6 views)
adarsh tiwari is offline   Reply With Quote

Old   February 4, 2014, 03:25
Default
  #6
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,133
Rep Power: 20
alexeym will become famous soon enoughalexeym will become famous soon enough
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 likes this.
alexeym is offline   Reply With Quote

Old   February 5, 2014, 07:23
Default
  #7
Member
 
adarsh tiwari's Avatar
 
adarsh tiwari
Join Date: Feb 2014
Location: Bangalore
Posts: 40
Rep Power: 3
adarsh tiwari is on a distinguished road
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
adarsh tiwari is offline   Reply With Quote

Old   February 5, 2014, 09:06
Default
  #8
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,133
Rep Power: 20
alexeym will become famous soon enoughalexeym will become famous soon enough
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

Code:
cellLimited
but

Code:
cellLimited Gauss linear 1;
Also you can just attach your files instead of posting them in the body of a message.
alexeym is offline   Reply With Quote

Old   February 7, 2014, 05:52
Default
  #9
Member
 
adarsh tiwari's Avatar
 
adarsh tiwari
Join Date: Feb 2014
Location: Bangalore
Posts: 40
Rep Power: 3
adarsh tiwari is on a distinguished road
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
adarsh tiwari is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 06:32.