CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Problem with inconsistent patches (AMI, interDyMFoam) (http://www.cfd-online.com/Forums/openfoam-solving/112522-problem-inconsistent-patches-ami-interdymfoam.html)

jrrygg January 30, 2013 10: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
    {
    }
}

)


maHein January 31, 2013 11:25

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 January 31, 2013 13:39

Quote:

Originally Posted by maHein (Post 405326)
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 February 1, 2013 04:23

Here are some pictures of my mesh:

Mesh consisting of a stationary (outer) domain with jet inlet, and a rotating (inner) domain with buckets
https://dl.dropbox.com/u/2820596/donut-mesh6.png

AMI1 / AMI2 patch
https://dl.dropbox.com/u/2820596/donut-mesh6-AMI1.png


All times are GMT -4. The time now is 10:26.