微波加热中热失控的一维数值模拟研究.pdfVIP

2005’全国微波毫米波会议论文集 微波加热中“热失控’’的一维数值模拟 刘长军 闫丽萍黄卡玛 四川I大学电子信息学院，成都(610064) 摘要：本文对微波加热陶瓷过程进行了数值模拟，使用时域有限差分的进行电磁场 分布的计算，耦合热传导方程计算微波加热过程中的温度变化。通过计算不同入射微波 功率对应的介质的稳态温度，对微波加热过程中的热失控现象进行了初步分析。 关键词：时域有限差分，Biot系数，热传导 微波加热技术在工业、化工、食品、材料等领域都用广泛的使用。微波加热过程中， 一些材料的介质特性随温度发生变化，导致一些特殊的宏观热现象。例如：在微波焊接 陶瓷、微波加热矿物、微波加热橡胶等过程中，出现材料温度持续升高直至烧毁的现象。 热失控现象阻碍了微波加热的产业化应用。 热失控是由于电磁波与物质的相互作用导致的一种非稳定的热现象。例如：在微波 加热介质的过程中，如果某个局部的热传递速率低于微波功率损耗速率，将会导致局部 温度的升高和介质特性的变化。如果介质特性的这种变化加速了介质局部对微波功率的 损耗，从而由于“正反馈”过程形成了热点，这就是一个局部的热失控现象。在一些条 件下，局部的热失控现象会体现为宏观的非稳定热现象。 使用时域有限差分算法耦合热传导方程对一维微波加热陶瓷介质中的热失控现象 进行数值模拟。在Biot系数较小的条件下，模拟了稳态温度与微波功率之间的不连续变 化。使用本方法可望对更为复杂的模型进行数值仿真，促进微波加热的研究。 1微波加热模型 微波加热过程中电磁场满足Maxwell方程，如方程(1)所示。 × l| Vl『E × II V t『H 况一扣．西一扭． 6E 其中。为电导率。在一些应用中使用TEM模式进行微波加热，可以适当简化为一维问题。 热传导方程为 (2) pq掣=K可02T(z)+易(z) 其中T为温度，p为介质的密度，cp为比热系数，K为热传导系数，Pd为微波耗散功率。 如果忽略介质的极化损耗，微波耗散功率可以用电场强度表示为 1040 2005’全国微波毫米波会议论文集 B(z)=去仃lE(z)12 (3) 假设陶瓷介质周围的温度为To，保持不变。热传导方程的边界条件为 (4) KiOr(z)=办[丁(z)一乃]慨[丁4(z)一万] 其中h为热对流常数，S和e分别为热辐射常数和表面辐射率。 实验研究发现，陶瓷的电导率随着温度的升高而增大，并且呈现出非线性的关系。 这种非线关系一般可以使用指数函数模型或者Arranius模型进行描述。如果采用指数函 数模型，陶瓷电导率与温度的关系表示为 巾)w冲卜孚I (5) 其中60为温度为T。时陶瓷介质的电导率，k为材料相关的常数。在陶瓷加热的热失控研 究中，通常使用Biot系数描述介质尺度与热参数之间的关系。Biot系数定义为 B：一lh f6) K 其中，为介质的尺