Merge branch 'issue/reflection' into 'development'

Normal of edges unitary

See merge request JorgeGonz/fpakc!8
This commit is contained in:
Jorge Gonzalez 2021-03-11 11:58:29 +00:00
commit 61885b4a37
9 changed files with 123 additions and 12 deletions

View file

@ -0,0 +1,51 @@
#Element Density(m^-3) Velocity (m/2) Temperature (K)
1 1683217500000000 -19386.6 0 0 500
2 1577004000000000 -19130.6 0 0 500
3 1483725500000000 -18673.4 0 0 500
4 1495368000000000 -18223.2 0 0 500
5 1601330000000000 -17804.8 0 0 500
6 1630851500000000 -17360.2 0 0 500
7 1569924000000000 -16986.4 0 0 500
8 1555674000000000 -16567.4 0 0 500
9 1501328000000000 -16181.2 0 0 500
10 1472318000000000 -15785.9 0 0 500
11 1458490000000000 -15374.6 0 0 500
12 1457809500000000 -15012.6 0 0 500
13 1474687000000000 -14477.2 0 0 500
14 1473073500000000 -14242.7 0 0 500
15 1511770500000000 -13848.8 0 0 500
16 1564946500000000 -13502.7 0 0 500
17 1533674000000000 -13157.7 0 0 500
18 1473108000000000 -12764.5 0 0 500
19 1496640000000000 -12488.3 0 0 500
20 1514440500000000 -12036.1 0 0 500
21 1494114500000000 -11665.3 0 0 500
22 1493046500000000 -11348.2 0 0 500
23 1500858500000000 -10977.3 0 0 500
24 1494480500000000 -10624.9 0 0 500
25 1542839500000000 -11111.7 0 0 500
26 1596599500000000 -9940.4 0 0 500
27 1650272000000000 -9566.4 0 0 500
28 1690817500000000 -9223.5 0 0 500
29 1724662000000000 -8917.5 0 0 500
30 1777677500000000 -8544.2 0 0 500
31 1809272500000000 -8129.9 0 0 500
32 1814150000000000 -7827.7 0 0 500
33 1838401500000000 -7458.4 0 0 500
34 1977935000000000 -7113.4 0 0 500
35 2093814500000000 -6707.0 0 0 500
36 2107688500000000 -6373.7 0 0 500
37 2244944000000000 -6024.8 0 0 500
38 2390213000000000 -5625.1 0 0 500
39 2484013500000000 -5223.9 0 0 500
40 2600636000000000 -4862.8 0 0 500
41 2751150500000000 -4640.7 0 0 500
42 2972491500000000 -4116.6 0 0 500
43 3282084500000000 -3661.9 0 0 500
44 3682392500000000 -3241.6 0 0 500
45 4103775500000000 -2772.5 0 0 500
46 4913489500000000 -2383.0 0 0 500
47 6405310500000000 -1911.6 0 0 500
48 9003078999999999 -1483.9 0 0 500
49 1.3388345E+016 -993.4 0 0 500
50 2.266645E+016 -592.5 0 0 500

View file

@ -0,0 +1,51 @@
#Element Density(m^-3) Velocity (m/2) Temperature (K)
1 0 0 0 0 11604
2 0 0 0 0 11604
3 0 0 0 0 11604
4 0 0 0 0 11604
5 0 0 0 0 11604
6 0 0 0 0 11604
7 0 0 0 0 11604
8 0 0 0 0 11604
9 0 0 0 0 11604
10 0 0 0 0 11604
11 0 0 0 0 11604
12 0 0 0 0 11604
13 0 0 0 0 11604
14 0 0 0 0 11604
15 0 0 0 0 11604
16 0 0 0 0 11604
17 0 0 0 0 11604
18 0 0 0 0 11604
19 0 0 0 0 11604
20 0 0 0 0 11604
21 0 0 0 0 11604
22 0 0 0 0 11604
23 0 0 0 0 11604
24 0 0 0 0 11604
25 0 0 0 0 11604
26 0 0 0 0 11604
27 0 0 0 0 11604
28 0 0 0 0 11604
29 0 0 0 0 11604
30 0 0 0 0 11604
31 0 0 0 0 11604
32 0 0 0 0 11604
33 0 0 0 0 11604
34 0 0 0 0 11604
35 0 0 0 0 11604
36 0 0 0 0 11604
37 0 0 0 0 11604
38 0 0 0 0 11604
39 0 0 0 0 11604
40 0 0 0 0 11604
41 0 0 0 0 11604
42 0 0 0 0 11604
43 529256500000 32881.47 0 0 11604
44 5035006500000 40868.50 0 0 11604
45 37296800000000 5264.45 0 0 11604
46 226267350000000 -12113.65 0 0 11604
47 928316300000000 -15834.22 0 0 11604
48 3119752000000000 -12139.49 0 0 11604
49 7748062000000000 -8297.72 0 0 11604
50 1.348286E+016 43.29 0 0 11604

View file

