《Results in Engineering》：RSM and FDM Role of Thermo-Diffusion and Heat source in MHD-Driven Heat and Mass transfer in Unsteady Casson Fluid flow over a porous vertical plate
编辑推荐：
本研究对非稳态辐射磁流体动力学（MHD）Casson流体在Soret效应和化学反应共同影响下的流动进行了数值研究，重点关注传热和传质过程。研究人员结合适当的初始和边界条件，建立了能量、动量和组分浓度的控制方程。通过使用有限差分隐式方法在MATLAB软件中求解偏
本研究对非稳态辐射磁流体动力学（MHD）Casson流体在Soret效应和化学反应共同影响下的流动进行了数值研究，重点关注传热和传质过程。研究人员结合适当的初始和边界条件，建立了能量、动量和组分浓度的控制方程。通过使用有限差分隐式方法在MATLAB软件中求解偏微分方程，扩展并验证了结果，并通过图形和表格研究了关键流动参数的影响。研究人员分析了不同物理参数对温度、速度和浓度分布的影响，以及不同参数对表面摩擦系数、Nusselt数和Sherwood数的影响。此外，利用响应面法（RSM）增强传输特性。研究人员获得了表面摩擦系数和Nusselt数的高度精确的二次回归模型，其决定系数（R2）分别大于99.95%。同时，化学反应通过消耗溶液降低了流体浓度，牛顿加热降低了表面传热速率。该研究强调了Casson流体流动特性在涉及受传热和传质影响的非牛顿流体的工程和工业应用中的重要性。
**论文解读文章**
**研究背景与问题**
在工程和工业应用中，磁流体动力学（MHD）流过多孔介质的现象广泛存在于天体物理、地球物理和工业设备中，其中磁场与导电流体相互作用，影响传热和传质。Casson流体模型能够有效描述血液、印刷油墨、巧克力、油漆等具有屈服应力的非牛顿流体的流变特性，比牛顿流体假设更贴近实际运输过程。然而，现有研究大多集中于单一或少数参数的影响，缺乏对非稳态Casson流体在MHD效应、热辐射、热扩散（Soret效应）、内热源、粘性耗散、化学反应以及多孔介质阻力等多因素耦合作用的全面理解。此外，传统数值分析常采用单因素方法，无法探究参数间的相互作用或优化操作条件。因此，需要建立一个综合数学模型，并结合统计优化方法，以深入分析非牛顿流体在复杂物理条件下的传热传质机理。
**研究内容与结论**
该研究数值探究了非稳态MHD Casson流体流经多孔垂直板的流动，整合了表面热流、可变吸力、热辐射、粘性耗散、焦耳热以及二阶化学反应等影响因素。通过有限差分法（FDM）求解控制方程，并利用响应面法（RSM）建立二次回归模型，优化传输特性。主要结论表明：较高的渗透率（K）和Grashof数（Gr, Gm）能提高流体速度；而Casson参数（β）和磁场参数（M）增大则会降低速度。温度随Prandtl数（Pr）、热源参数（Q）和辐射参数（R）增大而降低，但随Eckert数（Ec）增大而升高。浓度随Schmidt数（Sc）和化学反应参数（Kc）增大而下降，随Soret数（So）增大而上升。表面摩擦系数随Gr和Gm增大而增大，随M和β增大而减小；Nusselt数随Pr、Q、R增大而增大；Sherwood数随Sc增大而增大，随So增大而减小。该研究发表在《Results in Engineering》，为理解非牛顿流体在多孔介质中的耦合传输过程提供了理论依据，对材料加工、能源系统及地球物理应用具有重要意义。
**主要技术方法**
研究人员采用了有限差分法（FDM）对耦合的非线性偏微分方程组进行数值求解，利用MATLAB实现迭代计算。同时，应用响应面法（RSM）建立二次回归模型，以分析关键参数（磁场参数M、热源参数Q、Grashof数Gr）对表面摩擦系数、Nusselt数和Sherwood数的影响，并通过方差分析验证模型精度（决定系数R2>99.95%）。样本数据来源于FDM数值模拟结果，无特定临床队列来源。
**研究结果**
**1. Casson参数（β）对速度的影响**
图3显示，随着β增大，流体速度沿边界层持续下降。β增大意味着Casson流体塑性动态粘度μB增加，产生更大的流动阻力，导致速度降低。
**2. Grashof数（Gr）和修正Grashof数（Gm）对速度的影响**
图4和图5表明，Gr和Gm增大时，流体速度上升。表面冷却增强提高了Gr和Gm，增强了浮力效应，加速了流向多孔垂直板的流动。
**3. 磁场参数（M）对速度的影响**
图6显示，M增大使速度持续下降。磁场产生的Lorentz力阻碍流体运动，从而削弱动量边界层。
**4. Prandtl数（Pr）、热源参数（Q）和辐射参数（R）对温度的影响**
图7-9显示，Pr、Q和R增大时温度下降。Pr增大降低热扩散率，Q增大抽取更多能量，R增大增强辐射冷却，均导致热边界层变薄。
**5. Eckert数（Ec）对温度的影响**
图10显示，Ec增大使温度上升。更强的粘性耗散将动能转化为内能，提高了流体温度并增厚热边界层。
**6. Schmidt数（Sc）和化学反应参数（Kc）对浓度的影响**
图11和图12显示，Sc和Kc增大时浓度下降。Sc增大降低质量扩散率，Kc增大加速物种消耗，导致浓度边界层变薄。
**7. 渗透率参数（K）对速度的影响**
图13显示，K增大使速度逐渐上升。更高的渗透率降低多孔介质的阻力，促进流体运动。
**8. Soret数（So）对浓度的影响**
图14显示，So增大使浓度上升。Soret效应驱动质量流沿温度梯度方向移动，使更多溶质向高温区域聚集，增厚浓度边界层。
**9. 表面摩擦系数、Nusselt数和Sherwood数的变化**
表2展示了不同参数组合下的变化规律。表面摩擦系数随Gr、Gm、Kc增大而增大，随M、β增大而减小；Nusselt数随Pr、Q、R增大而增大；Sherwood数随Sc增大而增大，随So增大而减小。
**10. 响应面法（RSM）分析**
通过RSM建立了表面摩擦系数、Nusselt数和Sherwood数的二次回归模型。图15显示FDM数值结果与RSM预测值高度吻合，决定系数极高。图16和17的三维响应面和等高线图揭示了参数间的交互作用。图18的残差分析验证了模型的适用性，满足正态性和方差齐性假设。
**讨论与结论**
**讨论部分总结**
研究结果表明，各类参数对Casson流体流动、传热和传质的影响具有明确的物理机制。磁场和Casson参数通过增加阻力抑制流动，而浮力效应和渗透率则促进流动。热辐射和热源对温度分布起调节作用，化学反应和Soret效应分别降低和增强浓度。RSM成功建立了高精度回归模型，为快速预测工程响应提供了统计工具。
**结论部分翻译**
本研究数值探究了Casson流体通过嵌入多孔介质的垂直板的MHD流动。通过合适的无量纲参数将偏微分方程组简化，并利用有限差分法求解。以图形方式展示了不同参数下的流场分布。关键结论如下：
- 较高的渗透率（K）和Grashof数（Gr, Gm）通常能提高流体速度，因为它们促进浮力驱动流动并降低多孔介质阻力。
- 然而，随着Casson参数（β）和磁场参数（M）增大，速度下降，原因是屈服应力和Lorentz力限制流体运动。
- 随着热吸收参数（Q）和Prandtl数（Pr）增大，温度下降，因为增加的热吸收消耗更多能量，而高Pr导致热扩散减小，热边界层变薄。
- 随着化学反应参数（Kc）和Schmidt数（Sc）增大，浓度下降，归因于质量扩散率增大和溶质更快消耗。相反，浓度随Soret数（Sr）增大而升高，因为热扩散沿温度梯度增强质量传递。
- 表面摩擦系数随Kc、Gc和Gr增大而升高（壁面附近动量更强），但随M和β增大而降低（限制流动）。
- Sherwood数随Sc增大而升高，但随Sr降低（由于相反的质量和热扩散效应）；而Nusselt数随Q和Pr增大而升高，表明传热改善。
