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

Fo*(Bi)

于 2020-11-15 20:20:42 发布

阅读量1.6k

点赞数 1

分类专栏：数值传热文章标签： python

本文链接：https://blog.csdn.net/weixin_48615832/article/details/109708496

版权

本文介绍了一维稳态对流扩散问题的中心差分格式，结合陶文铨的数值传热学和李人宪的有限体积法基础，通过Python编程求解了一个边界条件固定的例子。问题中，流速u设为0.2，扩散系数Γ为0.1，物体长度为1，被划分为5个网格，左右边界温度分别为1和0。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

一维稳态对流扩散中心差分格式

参考书籍：陶文铨的数值传热学+李人宪的有限体积法基础
具体的理论可以参考此文章：有限体积法（5）——对流-扩散方程的离散
在这里插入图片描述

例题

一个长度L为1的规则物体，左边界温度恒为1，右边界温度恒为0，流速u=0.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 = 0.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(

最低0.47元/天解锁文章