readme.powbal

一维等离子体静电粒子模拟程序

                                                      October 2000

         INSTRUCTIONS FOR USING POWER BALANCE MODULE 
	             /Emi Kawamura/
  
  In order to use the power balance module in XPDC1, you must edit the 
  relevant makefile (makefile, makefile.par.p4, makefile.par.sh) so that 
  the line "DEFINES = -DPOW" is uncommented before compiling the code. 
  
  If you do this, then every POW_DT timesteps, XPDC1 will dump the power 
  and particle balance information into an ascii file with the suffix 
  ".dmp.pow". POW_DT is currently set to 10000 timesteps.  However, you 
  can change this by editing the line "#define POW_DT 10000" in the file 
  "pdc1.c".

  An example of the output from the code running on the reduced problem 
  (E_z = 100 V, Iin = 15.2 mA, height = 20 cm, R=1cm, p=2.83 mTorr)
  looks as follows:
  
   time = 1.680000e-05
   species=0, particle loss to wall =2.2049e+10
   species=1, particle loss to wall =2.19716e+10
   number of  ionization events =2.1928e+10
   species=0, power loss to wall = 0.0670558 W
   species=1, power loss to wall = 0.118345 W
   power loss to inelastic colls = 0.10719 W
   power loss to elastic colls = 6.18593e-05 W
   power loss to charge exchange colls = 0.00936686 W
   total power loss to colls = 0.116619 W
   total power loss to wall = 0.185401 W
   total power loss = 0.30202 W
   power Input = 0.300904 W
 
  Here, species 0 is the electron, and species 1 are the Ar+ ions.
  Note that the simulation has reached equilibrium since the
  particle loss to the wall for each species in POW_DT timesteps
  is equal to the number of ionization events in POW_DT timesteps.
  Also the total power loss equals the power input (Ez*Iz*height).
  
  The power balance module assumes a positive column discharge. So, you
  may want to suppress it by commenting it out of the relevant makefile 
  if you are simulating other types of discharges.