三维随机粗糙面与目标复合电磁散射的FDTD方法用时域-计算物理.PDF

[文章编号】 1001—246x(2007)05一0550．1l 三维随机粗糙面与目标复合电磁散射的FDTD方法 匡 磊，金亚秋 (复旦大学被散射与遥感信息教育部重点实验室．上海200433) [摘 要】 用时域有限差分方法(啪)研究三维周期性延拓的随机粗糙面与上方目标复合电磁散射，用周期 性延拓消除数值计算中截取有限大小粗糙面产生的边缘效应，讨论一个周期单元粗糙面的边长与其相关长度之 间的关系．给出在FDTD方法中向粗糙面加载入射波的方法，建立了粗糙面上单个三维目标的复合散射FDTD计 算模型．数值结果给出粗糙面与目标散射的近场分布．应用近远场变换得到全方位散射角的双站散射系数．比 较了三维与二维散射模型的区别．结果显示当粗糙面上放置目标时．其后向散射显著增强． [关键词】 三维粗糙面与目标；取站电磁散射；时域有限差分法；周期边界条件 【中围分类号】T№儿 【文献标识码】 A O引言 粗糙面与目标复合电磁散射研究在十分广泛的领域中有重要的应用，如复杂地海背景上舰船、飞行器和 其他目标的监测，地表下埋藏物探测等“。1．粗糙面散射研究的一些解析方法通常受其近似条件的约束，复 杂粗糙面与离散目标复合的散射计算存在方法上的困难，随着计算机与计算技术的发展，计算电磁学已成为 知数个数增加使得计算复杂度陡然增大．目前，已有一些工作研究粗糙面下方目标的复合散射模型． John一““运用准物理光学的近似方法和低精度的离散偶极子近似(DDA)方法分别对面散射与体散射进行预 处理，运用加速目标，粗糙面组合的矩阵方程快速收敛的迭代矩量法(IterativeMOM)计算了粗糙面上方有无 目标的后向同极化的差值散射．在求解目标与粗糙面相互耦合部分时，目标不能与粗糙面接触，否则计算会 产生奇异点． 本文提出周期性延拓随机粗糙面，用F啪方法只需计算一个周期单元的粗糙面，得到其近场分布，降 低三维粗糙面散射计算量和内存．然后，建立粗糙面上方放置单个三维目标的双站散射FDTD计算模型，得 到整个时域过程中三维目标和粗糙面相互作用的复合电磁散射．无论目标位于粗糙面上方还是贴在粗糙面 上，都可以数值地计算出其近场(总场和散射场)分布．最后运用近远场变换公式，得到全方位散射角度的随 机粗糙面和单个目标，粗糙面的双站散射系数． 1周期延拓粗糙面的近场散射模型 1．1粗糙面周期单元大小的选取 为模拟无限大粗糙面，文[14]用周期延拓面法代替在入射波上加载窗口函数法以克服有限粗糙面的边 缘效应，本文将其推广到三维情形下．把无限大粗糙面看成在x和y方向上具有周期延拓的结构，利用FDTD 方法，通过对一个周期单元粗糙面的计算，得到整个粗糙面的电磁散射，从而减少计算量和内存．如图1所 示，￡；和t分别为粗糙面在*和y方向上的周期长度． 【收稿日期]20嘶一05一10；[修回日期】2006一ll一∞ 【基金项目】国家重点基础研究计划项目(20Dl邙姗400)及国家自然科学基金6∞71050)资助项目 【作者简介】匡磊(1啪一)，男，安徽合肥，博士生．从事电磁场数值计算方法研究． 第5期 匡磊等：三维随机粗糙面与目标复合电磁散射的聊D方法 用Monte c耐。方法产生相关长度分别为f。和z，的 铂uss谱随机粗糙面，已知在二维粗糙面问题中，粗糙面 有限长度￡与相关长度f要满足”21 Z， (1) 工≥15 在三维问题中，Gauss谱密度函数及相关函数分别为 旷(疋，巧)=皇若争exp[一(足：圮+x；E)，4]，(2) c(，，y)：ex“一菩一事】，一 (3) -l ●f一 其中丘=2Ⅱ*，L，墨=2“y，L，分别为*和y方向上的空 间频率，^为均方根(nm)表面高度，设^=0．1m，f，=f， =z=o．2 m，用文[15]方法，随机实现20次粗糙面，平均 图l周期性延拓的三维粗糙面 后的随机粗糙面的础高度的相对误差和相关性的均方 3D Fig．1 r0。gh鲫血ce。iⅡ1p