jrrygg |
January 30, 2013 09:09 |
Problem with inconsistent patches (AMI, interDyMFoam)
Hi!
I have been working on simulating a case with a rotating Pelton turbine using the InterDyMFoam-solver and the AMI-interface in OpenFOAM 2.1.1.
When I try to run the case I get the following error-message:
Code:
Create time
Create mesh for time = 0
Selecting dynamicFvMesh solidBodyMotionFvMesh
Selecting solid-body motion function rotatingMotion
Applying solid body motion to entire mesh
Reading field p_rgh
Reading field alpha1
Reading field U
Reading/calculating face flux field phi
AMI: Creating addressing and weights between 7602 source faces and 7602 target faces
AMI: Patch source weights min/max/average = 1, 1.009881271, 1.000225935
AMI: Patch target weights min/max/average = 1, 1.009881271, 1.000225935
Reading transportProperties
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
sigmaEps 1.3;
}
Reading g
Calculating field g.h
PIMPLE: Operating solver in PISO mode
time step continuity errors : sum local = 3.598673146e-08, global = -3.598673146e-08, cumulative = -3.598673146e-08
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 1.848725019e-06, No Iterations 6
time step continuity errors : sum local = 6.652957099e-14, global = -1.355912807e-14, cumulative = -3.598674502e-08
Courant Number mean: 9.852053092e-07 max: 0.006905813148
Starting time loop
Interface Courant Number mean: 0 max: 0
Courant Number mean: 9.832388315e-07 max: 0.00689202909
deltaT = 1.196172249e-07
Time = 1.196172249e-07
solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 1.196172249e-07 transformation: ((0 0 0) (1 (0 3.961438105e-06 0)))
AMI: Creating addressing and weights between 7602 source faces and 7602 target faces
AMI: Patch source weights min/max/average = 0.9999103424, 1.009878216, 1.000226318
AMI: Patch target weights min/max/average = 0.9999695281, 1.009919616, 1.000226192
Execution time for mesh.update() = 0.73 s
time step continuity errors : sum local = 7.958082654e-14, global = -1.621905271e-14, cumulative = -3.598676123e-08
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 2.428583894e-06, No Iterations 5
time step continuity errors : sum local = 3.700829363e-10, global = 1.683202693e-11, cumulative = -3.596992921e-08
MULES: Solving for alpha1
Phase-1 volume fraction = 1.076158237e-08 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 2.152316475e-08 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 3.228474712e-08 Min(alpha1) = -2.818921773e-24 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 4.304632949e-08 Min(alpha1) = -4.654056954e-25 Max(alpha1) = 1
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 0.007821959321, No Iterations 6
time step continuity errors : sum local = 3.381447662e-10, global = -1.188843696e-10, cumulative = -3.608881358e-08
GAMG: Solving for p_rgh, Initial residual = 0.04888923893, Final residual = 0.0003682615198, No Iterations 2
time step continuity errors : sum local = 1.903051747e-09, global = -5.90183726e-10, cumulative = -3.66789973e-08
GAMGPCG: Solving for p_rgh, Initial residual = 0.009021835592, Final residual = 2.580695706e-09, No Iterations 5
time step continuity errors : sum local = 1.44550723e-14, global = -6.293997608e-17, cumulative = -3.667899737e-08
--> FOAM FATAL ERROR:
Cannot find face 471634 in any of the patches
12
(
jetinlet
jetwalls
BBC
MBC
TBC
statopening
rotopen
rotsym
statinlet
statsym
AMI1
AMI2
)
It seems your patches are not consistent with the mesh : internalFaces:471604 total number of faces:542821
From function label polyBoundaryMesh::whichPatch(const label faceIndex) const
in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 541.
FOAM aborting
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::polyBoundaryMesh::whichPatch(int) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#3 void Foam::fvMatrix<double>::setValuesFromList<Foam::UList>(Foam::UList<int> const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#4 Foam::fixedInternalValueFvPatchField<double>::manipulateMatrix(Foam::fvMatrix<double>&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#5 Foam::incompressible::RASModels::kEpsilon::correct() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#6
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/interDyMFoam"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/interDyMFoam"
Aborted (core dumped)
I think the problem might be since I manually removed some boundaries and faces generated by blockMesh from the boundary-file, but even if I add up the faces I removed, I can not find the mentioned face 471634.
Any advice on how I can solve this? I will be happy to provide additional information if this is not sufficient.
The boundary file BEFORE I deleted outside, inside, top and bottom:
Code:
17
(
outside
{
type patch;
nFaces 0;
startFace 479206;
}
inside
{
type patch;
nFaces 20;
startFace 479206;
}
top
{
type patch;
nFaces 2;
startFace 479226;
}
bottom
{
type wall;
nFaces 9;
startFace 479228;
}
jetinlet
{
type wall;
nFaces 25;
startFace 479237;
}
jetwalls
{
type wall;
nFaces 494;
startFace 479262;
}
BBC
{
type wall;
nFaces 9792;
startFace 479756;
}
MBC
{
type wall;
nFaces 9547;
startFace 489548;
}
TBC
{
type wall;
nFaces 9462;
startFace 499095;
}
statopening
{
type wall;
nFaces 600;
startFace 508557;
}
rotif
{
type wall;
nFaces 0;
startFace 509157;
}
rotopen
{
type wall;
nFaces 5822;
startFace 509157;
}
rotsym
{
type symmetryPlane;
nFaces 6347;
startFace 514979;
}
statinlet
{
type wall;
nFaces 9727;
startFace 521326;
}
statsym
{
type symmetryPlane;
nFaces 4166;
startFace 531053;
}
AMI1
{
type cyclicAMI;
nFaces 0;
startFace 535219;
neighbourPatch AMI2;
transform noOrdering;
surface
{
}
}
AMI2
{
type cyclicAMI;
nFaces 0;
startFace 535219;
neighbourPatch AMI1;
transform noOrdering;
surface
{
}
}
)
|