《计算物理(本科)》[第章].pptVIP

第九章 热传导方程的数值解法 一、传热过程：在t 时刻， dt 时间内通过dS截面的热量 §9.2 一维热传导方程的差分法 第三类边界条件 §9.3 二维热传导方程的差分法 二维热传导方程求解与一维的类似，在确定了差分格式和给出定解条件后，按时间序号分层计算，每一层由二维点阵组成，通常称为网格。 三、差分方法求解 建立差分格式 设空间步长为h，时间步长为τ。如图所示，将xoy平面均分成N×M的网格，并使 Nh=l Mh=s 对节点(i，j)，在k时刻(即t=kτ)的差分式 恒温热源边界：红色部分，即 {x=0, y ∈[M1h，M2h]} 归一化后高温源温度取“1”，即 §9.3 MATLAB实现热传导方程数值解法 实验八 热传导方程的差分解法 一、实验目的 1．了解一维热传导方程的差分解法； 2．掌握二维热传导方程的差分解法。 二、实验内容 1．建立一维热传导方程差分格式 2．建立二维热传导方程差分格式 3．在Matlab中，求解热传导方程函数 合肥工业大学电子科学与应用物理学院 * §9.1 热传导方程的引入 物理学中对热传导现象和扩散现象等物理过程的描述，都归结为的二阶偏微分方程，统称为热传导方程，也称为抛物线型方程。 本章以热传导现象进行介绍： 当物体内部各处温度不同时，就会有热量从温度较高处向温度较低处传递，这就是热传导现象。 其中K为热传导系数；　 是在ΔS面的法向上温度梯度。 对上式进行闭曲面积分，并用矢量积分定理得 其中V为闭曲面包括的体积；t1,t2表示t2~t1时间段的 积分； 三、闭体能量守恒：若体积v内无热源，热传导系数K是常数时，根据能量守恒Q=Q’得三维齐次热传导方程 。 二、体内耗热过程：设物体的比热容为c，密度为ρ，则v内温度变化所消耗的热量为 一、热传导方程：物体内无热源的一维热传导方程为 二、初始条件（t=0的温度分布情况）： 三、边界条件（物体边界的温度分布情况）： 第一类边界条件 ; 第二类边界条件 其中h1(t)、h2(t)、g1(t)、g2(t)为给定的函数，h1(t)≥0、h2(t)≥0且不同时为0。 并使 Nh=l Mτ=T 则有 t=kτ k=0,1,2,…,M (时间步序号) x=i h i=0,1,2,…,N (空间x方向的步序号) 四、差分方法求解 建立差分格式 通常分别取空间步长和时间步长为h和τ，均为常数；计算时的步序号空间用i表示，时间用k表示。 定义：一阶向前差商近似 一阶向后差商近似 二阶中心差商近似 用二阶中心差商作为二阶微商近似是 对时间的一阶差商近似是 由上两式，并令α=τλ/h2，得一维差分格式 其中N=[l/h],M=[T/τ ],“[ ]”表示取整。 定解条件 几何关系 把某一时刻计算的各点称为一层，计算过程是一层层地进行的，各点关系如图。由图可见，k+1时刻i点的值是由k时刻的i-1、i和i+1三点的值推算出来的。 t x k+1 k i+1 i i-1 i≠0,N 为什么? 可以证明，上式的差分格式收敛稳定条件为 α=τλ/h2≤1/2 五、归纳上述一维热传导的差分格式为 六、差分格式计算步骤 1.给定λ，l，h，α，T ； 2.计算N=[l/h]，M=[T/τ ] ，τ=αh2/λ； 计算边界值： 3.计算初始值： 4.差分格式的计算 一、热传导方程：均匀物体内部无热源的二维热传导方程为 二、初始条件 则有 t=kτ k=0,1,2,… (时间t序号) x=ih i=0,1,…,N (空间x序号) y=jh j=0,1,…,M (空间y序号) y x M N M2 M1 将上式代入二维热传导方程中，并令α=τλ/h2即可得二维热传导的差分格式 i=1,2,…,N-1； j=1,2,…,M-1。(注意： i≠0,N；j ≠0,M ) 几何关系 用上式的二维热传导的差分格式，计算k+1时刻点(i,j)处的温度值时，使用了k时刻点(i,j)及其相邻的四个点的温度值，如图所示相关各点之间的关系。 t x k+1 k i+1 i i-1 j+1 j-1 j 四、边界条件(如图所示的具体问题) 绝热边界：粉红色部分绝热壁，{x=0,y∈(0,M1h).a