含时Mazwell方程的辛方法分析.pdf

含时Maxwell方程的辛方法分析木 戴桂冬1一，唐贻发1’ 1中国科学院数学与系统科学研究院计算数学所，北京， 100080 2中国科学院研究生院，北京， 100080 擒要 本文给出了含时Maxwell方程的Hamilton表示；发现空间的有限元离散导 出的半离散方程是Poisson系统或带扰动的Pomson系统；证明了对该Foisson 系统隐含的辛结构． 自上个世纪60年代起，计算电磁学兴起并成为研究的热点．近些年，时 中包括有限差分u,-Jl盲-J区域(FDTD)方法【i，2】，时间区域有限元方法【3H6】， 法，它的优点是在时空处理上都比较便利，但是由于有限元方法在处理空间区 域尤其是不规则区域问题上具有优越性，因此后续的计算工作大多用有限元 处理空间．这方面工作的思路一般是：先利用磁场和电场的关系，消掉一个变 量，得到只关于电场E或磁场日的方程；然后先将时间进行离散，得到关于 空间的偏微分方程；最后用有限元方法求解． 们的方法属于FDTD方法，即时空都是用有限差分格式，只是时间上的离散 ‘本工作得到国家自然科学基金(No和中国科学院知识创新工程信息化建设专项”超级 计算环境建设与应用”】哽目科学计算平台建设应用课题(INFl05-SCF_,-02-20)资助． ’通讯作者，Email：tyfQtsec．cc．跹蕊． 程，比较理想地模拟了谐振腔中的共振波【11】．苏在她的博士论文中给出了 程的Hamilton表示形式，及其空间离散后半离散方程是否具有辛结构．我们 的目标是利用电场E或磁场日的耦合关系来建立对称的Hamilton方程，并 且希望有比较充分的理由在时间上能用辛方法离散． 本文指出在同一介质中介电系数和磁导率是常数的情况下，含时Maxwell 方程有Hamilton表示；发现MaweU方程空间的有限元离散导出的半离散方程 构． §2 Maxwell方程及其半离散形式 介质中的Maxwell方程为 雕_ (2．1) B为磁感应强度 其中变量对应的物理意义如下- E为电场强度(V·m以)， 流强度(A．m一2)，P表电荷密度(C·m．3)． 另外由介质的电磁性质得到以下两个关系式： ㈠‘去oE一+MP 仁2， 其中eo表示真空中的介电系数，／Jo表真空中的磁导率，P是介质的极化强 度，M是介质的磁化强度． {量三peE日 江3， 其中c是介质的介电常数，p是介质的磁导率，它们还可以通过eo和Po用相 对介电常数￡r和相对磁导率脚表示为：e=creo，p=坶瑚． 本文中我们只考虑线性、各向同性的介质中变化的电磁场问题． §2．1 Maxwell方程的Hamilton表示 裳+玑扣o 和(2．1)中的方程并不完全独立，因此要求解电场和磁场，我们只需用到含时 方程组(2．1)中的两个含时间偏导数的方程．利用本构关系(2．3)消掉E，D，B，日 中的两个变量，得到只关于E和B或者E和日的方程组． 比较简单的情形是介电常数e和磁导率弘都是与时间空间无关的常数， 这时，Maxwell方程组可以表示为无穷维Hamilton形式； 腾=壶V—J=警 I等=一V×E一等． x dy B-{B·J 何nm陋，引=上丢E·V×E+土2弘B．V 实质上，上面的方程能写成Hamiltion形式是因为具有某种对称性．将它 其中z=【E丁，B丁n冗=阿，jr]r，03表示三维零列向量，