有限差分法解热传导问题.ppt

有限差分法解决热传导问题,傅里叶定律：在导热现象中，单位时间内通过给定截面的热量，正比例于垂直于该界面方向上的温度变化率和截面面积，而热量传递的方向则与温度升高的方向相反。,基本概念,导热微分方程：,基本概念,定解条件：使微分方程获得适合某一问题的的解的附加条件。,基本概念,边界条件：,基本概念,NDSolveDux,t,t=Dux,t,x,2,ux,0=x*(1-x),u0,t=0,u1,t=0,u,x,0,1,t,0,0.3,NDSolveDux,t,t=Dux,t,x,2+x,ux,0=x*(1-x),u0,t=0,u1,t=0,u,x,0,1,t,0,0.3,NDSolveDux,t,t=Dux,t,x,2,ux,0=Sinx*Sinx,Derivative1,0u0,t=0,Derivative1,0uPi,t=0,u,x,0,Pi,t,0,1,NDSolveDux,t,t=Dux,t,x,2,ux,0=-x*(1-x),Derivative1,0u0,t=0,Derivative1,0u1,t=3-u1,t,u,x,0,1,t,0,0.3,传热问题的数学求解,传热问题的数学求解,泰勒级数展开法：,建立离散方程,建立离散方程,若取上面式右边的前三项，并将两式相加移项整理即二阶导数的中心差分： 同样可得：,建立离散方程,对于二维稳态导热问题，在直角坐标中，其导热 微分方程为： 其节点方程为：,热平衡法,建立离散方程,基本思想：对每个有限大小的控制容积应用能量守恒，从而获得温度场的代数方程组，它从基本物理现象和基本定律出发，不必事先建立控制方程，依据能量守恒和Fourier导热定律即可。 能量守恒：流入控制体的总热流量控制体内热源生成热 流出控制体的总热流量控制体内能的增量,建立离散方程,内部节点：,建立离散方程,17,1.边界节点离散方程的建立：,(1) 平直边界上的节点,18,(2) 外部角点,19,(3) 内部角点,写出所有内节点和边界节点的温度差分方程 n个未知节点温度，n个代数方程式：,GaussSeidel迭代,GaussSeidel迭代,GaussSeidel迭代,200,100,100,Tf=0 k=1W/(m*K) h=10W/(m2* ) (qw=h*(Tw-Tf),t1,t2,t3,t18,t17,t16,Ax=b,b=300,200,200,300,100,0,0,100,100,0,0,100,100,0,0,100,0,0 x=t1,t2,t3,.t18,偏微分方程工具箱,Step 1 “Draw模式”绘制平面有界区域 ，通过公式把Matlab系统提供的实体模型：矩形、圆、椭圆和多边形，组合起来，生成需要的平面区域. Step 2 “Boundary模式”定义边界，声明不同边界段的边界条件. Step 3 “PDE模式”定义偏微分方程，确定方程类型和方程系数c,a,f,d，根据具体情况，还可以在不同子区域声明不同系数. Step 4 “Mesh模式”网格化区域 ，可以控制自动生成网格的参数，对生成的网格进行多次细化，使网格分割更细更合理. Step 5 “Solve模式”解偏微分方程，对于椭圆型方程可以激活并控制非线性自适应解题器来处理非线性方程；对于抛物线型方程和双曲型方程，设置初始边界条件后可以求出给定时刻t的解；对于特征值问题，可以求出给定区间上的特征值.求解完成后，可以返回到Step 4，对网格进一步细化，进行再次求解. Step 6 “View模式