We studied a multi-orbital Hubbard model at half-filling for two and three orbitals per site on a two-site cluster via full exact diagonalization, in a wide range for the onsite repulsion U, from weak to strong coupling, and multiple ratios of the Hund coupling JH
to U. The hopping matrix elements among the orbitals were also varied extensively. At intermediate and large U, we mapped the results into a Heisenberg model. For two orbitals per site, the mapping is into a S=1 Heisenberg model where by symmetry both nearest-neighbor (S
i.S
j) and (S
i.S
j)2 are allowed, with respective couplings J1
and J2
. For the case of three orbitals per site, the mapping is into a S=3/2 Heisenberg model with (S
i.S
j), (S
i.S
j)2, and (S
i.S
j)3 terms, and respective couplings J1
, J2
, and J3
. The strength of these coupling constants in the Heisenberg models depend on the U, JH
, and hopping amplitudes of the underlying Hubbard model. Our study provides a first crude estimate to establish bounds on how large the ratios J2/J1
and J3/J1
can be. We show that those ratios appear rather limited and, as a crude guidance, we conclude that J2/J1
is less than 0.4 and J3/J1
is less than 0.2, establishing bounds on effective models for strongly correlated Hubbard systems. Moreover, the intermediate Hubbard U regime was found to be the most promising to enhance J2/J1
and J3/J1
.