In this paper, we study the dynamics of a reaction-diffusion Cellular Nonlinear Network (RD-CNN) nodel in which the reaction term is represented by Brusselator cell. We investigate the RD-CNN dynamics by means of describing function method. Comparison with classical results for Brusselator equation is provided. Then we introduce a new RD-CNN model with memristor coupling, for which the edge of chaos regime in the parameter space is determined. Numerical simulations are presented for obtaining dynamic patterns in the RD-CNN model with memristor coupling.