The momentum flux conservation equation includes the momentum transfer of two gases (heating gas and beam electron).

Through derivation the present thesis proves numerical value flux conservation of sub-domains interface boundary can make parallel numerical schemes conservative.

To relax the demand on computer memory and to raise the efficiency when the computational grid number is very large, the flux conservation multizone methods are applied.