From 95d5a1200d389375fa1341cefd8df9f9e500700a Mon Sep 17 00:00:00 2001 From: Jorge Gonzalez Date: Thu, 11 Mar 2021 12:54:51 +0100 Subject: [PATCH 1/2] Included a new cases of 1D Cathode emission based on plasma background previously computed. --- runs/1D_Cathode/Argon+_Background.dat | 51 +++++++++++++++++++ runs/1D_Cathode/Electron_Background.dat | 51 +++++++++++++++++++ .../{inputRad.json => inputRadEmission.json} | 25 +++++---- runs/1D_Cathode/mesh.geo | 2 +- 4 files changed, 117 insertions(+), 12 deletions(-) create mode 100644 runs/1D_Cathode/Argon+_Background.dat create mode 100644 runs/1D_Cathode/Electron_Background.dat rename runs/1D_Cathode/{inputRad.json => inputRadEmission.json} (72%) diff --git a/runs/1D_Cathode/Argon+_Background.dat b/runs/1D_Cathode/Argon+_Background.dat new file mode 100644 index 0000000..b09138f --- /dev/null +++ b/runs/1D_Cathode/Argon+_Background.dat @@ -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 diff --git a/runs/1D_Cathode/Electron_Background.dat b/runs/1D_Cathode/Electron_Background.dat new file mode 100644 index 0000000..e03cc49 --- /dev/null +++ b/runs/1D_Cathode/Electron_Background.dat @@ -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 diff --git a/runs/1D_Cathode/inputRad.json b/runs/1D_Cathode/inputRadEmission.json similarity index 72% rename from runs/1D_Cathode/inputRad.json rename to runs/1D_Cathode/inputRadEmission.json index 3e87739..ba3ce41 100644 --- a/runs/1D_Cathode/inputRad.json +++ b/runs/1D_Cathode/inputRadEmission.json @@ -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":{ diff --git a/runs/1D_Cathode/mesh.geo b/runs/1D_Cathode/mesh.geo index 657bb89..3f96cbb 100644 --- a/runs/1D_Cathode/mesh.geo +++ b/runs/1D_Cathode/mesh.geo @@ -1,4 +1,4 @@ -Lcell = 0.00004; +Lcell = 4e-5; x0 = 0.001; xf = x0 + 50.0*Lcell; From ff0c09242dfff42ed2536577c12f37d3bcb075b9 Mon Sep 17 00:00:00 2001 From: Jorge Gonzalez Date: Thu, 11 Mar 2021 12:56:15 +0100 Subject: [PATCH 2/2] Fixed an issue with normal vector not being norm unit, which was causing reflected particle velocities to increase. --- src/modules/mesh/1DCart/moduleMesh1DCart.f90 | 1 + src/modules/mesh/1DRad/moduleMesh1DRad.f90 | 1 + src/modules/mesh/2DCart/moduleMesh2DCart.f90 | 2 ++ src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 | 1 + src/modules/mesh/3DCart/moduleMesh3DCart.f90 | 1 + 5 files changed, 6 insertions(+) diff --git a/src/modules/mesh/1DCart/moduleMesh1DCart.f90 b/src/modules/mesh/1DCart/moduleMesh1DCart.f90 index 647d974..d45c7a0 100644 --- a/src/modules/mesh/1DCart/moduleMesh1DCart.f90 +++ b/src/modules/mesh/1DCart/moduleMesh1DCart.f90 @@ -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) diff --git a/src/modules/mesh/1DRad/moduleMesh1DRad.f90 b/src/modules/mesh/1DRad/moduleMesh1DRad.f90 index 3f67f15..1d2bdba 100644 --- a/src/modules/mesh/1DRad/moduleMesh1DRad.f90 +++ b/src/modules/mesh/1DRad/moduleMesh1DRad.f90 @@ -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) diff --git a/src/modules/mesh/2DCart/moduleMesh2DCart.f90 b/src/modules/mesh/2DCart/moduleMesh2DCart.f90 index cc2369a..ac22b28 100644 --- a/src/modules/mesh/2DCart/moduleMesh2DCart.f90 +++ b/src/modules/mesh/2DCart/moduleMesh2DCart.f90 @@ -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)) diff --git a/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 b/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 index 9303c1b..8ef5c6d 100644 --- a/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 +++ b/src/modules/mesh/2DCyl/moduleMesh2DCyl.f90 @@ -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)) diff --git a/src/modules/mesh/3DCart/moduleMesh3DCart.f90 b/src/modules/mesh/3DCart/moduleMesh3DCart.f90 index d01e820..1c713c9 100644 --- a/src/modules/mesh/3DCart/moduleMesh3DCart.f90 +++ b/src/modules/mesh/3DCart/moduleMesh3DCart.f90 @@ -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)