[Libmesh-users] Solving PDE with subdomains as BCs

[Viviana P. B. <vpa...@uc...>]

Hello,

I am trying to solve the diffusion equation with two types of Dirichlet
BCs. Let's say for example:
df/dt = lapl(f). f(0<x<0.5,y,t) = 1; f(0.5<x<1,y,t) = 0.
I was thinking of taking advantage of subdomains in which subdomain_id = 1
is the bulk. subdomain_id = 2 is surface1. subdomain_id=3 is surface 2.
Once the subdomains are read correctly, I could assemble the RHS
accordingly by using active_local_subdomain_elements_begin( subdomain_id ).
I was able to successfully assemble the corresponding terms using both type
of subdomains but I'm seeing that although it's solved simultaneously, the
solution of one subdomain is not being projected onto the other one. Is
this something I should project manually?

Thank you in advance.

Best,
Viviana.