- Thank you received: 0
The (projective) Lehrenfeld-Schölberl stabilization function
5 years 6 months ago #1627
by Greg
The (projective) Lehrenfeld-Schölberl stabilization function was created 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
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
5 years 6 months ago - 5 years 6 months ago #1628
by schruste
Replied by schruste on topic The (projective) Lehrenfeld-Schölberl stabilization function
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
. The key is the two flag
and
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
ngsolve.org/docu/latest/i-tutorials/unit-2.8-DG/DG.html
.
Best,
Christoph
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)
Code:
highest_order_dc
Code:
hide_highest_order_dc
Best,
Christoph
Attachments:
Last edit: 5 years 6 months ago by schruste.
5 years 6 months ago #1629
by Greg
Replied by Greg on topic The (projective) Lehrenfeld-Schölberl stabilization function
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
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