[Yage]

Dear All,

I got a problem when i running a case after 13.7s and the error message are as follow.

I dont know where to address the problem since there is no such file: /processor0/system/data

Anyone know the reason of the error and how to fix it? Your idea will be appreciated

Regards,
Yage

sixDoFRigidBodyMotion constraints converged in 1 iterations
GAMG:  Solving for cellDisplacementx, Initial residual = 0.000596631121592, Final residual = 8.80835806891e-07, No Iterations 1
GAMG:  Solving for cellDisplacementy, Initial residual = 0.00253622052072, Final residual = 3.59717966112e-06, No Iterations 1
Execution time for mesh.update() = 0.400000000001 s
time step continuity errors : sum local = 4.7209126193e-05, global = -2.81590345296e-14, cumulative = 1.33428814917e-11
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 8.98927227239e-08, No Iterations 31
GAMG:  Solving for pcorr, Initial residual = 0.0149491885051, Final residual = 8.39185354527e-08, No Iterations 18
time step continuity errors : sum local = 4.71626629007e-05, global = 9.40722598644e-16, cumulative = 1.33438222143e-11
MULES: Solving for alpha1
Liquid phase volume fraction = 0.497798489543  Min(alpha1) = -1.8248782399e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.497792870625  Min(alpha1) = -1.52706361513e-19  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.497787251904  Min(alpha1) = -8.64703258312e-20  Max(alpha1) = 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.00392932727907, Final residual = 2.87233874218e-10, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.0108536992625, Final residual = 1.26227260124e-10, No Iterations 4
GAMG:  Solving for p_rgh, Initial residual = 0.0347651461109, Final residual = 6.94583913428e-08, No Iterations 13
GAMG:  Solving for p_rgh, Initial residual = 0.00045207711816, Final residual = 8.32607727393e-08, No Iterations 9
time step continuity errors : sum local = 1.05245433716e-11, global = 8.80409595103e-13, cumulative = 1.42242318094e-11
GAMG:  Solving for p_rgh, Initial residual = 0.000912468782066, Final residual = 6.67286654929e-08, No Iterations 7
GAMG:  Solving for p_rgh, Initial residual = 1.36334179433e-05, Final residual = 9.57980711943e-08, No Iterations 3
time step continuity errors : sum local = 1.21228692274e-11, global = -5.68176717021e-15, cumulative = 1.42185500423e-11
GAMG:  Solving for p_rgh, Initial residual = 3.58461867604e-05, Final residual = 8.09445165244e-08, No Iterations 4
GAMG:  Solving for p_rgh, Initial residual = 6.20619996787e-07, Final residual = 5.33814232481e-09, No Iterations 5
time step continuity errors : sum local = 6.75521153489e-13, global = -2.39156450586e-14, cumulative = 1.41946343972e-11
ExecutionTime = 9295.69 s  ClockTime = 9574 s

Interface Courant Number mean: 0.000264695414423 max: 0.125884296441
Courant Number mean: 0.0119652986721 max: 0.501405722438
deltaT = 0.00017060922284
Time = 13.7385806798

