Numerical solution of 2D thermoacoustic problem