Hello, i want to make tests on benchmark Flow past the cylinder. In application Gmsh I generated .msh file for the problem, I generated also some refinement, but the grid is quite coarse. I used GridIn:readMsh, but the fem program throws following exception
----------------------------------------------------
Exception on processing:
--------------------------------------------------------
An error occurred in line <642> of file
</home/marek/deal.II/source/grid/grid_reordering.cc> in function
static void dealii::GridReordering<dim,
spacedim>::invert_all_cells_of_negative_grid(const
std::vector<dealii::Point<spacedim>, std::allocator<dealii::Point<spacedim> >
>&, std::vector<dealii::CellData<dim>, std::allocator<dealii::CellData<dim> >
>&) [with int dim = 2, int spacedim = 2]
The violated condition was:
n_negative_cells==0 || n_negative_cells==cells.size()
The name and call sequence of the exception was:
ExcInternalError()
Additional Information:
(none)--og!
----------------------------------------------------
I have thought, that Deal.II fully supports .msh format. Could anyone give me
an advice how to preprocess grids in gmsh in order to have correct numbering.
Alternatively, anyone having mesh in deal.ii compatible format for Flow past
the cylinder problem will be welcomed :-)
Could You give me an hint on generating non-standard grids (=not from
GridGenerator), which would be fully compatible with deal's GridIn:read_xxx?
Thank You in advance
Marek C
cylinder(3).msh
Description: Mesh model
lc = 0.02;
r = 0.1;
d = 2*r;
h = 2*d;
assym = 0.1*d;
depth = 0.5;
l1 = 2*d;
l2 = 20*d;
n = 1;
nl1 = n*8;
nl2 = n*88;
nc = n*8;
//nz = 16;
nz = 1;
Point(1) = {0,0,0,lc};
Point(2) = {0,r,0,lc};
Point(3) = {r,0,0,lc};
Point(4) = {-r,0,0,lc};
Point(5) = {-l1,0,0,lc};
Point(6) = {l2,0,0,lc};
Point(7) = {-l1,h,0,lc};
Point(8) = {l2,h,0,lc};
Point(9) = {2*r,0.0,0,lc};
Point(10) = {-2*r,0.0,0,lc};
Point(11) = {0.0,0.2,0,lc};
Point(12) = {0.0,h,0,lc};
Point(13) = {r*Sin(Pi/4),r*Sin(Pi/4),0,lc};
Point(14) = {-r*Sin(Pi/4),r*Sin(Pi/4),0,lc};
Point(15) = {2*r*Sin(Pi/4),2*r*Sin(Pi/4),0,lc};
Point(16) = {-2*r*Sin(Pi/4),2*r*Sin(Pi/4),0,lc};
Point(17) = {l2,2*r*Sin(Pi/4),0,lc};
Point(18) = {-l1,2*r*Sin(Pi/4),0,lc};
Point(19) = {2*r*Sin(Pi/4),h,0,lc};
Point(20) = {-2*r*Sin(Pi/4),h,0,lc};
Circle(1) = {15,1,11};
Circle(2) = {13,1,2};
Circle(3) = {2,1,14};
Circle(4) = {11,1,16};
Circle(5) = {16,1,10};
Circle(6) = {14,1,4};
Circle(7) = {13,1,3};
Circle(8) = {15,1,9};
Line(9) = {9,3};
Line(10) = {15,13};
Line(11) = {11,2};
Line(12) = {16,14};
Line(13) = {10,4};
Line(14) = {12,11};
Line(15) = {19,15};
Line(16) = {12,19};
Line(17) = {19,8};
Line(18) = {8,17};
Line(19) = {17,15};
Line(20) = {6,9};
Line(21) = {6,17};
Line(22) = {16,20};
Line(23) = {20,12};
Line(24) = {10,5};
Line(25) = {5,18};
Line(26) = {18,16};
Line(27) = {7,18};
Line(28) = {20,7};
Line Loop(29) = {17,18,19,-15};
Plane Surface(30) = {29};
Line Loop(31) = {21,19,8,-20};
Plane Surface(32) = {31};
Line Loop(33) = {10,7,-9,-8};
Plane Surface(34) = {33};
Line Loop(35) = {1,11,-2,-10};
Plane Surface(36) = {35};
Line Loop(37) = {16,15,1,-14};
Plane Surface(38) = {37};
Line Loop(39) = {4,22,23,14};
Plane Surface(40) = {39};
Line Loop(41) = {11,3,-12,-4};
Plane Surface(42) = {41};
Line Loop(43) = {6,-13,-5,12};
Plane Surface(44) = {43};
Line Loop(45) = {26,5,24,25};
Plane Surface(46) = {45};
Line Loop(47) = {28,27,26,22};
Plane Surface(48) = {47};
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
