Numerical modeling of surf zone dynamics under weakly plunging breakers with SPH method