feat: introduce normalized electron collision damping in J-equation

Author: alfkoehn

include/focal-struct.h

@@ -13,7 +13,8 @@ typedef struct gridConfiguration {      /*Saves the main physical parameters of
     double
         period,
         dx,dt,
-        ne_0, B0_value;
+        ne_0, B0_value,
+        nu0;
 } gridConfiguration;
 
 typedef struct saveData{                /*Variables related to simulation data saving*/

include/macros-grid.h

@@ -18,6 +18,7 @@
 #define B0_profileG(gridCfg)    gridCfg->B0_profile
 #define B0_valueG(gridCfg)      gridCfg->B0_value
 #define boundaryG(gridCfg)      gridCfg->sel_boundary
+#define nu0_G(gridCfg)          gridCfg->nu0
 
 #define NX                      NxG(gridCfg)            
 #define NY                      NyG(gridCfg)          
@@ -33,6 +34,7 @@
 #define B0_profile              B0_profileG(gridCfg)
 #define B0_value                B0_valueG(gridCfg)
 #define BOUNDARY                boundaryG(gridCfg)
+#define NU0                     nu0_G(gridCfg)
 
 /*Macros for save data*/
 #define projectPathSt(saveDCfg)         saveDCfg->projectPath

src/focal.c

@@ -67,26 +67,29 @@ int advance_J( gridConfiguration *gridCfg,
         for (jj=2 ; jj<NY-2 ; jj+=2) {
             for (kk=2 ; kk<NZ-2 ; kk+=2) {
                 // Jx: odd-even-even
-                J_B0[ii+1][jj  ][kk  ]    += + DT*(
-                        pow(2*M_PI,2) * n_e[(ii/2)][(jj/2)][(kk/2)] * EB_WAVE[ii+1][jj  ][kk  ]
-                        - 2*M_PI * ( J_B0[ii  ][jj+1][kk  ]*J_B0[ii+1][jj+1][kk  ]        // +Jy*B0z
-                                    -J_B0[ii  ][jj  ][kk+1]*J_B0[ii+1][jj  ][kk+1]        // -Jz*B0y
-                              )
-                        );
+                //J_B0[ii+1][jj  ][kk  ]    += + DT*(
+                J_B0[ii+1][jj  ][kk  ]  = (1. - NU0)*J_B0[ii+1][jj  ][kk  ] 
+                    + DT*( pow(2*M_PI,2) * n_e[(ii/2)][(jj/2)][(kk/2)] * EB_WAVE[ii+1][jj  ][kk  ]
+                          -2*M_PI * ( J_B0[ii  ][jj+1][kk  ]*J_B0[ii+1][jj+1][kk  ]        // +Jy*B0z
+                                      -J_B0[ii  ][jj  ][kk+1]*J_B0[ii+1][jj  ][kk+1]        // -Jz*B0y
+                                     )
+                         );
                 // Jy: even-odd-even
-                J_B0[ii  ][jj+1][kk  ]    += + DT*(
-                        pow(2*M_PI,2) * n_e[(ii/2)][(jj/2)][(kk/2)] * EB_WAVE[ii  ][jj+1][kk  ]
-                        -2*M_PI * (-J_B0[ii+1][jj  ][kk  ]*J_B0[ii+1][jj+1][kk  ]         // -Jx*B0z
-                                   +J_B0[ii  ][jj  ][kk+1]*J_B0[ii  ][jj+1][kk+1]         // +Jz*B0x
-                              )
-                        );
+                //J_B0[ii  ][jj+1][kk  ]    += + DT*(
+                J_B0[ii  ][jj+1][kk  ]  = (1. - NU0)*J_B0[ii  ][jj+1][kk  ] 
+                    + DT*( pow(2*M_PI,2) * n_e[(ii/2)][(jj/2)][(kk/2)] * EB_WAVE[ii  ][jj+1][kk  ]
+                          -2*M_PI * (-J_B0[ii+1][jj  ][kk  ]*J_B0[ii+1][jj+1][kk  ]         // -Jx*B0z
+                                     +J_B0[ii  ][jj  ][kk+1]*J_B0[ii  ][jj+1][kk+1]         // +Jz*B0x
+                                    )
+                         );
                 // Jz: even-even-odd
-                J_B0[ii  ][jj  ][kk+1]    += + DT*(
-                        pow(2*M_PI,2) * n_e[(ii/2)][(jj/2)][(kk/2)] * EB_WAVE[ii  ][jj  ][kk+1]
-                        -2*M_PI * ( J_B0[ii+1][jj  ][kk  ]*J_B0[ii+1][jj  ][kk+1]        // +Jx*B0y
-                                   -J_B0[ii  ][jj+1][kk  ]*J_B0[ii  ][jj+1][kk+1]        // -Jy*B0x
-                              )
-                        );
+                //J_B0[ii  ][jj  ][kk+1]    += + DT*(
+                J_B0[ii  ][jj  ][kk+1]  = (1. - NU0)*J_B0[ii  ][jj  ][kk+1]
+                    + DT*( pow(2*M_PI,2) * n_e[(ii/2)][(jj/2)][(kk/2)] * EB_WAVE[ii  ][jj  ][kk+1]
+                          -2*M_PI * ( J_B0[ii+1][jj  ][kk  ]*J_B0[ii+1][jj  ][kk+1]        // +Jx*B0y
+                                     -J_B0[ii  ][jj+1][kk  ]*J_B0[ii  ][jj+1][kk+1]        // -Jy*B0x
+                                    )
+                         );
             }
         }
     }

src/init_module.c

@@ -303,6 +303,7 @@ void print_systemConfiguration(gridConfiguration *gridCfg, beamAntennaConfigurat
     printf( "power detector position: %d\n", pwr_dect );
     printf( "Boundary condition set to '%d'\n", BOUNDARY );
     printf( "Courant number = %.2f. \n", DT/DX);
+    printf( "nu_0 (normalized electron collision frequency) = %f\n", NU0 );
 
 }//}}}