Domain dependent parameter

Hi All,

I was trying to impose domain dependent parameter. I have domain [0,2]\times[0,2] and a subdomain [0,1]\times[0,1]. I want to use nu = 20 in my subdomain (INCLUDING BOUNDARY WHERE x = 1 or y = 1) and use nu = 10 in the complement part of the subdomain. I am using:

real nup1 = 20;
real nus1 = 10;
func nus = (x >= 0 && x <= 1 && y >= 0 && y <= 1) ? nup1 : nus1;

But interestingly, some of the part of the boundary of the subdomain have values I expect to see in the complement domain. Please see the following mesh diagram.

This becomes more clear for another larger value parameter
lambda (see the following figure).

I thought this was interesting to share. I would greatly appreciate any feedback. Thank you.

Best,
Hemanta

Just use func nus = (x >= 0 && x <= 1+eps && y >= 0 && y <= 1+eps) ? nup1 : nus1;

That helps! Thank you! :slightly_smiling_face: