Test of examples

As all pushers and geometries have change, I'm testing all examples,
   adding a README.txt and a reference output for each case.
This commit is contained in:
Jorge Gonzalez 2022-12-10 10:37:14 +01:00
commit 2850f37551
15 changed files with 21137 additions and 12 deletions

Binary file not shown.

View file

@ -47,6 +47,7 @@
\newglossaryentry{git}{name={Git},description={Git is a distributed version-control system for tracking changes in a set of files}}
\newglossaryentry{gitlab}{name={GitLab},description={GitLab is a web-based lifecycle tool that provides a \Gls{git}-repository manager}}
\newglossaryentry{gmsh}{name={Gmsh},description={A three-dimensional finite element mesh generator with built-in pre- and post-processing facilities.}}
\newglossaryentry{gnuplot}{name={Gnuplot},description={A portable command-line driven graphing utility for Linux, OS/2, MS Windows, OSX, VMS, and many other platforms.}}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\bibliography{bibliography}
@ -287,7 +288,7 @@
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Gmsh}
Although \Gls{gmsh}\cite{gmshURL} is not required to compile and run \Gls{fpakc}, it is the default tool to generate finite element meshes and post-processing.
Right now, the only \acrshort{io} format available in \Gls{fpakc} is the v2.0 .msh format.
Right now, the only \acrshort{io} format available in \Gls{fpakc} is the v2.0 \lstinline|.msh| format.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Installation steps}
@ -340,7 +341,7 @@ make
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Mesh file}
\Gls{fpakc} accepts right now the version 2.0 of \Gls{gmsh} mesh format .msh in ASCII format.
\Gls{fpakc} accepts right now the version 2.0 of \Gls{gmsh} mesh format \lstinline|.msh| in ASCII format.
This file contains information about the nodes, edges and volumes that define the finite element mesh used by \Gls{fpakc} to scatter particle properties and compute the self-consistent electromagnetic field.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@ -743,6 +744,8 @@ make
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\chapter{Example runs}\label{ch:exampleRuns}
This chapter presents a description of the different example files distributed with \acrshort{fpakc}.
All examples in the repository have a \lstinline|README.txt| file and a reference output.
Plotting of the output is done with \Gls{gnuplot} or \Gls{gmsh}.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{1D Emissive Cathode (1D\_Cathode)}
@ -754,7 +757,7 @@ make
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{0D \ce{Ar}-\ce{Ar+} Elastic Collision (0D\_Argon)}
Test case to check the 0D geometry that includes the elastic collision between \ce{Ar} and \ce{Ar+}.
Initial states are read from the Argon\_Initial.dat and Argon+\_Initial.dat
Initial states are read from the \lstinline|Argon_Initial.dat| and \lstinline|Argon+_Initial.dat|.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{ALPHIE Grid system (ALPHIE\_Grid)}
@ -765,7 +768,7 @@ make
\section{Flow around cylinder (cylFlow)}
Simple case of neutral Argon flow around a cylinder in a 2D axial-symmetry geometry.
Elastic collisions between argon particles are included.
Two cases are presented here: one in which the same mesh (meshSingle.msh) for scattering and collisions is used (input.json) and a second one (inputDualMesh.json) in which a mesh is used for scattering (mesh.msh) and a second one is used only for collisions (meshColl.msh).
Two cases are presented here: one in which the same mesh (\lstinline|meshSingle.msh|) for scattering and collisions is used (\lstinline|input.json|) and a second one (\lstinline|inputDualMesh.json|) in which a mesh is used for scattering (\lstinline|mesh.msh|) and a second one is used only for collisions (\lstinline|meshColl.msh|).
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\printglossaries

View file

