[jrrygg]

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:

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:

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
    {
    }
}

)

[maHein]

I think there are some serious problems with your mesh, since none of your AMI patches contains any faces. Have you checked your mesh with paraview? Maybe you could post some images?

[jrrygg]

Originally Posted by maHein

I think there are some serious problems with your mesh, since none of your AMI patches contains any faces. Have you checked your mesh with paraview? Maybe you could post some images?

Thank you for your reply!

I also wondered about why it had 0 faces, but it seems to be the same in the propeller-tutorial which I have based my case on. When I look at the AMI-patches in paraView, however, it looks like they have faces. I will post some images tomorrow morning so it will be clearer.

[jrrygg]

Here are some pictures of my mesh:

Mesh consisting of a stationary (outer) domain with jet inlet, and a rotating (inner) domain with buckets


AMI1 / AMI2 patch