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

Benedikt Oswald benedikt.oswald at lspr.ch
Thu Jan 23 16:12:54 CET 2014 

ok, thanks, I will go on as suggested, greetings, Benedikt



On Jan 23, 2014, at 2:40 PM, Andreas Dedner <a.s.dedner at warwick.ac.uk> wrote:

> A final suggestion I could make:
> start a new project on the dune user wiki.
> Move the geometrygrid files from dune-grid/dune/grid/geometrygrid 
> to newproject/dune/grid/geometrygrid
> Now you should be able to test everything including dune-geometrygrid as
> before and everything can be easily copied back to dune-grid when completed.
> Andreas
> 
> On 23.01.2014 14:04, Benedikt Oswald wrote:
>> hello DUNE
>> 
>> so, reflecting the comments the approach to get curvilinear meshes
>> should be along those lines:
>> 
>> 1) extend the gush reader to read tetrahedral grids up to 3rd order
>>     gmsh at present works up to 3rd order, higher order mesh generation
>>     crashes
>> 
>> 2) use the grid factory facility to create the grid in dune
>> 
>> 3) implement the curvilinear geometry grid in
>> 
>> 	a) multilineargeometry.hh
>> 	b) geometry-grid
>> 
>> 
>> 4) since we would like to integrate this into Dune as soon as it is mature,
>>     hearing all the suggestions, where should we put the work ?
>> 
>>     keep it on our machines ? keep it in a project on the Dune user wiki ?
>> 
>>  Greetings, Benedikt
>> 
>>    
>> 
>> 
>> 
>> 
>> 
>> On Jan 23, 2014, at 11:48 AM, Oliver Sander <sander at igpm.rwth-aachen.de> wrote:
>> 
>>> Am 23.01.2014 10:25, schrieb Andreas Dedner:
>>>> What I would suggest would be to setup a project on the dune user wiki and
>>> A good idea, but while doing that I suggest to keep in mind that in the long
>>> run it may be preferable to merge your code into GeometryGrid proper.
>>> You are implementing a more general version of GeometryGrid; it is not obvious
>>> why it should be a completely separate implementation (except to get you started).
>>> 
>>>> giving it public access - that hopefully will give you some help from the
>>>> experienced dune developers if you get stuck.
>>>> I'm terrible with naming so will not make a suggestion but
>>>> CurvilinearGeometryGrid and simply CurvilinearGrid comes to mind.
>>> -1, I think this is too unspecific.  If you intend to implement polynomial
>>> geometries how about PolynomialGrid?  Or QuadraticGrid or SecondOrderGrid
>>> if you only want to do quadratic geometries?
>>> 
>>> Best,
>>> Oliver
>>> 
>>>> Best
>>>> Andreas
>>>> 
>>>> On 23.01.2014 10:13, Benedikt Oswald wrote:
>>>>> 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
>>>>> <mailto:benedikt.oswald at lspr.ch>
>>>>> --------------------------------------------------------------------------------------------------------------------------------------------------------------
>>>>> 
>>>>> On Jan 23, 2014, at 10:04 AM, Andreas Dedner <a.s.dedner at warwick.ac.uk
>>>>> <mailto: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
>>>>>>> <mailto: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 <mailto: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
>>>> 
>>> _______________________________________________
>>> 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
> 
> _______________________________________________
> 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/3d42f7d1/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/3d42f7d1/attachment.sig>

More information about the Dune mailing list