Smaller grid size and BC to ensure reaching a quasi-steady state after 80 ns in the domain

plasmaExpansion.f90

@@ -100,7 +100,7 @@ program plasmaExpansion
   ! Set domain boundaries (non-dimensional units)
   r0 =  10.0e-6_dp / L_ref
   rf =   2.0e-3_dp / L_ref
-  dr = 5.0e2_dp
+  dr =   1.0e-6_dp / L_ref
   nr = nint((rf - r0) / dr) + 1
   dr =      (rf - r0) / float(nr-1)
   allocate(r(1:nr))
@@ -280,7 +280,7 @@ program plasmaExpansion
     n_e = Zave * n_i
 
     ! Solve Poission (maximum number of iterations, break if convergence is reached before)
-    do k = 1, 1000
+    do k = 1, 2000
       ! Store previous value
       phi_old = phi
 
@@ -319,7 +319,7 @@ program plasmaExpansion
 
       ! Check if the solution has converged
       phiConv = maxval(abs(Res),1)
-      if (phiConv < 1.0e-4_dp) then
+      if (phiConv < 1.0e-6_dp) then
         exit
 
       end if