振荡矩形腔体中瞬态热流体动力学的分析研究：Oldroyd-B流体的Caputo分数阶导数方法

《Ain Shams Engineering Journal》：Analytical study of transient thermal fluid dynamics in an oscillatory rectangular enclosure: a Caputo fractional derivative approach for Oldroyd-B fluid

　　

随着分数阶导数的出现，分数阶微积分领域在众多领域产生了实质性的应用，包括控制系统、神经科学、电磁学、群体动力学和地震波传播。由于分数阶微积分的显著领域，本研究具体涉及基于Caputo分数阶导数（Caputo fractional derivative）的Oldroyd-B流体（Oldroyd-B fluid）在振荡加热矩形管道（oscillatory heated rectangular duct）内的瞬态自由对流流动行为及热机制。这项新颖的研究展示了Oldroyd-B流体在二维Caputo分数阶导数重要方法背景下的分数阶振荡特性。此外，通过两种高效技术（拉普拉斯变换（Laplace transform）和双有限傅里叶正弦变换（double finite Fourier sine transform））获得了分数阶速度和温度的精确解。在计算精确的分数阶结果后，准备了图形演示以便对相关参数下的温度（temperature）和速度（velocity）分布进行物理可视化。这项基于分数阶的研究表明，分数阶参数提升了温度和速度场。此外，两个流体参数（松弛时间参数（relaxation time parameter）和迟滞时间参数（retardation time parameter））也增强了速度分布。这些发现对工程应用具有实际意义，例如振荡换热器（oscillatory heat exchangers）、微流控冷却技术（microfluidic cooling technologies）和聚合物加工（polymer processing），这些领域需要对流动和传热进行精确调控。

论文解读文章

**研究背景、现存问题与研究动机**

分数阶微积分通过引入非整数阶导数，扩展了经典微积分的框架，在控制系统、神经科学、电磁学等领域展现出广泛应用。其中Caputo分数阶导数因其能够处理与经典导数相似的初始和边界条件，在描述黏弹性流体（viscoelastic fluid）的记忆依赖行为方面具有优势。Oldroyd-B流体作为速率型非牛顿流体，兼具松弛与迟滞特性，广泛应用于聚合物加工、生物流体模拟、润滑剂等领域。然而，先前关于分数阶Oldroyd-B流体流动与传热的研究多局限于简化几何，例如一维垂直板或平行表面间的稳态或非振荡情形，忽略了二维振荡边界效应及记忆效应对热量与动量耦合输运的综合影响。因此，研究人员开展本研究，旨在通过Caputo分数阶导数分析Oldroyd-B流体在振荡矩形腔体中的瞬态自由对流，以更真实地捕捉实际工程（如振荡换热器、微流控冷却）中的复杂热流体动力学行为。

**研究内容与结论**

研究人员建立了垂直振荡矩形管道（宽度?₀，高度h₀）内Oldroyd-B流体的瞬态自由对流模型，考虑重力引起的浮力项，并设定壁面正弦速度（U₀SinΩ?t?）和正弦温度（T?₀SinΩ?t?+T?₀）边界条件。通过无量纲化，得到含Caputo分数阶导数（阶数β₁、β₂）的动量与能量方程。采用Laplace变换与双有限Fourier正弦变换联合推导出速度和温度场的精确解析解（式(30)和式(43)），并通过将分数阶参数设为1（β₁=β₂=1）验证了与经典整数阶解的一致性。参数研究揭示了分数阶参数、Prandtl数、松弛/迟滞时间参数等对温度与速度分布的影响。研究发现：分数阶参数β₂从0.5增至0.9时，温度升高约270%–280%；Prandtl数从2.0增至4.0时，温度降低20%–30%；迟滞时间参数δ从0.2增至0.8时，温度增加15%–25%。速度场方面，β₁从0.7增至1.1使速度提升10%–22%，β₂从0.5增至0.9提升12%–25%；时间周期t使速度提升20%–35%；迟滞参数δ从0.2增至0.9使峰值速度提升15%–30%；频率参数Ω从2π增至12π使速度增加10%–25%；松弛参数ζ使速度提升18%–35%。这些结果表明，分数阶参数的增大降低了记忆效应，增强了传热与动量扩散；而松弛与迟滞参数通过改变流体的弹性响应促进了流动。论文发表在《Ain Shams Engineering Journal》。

**主要关键技术方法**

本研究采用的关键技术方法包括：(1) Caputo分数阶导数模型，用于描述黏弹性流体的记忆依赖行为；(2) 拉普拉斯变换（Laplace transform）与双有限傅里叶正弦变换（double finite Fourier sine transform）联合求解偏微分方程组，推导出速度和温度场的精确闭式解；(3) 无量纲化处理，将物理参数（如格拉晓夫数Gr、Prandtl数Pr）纳入控制方程。本研究为理论分析，无实验样本队列来源。

**研究结果**

- **温度场对分数阶参数β₂的响应**（图3）：随着β₂增大，系统记忆效应减弱，热扩散增强，温度分布显著升高（约270%–280%）。
- **温度场对Prandtl数Pr的响应**（图4）：Pr增大导致流体热导率降低，热量传导能力下降，温度分布降低20%–30%。
- **温度场对迟滞时间参数δ的响应**（图5）：δ增大使松弛时间减小，热能更快积累，温度升高15%–25%。
- **速度场对分数阶参数β₁的响应**（图6）：β₁增大降低黏弹性记忆，减少流动阻力，速度提升10%–22%。
- **速度场对分数阶参数β₂的响应**（图7）：β₂增大削弱迟滞效应，加速动量扩散，速度提升12%–25%。
- **速度场对时间周期t的响应**（图8）：t增大增强动量扩散与加速，速度提升20%–35%。
- **速度场对迟滞时间参数δ的响应**（图9）：δ增大减弱弹性恢复力，降低流动阻力，速度提升15%–30%。
- **速度场对频率参数Ω的响应**（图10）：振荡频率增大提高惯性力与能量交换，速度提升10%–25%。
- **速度场对松弛时间参数ζ的响应**（图11）：ζ增大使流体变形时间延长，动量传递增强，速度提升18%–35%。

**讨论总结**

本研究发现，分数阶参数、弛豫与迟滞参数、振荡频率等均可显著调控温度与速度分布。这些结果对优化含黏弹性流体的振荡热系统具有指导意义：例如，调节分数阶参数可有效控制热效率与流动稳定性；松弛与迟滞参数的调节可用于聚合物加工中流变特性的定制；频率参数的提升有助于加速混合与换热。该研究为工程中记忆依赖流体输运机制的设计提供了理论依据。

**研究结论翻译**

1. 温度场是分数阶参数的递增函数。较高的分数阶参数降低系统的热记忆，从而改善热扩散并放大热响应，导致温度升高。
2. 改进的迟滞时间参数增强了流体运动。
3. 温度随迟滞时间参数增大而升高，随Prandtl数增大而降低。
4. 分数阶参数增大时，流体运动随之加强。速度的增强归因于黏弹性记忆的减弱，这降低了动量传递的阻力并促进更大的流体加速度。
5. 时间周期的强烈影响显著提升了流动速度。
6. 流体速度是松弛时间参数和频率参数的递增函数。
7. 与普通导数相比，Caputo分数阶算子为本研究提供了更显著的结果。
8. 本模型基于二维充分发展流动等假设，忽略了黏性耗散、热辐射和内热生成，且未经过实验验证。未来工作可纳入更多物理因素，并通过实验测量或计算分析进行模型验证。