Perturbations of Zonal and Tesseral Harmonics on Frozen Orbits of Charged Satellites

The objective of this research is to address the issue of frozen orbits in charged satellites by incorporating geopotential zonal harmonics up to J6 and the initial tesseral harmonics. The employed model starts from the first normalized Hamiltonian to calculate specific sets of long-term frozen orbits for charged satellites. To explore the frozen orbits acquired, a MATHEMATICA CODE is developed. The investigation encompasses extensive variations in orbit altitudes by employing the orbital inclination and argument of periapsis as freezing parameters. The determined ranges ensuring frozen orbits are derived from the generated figures. Three-dimensional presentations illustrating the freezing inclination in relation to eccentricity, argument of periapsis, and semi-major axis parameters are presented. Additional three-dimensional representations of the phase space for the eccentricity vector and its projection onto the nonsingular plane are examined. In all investigated scenarios, the impacts of electromagnetic (EM) field perturbations on the freezing parameters of a charged satellite are demonstrated.