Modeling a 2 sub-domain diffusion PDE with concentration discontinuity

Hi all, I'm modeling mass transfer between liquid and gas. and in the interface between two part (x=L), i know that i have a fixed value of the concentration, due to liquid-vapor equilibrium. thus, from the left side of x=L (boundary of my left subdomain,the gas), i know the value of the concentration.

It also causes discontinuity in the value of C - as the concentration in the gas phase in much lower. it's like there is a "jump" in a water concentration right above a water cup bewteen the concentration in the cup and the humidity (=gaseous water) above it - but there is a flux between those two (hope this example helped and did not confuse)

can anyone know how can i solve this one? (a code example will be helpful)

Thanks, Bar

Answer 1

[Giuseppe Pontrelli]

Hi Bar, I have a very similar problem, and I am convinced the problem should be split in two sub-problems to be solved in iterative way: 1) solve one domain (with the flux or concentration at the boundary coming from the other) 2) solve the other domain (with the flux or concentration at the boundary coming from the first) and so on, till convergence, hopefully. if you solve in one domain, your solution is supposed to be globally continuous, that is not!

I am trying to implement this for my case, keep in contact, then. Giuseppe