@ -2,7 +2,7 @@
{
"output": {
"path": "./runs/1D_Cathode/",
"triggerOutput": 1000,
"triggerOutput": 100,
"cpuTime": false,
"numColl": false,
"EMField": true,
@ -39,7 +39,7 @@
],
"solver": {
"tau": [1.0e-11],
"finalTime": 3.0e-7,
"finalTime": 1.0e-8,
"pusher": ["Electrostatic"],
"EMSolver": "Electrostatic"
},

View file

@ -2,7 +2,7 @@
{
"output": {
"path": "./runs/1D_Cathode/",
"triggerOutput": 1000,
"triggerOutput": 100,
"cpuTime": false,
"numColl": false,
"EMField": true,
@ -15,7 +15,7 @@
},
"geometry": {
"dimension": 1,
"type": "Cart",
"type": "Rad",
"meshType": "gmsh2",
"meshFile": "mesh.msh"
},
@ -39,7 +39,7 @@
],
"solver": {
"tau": [1.0e-11],
"finalTime": 3.0e-7,
"finalTime": 1.0e-8,
"pusher": ["Electrostatic"],
"EMSolver": "Electrostatic"
},

View file

@ -0,0 +1,185 @@
$MeshFormat
2.2 0 8
$EndMeshFormat
$NodeData
1
"Potential (V)"
1
1.0000000000000000E-008
3
1000
1
51
1 -7.2180834934653183E-012
2 -1.6164415181357834
3 -0.36980461527280883
4 -0.57492573345192088
5 -0.71403100967357258
6 -0.82100545766928923
7 -0.90752842743341133
8 -0.97972921034487059
9 -1.0411791818607907
10 -1.0939318279574932
11 -1.1399404313046071
12 -1.1807139724083482
13 -1.2174163304531369
14 -1.2504990809463861
15 -1.2803540949915220
16 -1.3073500616819811
17 -1.3316282169084139
18 -1.3537385275137412
19 -1.3739696197563589
20 -1.3927592194128697
21 -1.4103076231740377
22 -1.4266496533259061
23 -1.4418485519615323
24 -1.4562354448986332
25 -1.4697852468324024
26 -1.4825221356246889
27 -1.4946048104025964
28 -1.5059524644639768
29 -1.5166545901603037
30 -1.5269837565446036
31 -1.5366931422515575
32 -1.5457467079858906
33 -1.5541298965687140
34 -1.5618551829370724
35 -1.5689172611267244
36 -1.5752150587183313
37 -1.5808462624732060
38 -1.5860612571248669
39 -1.5910404498436017
40 -1.5957279636246224
41 -1.5998828959091791
42 -1.6032394878414955
43 -1.6059267946667977
44 -1.6082423003097928
45 -1.6102897646308323
46 -1.6120678175170864
47 -1.6134998424675298
48 -1.6146588662864141
49 -1.6155317189248377
50 -1.6160774066120194
51 -1.6163242439958503
$EndNodeData
$ElementData
1
"Electric Field (V m^-1)"
1
1.0000000000000000E-008
3
1000
3
50
3 9245.1153816612405 0.0000000000000000 0.0000000000000000
4 5128.0279544843424 0.0000000000000000 0.0000000000000000
5 3477.6319055457411 0.0000000000000000 0.0000000000000000
6 2674.3611998963270 0.0000000000000000 0.0000000000000000
7 2163.0742441058046 0.0000000000000000 0.0000000000000000
8 1805.0195727920382 0.0000000000000000 0.0000000000000000
9 1536.2492879082888 0.0000000000000000 0.0000000000000000
10 1318.8161524257080 0.0000000000000000 0.0000000000000000
11 1150.2150836799428 0.0000000000000000 0.0000000000000000
12 1019.3385275948289 0.0000000000000000 0.0000000000000000
13 917.55895112089058 0.0000000000000000 0.0000000000000000
14 827.06876233228445 0.0000000000000000 0.0000000000000000
15 746.37535112673595 0.0000000000000000 0.0000000000000000
16 674.89916726183060 0.0000000000000000 0.0000000000000000
17 606.95388066159535 0.0000000000000000 0.0000000000000000
18 552.75776513387154 0.0000000000000000 0.0000000000000000
19 505.77730606825133 0.0000000000000000 0.0000000000000000
20 469.73999141591514 0.0000000000000000 0.0000000000000000
21 438.71009403214083 0.0000000000000000 0.0000000000000000
22 408.55075379944572 0.0000000000000000 0.0000000000000000
23 379.97246589319053 0.0000000000000000 0.0000000000000000
24 359.67232342854714 0.0000000000000000 0.0000000000000000
25 338.74504834466700 0.0000000000000000 0.0000000000000000
26 318.42221980693205 0.0000000000000000 0.0000000000000000
27 302.06686944643235 0.0000000000000000 0.0000000000000000
28 283.69135153333235 0.0000000000000000 0.0000000000000000
29 267.55314240706213 0.0000000000000000 0.0000000000000000
30 258.22915960642655 0.0000000000000000 0.0000000000000000
31 242.73464267283171 0.0000000000000000 0.0000000000000000
32 226.33914335739505 0.0000000000000000 0.0000000000000000
33 209.57971456971421 0.0000000000000000 0.0000000000000000
34 193.13215920815529 0.0000000000000000 0.0000000000000000
35 176.55195474056683 0.0000000000000000 0.0000000000000000
36 157.44493979054459 0.0000000000000000 0.0000000000000000
37 140.78009387281193 0.0000000000000000 0.0000000000000000
38 130.37486629147924 0.0000000000000000 0.0000000000000000
39 124.47981796785177 0.0000000000000000 0.0000000000000000
40 117.18784452503036 0.0000000000000000 0.0000000000000000
41 103.87330711348883 0.0000000000000000 0.0000000000000000
42 83.914798307876424 0.0000000000000000 0.0000000000000000
43 67.182670633004349 0.0000000000000000 0.0000000000000000
44 57.887641075013939 0.0000000000000000 0.0000000000000000
45 51.186608025777616 0.0000000000000000 0.0000000000000000
46 44.451322156171258 0.0000000000000000 0.0000000000000000
47 35.800623760933746 0.0000000000000000 0.0000000000000000
48 28.975595471990374 0.0000000000000000 0.0000000000000000
49 21.821315960497699 0.0000000000000000 0.0000000000000000
50 13.642192179482622 0.0000000000000000 0.0000000000000000
51 6.1709345957464548 0.0000000000000000 0.0000000000000000
52 2.9318534983186093 0.0000000000000000 0.0000000000000000
$EndElementData
$NodeData
1
"Magnetic Field (T)"
1
1.0000000000000000E-008
3
1000
3
51
1 0.0000000000000000 0.0000000000000000 0.0000000000000000
2 0.0000000000000000 0.0000000000000000 0.0000000000000000
3 0.0000000000000000 0.0000000000000000 0.0000000000000000
4 0.0000000000000000 0.0000000000000000 0.0000000000000000
5 0.0000000000000000 0.0000000000000000 0.0000000000000000
6 0.0000000000000000 0.0000000000000000 0.0000000000000000
7 0.0000000000000000 0.0000000000000000 0.0000000000000000
8 0.0000000000000000 0.0000000000000000 0.0000000000000000
9 0.0000000000000000 0.0000000000000000 0.0000000000000000
10 0.0000000000000000 0.0000000000000000 0.0000000000000000
11 0.0000000000000000 0.0000000000000000 0.0000000000000000
12 0.0000000000000000 0.0000000000000000 0.0000000000000000
13 0.0000000000000000 0.0000000000000000 0.0000000000000000
14 0.0000000000000000 0.0000000000000000 0.0000000000000000
15 0.0000000000000000 0.0000000000000000 0.0000000000000000
16 0.0000000000000000 0.0000000000000000 0.0000000000000000
17 0.0000000000000000 0.0000000000000000 0.0000000000000000
18 0.0000000000000000 0.0000000000000000 0.0000000000000000
19 0.0000000000000000 0.0000000000000000 0.0000000000000000
20 0.0000000000000000 0.0000000000000000 0.0000000000000000
21 0.0000000000000000 0.0000000000000000 0.0000000000000000
22 0.0000000000000000 0.0000000000000000 0.0000000000000000
23 0.0000000000000000 0.0000000000000000 0.0000000000000000
24 0.0000000000000000 0.0000000000000000 0.0000000000000000
25 0.0000000000000000 0.0000000000000000 0.0000000000000000
26 0.0000000000000000 0.0000000000000000 0.0000000000000000
27 0.0000000000000000 0.0000000000000000 0.0000000000000000
28 0.0000000000000000 0.0000000000000000 0.0000000000000000
29 0.0000000000000000 0.0000000000000000 0.0000000000000000
30 0.0000000000000000 0.0000000000000000 0.0000000000000000
31 0.0000000000000000 0.0000000000000000 0.0000000000000000
32 0.0000000000000000 0.0000000000000000 0.0000000000000000
33 0.0000000000000000 0.0000000000000000 0.0000000000000000
34 0.0000000000000000 0.0000000000000000 0.0000000000000000
35 0.0000000000000000 0.0000000000000000 0.0000000000000000
36 0.0000000000000000 0.0000000000000000 0.0000000000000000
37 0.0000000000000000 0.0000000000000000 0.0000000000000000
38 0.0000000000000000 0.0000000000000000 0.0000000000000000
39 0.0000000000000000 0.0000000000000000 0.0000000000000000
40 0.0000000000000000 0.0000000000000000 0.0000000000000000
41 0.0000000000000000 0.0000000000000000 0.0000000000000000
42 0.0000000000000000 0.0000000000000000 0.0000000000000000
43 0.0000000000000000 0.0000000000000000 0.0000000000000000
44 0.0000000000000000 0.0000000000000000 0.0000000000000000
45 0.0000000000000000 0.0000000000000000 0.0000000000000000
46 0.0000000000000000 0.0000000000000000 0.0000000000000000
47 0.0000000000000000 0.0000000000000000 0.0000000000000000
48 0.0000000000000000 0.0000000000000000 0.0000000000000000
49 0.0000000000000000 0.0000000000000000 0.0000000000000000
50 0.0000000000000000 0.0000000000000000 0.0000000000000000
51 0.0000000000000000 0.0000000000000000 0.0000000000000000
$EndNodeData

View file

@ -0,0 +1,247 @@
$MeshFormat
2.2 0 8
$EndMeshFormat
$NodeData
1
"Electron density (m^-3)"
1
1.0000000000000000E-008
3
1000
1
51
1 1.663355E+017
2 1.635218E+014
3 6.971311E+016
4 2.630990E+016
5 1.405426E+016
6 9.465391E+015
7 6.637877E+015
8 4.988303E+015
9 4.189069E+015
10 3.156337E+015
11 2.494143E+015
12 1.855358E+015
13 1.740443E+015
14 1.546003E+015
15 1.316731E+015
16 1.372680E+015
17 9.737899E+014
18 9.388176E+014
19 6.512445E+014
20 5.832764E+014
21 5.692614E+014
22 5.935914E+014
23 3.292208E+014
24 4.143419E+014
25 4.100487E+014
26 2.728780E+014
27 3.714897E+014
28 3.455677E+014
29 9.442239E+013
30 3.445456E+014
31 3.018612E+014
32 3.256579E+014
33 3.148337E+014
34 2.986174E+014
35 3.894979E+014
36 3.315681E+014
37 1.927240E+014
38 8.916509E+013
39 1.257294E+014
40 2.430101E+014
41 4.270402E+014
42 3.345209E+014
43 1.510759E+014
44 1.256626E+014
45 1.136905E+014
46 1.909151E+014
47 1.133460E+014
48 1.373180E+014
49 1.567058E+014
50 1.725501E+014
51 8.718697E+012
$EndNodeData
$NodeData
1
"Electron velocity (m s^-1)"
1
1.0000000000000000E-008
3
1000
3
51
1 6.822573E+002 7.341184E+002 -1.294679E+003
2 2.275900E+005 -1.156323E+005 -3.994264E+004
3 6.011992E+003 -2.085702E+003 -6.955356E+003
4 2.674256E+003 -1.413909E+003 1.028663E+003
5 -6.262168E+003 -5.034055E+003 1.381887E+004
6 1.057957E+003 8.436700E+002 -1.017065E+004
7 1.182090E+004 1.797951E+004 6.390618E+003
8 1.868273E+004 1.131158E+004 1.107675E+004
9 1.821590E+004 5.287286E+003 9.958486E+003
10 5.545844E+003 -1.134700E+003 -1.747390E+004
11 4.595683E+004 1.909957E+004 4.055354E+004
12 5.242670E+004 -1.209486E+003 1.046959E+004
13 -6.508610E+003 4.362341E+004 -8.868030E+003
14 -5.779804E+003 7.888582E+003 -1.190795E+004
15 -3.256102E+003 -2.423327E+004 1.680404E+004
16 2.026515E+004 5.900866E+003 4.573129E+004
17 1.881938E+004 -5.907004E+004 4.787712E+003
18 2.471440E+004 3.126628E+004 3.276285E+004
19 1.927866E+004 -1.774386E+004 1.728677E+004
20 4.357578E+004 -1.754746E+004 2.441953E+004
21 -2.293113E+004 5.501146E+003 3.536741E+004
22 2.289073E+004 2.308340E+004 5.214795E+004
23 1.380393E+005 8.134296E+004 6.839466E+004
24 1.059159E+005 7.461802E+003 -9.633646E+003
25 1.093199E+004 -1.991657E+004 -1.139848E+005
26 2.642203E+004 -2.481740E+004 3.143777E+004
27 9.005040E+004 8.070628E+003 -3.359234E+004
28 1.241325E+005 5.096109E+004 7.419358E+002
29 1.346489E+005 -1.573873E+005 6.203455E+003
30 5.610495E+004 -3.368553E+004 1.507073E+005
31 3.535163E+004 1.988632E+004 6.165400E+004
32 9.376597E+004 1.073849E+004 1.451017E+004
33 1.084345E+005 3.421118E+004 2.083726E+004
34 7.746595E+004 3.203130E+004 6.852780E+003
35 1.133766E+005 2.337750E+003 -4.867282E+004
36 1.025276E+005 1.014552E+005 -8.405875E+004
37 7.241202E+004 2.972254E+002 1.876103E+004
38 4.204461E+004 9.161183E+004 -8.625228E+004
39 2.000094E+004 1.616608E+005 -4.124777E+004
40 1.065082E+005 1.060744E+005 -9.625557E+004
41 2.353957E+005 -3.660256E+003 -6.119245E+003
42 2.043649E+005 1.172239E+003 -2.283045E+004
43 1.632956E+005 -3.766183E+004 -5.714495E+004
44 2.139853E+005 -1.091968E+005 -7.473445E+004
45 6.523701E+004 -5.437454E+004 2.092413E+005
46 6.804248E+004 2.190917E+004 -3.051095E+004
47 1.247063E+005 9.312863E+004 -2.110944E+004
48 1.488446E+005 4.703577E+003 1.162251E+005
49 1.444633E+005 2.511461E+004 7.188944E+004
50 1.343558E+005 -1.246722E+004 2.197029E+005
51 1.789699E+005 -1.338478E+005 -5.783810E+004
$EndNodeData
$NodeData
1
"Electron pressure (Pa)"
1
1.0000000000000000E-008
3
1000
1
51
1 5.385580E-003
2 1.002571E-006
3 2.356496E-003
4 8.962727E-004
5 4.794124E-004
6 3.503273E-004
7 2.442948E-004
8 1.546660E-004
9 1.322598E-004
10 9.957574E-005
11 8.063749E-005
12 5.468381E-005
13 5.120468E-005
14 4.099773E-005
15 4.147335E-005
16 4.439409E-005
17 3.060211E-005
18 3.239301E-005
19 1.928871E-005
20 1.720661E-005
21 1.615443E-005
22 2.301910E-005
23 9.019928E-006
24 8.246059E-006
25 9.022171E-006
26 5.090623E-006
27 1.430550E-005
28 1.078930E-005
29 2.237441E-006
30 8.720790E-006
31 7.185960E-006
32 1.040168E-005
33 8.947370E-006
34 8.505472E-006
35 1.011847E-005
36 1.009116E-005
37 6.487404E-006
38 9.794017E-007
39 1.137051E-006
40 4.115733E-006
41 1.110521E-005
42 5.195716E-006
43 2.249886E-006
44 2.702675E-006
45 1.909427E-006
46 4.939937E-006
47 1.983386E-006
48 2.345377E-006
49 4.284590E-006
50 4.071858E-006
51 4.700991E-008
$EndNodeData
$NodeData
1
"Electron temperature (K)"
1
1.0000000000000000E-008
3
1000
1
51
1 2.345117E+003
2 4.440750E+002
3 2.448325E+003
4 2.467391E+003
5 2.470690E+003
6 2.680725E+003
7 2.665642E+003
8 2.245737E+003
9 2.286795E+003
10 2.285004E+003
11 2.341707E+003
12 2.134754E+003
13 2.130918E+003
14 1.920730E+003
15 2.281334E+003
16 2.342463E+003
17 2.276161E+003
18 2.499120E+003
19 2.145241E+003
20 2.136672E+003
21 2.055403E+003
22 2.808779E+003
23 1.984416E+003
24 1.441466E+003
25 1.593648E+003
26 1.351199E+003
27 2.789158E+003
28 2.261399E+003
29 1.716301E+003
30 1.833268E+003
31 1.724227E+003
32 2.313443E+003
33 2.058406E+003
34 2.063004E+003
35 1.881598E+003
36 2.204374E+003
37 2.438102E+003
38 7.955782E+002
39 6.550284E+002
40 1.226704E+003
41 1.883539E+003
42 1.124965E+003
43 1.078654E+003
44 1.557774E+003
45 1.216454E+003
46 1.874123E+003
47 1.267412E+003
48 1.237093E+003
49 1.980347E+003
50 1.709206E+003
51 3.905303E+002
$EndNodeData

View file

@ -0,0 +1,185 @@
$MeshFormat
2.2 0 8
$EndMeshFormat
$NodeData
1
"Potential (V)"
1
1.0000000000000000E-008
3
1000
1
51
1 -2.0125602360500233E-012
2 -2.1476540511245403
3 -0.92254600934286035
4 -1.1693722452494442
5 -1.3055892978798838
6 -1.4090596211031350
7 -1.4933143808451328
8 -1.5638362561671308
9 -1.6231896643635058
10 -1.6740175110100952
11 -1.7185366931231973
12 -1.7583784348355216
13 -1.7946512154913346
14 -1.8278278087327231
15 -1.8580164245432433
16 -1.8852831251773992
17 -1.9099218127680115
18 -1.9323133093677307
19 -1.9527025003199079
20 -1.9712632206620335
21 -1.9881621423695754
22 -2.0035781018067564
23 -2.0176479044718039
24 -2.0305087603061467
25 -2.0422748378993707
26 -2.0530348175074842
27 -2.0628746022811000
28 -2.0718874967016756
29 -2.0801381491382527
30 -2.0876833658364853
31 -2.0945891263065528
32 -2.1008985562257543
33 -2.1066381420356524
34 -2.1118556546612779
35 -2.1165897969031175
36 -2.1208699819402956
37 -2.1247303076336199
38 -2.1281989125689238
39 -2.1312983972494055
40 -2.1340576795392843
41 -2.1364955581097296
42 -2.1386339190704353
43 -2.1405011196666606
44 -2.1421213004925956
45 -2.1435109859394803
46 -2.1446785599633511
47 -2.1456400594108116
48 -2.1464001968748745
49 -2.1469649987104309
50 -2.1473539773488697
51 -2.1475805166422961
$EndNodeData
$ElementData
1
"Electric Field (V m^-1)"
1
1.0000000000000000E-008
3
1000
3
50
3 23063.650233574754 0.0000000000000000 0.0000000000000000
4 6170.6558976724700 0.0000000000000000 0.0000000000000000
5 3405.4263157653404 0.0000000000000000 0.0000000000000000
6 2586.7580805845814 0.0000000000000000 0.0000000000000000
7 2106.3689935526277 0.0000000000000000 0.0000000000000000
8 1763.0468830553750 0.0000000000000000 0.0000000000000000
9 1483.8352049193140 0.0000000000000000 0.0000000000000000
10 1270.6961661725816 0.0000000000000000 0.0000000000000000
11 1112.9795528295804 0.0000000000000000 0.0000000000000000
12 996.04354280937548 0.0000000000000000 0.0000000000000000
13 906.81951639648162 0.0000000000000000 0.0000000000000000
14 829.41483103577525 0.0000000000000000 0.0000000000000000
15 754.71539526132392 0.0000000000000000 0.0000000000000000
16 681.66751585425504 0.0000000000000000 0.0000000000000000
17 615.96718976609407 0.0000000000000000 0.0000000000000000
18 559.78741499368221 0.0000000000000000 0.0000000000000000
19 509.72977380725575 0.0000000000000000 0.0000000000000000
20 464.01800855624742 0.0000000000000000 0.0000000000000000
21 422.47304269137754 0.0000000000000000 0.0000000000000000
22 385.39898593210717 0.0000000000000000 0.0000000000000000
23 351.74506662853850 0.0000000000000000 0.0000000000000000
24 321.52139585948811 0.0000000000000000 0.0000000000000000
25 294.15193983097186 0.0000000000000000 0.0000000000000000
26 268.99949020263892 0.0000000000000000 0.0000000000000000
27 245.99461933938235 0.0000000000000000 0.0000000000000000
28 225.32236051344728 0.0000000000000000 0.0000000000000000
29 206.26631091357709 0.0000000000000000 0.0000000000000000
30 188.63041745502420 0.0000000000000000 0.0000000000000000
31 172.64401175097743 0.0000000000000000 0.0000000000000000
32 157.73574797938116 0.0000000000000000 0.0000000000000000
33 143.48964524685633 0.0000000000000000 0.0000000000000000
34 130.43781564009393 0.0000000000000000 0.0000000000000000
35 118.35355604550293 0.0000000000000000 0.0000000000000000
36 107.00462592970442 0.0000000000000000 0.0000000000000000
37 96.508142333755885 0.0000000000000000 0.0000000000000000
38 86.715123382567057 0.0000000000000000 0.0000000000000000
39 77.487117011718453 0.0000000000000000 0.0000000000000000
40 68.982057246684732 0.0000000000000000 0.0000000000000000
41 60.946964260881018 0.0000000000000000 0.0000000000000000
42 53.459024017616187 0.0000000000000000 0.0000000000000000
43 46.680014905952270 0.0000000000000000 0.0000000000000000
44 40.504520648464592 0.0000000000000000 0.0000000000000000
45 34.742136171983248 0.0000000000000000 0.0000000000000000
46 29.189350596637517 0.0000000000000000 0.0000000000000000
47 24.037486186421319 0.0000000000000000 0.0000000000000000
48 19.003436601487543 0.0000000000000000 0.0000000000000000
49 14.120045888854746 0.0000000000000000 0.0000000000000000
50 9.7244659609280415 0.0000000000000000 0.0000000000000000
51 5.6634823356386530 0.0000000000000000 0.0000000000000000
52 1.8383620561016323 0.0000000000000000 0.0000000000000000
$EndElementData
$NodeData
1
"Magnetic Field (T)"
1
1.0000000000000000E-008
3
1000
3
51
1 0.0000000000000000 0.0000000000000000 0.0000000000000000
2 0.0000000000000000 0.0000000000000000 0.0000000000000000
3 0.0000000000000000 0.0000000000000000 0.0000000000000000
4 0.0000000000000000 0.0000000000000000 0.0000000000000000
5 0.0000000000000000 0.0000000000000000 0.0000000000000000
6 0.0000000000000000 0.0000000000000000 0.0000000000000000
7 0.0000000000000000 0.0000000000000000 0.0000000000000000
8 0.0000000000000000 0.0000000000000000 0.0000000000000000
9 0.0000000000000000 0.0000000000000000 0.0000000000000000
10 0.0000000000000000 0.0000000000000000 0.0000000000000000
11 0.0000000000000000 0.0000000000000000 0.0000000000000000
12 0.0000000000000000 0.0000000000000000 0.0000000000000000
13 0.0000000000000000 0.0000000000000000 0.0000000000000000
14 0.0000000000000000 0.0000000000000000 0.0000000000000000
15 0.0000000000000000 0.0000000000000000 0.0000000000000000
16 0.0000000000000000 0.0000000000000000 0.0000000000000000
17 0.0000000000000000 0.0000000000000000 0.0000000000000000
18 0.0000000000000000 0.0000000000000000 0.0000000000000000
19 0.0000000000000000 0.0000000000000000 0.0000000000000000
20 0.0000000000000000 0.0000000000000000 0.0000000000000000
21 0.0000000000000000 0.0000000000000000 0.0000000000000000
22 0.0000000000000000 0.0000000000000000 0.0000000000000000
23 0.0000000000000000 0.0000000000000000 0.0000000000000000
24 0.0000000000000000 0.0000000000000000 0.0000000000000000
25 0.0000000000000000 0.0000000000000000 0.0000000000000000
26 0.0000000000000000 0.0000000000000000 0.0000000000000000
27 0.0000000000000000 0.0000000000000000 0.0000000000000000
28 0.0000000000000000 0.0000000000000000 0.0000000000000000
29 0.0000000000000000 0.0000000000000000 0.0000000000000000
30 0.0000000000000000 0.0000000000000000 0.0000000000000000
31 0.0000000000000000 0.0000000000000000 0.0000000000000000
32 0.0000000000000000 0.0000000000000000 0.0000000000000000
33 0.0000000000000000 0.0000000000000000 0.0000000000000000
34 0.0000000000000000 0.0000000000000000 0.0000000000000000
35 0.0000000000000000 0.0000000000000000 0.0000000000000000
36 0.0000000000000000 0.0000000000000000 0.0000000000000000
37 0.0000000000000000 0.0000000000000000 0.0000000000000000
38 0.0000000000000000 0.0000000000000000 0.0000000000000000
39 0.0000000000000000 0.0000000000000000 0.0000000000000000
40 0.0000000000000000 0.0000000000000000 0.0000000000000000
41 0.0000000000000000 0.0000000000000000 0.0000000000000000
42 0.0000000000000000 0.0000000000000000 0.0000000000000000
43 0.0000000000000000 0.0000000000000000 0.0000000000000000
44 0.0000000000000000 0.0000000000000000 0.0000000000000000
45 0.0000000000000000 0.0000000000000000 0.0000000000000000
46 0.0000000000000000 0.0000000000000000 0.0000000000000000
47 0.0000000000000000 0.0000000000000000 0.0000000000000000
48 0.0000000000000000 0.0000000000000000 0.0000000000000000
49 0.0000000000000000 0.0000000000000000 0.0000000000000000
50 0.0000000000000000 0.0000000000000000 0.0000000000000000
51 0.0000000000000000 0.0000000000000000 0.0000000000000000
$EndNodeData

View file

@ -0,0 +1,247 @@
$MeshFormat
2.2 0 8
$EndMeshFormat
$NodeData
1
"Electron density (m^-3)"
1
1.0000000000000000E-008
3
1000
1
51
1 6.926728E+018
2 4.338487E+014
3 1.201033E+018
4 1.505055E+017
5 7.986344E+016
6 4.403063E+016
7 3.228092E+016
8 2.775434E+016
9 2.012556E+016
10 1.435034E+016
11 1.012383E+016
12 7.086783E+015
13 6.226656E+015
14 6.362452E+015
15 6.664819E+015
16 5.946772E+015
17 4.875712E+015
18 4.419041E+015
19 4.092162E+015
20 3.771210E+015
21 3.309077E+015
22 3.071732E+015
23 2.732966E+015
24 2.484798E+015
25 2.313440E+015
26 2.142727E+015
27 1.892093E+015
28 1.776155E+015
29 1.672823E+015
30 1.493186E+015
31 1.407724E+015
32 1.403006E+015
33 1.262097E+015
34 1.182309E+015
35 1.134819E+015
36 1.047928E+015
37 9.865986E+014
38 9.544875E+014
39 8.680888E+014
40 8.440321E+014
41 7.961164E+014
42 7.174005E+014
43 6.589875E+014
44 6.158189E+014
45 6.203645E+014
46 5.604593E+014
47 5.677392E+014
48 5.696204E+014
49 5.003183E+014
50 4.735248E+014
51 4.522101E+014
$EndNodeData
$NodeData
1
"Electron velocity (m s^-1)"
1
1.0000000000000000E-008
3
1000
3
51
1 6.755624E+003 2.474633E+002 2.390054E+002
2 2.054354E+005 -5.163376E+003 -9.846850E+003
3 1.851665E+004 5.899050E+001 2.941097E+002
4 5.874981E+003 -2.544408E+001 3.338329E+002
5 5.745724E+003 -2.992103E+002 -1.562428E+002
6 9.042765E+003 5.755962E+002 -7.509840E+002
7 1.189973E+004 1.251817E+003 -1.278355E+003
8 1.171915E+004 2.483533E+002 -1.798212E+003
9 1.539467E+004 -4.794788E+002 -1.578985E+002
10 1.749366E+004 -2.133315E+003 2.884641E+003
11 3.073569E+004 -1.813877E+003 2.015157E+002
12 4.338021E+004 4.080113E+003 -1.903833E+002
13 4.108799E+004 1.245391E+003 -1.163594E+003
14 3.541737E+004 2.090477E+003 2.419701E+002
15 3.427041E+004 2.013895E+003 -2.234277E+003
16 4.230007E+004 -2.398689E+003 1.371723E+002
17 4.749322E+004 1.597246E+003 -3.628288E+003
18 4.997666E+004 1.160104E+003 -8.401792E+003
19 5.195058E+004 -8.364161E+002 -1.802101E+003
20 5.591312E+004 -3.939060E+003 2.388800E+003
21 6.333098E+004 -6.945454E+002 2.258530E+002
22 6.525072E+004 -3.210740E+003 1.801044E+003
23 7.399242E+004 -1.986068E+003 2.887907E+002
24 7.895156E+004 3.349907E+003 -2.339613E+003
25 8.100325E+004 3.214988E+003 9.916669E+002
26 8.973728E+004 7.594115E+002 -3.305504E+003
27 9.309609E+004 1.417677E+003 -1.436153E+002
28 9.726401E+004 -3.383752E+003 2.393699E+003
29 1.004495E+005 -2.281003E+003 3.209291E+003
30 1.107262E+005 -1.404186E+003 7.470032E+002
31 1.202574E+005 -2.703508E+003 4.930279E+003
32 1.224420E+005 -3.973287E+002 -7.416719E+003
33 1.312544E+005 -5.905453E+003 -9.093921E+002
34 1.347764E+005 -3.491930E+003 -4.010455E+003
35 1.376741E+005 1.357565E+003 2.096097E+003
36 1.419675E+005 1.982650E+003 6.535850E+002
37 1.437177E+005 9.484341E+002 -1.187650E+004
38 1.468643E+005 2.894937E+003 -5.271059E+002
39 1.559764E+005 2.200720E+003 -6.344991E+003
40 1.599147E+005 -2.663692E+003 5.719357E+003
41 1.567097E+005 -7.359025E+000 -4.033384E+003
42 1.690415E+005 -6.956095E+002 3.022908E+003
43 1.814588E+005 -9.615306E+001 6.917210E+003
44 1.863152E+005 -7.990127E+002 1.523384E+004
45 1.886311E+005 4.934329E+003 4.438134E+002
46 1.938832E+005 -3.121504E+002 -3.067774E+003
47 1.900862E+005 -1.160531E+003 -4.620352E+003
48 1.923784E+005 -3.210640E+003 2.151672E+003
49 2.043583E+005 2.368308E+003 1.749902E+002
50 2.054704E+005 4.364283E+003 1.393513E+003
51 2.125792E+005 -4.749345E+003 -2.018864E+003
$EndNodeData
$NodeData
1
"Electron pressure (Pa)"
1
1.0000000000000000E-008
3
1000
1
51
1 2.252567E-001
2 7.546705E-006
3 4.045878E-002
4 5.008147E-003
5 2.551765E-003
6 1.488666E-003
7 1.021884E-003
8 7.937038E-004
9 5.788516E-004
10 4.290930E-004
11 3.263161E-004
12 2.445814E-004
13 2.004129E-004
14 1.799644E-004
15 1.711394E-004
16 1.505493E-004
17 1.238843E-004
18 1.124456E-004
19 9.851639E-005
20 8.889251E-005
21 7.626551E-005
22 6.799265E-005
23 5.942787E-005
24 5.395244E-005
25 4.992406E-005
26 4.474483E-005
27 4.103992E-005
28 3.826699E-005
29 3.404848E-005
30 3.038250E-005
31 2.882086E-005
32 2.844685E-005
33 2.486692E-005
34 2.305064E-005
35 2.204186E-005
36 2.045655E-005
37 1.912248E-005
38 1.817270E-005
39 1.663693E-005
40 1.616886E-005
41 1.470869E-005
42 1.249547E-005
43 1.189200E-005
44 1.112158E-005
45 1.148139E-005
46 1.015354E-005
47 9.973794E-006
48 9.467041E-006
49 8.605269E-006
50 8.068109E-006
51 8.336811E-006
$EndNodeData
$NodeData
1
"Electron temperature (K)"
1
1.0000000000000000E-008
3
1000
1
51
1 2.355409E+003
2 1.259899E+003
3 2.439915E+003
4 2.410136E+003
5 2.314246E+003
6 2.448833E+003
7 2.292833E+003
8 2.071307E+003
9 2.083225E+003
10 2.165738E+003
11 2.334589E+003
12 2.499719E+003
13 2.331244E+003
14 2.048702E+003
15 1.859852E+003
16 1.833641E+003
17 1.840328E+003
18 1.843025E+003
19 1.743703E+003
20 1.707266E+003
21 1.669315E+003
22 1.603229E+003
23 1.574972E+003
24 1.572667E+003
25 1.563034E+003
26 1.512492E+003
27 1.571017E+003
28 1.560488E+003
29 1.474228E+003
30 1.473759E+003
31 1.482881E+003
32 1.468559E+003
33 1.427073E+003
34 1.412111E+003
35 1.406820E+003
36 1.413897E+003
37 1.403850E+003
38 1.379006E+003
39 1.388116E+003
40 1.387513E+003
41 1.338179E+003
42 1.261559E+003
43 1.307058E+003
44 1.308069E+003
45 1.340493E+003
46 1.312170E+003
47 1.272414E+003
48 1.203776E+003
49 1.245762E+003
50 1.234087E+003
51 1.335293E+003
$EndNodeData

7
runs/cylFlow/README.txt Normal file
View file

@ -0,0 +1,7 @@
Neutral flow around a cylinder with elastic collisions.
Case with the same mesh for particle weighting and collisions and another case using the dual mesh mode, with a more refined mesh around the cylinder.
This demostrates the capacity of fpakc to work with independent meshes for particles and collisions.
CPU time is outputed.
Gnuplot scripts for plotting the CPU time are provided.

View file

@ -43,9 +43,8 @@
},
"solver": {
"tau": [5.0e-7],
"finalTime": 1.0e-3,
"finalTime": 5.0e-4,
"pusher": ["Neutral"],
"WeightingScheme": "Volume"
},
"interactions": {
"folderCollisions": "./data/collisions/",

View file

@ -44,7 +44,7 @@
},
"solver": {
"tau": [5.0e-7],
"finalTime": 1.0e-3,
"finalTime": 5.0e-4,
"pusher": ["Neutral"],
"WeightingScheme": "Volume"
},

File diff suppressed because it is too large Load diff

File diff suppressed because it is too large Load diff

File diff suppressed because it is too large Load diff

File diff suppressed because it is too large Load diff