微波加热物体的数值仿真.ppt

微波加热物体的数值仿真,2011年6月23日,工作内容,二,物体在非正弦波激励下升温过程的计算,一,物体在正弦波激励下升温过程的计算,三,半导体器件发热的计算,四,NND有限元求解一维欧拉方程,目的: 解决电磁场、热场耦合问题（电磁场耗散能量使物体升温，温升改变物体的介电常数等性质，从而再影响电磁场分布）；给出物体的时域温度变化过程。 核心问题: 电磁场变化快（ps 数量级）与热场变化相对慢( ms-s数量级）的交互影响耦合问题。,一,物体在正弦波激励下升温过程的计算,问题描述,电磁场迭代达到相对稳态之后，利用热损耗公式计算出电磁损耗，并将其作为热源代入热模型进行热问题求解，若干时间步后更新物性参数（介电系数、电导率等），重新计算电磁场，如此循环直至预定加热时间。基本流程如下：,计算方法,算例1：微波加热水的稳态温度分布,C. A. Vriezinga, Thermal profiles and thermal runaway in microwave heated slabs, J. Appl. Phys. vol.85,no.7,pp.37743779,April.1999.,被加热物体是置于空气中无限大 水片,入射波为频率2.45GHz的 平面波。,算例验证,Tanmay Basak , Badri S. Rao，Theoretical analysis on pulsed microwave heating of pork meat supported on ceramic plate, Meat Science,86(2010),pp780-793,Aug,2010.,算例2：微波加热猪肉,猪肉厚度1.6cm，右边是 0.2cm厚的陶瓷片。两侧 入射波均为1.5W/cm2 的 平面波,频率2.45GHz。 加热时间194s。,算例验证,加热过程中猪肉上温度差的变化,Tanmay Basak , Badri S. Rao,Theoretical analysis on pulsed microwave heating of pork meat supported on ceramic plate, Meat Science,86(2010),pp780-793,Aug,2010.,算例2：微波加热猪肉,算例验证,P. Rattanadecho,N. Suwannapum,Interactions Between Electromagnetic and Thermal Fields in Microwave Heating of Hardened Type I-Cement Paste Using a Rectangular Waveguide,Journal of Heat Transfer,Vol.131,pp.082101-1-082101-12,Aug.2009.,算例验证,算例3：微波加热水泥浆,文献结果,FDTD-HTE,加热20秒后水泥浆纵切面温度分布,入射波为2.45GHz,1000W,材料尺寸 为(xyz)110mmx55mmx80mm.加热 时间20s.,前面介绍了正弦电磁场作用下物体热效应的分析方法,当电磁场不是正弦形式时,这种方法并不适用.这时需要用电磁场的瞬时耗散功率分析热效应.,二,物体在非正弦波激励下升温过程的计算,计算流程,脉冲作用下互联线最高温度变化,史炎冰,多层高密度互连在静电脉冲作用下的时域非线性有限元电热耦合模拟,上海交通大学硕士论文,2009.,互联线温度变化,算例验证,90nm工艺三层互联线. 脉冲作用与顶层中间的互 联线上.脉冲形式如下:,三,半导体器件发热的计算,与普通物体温度场计算不同,半导体器件温度场计算有 以下几个特点. 1,热源一般根据器件内部电流,电压直接计算. 2,当只关心稳态温度时,可以只在空间上迭代求解.,上图为SOI MOSFET的三维模型. 热源位于栅极下方,P=1mW/um.,算例验证,算例1：SOI MOSFET的稳态温度分布,Kun Zhang, M.C.Cheng,Thermal Circuit for SOI MOSFET Structure Accounting for Nonisothermal Effects, IEEE Trans. Electron Devices, vol. 57, no. 11, pp.2838-2847,Nov. 2010.,算例验证,算例1：SOI MOSFET的稳态温度分布,xoy面温度分布-FEM,xoy面温度分布- FDTD-THE,Kun Zhang, M.C.Cheng,Thermal Circuit for SOI MOSFET Structure Accounting for Nonisothermal Effects, IEEE Trans. Electron Devices, vol. 57, no. 11, pp.2838-2847,Nov. 2010.,使用瞬时功率计算方法计算正弦波加热问题,如果采用瞬时损耗功率直接计算正弦波加热问题,由于实际加热时间约为电磁场和热场交替计算时间步的1012倍,计算耗时长到无法接受. 采用alpha因子,人为加快热场的变化过程后,可以将实际加热时间缩短 (算例中可以缩短10-9s数量级),从而大大缩短计算时间.,alpha=8.5e9时的升温曲线,算例验证,算例1：使用瞬时功率计算方法计算正弦波加热问题,Hx,Hz,矩形波导TE10入射，Pin=1000W, f=2.45GHz, 加热甲基纤维素，复介 电常数虚部：,S. Curet, O. Rouaud, L. Boillereaux, Microwave tempering and heating in a single-mode cavity: Numerical and experimental investigations, Chemical Engineering and Processing 47 (2008) 16561665.,alpha=8.5e8时的升温曲线,算例验证,算例1：使用瞬时功率计算方法计算正弦波加热问题,S. Curet, O. Rouaud, L. Boillereaux, Microwave tempering and heating in a single-mode cavity: Numerical and experimental investigations, Chemical Engineering and Processing 47 (2008) 16561665.,一维欧拉方程,NND有限元解一维欧拉方程,前面三个方程可以记为,记,NND有限元解一维欧拉方程,则各个方程可以离散为（以第一个质量守恒方程为例）,表示为矩阵形式：,NND有限元解一维欧拉方程,对线积分（三维情况是面积分）特殊处理1,以第K单元为例,K,每个单元端点处的f都可以根据流场的特征量分解为正负通量之和。正通量反映上游对下游的影响，负通量反映下游对上游的影响。所以取单元左侧（上游）的正通量和单元右侧（下游）的负通量之和作为f在单元的上的分布。,1王岩，计算流体力学中有限元的强间断处理，空气动力学学报，1994年3月，12卷1期。,于是线积分项为,NND有限元解一维欧拉方程,Steger-Warming通量分裂法,Jacobi矩阵：,A的特征值为：,则Steger-Warming通量分裂为,算例：一维sod问题,下图表示一个细长的激波管。t=0时，管子中心有一隔板，此时激波管中压力，密度速度的分布如图所示。,隔板,t=0时刻，突然撤去隔板，则在t=0.5时刻，激波管中密度，压力， 速度的分布分别为：,X=-0.5,X=0.5,贺立新，间断Galerkin有限元方法及其与有限体积混合计算方法研究，中国空气动力 研究与发展中心研究生部博士论文，2008年8月。,算例：一维sod问题,算例：一维sod问题,算例：一维sod问题,不连续Galerkin有限元方法(DGM),DGM的优点,DGM与传统有限元方法的本质区别：传统有限元建立在整个求解区域加权余量 为0的基础上，因此需要形成并求解全局稀疏矩阵，DGM则建立在单元加权余量为0 的基础上，只需求解单元矩阵。因此大大提高的技术速度。 DGM中不同单元之间仅通过数值通量联系，各个单元可以根据需要采用不同类型、 不同阶数的基函数，有利于采用自适应网格，且高度并行。 DGM不再要求单元边界连续，单元边界左右的值可以不同，不同单元所对应的同 一节点值可以不同，这在计算含有间断解的流体力学中有天然优势。,由于不同单元交界面处的场值可以不同，需要分别记录，相对传统有限元消耗 内存更多。,1龚小权，欧拉方程间断Galerkin有限元方法初步研究，中国空气动力研究与发展中心研究生部硕士论文,2009年。,DGM的缺点,不连续Galerkin有限元方法(DGM),三维NS方程可记为如下形式,在单元 上对其进行Galerkin测试,其中,是单元,的表面,表示面,上的数值通量,不连续Galerkin有限元方法(DGM),1李宏，高分辨率间断有限元方法，计算物理，第21卷第4期，2004年7月。,求解 的方法有Godunov, LaxFriedrichs, Roe,EngquistOsher, 和HLLE等,但常用的是Lax-Friedrichs方法1,2.,2 F. Bassi, S. Rebay,A High-Or