In this paper, a new formulation based on the method of fundamental solutions for two/three-dimensional steady-state heat conduction problems involving internal curved line/surface heat sources is presented. Arbitrary shapes and non-uniform intensities of the curved heat sources can be modeled by an assemblage of several parts with quadratic variations. The presented mesh-free modeling does not require any internal points as in domain methods. Four numerical examples are studied to verify the validity and efficiency of the proposed method. Our analyses have shown that the presented mesh-free formulation is very efficient in comparison with conventional boundary or domain solution techniques.