The study of electromagnetic scattering in optics and nanophotonics is crucial for understanding complex nanostructures and optical devices. However, numerical analysis of scattering spectra, even for simple shape nanocavities, presents significant computational challenges. This paper focuses on simplifying the simulation process by utilizing system symmetries and the separation of variables, which reduces the dimension of the problem. Specifically, we present a practical guide to efficiently simulate linear and nonlinear scattering problems in COMSOL Multiphysics for axisymmetric objects. This includes computing the scattering cross-section, multipolar decomposition, optical forces, and second harmonic generation (SHG). We show that by reducing the calculation time we are able to analyze the generation spectra of SH under varying geometrical parameters of nanoantennas across a broad range. Additionally, we study the increase in SHG within a high-Q state at the band edge in a chain of dielectric discs with lengths ranging from 5 to 20 particles. We also accompany the provided guide with the ready-to-run COMSOL models.