There is an inherent problem: N-P uses electroneutrality condition which you try to overrule by P-B. In an older version of Comsol there was physics "N-P without electroneutrality" (or something like that). It has been removed - maybe Comsol developers realized that it did not work properly?
I have used P-B according to Gouy-Chapman theory, i.e.
Nabla^2 phi = - SUM{c_i^b*exp(-(z_i*F/RT)*phi)}
In a physiological solution the electric double layer is very thin, ca. 1 nm, so do you really need that?
br
Lasse
I have used P-B according to Gouy-Chapman theory, i.e.
Nabla^2 phi = - SUM{c_i^b*exp(-(z_i*F/RT)*phi)}
In a physiological solution the electric double layer is very thin, ca. 1 nm, so do you really need that?
br
Lasse