A time-dependent two-dimensional photochemical-dynamical coupling gravity wave model of sodium layer is developed, which combines the sodium photochemical theory, a time-dependent two-dimensional atmospheric photochemical model, a two-dimensional gravity wave model, and the International Reference Ionosphere model (IRI-95)with the diabatic process induced by photochemical reactions and the transport of chemical species by gravity waves included. The pseudospectral method is used in the horizontal direction, the finite difference approximations are used in vertical direction z and time t. And FICE method is used to solve the model. The simulation results indicate that intense perturbations of the sodium layer can be induced by the propagation of gravity waves. The results are consistent with the observations.