Entropy generation for a paramagnetic fluid in a square enclosure with thermomagnetic convection is numerically investigated under the influence of a magnetic quadrupole field. The magnetic field is calculated using the scalar magnetic potential approach. The finite-volume method is applied to solve the coupled equation for flow, energy, and entropy generation. Simulations are conducted to obtain streamlines, isotherms, Nusselt numbers, entropy generation, and the Bejan number for various magnetic forces (1 ≤ γ ≤ 100) and Rayleigh numbers (10 4 ≤ Ra ≤ 10 6 ). In the absence of gravity, the total entropy generation increases with the increasing magnetic field number, but the average Bejan number decreases. In the gravitational field, the total entropy generation respects the insensitive trend to the change of the magnetic force for low Rayleigh numbers, while it changes significantly for high Rayleigh numbers. When the magnetic field enhances, the share of viscous dissipation in energy losses keeps growing.