We consider heat transfer in an infinite two-dimensional square harmonic scalar lattice lying in a viscous environment and subjected to a heat source. The basic equations for the particles of the lattice are stated in the form of a system of stochastic ordinary differential equations. We perform a continualization procedure and derive an infinite system of linear partial differential equations for covariance variables. The most important results of the paper are the deterministic differential-difference equation describing non-stationary heat propagation in the lattice and the analytical formula in the integral form for its steady-state solution describing kinetic temperature distribution caused by a point heat source of a constant intensity. The comparison between numerical solution of stochastic equations and obtained analytical solution demonstrates a very good agreement everywhere except for the main diagonals of the lattice (with respect to the point source position), where the analytical solution is singular.Keywords ballistic heat transfer · 2D harmonic scalar lattice · kinetic temperature 1 Introduction At the macroscale, Fourier's law of heat conduction is widely and successfully used to describe heat transfer processes. However, recent experimental observations demonstrate that Fourier's law is violated at the microscale and nanoscale, in particular, in low-dimensional nanostructures [1][2][3][4][5][6], where the ballistic heat transfer is realized. The anomalous heat transfer also may be related with the spontaneous emergence of long-range correlations; the latter is typical for momentum-conserving systems [7,8]. The simplest theoretical approach to describe the ballistic heat propagation is to use harmonic lattice models. In some cases such models allow one to obtain the analytical description of thermomechanical processes in solids [9][10][11][12][13][14][15]. In the literature, the problems concerning heat transfer in harmonic lattices are mostly considered