用于备向异性材料电磁问题中的一种新时域有限差分算法.pdfVIP

电工数学进晨 =a扛直垒殛坠Q签垫幽￡丛鲴墅g型竖 !塑 用于各向异性材料电磁问题中的一种新时域有限差分算法 李永明俞集辉 ‘(重庆大学电气工程学劈已重庆400044) 摘要 本文针对时麓青最差分法在传统的‰氏一格壹■求髂向异性材料电硅同曩中一近慨通多的情 况，采用了利用一种新的_I掐变量法．并利用壹舫程的对■性推导了这种膏幡■式下塌量的计算格 式． 美t膏t时麓：一—t分涪Y∞氏罔格对儡性 1前言 时域有限差分法首先由kY∞作为一种求解麦克斯韦旋度方程的数值计算方法而提出的， 它的特点是直接从麦克斯韦方程组中的旋度方程出发，采用所谓的Y∞氏同格结构(如圈1)， 在空间和时间上进行离散，然后用差分来近似瞢代麦克斯韦方程的旋度方程．但这种F¨D法 应用封计算各向异性材料的电磁过程同惠中时，其计算公式表述十分繁杂，近戗取得过多，这 无形中兢降低了计算精度．为此，我们引入一种囊的差分嘲格模型，并将其朗翻鲁向异性材 料中，得到一种新的、更好的有限时域差分算法． 面 ■I k氏t骨一格蔓■ ●2奉土所tt舟一幡童t 国衷自然科学基金赍助项目及博士学科点基金赍助项目 ．49 ‘ 2差分网格模型 要差分麦克斯韦旋度方程，建立正确的差分算法，首先必须建立适当的物理模型，以保证 电场和磁场在时空上是相关的。 考虑图2的离散网格模式，我们把电场的各分量置于网格的整同格点上，磁场的各分量置 于网格的半网格点上。这样的电磁场空间配置使得电场分量和磁场分量在空间中是相互环绕的， 因而符合电磁场的基本规律一Faraday电磁感应定律和^呻ere环流定律，亦即麦克斯韦方程组 的基本要求。下面我们根据这种差分网格模式来推导各向异性媒质材料中的场量计算公式． 3．应用于各向异性媒质中的新FD_TD法 ． 对于一般任意各向异性线性媒质，基本方程组的徽分形式是： ‘vx茸：一尝一歹． (1I) v。厅：罂+，． (1” 。 口f 上，厶分另唾表示电流密度，礁流密度． 其中，本构关系用矩阵形式的媒质参数表示为(s为鞋阻率矩阵) 西=蠢。Z=荔，豆=面，jlI=荫 利用上述的电磁参量表达式，展开基本方程组中的式(1a)并整理得 堡一生 ll az ay (2) ，