Natural convection of liquid metal in a cubical enclosure under an external magnetic field was investigated by three-dimensional numerical analyses. The system parameters were Ra = 105 and 106, Pr = 0.025, and Ha = 0–1000. One vertical wall of the cubical enclosure was heated, and the opposing vertical wall was cooled, both isothermally; the other four walls were thermally insulated. A uniform horizontal magnetic field was applied parallel to the heated and cooled walls. At Ra = 105 and Ha = 50, the average Nusselt number on the heated wall attained almost the maximum value and was greater than that at Ha = 0. The velocity vectors along the vertical walls, and those along the horizontal planes, were rectified in a two-dimensional way at Ha = 50 or over, and the average Nusselt number decreased gradually for higher values of the Hartmann number. Similar characteristics were obtained at Ra = 106. The agreement with our earlier experiments was moderately good.