October 11, 2021, 16:54
|
How to solve the error message "Foam::error::printStack(Foam::Ostream&)"
|
#1
|
|
New Member
F
Join Date: Oct 2021
Posts: 10
Rep Power: 5
|
Good day,
I am new to OpenFoam, so probably my mistakes are going to be obvious to some of you.
I am trying to simulate an external flow of water around a cylinder (2d).
I copied the tutorial Airfoil2d (simpleFoam) folder. I then imported a Fluent Mesh.
I changed just a couple of parameters here and there, but when I run the solver I get the below error message
Could anybody give me some help please?
Thanks a lot
- Control Dict
Code:
application simpleFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 500;
deltaT 1;
writeControl timeStep;
writeInterval 50;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression off;
timeFormat general;
timePrecision 6;
runTimeModifiable true;
functions
{
momErr
{
type momentumError;
libs (fieldFunctionObjects);
executeControl writeTime;
writeControl writeTime;
}
contErr
{
type div;
libs (fieldFunctionObjects);
field phi;
executeControl writeTime;
writeControl writeTime;
}
turbulenceFields1
{
type turbulenceFields;
libs (fieldFunctionObjects);
fields
(
k
epsilon
nut
nuEff
R
devReff
L
I
);
executeControl writeTime;
writeControl writeTime;
}
mag1
{
type mag;
libs (fieldFunctionObjects);
field turbulenceProperties:R;
result magR;
executeControl writeTime;
writeControl writeTime;
}
}
- Turbolence Properties
Code:
simulationType RAS;
RAS
{
RASModel SpalartAllmaras;
turbulence on;
printCoeffs on;
}
- Transport Properties
Code:
transportModel Newtonian;
nu 1e-06;
- U
Code:
internalField uniform (1.5 0.01 0);
boundaryField
{
inlet
{
type freestreamVelocity;
freestreamValue $internalField;
}
outlet
{
type freestreamVelocity;
freestreamValue $internalField;
}
walls
{
type noSlip;
}
frontandback
{
type empty;
}
}
- p
i
Code:
nternalField uniform 0;
boundaryField
{
inlet
{
type freestreamPressure;
freestreamValue $internalField;
}
outlet
{
type freestreamPressure;
freestreamValue $internalField;
}
walls
{
type zeroGradient;
}
frontandback
{
type empty;
}
}
- nut
Code:
internalField uniform 0;
internalField uniform 0.14;
boundaryField
{
inlet
{
type freestream;
freestreamValue uniform 0.14;
}
outlet
{
type freestream;
freestreamValue uniform 0.14;
}
walls
{
type nutUSpaldingWallFunction;
value uniform 0;
}
frontandback
{
type empty;
}
}
- nuTIlda
Code:
internalField uniform 0.14;
boundaryField
{
inlet
{
type freestream;
freestreamValue uniform 0.14;
}
outlet
{
type freestream;
freestreamValue uniform 0.14;
}
walls
{
type fixedValue;
value uniform 0;
}
frontandback
{
type empty;
}
}
- Boundary of the polymesh
Code:
walls
{
type wall;
physicalType wall;
nFaces 28;
startFace 43092;
}
inlet
{
type patch;
physicalType inlet;
nFaces 483;
startFace 43120;
}
outlet
{
type patch;
physicalType outlet;
nFaces 117;
startFace 43603;
}
frontandback
{
type empty;
physicalType empty;
nFaces 43406;
startFace 43720;
}
- fvSchemes
Code:
ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default none;
div(phi,U) bounded Gauss linearUpwind grad(U);
div(phi,nuTilda) bounded Gauss linearUpwind grad(nuTilda);
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
wallDist
{
method meshWave;
}
- fvSolution
Code:
solvers
{
p
{
solver GAMG;
tolerance 1e-06;
relTol 0.1;
smoother GaussSeidel;
}
U
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-08;
relTol 0.1;
}
nuTilda
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-08;
relTol 0.1;
}
}
SIMPLE
{
nNonOrthogonalCorrectors 0;
residualControl
{
p 1e-5;
U 1e-5;
nuTilda 1e-5;
}
}
relaxationFactors
{
fields
{
p 0.8;
}
equations
{
U 0.7;
nuTilda 0.7;
}
}
- Error Message
Code:
smoothSolver: Solving for Ux, Initial residual = 0.504513, Final residual = 0.0458849, No Iterations 52
smoothSolver: Solving for Uy, Initial residual = 0.506237, Final residual = 0.0415468, No Iterations 52
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib/x86_64-linux-gnu/libpthread.so.0
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::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 ??:?
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7 ? at ??:?
#8 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
#9 ? in /usr/lib/openfoam/openfoam2106/platforms/linux64GccDPInt32Opt/bin/simpleFoam
#10 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#11 ? in /usr/lib/openfoam/openfoam2106/platforms/linux64GccDPInt32Opt/bin/simpleFoam
Floating point exception (core dumped)
|
|
|
Similar Threads
Linear Mode
