随机粗糙面散射的微扰法及有限元法研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
本论文采用微扰法(small perturbation method,SPM)和有限元法(finite elementmethod,FEM)深入系统地研究了随机粗糙面的散射问题以及相关问题。首先总结归纳了现有的应用最广泛的三种锥形波,并讨论了这三种锥形波的锥型因子g的合理取值范围;随后将Thorsos锥形波引入SPM中,并结合Monte Carlo方法生成了混合方法—混合SPM,并利用混合SPM模拟了一维导体粗糙面的近场散射特性;研究了两种经典的用于减小刚度矩阵带宽和外形的结点重编码算法,并基于经典的结点重编码算法—GPS(Gibbs-Poole-Stockmeyer)算法提出了一种新的结点重编码算法—GGPS(Generalized GPS)算法;基于GPS算法中的伪直径端点寻找算法提出了一种新的伪直径端点寻找算法—WDR(width-depth-ratio)算法,用WDR算法替换GPS算法中的伪直径端点寻找算法得到改进版的GPS算法—GPS-WDR算法;将GPS-WDR算法引入FEM/ABC(absorbing boundary condition,ABC)模拟了一维导体高斯随机粗糙面的散射,并讨论了其远场特性。论文主要工作如下:
     1.介绍了随机粗糙面建模涉及的理论和描述粗糙表面的统计参量—均方根高度h和相关长度l,介绍了随机粗糙面建模的的蒙特卡罗方法,生成了一维和二维高斯随机粗糙面,并讨论了粗糙面表面的高度起伏与均方根高度h和相关长度l之间的关系;归纳总结了现有的应用最广泛的三种锥形波—Thorsos锥形波(空间域锥形波)、谱域锥形波和角域锥形波,讨论了三种锥形波的锥型因子g的合理取值范围,给出了Thorsos锥形波和谱域锥形波入射下双站散射系数的推导过程。
     2.介绍了积分方程法、SPM和Monte Carlo法,提出一种用来计算Thorsos锥形波入射下一维导体高斯随机粗糙面的散射场的新的快速模拟算法—混合SPM,推导出了TE锥形波入射下的电场散射波公式。采用高斯随机粗糙面模拟微粗糙面,讨论了方法中涉及到的近场表达式的积分方程中的谱域变量k x的有效范围。通过数值计算得到了TE极化散射电场的表面电场结果。给出了垂直入射下混合SPM的近似公式,具有更简单的形式,并且数值精确度与矩量法(Method of Moments,MOM)比较后得到了验证。相干波和非相干波随着均方根高度h的变化而变化的规律以及与此相关的物理解释也得到了讨论。
     3.研究和讨论了电磁场有限元结点编码优化的两种典型算法:改进的AD(Akhras和Dhatt)算法和GPS算法,给出了实施两种算法的具体措施。以球冠型共形天线模型的网格结点数据作为实际算例,比较了两种算法对该类问题优化效果的优劣。在GPS方法的基础上,提出了两种新的结点重编码方法—GGPS法和GPS-WDR法,给出了算法的理论依据和详细过程,并通过有限元实例验证了GGPS方法和GPS-WDR方法的有效性,提高了FEM刚度矩阵方程的直接法求解速度。
     4.将GPS-WDR引入FEM/ABC并应用于模拟一维导体高斯随机粗糙面的散射,有限元的网格模型通过商业软件ANSYS的脚本语言APDL实现,所得到的双站散射系数的数值结果和MOM得到的结果高度吻合,证明了FEM/ABC的准确性,比较了单纯用GPS算法优化刚度矩阵带宽和外形后的矩阵计算时间和利用GPS-WDR优化刚度矩阵后的矩阵计算时间,时间结果显示了后者的时间优越性,讨论了均方根高度h和相关长度l的变化对双站散射系数的影响,并且给出了结果的合理的物理解释。
