A distributed nonlinear mathematical model for investigation of regenerative heat exchangers of both a continuous and periodic operation is described in the paper. The non-iterative numerical integration scheme for conjugate unsteady heat exchange problem of one dimensional flows and two dimensional matrix wall conductivity is developed. Case study of a regenerative heat exchanger with a rotary ceramic matrix is presented. The range of optimum rotation rates of the regenerator providing the greatest calorific efficiency is determined.