[Dune] how to integrate curvilinear tetrahedral elements into DUNE/Alugrid ?

[Benedikt Oswald]

Hello Andreas, thanks.

Then, I shall do so!

However, not being an expert in Dune internals programming,
where would I start ?

would the extension of dune / grid / geometry grid / coordfunction be the right place ?

And, more technically speaking, should copy the geometry grid directory into my
own source tree, implement the things there and, once it works, submit back to Dune?

What about naming the new implementation ?

Greetings, Benedikt




--------------------------------------------------------------------------------------------------------------------------------------------------------------
Dr. sc. techn. Benedikt Oswald - first engineer - LSPR AG - phone - +41 43 366 90 74
Technoparkstrasse 1, CH-8005 Zürich, benedikt.oswald at lspr.ch
--------------------------------------------------------------------------------------------------------------------------------------------------------------

On Jan 23, 2014, at 10:04 AM, Andreas Dedner <a.s.dedner at warwick.ac.uk> wrote:

> Hi Benedikt.
> You would need to extend GeometryGrid (something many of us would
> really appreciate). Basically GeometryGrid replaces the geometries of the host grid
> using the GenericGeometries which are only "linear". This part of GeometryGrid
> would need to be replaced with something higher order, e.g., some curvilinear approximation.
> Best
> Andreas
> 
> On 23.01.2014 08:32, Benedikt Oswald wrote:
>> Dear Dune
>> 
>> as it happens we are forced to consider drastic steps to increase the efficiency of our electromagnetic solvers,
>> especially in the frequency domain.
>> 
>> Recently, Peter Bastian suggested we use curvilinear elements in order to model our curved geometry and domain
>> boundaries. In fact, we believe that this will help us enormously since the element count will be much reduced.
>> 
>> Peter also suggested we use the geometry grid on top of ALUGrid in order to achieve this.
>> 
>> Now, my concrete questions is: given an existing curvilinear mesh (e.g. from gush), how would we do this ?
>> 
>> I read that geometry grid manages the vertices of the underlying grid, but, a curvilinear tetrahedral element has
>> vertices not only in the corners but also on the edges, faces and inside the volume.
>> Perhaps my question is trivial, nevertheless, I'd appreciate your suggestions.
>> 
>> I should say that gmsh is capable of generating quite nice curvilinear tetrahedral meshes, at least up to 3rd order.
>> 
>> Greetings and thanks in advance, Benedikt
>> 
>> 
>> 
>> --------------------------------------------------------------------------------------------------------------------------------------------------------------
>> Dr. sc. techn. Benedikt Oswald - first engineer - LSPR AG - phone - +41 43 366 90 74
>> Technoparkstrasse 1, CH-8005 Zürich, benedikt.oswald at lspr.ch
>> --------------------------------------------------------------------------------------------------------------------------------------------------------------
>> 
>> 
>> 
>> _______________________________________________
>> Dune mailing list
>> Dune at dune-project.org
>> http://lists.dune-project.org/mailman/listinfo/dune
> 
> _______________________________________________
> Dune mailing list
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20140123/4b387f45/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 841 bytes
Desc: Message signed with OpenPGP using GPGMail
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20140123/4b387f45/attachment.sig>