[Gmsh] High-order visualization for hp-FEM solutions

[Beriot, Hadrien]

Dear Nicolas,



Thank you for your help. I followed your advice, and projected on the Lagrange basis before generating a results file for gmsh. It worked like a charm, I have been able to interface easily with gmsh and can now visualize hp-FEM solutions up to order 10 on very coarse meshes (see some nice snapshots below). It looks awesome, thanks.

Kind regards,

Hadrien


A duct mode, coarse T4 mesh, p=8, visu. order 8
[cid:image001.png at 01D2AD70.392FEBF0]

Diffraction of a plane wave by a sphere, coarse T10 mesh, visualization order 6, computational order 6
[cid:image003.png at 01D2AD72.61A00640]







-----Original Message-----
From: gmsh [mailto:gmsh-bounces at ace20.montefiore.ulg.ac.be] On Behalf Of Marsic, Nicolas
Sent: samedi 1 avril 2017 12:56
To: gmsh at onelab.info
Subject: Re: [Gmsh] High-order visualization for hp-FEM solutions



Dear Hadrien,



The current solution to visualize higher-order fields, is to project it into a Lagrange basis (the one used by gmsh for the high-order visu).

A better solution is perhaps to allow the user to enter his own basis functions, defining thus an InterpolationSheme per orientation, as you suggested.

But this is not possible yet...

The 'projection-on-a-Lagrange-basis' solution has the good point that only one well defined basis is used, and it is up to the user to accommodate to the interface...



If you want an example, you can have a look at the small_fem solver:

http://gitlab.onelab.info/gmsh/small_fem

It uses the 'projection-on-a-Lagrange-basis' solution.

The class postprocessing/FEMSolution (FEMSolution.cpp; FEMSolution.h and

FEMSolutionInclusion.h)

could be a good starting point if you want to get some inspiration :-).

By the way, I think that small_fem uses the same basis functions as hp-fem ;-).



Regarding the CG vs DG, the gmsh Lagrange basis function are defined in a CG way.

So you will end up with correct solution.

And finally, if you want a high-order vector solution, you have to project the 3 components of the vector field.



I hope this helps.

Best,

Nicolas.



On 28/03/17 12:15, Beriot, Hadrien wrote:

>

>

> Hello there,

>

>

>

> I have a question regarding the visualization of high-order FEM

> solutions. I would like to use this feature for visualizing highly

> oscillatory solutions arising from high-order p-FEM solutions (for

> instance used in frequency domain electromagnetics and acoustics)?

> There, by contrast to high-order DG, the difficulty I see is that each

> element has a unique edge/face orientation, which results pretty much

> in a different basis definition for each element.

>

>

>

> Considering the tetrahedral element for instance, there are usually 6

> different possible orientations per face, which yields a total of 1296

> different possible basis definitions. In my current understanding,

> this would require defining 6^4=1296 different "$InterpolationScheme"

> entities in the input file for gmsh (for the edges shape functions, a

> workaround can be found by simply multiplying the shape function

> contribution by +1 or -1 depending on the local edge orientation).

>

>

>

> Before I start the actual implementation to make some tests, I would

> like to know whether

>

>

>

> 1.       Has anyone tried the high-order visualization feature for

> viewing 3D p-FEM solutions (integrated Legendre basis)?

>

> 2.       Can the current implementation cope with a very large number of

> interpolation schemes (approx. one per element)? Would it still

> perform correctly on large models?

>

>

>

> Thanks for your help!

>

> Kind regards,

>

>

>

> Hadrien

>

>

>

>

>

>

>

>

>

> _______________________________________________

> gmsh mailing list

> gmsh at onelab.info<mailto:gmsh at onelab.info>

> http://onelab.info/mailman/listinfo/gmsh

>

_______________________________________________

gmsh mailing list

gmsh at onelab.info<mailto:gmsh at onelab.info>

http://onelab.info/mailman/listinfo/gmsh
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/gmsh/attachments/20170404/e0840651/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.png
Type: image/png
Size: 339782 bytes
Desc: image001.png
URL: <http://onelab.info/pipermail/gmsh/attachments/20170404/e0840651/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image002.png
Type: image/png
Size: 342945 bytes
Desc: image002.png
URL: <http://onelab.info/pipermail/gmsh/attachments/20170404/e0840651/attachment-0003.png>