This dissertataion is mainly about the use of small perturbation method (SPM) andFEM in the electromagnetic scattering from random rough surface and issues related tothis topic. Three mostly used forms of tapered waves are discussed. The range oftapered parameter g is discussed, and proper ranges for all the three tapered waves aregiven. The tapered wave is introduced into the classical SPM to simulate the scatteringfrom1-D random rough surface. Two classical nodes renumbering algorithms for thereduction of matrix bandwidth and profile are studied and discussed. A novel nodesrenumbering algorithm-GPS (Gibbs-Poole-Stockmeyer) algorithm is proposed basedon GPS algorithm. A novel pseudo-peripheral nodes finding algorithm-WDR(width-depth-ratio) algorithm is proposed based on the pseudo-peripheral nodes findingalgorithm in the GPS algorithm. The WDR algorithm replaces the pseudo-peripheralnodes finding algorithm in GPS algorithm and an improved GPSalgorithm-GPS-WDR algorithm is developed. The GPS-WDR algorithm is introducedto FEM/ABC to simulate the scattering from1-D conducting Gaussian random roughsurface, and the far field results are discussed. The main work of the dissertation is:
     1. The basic theories and the parameters to describe the random rough surface isintroduced-root mean square (rms) height h and correlation length l aremainly introduced. The random rough surface generating method-MonteCarlo method is introduced, and the1-D and2-D random rough surface aregenerated. The relationship between the profile of the rough surface and rmsheight h are discussed, as well as the relationship between profile andcorrelation length l. Three mainly tapered waves–Thorsos tapered wave,spectral domain tapered wave and the angular domain tapered wave aresummarized and introduced. The proper ranges of tapered parameter g for allthe three tapered waves are discussed and given. The derivations of bistaticscattering coefficient under Thorsos tapered wave incident are given.
     2. The integral equation method, SPM and Monte Carlo method are introduced.A novel analytical method-hybrid SPM is proposed to simulate the scatteringfrom1-D conducting Gaussian random rough surface. The equations for thescattering electric field are derived. The range for the spectral domain variablek xis discussed. The surface field for the scattering field under tapered wave incident is given by hybrid SPM. The approximation equations for hybridSPM with normal incidence are given. The approximation equations havesimpler forms and the accuracy is testified. The varying of the coherent andincoherent waves with the changes of rms height h has been discussed, andthe physical explanation has been given, too.
     3. The two classical nodes renumbering algorithms-the improved AD algorithmand the GPS algorithm are studied and discussed. The steps of the improvedAD algorithm and GPS algorithm are given. Taking the spherical crownconformal antenna as the model to testify both the algorithms, the results showtheir advantages and disadvantages in the application as the nodesrenumbering algorithm. Two novel nodes renumbering algorithms areproposed-GGPS algorithm and GGPS-WDR algorithm, and the use of thesetwo algorithms can accelerate the speed of solving the stiffness equations inFEM.
     4. The GPS-WDR is introduced to FEM/ABC, and the combined method is usedto simulate the scattering from a1-D conducting Gaussian random roughsurface. The bistatic scattering coefficient is obtained, and the results arecompared with the bistatic scattering coefficient obtained by MOM to testifyits accuracy. The effect of the varying rms height and the correlation length tobistatic scattering coefficient is discussed. The proper physical interpretation isgiven.
