Message by Richard Marchand:
I am interested in finding out how many particles were used in a simulation (for all species treated kinetically). An average number for each species would be fine. A plot of N_i as a function of time would be interesting as well. How can one obtain this information?
P.S.: I don’t know how this happened, but I posted this question a few minutes ago, and it showed as being from Julien instead of me.
Message by Jean-Francois Roussel:
This is printed in the log (on the screen and in file SpisNum.log in yor spistmp directory).
Typical form is:
Monitoring numerics of simulation at time 0.75082004
PICVolDistrib ions1 : 84862 particles, last numerical time step used = 1.1759414E-6s (max allowed = 1.7365647E-6s), numerical speed up = 1.0
PICVolDistrib photoElec: 17715 particles, last numerical time step used = 3.0E-5s (max allowed = 5.0E19s), numerical speed up = 1.0
…
I don’t think you can get it as number_versus_time plot (you can still take a look at the “2D plots”)
Message by Richard Marchand:
Thanks for the reply.
I suppose one could get the number of ions as a function of time by writing a little program that would parse the .log file and extract the numbers at the various timesteps.
I didn’t find the numbers of electrons in the log file. I guess those are roughly equal to the numbers of ions.
Message by Richard Marchand:
Thanks for the reply.
I suppose one could get the number of ions as a function of time by writing a little program that would parse the .log file and extract the numbers at the various timesteps.
I didn’t find the numbers of electrons in the log file. I guess those are roughly equal to the numbers of ions.
Message by Jean-Francois Roussel:
Parsing the log file is rathing indirect but possible…
If you don’t get the number of electron, you probably don’t model them through Monte Carlo (PIC) but analytically (Boltzmann).
Message by Richard Marchand:
For the electron distribution, I have “PICVolDistrib”. Shouldn’t that produce PIC electrons?
Message by Jean-Francois Roussel:
It should.
Could you please post (part of) the screen log, including an section where you have somthing about “Monitoring numerics”, which would look like:
Monitoring numerics of simulation at time 0.011963652
PICVolDistrib ions1 : 252542 particles, last numerical time step used = 3.1921337E-7s (max allowed = 1.7365647E-6s), numerical speed up = 1.0
PICVolDistrib photoElec: 2857 particles, last numerical time step used = 3.0E-5s (max allowed = 5.0E19s), numerical speed up = 1.0
Message by Richard Marchand:
Got it.
I’m afraid I had missed it when I first read the SpisNum.log file.
Thanks. Richard.