# Problem with inconsistent patches (AMI, interDyMFoam)

 Register Blogs Members List Search Today's Posts Mark Forums Read

 January 30, 2013, 10:09 Problem with inconsistent patches (AMI, interDyMFoam) #1 Member   Anon Join Date: Oct 2012 Posts: 33 Rep Power: 7 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::setValuesFromList(Foam::UList const&, Foam::UList const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so" #4 Foam::fixedInternalValueFvPatchField::manipulateMatrix(Foam::fvMatrix&) 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 { } } )```

 January 31, 2013, 11:25 #2 Member   Join Date: Jun 2012 Posts: 73 Rep Power: 7 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?

January 31, 2013, 13:39
#3
Member

Anon
Join Date: Oct 2012
Posts: 33
Rep Power: 7
Quote:
 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?

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.

 February 1, 2013, 04:23 #4 Member   Anon Join Date: Oct 2012 Posts: 33 Rep Power: 7 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

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Leech OpenFOAM Running, Solving & CFD 10 March 29, 2012 15:24 ebah6 OpenFOAM Running, Solving & CFD 5 March 26, 2012 18:06 jignesh_thaker2007 OpenFOAM 2 March 18, 2012 22:11 jiejie OpenFOAM Running, Solving & CFD 0 March 7, 2012 00:59 anmartin OpenFOAM Running, Solving & CFD 0 July 20, 2010 13:21

All times are GMT -4. The time now is 04:50.