THE method of moment (MOM) is widely used to extract. Fast Green Function Evaluation for Method of Moment. arxiv: v1 [cs.

Transcription

1 1 Fas Green Funcon Evaluaon for Mehod of Momen Shunchuan Yang, Member, IEEE, Dongln Su, Member, IEEE arxv: v1 [cs.ce] 14 Jan 2019 Absrac In hs leer, an approach o accelerae he marx fllng n mehod of momen (MOM) s proposed. Based on he fac ha he Green funcon s dependen on he Eucldean dsance beween he source and he observaon pons, we consruced an effcen adapve one-dmensonal nerpolaon approach o fas calculae he Green funcon and Exp funcon values. In he proposed mehod, several adapve nerpolaon ables are consruced based on he maxmum and mnmum dsance beween any wo negraon pons wh local refnemen near zero funcon values o oban relavely unform error rae n he overall compuaonal doman. An effcen approach o oban he samplng pons used n he nerpolaon phase s carefully desgned. Then, any funcon values could be effcenly calculaed hrough a lnear nerpolaon mehod for Exp and a Lagrange polynomal nerpolaon mehod for he Green funcon. In addon, he error bound of he proposed mehod s rgorously nvesgaed. The proposed mehod could be que easly negraed no he avalable MOM codes for dfferen negraon equaon (IE) formulaons, lke he elecrcal feld negral equaon (EFIE), magnec feld negral equaon (MFIE), he combned feld negral equaon (CFIE), he Poggo- Mller-Chan-Harrngon-Wu-Tsa(PMCHWT), wh few effors. Comprehensve numercal expermens valdae s accuracy and effcency hrough several IE formulaons. Resuls show ha over 20% effcency mprovemen could be acheved whou sacrfcng he accuracy. Index Terms Green funcon, nerpolaon, local refnemen, mehod of momen I. INTRODUCTION THE mehod of momen (MOM) s wdely used o exrac compose parameers of large-scale negraed crcus (ICs) [1], solve elecrcally large, mulscale scaerng problems [2] due o s unknowns resdng on he nerface of dfferen homogenous meda. Therefore, he coun of unknowns s much smaller han ha of oher dfferenal equaon based mehods, lke he fne-dfference me-doman (FDTD) mehod [3], he fne elemen mehod (FEM) [4], whch requre volumerc dscrezaon. However, a dense and full mpedance marx n MOM s consruced. Invesgaons found ha marx fllng n MOM s que me-consumng and s one of he bolenecks o prohb MOM from solvng praccal problems, whch s drecly relaed o calculae he Green funcon, Exp( jkr)/r, Manuscrp receved xxx; revsed xxx. Ths work was suppored n par by he Naonal Naural Scence Foundaon of Chna hrough Gran , Gran , he 13h Fve- Year Equpmen Pre-research Fund hrough Gran and Bejng Naural Scence Foundaon hrough Gran All auhors are wh he School of Elecronc and Informaon Engneerng, Behang Unversy, Bejng, 10083, Chna (e-mal: scyang@buaa.edu.cn, sdl@buaa.edu.cn). and Exp( jkr). When consrucng he mpedance marx, one could drecly evaluae hose analycal expressons. However, s repored ha s me-consumng o calculae Exp funcon, especally for complex numbers [5], whch s he exac scenaro n MOM for he me-harmonc elecromagnec smulaons. In hs paper, we show ha drec evaluaon of hose analycal formulaons s no opmal snce a large number of me-consumng and redundan evaluaon of Exp funcon may exs. We proposed an effcen adapve nerpolaon approach o fas evaluae Exp( jkr)/r and Exp( jkr). Through careful desgn of nerpolaon ables and an effcen approach o rereve he samplng pons, we could effcenly calculae Exp( jkr)/r and Exp( jkr) values wh conrollable accuracy and sgnfcan performance mprovemen. Ths paper s organzed as follows. In Secon II, some key observaons relaed o he negral denes used n MOM are presened and he proposed approach s llusraed n deals. In Secon III, s error bound s rgorously analyzed. In Secon IV, s accuracy and effcency are comprehensvely nvesgaed hrough numerous expermens. A las, we draw some conclusons n Secon V. II. METHOD A. Observaons of Inegral Idenes Used n MOM To effcenly consruc he mpedance marx n MOM, several ypes of negrand are nvolved. In hs paper, sngular cancellaon echnque along wh angular ransformaon [6] s used o demosrae he proposed mehod. However, could be easly exended o oher negraon echnques whou any dffcules. The followng four negrals are requred o be evaluaed [6], I a = S R dr = 1 3 h ( k =1 x h 2 e jkd) dx (1) x2 I b = ρ e jkr S R dr = 1 3 û dx (2) k =1 x (1jkR) I c = dr S R 3 ( ) h e jkr = F h 2 x2 R e jkd dx (3) d =1 x

2 I d = = F ρ e jkr (1jkR) dr S R 3 =1 x û dx (4) R where k s he wavenumber, d s he normal dsance from he evaluaon pon o he source rangle plane, h s he dsance from he projecon pon of he evaluaon pon on he source plane o he negraon lne, x s he angular coordnae a he negral lne. To effcenly evaluae (1), we obvously need o fas calculae hree exponenal funcons Value Real/Imag par of Exp(-jkr) Real Imag Value Real/Imag par of Exp(-jkr)/r Fg. 1. Real and magnary pars of Exp( jkr) and Exp( jkr)/r. Real Imag 2, e jkr R, e jkd. (5) d I s nrnsc o evaluae (1) hrough (5) snce hey are analycal expressons and we could oban he funcon values easly. However, as saed n [5], s que me-consumng for exponenal funcon evaluaon snce numerous numercal operaons are requred o oban convergence resuls. Unforunaely, even worse reducan calculaons of hose funcons may exs when we fll he mpedance marx n MOM. To beer undersand hs problem, we could roughly esmae he coun of funcon calls o consruc he mpedance marx as follows. To oban he quanave analyss, we make he followng assumpons. The srucures are frs dscrezed no N planar rangles. n-pon numercal negraon a each edge of he source rangle o evaluae he nner surface negraon s used and hree lnes overall for each rangle are requred. m-pon numercal negraon on average for he ouer surface negraon s consdered. Therefore, he oal coun of exponenal funcon calls approxmaely equal o 3mnN 2 o generae he mpedance marx. I s easy o see ha for a large number N he coun s exremely large because depends on he square of N. When acceleraon echnques, lke MLFMP [2], pfft [1] and AIM [7], are used, he oal coun of funcon calls s approxmaely reduced o 3mnN C, where C s he averaged coun of rangles n he near feld regon and we only consder he near-feld neracon calculaon. In hose scenaros, n whch he elecrcally large or mulscale srucures are ofen nvolved, N self s an exremely large number. Therefore, we sll make a large number of funcon calls o evaluae (5) whch s nevably me-consumng. Furher carefully nvesgaons could fnd ha reducan funcon calls possbly exs n he drec evaluaon mehod snce (5) are funcons of dsance beween he feld pon or s projec pon o he source pon. In he followng subsecon, an adapve nerpolaon mehod along wh an effcen approach o rereve samplng pons s proposed o fas calculae (5). B. Effcen Evaluaons of Exp Funcons To effcen evaluae values for he expoenal ype funcons wh complex numbers, n he proposed mehod we should frs deermne samplng nervals. As shown n Fg. 1 and, real and magnary pars of Exp( jkr) and Exp( jkr)/r show dfferen oscllang behavors, whch mples ha specal aenon mus be pad o selecng approprae samplng pons. Exp( jkr) s a perodc funcon whou any sngulary near he orgn. Therefore, we could use a relavely unform samplng nerval and a lnear nerpolaon mehod o oban accurae nerpolaed values. For Exp( jkr)/r, we sar o sample a he mnmum dsance R mn near he orgn snce here s a sngulary pon where shows large funcon varaon. Large nerpolaon error could occur f he lnear nerpolaon mehod s sll used. In he proposed mehod, we selec he Langrage polynomal mehod o nerpolae Exp( jkr)/r. We frs unformly sample he funcon usng he nervals defned as follows = λ εr µ r λ εr µ r. (6) where λ 0, ǫ r, µ r denoe he wavelengh n free space, relave permvy and relave permeably of resdng homogenous medum, respecvely. As saed n (6), o oban accepable accuracy of funcon values, we use samplng pons per wavelengh due o he error possbly accumulang n he negraon phase when we evaluae double surface negrals o consruc he mpedance marx. Snce we only requre o consruc one-dmensonal nerpolaon, he compuaonal resource consumpons n erms of memory usage and CPU me could be gnored compared wh hose n he smulaon even f we use such a dense samplng rae. However, varous numercal expermens reveal ha que large relave nerpolaon error wll occur when he desred funcon values are small near he zero funcon values. We proposed he followng sraegy o keep he error unformly dsrbued n he whole compuaonal doman. Local refnemen samplng s used near he zero funcon values. When samplng locaon s near zero funcon values, a local refnemen o use smaller samplng nerval s used. Zero funcon values s forced o be sampled. If we mss zero funcon values n, we drecly sample zero funcon values when consrucng he nerpolaon ables. One mus pay specal aenon o avodng ha wo samplng pons are oo close, whch could keep he hash array nroduced n he nex subsecon wh reasonable dmenson. Therefore, he proposed mehod could effcenly ulze he memory and consruc an opmal nerpolaon ables. I s a banlance beween he accuracy and compuaonal cos. Wh he defnon of he samplng pons, several ables are effcenly consruced hrough (5) whou any redundan

3 3 calculaons. Therefore, we could sgnfcanly mprove he effcency of marx fllng n MOM. Varous nerpolaon mehods, lke he lnear nerpolaon mehod and he Lagrange nerpolaon mehod [9], are avalable o effcenly nerpolae he connuous funcon values for Exp( jkr) and Exp( jkr)/r. They are wdely used n he mullayer Green funcon evaluaons [7][10] snce s spaal analycal expresson s unavalable. We proposed o use he lnear nerpolaon mehod o nerpolae Exp( jkr) and he Lagrange polynomal mehod for Exp( jkr)/r. Therefore, he nerpolaed funcon values could be expressed as n f (r) = f (r j )P n,j, (7) where r s he desred pon locaon, r j s he samplng locaon, n s he number of samplng pons n he local nerpolaon, P n,j s he lnear nerpolaon funcon or he Lagrange nerpolaon polynomal of degree n. C. Effcen Searchng Elemens n he Samplng Table The dmenson of samplng ables s usually que large for elecrcally large or mulscale srucures. If we navely search elemens one-by-one n he samplng ables requred by he nerpolaon mehods, he performance wll be severely deeroraed and even worse han ha of orgnal analycal mehod. We proposed he followng effcen approach o mgae he problem. The hash-able mehod s desgned o fas oban he desred elemens n a large conaner [11]. We use he hash mehod o map he dsance r o he ndex of samplng pon n he samplng ables. The followng many-o-one hash funcon s proposed [ ] R Rmn (8) R mn where R s he locaon of desred funcon values, R mn s he mnmum samplng locaon, R mn s he mnmum samplng nerval, [ ] denoes he neger par of he correspondng floa number. I could be found ha once he desred r s obaned we could fas locae he samplng pon n he ables hrough (8). On he mplemenaon aspec, one could smply consruc a one dmensonal array wh he dmenson of [R max R mn ]/( R mn ), where R max s he maxmum samplng locaon. We frs fll he array wh he values of samplng locaons usng (8) and hen oher posons are recursvely flled by s prevous elemens sarng a he lower bound of he hash array. number could be easly derved accordng o he Mean Value heorem smlar n [9] as follows. max f n1 (ξ r ) max f (r) f n(r) max ω r [a,b] n1(r), r [a,b] r [a,b] (n1)! (9) where n ω n1 (r) = (r r j ), (10) f(r) s he complex funcon value, f n (r) s he nerpolaed funcon value n [a,b], f n1 s he n-h dervave of f. For he unform samplng pons, denoed as r = a (b a), = 0,1,...,n, (11) n we could arrve a he analycal error bound wh some mahemacal manpulaons based on (9) as follows ( ) n1 max f n1 (ξ r ) b a max f (r) f r [a,b] n(r). r [a,b] n 4(n1) (12) I s obvous ha he error decreases as we ncrease he nerpolaon order or use smaller nerpolaon nervals. Ths s he reason why we locally refne he nervals near zero funcon values o oban unform relave error bound n he whole smulaon doman. Therefore, o balance memory cos and accuracy, we use he unform samplng far away from zero funcon values and non-unform samplng near zero funcon values n he proposed mehod. B. Error Esmaes for Lagrange Inerpolaon I s easy o derve he error bound for he Lagrange nerpolaon hrough he Generalzed Rolle s heorem [9] as follows max f (r) f n(r) r [a,b] n f n1 (ξ r ) r r j. (13) (n1)! The general remarks for (13) are he same as ha of (12). The Lagrange nerpolaon shows beer performance compared wh he lnear nerpolaon for hghly oscllang kernel. However, he Lagrange mehod s more me-consumng compared wh he lnear nerpolaon. As saed before, we proposed he lnear nerpolaon for Exp( jkr) and he Lagrange nerpolaon mehod for Exp( jkr)/r o balance he accuracy and effcency. III. ERROR ESTIMATES In hs sesson, we wll analyze he error bound for he proposed mehod wh he lnear and Lagrange polynomal nerpolaon mehod for Exp( jkr) and Exp( jkr)/r. A. Error Esmaes for Lnear Inerpolaon Snce Exp( jkr) and Exp( jkr)/r are nfne smooh, he n-h dervave of hese wo funcons are always avalable. The error bound for he polynomal nerpolaons for complex IV. NUMERICAL EXPERIMENTS AND DISCUSSION A. Accuracy of he proposed mehod We frs valdae he accuracy of he proposed mehod wh an nerpolaon problem. The confguraons are se up as, R max = 1 m, correspondng o one wavelengh a 300 MHz, R mn = m, he samplng nerval far away from zero funcon value s, local refnemen samplng nerval s /2, and local refnemen range s 4. We nerpolaed 10 4 funcon values unformly dsrbued n [ ]. The 2-pon lnear

4 4 Relave Error Relave Error Max. error whou local refnemen /r /r Real par of Exp(-jkr) Real par of Exp(-jkr) Max. error whou local refnemen (c) /r Real par of Exp(-jkr) Relave Error Relave Error Max. error wh local refnemen Max. error wh local refnemen /r Real par of Exp(-jkr) Fg. 2. and denoe he maxmum relave error of real and magnary par of Exp( jkr) and Exp( jkr)/r obaned from lnear nerpolaon mehod, (c) and (d) from Lagrange polynomal mehod. mehod and Lagrange polynomnal mehod of order 3 are used n our smulaons. The maxmum relave error s defned as ( ) f (r) f (r) E max = max. (14) f (r) I s easy o see ha n Fg. 2 and, (c) and (d), he proposed local refnemen mehod could approxmaely mprove he accuracy up o wo order for boh he real and magnary pars of Exp( jkr) and Exp( jkr)/r. Compared Fg. 2 and wh (c) and (d), he lnear mehod shows smlar performance n erms of he maxmum relave error for Exp( jkr) as ha of Lagrange polynomal mehod. Snce he lnear mehod has hgher effcency han ha of Lagrange, we use he lnear mehod o nerpolae Exp( jkr) o acheve beer effcency n he proposed mehod. As shown n Fg. 2, he Lagrange polynomal mehod shows sgnfcan beer performance han he lnear mehod for Exp( jkr)/r. Therefore, we use he Lagrange mehod for he Green funcon. B. Performance of he Proposed Mehod Incorporaed n MOM based on Varous Formulaons To comprehensvely es he performance of he proposed mehod ncorporaed n MOM, hree scenaros, a sphere wh he rad of 0.5 wavelengh s consdered for varous negral formulaons ncludng he EFIE, he MFIE, he CFIE and he PMCHWT. For he EFIE, MFIE, CFIE, he sphere s perfecly elecrcal conducor (PEC) and homogenous delecrc srucure wh ǫ r = 4 for he PMCHWT. A x-polarzed plane wave ncdens along z axs. All numercal expermens are carred ou on a worksaon wh an nel GHz CPU and 32 G memory and all n-house codes are runnng wh a hread whou explorng any he parallel compuaon o make far comparsons. (d) (c) Fg. 3. The relave error of real and magnary par of elemens n mpedance marx hrough EFIE and MFIE obaned from he proposed mehod. TABLE I COMPARISON OF TIME COST OF DIFFERENT FORMULATIONS BETWEEN THE ANALYTICAL METHOD AND THE PROPOSED METHOD Merc EFIE (d) MFIE Analycal Proposed Analycal Proposed Couns of unknowns Marx fllng me [s] Overall me [s] In Fg. 3, we plo he relave error of all elemens n he mpedance marx obaned from he EFIE nvolvedl operaor and he MFIE nvolved K operaor [12]. Oher formulaons n essence are varous combnaon of he EFIE and MFIE, herefore, he performance of he proposed mehod n oher formulaons are smlar. As shown n Fg. 3 and, (c) and (d), he relave error of boh he real and magnary pars for he EFIE reaches 10 5 and For he MFIE, he relave error s furher decreased o 10 6 and Tha mples he proposed mehod could provde an excellen level of elemenwse accuracy. The L 2 error of RCS for boh EFIE and MFIE could reach Those resuls show excellen performance n erms of accuracy for he proposed mehod. Table I shows he me cos for he EFIE and he MFIE ncorporaed wh he analycal mehod and he proposed mehod. For he EFIE, marx fllng n he analycal mehod akes 85s. However, when ncorporaed wh he proposed mehod, marx fllng only akes 61s wh L 2 error of RCS I shows 27% and 24 % effcency mprovemen for he EFIE and he MFIE formulaons, respecvely, whou sacrfcng accuracy. We perform a number of numercal expermens wh he EFIE, MFIE, CFIE and PMCHWT formulaons for dfferen

5 5 srucures. Due o lmed page space, we only showed par of resuls. However, he proposed mehod shows que good performance of boh he accuracy and effcency, whch mply ha he proposed mehod s general, robus and effcen o accelerae marx fllng n MOM. V. CONCLUSIONS In hs paper, we proposed an effcen nerpolaon approach o accelerae marx fllng n MOM. In he proposed mehod, several nerpolaon ables are consruced wh carefully consderaon of local refnemen near he zero funcon values. The lnear mehod, Lagrange polynomal are used o nerpolae Exp( jkr) and Exp( jkr)/r, respecvely. An effcen approach o oban samplng pons are also presened. Compared wh he analycal mehod, shows over 20% effcency mprovemen whou sacrfcng accuracy. In all, he proposed mehod s smple, robus and effcen o accelerae marx fllng for MOM and applcable for almos all SIE formulaons and s que easy o negrae he avalable MOM codes wh only a few modfcaons. REFERENCES [1] Z. Zhu, B. Song, and J. K. Whe, Algorhms n FasImp: a fas and wde-band mpedance exracon program for complcaed 3-D geomeres, IEEE Trans. Compu.-Aded Desgn Inegr. Crcus Sys., vol. 24, no. 7, pp , Jul., [2] J. Song, C. Lu, W. Chew, Mullevel fas mulpole algorhm for elecromagnec scaerng by large complex objecs, IEEE Trans. Anennas Propag., vol. 45, no. 10, pp , Oc., [3] A. Taflove and S. Hagness, Compuaonal Elecrodynamcs: he Fne- Dfference Tme-Doman Mehod. Arech house, [4] J. Jn, The fne elemen mehod n elecromagnecs. New York, NY, USA:Wley, [5] N. Schraudolph, A fas, compac approxmaon of he exponenal funcon, Neural Compuaon, vol. 11, no. 4, pp , May, [6] Z. G. Qan, W. C. Chew, and R. Suaya, Generalzed Impedance Boundary Condon for Conducor Modelng n Surface Inegral Equaon, IEEE Trans. Mcrow. Theory Techn., vol. 55, no. 11, pp , Nov [7] E. Bleszynsk, M. Bleszynskand, and T. Jaroszewcz, AIM: Adapve negral mehod for solvng large-scale elecromagnec scaerng and radaon problems, Rado Sc., vol. 31, no. 5, pp , Sep, [8] K. Yang and A. Ylmaz, A hree-dmensonal adapve negral mehod for scaerng from srucures embedded n layered meda, IEEE Trans. Geosc. Remoe Sens., vol. 50, no. 4, pp , Apr [9] R. Kress, Numercal analyss, Sprnger, [10] Y. Lu, L. L, T. Yeo, M. Leong, Applcaon of DCIM o MPIE-MoM analyss of 3D PEC objecs n mullayered meda, IEEE Trans. Anennas Propag., vol. 50, no. 2, pp , [11] A. Aho, D. Jeffrey, Daa srucures and algorhms. Pearson, [12] W. Gbson, The mehod of momens n elecromagnecs. CRC Press, 2007.