请教如何编写一个source／sink边界条件

[hurricane007]

最近想做一个简化的模拟，需要一种边界条件，在这个边界上，某个组分的流量等于min(最大可能的流量，一个常数）。考虑准稳态情况下，菲克扩散定律 J =-D gradient (C)*n. 在使得边界网格内这个组分不为零的情况下，从这个网格向边界扩散的流量应该小于这个网格接受其他网格扩散的量。
那么问题来了…… 这个边界上的网格，除了边界面，其他面上的扩散流怎么求……求助一下大神们……