We use 2D (axisymmetric) and 3D hydrodynamic simulations to study Bondi-Hoyle-Lyttleton (BHL) accretion with and without transverse upstream gradients. We mainly focus on the regime of high (upstream) Mach number, weak upstream gradients and small accretor size, which is relevant to neutron star (NS) accretion in wind-fed Supergiant X-ray binaries (SgXBs). We present a systematic exploration of the flow in this regime. When there are no upstream gradients, the flow is always stable regardless of accretor size or Mach number. For finite upstream gradients, there are three main types of behavior: stable flow (small upstream gradient), turbulent unstable flow without a disk (intermediate upstream gradient), and turbulent flow with a disk-like structure (relatively large upstream gradient). When the accretion flow is turbulent, the accretion rate decreases non-convergently as the accretor size decreases. The flow is more prone to instability and the disk is less likely to form than previously expected; the parameters of most observed SgXBs place them in the regime of a turbulent, diskless accretion flow. Among the SgXBs with relatively well-determined parameters, we find OAO 1657-415 to be the only one that is likely to host a persistent disk (or disk-like structure); this finding is consistent with observations. c 0000 The Authors arXiv:1907.06108v1 [astro-ph.HE] 13 Jul 2019 R , v Φ ) is evaluated in the rotating frame in which the binary is fixed. D = a b − Rc is the separation between the NS and the surface of the companion (assuming circular orbit). ρ,v are the transverse gradients of density and velocity for single-star wind profile at the NS, and physically correspond to the fractional change of ρ, v per Ra (see text for detailed definition). Ra/R H characterizes the importance of the companion's gravity (with R H being the Hill sphere radius), and Ω b ta characterizes the importance of Coriolis force. The last five parameters show how significantly the system differs from axisymmetric BHL accretion.