In this thesis we have focussed on two types of axially symmetric ion trap mass analysers viz., the quadrupole ion trap mass analyser and the toroidal ion trap mass analyser. We have undertaken three numerical studies in this thesis, one study is on the quadrupole ion trap mass analysers and two studies are on the toroidal ion trap mass analysers. The first study is related to improvement of the sensitivity of quadrupole ion trap mass analysers operated in the resonance ejection mode. In the second study we have discussed methods to determine the multipole coefficients in the toroidal ion trap mass analysers. The third study investigates the stability of ions in the toroidal ion trap mass analysers. The first study presents a technique to cause unidirectional ion ejection in a quadrupole ion trap mass spectrometer operated in the resonance ejection mode. In this technique a modified auxiliary dipolar excitation signal is applied to the endcap electrodes. This modified signal is a linear combination of two signals. The first signal is the nominal dipolar excitation signal which is applied across the endcap electrodes and the second signal is the second harmonic of the first signal, the amplitude of the second harmonic being larger than that of the fundamental. We have investigated the effect of the following parameters on achieving unidirectional ion ejection: primary signal amplitude, ratio of amplitude of second harmonic to that of primary signal amplitude, different operating points, different scan rates, different mass to charge ratios and different damping constants. In all these simulations unidirectional ejection of destabilized ions has been successfully achieved. The second study presents methods to determine multipole coefficients for describing the potential in toroidal ion trap mass analysers. Three different methods have been presented to compute the toroidal multipole coefficients. The first method uses a least square fit and is useful when we have ability to compute potential at a set of points in the trapping region. In the second method we use the Discrete Fourier Transform of potentials on a circle in the trapping region. The third method uses surface charge distribution obtained from the Boundary Element Method to compute these coefficients. Using these multipole coefficients we have presented (1) equations of ion motion in toroidal ion traps (2) the Mathieu parameters in terms of multipole coefficients and (3) the secular frequency of ion motion in these traps. It has been shown that the secular frequency obtained from our method has a good match with that obtained from numerical trajectory simulation. The third study presents stability of ions in practical toroidal ion trap mass analysers. Here we have taken up for investigation four geometries with apertures and truncation of electrodes. The stability is obtained in UDC-VRF plane and later this is converted into A-Q plane on the Mathieu stability plot. Though the plots in terms of Mathieu parameters for these structures are qualitatively similar to the corresponding plot of linear ion trap mass analysers, there is a significant difference. The stability plots of these have regions of nonlinear resonances where ion motion is unstable. These resonances have been briefly investigated and it is proposed that they occur on account of hexapole and octopole contributions to the field in these toroidal ion traps.