一维热传导问题与C++描述

1一维热传导问题与C++描述

1.1理论推导

扩散问题的通用方程可以写为：

$$\begin{equation} \nabla\nabla\Gamma +S = 0\tag{1} \end{equation}$$

一维热传导问题可以讲上式子简化为来：

$$\begin{equation} \frac{d}{dx}(\Gamma\frac{d\phi}{dx}) +S = 0\tag{2} \end{equation}$$
将其在网格上离散可以得到下图的控制体单元表示：

将1.2式在网格上积分并利用Gauss定理就可以得到：

$$\begin{equation} \int_{{\Delta}V}\frac{d}{dx}(\Gamma\frac{d\phi}{dx})dV+\int_{{\Delta}V}SdV =({\Gamma}A\frac{d\phi}{dx})_{e}-({\Gamma}A\frac{d\phi}{dx})_w+\bar S{\Delta}V=0\tag{3} \end{equation}$$
${\Gamma}_{w}$与 ${\Gamma}_{e}$可以用线性插值的方法得到。1.3再在w,e面上的值即可以写为