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