Objectives. To develop mathematical model representations of the energy effect in non-cylindrical domains having a thermally insulated moving boundary; to introduce a new boundary condition for thermal insulation of a moving boundary both for locally equilibrium heat transfer processes in the framework of classical Fourier phenomenology, as well as for more complex locally non-equilibrium processes in the framework of Maxwell–Cattaneo–Lykov–Vernott phenomenology, taking into account the finite rate of heat propagation into analytical thermophysics and applied thermomechanics; to consider an applied problem of analytical thermophysics according to the theory of thermal shock for a domain with a moving thermally insulated boundary free from external and internal influences; to obtain an exact analytical solution of the formulated mathematical models for hyperbolic type equations; to investigate the solutions obtained using a computational experiment at various values of the parameters included in it; to describe the wave nature of the kinetics of the processes under consideration.Methods. Methods and theorems of operational calculus, Riemann–Mellin contour integrals are used in calculating the originals of complex images with two branch points. A new mathematical apparatus for the equivalence of functional constructions for the originals of the obtained operational solutions, which considers the computational difficulties in finding analytical solutions to boundary value problems for equations of hyperbolic type in the domain with a moving boundary, is developed.Results. Developed mathematical models of locally nonequilibrium heat transfer and the theory of thermal shock for equations of hyperbolic type in a domain with a moving thermally insulated boundary are presented. It is shown that, despite the absence of external and internal sources of heat, the presence of a thermally insulated moving boundary leads to the appearance of a temperature gradient in the domain and, consequently, to the appearance of a temperature field and corresponding thermoelastic stresses in the domain, which have a wave character. A stochastic analysis of this energy effect forms the basis for a proposed transition of the kinetic energy of a moving thermally insulated boundary into the thermal energy of the domain. The presented model representations of the indicated effect confirmed the stated assumption.Conclusions. Mathematical models for locally nonequilibrium heat transfer processes and the theory of thermal stresses are developed and investigated on the basis of constitutive relations of the theory of thermal shock for equations of hyperbolic type in a domain with a thermally isolated moving boundary. A numerical experiment is presented to demonstrate the possibility of transiting from one form of analytical solution of a thermophysical problem to another equivalent form of a new type. The described energy effect manifests itself both for parabolic type equations based on the classical Fourier phenomenology, as well as for hyperbolic type equations based on the generalized Maxwell–Cattaneo–Lykov–Vernott phenomenology.