Hi Thomas,
thanks. the result looks nice. thats exactly what i´ve want. i hope i can transfer the setting on a mesh made by shm (different number of faces on the opposite sides)... what steps did you take? 1) blockMesh 2) faceSet 3) setsToZones -noFlipMap 4) change boundary-file best regards |
Quote:
Hello FabOr, for the second step, I used setSet -batch system/setSetBatch (just for the sake of completeness, faceSet should also work (I think ;)) ), all other steps are ok. Best regards, -Thomas |
Hi Thomas,
it works also with shm-meshes:) Thank you very much! Best regards |
Hi Thomas,
Thank you very much for the case file and it helped me to understand ggi well. Can you suggest me an idea of using the non-matching grids(think-cyclicGgi) and a rotating solid(ggi) together in the same case? |
further explanation needed
Quote:
sorry, I don't understand what you mean, can you explain it furhter and/or supply a sketch? Best regards, -Thomas |
Quote:
Dear Thomas, I am using OF 1.5-dev (the last version released). I am simulating a blade of a radial pump. I am using cyclicGgi patch, with the bridgeOverlap=false option. But when I run the MRFSimpleFoam solver, the message "Found uncovered faces for GGI interface Per1/Per2 while the bridgeOverlap option is not set in the boundary file." appears. I personally think that i would have to set this option true, because the cyclic patches do not overlap with themselves. What can You kindly suggest me to this purpose? Best regards, Claudio However, when I run |
Hi Claudio,
meanwhile there is OF-1.6-ext which is also capable of GGI interpolation and is the successor of 1.5-dev, so the things stated above are a little bit outdated. Regarding your problem, if the patches do not geometrically match, then set bridgeOverlap to true. However, for a radial pump, I'm not sure why they shouldn't do so. If you created your mesh with sHM, than there might be problems because sHM edges are not following strictly the stl geometry, and if this was the case, you'd have to be clear about the effects of setting bridgeOverlap to true to get the simulation running. Regards, -Thomas |
Quote:
Does rotationAngle mean the angle between the two patches which are setted as cyclicGgi? How to decide the angle is positive or negative? Because in the 1.6-ext tutorial, the case "turboPassangeRotating", the rotationAngle is 30 degree for one patch and -30 degree for the couple patch..I don't know why... Thank you so much! Aqua |
I see it's quite a past post, but still felt relevant for my problem.
I am trying to run overlapGgi and getting the error as: --> FOAM FATAL ERROR: Found uncovered faces for GGI interface innercylinder/outercylinder This is an unrecoverable error. Aborting. boundary blade_1 { type wall; nFaces 46561; startFace 5183059; } blade_2 { type wall; nFaces 39728; startFace 5229620; } innercylinder { type overlapGgi; nFaces 41878; startFace 5269348; rotationAxis (0 0 1); nCopies 30; shadowPatch outercylinder; zone innercylinder_faces; } outercylinder { type overlapGgi; nFaces 41880; startFace 5326982; rotationAxis (0 0 1); nCopies 30; shadowPatch innercylinder; zone outercylinder_faces; } farfield { type patch; nFaces 15756; startFace 5311226; } setBatch faceSet innercylinder_faces new patchToFace innercylinder faceSet outercylinder_faces new patchToFace outercylinder cellSet rotor new cellToCell cellRegion0 quit dynamicMeshDict dynamicFvMeshLib "libtopoChangerFvMesh.so"; dynamicFvMesh turboFvMesh;//mixerGgiFvMesh; turboFvMeshCoeffs//mixerGgiFvMeshCoeffs { coordinateSystem { type cylindrical; origin (0 0 0); axis (0 0 1); direction (1 0 0); } rpm { cellRegion0 360; cellRegion1 0.0; } slider { insidecylinder_faces 360; } } MRFZones cellRegion0//rotor { nonRotatingPatches (farfield); origin origin [0 1 0 0 0 0 0] (0 0 0); axis axis [0 0 0 0 0 0 0] (0 0 1); omega omega [0 0 -1 0 0 0 0] 360; } Command regionCellSets setSet -batch setBatch setsToZones -noFlipMap solverName Can any one tell me how to resolve this error. The figure shows blade1, blade2 and innercylinder. Thanks a lot. |
All times are GMT -4. The time now is 10:25. |