Forum Message

 

 

We have moved the forum to https://forum.ngsolve.org . This is an archived version of the topics until 05/05/23. All the topics were moved to the new forum and conversations can be continued there. This forum is just kept as legacy to not invalidate old links. If you want to continue a conversation just look for the topic in the new forum.

Notice

The forum is in read only mode.

The (projective) Lehrenfeld-Schölberl stabilization function

More
4 years 11 months ago #1627 by Greg
Hi,

I'm trying to implement an HDG method with reduced stabilization for second-order diffusion problems. Precisely, I want to use the (projective) Lehrenfeld-Schöberl stabilization function.

Is this function already available in NGSolve, or should I create it?

Regards,
Gregory
More
4 years 11 months ago - 4 years 11 months ago #1628 by schruste
Hi Gregory,

It's available in NGSolve. Attached is a minimal example. You may also want to have a look at the documentation of FacetFESpace. To do so type
Code:
help(FacetFESpace)
. The key is the two flag
Code:
highest_order_dc
and
Code:
hide_highest_order_dc
which do the projected-jumps trick as explained in 2.2.2. in this (preprint version of the ) paper: arxiv.org/pdf/1508.04245.pdf . For the even more basic stuff about HDG in NGSolve I refer you to the i-tutorial 2.8, see

File Attachment:

File Name: hdg.py
File Size:1 KB
ngsolve.org/docu/latest/i-tutorials/unit-2.8-DG/DG.html .

Best,
Christoph
Attachments:
Last edit: 4 years 11 months ago by schruste.
More
4 years 11 months ago #1629 by Greg
Hi Christoph,

Thank you for your response.
I follow the same steps for the FESpace as in the attached script. In the function, I have set "highest_order_dc = true" and "hide_highest_order_dc = True".
Now it works fine.


Thanks,
Gregory
Time to create page: 0.124 seconds