一维常物性无源对流导热微分方程数值解(大概)

最新推荐文章于 2022-07-02 19:52:37 发布

qq_41640870

最新推荐文章于 2022-07-02 19:52:37 发布

阅读量1.7k

点赞数 2

本文链接：https://blog.csdn.net/qq_41640870/article/details/104650455

版权

这篇博客详细介绍了计算传热学中的一维常物性无源对流导热微分方程的数值解方法，包括Taylor级数展开法、控制容积积分法、Gauss-Seidel迭代法和TDMA方法。通过对网格划分和离散方程的讨论，展示了如何求解和验证数值解的网格独立性。

摘要由CSDN通过智能技术生成

计算传热学第一次大作业

1 Taylor级数展开法

1.1 网格划分

Taylor级数展开发的网格划分使用外点法
在这里插入图片描述

1.2 离散方程表达式

$\rho u\frac{d\phi}{d x}=\Gamma \frac{d^{2}\phi}{dx^{2}}\tag{1}$
将 $\phi$ 分别在 $i$ +1点和 $i$ -1点在 $i$ 点展开
$\phi(i+1) = \phi(i)+\left.\frac{d\phi}{dx}\right|_{i}\delta x+\left.\frac{d^{2}\phi}{dx^{2}}\right|_{i}\frac{\delta x^{2}}{2!}+…… \tag{2}$
$\phi(i-1) = \phi(i)-\left.\frac{d\phi}{dx}\right|_{i}\delta x+\left.\frac{d^{2}\phi}{dx^{2}}\right|_{i}\frac{\delta x^{2}}{2!}+…… \tag{3}$
联立(2)、(3)式得
$\left.\frac{d\phi}{dx}\right|_{i}=\frac{\phi_{i+1}-\phi_{i-1}}{2\delta x},o(\delta x^{2})\tag{4}$
同样的，联立(2)、(3)式得
$\left.\frac{d^{2}\phi}{dx^{2}}\right|_{i}=\frac{\phi_{i+1}-2\phi_{i}+\phi_{i-1}}{\delta x^{2}},o(\delta x^{2})\tag{5}$
将(4)、(5)式带入(1)式
$\rho u \frac{\phi_{i+1}-\phi_{i-1}}{2\delta x}=\Gamma \frac{\phi_{i+1}-2\phi_{i}+\phi_{i-1}}{\delta x^{2}}\tag{6}$