python解zuobiaoxi方程_一维稳态有源导热——从方程到代码

1. 问题描述

如图，厚度为

的无限大平板，导热系数

，板内有均匀内热源

，表面A温度保持在

，表面B温度保持在

。求板内厚度方向温度分布。

2. 问题分析

一维稳态无源导热控制方程如下：

对于上述方程，利用有限体积法来进行离散求解。

离散求解PED方程的步骤为：

1. 离散控制域（网格划分）
2. 在每一个控制体上离散控制方程
3. 插值得到界面值，完成单元离散方程
4. 组装单元控制方程，形成整体控制方程组（Ax=b）
5. 求解代数方程组（直接，迭代）
6. 得到离散场变量

下面分别按照上面给出的六个步骤进行分析。

3. 问题求解

3.1 域离散

采用均匀网格将棒沿着长度方向离散，如图所示：

为了编写程序，我们需要容器来储存网格数据。对于一维问题，我们需要存储离散节点的位置坐标，单元的体积（在本例中为单元的长度乘以截面积）。

代码实现

import 

3.2 方程离散

在如图的控制体上，对控制方程进行积分得到：

上式中间部分继续展开&#