有耗色散介质中的FDTD方法..doc

有耗色散介质中的FDTD方法 李庆伟(核工业西南勘察设计研究院有限公司 四川成都 610059) 摘 要：讨论了时域有限差分(FDTD) 法在地下有耗色散介质中的应用。用数值模拟地下有耗色散目标的瞬态响应,并对结果进行了分析. 从以上工作可看出该方法在有耗色散介质中是有效可行的。 关键词：时域有限差分 有耗色散介质 Debye方程 吸收边界条件 STWBC边界 作者简介：李庆伟 (1975.4-)，男，工程师，当前主要从事物探及工程检测工作 一 引言 用时域有限差分法（FDTD）分析与频率无关的非色散的复杂电磁场问题具有显著的优越性.对单色稳定电磁波的问题，媒质的色性质不产生影响，对窄频带的问题也可以近似地忽略媒质的色散特性，但对包含宽频谱的瞬态电磁场问题媒质的色散性质起着作用．uebbers等-1990）和直接微分法。由于很多色散媒质的极化率都无法表示为指数形式，因此无法应用这种递推卷积方法，另外，即使极化率可以表示成指数形式，但每种色散媒质都必须单独推导自己的FDTD方程.显然这种递推法的通用性不强，使用极其不便.直接微分法则依据本构关系的时域微分方程，建立辅助微分方程，最终通过求解关于矢量H,D,E的时域差分方程组得到场解。本文拟从直接微分法出发探讨FDTD法在有耗色散介质中的应用。 二 实现方法 2.1 基本方程式 无源有耗(不考虑磁损耗，即)色散媒质，时域内的Maxwell方程为 （1） 其中，E-电场强度(V/m)；H-磁场强度(A/m)；B-磁通密度(Wb/m2)，D-电通密度(C/m2)；σ-电导率(S/m);μ-磁导率（H/m），这里设磁导率与频率无关，即。 以二维TM波为例，采用Yee氏网络，由(1)式中第一式可得到H各分量在t=n△t时刻随时间推进的FDTD公式： （2） 其中，为每一计算区域的μ值，△s为等距空间步长。 将(1)式中第二式在t=(n+1/2)△t时刻做差分离散，可得 即 （3） 按照Yee氏网络，可得Dz分量的FDTD公式： （4） 2.2 转换公式 下面通过直接微分法实现由D到E的转换[1]： 设介质为一阶Debye型有耗色散介质，可以写成以下形式 （5） 式中，为复介电系数，为真空介电常数，分别为静态和无限频率的相对介电常数，σ为电导率,τ0为驰豫时间常数。 将（5）式代入色散介质频域本构关系，可得 （6） 即 （7） 利用频域到时域的算子转换关系，可得（7）式的时域方程为 （8） 将（8）式在t=(n+1/2)△t时刻离散，可得 （9） 即 （10） 由此可以得到时域中由D到E转换的FDTD公式 （11） 其中， ，，， 在以上各式中，对二维TM波，考虑解的稳定性及FDTD法的收敛条件，可取 综上所述，有耗色散介质中FDTD随时间推进的的步骤为，通过式（2）、（4）、（11），依次计算出各分量。 三 吸收边界条件 3.1 吸收边界条件[2] 用FDTD求解电磁场问题时假定空间是无限大的，受计算机存储量和计算时间的限制, FDTD的计算只能在有限区域进行。为了使框定的有限空间与无限空间等效，在计算区域的截断边界处必须给出吸收边界条件.使得向边界面行进的波在边界处无明显的反射现象。吸收边界条件的选取与合理设置是FDTD 法计算中必须考虑的重要内容. 吸收边界条件从开始简单的插值边界，到后来较常用的Engqusit-Majda的吸收边界条件(Engqusit and majda,1977),Mur吸收边界(Mur,1981)，Lao吸收边界条件(Liao et al.,1984)以及近年发展的PML完全匹配层（Berenger,1994）和UPML各向异性介质完全匹配层(Sacks 1995,Genney 1996).其吸收效果越来越好。 3.2 驻波-行波吸收边界条件(STWBC) [3] PML边界体现了更大的优越性,其缺点在于复杂性和适应性,如果是研究有耗问题或色散问题,PML 的应用困难将会大大增加. Tan 于2001 年提出的驻波-行波吸收边界条件(STWBC)，由于该边界比单向波边界所需计算空间小,数值稳定性好,同时不象PML边界需要进行场量分离和附加额外的吸收层,因此计算效率较高. 驻波-行波边界是在计算区域边界处加理想导体,根据电磁场理论,波到达边界面将发生全反射,若是理想导电(磁)壁,则切向电(磁) 场为零,切向磁(电) 场是入射场的两倍,这是严格正确的.同时,反射场将往回传播,在区域内部形成驻波,随时间推移,驻