CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Problem with inconsistent patches (AMI, interDyMFoam) (https://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?

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 05:56.