Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertorOpen Access