3D numerical modeling of YSO accretion shocks