|
[Sponsors] |
Problem with inconsistent patches (AMI, interDyMFoam) |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
January 30, 2013, 10:09 |
Problem with inconsistent patches (AMI, interDyMFoam)
|
#1 |
Member
Anon
Join Date: Oct 2012
Posts: 33
Rep Power: 14 |
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) 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: 76
Rep Power: 14 |
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: 14 |
Quote:
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: 14 |
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 | Search this Thread |
Display Modes | |
|
|
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 |