[Dune] parallel DUNE grid for implementing a spherical shell (UGGrid)

[Christian Engwer]

On Mon, Nov 05, 2012 at 04:45:24PM +0000, Eike Mueller wrote:
> Dear dune-list,
> 
> can anyone help me with the following question about UG grid:
> 
> I want to create a parallel grid representing a spherical shell, with one level in the radial direction only. I have tried this with UGGrid and the GridFactory class. To create a cubed sphere grid, I created 2x8 vertices, describing two nested cubes with factory.insertVertex(). I then create the corresponding 12 elements with factory.insertElement() and describe the projection onto the sphere by using factory.insertBoundarySegment() for each of the 2x12 faces (12 inner faces and 12 outer faces). Before doing any refinement I get the macro grid in the attached file cubed sphere_macro.vtu.
> If I then call globalRefine(), the grid gets refined in all three dimensions and I end up with the attached refined grid cubed sphere_refined.vtu.
> 
> Now my two questions are as follows:
> 
> (1) How can I ensure that the grid is not refined in the radial direction, i.e. I always only end up with one layer of cells in the radial direction? I guess I have to run over the grid and mark cells for refinement, but is this documented anywhere in more detail?

there are some hidden features in ug to do anisotropic refinement, but
I don't think they are exported to dune and they will turn off all
closure computations in ug, so you might end up with non-conforming
grids.

> (2) The parallel load balancing does not seem to work. If I run on 24 processes and I call globalRefine(), followed by loadBalance(), some processes end up with no part of the grid. I have checked that I use the parallel version of UG grid.

load-balancing with ug is still very much work in progress... sorry.

Christian

> 
> I attach the code I used to generate the grid, as well as the .vtu files mentioned above.
> 
> Thanks a lot for any help,
> 
> Eike
> 
> PS: Alternatively, I would also be equaully happy if I can create a PARALLEL 2d surface grid for the surface of the sphere (again, either cubed sphere or a icosahedral grid). However, as far as I can see it, there is no parallel implementation for this, we've tried 2d ALUGRID, parallelised with the metagrid approach, but this does not scale beyond a few hundred processors as the entire grid is stored on each processor and my code runs out of memory.
> 




> _______________________________________________
> Dune mailing list
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune


-- 
Prof. Dr. Christian Engwer 
Institut für Numerische und Angewandte Mathematik
Fachbereich Mathematik und Informatik der Universität Münster
Einsteinstrasse 62
48149 Münster

E-Mail	christian.engwer at uni-muenster.de
Telefon	+49 251 83-35067
FAX		+49 251 83-32729