关于非均一扩散系数的扩散项离散问题
-
采用二阶精度有限体积法,扩散系数是关于空间位置的函数,面心扩散系数直接给定扩散系数的解析值,扩散项用中心差分格式,时间项用一阶欧拉隐式,计算结果是变量在全域的总积分值会逐渐减少,令扩散系数为常数时计算结果是守恒的。
这种扩散系数应该怎么离散啊。。。
大佬们指点下 -
已经有$\alpha$解析解了还需要离散么?
守恒是什么守恒?全局的$\phi$守恒?
-
单元体积乘单元变量值(在全计算域所有单元求和),这个积分值和初始场的积分值相等,既是这里指的守恒
-
如果不考虑时间项
\begin{equation}
\int^{f+1}_ {f} \frac{\partial}{\partial x}\left(\alpha \frac{\partial \phi}{\partial x}\right)\mathrm{d}x=\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f}
\end{equation}
把后面两项$\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_{f}$把所有的面连起来,看起来是守恒的,剩下的只有边界条件