▎ 摘　　要

Rigorous mathematical models of nonlinear terahertz (THz) graphene-based devices are developed by solving nonlinear Maxwell's equations with boundary conditions simultaneously with a model of graphene surface conductivity as a nonlinear function on the time-varying electric field. The numerical approach is based on the projection method. Nonlinear phenomena, related to the parametric instability due excitations by the incident pumping TEM-wave in the multilayer graphene-dielectric nanostructures, are investigated numerically at THz frequencies.