GGI and Topological changes - checkMesh fails
Hi foamers,
I am working to combine the GGI interface and topological changes of the mesh in OpenFOAM 1.6 ext. In the video you can see my test case, a cylinder where the right (or front) part is growing (translation) in the negative z-direction whereas the left (or back) part is rotating. https://dl.dropbox.com/u/86131187/GgiTopoChange.avi Calculation was done with the transientSimpleDyMFoam solver on a single processor, at first sight the results look ok. Post processing works fine with with the paraFoam utility. But running checkMesh gives the following error Code:
--> FOAM FATAL ERROR: Code:
/*---------------------------------------------------------------------------*\ Code:
... faceZone at time step 0 Code:
... Someone out there who can shed some light on this issue? Best regards, Klaus |
is it a bug...
some update to my post above
running checkMesh utility for a single time step works perfect Code:
checkMesh -time 0 Code:
checkMesh -time 0.005 Today i tried to find out what's going wrong. In the file polyMeshIO.C the line Code:
boundary_.calcGeometry(); Code:
... Code:
... Question to the experts: is it allowed to do that or are there any side effects which i should be aware of? Comments would be fine. Best regards, Klaus |
Hmm... looks like the zones were written to disk with addressing prior to topo-changes. This may not be an actual problem, since you say that the solution looks fine, but an issue where inconsistent addressing is written to disk. Still needs to be fixed though...
|
Quote:
t=0; Code:
... screw Code:
... screw The problem arises when the mesh is read from disk. Best, Klaus |
Dear foamer
Did you involve with this error? Quote:
Quote: // Create neighbouring face centres using interpolation if (master()) { const label shadowID = shadowIndex(); // Get the transformed and interpolated shadow face cell centers reconFaceCellCentresPtr_ = new vectorField ( interpolate ( boundaryMesh()[shadowID].faceCellCentres() - boundaryMesh()[shadowID].faceCentres() ) + faceCentres() ); } else { FatalErrorIn("void ggiPolyPatch::calcReconFaceCellCentres() const") << "Attempting to create reconFaceCellCentres on a shadow" << abort(FatalError); } Do you know, why this error appears? when I use cyclicGgi for top and bottom of rectangular domain, the problem appear when running simulation. |
All times are GMT -4. The time now is 15:14. |