实现 Robin 边界条件(有限差分)

计算科学 pde 边界条件 数字 平流扩散
2021-12-21 07:59:34

我有兴趣将 Robin 边界条件应用于 1D 中的对流扩散问题。

在以下系统中,

Ct=D2Cx2vCx

Ct=x(DCxvC)
实现无通量边界条件,通量 在左右边界处设置为零。
N=DCxvC

我想知道通量中项的符号是否会在左右边界处变化。

根据wiki中给出的以下描述

在此处输入图像描述

在此处输入图像描述

有人可以解释上述方法是否是正确的实施方式吗?

1个回答

我认为第二个边界条件方程是不正确的。第一个应该适合两端。按照你的符号,域中的质量通量应该是: 包括边界点在内的所有地方。

N=vCDCx

请记住,如果质量通量的方向为正,则边界条件值的符号为正,反之亦然。