The paper is an overview of the main contributions of a Bulgarian team of researchers to the problem of finding the possible structures and waves in the open nonlinear heat conducting medium, described by a reaction-diffusion equation. Being posed and actively worked out by the Russian school of A. A. Samarskii and S.P. Kurdyumov since the seventies of the last century, this problem still contains open and challenging questions.of A.A. Samarskii and M.I. Sobol in 1963 [51]. The problem of localization for quasilinear equations with a source is posed by S.P. Kurdyumov [42] in 1974. The works of I.M. Gelfand, A.S. Kalashnikov, the scientists of the school of A.A. Samarskii and S.P. Kurdyumov are devoted to the challenging physical and mathematical problems, related with this model and its generalizations. Among them are: localization in space of the process of burning, different types of blow-up, arising of structurestraveling and standing waves, complex structures with varying degrees of symmetry. The combination of the computational experiment with the progress in the qualitative and analytical methods of the theory of ordinary and partial differential equations, the Lie and the Lie-Bäcklund group theory, has been crucial for the success of these investigations. The book [50] contains many of these results, achieved to 1986, in the review [34] there are citations of later works.A special part of these investigations is devoted to finding and studying different kinds of selfsimilar and invariant solutions of equation (1) with power nonlinearities :This choice is suggested by the following reasoning.First, such temperature dependencies are usual for many real processes [5], [54], [57]. For example, when σ i = σ = 2.5, β ≤ 5.2, equation (1) describes thermo-nuclear combustion in plasma in the case of electron heat-conductivity; the parameters σ = 0, 2 ≤ β ≤ 3 correspond to the models of autocatalytic processes with diffusion in the chemical reactors; σ≈6.5 corresponds to the radiation heat-conductivity of the high-temperature plasma in the stars, and so on.Second, it is shown in [25], that in the class of power functions the symmetry of equation (1) is maximal in some sense -the equation admits a rich variety of invariant solutions. In general, almost all of the dissipative structures known so far are invariant or partially invariant solutions of nonlinear equations. The investigations of the dissipative structures provide reasons to believe that the invariant solutions describe the attractors of the dissipative structures' evolution and thus they characterize important internal properties of the nonlinear dissipative medium.Third, this rich set of invariant solutions of equation (1) with power nonlinearities is necessary for the successful application of the methods for investigating the same equation in the case of more general dependencies k i (u), Q(u). By using the methods of operator comparison [29] and stationary states [32] it is possible to analyze the properties of the solutions (such as localization, blo...