The (projective) Lehrenfeld-Schölberl stabilization function

More
5 years 6 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
5 years 6 months ago - 5 years 6 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: 5 years 6 months ago by schruste.
More
5 years 6 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.096 seconds