@ -6,12 +6,12 @@
"cpuTime": false,
"numColl": false,
"EMField": true,
"folder": "Radial"
"folder": "Radial_Emision"
},
"reference": {
"density": 1.0e16,
"mass": 9.109e-31,
"temperature": 2500.0
"temperature": 11604.0
},
"geometry": {
"type": "1DRad",
@ -33,22 +33,25 @@
]}
],
"boundaryEM": [
{"name": "Cathode", "type": "dirichlet", "potential": -10.0, "physicalSurface": 1},
{"name": "Infinite", "type": "dirichlet", "potential": 0.0, "physicalSurface": 2}
{"name": "Cathode", "type": "dirichlet", "potential": 0.0, "physicalSurface": 1}
],
"inject": [
{"name": "Cathode Electron", "species": "Electron", "flow": 1.0e-3, "units": "A", "v": 27500.0, "T": [2500.0, 2500.0, 2500.0],
"velDist": ["Maxwellian", "Maxwellian", "Maxwellian"], "n": [ 1, 0, 0], "physicalSurface": 1},
{"name": "Plasma Inf Ar+", "species": "Argon+", "flow": 1.0e-6, "units": "A", "v": 500.0, "T": [500.0, 500.0, 500.0],
{"name": "Plasma Inf Ar+", "species": "Argon+", "flow": 1.00e-6, "units": "A", "v": 300.0, "T": [ 500.0, 500.0, 500.0],
"velDist": ["Maxwellian", "Maxwellian", "Maxwellian"], "n": [-1, 0, 0], "physicalSurface": 2},
{"name": "Plasma Inf e", "species": "Electron", "flow": 1.0e-6, "units": "A", "v": 500.0, "T": [500.0, 500.0, 500.0],
"velDist": ["Maxwellian", "Maxwellian", "Maxwellian"], "n": [-1, 0, 0], "physicalSurface": 2}
{"name": "Plasma Inf e", "species": "Electron", "flow": 2.64e-4, "units": "A", "v": 0.0, "T": [11604.0, 11604.0, 11604.0],
"velDist": ["Maxwellian", "Maxwellian", "Maxwellian"], "n": [-1, 0, 0], "physicalSurface": 2},
{"name": "Plasma Cat e", "species": "Electron", "flow": 2.64e-2, "units": "A", "v": 180000.0, "T": [ 2300.0, 2300.0, 2300.0],
"velDist": ["Maxwellian", "Maxwellian", "Maxwellian"], "n": [ 1, 0, 0], "physicalSurface": 1}
],
"case": {
"tau": [1.0e-11, 1.0e-11],
"time": 1.0e-5,
"time": 4.0e-6,
"pusher": ["1DRadCharged", "1DRadCharged"],
"EMSolver": "Electrostatic"
"EMSolver": "Electrostatic",
"initial": [
{"speciesName": "Argon+", "initialState": "Argon+_Background.dat"},
{"speciesName": "Electron", "initialState": "Electron_Background.dat"}
]
},
"parallel": {
"OpenMP":{

View file

@ -1,4 +1,4 @@
Lcell = 0.00004;
Lcell = 4e-5;
x0 = 0.001;
xf = x0 + 50.0*Lcell;

View file

@ -145,6 +145,7 @@ MODULE moduleMesh1DCart
self%x = r1(1)
self%normal = (/ 1.D0, 0.D0, 0.D0 /)
self%normal = self%normal/NORM2(self%normal)
!Boundary index
self%boundary => boundary(bt)

View file

@ -146,6 +146,7 @@ MODULE moduleMesh1DRad
self%r = r1(1)
self%normal = (/ 1.D0, 0.D0, 0.D0 /)
self%normal = self%normal/NORM2(self%normal)
!Boundary index
self%boundary => boundary(bt)

View file

@ -193,6 +193,8 @@ MODULE moduleMesh2DCart
self%normal = (/ self%y(2)-self%y(1), &
self%x(2)-self%x(1), &
0.D0 /)
self%normal = self%normal/NORM2(self%normal)
!Boundary index
self%boundary => boundary(bt)
ALLOCATE(self%fboundary(1:nSpecies))

View file

@ -194,6 +194,7 @@ MODULE moduleMesh2DCyl
self%normal = (/ self%r(2)-self%r(1), &
self%z(2)-self%z(1), &
0.D0 /)
self%normal = self%normal/NORM2(self%normal)
!Boundary index
self%boundary => boundary(bt)
ALLOCATE(self%fboundary(1:nSpecies))

View file

@ -159,6 +159,7 @@ MODULE moduleMesh3DCart
self%normal = (/ (self%y(2)-self%y(1))*(self%z(3)-self%z(1)) - (self%z(2)-self%z(1))*(self%y(3)-self%y(1)), &
(self%x(2)-self%x(1))*(self%z(3)-self%z(1)) - (self%z(2)-self%z(1))*(self%x(3)-self%x(1)), &
(self%x(2)-self%x(1))*(self%y(3)-self%y(1)) - (self%z(2)-self%z(1))*(self%y(3)-self%y(1)) /)
self%normal = self%normal/NORM2(self%normal)
!Boundary index
self%boundary => boundary(bt)