-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathexample1.log
87 lines (80 loc) · 4.84 KB
/
example1.log
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
./example1.out ex.dat
"ipw.txt" No.00 -- plane wave --
wave length of incident beam in vacuum : 1.06
refractive index of surrounding : 1
incident beam power per unit area : 0.001
x-component of polarization coefficient : 1+ 0i
y-component of polarization coefficient : 0+ 0i
x-component of translation vector : 0
y-component of translation vector : 0
z-component of translation vector : 0
rotation parameter theta [rad]: 0
rotation parameter phi [rad]: 0
---- sphere data ( msphr.txt ) ----
number of spheres : 1
Sphere ID 0
radius of sphere : 2.5
refractive index of sphere : 1.33+ 5e-06i
x-coordinate of sphere center : 0
y-coordinate of sphere center : 0
z-coordinate of sphere center : 0
basic sampling number on sphere surface : 32
division number for sphere surface (per PI): 4
limit of order number l : 40
"ipw.txt" No.00 -- plane wave, MKSA system --
wave length of incident beam in vacuum [m]: 1.06e-06
refractive index of surrounding : 1
incident beam power per unit area [W/m^2]: 1000000
x-component of polarization coefficient : 1+ 0i
y-component of polarization coefficient : 0+ 0i
x-component of translation vector [m]: 0
y-component of translation vector [m]: 0
z-component of translation vector [m]: 0
rotation parameter theta [rad]: 0
rotation parameter phi [rad]: 0
---- sphere data ( msphr.txt ), MKSA system ----
number of spheres : 1
Sphere ID 0
radius of sphere [m]: 2.5e-06
refractive index of sphere : 1.33+ 5e-06i
x-coordinate of sphere center [m]: 0
y-coordinate of sphere center [m]: 0
z-coordinate of sphere center [m]: 0
basic sampling number on sphere surface : 32
division number for sphere surface (per PI): 4
limit of order number l : 40
Electromagnetic field at r=( 0, 0,-1.5 )
Ex = 2.40261009451906e-02 -2.36249997098309e-02 I (= 1.47468268462630e+04 -1.45006374841536e+04 I [V/m](MKSA))
Ey = 6.44419191875127e-06 +4.58319376006600e-06 I (= 3.95533934560177e+00 +2.81308918732798e+00 I [V/m](MKSA))
Ez = 1.39150611378504e-09 +3.58023106936773e-09 I (= 8.54083638552759e-04 +2.19748713159991e-03 I [V/m](MKSA))
Hx = -1.97403610922506e-06 +3.04514939414302e-06 I (=-3.21617601699252e-03 +4.96127522887451e-03 I [A/m](MKSA))
Hy = 5.41251335953351e-02 -3.93681179901579e-02 I (= 8.81827620945443e+01 -6.41400612290636e+01 I [A/m](MKSA))
Hz = -1.93622238168653e-09 -7.14579821044568e-09 I (=-3.15456842883638e-06 -1.16422109602268e-05 I [A/m](MKSA))
Real electromagnetic field at t=0.001
Ex = 2.38856417074318e-02 (= 1.46606152689907e+04 [V/m](MKSA))
Ey = 6.47124558322777e-06 (= 3.97194444130590e+00 [V/m](MKSA))
Ez = 1.41270348276209e-09 (= 8.67094236094743e-04 [V/m](MKSA))
Hx = -1.95595131104183e-06 (=-3.18671156397809e-03 [A/m](MKSA))
Hy = 5.38908282730435e-02 (= 8.78010227967237e+01 [A/m](MKSA))
Hz = -1.97854507543901e-09 (=-3.22352219922849e-06 [A/m](MSKA))
Radiation force and torque
Mie coefficients
sphere id 0, F=( -1.3034956917616e-18,-5.7808764512687e-19, 0.012210628168911 )
0, N=( 6.7323751696673e-22, 4.9160192935406e-22, 2.6183443877571e-19 )
0, F=( -4.3479936101715e-30,-1.9282928229197e-30, 4.0730271369639e-14 ) [ N ](MKSA)
0, N=( 2.2456786320046e-39, 1.6398075276265e-39, 8.7338567661934e-37 ) [N m](MKSA)
Surface integral of maxwell stress tensor (verification)
sphere id 0, F=( 3.0988700375299e-18, 2.0267399456991e-18, 0.012210628168911 )
0, N=( -3.0077014991367e-18, 3.2046644526419e-17,-2.0129599180252e-18 )
0, F=( 1.0336717802053e-29, 6.7604767618907e-30, 4.0730271369639e-14 ) [ N ](MKSA)
0, N=( -1.0032612291857e-35, 1.0689609985592e-34,-6.7145115372622e-36 ) [N m](MKSA)
Absorbed energy
Mie coefficients
sphere id 0, P= 7.2644680605785e-06
0, P= 7.2644680605785e-09 [W]
Surface integral of Poynting vector (verification)
sphere id 0, P= 7.2644680604423e-06
0, P= 7.2644680604423e-09 [W]
Volume integral of loss by Joulian heat (verification)
sphere id 0, P= 7.2644680605787e-06
0, P= 7.2644680605787e-09 [W]