如何拿COMSOL使用Poisson-Nernst-Planck模型仿真流动电势的产生呢？

我想这类问题最简单的物理情境应该是通过压力驱动的方式，使溶液进入直圆管中，固壁接触到电解质溶液后带电，带电表面附近产生双电层，扩散层中富集的反离子随溶液朝下游流动，此过程对应流动电流Is的产生，且通道出口处积聚反离子（“反”相对于固表电荷的电性而言）的净电荷，累积的电荷使通道两端产生电位差，即流动电势ψs，对应电场为“流动诱导电场”Es，阻碍反离子朝下游的电迁移，形成与流动电流相反的传导电流Ic，出口处净电荷不断累加，流动电势ψs越来越大，很快传导电流Ic增大至与流动电流Is相等，即达到稳态，出口处不再积聚更多净电荷。

微通道中流体流动和离子传质的瞬态时间非常短(通常在毫秒或更小的数量级)，可以忽略不计。因此可求解稳态下电场中的Possion方程，流场中的Navier-Stokes方程和离子浓度场的Nernst-Plank方程。采用es静电场接口、spf层流接口和tds稀物质传递接口全耦合求解电场、速度场和浓度场，壁面ζ电势、压力梯度、壁面无滑移和入出口的浓度作为已知边界条件。仿真过程中遇到的很大一个困难就是流动电势ψs的求解，因为Is+Ic=0,即∫ρudS+σEs·S=0，而出口处ψs=Es·l（此处l为通道长度），若我通过这个电流平衡条件表征出流动电势ψs，以此作为静电场中出口处的边界条件，求解会出现不收敛等问题，该怎么去解决呢？抑或这个求解思想本身就是错误的吗？那么真正正确的求解方式就该是什么呢？感谢您拨冗阅读，为我拨开眼前的重重迷雾。