引文
[1] Yang, T. Numerical studies of rough surface scattering models. Dissertation.Washington State University.1993.
    [2]陈颖.小粗糙度表面上方物体的散射. Dissertation.复旦大学.2008.
    [3]李娟.粗糙面及其与目标复合电磁散射的FDTD方法研究. Dissertation.西安电子科技大学.2010.
    [4]柴草.矩量法及其加速算法研究一维粗糙面与目标复合电磁散射. Dissertation.西安电子科技大学.
    [5] Kuo, C. Full wave two-dimensional modeling of scattering and inverse scatteringfor layered rough surfaces with buried objects. ProQuest2008.
    [6] Johnson, J.T.andBurkholder, R.J. A study of scattering from an object below arough surface. Geoscience and Remote Sensing, IEEE Transactions on.2004,42(1).59.
    [7] Johnson, J. A numerical study of scattering from an object above a rough surface.Antennas and Propagation, IEEE Transactions on.2002,50(10).1361.
    [8] Tsang, L., Kong, J., Ding, K.andAo, C. Scattering of Electromagnetic Waves,Numerical Simulations.2001,
    [9] Yaqiu, J., Peng, L.andHongxia, Y. Theory and Method of Numerical Simulation ofComposite Scattering from the Objects and Randomly Rough Surface. Beijing:Science Press2008.
    [10]Rice, S.O. Reflection of electromagnetic waves from slightly rough surfaces.Communications on pure and applied mathematics.1951,4(2‐3).351.
    [11]Beckmann, P.andSpizzichino, A. The scattering ofelectromagnetic wavesfromrough surfaces. New York: Macmillin.1963,218.
    [12]陶爱华.粗糙面散射及地球极区中间层顶PMSEs现象研究. Dissertation.中国科学院研究生院(电子学研究所).2007.
    [13]叶红霞.随机粗糙面与目标复合电磁散射的数值计算方法. Dissertation.复旦大学.2007.
    [14]张景剑.海面上方目标复合电磁散射的实验测量与理论研究. Dissertation.西安电子科技大学.2010.
    [15]Thorsos, E. The validity of the Kirchhoff approximation for rough surfacescattering using a Gaussian roughness spectrum. J. Acoust. Soc. Am.1988,83(1).78.
    [16]任新成.粗糙面电磁散射及其与目标的复合散射研究. Dissertation.西安电子科技大学.2008.
    [17]王蕊.粗糙面及其与目标复合电磁散射中的相关问题研究. Dissertation.西安电子科技大学.2009.
    [18]刘长学.超大规模稀疏矩阵计算方法.上海科学技术出版社1991.
    [19]Cuthill, E.andMcKee, J., Reducing the bandwidth of sparse symmetric matrices,Proceedings of the196924th national conference,1969, pp.157.
    [20]George, A.andLiu, J.W.H. An Implementation of a Pseudoperipheral Node Finder.ACM Trans. Math. Softw.1979,5(3).284.
    [21]Gibbs, N.E., Poole Jr, W.G.andStockmeyer, P.K. A comparison of several bandwidthand profile reduction algorithms. ACM Transactions on Mathematical Software.1976,2(4).322.
    [22]Gibbs, N.E., Poole Jr, W.G.andStockmeyer, P.K. An Algorithm for Reducing theBandwidth and Profile of a Sparse Matrix. SIAM Journal on Numerical Analysis.1976,13(2).236.
    [23]Pinana, E., Plana, I., Campos, V.andMarti, R. GRASP and path relinking for thematrix bandwidth minimization. European Journal of Operational Research.2004,153(1).200.
    [24]Souza, L.T.andMurray, D.W. An alternative pseudoperipheral node finder forresequencing schemes. International Journal for Numerical Methods in Engineering.1993,36(19).3351.
    [25]Lewis, J.G. Implementation of the Gibbs-Poole-Stockmeyer and Gibbs-Kingalgorithms. ACM Transactions on Mathematical Software.1982,8(2).180.
    [26]王新荣,梁迎春and王细洋.有限元网格结点优序标号的WGPS法.哈尔滨工业大学学报.1992,24(4).99.
    [27]Beckman, P.andSpizzichino, A. The Scattering of Electromagnetic Waves FromRough Surfaces (Artech House, Norwood, MA,1987). Author Affiliations MichaelGreenstein The author is with Hewlett Packard Laboratories.3500(
    [28]Ogilvy, J.A.andMerklinger, H.M. Theory of wave scattering from random roughsurfaces. The Journal of the Acoustical Society of America.1991,90(3382.
    [29]Fung, A.K.andChen, K.S. Microwave scattering and emission models and theirapplications.1994,
    [30]Leung Tsang, Jin Au KongandKung-Hau. Scattering of Electromagnetic Waves:Theories and Applications. John Wiley&Sons Inc.2000.
    [31]Wagner, R.L., Song, J.andChew, W.C. Monte Carlo simulation of electromagneticscattering from two-dimensional random rough surfaces. Antennas and Propagation,IEEE Transactions on.1997,45(2).235.
    [32]Kuang, L.andJin, Y.Q. Bistatic scattering from a three-dimensional object over arandomly rough surface using the FDTD algorithm. Antennas and Propagation,IEEE Transactions on.2007,55(8).2302.
    [33]Lou, S., Tsang, L., Chan, C.andIshimaru, A. Monte Carlo simulations of scatteringof waves by a random rough surface with the finite element method and the finitedifference method. Microwave and Optical Technology Letters.1990,3(5).150.
    [34]Chih-hao, K.andMoghaddam, M., Electromagnetic scattering from a buriedcylinder in layered media with rough interfaces, Antennas and Propagation SocietyInternational Symposium2006, IEEE,2006, pp.653.
    [35]Kuo, C.andMoghaddam, M. Scattering from multilayer rough surfaces based on theextended boundary condition method and truncated singular value decomposition.Ieee Transactions on Antennas and Propagation.2006,54(10).2917.
    [36]郭立新,王蕊and吴振森.随机粗糙面散射的基本理论与方法.北京:科学出版社2010.
    [37]Tsang, L., Kong, J., Ding, K.andAo, C. Scattering of Electromagnetic Waves,Numerical Simulations. John Wiley&Sons Inc.2001.
    [38]Tabatabaeenejad, S.A. Forward and inverse models of electromagnetic scatteringfrom layered media with rough interfaces. ProQuest2008.
    [39]Kim, H.andJohnson, J.T., Radar Image studies of scattering from random roughsurfaces, IEEE,2001, pp.351.
    [40]Kapp, D.A.andBrown, G.S. A new numerical method for rough-surface scatteringcalculations. Antennas and Propagation, IEEE Transactions on.1996,44(5).711.
    [41]Toporkov, J. Study of electromagnetic scattering from randomly rough ocean-likesurfaces using integral-equation-based numerical technique. Dissertation. VirginiaPolytechnic Institute and State University.1998.
    [42]薛谦忠.粗糙面波束散射研究及其应用. Dissertation.西安电子科技大学.2000.
    [43]Rui, W., Li-Xin, G.andJun, M. Electromagnetic scattering from two-layer roughinterfaces in the Kirchhoff approximation. Chinese Physics B.2009,18(3422.
    [44]Jun Ma, Lixin GuoandCheng, X. Unification of the Kirchhoff approximation andthe method of moment for optical wave scattering from the lossy dielectricGaussian random rough surface. Chinese Optics Letters.2009,7(3).
    [45]Chan, C.H., Tsang, L.andLi, Q. Monte Carlo Simulations of Large-ScaleOne-Dimensional Random Rough-Surface Scattering at Near-Grazing Incidence:Penetrable Case. IEEE Transactions on Antennas and Propagation.1998,46(1).142.
    [46]Johnson, J., Shin, R., Kong, J., Tsang, L.andPak, K. A numerical study of oceanpolarimetric thermal emission. Geoscience and Remote Sensing, IEEE Transactionson.1999,37(1).8.
    [47]Hastings, F.D., Schneider, J.B.andBroschat, S.L. A Monte-Carlo FDTD Techniquefor Rough Surface Scattering. IEEE Transactions on Antennas and Propagation.1995,43(11).1183.
    [48]El-Shenawee, M., Rappaport, C., Miller, E.andSilevitch, M. Three-dimensionalsubsurface analysis of electromagnetic scattering from penetrable/PEC objectsburied under rough surfaces: Use of the steepest descent fast multipole method.Geoscience and Remote Sensing, IEEE Transactions on.2001,39(6).1174.
    [49]叶红霞and金亚秋.随机粗糙面上介质目标差场散射的快速计算方法.中国科学: G.2006,36(002).142.
    [50]Tsang, L., Kong, J., Ding, K.andAo, C. Scattering of electromagnetic waves,Theories and Application. John Wiley&Sons Inc.2000.
    [51]Fuks, I.M. Wave diffraction by a rough boundary of an arbitrary plane-layeredmedium. Antennas and Propagation, IEEE Transactions on.2001,49(4).630.
    [52]Tabatabaeenejad, A.andMoghaddam, M. Bistatic scattering from three-dimensionallayered rough surfaces. Geoscience and Remote Sensing, IEEE Transactions on.2006,44(8).2102.
    [53]Ye, H.andJin, Y. Parameterization of the tapered incident wave for numericalsimulation of electromagnetic scattering from rough surface. Antennas andPropagation, IEEE Transactions on.2005,53(3).1234.
    [54]Braunisch, H., Zhang, Y., Ao, C., Shih, S., Yang, Y., Ding, K., Kong, J.andTsang, L.Tapered wave with dominant polarization state for all angles of incidence. Antennasand Propagation, IEEE Transactions on.2000,48(7).1086.
    [55]Kuo, C.H.andMoghaddam, M. Electromagnetic scattering from a buried cylinder inlayered media with rough interfaces. Antennas and Propagation, IEEE Transactionson.2006,54(8).2392.
    [56]Kuo, C.H.andMoghaddam, M. Scattering from multilayer rough surfaces based onthe extended boundary condition method and truncated singular valuedecomposition. Antennas and Propagation, IEEE Transactions on.2006,54(10).2917.
    [57]魏宗舒.概率论和数理统计教程.高等教育出版社2008.
    [58]Papoulis, A. Probability&statistics. Prentice-Hall1990.
    [59]Ishimaru, A. Wave propagation and scattering in random media. Wiley-IEEE Press1999.
    [60]丁锐.随机介质中波散射与传播的随机泛函分析. Dissertation.复旦大学.2010.
    [61]Ishimaru, A. Theory and application of wave propagation and scattering in randommedia. Proceedings of the IEEE.1977,65(7).1030.
    [62]Millar, R. The Rayleigh hypothesis and a related least-squares solution to scatteringproblems for periodic surfaces and other scatterers. Radio Science.1973,8(8,9).785.
    [63]Van Den Berg, P.andFokkema, J. The Rayleigh hypothesis in the theory ofdiffraction by a perturbation in a plane surface. Radio Science.1980,15(4).723.
    [64]Oppenheim, A.V., Schafer, R.W.andBuck, J.R. Discrete-time signal processing.Prentice hall Englewood Cliffs, NJ:1989.
    [65]金建铭,电磁场有限元分析,西安电子科技大学出版社,2001.
    [66]贾建军and彭颖红.三种基于图论的有限元结点编号优化算法.机械科学与技术(西安).1998,17(005).725.
    [67]刘长学.超大规模稀疏矩阵计算方法.上海:上海科学技术出版杜.1991,
    [68]Cuthill, E.andMcKee, J. Reducing the bandwidth of sparse symmetric matrices.Proceedings of the196924th national conference.1969,157.
    [69]George, J.A. Computer implementation of the finite element method. Dissertation.Dept. of Computer Science, Stanford University.1971.
    [70]Gibbs, N.E., Poole Jr, W.G.andStockmeyer, P.K. A Comparison of SeveralBandwidth and Profile Reduction Algorithms. ACM Transactions on MathematicalSoftware (TOMS).1976,2(4).322.
    [71]Luo, J.C. Algorithms for reducing the bandwidth and profile of a sparse matrix.Computers and Structures.1992,44(3).535.
    [72]Barnard, S.T., Pothen, A.andSimon, H. A spectral algorithm for envelope reductionof sparse matrices. Numerical Linear Algebra with Applications.1995,2(4).317.
    [73]黄志超,包忠诩and周天瑞.板料成形数值模拟中网格结点编号优化算法研究.金属成形工艺.2003,21(004).16.
    [74]George, A.andLiu, J.W.H. An Implementation of a Pseudoperipheral Node Finder.ACM Transactions on Mathematical Software (TOMS).1979,5(3).284.
    [75]欧阳兴,陈中奎and施法中.有限元网格结点编号.北京航空航天大学学报.2002,28(3).339.
    [76]Boutora, Y., Takorabet, N., Ibtiouen, R.andMezani, S. A New Method forMinimizing the Bandwidth and Profile of Square Matrices for Triangular FiniteElements Mesh. IEEE TRANSACTIONS ON MAGNETICS MAG.2007,43(4).1513.
    [77]姜涛,王安麟and朱灯林.有限元結点编号的综合带宽优化算法.机械设计.2005,22(11).3.
    [78]Crane Jr, H.L., Gibbs, N.E., Poole Jr, W.G.andStockmeyer, P.K. Algorithm508:matrix bandwidth and profile reduction [F1]. ACM Transactions on MathematicalSoftware.1976,2(4).375.
    [79]黄志超,包忠诩and周天瑞.有限元节点编号优化.南昌大学学报.2004,28(003).281.
    [80]任华.板料成形模拟中的曲面网格划分和节点编号优化的研究. Dissertation.南昌大学.2005.
    [81]贾建军and彭颖红.三种基于图论的有限元结点编号优化算法.机械科学与技术.1998,17(5).725.
    [82]王萍and吴海容.电磁场有限元节点编号优化的两种方法.电机与控制学报.1997,1(002).101.
    [83]Gibbs, N. Algorithm509: A Hybrid Profile Reduction Algorithm [F1]. ACMTransactions on Mathematical Software.1976,2(4).378.
    [84]Courant, R. Variational methods for the solution of problems of equilibrium andvibrations. LECTURE NOTES IN PURE AND APPLIED MATHEMATICS.1994,1.
    [85]Levy, S. Structural analysis and influence coefficients for delta wings. J. Aero. Sci.1953,20(7).449.
    [86]Silvester, P. Finite element solution of homogeneous waveguide problems. AltaFrequenza.1969,38(313.
    [87]Coccioli, R., Itoh, T., Pelosi, G.andSilvester, P. Finite-element methods inmicrowaves: A selected bibliography. Antennas and Propagation Magazine, IEEE.1996,38(6).34.
    [88]Jin, J.andJin, J.M. The finite element method in electromagnetics. Wiley New York1993.
    [89]Wilcox, C.H. An expansion theorem for electromagnetic fields. Communications onpure and applied mathematics.1956,9(2).115.
    [90]Mittra, R.andRamahi, O. Absorbing boundary conditions for the direct solution ofpartial differential equations arising in electromagnetic scattering problems.Progress In Electromagnetics Research.1990,2(133.
    [91]Engquist, B.andMajda, A. Absorbing boundary conditions for numerical simulationof waves. Proceedings of the National Academy of Sciences.1977,74(5).1765.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700