We numerically investigate the interaction between two spherical particles in a nematic liquid crystal mediated by elastic distortions in the orientational order. We pay attention to the cases where two particles with equal radii R0 impose rigid normal anchoring on their surfaces and carry a pointlike topological defect referred to as a hyperbolic hedgehog. To describe the geometry of our system, we use bispherical coordinates, which prove useful in the implementation of boundary conditions at the particle surfaces and at infinity. We adopt the Landau-de Gennes continuum theory in terms of a second-rank tensor order parameter Q(ij) for the description of the orientational order of a nematic liquid crystal. We also utilize an adaptive mesh refinement scheme that has proven to be an efficient way of dealing with topological defects whose core size is much smaller than the particle size. When the two "dipoles," composed of a particle and a hyperbolic hedgehog, are in parallel directions, the two-particle interaction potential is attractive for large interparticle distances D and proportional to D-3 as expected from the form of the dipole-dipole interaction, until the well-defined potential minimum at D approximately 2.46 R0 is reached. For the antiparallel configuration with no hedgehogs between the two particles, the interaction potential is repulsive and behaves as D-2 for D less than or approximately equal 10R0, which is stronger than the dipole-dipole repulsion (approximately D-3 ) expected theoretically as an asymptotic behavior for large D.