一维稳态对流扩散问题，无源项，一阶迎风差分格式的python程序

一维稳态对流扩散无源项一阶迎风差分格式

参考书籍：陶文铨的数值传热学+李人宪的有限体积法基础
具体的理论可以参考此文章：有限体积法（7）——迎风格式

例题

一个长度L为1的规则物体，左边界温度恒为1，右边界温度恒为0，流速u=2，扩散系数Γ为0.1，密度 ρ为1。
本题划分5个网格。

代码：

#一维稳态对流扩散问题
#无源项，一阶迎风差分格式
from numpy import *
from numpy.linalg import *
import matplotlib.pyplot as plt
L = 1  #物体总长度
n = 5  #划分网格数
I = 0.1  #扩散系数
p = 1  #密度
u = 2  #速度
#边界条件
tA = 1  #左边界A点，定壁温，为1
tB = 0  #右边界B点，定壁温，为0

dx = L/n  #控制体宽度
De, Dw = I/dx, I/dx
DA, DB = 2*I/dx, 2*I/dx
Fe, Fw, FA, FB = p*u, p*u, p*u, p*u

A = zeros((n,n)) #创建零数组
B = zeros((n,1))

if u