基于非结构Cartesian网格电磁散射场计算.pdfVIP

清华大学学报(自然科学版)2007年第47眷第1l期 33／36 !墅塑!!!Q：!15± cN 11—2223／NJTsi“gh岫U面v(scl＆Tcch)t2007，V01．47，No．11 Z068_2071 基于非结构cartesian网格的电磁散射场计算 张兆1，沈孟育1，张涵信2 (1．请华大学航天航空学院，北京1000HI2．中国空气动力研究与发屉中心，绵阳621000) 摘要：计算洼体力学中的网格技术和数值方法对于计算 时域中的电磁散射场求解是计算电磁学 电磁学的应用具有借鉴意义．该文结台计算流体力学中的披 (cEM)中十分重要的研究领域。时域电磁场的计算 值方法．建立了基于非结构carteajan阿椿上的时蛾有限体 主要有有限差分和有限体积2种方法“】。Yeeoo提出 积法．用于末解电磁散射问题．为保证这种方法在时问和空 的时域有限差分法，采用二阶中心差分格式在电磁 何上具有二阶精度，在时问离散上采用二步Runge—Kutta方 场矢量时空交错的cartesian阿格上求解，此后 法，空J可高教采用非结构坤二阶NND格式。通过末解时域 HoIland01进一步研究了非正交曲线坐标系中的有 中的Maxwen方程组．计算了二维翼型和三维秉波体外形 限差分法。由于Maxwell方程组和流体力学中 等典型完垒导电散射悻的雷遮截面．甘算i占皋表明：该方法 计算精度高．适用于复杂外形的电磁散射问题的隶解． 关键词：计算电硅学l非结构cane8hn网格。曹迭藏面I时 限体积方法应用于电磁散射和辐射问题的计算， 域有限体积法INND格武 Sha“g嘲还把高阶紧致格式应用于求解电磁传播 中圉分类号：0441 文献标识码；A 问题。 文章编号：1000—0054(2007)11—206B_04 技术。吲移植到电磁散射场求解上来。为了克服 ComputatioⅡof scattering electromagneticcartesian网格在壁面边界带来的阶梯效应，在壁面 fieldbasedonunstructuredCartesian grids内部添加虚单元采用镜象边界条件。远场边界上采 刮^NG 用特征方向上的波动方程作为吸收边界条件。方程 Zh∞l，SHENM鲫gyul．ZH^N0№nxi巾 (1．鼬捌0，如瑚哗“nIn曲岫u耐Ⅷ竹， 在空间上采用cFD中的非结构NND格式，时间上 酬●qJ∞O¨，C‰， 的推进采用二阶RuⅡge—Kutta法。 2．CMn^州y山∞R—岫●州Ik代l叩Ⅲ阻t(：呲悻， MhⅡp“g6210∞-伽皿) l 非结构的cartesian网格 ^蛔h‘-n．AEIik—volumetime一蜘aInm吐h甜b蚰edon u珊trIlctu捌C缸幅hn旷ldsw衄de垤lopedu3i雌∞m叫tatb鼬l num t0 中，这种网格对于复杂外形具有良好的适应型。一般 mHhodBsolvethc dy蚰皿cB 8cBt№rmgp10bl锄of coTImuta