[论文解读] An exact axisymmetric spiral solution of incompressible 3D Euler equations

Liang Sun|arXiv (Cornell University)|Jan 29, 2011

Computational Fluid Dynamics and Aerodynamics参考文献 5被引用 5

一句话总结

本文通过分离变量法，为三维不可压缩欧拉方程提出了一类通用的精确轴对称螺旋解，得到了有界且具有物理意义的涡旋结构，如多层伞状涡旋和多抛物面涡旋。该解能够捕捉伯努利面上的螺旋流体轨迹，并提供具有三个可调参数的台风景、龙卷风和中尺度涡旋的显式模型，其参数可匹配最大风速、螺旋角等可观测流动特征。

ABSTRACT

Spiral structure is one of the most common structures in the nature flows. A general steady spiral solution of incompressible inviscid axisymmetric flow was obtained analytically by applying separation of variables to the 3D Euler equations. The solution, depending on 3 parameters, describes the spiral path of the fluid material element on the Bernoulli surface, whereas some new exact solutions were obtained to be bounded within the whole region. The first one is a continued typhoon-like vortex solution, where there are two intrinsic length scales. One is the radius of maximum circular velocity $r_m$, the other is the radius of the vortex kernel $r_k=\sqrt{2}r_m$. The second one is a multi-planar solution, periodically in $z$-coordinate. Within each layer, the solution is a umbrella vortex similar to the first one. The third one is also a multi-planar solution in $z$-coordinate. In each layer, it is a combination of two independent solutions like the Rankine vortex, which is also finite but discontinued for either vertical or horizontal velocity. The fourth one is a multi-paraboloid vortex solution finite for $z$-coordinate but infinite for $r$-coordinate. Besides, some classical simple solutions (Rankine vortex, Batchelor vortex, Hill spherical vortex, etc.) were also obtained. The above explicit solutions can be applied to study the radial structure of the typhoon, tornados and mesoscale eddies. Both the solutions and approaches used here could also be applied to other complex flows by the Navier-Stokes equations.

研究动机与目标

• 在轴对称、定常、无粘流动条件下，推导三维不可压缩欧拉方程的精确解析解。
• 通过强制径向、角向和轴向速度分量之间的耦合，克服先前解存在无界或解耦的局限性。
• 构建物理上合理、有限且有界的涡旋解，以捕捉热带气旋和龙卷风中观测到的自然螺旋流动模式。
• 提出一种系统化方法，利用分离变量法生成具有可调参数的新精确解，使其匹配真实世界涡旋特征。

提出的方法

• 在柱坐标系 (r, θ, z) 下对三维不可压缩欧拉方程应用分离变量法，假设轴对称性和定常流动。
• 将速度分量 (Vr, Vθ, Vz) 表示为仅与 r 和 z 有关的函数，其中 Vθ 与 Ψ/r 成正比，Ψ 为流函数。
• 推导广义贝特拉米流条件，即涡度与速度平行，从而将非线性欧拉方程简化为可解系统。
• 引入径向函数 R(r) 和轴向函数 H(z)，其中伯努利函数 B = −λR(r)H(z) 定义流动结构。
• 利用质量守恒和动量方程，推导出 R(r) 的二阶常微分方程，在特定参数约束下可解。
• 将不同径向区域（如内层和外层）的解组合，构造在界面半径处导数连续的分段连续解。

实验结果

研究问题

• RQ1 在定常、无粘、轴对称条件下，能否为三维不可压缩欧拉方程导出一类通用的精确轴对称螺旋解？
• RQ2 哪些物理约束和参数依赖关系可产生有界、有限的涡旋解，使其与真实大气和海洋涡旋相似？
• RQ3 如何将分离变量法应用于耦合的非线性欧拉方程，以生成在伯努利面上具有螺旋轨迹的显式解？
• RQ4 新涡旋解的结构和运动学特性是什么？与经典模型（如兰金或巴托勒涡旋）相比有何异同？
• RQ5 所导出的解能否用于建模和解释台风与龙卷风中的径向速度分布和角动量分布？

主要发现

• 为三维不可压缩欧拉方程导出了通用的精确螺旋解，描述流体质点在伯努利面上沿螺旋路径运动。
• 获得了三种新的有界涡旋解：一种具有两个本征长度尺度（rm 和 rk = √2 rm）的多平面伞状涡旋；一种结合两个兰金型涡旋的多平面涡旋；以及一种在 z 方向有界但在 r 方向无界的多抛物面涡旋。
• 该解包含兰金涡旋、巴托勒涡旋和希尔球形涡旋作为特例，验证了其与已知经典解的一致性。
• 角动量表达式为 m = rVθ = −λR(r)H(z)，R(r) 的径向依赖关系允许对最大风速出现在 rm 处的观测涡旋结构进行建模。
• 通过调节三个参数（rm：最大风速半径，Vθm：最大切向速度，k/λ：螺旋角），模型可拟合真实观测数据，从而实现对台风和龙卷风的流动结构分类。
• 该方法可扩展至非定常粘性流和其它复杂系统，如 f 平面上的地球流或磁流体动力学流。