sixDoFRigidBodyMotion constraints converged in 1 iterations
GAMG:  Solving for cellDisplacementx, Initial residual = 0.000592454403998, Final residual = 8.73829434669e-07, No Iterations 1
GAMG:  Solving for cellDisplacementy, Initial residual = 0.00252591413159, Final residual = 3.59621763286e-06, No Iterations 1
Execution time for mesh.update() = 0.189999999999 s
time step continuity errors : sum local = 4.69013739192e-05, global = -2.37826724602e-14, cumulative = 1.41708517247e-11
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 4.00242832433e-08, No Iterations 33
GAMG:  Solving for pcorr, Initial residual = 0.0064818094357, Final residual = 8.86289596269e-08, No Iterations 19
time step continuity errors : sum local = 4.68144541874e-05, global = 8.26885557342e-16, cumulative = 1.41716786103e-11
MULES: Solving for alpha1
Liquid phase volume fraction = 0.497801047471  Min(alpha1) = -1.73699833221e-19  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.497795474835  Min(alpha1) = -2.43657115976e-19  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.497789902393  Min(alpha1) = -3.75155380006e-20  Max(alpha1) = 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.00392628959014, Final residual = 2.72994517484e-10, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.0108589761566, Final residual = 1.33756418069e-10, No Iterations 4
GAMG:  Solving for p_rgh, Initial residual = 0.0347461863115, Final residual = 1.8200314363e+36, No Iterations 1000
GAMG:  Solving for p_rgh, Initial residual = 0.905876455385, Final residual = 6.4145861131e+39, No Iterations 1000
time step continuity errors : sum local = 1.61815689411e+72, global = 6.56972143739e+67, cumulative = 6.56972143739e+67
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 8.73287810527e-08, No Iterations 21
GAMG:  Solving for p_rgh, Initial residual = 0.0170331773648, Final residual = 6.86355039561e-08, No Iterations 17
time step continuity errors : sum local = 1.29814723091e+124, global = -6.41295148562e+122, cumulative = -6.41295148562e+122
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = nan, No Iterations 1000
GAMG:  Solving for p_rgh, Initial residual = nan, Final residual = nan, No Iterations 1000
[3] [0] [1] 
[1] 
[1] --> FOAM FATAL IO ERROR: 
[1] wrong token type - expected Scalar, found on line 7 the word 'nan'
[1] 
[1] file: /home/yage/waveDyMFoam/waveDyMFlume4/processor1/system/data::solverPerformance::p_rgh at line 7[2] 

[0] 
[0] --> FOAM FATAL IO ERROR: 
[0] wrong token type - expected Scalar, found on line 7 the word 'nan'
[0] 
[0] file: /home/yage/waveDyMFoam/waveDyMFlume4/processor0/system/data::solverPerformance::p_rgh at line 7.
[0] 
[0]     From function operator>>(Istream&, Scalar&)
[0]     in file lnInclude/Scalar.C at line 91.
[0] 
FOAM parallel run exiting
[0] 
.
[1] 
[1]     From function operator>>(Istream&, Scalar&)
[1]     in file lnInclude/Scalar.C at line 91.
[1] 
FOAM parallel run exiting
[1] 
[2] 
[2] --> FOAM FATAL IO ERROR: 
[2] wrong token type - expected Scalar, found on line 7 the word 'nan'
[2] 
[2] file: /home/yage/waveDyMFoam/waveDyMFlume4/processor2/system/data::solverPerformance::p_rgh at line 7.
[2] 
[2]     From function operator>>(Istream&, Scalar&)
[2]     in file lnInclude/Scalar.C at line 91.
[2] 
FOAM parallel run exiting
[2] 

[3] 
[3] --> FOAM FATAL IO ERROR: 
[3] wrong token type - expected Scalar, found on line 7 the word 'nan'
[3] 
[3] file: /home/yage/waveDyMFoam/waveDyMFlume4/processor3/system/data::solverPerformance::p_rgh at line 7.
[3] 
[3]     From function operator>>(Istream&, Scalar&)
[3]     in file lnInclude/Scalar.C at line 91.
[3] 
FOAM parallel run exiting
[3] 
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD 
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.

[alexeym]

Hi,

your solution is diverging:

time step continuity errors : sum local = 1.61815689411e+72, global = 6.56972143739e+67, cumulative = 6.56972143739e+67
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 8.73287810527e-08, No Iterations 21
GAMG:  Solving for p_rgh, Initial residual = 0.0170331773648, Final residual = 6.86355039561e-08, No Iterations 17
time step continuity errors : sum local = 1.29814723091e+124, global = -6.41295148562e+122, cumulative = -6.41295148562e+122
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = nan, No Iterations 1000
GAMG:  Solving for p_rgh, Initial residual = nan, Final residual = nan, No Iterations 1000

Check boundary conditions. Check your mesh (if sixDoFRigidBodyMotion distorts mesh, the problem can be here).

[Yage]

