Accurate and robust spatial orders are ubiquitous in living systems. In 1952, Turing proposed a general mechanism for pattern formation exemplified by a reaction–diffusion model with two chemical species in a large system. However, in small biological systems such as a cell, the existence of multiple Turing patterns and strong noise can lower the spatial order. Recently, a modified reaction–diffusion model with an additional chemical species is shown to stabilize the Turing pattern. Here, we study non-equilibrium thermodynamics of this three-species reaction–diffusion model to understand the relationship between energy cost and the performance of self-positioning. By using computational and analytical approaches, we show that beyond the onset of pattern formation the positioning error decreases as energy dissipation increases. In a finite system, we find that a specific Turing pattern exists only within a finite range of total molecule number. Energy dissipation broadens this range, which enhances the robustness of Turing pattern against molecule number fluctuations in living cells. The generality of these results is verified in a realistic model of the Muk system underlying DNA segregation in Escherichia coli, and testable predictions are made for the dependence of the accuracy and robustness of the spatial pattern on the ATP/ADP ratio.