近期，中国科学院上海天文台研究人员通过理论推导，深入分析了处于旋转自引力平衡的椭球Boussinesq流体中的热不稳定性问题，首次得到快速旋转椭球形天体内部流体热不稳定性的判据，这一理论不仅可以直接用于理解木星与土星这样显著非球形行星内部的对流动力学，甚至可能被用于研究黑洞吸积盘等极端扁平的旋转流体系统。相关研究成果日前发表在流体力学期刊Physical Review Fluids上，并被美国物理学会选为媒体推荐成果。

热对流运动是天体内部普遍存在的动力学现象，对能量传递、物质输运、磁场产生和演化路径起到支配性的作用。如果一个天体没有自转，可以简单地依据实际温度梯度和绝热温度梯度之间的比较来判断该天体内部是否发生热对流。但实际上，行星与恒星总是自转的，科里奥利力（即地转偏向力）破坏了对流方程的对称性，导致动力系统的临界特性完全改变。

几十年以来，众多学者系统性地研究了球形近似条件下流体从传导态到对流态复杂的线性分岔性质。迄今为止，以上所有研究都忽略了天体由于自转产生的离心力而形成的非球形形状，因此这一系列基于球形近似条件所获得的研究结论的合理性和有效性从未得到任何理论或数值分析的检验。但随着人们对木星、土星和α Eridani等快速自转天体观测和研究的不断加深，椭球形天体内部全局热不稳定性判据变得愈发必要。

在快速旋转非球形稳定分层模型的基础上，中国科学院上海天文台行星物理与磁流体力学课题组孔大力研究员指导澳门科技大学月球与行星科学国家重点实验室博士研究生李文博，运用全局渐进展开的分析手段深入分析了处于旋转自引力平衡的椭球Boussinesq流体中的热不稳定性问题，首次获得了全局热对流线性临界模式和临界参数的解析解，并系统地探讨了非球形形状与对流动力学分岔性质之间的联系。

图1 以木星的扁率情况为例，展示了全局热不稳定性发生时的临界对流模式。

“研究结果证明，对于像木星这样的快速自转且扁率较大的行星，非球形旋转对流的临界参数与球形近似下的结果相比，会出现大幅度变化。”孔大力研究员表示，“如果严格采用与天体自转相自洽的非球形模型，很多快速自转行星与恒星内部的对流输运效率会与前人的预期相差很大。”他还指出，这一研究成果对天体热演化、磁场发电机动力学等问题的影响将在后续同一系列的学术论文中深入探讨。

*文汇独家稿件，转载请注明出处。