here is the CheckMesh

Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 151992
internal points: 0
faces: 301932
internal faces: 149940
cells: 75312
boundary patches: 7
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 75312
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
inlet 336 674 ok (non-closed singly connected)
bottom 228 458 ok (non-closed singly connected)
body 240 480 ok (non-closed singly connected)
outlet 336 674 ok (non-closed singly connected)
atmosphere 228 458 ok (non-closed singly connected)
frontBack 0 0 ok (empty)
defaultFaces 150624 151992 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-4.2 -2.1 0) (4.2 2.1 0.1)
Mesh (non-empty, non-wedge) directions (1 1 0)
Mesh (non-empty) directions (1 1 0)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (-1.19175836492e-19 1.72804962913e-19 -1.47987190843e-15) OK.
Max cell openness = 1.66428773226e-16 OK.
Max aspect ratio = 10.0202965596 OK.
Minumum face area = 0.000152083298349. Maximum face area = 0.0083164539. Face area magnitudes OK.
Min volume = 1.52083298349e-05. Max volume = 0.000346518912777. Total volume = 3.503. Cell volumes OK.
Mesh non-orthogonality Max: 0 average: 0
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 1.6438385609e-10 OK.
Coupled point location match (average 0) OK.

Mesh OK.

Time = 0.1


Checking geometry...
Overall domain bounding box (-4.2 -2.1 -1.38777878078e-16) (4.2 2.1 0.1)
Mesh (non-empty, non-wedge) directions (1 1 0)
Mesh (non-empty) directions (1 1 0)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (-1.87770967813e-18 1.71017541375e-18 5.24397515732e-16) OK.
Max cell openness = 2.5595914843e-16 OK.
Max aspect ratio = 10.0205754735 OK.
Minumum face area = 0.000152078920886. Maximum face area = 0.00831724976309. Face area magnitudes OK.
Min volume = 1.52078920886e-05. Max volume = 0.000346552070955. Total volume = 3.503. Cell volumes OK.
Mesh non-orthogonality Max: 0.00276820289483 average: 0.00142697429581
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.000115423565825 OK.
Coupled point location match (average 0) OK.

Mesh OK.

[alexeym]

OK, your initial mesh non-orthogonality is 0, but at 0.1 seconds it's already 0.00276820289483, do you know the value at 13.7 s? Also mesh is just one of possible reasons for solution divergence, what are your boundary conditions?

[Yage]

Hi Alexeym,

Thanks for your reply, is there any experience for the boundary condition for a floating object on surface wave?

Regards,
Yage

[Yage]

the boundary is as followed:
alpha1:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5-dev |
| \\ / A nd | Revision: 1740 |
| \\/ M anipulation | Web: http://www.OpenFOAM.org |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object alpha.org;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type waveAlpha;
refValue uniform 0;
refGrad uniform 0;
valueFraction uniform 1;
value uniform 0;
}
bottom
{
type zeroGradient;
}
body
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
atmosphere
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
frontBack
{
type empty;
}
}

p_rgh:

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type zeroGradient;
}
bottom
{
type zeroGradient;
}
body
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
atmosphere
{
type totalPressure;
U U;
phi phi;
rho none;
psi none;
gamma 1;
p0 uniform 0;
value uniform 0;
}
frontBack
{
type empty;
}
}

U:

FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U.org;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform ( 0.0 0.0 0.0);

boundaryField
{
inlet
{
type waveVelocity;
refValue uniform ( 0 0 0 );
refGradient uniform ( 0 0 0 );
valueFraction uniform 1;
value uniform ( 0 0 0 );
}
bottom
{
type fixedValue;
value uniform ( 0 0 0 );
}
body
{
type fixedValue;
value uniform ( 0 0 0 );
}
outlet
{
type fixedValue;
value uniform ( 0 0 0 );
}
atmosphere
{
type pressureInletOutletVelocity;
value uniform ( 0 0 0 );
}
frontBack
{
type empty;
}
}

[Yage]

for pointdisplacement:

FoamFile
{
version 2.0;
format ascii;
class pointVectorField;
location "0.01";
object pointDisplacement;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 0 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
inlet
{
type fixedValue;
value uniform (0 0 0);
}
bottom
{
type fixedValue;
value uniform (0 0 0);
}
outlet
{
type fixedValue;
value uniform (0 0 0);
}
atmosphere
{
type fixedValue;
value uniform (0 0 0);
}

body
{
// type sixDoFRigidBodyTestDisplacement;

type sixDoFRigidBodyDisplacement;
mass 15;
centreOfMass (0.0 0 0.05);
momentOfInertia (1e10 1e10 4.05);

rhoInf 1; // needed only for solvers solving for kinematic pressure
velocity (0 0 0);
acceleration (0 0 0);
angularMomentum (0 0 0);
torque (0 0 0);
rhoName rhoInf;
rhoInf 1;

constraints
{
maxIterations 1000;
fixedPlane1
{
sixDoFRigidBodyMotionConstraint fixedPlane;
tolerance 1e-06;
relaxationFactor 0.7;
fixedPlaneCoeffs
{
refPoint ( -2.0 -0.5 0.1 )( -2.1 -0.5 0.1 )( -2.0 -0.6 0.1 );
normal (0 0 1);
}
}
}


value uniform (0 0 0);
}

frontBack
{
type empty;
}
}

