关于代码:我有一个代码,模拟浓度从advection-diffusion-reaction PDE在二维空间(X, Y)。获得的解决方案是使用全隐式有限差分方法,包括模拟媒体的功能与空间不同渗透率和反应常数(通过逆风调和平均数)。我已经能够测试的代码同质媒体和它将正常工作。这是一个解下列方程:$ ${对齐}\ \开始压裂{\部分C}{\部分t} + \微分算符。\离开(v C - D \微分算符{C} \右)= \αC \{对齐}$ $

问题:代码提供了现实的解决方案与均匀介质层的渗透率,如下所示。然而,一旦中开始得到更多的异构它开始把不切实际的结果(即负浓度随时间和大波动)。任何异质媒体猜测为什么我变得不切实际的结果(在类似的范围值均匀分层媒体)?

更新:我想我已经找到这个问题的负面价值观和浓度剧烈波动。我认为,问题不在于数值模型,但由于物理值。虽然我已经回到了现实的参数值的我变得不切实际的浓度,我必须进一步调整他们(减少对流速度)来得到真实的结果