We solve numerically for the first time the two-fluid Hall–Vinen–Bekarevich–Khalatnikov (HVBK) equations for an He-II-like superfluid contained in a differentially rotating spherical shell, generalizing previous simulations of viscous spherical Couette flow (SCF) and superfluid Taylor–Couette flow. The simulations are conducted for Reynolds numbers in the range 1 × 102≤Re≤3 × 104, rotational shear 0.1≤ΔΩ/Ω≤0.3, and dimensionless gap widths 0.2≤δ≤0.5. The system tends towards a stationary but unsteady state, where the torque oscillates persistently, with amplitude and period determined by δ and ΔΩ/Ω. In axisymmetric superfluid SCF, the number of meridional circulation cells multiplies as Re increases, and their shapes become more complex, especially in the superfluid component, with multiple secondary cells arising for Re > 103. The torque exerted by the normal component is approximately three times greater in a superfluid with anisotropic Hall–Vinen (HV) mutual friction than in a classical viscous fluid or a superfluid with isotropic Gorter–Mellink (GM) mutual friction. HV mutual friction also tends to ‘pinch’ meridional circulation cells more than GM mutual friction. The boundary condition on the superfluid component, whether no slip or perfect slip, does not affect the large-scale structure of the flow appreciably, but it does alter the cores of the circulation cells, especially at lower Re. As Re increases, and after initial transients die away, the mutual friction force dominates the vortex tension, and the streamlines of the superfluid and normal fluid components increasingly resemble each other. In non-axisymmetric superfluid SCF, three-dimensional vortex structures are classified according to topological invariants. For misaligned spheres, the flow is focal throughout most of its volume, except for thread-like zones where it is strain-dominated near the equator (inviscid component) and poles (viscous component). A wedge-shaped isosurface of vorticity rotates around the equator at roughly the rotation period. For a freely precessing outer sphere, the flow is equally strain- and vorticity-dominated throughout its volume. Unstable focus/contracting points are slightly more common than stable node/saddle/saddle points in the viscous component, but not in the inviscid component. Isosurfaces of positive and negative vorticity form interlocking poloidal ribbons (viscous component) or toroidal tongues (inviscid component) which attach and detach at roughly the rotation period.