and mesh at 13.7

Time = 13.7


Checking geometry...
Overall domain bounding box (-4.2 -2.10000000001 -1.98919783423e-10) (4.20000000013 2.1 0.100000000167)
Mesh (non-empty, non-wedge) directions (1 1 0)
Mesh (non-empty) directions (1 1 0)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (2.4257806043e-18 -1.38551600356e-17 -3.25594344088e-16) OK.
Max cell openness = 2.28374555038e-16 OK.
Max aspect ratio = 10.66552642 OK.
Minumum face area = 0.000114842446218. Maximum face area = 0.00899173923168. Face area magnitudes OK.
Min volume = 1.14842446219e-05. Max volume = 0.00036591130215. Total volume = 3.503. Cell volumes OK.
Mesh non-orthogonality Max: 72.6204018364 average: 3.54511388904
*Number of severely non-orthogonal faces: 2.
Non-orthogonality check OK.
<<Writing 2 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 0.961251004341 OK.
Coupled point location match (average 0) OK.

Mesh OK.

[Yage]

And in addition,
the object is allowed to do heave, surge and pitch motion.

[alexeym]

Hi,

well, now you can also post your fvSchemes, as at 13.7 you've got rather non-orthogonal mesh.

Also I can suppose that at 13.7 you've got whole domain filled with water and you've got problems at the atmosphere patch.

In fact there are lots of possibilities, it'll be easier for everybody if you just describe your problem, attach archive with your case files, it'll allow people not to post questions here but look for the answer in the case files.

[Yage]

Hi alexeym,

OK, the problem is a 2D case and that is a object floating on the free surface. the wave is generated by waveDyMFoam package. the scheme and solution file is attached below

[Yage]

here are the fvschemes and fvsolution

FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{

default Euler;
}

gradSchemes
{

default Gauss linear;
grad(U) Gauss linear;
grad(alpha1) Gauss linear;
}

divSchemes
{

div(rho*phi,U) Gauss limitedLinearV 1;

div((muEff*dev(T(grad(U))))) Gauss linear;
div(phi,alpha) Gauss MUSCL;
div(phirb,alpha) Gauss interfaceCompression;
}

laplacianSchemes
{
default Gauss linear corrected;

}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p_rgh;
pcorr;
gamma;
}

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

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

solvers
{
cellDisplacement
{
solver GAMG;
tolerance 1e-5;
relTol 0;
smoother GaussSeidel;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

pcorr GAMG
{
tolerance 1e-7;
relTol 0.0;

smoother DIC;//GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;

cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
};

p_rgh GAMG
{
tolerance 1e-7;
relTol 0.0;

smoother DIC;//GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;

cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
};

p_rghFinal GAMG
{
tolerance 1e-8;
relTol 0.0;

smoother DIC;//GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;

cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
};


U PBiCG
{
preconditioner DILU;
tolerance 1e-09;
relTol 0;
};

UFinal PBiCG
{
preconditioner DILU;
tolerance 1e-09;
relTol 0;
};


gamma PBiCG
{
preconditioner DILU;
tolerance 1e-07;
relTol 0;
};
}


PIMPLE
{
pdRefCell 0;
pdRefValue 0;
momentumPredictor yes;
nOuterCorrectors 1;
nCorrectors 3;
nNonOrthogonalCorrectors 1;
nAlphaCorr 1;
nAlphaSubCycles 3;
cAlpha 1;
}

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

[Yage]

Screenshot from 2014-10-24 11:16:55.jpg

and this is the screenshot for t=13.7

[muth]

Mr.YageEven i face the same error. Did you solved this error?

[foamiste]

Mr. Muth did you solved this error? because I have the same problem even everything seems OK