A thermoacoustic sound generation model, based on the classical balance equations of the continuum mechanics, is here developed for the cylindrical and the spherical thermoacoustic wave generation. In both geometries, the model considers an arbitrary multilayered structure, where each layer can be fluid or solid and it is characterized by the fully coupled thermo-visco-acoustic response. It means that the viscous behavior and the thermal conduction are considered in each layer. The model is based on a unified representation of cylindrical or spherical thermoacoustic waves, which is valid for both fluid and solid phases. Thanks to the continuity of temperature, particle velocity, normal stress and heat flux between adjacent layers, the model can be implemented by means of a versatile matrix approach, allowing flexible analysis and design of cylindrical or spherical thermophones. Any thermoacoustic variable can be determined at any position, any frequency and for any input power. The results are compared with the models already existing in the literature and the underlying physics is thoroughly discussed. The analysis is focused on the better understanding of the thermoacoustic generation with application to the state of the art of the thermophone technology.