The model you described is known as the transverse field quantum Ising model. Yes, you are correct that for one phase ($B\gg J$), we have a paramagnetic phase, and for ($B\ll J$), we have a ferromagnetic phase. The quantum phase transition boundary is $B=J$, which can be shown as from Karmer's Duality relation.
What happens when $B\sim J$ or when we reach near criticality. Pretty much everything that happens for our ole classical ising model, like correlation length diverges. The traditional 2D ising model and Transverse field quantum Ising model belong to the same universality class.
What else happens? "Bulk gap" closes here, which has an interesting aspect related to the topological superconductor.
But let's come back to point "universality class." When I say two systems belong to the same universality class, it means that all critical exponents for two systems are identical. In classical Ising model, the magnetisation or the order parameter below criticality goes as
$$M_z\sim (T-T_c)^{-\beta}$$
The same thing happens here.
In classical Ising model, the phase transtion is driven by thermal fluctuation. But, here, the transitions are driven by quantum fluctuations. So treatment is a different form of standard statistical mechanics treatment.
An important feature of this setup is that, in a quantum sense, the
spin projection along the x axis and the spin
projection along the z axis are not commuting
observable quantities. That is, they cannot both be observed
simultaneously. This means classical statistical mechanics cannot
describe this model, and a quantum treatment is needed.
You can find treatment in any standard Quantum Phase Transitions textbooks