 |
 |
 |
|
|
|
[Sponsors] | ||||
PimpleDyMFoam errors- Foam::error::printStack(Foam::Ostream&) at |
 |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 13, 2018, 19:57
|
PimpleDyMFoam errors- Foam::error::printStack(Foam::Ostream&) at
|
#1 |
|
Member
Join Date: Jan 2018
Location: Malaysia
Posts: 58
Rep Power: 8  |
Hi, i am doing a simulation of a small box moving within a rectangular bigger box to see it's air flow and pressure. i constructed my geo and mesh using blockMesh. I use k epsilon RAS model. The small box movement is solid body motion at Uy=20 m/s.
When running it in pimpeDyMFoam, the simulation stopped and errors were reported at the 5th time step. I tried adjusting the timeStep, BC, IC, fvSchemes and fvSolutions but the furthest it can go is this t=0.04s  Here are some attachment on my case. I will be very grateful if can get some suggestions and advice to solve this problem. Thanks a lot everyone. Code:
Create time
Create mesh for time = 0
Selecting dynamicFvMesh solidBodyMotionFvMesh
Selecting solid-body motion function linearMotion
Applying solid body motion to cellZone carMovingCZ
PIMPLE: no residual control data found. Calculations will employ 2 corrector loops
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 kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
}
No MRF models present
Reading/calculating face velocity Uf
No finite volume options present
Courant Number mean: 0.009848378 max: 0.01
Starting time loop
Courant Number mean: 0.009848378 max: 0.01
Time = 0.001
PIMPLE: iteration 1
GAMG: Solving for p, Initial residual = 1, Final residual = 0.003925482, No Iterations 5
time step continuity errors : sum local = 8.057605e-06, global = -3.441193e-06, cumulative = -3.441193e-06
GAMG: Solving for p, Initial residual = 0.00359735, Final residual = 0.00359735, No Iterations 0
time step continuity errors : sum local = 0.007444914, global = -3.441193e-06, cumulative = -6.882386e-06
PIMPLE: iteration 2
GAMG: Solving for p, Initial residual = 0.002297544, Final residual = 0.002297544, No Iterations 0
time step continuity errors : sum local = 0.004748193, global = -3.408226e-06, cumulative = -1.029061e-05
GAMG: Solving for p, Initial residual = 0.002297544, Final residual = 7.776222e-07, No Iterations 6
time step continuity errors : sum local = 1.607064e-06, global = -9.355428e-07, cumulative = -1.122615e-05
smoothSolver: Solving for epsilon, Initial residual = 0.01557478, Final residual = 1.028783e-10, No Iterations 2
smoothSolver: Solving for k, Initial residual = 1, Final residual = 4.677278e-08, No Iterations 5
ExecutionTime = 2.7 s ClockTime = 3 s
Courant Number mean: 0.04662531 max: 1.072731
Time = 0.002
PIMPLE: iteration 1
GAMG: Solving for p, Initial residual = 9.549607e-06, Final residual = 9.549607e-06, No Iterations 0
time step continuity errors : sum local = 0.0020843, global = -9.612726e-07, cumulative = -1.218743e-05
GAMG: Solving for p, Initial residual = 9.549607e-06, Final residual = 9.549607e-06, No Iterations 0
time step continuity errors : sum local = 0.0020843, global = -9.612726e-07, cumulative = -1.31487e-05
PIMPLE: iteration 2
GAMG: Solving for p, Initial residual = 9.548678e-06, Final residual = 9.548678e-06, No Iterations 0
time step continuity errors : sum local = 0.002084102, global = -9.610203e-07, cumulative = -1.410972e-05
GAMG: Solving for p, Initial residual = 9.548678e-06, Final residual = 4.039395e-07, No Iterations 3
time step continuity errors : sum local = 8.816414e-05, global = -5.810901e-07, cumulative = -1.469081e-05
smoothSolver: Solving for epsilon, Initial residual = 3.758072e-10, Final residual = 3.758072e-10, No Iterations 0
smoothSolver: Solving for k, Initial residual = 0.9701054, Final residual = 0.0002894628, No Iterations 1000
ExecutionTime = 8.02 s ClockTime = 8 s
Courant Number mean: 0.04708371 max: 3.388358
Time = 0.003
PIMPLE: iteration 1
GAMG: Solving for p, Initial residual = 0.7659974, Final residual = 0.008425877, No Iterations 3
time step continuity errors : sum local = 0.1417327, global = -1.503157e-07, cumulative = -1.484113e-05
GAMG: Solving for p, Initial residual = 0.4797737, Final residual = 0.005124784, No Iterations 2
time step continuity errors : sum local = 0.05755669, global = -6.746554e-08, cumulative = -1.490859e-05
PIMPLE: iteration 2
GAMG: Solving for p, Initial residual = 0.02166699, Final residual = 0.0005683549, No Iterations 1
time step continuity errors : sum local = 0.008305296, global = -5.039098e-08, cumulative = -1.495898e-05
GAMG: Solving for p, Initial residual = 0.0008800453, Final residual = 6.720886e-07, No Iterations 6
time step continuity errors : sum local = 9.477744e-06, global = 2.622759e-08, cumulative = -1.493276e-05
smoothSolver: Solving for epsilon, Initial residual = 3.418658e-07, Final residual = 3.418658e-07, No Iterations 0
smoothSolver: Solving for k, Initial residual = 0.7531732, Final residual = 9.384997e-07, No Iterations 82
ExecutionTime = 10.4 s ClockTime = 11 s
Courant Number mean: 0.05138682 max: 8.541179
Time = 0.004
PIMPLE: iteration 1
GAMG: Solving for p, Initial residual = 0.0117465, Final residual = 0.00123785, No Iterations 1
time step continuity errors : sum local = 0.007715017, global = 4.512535e-08, cumulative = -1.488763e-05
GAMG: Solving for p, Initial residual = 0.251486, Final residual = 0.007505026, No Iterations 1
time step continuity errors : sum local = 0.02218468, global = 7.138906e-08, cumulative = -1.481624e-05
PIMPLE: iteration 2
GAMG: Solving for p, Initial residual = 0.008602272, Final residual = 0.008602272, No Iterations 0
time step continuity errors : sum local = 0.02303742, global = 7.094134e-08, cumulative = -1.47453e-05
GAMG: Solving for p, Initial residual = 0.008602272, Final residual = 6.099818e-07, No Iterations 9
time step continuity errors : sum local = 1.633569e-06, global = 1.669112e-07, cumulative = -1.457839e-05
smoothSolver: Solving for epsilon, Initial residual = 4.358505e-09, Final residual = 4.358505e-09, No Iterations 0
smoothSolver: Solving for k, Initial residual = 0.2787908, Final residual = 8.078052e-07, No Iterations 29
ExecutionTime = 12.57 s ClockTime = 13 s
Courant Number mean: 0.04994512 max: 5.830545
Time = 0.005
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::surfaceInterpolation::makeWeights() const at ??:?
#4 Foam::surfaceInterpolation::weights() const at ??:?
#5 Foam::surfaceInterpolation::makeDeltaCoeffs() const at ??:?
#6 Foam::surfaceInterpolation::deltaCoeffs() const at ??:?
#7 Foam::fvPatch::deltaCoeffs() const at ??:?
#8 Foam::mixedFvPatchField<Foam::Vector<double> >::evaluate(Foam::UPstream::commsTypes) at ??:?
#9 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() in "/opt/openfoam4/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"
#10 Foam::solidBodyMotionFvMesh::update() at ??:?
#11 ? in "/opt/openfoam4/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? in "/opt/openfoam4/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"
Floating point exception (core dumped)
Code:
Time = 0
Mesh stats
points: 77950
faces: 213712
internal faces: 193880
cells: 67932
faces per cell: 6
boundary patches: 5
point zones: 0
face zones: 0
cell zones: 5
Overall number of cells of each type:
hexahedra: 67932
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: 2
The mesh has multiple regions which are not connected by any face.
<<Writing region information to "0/cellToRegion"
<<Writing region 0 with 2364 cells to cellSet region0
<<Writing region 1 with 65568 cells to cellSet region1
Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
car 2608 2610 ok (closed singly connected)
shaftTop 744 800 ok (non-closed singly connected)
shaftSide 11000 11110 ok (non-closed singly connected)
shaftBase 744 800 ok (non-closed singly connected)
oldInternalFaces 4736 4740 ok (closed singly connected)
Checking geometry...
Overall domain bounding box (0 0 0) (2.4 10 3.1)
Mesh has 0 geometric (non-empty/wedge) directions (0 0 0)
Mesh has 0 solution (non-empty) directions (0 0 0)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (1.2036e-17 5.109177e-17 1.36927e-16) OK.
Max cell openness = 1.734723e-16 OK.
Max aspect ratio = -1 OK.
Minimum face area = 0.01. Maximum face area = 0.01. Face area magnitudes OK.
Min volume = 0.001. Max volume = 0.001. Total volume = 67.932. Cell volumes OK.
Mesh non-orthogonality Max: 0 average: 0
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 1.065816e-13 OK.
Coupled point location match (average 0) OK.
Mesh OK.
Time = 0.001
Checking geometry...
Overall domain bounding box (0 0 0) (2.4 10 3.1)
Mesh has 0 geometric (non-empty/wedge) directions (0 0 0)
Mesh has 0 solution (non-empty) directions (0 0 0)
***Number of edges not aligned with or perpendicular to non-empty directions: 2746
<<Writing 5084 points on non-aligned edges to set nonAlignedEdges
Boundary openness (1.185061e-17 5.108442e-17 1.346771e-16) OK.
Max cell openness = 1.971277e-16 OK.
Max aspect ratio = -1 OK.
Minimum face area = 0.008. Maximum face area = 0.012. Face area magnitudes OK.
Min volume = 0.0008. Max volume = 0.0012. Total volume = 67.932. Cell volumes OK.
Mesh non-orthogonality Max: 11.30993 average: 1.494921
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.1818182 OK.
Coupled point location match (average 0) OK.
Failed 1 mesh checks.
Time = 0.002
Checking geometry...
Overall domain bounding box (0 0 0) (2.4 10 3.1)
Mesh has 0 geometric (non-empty/wedge) directions (0 0 0)
Mesh has 0 solution (non-empty) directions (0 0 0)
***Number of edges not aligned with or perpendicular to non-empty directions: 2746
<<Writing 5084 points on non-aligned edges to set nonAlignedEdges
Boundary openness (3.012562e-17 5.106319e-17 1.502784e-16) OK.
Max cell openness = 1.734723e-16 OK.
Max aspect ratio = -1 OK.
Minimum face area = 0.006. Maximum face area = 0.014. Face area magnitudes OK.
Min volume = 0.0006. Max volume = 0.0014. Total volume = 67.932. Cell volumes OK.
Mesh non-orthogonality Max: 21.80141 average: 2.90102
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.362069 OK.
Coupled point location match (average 0) OK.
Failed 1 mesh checks.
Time = 0.003
Checking geometry...
Overall domain bounding box (0 0 0) (2.4 10 3.1)
Mesh has 0 geometric (non-empty/wedge) directions (0 0 0)
Mesh has 0 solution (non-empty) directions (0 0 0)
***Number of edges not aligned with or perpendicular to non-empty directions: 2746
<<Writing 5084 points on non-aligned edges to set nonAlignedEdges
Boundary openness (2.959066e-17 5.103013e-17 1.497705e-16) OK.
Max cell openness = 1.734723e-16 OK.
Max aspect ratio = -1 OK.
Minimum face area = 0.004. Maximum face area = 0.016. Face area magnitudes OK.
Min volume = 0.0004. Max volume = 0.0016. Total volume = 67.932. Cell volumes OK.
Mesh non-orthogonality Max: 30.96376 average: 4.162129
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.75 OK.
Coupled point location match (average 0) OK.
Failed 1 mesh checks.
Time = 0.004
Checking geometry...
Overall domain bounding box (0 0 0) (2.4 10 3.1)
Mesh has 0 geometric (non-empty/wedge) directions (0 0 0)
Mesh has 0 solution (non-empty) directions (0 0 0)
***Number of edges not aligned with or perpendicular to non-empty directions: 2746
<<Writing 5084 points on non-aligned edges to set nonAlignedEdges
Boundary openness (2.957479e-17 5.09877e-17 1.429332e-16) OK.
Max cell openness = 1.734723e-16 OK.
Max aspect ratio = -1 OK.
Minimum face area = 0.002. Maximum face area = 0.018. Face area magnitudes OK.
Min volume = 0.0002. Max volume = 0.0018. Total volume = 67.932. Cell volumes OK.
Mesh non-orthogonality Max: 38.65981 average: 5.26124
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 1 OK.
Coupled point location match (average 0) OK.
Failed 1 mesh checks.
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: plus |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver GAMG;
tolerance 1e-02;
relTol 0;
smoother DICGaussSeidel;
cacheAgglomeration no;
maxIter 50;
}
pFinal
{
$p;
tolerance 1e-06;
relTol 0;
}
"pcorr.*"
{
$p
tolerance 1e-05;
relTol 0;
}
"(U|k|epsilon|nuTilda)"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-06;
relTol 0.1;
}
"(U|k|epsilon)Final"
{
$epsilon;
tolerance 1e-06;
relTol 0;
}
cellMotionUy
{
solver PCG;
preconditioner DIC;
tolerance 1e-08;
relTol 0;
}
}
PIMPLE
{
correctPhi no;
nOuterCorrectors 2;
nCorrectors 2;
nNonOrthogonalCorrectors 0;
}
residualControl
{
p 1e-06;
"(U|k|epsilon|nuTilda)" 1e-06;
}
relaxationFactors
{
"(U|k|epsilon).*" 1;
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: plus |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) cellLimited Gauss linear 1;
}
divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear limited corrected 0.5;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default limited corrected 0.5;
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: plus |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application pimpleDyMFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 0.01;
deltaT 0.001;
writeControl timeStep;
writeInterval 1;
purgeWrite 0;
writeFormat ascii;
writePrecision 7;
writeCompression off;
timeFormat general;
timePrecision 6;
runTimeModifiable true;
adjustTimeStep no;
maxCo 0.2;
functions
{
}
// ************************************************************************* //
Last edited by jiahui_93; February 13, 2018 at 20:58. Reason: attach files |
|
|
|
|
|
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| Building OpenFOAM1.7.0 from source | ata | OpenFOAM Installation | 46 | March 6, 2022 14:21 |
| How to write k and epsilon before the abnormal end | xiuying | OpenFOAM Running, Solving & CFD | 8 | August 27, 2013 16:33 |
| Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 | bookie56 | OpenFOAM Installation | 8 | August 13, 2011 05:03 |
| IcoFoam parallel woes | msrinath80 | OpenFOAM Running, Solving & CFD | 9 | July 22, 2007 03:58 |
| Could anybody help me see this error and give help | liugx212 | OpenFOAM Running, Solving & CFD | 3 | January 4, 2006 19:07 |
Linear Mode
