[版主推荐]從二面角RCS看FEKO MOM and PO 物理光學法計算整理

這個是之前FEKO的一個學習檔，在論壇上有整個 PDF (可以於資料庫中找到) Z}>+!Z  
eGI&4JgJ.  
關於這兩章圖片應該可以在更深入的討論，請各位大俠提供看法.. ::Pf\Lb>  
J*qo3aJjE  
9YwS"~Q =w  

图片:mom_po01.jpg

^>h2.AJ  
2VkA!o4nP  

图片:mom_po02.jpg

新的 RUNTIME )BLoj:gYn  
pGGx.&5#82  
Intel(R) Xeon(R) CPU           X5482  @ 3.20GHz; GenuineIntel family 6 model 7 stepping 6 brand id 0 hKW!kA=gZ  
2 physical CPUs with a total of 8 processors found (multi-core CPUs with max. 4 cores per physical CPU) ._z[T@!9  
Number of parallel processes:    8 on    2 CPUs dbLxm!;(  
4lfJc9J  
max. triangles:MAXNDR   =    25938 Q CB~x2C  
K1>(Fs$  
MOM         0.601hr(36.06 min) n;Etn!4M  
kXbdR  
MOM+ MLFMM  0.208hr(12.48 min) A?829<  
S=~8nr/V  
MOM+ PO_multiple_reflections_10  8.846hr(8 hr 50.76 min) LlG~aGhel  
v;@-bED(Qs  
T3 k#6N.  
XS/5y(W  
                    SUMMARY OF REQUIRED TIMES IN SECONDS >w=xGb7  
da&f0m U  
                                            CPU-time      runtime ' ! ls"qo  
Reading and constructing the geometry         0.390        0.390 %S<( z5  
Checking the geometry                         0.219        0.219 BV HO_  
Initialisation of the Greens function         0.000        0.000 ~cEr<mzR  
Calcul. of coupling for PO/Fock               0.000        0.000 F I\V6\B/  
Calcul. of matrix elements                  383.359      383.359 "A?_)=zZ  
Calcul. of right-hand side vector             1.032        1.028 o(@F37r{?  
Preconditioning system of linear eqns.       59.358       59.360 9#@CmiIhy  
Solution of the system of linear eqns.     1717.345     1717.347 3 e9fziQ~  
Determination of surface currents             0.000        0.000 D`^9 u K  
Calcul. of impedances/powers/losses           0.091        0.092 ba  
Calcul. of averaged SAR values                0.000        0.000 2=l!b/m  
Calcul. of power ideal receiving ant.         0.000        0.000 &Wup 7  
Calcul. of cable coupling                     0.000        0.000 L~FE;*>7  
Calcul. of electric near field                0.000        0.000 @*c) s_  
Calcul. of magnetic near field                0.000        0.000 XP-4=0zd  
Calcul. of far field                          0.375        0.375 +2ih!$T;7>  
                            other              0.237        0.236 Sm%MoFf  
                                         -----------  ----------- H;n(qBSB  
                         total times:       2162.406     2162.406 \;A\ vQ[  
               (total times in hours:          0.601        0.601) M^^u{);q  
~:srm#IX  
Specified CPU-times are referring to the master process only {XNu4d9w(  
Sum of the CPU-times of all processes:    17299.251 seconds (    4.805 hours) Kv:ih=?  
On average per process:                    2162.406 seconds (    0.601 hours) Ma3H
n  
+>5 "fs$Y  
Peak memory usage during the whole solution:    1.413 GByte jP";ll|c  
(refers to the master process only) gZs8BKO  
Sum of the peak memory of all processes:       11.236 GByte c!%:f^7g  
On average per process:                         1.405 GByte cNG6 A4  
Aq{m42EAj  
MOM+ MLFMM  0.208hr(12.48 min) 1xo<V5  
c};Qr@v
po  
N7xkkAS{  
                    SUMMARY OF REQUIRED TIMES IN SECONDS hD/bO  
%j=,c{`Q  
                                            CPU-time      runtime $V8B =k~  
Reading and constructing the geometry         0.454        0.453 xf@D<}~1  
Checking the geometry                         0.203        0.203 :8aIj_qds  
Initialisation of the Greens function         0.000        0.000 ?D6rFUs9;  
Calcul. of coupling for PO/Fock               0.000        0.000 ;G3{ e  
Calcul. of the FMM transfer function          0.500        0.500 |7I.DBjR;  
Fourier transform of FMM basis funct.         0.250        0.250 w`Dzk.2  
Calcul. of matrix elements                   12.360       12.359 u[oYVpe)IG  
Calcul. of right-hand side vector             1.145        1.137 #+eV5%Si  
Preconditioning system of linear eqns.        7.297        7.297 !2N#H~{  
Solution of the system of linear eqns.      720.079      720.087 -~h2^Oez  
Determination of surface currents             0.000        0.000 Q=9
S?p M  
Calcul. of impedances/powers/losses           0.077        0.079 ;
D[b25  
Calcul. of averaged SAR values                0.000        0.000 {hM*h(W~3  
Calcul. of power ideal receiving ant.         0.000        0.000 mYk5f_}  
Calcul. of cable coupling                     0.000        0.000 T`=N^Ca1!`  
Calcul. of electric near field                0.000        0.000 |C S[>0mV!  
Calcul. of magnetic near field                0.000        0.000 @]HV:7<q  
Calcul. of far field                          7.636        7.638 2'J.$ h3  
                            other              0.406        0.403 |
[TH ~o  
                                         -----------  ----------- $sO}l  
                         total times:        750.407      750.406 7j&l2Z  
               (total times in hours:          0.208        0.208) >}/"gx  
&w3LMOT  
Specified CPU-times are referring to the master process only YAZ=-@]`\  
Sum of the CPU-times of all processes:     6003.250 seconds (    1.668 hours) "FaG5X(  
On average per process:                     750.406 seconds (    0.208 hours) RS/%uxS?  
Nu{RF  
Peak memory usage during the whole solution:   62.039 MByte +Z[%+x92  
(refers to the master process only) ]IJ.}  
Sum of the peak memory of all processes:      481.675 MByte \o@b5z]e  
On average per process:                        60.209 MByte ?lU]J]  
s/Wg^(&M  
N_r*Ig  
DRVvW6s  
MOM+ PO_multiple_reflections_10  8.846hr(8 hr 50.76 min) ][b
z5aV  
$(}kau  
  Physical Optics approximation or the Fock-theory is used for i\rI j0+  
  triangles with the following labels: v.v3HB8p  
j&Ayk*  
        Label                 PO/Fock  NOSHADE  MAXREF ~eL7=G@{  
              Cuboid1.Face1     PO        0         10 ^LT9t2  
              Cuboid1.Face2     PO        0         10 M)EUR0>8  
a@,tf'Sr  
  Meaning of NOSHADE: sa+:c{  
   0 = Full ray tracing is performed, each side of the triangle Y@u{73H  
       can be illuminated. BoIe<{X(9  
   1 = No ray tracing, it is assumed that the triangle is directly K"b vUH  
       illuminated by the sources and by the whole MoM region. D
l/UZ@8pl  
       The side is determined automatically. ^m8\fCA*  
   2 = Ray tracing is performed, however, the triangle can only miHW1h[=  
       be illuminated on the side from which the normal vector N &p=4  
       is pointing away, not on the back. \{+7`4g  
xyz-T1ib  
>e-XZ2>Sj  
                    SUMMARY OF REQUIRED TIMES IN SECONDS L*h X_8J  
1xq1te)  
                                            CPU-time      runtime Zml9ndzT  
Reading and constructing the geometry         0.453        0.453 Q:VD2<2  
Checking the geometry                         0.219        0.219 <WBGPzVZE  
Initialisation of the Greens function         0.000        0.000 YQX>)'  
Calcul. of coupling for PO/Fock           31843.845    31843.846 B?Y%y@.  
Calcul. of matrix elements                    0.000        0.000 W.:kE|a.g  
Calcul. of right-hand side vector             0.000        0.000 hY'"^?OP  
Preconditioning system of linear eqns.        0.000        0.000 w6Tb<ja  
Solution of the system of linear eqns.        0.000        0.000 WM=kr$/3  
Determination of surface currents             0.000        0.000 Punbw\9!d,  
Calcul. of impedances/powers/losses           0.092        0.092 J(/ eR,ak  
Calcul. of averaged SAR values                0.000        0.000 '}4[m>/  
Calcul. of power ideal receiving ant.         0.000        0.000 ",8h>eEWK  
Calcul. of cable coupling                     0.000        0.000 6x/ X8zu  
Calcul. of electric near field                0.000        0.000 +S%@/q  
Calcul. of magnetic near field                0.000        0.000 p7H*Ff`  
Calcul. of far field                          0.466        0.468 5I(` s#O  
                            other              0.644        0.641 n7<<}wcV  
                                         -----------  ----------- bJB:]vs$  
                         total times:      31845.719    31845.719 s9`T%pg  
               (total times in hours:          8.846        8.846) s[s6E`Q  
;co{bk|rj  
Specified CPU-times are referring to the master process only YBS]JCO  
Sum of the CPU-times of all processes:   254765.749 seconds (   70.768 hours) YN>#zr+~  
On average per process:                   31845.719 seconds (    8.846 hours) ]P<&CEk  
o~>p=5t  
Peak memory usage during the whole solution:   12.974 MByte ( v ~/glf  
(refers to the master process only) [S</QS!  
Sum of the peak memory of all processes:      103.782 MByte 9CBB,  
On average per process:                        12.973 MByte

详细的报告。 3UN Jj&-`  
!&'xkw`  
1、MOM的时间主要用于构建矩阵和求解方程组。 >2g CM  
2、MLFMM的时间主要是求解方程组。同时发现： #sqD
Z]\B  
Calcul. of far field                          7.636        7.638 0U/:Tpyr  
计算远场比MOM慢。 3&?Tc|F+  
3、PO计算的时间全耗费在计算PO电流了，正因为多次反射导致了确定PO电流需要先考虑一次、再考虑二次激励来得到。 1;:2
=8  
[#R%jLEJ2  
结论：PO的多次反射有点鸡肋，不适合处理两个面都很大的情况。MLFMM的确挺强的。 |g'sRTKJ  
*10e)rzM  
nW drVT$  
楼主要是再做个多进程的并行对比，测试一下不同算法的加速比就更好了！

在略入射时 尤其高频会发生一些问题，比如遮蔽问题等等 5，85恰好是这个角度，所以才会那样

PO 2次 uit-Q5@~  
j#e.rNG  
total times in hours:          8.846hr kP)o=\|W{z  
,0Zn hS)kq  
  Physical Optics approximation or the Fock-theory is used for 'Y[A'.*}4  
  triangles with the following labels: B/=q_.1F>  
x~;EH6$5'/  
        Label                 PO/Fock  NOSHADE  MAXREF Q$)|/Y))  
              Cuboid1.Face1     PO        0         10 ,GX~s5S8  
              Cuboid1.Face2     PO        0         10 <8)s  
Fd[h9 G  
  Meaning of NOSHADE: r.W,-%=bL  
   0 = Full ray tracing is performed, each side of the triangle B~>cNj<  
       can be illuminated. R9l7CJM@  
   1 = No ray tracing, it is assumed that the triangle is directly B#aH
\$_U  
       illuminated by the sources and by the whole MoM region. EyPJvs  
       The side is determined automatically. |2@en=EYk  
   2 = Ray tracing is performed, however, the triangle can only 8A:^K:Q  
       be illuminated on the side from which the normal vector 9"aFS=><  
       is pointing away, not on the back. b%;59^4AjD  
bw<w u}ED  
                    SUMMARY OF REQUIRED TIMES IN SECONDS >C3NtGvy  
DA
nb.0  
                                            CPU-time      runtime DvX3/z#T  
Reading and constructing the geometry         0.453        0.453 )_v\{N  
Checking the geometry                         0.219        0.219 ^W"Q(sh  
Initialisation of the Greens function         0.000        0.000 fT3*>^Uv  
Calcul. of coupling for PO/Fock           31843.845    31843.846 2Nn1-wdhb  
Calcul. of matrix elements                    0.000        0.000
cH;TnuX  
Calcul. of right-hand side vector             0.000        0.000 W3/ 7BW`  
Preconditioning system of linear eqns.        0.000        0.000 z8[H:W#G  
Solution of the system of linear eqns.        0.000        0.000 ChTXvk
dH  
Determination of surface currents             0.000        0.000 tQ~<i %;  
Calcul. of impedances/powers/losses           0.092        0.092 ]&:b<]K3  
Calcul. of averaged SAR values                0.000        0.000 90T%T2K  
Calcul. of power ideal receiving ant.         0.000        0.000 _~[?>cF%  
Calcul. of cable coupling                     0.000        0.000 {KDgK  
Calcul. of electric near field                0.000        0.000 ^$IZLM?E~  
Calcul. of magnetic near field                0.000        0.000 Q)S>VDLA  
Calcul. of far field                          0.466        0.468 _E6}XNS  
                            other              0.644        0.641 V-_/(xt*  
                                         -----------  -----------
h4anr7g{  
                         total times:      31845.719    31845.719 8rwYNb.P  
               (total times in hours:          8.846        8.846) v'@b.R,  
Mjj}E >&  
Specified CPU-times are referring to the master process only kwHqv
O!G  
Sum of the CPU-times of all processes:   254765.749 seconds (   70.768 hours) x48'1&m  
On average per process:                   31845.719 seconds (    8.846 hours) gKN}Of@^1  
`4%;qLxngP  
Peak memory usage during the whole solution:   12.974 MByte 5_)@B]~nM  
(refers to the master process only) ISS\uj63M  
Sum of the peak memory of all processes:      103.782 MByte OB6J.dF[%  
On average per process:                        12.973 MByte

谁有二面角的算例，能否上传一下，我按材料上设置，每次算出来和测试值都不靠谱

是啊。传个模型上来阿！.cfx文件就行

算的模型对大呀，频率多高呀

引用第7楼流浪的老兵于2008-09-10 09:02发表的  :
s1N?/>lmB  
楼主的报告非常详细，学习了。
有个问题不明白， 请指点：
我对PO计算二次反射步骤的理解是： w'Tq3-%V  
....... -~{c u47_  

5z}w}zdg  
U(>4s]O6  
lSwcL  
認同您的說法..

引用第15楼信青蛙于2008-12-16 16:38发表的  : (Fq5IGs  
谁有二面角的算例，能否上传一下，我按材料上设置，每次算出来和测试值都不靠谱 K8n4oz#z  

t*z~5_/  
/kL$4CA  
3~,d+P  
http://bbs.rfeda.cn/read-htm-tid-14096.html