A three-dimensional fluid model is developed to investigate the radio-frequency inductively coupled H2 plasma in a reactor with a rectangular expansion chamber and a cylindrical driver chamber, for neutral beam injection system in CFETR. In this model, the electron effective collision frequency and the ion mobility at high E-fields are employed, for accurate simulation of discharges at low pressures (0.3 Pa–2 Pa) and high powers (40 kW–100 kW). The results indicate that when the high E-field ion mobility is taken into account, the electron density is about four times higher than the value in the low E-field case. In addition, the influences of the magnetic field, pressure and power on the electron density and electron temperature are demonstrated. It is found that the electron density and electron temperature in the xz-plane along permanent magnet side become much more asymmetric when magnetic field enhances. However, the plasma parameters in the yz-plane without permanent magnet side are symmetric no matter the magnetic field is applied or not. Besides, the maximum of the electron density first increases and then decreases with magnetic field, while the electron temperature at the bottom of the expansion region first decreases and then almost keeps constant. As the pressure increases from 0.3 Pa to 2 Pa, the electron density becomes higher, with the maximum moving upwards to the driver region, and the symmetry of the electron temperature in the xz-plane becomes much better. As power increases, the electron density rises, whereas the spatial distribution is similar. It can be summarized that the magnetic field and gas pressure have great influence on the symmetry of the plasma parameters, while the power only has little effect.