旋转Rayleigh-Benard问题的Lorenz模型及数值模拟.docxVIP

旋转Rayleigh-Benard问题的Lorenz模型及数值模拟

一、物理背景与模型推导

旋转Rayleigh-Benard（R-RB）对流是流体力学中典型的非线性系统，指在旋转坐标系下，流体层受垂直温度梯度驱动形成的对流现象，广泛存在于大气环流、地幔对流等自然过程中。当旋转角速度Ω与温度梯度引起的浮力达到特定平衡时，系统会从稳态过渡到混沌态，其动力学行为可通过简化的Lorenz模型描述。

模型推导基于以下假设：

流体为不可压缩牛顿流体，满足Boussinesq近似；

旋转沿垂直方向（z轴），科里奥利力为主要旋转效应；

对流结构呈轴对称，用柱坐标（r,φ,z）下的径向速度u、轴向速度w、温度扰动θ表征；

忽略高次非线性项，保留关键耦合项（速度-温度耦合、科里奥利力-速度耦合）。

二、旋转Rayleigh-Benard系统的Lorenz模型

在传统Lorenz模型基础上引入旋转项（科里奥利力修正），最终无量纲化方程为：

\begin{cases}\dot{x}=\sigma(y-x)+\alphaz\\\dot{y}=rx-y-xz\\\dot{z}=xy-bz\end{cases}

1.变量定义

x：表征径向对流速度的振幅；

y：表征温度扰动的振幅；

z：表征旋转效应与轴向流动的耦合项；

t：无量纲时间（与热扩散时间尺度相关）。

2.参数物理意义

参数

表达式

物理意义

典型取值

σ

ν/κ

普朗特数（动量扩散率/热扩散率）

10（流体如空气）

r

Ra/Ra_c

相对瑞利数（实际Ra/临界Ra）

28（混沌阈值）

α

Ω2d?/νκ

旋转数（Ω：角速度，d：流体层厚度）

0~5（弱旋转→强旋转）

b

4/(1+Pr)

几何修正系数

8/3（Pr=0.7时）

三、数值模拟方案

1.数值方法选择

采用四阶龙格-库塔（RK4）方法求解常微分方程组，该方法精度高（局部截断误差O(h?)），且能稳定捕捉混沌系统的分岔与吸引子结构，迭代公式为：

\begin{cases}k_1=f(t_n,X_n)\\k_2=f(t_n+h/2,X_n+hk_1/2)\\k_3=f(t_n+h/2,X_n+hk_2/2)\\k_4=f(t_n+h,X_n+hk_3)\\X_{n+1}=X_n+\frac{h}{6}(k_1+2k_2+2k_3+k_4)\end{cases}

其中X=(x,y,z)^T，h为时间步长（取h=0.01以平衡精度与效率）。

2.初始条件与边界

初始扰动：X_0=(1,0,0)（微小速度扰动触发对流）；

模拟时长：t∈[0,1000]（前500步为暂态，后500步用于分析稳态行为）。

3.关键分析指标

时间序列：x(t)、y(t)曲线，观察振幅波动与周期性；

相空间吸引子：x-y、x-z、y-z投影图，判断系统稳态（定点/周期/混沌）；

Lyapunov指数（LE）：计算最大LE，LE0时系统呈混沌态。

四、模拟结果与分析

以典型参数（σ=10，b=8/3，r=28）为例，分析不同旋转强度α的影响：

1.弱旋转（α=0.5）

时间序列：x(t)呈非周期波动，振幅在[-20,20]间不规则变化；

相空间：形成类似传统Lorenz吸引子的“蝴蝶翼”结构，但翼面更紧凑（图1）；

LE：最大LE≈0.92，系统处于混沌态，旋转抑制了部分无序性。

2.强旋转（α=4.0）

时间序列：x(t)振幅减小（[-5,5]），出现准周期振荡（图2）；

相空间：吸引子从混沌态过渡到2维环面（周期2运动）；

LE：最大LE≈0，系统进入周期态，旋转主导了流体运动的有序性。

3.临界旋转（α≈2.3）

系统发生倍周期分岔：随α增大，混沌吸引子逐渐破碎，最终坍缩为定点吸引子（稳态对流），对应物理上“旋转抑制对流混沌”的现象。

五、结论

旋转强度α是R-RB系统动力学行为的关键控制参数：α较小时，系统呈混沌态；α超过临界值后，旋转抑制非线性扰动，系统过渡到周期态或稳态；

修正后的Lorenz模型能有效描述旋转效应，其相空间结构与Lyapunov指数可定量反映对流从有序到混沌的转变；

数值模拟表明，RK4方法能准确捕捉系统分岔过程，为大气、地球物理等领域的旋转对流研究提供简化模型支持。

（注：文中图表需通过编程实现，可使用Python（scipy.integrate.solve_ivp+matplotlib）或MATLAB（ode45+