扩散方程的数值解法

扩散方程的数值解法

一维稳态导热的离散

​ 一维稳态导热问题的数学描述可表示为：
$$\frac{1}{A(x)}\frac{d}{dx}[\lambda A(x)\frac{dT}{dx}]+S=0\text{\hspace{1em}(1)}$$
​ 其中$S$表示为源项，假设源项可表示为温度的线性函数，即$S=S_C+S_P{T}_P$

​ 对方程$(1)$式两边同乘$A(x)$ ，并对控制体$P$ 积分，可得
$$T_P(\frac{{\lambda }_e{A}_e}{(\delta x{)}_e}+\frac{{\lambda }_w{A}_w}{(\delta x{)}_w}-S_P{A}_P\Delta x)=T_W(\frac{{\lambda }_w{A}_w}{(\delta x{)}_w})+T_E(\frac{{\lambda }_e{A}_e}{(\delta x{)}_e})+S_C{A}_P\Delta x\text{\hspace{1em}(2)}$$
​ 最终可将其简化为：
$$a_P{T}_P=a_W{T}_W+a_E{T}_E+b\text{\hspace{1em}(3)}$$

​ $(3)$式即为一维稳态导热问题的最终离散形式，这时的问题关键在于表示控制体上的导热系数${\lambda }_w,{\lambda }_e$，确定界面上的导热系数，方法包括算术平均(arithmetic mean) 和调和平均( harmonic mean)两种，经验证，调和平均法更接近物理事实。

一维非稳态导热的离散

​ 一维非稳态导热问题的数学描述为：
KaTeX parse error: Undefined control sequence: \part at position 15: \rho c \frac{\̲p̲a̲r̲t̲ ̲T}{\part \tau} …

​ 对控制体在$[\tau ,\tau +\Delta \tau ]$其进行积分，
( ρ c ) P A P Δ x ( T P τ + Δ τ − T P τ ) = ∫ τ τ + Δ τ [ λ e A e ( T E − T P ) ( δ x ) e − λ w A w ( T P − T W ) ( δ x ) w ] d τ + ∫ τ τ + Δ τ ( S C + S P T P ) A P Δ x d τ (5) (\rho c)_P A_P \Delta x (T_P^{\tau + \Delta \tau} - T_P^\tau) = \int_\tau^{\tau + \Delta \tau} [\frac{\lambda_e A_e (T_E - T_P)}{(\delta x )_e} - \frac{\lambda_w A_w(T_P - T_W)}{(\delta x)_w}]d\tau \\ + \int_\tau^{\tau + \Delta \tau}(S_C + S_PT_P) A_P \Delta x d\tau \tag{5} (ρc)P​AP​Δx(TPτ+Δτ​−TPτ​)=∫ττ+Δτ​[(δx)e​λe​Ae​(TE​−TP​)​−(δx)w​λ