Reactivities were inaccurate. In Figure 1 below, more accurate reactivities are taken into account.

It was supposed there would be permanent losses by charge exchanges on neutrals. This behavior exists at the beginning of the operation but disappears progressively as the fusion reactor wall is not in contact with the gas, so charge exchanges will not be considered. Therefore, the variables $\gamma cep=\gamma ceT=\gamma ceHe3=\gamma ceHe4=\gamma ceD=0$ are forced to 0.

For a fusion radius of 4.5 m, the results are the following:

These low mechanical gains Q greatly limit the advantage of this type of D-D reactor.

The modified executable program with its Delphi 6 source can be downloaded from this direct link: http://f6cte.free.fr/D_D_reactor_model_V_1_1.zip or on the Zenodo WEB open repository by searching with the title of this article.