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.

DG method : how to construct a weighted mean operator on an interior interface

More
4 years 8 months ago #1773 by vfont
Hi,
I need to solve the heterogeneous diffusion equation with the DG SWIP method. To this aim, I need to construct a new weighted mean operator on interior interfaces where the weighting coefficients depend on the (constant) value of the diffusivity defined on both sides of the interface. How to access the coefficient value of the neighboring element? Can I access it via the .Other() command like this:

k = CoefficientFunction(....)
w1=k/(k+k.Other())
w2=k.Other()/(k+k.Other())
wmean_gradu = w1*grad(u)+w2*grad(u.Other())

Any idea?
Regards
Vincent
Time to create page: 0.151 seconds