The scaling law for the error field (EF) penetration threshold is predicted numerically based on non-linear single-fluid and two-fluid modeling using the TM1 code. The simulated penetration threshold of radial magnetic field br at the plasma edge is scaled to the electron density ne, temperature Te, viscous time , toroidal field Bt and the natural frequency ω in the form of by scanning these parameters separately. Here, αn, αT, , αB and are the scaling coefficients on ne, Te, , Bt and ω, respectively. Single-fluid modeling shows that the 3/2 EF threshold scales as , which is similar with the analytical scaling law in both the Rutherford and visco-resistive regimes. However, two-fluid modeling shows that the scaling law differs significantly in particular regarding the dependence on plasma rotation. In detail, the scaling coefficient αn on density decreases from 0.67 to 0.56 and αT on temperature decreases from 0.67 to 0.32, while on viscous time is around -0.45 and αB on toroidal field decreases slightly from -1.15 to -1, when the ratio between plasma rotation frequency ωE and diamagnetic drift frequency ω*e varies from 0 to 10. Scans of the plasma rotation reveals that the penetration threshold linearly depends on the perpendicular electron flow frequency (or natural frequency) , and there is a minimum in the required field amplitude when . In addition, the enduring mystery of non-zero penetration threshold at zero plasma natural frequency in EF experiments is resolved by two-fluid simulations. We find that the very small island and smooth bifurcation in EF penetration near zero frequency is hard to detect in the experiment, leading to a finite penetration threshold within the capability of the experimental measurements.
Nonlinear error-field penetration in low density ohmically heated tokamak plasmas