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

第二章 材料科学与工程中典型物理场的数值模拟 2.1 温度场的计算 材料科学与工程领域存在大量传热问题 金属材料的热加工 高分子材料的成型 陶瓷材料的成型、烧结 …… 热量从一个物体传至另一个物体，或由同一物体的这一部分传至另一部分的过程称为传热或换热。 只要两物体间或同一物体内部存在温度差，就会发生热量的传递，而且热量总是从高温部分向低温部分传递。 按物体温度是否随时间的变化而变化，热量传递过程可分为稳态过程（又称定常过程）与非稳态过程（又称非定常过程）两大类。 凡是物体中各点温度不随时间改变而变化的热传递过程称为稳态热传递过程，反之称为非稳态热传递过程。 物体传热的三种基本形式 导热 对流 辐射 物体各部分之间不发生相对位移，依靠分子、原子及自由电子等微观粒子的热运动进行的热量传递称为导热。 导热现象的规律用傅立叶定律描述。 一.导热方程 1. 傅立叶（Fourier）定律 如图 3-1 所示的两个表面分别维持均匀恒定温度的平板，对于x方向上任意一个厚度的微元层来说，单位时间内通过该层的导热热量Q（W）与当地的温度变化率及平板面积A成正比，即: 3. 对流 对流是指流体各部分之间发生相对位移，冷热流体相互掺混所引起的热量传递方式。 对流仅发生在流体中，而且必然伴随着导热。 牛顿对流换热计算公式：对流换热所传递的热流量与流体和固体表面间的温度差以及两者的接触面积成正比。对流换热的基本计算式是牛顿冷却公式： 流体被加热时：q＝α（tw-tf） 流体被冷却时： q＝α（tf-tw） 式中， α 称为对流换热系数，简称换热系数 4. 热辐射 物体通过电磁波来传递能量的方式称为辐射。 由于热的原因发出辐射能的现象称为热辐射。自然界中各个物体都不停地向空间发出热辐射，同时又不断地吸收其他物体发出的热辐射。 辐射与吸收过程的综合就造成了以辐射方式进行的物体间能量传递－辐射换热。 黑体单位时间内发射出的辐射热流量用斯蒂芬一玻尔兹曼定律表示： Q＝Fσ0T4 T为黑体的绝对温度；σ0为黑体的辐射常数。 实际物体辐射热流量的计算采用斯蒂芬玻尔兹曼定律的经验修正形式： Q＝εFσ0T4 ε称为该物体的黑度（又称发射率）。 二、导热微分方程 对于一维导热问题，直接对傅立叶定律的表达式进行积分就可获得答案。 多维导热问题的情况则较为复杂，虽然傅立叶定律仍然可以应用，但是还必须解决不同坐标方向间导热公式的互相联系问题。这时，导热问题的数学描写必须针对从物体中分割出来的微元平行六面体进行分析才能得到。这种数学描写称为导热微分方程式。 它的建立，除了依靠傅立叶定律之外，还以能量守恒定律为基础。 导热的微分方程的建立 为了便于推导，把讨论对象局限于常物性（即物性λ，ρ，Cp都是常量）的各向同性的材料。 对如下图所示的微元平行六面体进行导热分析，按照能量守恒定律，微元体的热平衡式可以表示为： （导入微元体的总热流量）＋ （微元体内热源的生成热）＝ （微元体内能的增量）＋（导出微元体的总热流量） 导入及导出微元体的总热流量可用傅立叶定律推出。任意方向的热流量总可以分解成为三个坐标轴方向的分热流量，上图中标出了微元的分热流量。 根据傅立叶定律，通过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）第三类边界条件 四、导热问题的数值解法 数值解法是以离散数学为基础，以计算机为工具的一种求解方法。常用的数值解法有有限差分法和有限元法。 有限差分法是把原来求解物体内随空间和时间连续分布的温度问题，转化为求在时间领域内和空间领域内有限个离散点的温度值问题，再用这些离散点上的温度值去逼近连续的温度分布。 它以差分代替微分，差商代替微商，建立以节点温度为未知