02材料科学与工程中典型物理场的数值模拟.pptVIP

第二章 材料科学与工程中典型物理场的数值模拟;2.1 温度场的计算;热量从一个物体传至另一个物体，或由同一物体的这一部分传至另一部分的过程称为传热或换热。 只要两物体间或同一物体内部存在温度差，就会发生热量的传递，而且热量总是从高温部分向低温部分传递。;按物体温度是否随时间的变化而变化，热量传递过程可分为稳态过程（又称定常过程）与非稳态过程（又称非定常过程）两大类。 凡是物体中各点温度不随时间改变而变化的热传递过程称为稳态热传递过程，反之称为非稳态热传递过程。;物体传热的三种基本形式 导热 对流 辐射;物体各部分之间不发生相对位移，依靠分子、原子及自由电子等微观粒子的热运动进行的热量传递称为导热。 导热现象的规律用傅立叶定律描述。;一.导热方程 ;2. 热流密度（单位面积热流量） 单位时间内通过单位面积的热量称为热流密度，记为 q （W/m2） 。;3. 对流;牛顿对流换热计算公式：对流换热所传递的热流量与流体和固体表面间的温度差以及两者的接触面积成正比。对流换热的基本计算式是牛顿冷却公式： 流体被加热时：q＝α（tw-tf） 流体被冷却时： q＝α（tf-tw） 式中， α 称为对流换热系数，简称换热系数;4. 热辐射;黑体单位时间内发射出的辐射热流量用斯蒂芬一玻尔兹曼定律表示： Q＝Fσ0T4 T为黑体的绝对温度；σ0为黑体的辐射常数。 实际物体辐射热流量的计算采用斯蒂芬玻尔兹曼定律的经验修正形式： Q＝εFσ0T4 ε称为该物体的黑度（又称发射率）。;二、导热微分方程;导热的微分方程的建立;导入及导出微元体的总热流量可用傅立叶定律推出。任意方向的热流量总可以分解成为三个坐标轴方向的分热流量，上图中标出了微元的分热流量。 根据傅立叶定律，通过x，y，z三个表面导入微元体的热量可直接写出：;同理，通过x+dx，y+dy，z+dz导出微元的热流量为：;微元体内的增量为： 微元体内的生成热为： 由上述的式子可以得到：;一般的三维导热微分方程（非稳态）;等式左边：体元升温需要的热量； 等式右边：前三项分别是x，y，z方向流入体元的热量；最后一项是体元内???源产生的热量 微分方程的物理意义： 体元升温所需的热量应等于流入体元的热量与体元内产生热量的总和。;二维热传导微分方程;定解条件是指使导热微分方程获得适合某一特定导热问题的求解的附加条件。;三、初始条件与边界条件 1、初始条件：所求问题的初始温度场，也就是零时刻的温度场。 均匀温度场：T|t=0=T0 温度场不均匀： T|t=0=T0(x, y, z);2、边界条件 （1）第一类边界条件 T=Tw T=Tw(x, y, z，t);（2）第二类边界条件 q= q= （3）第三类边界条件;四、导热问题的数值解法;*泰勒级数展开法;向前差商： 向后差商： 中心差商：;将上两式相加可得;例1：以一维稳态导热为例来说明有限差分方法的应用。设有一炉墙，厚度为δ ，已知内壁温t0＝500℃，外壁温度为t4＝100℃，求厚度方向上的温度分布。 解：该问题可用下式求解： ;t=c1x+c2 x= 0, t=c2=500 x= δ,t=c1δ +c2=100 c1=100-500/ δ=-400/ δ 即得：t=c1x+c2= （-400x/ δ）＋500 因此在特定点上的值为：;采用有限差分求解： 如图 所示，把需 要求解的空间区域 0 ～ δ 以间距Δx（称为步长）划分m等分，这些等分线称为网格线。;为了计算节点温度，需要建立节点的差分方程。 将差分方程代入得：ti+1-2ti+ti-1=0;设m＝4，则： t0＝500 t2－2t1+t0＝0 t3－2t2+t1＝0 t4－2t3+t2＝0 t4＝0 利用高斯－赛德尔法可求解得到各节点温度值。;平面温度场的有限差分求解;（2）差分方程的建立;?根据导热问题的控制方程 ( 导热微分方程 );边界条件的差分格式： 1）对流传热边界条件 2）热流边界条件： 3）绝热边界条件： 4）给定温度边界条件：;方程组：;例2：设有如图所示的温度区域，边界点温度已知，且区域内无内热源，求内部节点温度。 已知T1＝ T2 ＝T3＝ T4＝100℃， T5＝860℃， T6＝400℃，T7＝ 300℃，T8＝200℃。;解：根据;雅可比叠代法：;非稳态导热问题;考察一个无内热源，导热系数为常数，宽、高无限，厚为L的平板的一维非稳态导热问题。当无内热源，λ，ρ，cp 均为常数时，