Discussion Forum

I have done what you asked by using two physics nodes, one for each domain. Let's say that in domain 1 you have concentration C1 and in domain concentration C2 (they are the same molecules but in different domains). At the internal boundary, in domain 1 define

C1 = P*C2

and in domain 2

C2 = C1/P

and P is the partition coefficient between the domains.

For closer details, you can have a look at

dx.doi.org/10.1016/j.mbs.2014.03.014

BR
Lasse

Hi Lasse,

thanks for the reply. How about define flux discontinuity when have two physics nodes? Shall I use the same way you mentioned below for defining flux discontinuity or I can use identity pairs for it?

regards,
Hansong

I do not understand these identity pairs, never used them.

Flux discontinuity is formed if you define different fluxes in domain 1 and 2 at the boundary. But as you certainly understand, flux discontinuity means accumulation or depletion of material at the interface.

BR
Lasse

As what you said above, I use two physics nodes, one for each domain. And I set the concentration and flux are continuous across the internal boundary, but now it is not converged. Do you know any idea why? The difference is just that when I use 1 physics node defined on two domain, by default the concentration and the flux are continuous and it is converged.

thanks,
Hansong

I need to see the model file, can you submit to me privately?

Lasse

Hi Lasse,

I have figured out my problem already. It is just the sign of flux at the internal boundary which cause not convergence.

thanks,
Hansong