jamesr wrote:
The link to the PDF is fixed now, sorry for the hang up.
I realise the paper is still work-in-progress, and the model as a whole seems a pretty solid 1D scheme. But one thing struck me, you seem to be assuming the ions are singly charged, so the quasi-neutrality condition is simply n = n_e = n_i, and subsequently the current is J_z=qn(v_e – v_i). Would it not be easier to set the ion species charge as q=Ze, so
n=n_e = Zn_i and J_z=en(v_e – Zv_i). The step up from a deuterium plasma to one with a higher Z_effective is then carried through in all the equations.
You can say these simulations are for deuterium so Z=1, but having it explicitly in the equations always seems clearer to me.
You make a valid point, but the model changes to do p-B11 will actually have to be a lot more complicated. We are now contemplating a model with many different species (still in context of deuterium) which clearly would be different with decaborane.