If only there was a Wikipedia article - https://en.wikipedia.org/wiki/Saturn’s_hexagon

A hypothesis developed at Oxford University is the hexagon forms where there is a steep latitudinal gradient in the speed of the atmospheric winds in Saturn’s atmosphere.[22] Similar regular shapes were created in a laboratory when a circular tank of liquid was rotated at different speeds at its centre and periphery. The most common shape was six sided, but shapes with three to eight sides were also produced. The shapes form in an area of turbulent flowbetween the two different rotating fluid bodies with dissimilar speeds.[22][23]A number of stable vortices of similar size form on the slower (south) side of the fluid boundary, and these interact with each other to space themselves out evenly around the perimeter. The presence of the vortices influences the boundary to move northward where each is present, and this gives rise to the polygon effect.[23] Polygons do not form at wind boundaries unless the speed differential and viscosity parameters are within certain margins and thus absent at other likely places, such as Saturn’s south pole or the poles of Jupiter.

Other researchers claim that lab studies exhibit vortex streets, a series of spiraling vortices not observed in Saturn’s hexagon. Simulations show a shallow, slow, localized meandering jetstream in the same direction as Saturn’s prevailing clouds are able to match the observed behaviors of Saturn’s hexagon with the same boundary stability.[24]

Developing barotropic instability of Saturn’s North Polar hexagonal circumpolar jet (Jet) plus North Polar vortex (NPV) system produces a long-living structure akin to the observed hexagon, which is not the case of the Jet-only system, which was studied in this context in a number of papers in literature. The NPV, thus, plays a decisive dynamical role to stabilize hexagon jets. The influence of moist convection, which was recently suggested to be at the origin of Saturn’s NPV system in the literature, is investigated in the framework of the barotropic rotating shallow water model and does not alter the conclusions.[25]

A 2020 mathematical study at the California Institute of Technology found that a stable geometric arrangement of the polygons can occur on any planet when a storm is surrounded by a ring of winds turning in the opposite direction to the storms itself, called an anticyclonic ring, or anticyclonic shielding.[26][27]Such shielding creates a vorticity gradient in the background of a neighbor cyclone, causing mutual rejection between the cyclones (similar to the effect of beta-drift). Although apparently shielded, the polar cyclone on Saturn cannot hold a polygonal pattern of circumpolar cyclones such as Jupiter’s due to the bigger size and slower wind speed of Saturn’s polar cyclone, so the side-adjacent vortices and deep barotropic instability (Cassini’s wind speed measurements preclude shallower barotropic instability at least at the time of the Cassini encounter), or possibly baroclinic instabilities remain as the most viable explanations for Saturn’s sustained hexagon.[28]