-
Notifications
You must be signed in to change notification settings - Fork 0
/
example1.log
96 lines (89 loc) · 5.15 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
88
89
90
91
92
93
94
95
96
./example1.out ex.dat
"ipw.txt" No.00 -- plane wave --
wave length of incident beam in vacuum : 1.064
refractive index of surrounding : 1.33
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 ---------
sphere id 0
basic sampling number on sphere surface : 32
division number of sphere surface (per PI): 4
limit of order number l : 40
x-coordinate of sphere center : 0
y-coordinate of sphere center : 0
z-coordinate of sphere center : 0
layer id 0
radius of layer : 0.5
refractive index of layer : 1+ 0I
layer id 1
radius of layer : 1.5
refractive index of layer : 1.6+ 1e-06I
layer id 2
radius of layer : 2.5
refractive index of layer : 1.5+ 0.1I
"ipw.txt" No.00 -- plane wave, MKSA system --
wave length of incident beam in vacuum [m]: 1.064e-06
refractive index of surrounding : 1.33
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, MKSA system ---------
sphere id 0
basic sampling number on sphere surface : 32
division number of sphere surface (per PI): 4
limit of order number l : 40
x-coordinate of sphere center [m]: 0
y-coordinate of sphere center [m]: 0
z-coordinate of sphere center [m]: 0
layer id 0
radius of layer [m]: 5e-07
refractive index of layer : 1+ 0I
layer id 1
radius of layer [m]: 1.5e-06
refractive index of layer : 1.6+ 1e-06I
layer id 2
radius of layer [m]: 2.5e-06
refractive index of layer : 1.5+ 0.1I
Electromagnetic field at r=( 0.5, 0,-1 )
Ex = 2.06142583952747e-02 +1.46453556424817e-03 I (= 1.26526938271310e+04 +8.98907917911043e+02 I [V/m](MKSA))
Ey = 2.14900518082931e-16 +1.68911738263534e-16 I (= 1.31902414651911e-10 +1.03675255596237e-10 I [V/m](MKSA))
Ez = 1.11869668214660e-03 -1.11656045808259e-03 I (= 6.86637682191504e+02 -6.85326502884943e+02 I [V/m](MKSA))
Hx = 7.46360281005306e-17 -1.43360490406160e-16 I (= 1.21599905117612e-13 -2.33568458486597e-13 I [A/m](MKSA))
Hy = 3.30334352405617e-02 +7.04875501034219e-03 I (= 5.38193509647992e+01 +1.14841044234079e+01 I [A/m](MKSA))
Hz = 2.87872059778360e-16 -3.11430636027874e-16 I (= 4.69012299367136e-13 -5.07394843421977e-13 I [A/m](MKSA))
Real electromagnetic field at t=0.001
Ex = 2.06225473638429e-02 (= 1.26577814601385e+04 [V/m](MKSA))
Ey = 2.15894231226761e-16 (= 1.32512339487420e-10 [V/m](MKSA))
Ez = 1.11208364702845e-03 (= 6.82578709658340e+02 [V/m](MKSA))
Hx = 7.37881522293354e-17 (= 1.20218512938623e-13 [A/m](MKSA))
Hy = 3.30744836857694e-02 (= 5.38862286801535e+01 [A/m](MKSA))
Hz = 2.86027975588111e-16 (= 4.66007846045194e-13 [A/m](MKSA))
Radiation force and torque
Mie coefficient
sphere id 0, F=( -1.9285707478781e-19,-2.3363916771048e-19, 0.028806215438706 )
0, N=( 1.0642723522147e-18, 7.5150891511086e-18,-4.160776480023e-18 )
0, F=( -6.4330195654158e-31,-7.7933637580195e-31, 9.6087191888948e-14 ) [ N ](MKSA)
0, N=( 3.5500304421089e-36, 2.506763913023e-35,-1.3878856418806e-35 ) [N m](MKSA)
Surface integral
sphere id 0, F=( 8.5960609256511e-21,-2.8893399227678e-19, 0.028806215438705 )
0, N=( 1.4785104850931e-18, 2.8951586137152e-17,-3.0861267648354e-18 )
0, F=( 2.867337284933e-32,-9.6378005705794e-31, 9.6087191888947e-14 ) [ N ](MKSA)
0, N=( 4.9317801220106e-36, 9.6572096343905e-35,-1.0294210819791e-35 ) [N m](MKSA)
Absorbed energy
Mie coefficients
sphere id 0, P= 0.021096112253686
0, P= 2.1096112253686e-05 [W]
Surface integral of Poynting vector (verification)
sphere id 0, P= 0.021096112253686
0, P= 2.1096112253686e-05 [W]