相间传质问题

zy-nwu

对于一个对流扩散方程，时间上每求解一步，求解出来 得到一条浓度曲线C（z），用此浓度分布曲线与预先知道的溶解度曲线C1(z)做比较，在z上一点，If C>C1, If C<C1， C>C1,则此点的浓度度为C1，也就是说最大浓度为C1，不能让C大于C1（C-C1部分被转化为其他相中的M, M不参与对流扩散）；然后继续求解下一个时间步的浓度分布。请问在此过程中，为了保持质量方程的守恒，我是不是应该修改源项R为：R=-（C-C1）*(C >C1)   但加上（C>C1）这一限制后，求解就会报错

与Comsol自带的碳沉淀（carbon deposition）例子有些像，但碳沉淀例子是只要在反应区就会有物质转化，转化率定为R；而我这个问题是在区间满足大于溶解度C1的部分才转化。


或者是我的修改方法压根错了，请帮忙指正！谢谢大家！