A Complete Example for
...
Code Block |
---|
*************************************************************
*** Omega3P V8.0.0 07/01/2009 $ ***
-------------------------------------------------------------
Copyright 2009, Advanced Computations Department
SLAC National Accelerator Laboratory
*************************************************************
Read Mesh: gun-v4.ncdf
Partitioning Method: parmetis
Setup: Max: 0.007, Sum: 0.015, Balance: 1.002
Matching: Max: 0.008, Sum: 0.015, Balance: 1.001
Contraction: Max: 0.011, Sum: 0.022, Balance: 1.001
InitPart: Max: 0.001, Sum: 0.001, Balance: 1.000
Project: Max: 0.000, Sum: 0.001, Balance: 1.000
Initialize: Max: 0.002, Sum: 0.004, Balance: 1.002
K-way: Max: 0.006, Sum: 0.011, Balance: 1.000
Remap: Max: 0.000, Sum: 0.000, Balance: 1.000
Total: Max: 0.035, Sum: 0.070, Balance: 1.000
***********************************************************
* Total Number of Elements read: 24868
* Total Number of Elements used: 24868
* Total Number of DOFs: 152902
***********************************************************
Total Volume of the strucutre is : 0.0001593555890525991
Number of Grad DOFs: 27613
**********************************************************
ARPACK Loop:
Shift = 439.2566356039645
**********************************************************
factorizing the matrix using MUMPS ...
Using ParMETIS for ordering...
Analysis step: 0.761027 seconds
Maximal per-core estimated memory 451 MB
Aggregated estimated memory 875 MB
Maximal per-core estimated memory if OOC 184 MB
Aggregated estimated memory if OOC 365 MB
Factorization step: 39.184121 seconds
ncv=6 nev=2
Number of converged eigenpairs = 2
eigenvalue: 2.726061596672462e+03 Frequency: 2.491200411267457e+09 Residual: 2.54e-11
eigenvalue: 3.545447315320415e+03 Frequency: 2.841033486606133e+09 Residual: 1.25e-09
COMMIT MODE: 0 FREQ = 2491200411.267457 k= 52.21169980638881 norm(v[0]) = 31.04704308504932
COMMIT MODE: 1 FREQ = 2841033486.606133 k= 59.54365890101494 norm(v[1]) = 33.84520753989371
Number of TriSolve: 39. Average time for one TriSolv: 0.287939
Computed Total Energy (normalized by Epsilon0/2): 0.999999999999992
Computed Total Energy (normalized by Epsilon0/2): 0.9999999999999895
|
Some Explanations
The following lines show the number of elements and number of DOFs in the computation:
Code Block |
---|
***********************************************************
* Total Number of Elements read: 24868
* Total Number of Elements used: 24868
* Total Number of DOFs: 152902
***********************************************************
|
The following lines show the memory usage:
Code Block |
---|
Maximal per-core estimated memory 451 MB
Aggregated estimated memory 875 MB
Maximal per-core estimated memory if OOC 184 MB
Aggregated estimated memory if OOC 365 MB
|
If the Maximal per-core memory is larger than what is available, user should either increase the number of cores in the computation or use other options such as Out-of-core (OOC) solver. Note that it may not be good to use excessively large number of cores in the computation (in fact, it may hurt performance or it may fail to get results if doing so). It is often a good idea to have a few thousand elements per core in the parallel computations.
The following lines show the resulting eigen frequencies.
Code Block |
---|
Number of converged eigenpairs = 2
eigenvalue: 2.726061596672462e+03 Frequency: 2.491200411267457e+09 Residual: 2.54e-11
eigenvalue: 3.545447315320415e+03 Frequency: 2.841033486606133e+09 Residual: 1.25e-09
|
multipacting field level scan
Code Block |
---|
h3. 1. the block example for field level scan in the input file:
FieldScales:
{
Type: FieldGradient
ScanToken: 1 // 1: scan, 0: no scan
Minimum: 1.44e+6
Maximum: 1.7e+6
Interval: 0.010e+06
}
h3. 2. Information related all particles with resonant trajectories are in the file ./job/OUTPUT/resonantparticles, the format of resonantparticles shows below:
Field Level ID X Y Z Energy Phase RF Impact_num
16000000 4310 -0.00079 0.100003 -0.00017 10.58931 1.48E-08 19.41 39
16000000 4310 -0.00079 0.100003 0.000168 10.50246 1.52E-08 19.91 40
16000000 4440 -0.00079 0.100003 -0.00017 10.58718 1.48E-08 19.41 39
16000000 4440 -0.00079 0.100003 0.000168 10.50515 1.52E-08 19.91 40
...
17000000 4500 -0.00078 0.100008 -0.00017 15.45694 1.48E-08 19.39 39
17000000 4500 -0.00078 0.100008 0.000175 15.46824 1.51E-08 19.89 40
17000000 4630 -0.00078 0.100008 -0.00017 15.43081 1.48E-08 19.39 39
17000000 4630 -0.00078 0.100008 0.000175 15.38785 1.51E-08 19.89 40
17000000 4760 -0.00078 0.100008 -0.00017 15.36465 1.48E-08 19.39 39
17000000 4760 -0.00078 0.100008 0.000175 15.45242 1.51E-08 19.89 40
...
18000000 10901 -0.00078 0.100005 -0.00017 20.42707 1.47E-08 19.37 38
18000000 10901 -0.00078 0.100006 0.000178 20.46064 1.51E-08 19.87 39
18000000 11021 -0.00078 0.100005 -0.00018 20.49802 1.47E-08 19.37 38
18000000 11021 -0.00078 0.100006 0.000178 20.46099 1.51E-08 19.87 39
18000000 11141 -0.00078 0.100006 -0.00018 20.42072 1.47E-08 19.37 38
18000000 11141 -0.00078 0.100005 0.000177 20.4386 1.51E-08 19.87 39
19000000 3460 -0.00078 0.100016 -0.00017 24.09866 1.47E-08 19.36 39
...
Explanation:
First colum: Field level;
Second colum: particle ID;
Third colum: the x coordinator of particle's impact position;
Fourth colum: the y coordinator of particle's impact position;
Fifth colum: the z coordinator of particle's impact position;
Sixth colum: particle's impact energy in eV;
Seventh colum: particle's impact phase;
Eighth colum: particle's impact time in RF cycle;
Ninth colum: particle's impact number;
For example, from first two lines in the above example output data, we know that at field level 16MV/m, the particle's (ID 4310)last two impacts' (39, 40) impact energies are 10.58931 and 10.50246, we also know their corresponding impact positions. For details of each resonant particle's information, see the file in ./job/OUTPUT/MPParticles\*
h3. 3. Information related all particles with resonant trajectories at a single field level is in the file ./job/OUTPUT/MPParticles*, the format of the file shows below:
Particle ID 47001 total 51 impacts
0 -1.21169903e-01 2.59999000e-01 -1.10533933e-01 2.00000000e+00 4.60000000e-01
1 -1.20286001e-01 2.60014656e-01 -1.10919628e-01 1.65458514e+02 1.23000000e+00
2 -1.21815314e-01 2.60075499e-01 -1.12468337e-01 2.00254221e+02 2.16000000e+00
3 -1.17494040e-01 2.60168428e-01 -1.09433448e-01 5.31140075e+02 3.33000000e+00
4 -1.20069099e-01 2.60014500e-01 -1.11640855e-01 7.07904283e+01 4.15000000e+00
5 -1.21872443e-01 2.60055470e-01 -1.13187720e-01 1.54175170e+02 5.13000000e+00
6 -1.15927695e-01 2.60106057e-01 -1.08433821e-01 5.91315164e+02 6.37000000e+00
7 -1.18458881e-01 2.60129450e-01 -1.10663427e-01 5.39523370e+01 7.16000000e+00
8 -1.20625111e-01 2.60152822e-01 -1.12498558e-01 7.83412829e+01 8.13000000e+00
9 -1.20982466e-01 2.60050997e-01 -1.13368367e-01 2.15465202e+02 9.18000000e+00
10 -1.20145835e-01 2.60199572e-01 -1.13260591e-01 3.81278009e+02 1.02300000e+01
11 -1.21081954e-01 2.60209524e-01 -1.14354168e-01 2.96273552e+02 1.11800000e+01
12 -1.18728446e-01 2.60127100e-01 -1.13069714e-01 5.30929865e+02 1.22700000e+01
13 -1.20875235e-01 2.60133494e-01 -1.15022549e-01 1.49115895e+02 1.31400000e+01
14 -1.18077892e-01 2.60001701e-01 -1.13265186e-01 4.78813806e+02 1.42800000e+01
15 -1.20298504e-01 2.60145823e-01 -1.15279976e-01 1.19531080e+02 1.51400000e+01
16 -1.19283032e-01 2.60261229e-01 -1.14922126e-01 3.71941381e+02 1.62300000e+01
17 -1.20626039e-01 2.60137910e-01 -1.16301477e-01 2.86617076e+02 1.71600000e+01
18 -1.26545724e-01 2.56247874e-01 -1.19998098e-01 2.94209998e+01 1.79000000e+01
19 -1.13254396e-01 2.60410214e-01 -1.11854058e-01 9.17291745e+02 1.83100000e+01
20 -1.16216727e-01 2.60006064e-01 -1.14613063e-01 6.06629977e+01 1.91500000e+01
21 -1.18181022e-01 2.60141588e-01 -1.16514923e-01 7.16122961e+01 2.01300000e+01
22 -1.19455437e-01 2.60197068e-01 -1.17765165e-01 1.78624164e+02 2.11400000e+01
23 -1.16661134e-01 2.60219134e-01 -1.15335359e-01 5.55110182e+02 2.22900000e+01
24 -1.19013362e-01 2.60005150e-01 -1.17530757e-01 1.19042630e+02 2.31300000e+01
25 -1.18060882e-01 2.60144940e-01 -1.16745516e-01 4.04362017e+02 2.42300000e+01
26 -1.19495663e-01 2.60121237e-01 -1.18219764e-01 2.81250204e+02 2.51500000e+01
27 -1.24685415e-01 2.55970645e-01 -1.22821171e-01 4.42499046e+01 2.58900000e+01
28 -1.12511919e-01 2.60104032e-01 -1.12167073e-01 1.01041046e+03 2.63000000e+01
29 -1.15350471e-01 2.60050797e-01 -1.14952432e-01 6.19635677e+01 2.71500000e+01
30 -1.17335755e-01 2.60099282e-01 -1.16863940e-01 6.37315316e+01 2.81300000e+01
31 -1.18941468e-01 2.60010636e-01 -1.18486082e-01 1.41659527e+02 2.91200000e+01
32 -1.15912019e-01 2.60229438e-01 -1.15483606e-01 5.45095116e+02 3.03000000e+01
33 -1.18120098e-01 2.60041808e-01 -1.17699454e-01 1.01173353e+02 3.11400000e+01
34 -1.18760197e-01 2.60142979e-01 -1.18302103e-01 2.97733684e+02 3.21700000e+01
35 -1.17946438e-01 2.60144149e-01 -1.17407052e-01 5.20894363e+02 3.32400000e+01
36 -1.18972829e-01 2.60056485e-01 -1.18449114e-01 3.43538106e+02 3.41700000e+01
37 -1.18249097e-01 2.60303235e-01 -1.17635259e-01 5.53380420e+02 3.52400000e+01
38 -1.18519568e-01 2.60211948e-01 -1.17860002e-01 4.44891396e+02 3.62100000e+01
39 -1.18466696e-01 2.60286228e-01 -1.17725106e-01 4.75773402e+02 3.72200000e+01
40 -1.18702774e-01 2.60114358e-01 -1.17894782e-01 4.67737859e+02 3.82100000e+01
41 -1.18561070e-01 2.60127665e-01 -1.17680567e-01 4.92410647e+02 3.92200000e+01
42 -1.18683590e-01 2.60201748e-01 -1.17790259e-01 4.54332311e+02 4.02100000e+01
43 -1.18590345e-01 2.60203845e-01 -1.17668108e-01 4.88547811e+02 4.12200000e+01
44 -1.18774894e-01 2.60157822e-01 -1.17832953e-01 4.67193251e+02 4.22100000e+01
45 -1.18659842e-01 2.60089707e-01 -1.17682719e-01 5.00522149e+02 4.32200000e+01
46 -1.18760866e-01 2.60163027e-01 -1.17781515e-01 4.61325010e+02 4.42100000e+01
47 -1.18652996e-01 2.60206130e-01 -1.17666973e-01 4.91540232e+02 4.52200000e+01
48 -1.18819730e-01 2.60072496e-01 -1.17824145e-01 4.75822042e+02 4.62100000e+01
49 -1.18659769e-01 2.60214521e-01 -1.17666950e-01 4.91700497e+02 4.72200000e+01
50 -1.18825841e-01 2.60058104e-01 -1.17823183e-01 4.78194137e+02 4.82100000e+01
Particle ID 42001 total 51 impacts
0 -1.21169903e-01 2.59999000e-01 -1.10533933e-01 2.00000000e+00 4.10000000e-01
1 -1.19717870e-01 2.60078718e-01 -1.10396950e-01 2.29748826e+02 1.27000000e+00
2 -1.21626815e-01 2.60171944e-01 -1.12180022e-01 1.55856256e+02 2.16000000e+00
3 -1.18380325e-01 2.60108270e-01 -1.10149169e-01 4.49658586e+02 3.31000000e+00
4 -1.20812024e-01 2.60100620e-01 -1.12234595e-01 9.30324517e+01 4.15000000e+00
5 -1.21070594e-01 2.60130196e-01 -1.12987526e-01 2.49027462e+02 5.19000000e+00
6 -1.20041714e-01 2.60203443e-01 -1.12791026e-01 4.05999735e+02 6.24000000e+00
7 -1.21337140e-01 2.60171838e-01 -1.14070497e-01 2.58355081e+02 7.17000000e+00
8 -1.18376238e-01 2.60032235e-01 -1.12356484e-01 5.66913680e+02 8.28000000e+00
9 -1.20610667e-01 2.60059309e-01 -1.14360934e-01 1.13629350e+02 9.14000000e+00
10 -1.20088283e-01 2.60008309e-01 -1.14447177e-01 3.43602211e+02 1.02100000e+01
11 -1.20769554e-01 2.60066577e-01 -1.15233935e-01 3.28152116e+02 1.11800000e+01
12 -1.19073090e-01 2.60106117e-01 -1.14365911e-01 5.12065738e+02 1.22500000e+01
13 -1.20685281e-01 2.60064625e-01 -1.15835253e-01 2.13176558e+02 1.31500000e+01
14 -1.18277477e-01 2.60087515e-01 -1.14257214e-01 5.43008579e+02 1.42700000e+01
15 -1.20476288e-01 2.60037643e-01 -1.16257849e-01 1.54201757e+02 1.51300000e+01
16 -1.17388872e-01 2.60165414e-01 -1.14039364e-01 5.26629438e+02 1.62900000e+01
17 -1.19734173e-01 2.60183272e-01 -1.16151330e-01 1.15992116e+02 1.71400000e+01
18 -1.19000807e-01 2.60228824e-01 -1.15888864e-01 3.61963392e+02 1.82200000e+01
19 -1.20261893e-01 2.60015653e-01 -1.17098927e-01 3.21737896e+02 1.91600000e+01
20 -1.25827427e-01 2.56157759e-01 -1.21187823e-01 4.54170936e+01 1.98900000e+01
21 -1.13109491e-01 2.60288425e-01 -1.12227070e-01 9.48625096e+02 2.03000000e+01
22 -1.15987437e-01 2.60022222e-01 -1.14970732e-01 6.14881131e+01 2.11500000e+01
23 -1.17935509e-01 2.60136698e-01 -1.16836908e-01 7.20874106e+01 2.21300000e+01
24 -1.19248827e-01 2.60019236e-01 -1.18281619e-01 1.77513855e+02 2.31300000e+01
25 -1.16472021e-01 2.60233169e-01 -1.15657858e-01 5.66740418e+02 2.42900000e+01
26 -1.18787381e-01 2.60020176e-01 -1.17887913e-01 1.23569898e+02 2.51300000e+01
27 -1.17607739e-01 2.60150432e-01 -1.16757938e-01 4.19837359e+02 2.62400000e+01
28 -1.19283114e-01 2.60115125e-01 -1.18495326e-01 2.39178714e+02 2.71400000e+01
29 -1.24366967e-01 2.55948066e-01 -1.23175697e-01 2.92332483e+01 2.79000000e+01
30 -1.12942763e-01 2.60015484e-01 -1.12666143e-01 9.64364448e+02 2.83000000e+01
31 -1.15612580e-01 2.60008597e-01 -1.15288293e-01 6.06920203e+01 2.91500000e+01
32 -1.17568462e-01 2.60126685e-01 -1.17163137e-01 7.07355384e+01 3.01300000e+01
33 -1.18930168e-01 2.60136415e-01 -1.18641507e-01 1.79147123e+02 3.11300000e+01
34 -1.16107438e-01 2.60013381e-01 -1.15898158e-01 5.97637027e+02 3.22900000e+01
35 -1.18216705e-01 2.60153276e-01 -1.18029486e-01 1.13940087e+02 3.31400000e+01
36 -1.17987820e-01 2.60235956e-01 -1.17675404e-01 3.56535906e+02 3.42100000e+01
37 -1.18790059e-01 2.60165498e-01 -1.18428150e-01 3.73035619e+02 3.51800000e+01
38 -1.18359891e-01 2.60024533e-01 -1.17914088e-01 5.46924437e+02 3.62300000e+01
39 -1.18445276e-01 2.60216168e-01 -1.17952010e-01 4.44574962e+02 3.72100000e+01
40 -1.18365373e-01 2.60197599e-01 -1.17731690e-01 4.80906956e+02 3.82200000e+01
41 -1.18746611e-01 2.60042218e-01 -1.18062952e-01 4.42132549e+02 3.92000000e+01
42 -1.18550546e-01 2.60066400e-01 -1.17759540e-01 4.99081539e+02 4.02200000e+01
43 -1.18635329e-01 2.60207745e-01 -1.17827420e-01 4.53189934e+02 4.12100000e+01
44 -1.18539204e-01 2.60217526e-01 -1.17692910e-01 4.86772805e+02 4.22200000e+01
45 -1.18733981e-01 2.60163207e-01 -1.17853133e-01 4.65522750e+02 4.32100000e+01
46 -1.18614947e-01 2.60074546e-01 -1.17686311e-01 4.98291037e+02 4.42200000e+01
47 -1.18708470e-01 2.60208878e-01 -1.17774663e-01 4.54374656e+02 4.52100000e+01
48 -1.18611960e-01 2.60157058e-01 -1.17650795e-01 4.90332007e+02 4.62200000e+01
49 -1.18760267e-01 2.60185775e-01 -1.17791781e-01 4.61096091e+02 4.72100000e+01
50 -1.18662943e-01 2.60149916e-01 -1.17668177e-01 4.97097101e+02 4.82200000e+01
...
Explanation: This file records all the resonant particles' all impact information at this field level. For example, for particle ID 47001, it has 51 impacts,
Impact number impact position (x, y, z) impact energy(eV) impact time (in RF cycle)
0 -1.21169903e-01 2.59999000e-01 -1.10533933e-01 2.00000000e+00 4.60000000e-01
1 -1.20286001e-01 2.60014656e-01 -1.10919628e-01 1.65458514e+02 1.23000000e+00
2 -1.21815314e-01 2.60075499e-01 -1.12468337e-01 2.00254221e+02 2.16000000e+00
...
48 -1.18819730e-01 2.60072496e-01 -1.17824145e-01 4.75822042e+02 4.62100000e+01
49 -1.18659769e-01 2.60214521e-01 -1.17666950e-01 4.91700497e+02 4.72200000e+01
50 -1.18825841e-01 2.60058104e-01 -1.17823183e-01 4.78194137e+02 4.82100000e+01
According to above information, we have the conclusion: at this field level, there are resonant trajectories: one point first order and the impact energies is around 470eV \~500eV
|
A Complete Example for particles trajectories at a given field level
Code Block |
---|
h3. 1. the block example for particles trajectories in the input file:
FieldScales:
{
Type: FieldGradient
Scale: 0.15e+06 //field scale for particle trajectory
}
By default, the code will save all particles locations at each time step, there are other options for only save information at some time steps, see the examples below:
ParticlesTrajectories: // record particles' trajectory, only for running single
//field case
{ ParticleFile: p // file name
Skip: 10 // write file each 10 steps
Start: 10 // start time step for writing file
Stop: 100000 // stop time step for writing file
}
h3. 2. Output files saved in the directories ./job/PARTICLE/*/p*_ts*, the file is in ncdf format, can be viewed through paraview. The particle's ID, location, momentum... is saved in these files.
|
...