We investigate accretion onto a central star, with the size, rotation rate, and magnetic dipole of a young stellar object, to study the flow pattern (velocity and density) of the fluid within and outside of the disc. We perform resistive magnetohydrodynamic (MHD) simulations of thin discs, varying the parameters such as the stellar rotation rate and (anomalous) coefficients of viscosity and resistivity in the disc. To provide a benchmark for the results and to compare them with known analytic results, we also perform purely hydrodynamic (HD) simulations for the same problem. Although obtained for different situations with differing inner boundary condition, the disc structure in the HD simulations closely follows the analytic solution of Kluzniak and Kita - in particular, a region of 'mid-plane' backflow exists in the right range of radii, depending on the viscosity parameter. In the MHD solutions, whenever the magnetic Prandtl number does not exceed a certain critical value, the mid-plane backflow exists throughout the accretion disc, extending all the way down to the foot point of the accretion funnel flow where the disc transitions to a magnetic funnel flow. For values of the magnetic Prandtl number close to the critical value the backflow and the inner disc undergo a quasi-periodic radial oscillation, otherwise the backflow is steady, as is the disc solution. From our results, supplemented by our reading of the literature, we suggest that mid-plane backflow is a real, physical, and not only numerical feature of at least some accretion discs.