The heat losses from pre-insulated double-pipe district heating (DH) systems buried in a homogeneous soil are studied numerically. The study is conducted using the diameter of the pipes and their distance, the size of the insulation, the thermal conductivity ratio between the insulation and the soil, as well as the burial depth of the double-pipe system, as controlling parameters. A computational code based on a control-volume formulation of the finite-difference method has been developed using the open-source framework OpenFOAM with the purpose to compute the heat transfer rate across adjacent solid regions with different thermophysical properties. The main scopes of the study are: (a) to investigate in what measure the geometry and the relative position of the warm and cold pipes, as well as the temperature imbalance, the burial depth and the physical properties of both the insulation and the soil, affect the heat losses; (b) to analyze the existence of an optimal configuration of the DH system by the thermal resistance enhancement viewpoint; and (c) to develop accurate correlating equations for the evaluation of the thermal resistance existing between each pipe and its surroundings, useful for practical thermal engineering applications.