下载文档原格式
一维热传导方程的差分法一维热传导方程是指在一维空间中，描述材料内部温度分布随时间的变化过程的方程式。
可以表示为：$$\frac{\partial u(x, t)}{\partial t} = \alpha \frac{\partial^2 u(x,t)}{\partial x^2}$$其中，$u(x, t)$ 是空间坐标为 $x$，时间为 $t$ 时的温度，$\alpha$ 是热扩散系数。
在边界条件确定的情况下，可以得到一维热传导方程的解。
然而，在实际应用中，解析解并不总是容易或可行的，因此需要使用数值方法进行近似求解。
常用的数值方法有有限差分法、有限元法和谱方法等，其中有限差分法是最为简单、易于实现的方法之一。
有限差分法的基本思想是将连续的空间区间离散化为若干个节点，将时间轴离散化为若干个时间步。
在空间和时间轴两个方向上，分别对热传导方程进行差分，得到离散的差分方程组，从而可以求得数值解。
在一维热传导方程的差分过程中，我们首先需要将空间区间 $(0, L)$ 划分为 $N$ 个等间距的节点，每个节点间距为 $h = \frac{L}{N}$。
我们使用 $u_i^n$ 表示节点$i$ 在时间步 $n$ 时的温度，其中 $i = 0, 1, ..., N$，$n = 0, 1, ..., M$。
接下来，我们对一维热传导方程进行中心差分，得到：$$\frac{u_i^{n+1} - u_i^n}{\Delta t} = \alpha \frac{u_{i+1}^n - 2u_i^n + u_{i-1}^n}{h^2}$$其中，$\Delta t$ 是时间步长。
可以将上式改写为：其中，$r = \frac{\alpha \Delta t}{h^2}$。
由于在一维空间中，有两个边界，因此需要对边界进行特殊处理。
常见的边界处理方式有三种：1. 固定边界：将边界上的温度固定为某一值；2. 自然边界：假设边界处热通量为零，从而根据傅里叶定律可以求得边界处温度的梯度值，从而推算出边界的温度值；3. 第二类边界：将边界节点的温度根据边界条件与内部的节点做差分，从而计算出边界节点的温度值。
分类号密级UDC学位论文一维辐射传递方程的谱方法求解作者姓名：张大伟指导教师：李本文 教授东北大学材料电磁过程研究教育部重点实验室申请学位级别：硕士学科类别：工学学科专业名称：工程热物理论文提交日期：2008年2月22日论文答辩日期：2008年2月25日学位授予日期：2008年3月 答辩委员会主席：王恩刚 教授评阅人：陈海耿（东北大学教授）聂宇宏（江苏科技大学副教授）东北大学2008年2月A Dissertation in Engineering ThermophysicsSPECTRAL METHOD FOR SOLVINGONE-DIMENSIONAL RADIATIVETRANSFER EQUATIONby Zhang DaweiSupervisor: Professor Li BenwenNortheastern UniversityFebruary 2008独创性声明本人声明，所呈交的学位论文是在导师的指导下完成的。
论文中取得的研究成果除加以标注和致谢的地方外，不包含其他人己经发表或撰写过的研究成果，也不包括本人为获得其他学位而使用过的材料。
与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示谢意。
学位论文作者签名：日期：学位论文版权使用授权书本学位论文作者和指导教师完全了解东北大学有关保留、使用学位论文的规定：即学校有权保留并向国家有关部门或机构送交论文的复印件和磁盘，允许论文被查阅和借阅。
本人同意东北大学可以将学位论文的全部或部分内容编入有关数据库进行检索、交流。
（如作者和导师不同意网上交流，请在下方签名；否则视为同意。
）学位论文作者签名：导师签名：签字日期：签字日期：一维辐射传递方程的谱方法求解摘要对于辐射传递方程的求解已经发展了很多方法，区域法、热通量法、蒙特卡洛法、离散坐标法等。
综合看来，基于微分形式辐射传递方程离散的方法具有很好的适应性，上述传统方法在求解辐射传递方程时各有优点，但也有其本身不可克服的缺点。
发展一种高精度、高效率、简便的数值求解辐射传递方程的方法仍然必要。
谱方法可以很好的克服低阶方法的不足，能高效地得到高精度的计算结果，而其所用的计算方法简便。
本文利用谱方法解决半透明介质中的辐射换热问题，主要工作分两部分：1．利用谱方法求解一维直角坐标系下的辐射传递方程，分别对不考虑介质散射和考虑介质散射两种情况进行计算，研究了该方法的数值稳定性，并与精确解或其他方法的数值解进行比较。
结果表明，谱方法适用于一维直角坐标下辐射传递方程的求解，而且相对于其他方法具有较高的精度。
2．利用谱方法解决一维球面坐标系下的辐射换热问题以及一维球面坐标系下辐射与导热耦合问题。
通过对谱方法求解结果与文献的结果进行比较和误差分析，验证了谱方法能简便、高效、高精度的解决一维球面坐标系下辐射问题以及辐射与导热相耦合的问题。
关键词辐射换热；导热；Chebyshev配置点谱方法；半透明介质Spectral Method for Solving One Dimensional RadiativeTransfer EquationAbstractMany numerical methods have been developed for solving radiative heat transfer problems in semitransparent media, such as, zonal method, heat flux method, Monte-Carlo method, and discrete ordinates method, etc. By comprehensive comparison, the methods based on discretization of the differential form of radiative transfer equation possess very good flexibility. However, there are still some disadvantages existing in the former methods. It is necessary to develop a high accurate, efficient and simple numerical method to solve radiative transfer equations.Spectral methods can somewhat overcome the drawbacks of low order methods, and can provide accurate results, while their solution procedures of spectral methods are something efficient and simple. The author develops and studies the characteristics of spectral methods to solve radiative transfer equation in semitransparent media. The present research contains two parts:1. The spectral method is used to solve radiative transfer equation in one-dimensional Cartesian coordinate systems. The numerical stability of spectral methods is investigated for solving radiative heat transfer under the both cases of non-scattering media and scattering media. The present numerical solutions are compared with those of references by exact or other numerical methods. The results show that the spectral methods have good accuracy and efficiency for solving one-dimensional radiative heat transfer problems.2. Spectral methods are presented to solve radiative heat transfer equations and coupled conduction-radiation problems in spherical media. An error comparison is made between present results and those from references. The final conclusions can be obtained that: the spectral methods can be used to solve radiative heat transfer alone, and can also to solve the coupled conduction-radiation problems in one-dimensional spherical media system efficiently.Keywords: radiative heat transfer; conduction; Chebyshev collocation spectral method;semitransparent media目录独创性声明 (I)摘要 (I)Abstract (I)目录 (I)第1章 绪论 (1)1.1研究背景及意义 (1)1.2辐射换热计算模型简介 (1)1.2.1区域法（Zone method） (1)1.2.2热通量法（Heat Flux method） (1)1.2.3蒙特卡洛法（Monte-Carlo method） (1)1.2.4离散坐标法（Discrete Ordinates method） (1)1.3谱方法介绍 (1)1.3.1谱方法的研究现状 (1)1.3.2谱方法的特点 (1)1.3.3谱方法的分类 (1)1.3.4谱方法的优缺点 (1)1.4 Chebyshev配置点谱方法简介 (1)1.4.1 Chebyshev 多项式 (1)1.4.2 Chebyshev 多项式的性质 (1)1.4.3几种Chebyshev配置点 (1)1.5本文的研究内容及意义 (1)第2章 一维辐射传递方程的谱方法求解 (1)2.1 Chebyshev配置点谱方法求解辐射传递方程的准备工作 (1)2.1.1配置点的选择 (1)2.1.2配置点均匀化处理 (1)2.1.3计算区域的转化 (1)2.1.4系数矩阵的确定 (1)2.1.5边界条件的处理 (1)2.2辐射传递方程及其特性 (1)2.3辐射传递方程的Chebyshev配置点谱方法求解方案 (1)2.4几个算例的验证 (1)2.4.1不考虑介质散射的辐射换热 (1)2.4.2考虑介质散射的辐射换热 (1)2.5本章小结 (1)第3章 一维辐射传递方程与导热方程耦合的谱方法求解 (1)3.1几种常见误差 (1)3.2导热方程及其特性 (1)3.3一维导热方程与辐射传递方程耦合的配置点谱方法求解 (1)3.3.1一维辐射传递方程配置点谱方法求解的准备工作 (1)3.3.2无量纲各物理量的设定 (1)3.3.3无量纲导热方程Chebyshev配置点谱方法求解方案 (1)3.3.4无量纲一维辐射传递方程的Chebyshev配置点谱方法求解 (1)3.4计算实现及求解步骤 (1)3.5几个算例的验证 (1)3.5.1算例一：只考虑辐射换热问题 (1)3.5.2算例二：导热与辐射换热耦合问题 (1)3.6本章小结 (1)第4章 结论 (48)参考文献 (49)致谢 (1)个人简历 (1)第1章绪论1.1研究背景及意义热量传递有三种基本方式：导热、对流和热辐射。
在高温状态进行的热量传递中，辐射换热处于主导地位，如工程用燃烧设备中，由于火焰对周围壁面和介质有强烈的传热，燃烧温度一般较高，能量传递方式主要以辐射传递为主，此时辐射换热约占总换热量的70%-90%[1]。
