CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Problem with inconsistent patches (AMI, interDyMFoam)

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 30, 2013, 10:09
Default Problem with inconsistent patches (AMI, interDyMFoam)
  #1
Member
 
Anon
Join Date: Oct 2012
Posts: 33
Rep Power: 14
jrrygg is on a distinguished road
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
    {
    }
}

)
jrrygg is offline   Reply With Quote

Old   January 31, 2013, 11:25
Default
  #2
Member
 
Join Date: Jun 2012
Posts: 76
Rep Power: 14
maHein is on a distinguished road
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?
maHein is offline   Reply With Quote

Old   January 31, 2013, 13:39
Default
  #3
Member
 
Anon
Join Date: Oct 2012
Posts: 33
Rep Power: 14
jrrygg is on a distinguished road
Quote:
Originally Posted by maHein View Post
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 is offline   Reply With Quote

Old   February 1, 2013, 04:23
Default
  #4
Member
 
Anon
Join Date: Oct 2012
Posts: 33
Rep Power: 14
jrrygg is on a distinguished road
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
jrrygg is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with FloatingObject Leech OpenFOAM Running, Solving & CFD 10 March 29, 2012 16:24
openFoam 2.1.x: creating AMI patches from a cylinder. ebah6 OpenFOAM Running, Solving & CFD 5 March 26, 2012 19:06
problem occurs in interDyMFoam solver jignesh_thaker2007 OpenFOAM 2 March 18, 2012 22:11
Lagrangain tracking with sliding mesh with AMI patches jiejie OpenFOAM Running, Solving & CFD 0 March 7, 2012 00:59
error using interDyMFoam with kOmegaSST to simulate sloshing anmartin OpenFOAM Running, Solving & CFD 0 July 20, 2010 14:21


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