The scaling law for the m/n = 3/2 error field (EF) penetration threshold is predicted numerically based on nonlinear single-fluid and two-fluid modeling using the TM1 code. The simulated penetration threshold of radial magnetic field b r at the plasma edge is scaled to the electron density n e , temperature T e , viscous time τ µ , toroidal field B t and the natural frequency ω in the form of b r /B t ∝ n αn e T α T e τ αµ µ B α B t ω αω by scanning these parameters separately. Here, α n , α T , α µ , α B and α ω are the scaling coefficients on n e , T e , τ µ , B t and ω, respectively. Singlefluid modeling shows that the 3/2 EF threshold scales as b r /B t ∝ n 0.56 e T 0.6 e τ −0.59 µ B −1.15 t ω, 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 |ω E /ω * e | 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 rotation frequency (or natural frequency) ω ⊥e = ω E + ω * e , and there is a minimum in the required field amplitude when ω ⊥e ∼ 0. 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.