这篇博客详细介绍了计算传热学中的一维常物性无源对流导热微分方程的数值解方法，包括Taylor级数展开法、控制容积积分法、Gauss-Seidel迭代法和TDMA方法。通过对网格划分和离散方程的讨论，展示了如何求解和验证数值解的网格独立性。
计算传热学第一次大作业
1 Taylor级数展开法
1.1 网格划分
Taylor级数展开发的网格划分使用外点法
1.2 离散方程表达式
ρ u d ϕ d x = Γ d 2 ϕ d x 2 (1) \rho u\frac{d\phi}{d x}=\Gamma \frac{d^{2}\phi}{dx^{2}}\tag{1} ρudxdϕ=Γdx2d2ϕ(1)
将 ϕ \phi ϕ分别在 i i i+1点和 i i i-1点在 i i i点展开
ϕ ( i + 1 ) = ϕ ( i ) + d ϕ d x ∣ i δ x + d 2 ϕ d x 2 ∣ i δ x 2 2 ! + … … (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{2} ϕ(i+1)=ϕ(i)+dxdϕ∣∣∣∣iδx+dx2d2ϕ∣∣∣∣i2!δx2+……(2)
ϕ ( i − 1 ) = ϕ ( i ) − d ϕ d x ∣ i δ x + d 2 ϕ d x 2 ∣ i δ x 2 2 ! + … … (3) \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} ϕ(i−1)=ϕ(i)−dxdϕ∣∣∣∣iδx+dx2d2ϕ∣∣∣∣i2!δx2+……(3)
联立(2)、(3)式得
d ϕ d x ∣ i = ϕ i + 1 − ϕ i − 1 2 δ x , o ( δ x 2 ) (4) \left.\frac{d\phi}{dx}\right|_{i}=\frac{\phi_{i+1}-\phi_{i-1}}{2\delta x},o(\delta x^{2})\tag{4} dxdϕ∣∣∣∣i=2δxϕi+1−ϕi−1,o(δx2)(4)
同样的，联立(2)、(3)式得
d 2 ϕ d x 2 ∣ i = ϕ i + 1 − 2 ϕ i + ϕ i − 1 δ x 2 , o ( δ x 2 ) (5) \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} dx2d2ϕ∣∣∣∣i=δx2ϕi+1−2ϕi+ϕi−1,o(δx2)(5)
将(4)、(5)式带入(1)式
ρ u ϕ i + 1 − ϕ i − 1 2 δ x = Γ ϕ i + 1 − 2 ϕ i + ϕ i − 1 δ x 2 (6) \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} ρu2δxϕi+1−ϕi−1=Γδx2ϕi+1−2ϕi+ϕi−1(6)
最低0.47元/天 解锁文章
扫一扫
1325
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
