[论文解读] Shifted polynomials in a convection problem

Ioana Dragomirescu|ArXiv.org|Sep 14, 2007

Differential Equations and Numerical Methods参考文献 18被引用 3

一句话总结

本论文采用基于移位勒让德多项式和切比雪夫多项式的伽辽金法，研究了在垂直方向变化的重力场下，水平流体层从下加热时的线性稳定性。结果表明，重力场减小会扩大稳定域，非均匀重力场下瑞利数临界值显著提高，且移位切比雪夫多项式方法在求解该特征值问题方面比勒让德多项式更有效。

ABSTRACT

The onset of convection in a horizontal layer of fluid heated from below in the presence of a gravity field varying across the layer is investigated. The eigenvalue problem governing the linear stability of the mechanical equilibria of the fluid layer in the case of free boundaries is solved using a Galerkin method based on shifted polynomials (Legendre and Chebyshev polynomials).

研究动机与目标

• 分析垂直方向变化的重力场对水平流体层热对流起始的影响。
• 采用基于正交多项式的谱方法求解自由滑移边界条件下线性稳定性特征值问题。
• 比较移位勒让德多项式与移位切比雪夫多项式在逼近临界瑞利数方面的精度与效率。
• 评估非均匀重力场（特别是 h(z) = -z, -z², 和 z²-2z）对稳定域的影响。
• 通过经典结果和先前的三角函数伽辽金方法验证该方法的有效性。

提出的方法

• 研究构建了由变量系数微分方程组控制的速度与温度扰动的两点边值问题，其系数因重力变化 H(z) = 1 + εh(z) 而变化。
• 采用移位正交多项式（勒让德与切比雪夫）作为基函数，应用伽辽金法，这些基函数满足齐次边界条件。
• 引入新变量 Ψ = (D² - a²)W，将四阶系统转化为关于 U = (W, Ψ, Θ) 的三个耦合的一阶常微分方程组。
• 试函数由正交多项式 Q_i(z) 的积分构造，构成 H₀¹(0,1) 中的完备集，且满足零边界条件。
• 推导出特征值问题的弱形式，得到涉及基函数内积的矩阵形式的特征方程。
• 针对 n = 4 个基函数进行数值求解，并通过比较不同 ε 和 a 值下的结果评估收敛性。

实验结果

研究问题

• RQ1 非均匀重力场（特别是 H(z) = 1 + εh(z)）如何影响具有自由边界的流体层中对流起始的临界瑞利数？
• RQ2 在该对流设置中，采用移位勒让德多项式还是移位切比雪夫多项式作为正交多项式基，能提供更精确且高效的数值逼近？
• RQ3 与经典均匀重力场情况相比，重力场减小（如 h(z) = -z, -z²）在多大程度上稳定了流体层？
• RQ4 基于多项式基的伽辽金方法所得结果与经典解及先前的三角函数伽辽金方法结果相比如何？
• RQ5 在考虑的变重力配置与自由边界条件下，稳定性交换原理是否成立？

主要发现

• 当 ε = 0.75 且 h(z) = z² - 2z 时，临界瑞利数达到 1251.178（SCP）和 1276.05（SLP），相比经典情况（R² ≈ 657.51）显著提高，表明稳定性明显增强。
• 当 h(z) = -z 时，瑞利数从 ε = 0 时的 657.51 上升至 ε = 0.33 时的 787.36，表明重力减小可增强稳定性。
• 当 h(z) = -z² 时，瑞利数从 ε = 0 时的 657.51 上升至 ε = 0.75 时的 993.51，证实了类似的稳定趋势。
• 移位切比雪夫多项式方法始终产生更接近经典值的结果，并表现出优于勒让德方法的收敛特性。
• 尽管勒让德多项式在伽辽金方法中应用广泛，但本研究发现基于移位切比雪夫多项式的数值方法更为有效。
• 当 ε = 0.0 时，数值结果与经典值一致（R² ≈ 657.51），验证了该方法的准确性。