含交叉扩散的生物对流Casson纳米流体流动中热辐射影响的三维磁流体动力学分析

《Results in Chemistry》：Impact of thermal radiation in bioconvection Casson nanofluid flow with cross-diffusion

　　

研究人员研究了Darcy多孔介质中拉伸薄板上的稳态三维磁流体动力学（MHD）生物对流Casson纳米流体流动。该模型综合了Hall电流修正的磁力、非线性Rosseland辐射、黏性耗散、Soret–Dufour交叉扩散、Arrhenius活化能以及速度、热、溶质和微生物滑移边界条件。通过带Newton迭代的自适应Runge–Kutta–Fehlberg（RKF45）打靶法求解相似变换后的边值问题。研究人员考察了关键参数对速度、温度、浓度和微生物分布的影响。结果表明：增大Casson参数会提高流体速度和温度；增强磁场会减慢流动，但因Joule热效应而使温度升高；渗透率参数增大使速度上升但温度下降；Prandtl数增大使温度降低。浓度随Schmidt数增大而减小，但随Soret效应或化学反应的增强而上升。最后，Lewis数和Peclet数增大均会降低流动中的微生物密度。主要定量趋势如下：当Prandtl数从0.7增至2时，局部Nusselt数从0.32464升至0.55131（约70%）；辐射参数从0.1增至2时，局部Nusselt数从0.42871降至0.19450（约55%）；此外，Schmidt数从0.1增至0.5时，Sherwood数从0.28890升至0.59234。研究结果阐明了流变学、电磁阻尼、辐射、交叉扩散与生物对流如何在热管理、微流控、反应加工和生物工程装置中被协同调控。

研究背景方面，非牛顿纳米流体在拉伸表面上的热特性表征是换热器、电子冷却、聚合物加工和生物医学装置等领域的重要研究方向。Casson流体具有屈服应力和剪切变稀特性，适合模拟血液流变、聚合物熔体和工业悬浮液。三维流动构型比二维理想化更接近实际工业热过程，而实际热应用往往涉及多方向空间变化，需要更精确的数学框架。同时，热辐射在高温或光学参与介质中显著影响热量传递，黏性耗散在高速度梯度流动中会产生额外热源，磁场效应可提供非接触流动控制，Soret–Dufour交叉扩散和Arrhenius活化能化学反应会进一步增强热场与浓度场的耦合，旋转微生物引起的生物对流则通过定向游泳产生密度梯度并驱动对流运动。然而，已有三维Casson流动研究通常只处理上述机制的有限子集；参考构型虽然包含三维辐射磁流体动力学Casson流动和化学反应，但未同时纳入Hall电流修正磁力、旋转微生物生物对流、Arrhenius活化能、Soret–Dufour耦合以及四种广义滑移条件。因此，该研究旨在统一的二维边界层模型中定量揭示这些耦合机制如何共同改变动量、热量、溶质和微生物的输运率，为先进热管理系统和生物工程装置提供设计依据。

研究人员构建了稳态三维磁流体动力学生物对流Casson纳米流体流过拉伸薄板的数学模型，该模型位于Darcy多孔介质中，并综合了Hall电流修正磁场力、非线性Rosseland辐射、黏性耗散、Soret–Dufour交叉扩散、Arrhenius活化能驱动的化学反应，以及速度、热、溶质和微生物四种滑移边界条件。通过相似变换将控制方程化为耦合非线性常微分方程组，再采用数值方法求解，并对速度、温度、浓度和微生物分布进行系统的单参数分析。研究得到的关键定量关系表明，流变参数、电磁阻尼、辐射、交叉扩散和生物对流可以被联合调节，从而用于热管理、微流控、反应加工和生物工程装置的性能优化。

在关键的技术方法方面，研究人员采用相似变换将原始偏微分方程组化为耦合非线性常微分方程组；使用四阶自适应Runge–Kutta–Fehlberg（RKF45）打靶法并结合Newton–Raphson迭代求解，迭代容差为ε=10^?6，计算域截断至η_∞=10。研究人员以二维无磁场、无Hall效应、无多孔介质、无滑移且无交叉扩散的极限情形验证数值代码，数值壁面剪切结果与解析解一致，并进行了远场残差、域截断和参数连续性检查。随后在基线参数基础上逐次改变单个参数，得到速度、温度、浓度和微生物分布及壁面传递量。本研究为理论数值研究，不含实验样本队列。

研究结果部分，4.1速度分布方面，通过单参数扫描发现，Casson参数β增大时，主速度f′(η)明显上升，原因是屈服应力效应减弱，流动更易维持；Hall修正磁参数M₁增大时，两个面内速度分量均下降，源于增强的Lorentz阻力；渗透率参数K增大时，Darcy阻力减小，速度总体增强。4.2温度分布方面，β增大使温度剖面下移，热边界层厚度减小；M₁增大使温度升高，归因于Joule加热；K增大使温度明显下降；Prandtl数Pr从0.71增至2.00时，在η=2处温度从约0.68降至约0.12，降幅约82%；辐射参数Rd增大使热边界层增厚、温度剖面上升，因为有效热扩散率增大；Eckert数Ec增大通过黏性耗散使温度升高；Dufour数Du增大使温度剖面下移，壁面温度梯度?θ′(0)增大，强化换热。4.3浓度分布方面，Schmidt数Sc增大使浓度边界层变薄、浓度下降；Soret数Sr增大使浓度剖面略微上移；微生物Lewis数Lb增大使微生物边界层急剧变薄，壁面微生物密度梯度增强；生物对流Peclet数Pe增大使微生物密度降低，因为更强平流将微生物带离表面；化学反应的增强也会使浓度上升。4.4拉伸速率比影响方面，λ=b/a改变y方向壁面拉伸强度，通过耦合速度场间接影响主速度f′(η)和温度分布，它不是滑移参数效应。4.5壁面传递量的定量结果方面，局部Nusselt数Nu_x随Pr增大而显著升高，随Rd增大而降低；Du增强Nu_x，Ec降低Nu_x。Sherwood数Sh_x随Sc增大而明显升高，随Sr增大而略有下降。微生物密度数Nn_x随Lb和Pe增大而增加。磁参数M从0增至1.5时，壁面摩擦系数绝对值|C_f,x|增加约19%。

讨论部分指出，无量纲参数对输运现象的影响具有层级性：动量输运主要受Casson参数、磁场和孔隙率控制；热输运主要由Pr和Rd控制，Ec在高值时变得关键；浓度输运主要取决于Sc；生物对流受Lb和Pe控制；Soret–Dufour交叉扩散和滑移条件属于次级但不可忽略的修正。这种层级性有助于设计者优先调控主导参数，同时将次级效应作为扰动处理。研究还强调辐射虽然增强体相热输运，但会削弱壁面热传递；黏性耗散可能造成热过冲；在高梯度多组分系统和生物反应器中，交叉扩散与生物对流的耦合效应不能被忽略。

最后，研究结论部分可概括为：该研究对拉伸薄板上的三维磁流体动力学生物对流Casson纳米流体流动进行了全面热分析，模型包含热辐射、黏性耗散、Soret–Dufour交叉扩散、Arrhenius活化能驱动化学反应、Hall电流和多滑移边界条件，并采用经过验证的RKF45打靶法进行数值求解。主要结论是：β增大使流体速度显著增强；M₁增大使两速度分量降低；K增大通过降低流动阻力使速度增强；温度随β、M₁、Ec和Rd的增大而升高，随K、Pr和Du的增大而降低；浓度随Sc增大而显著下降，随Sr和化学反应参数增大而明显上升；Lb和Pe增大会降低微生物密度。未来工作应考虑非定常强迫、变热物性、非Fourier热输运、纳米颗粒团聚以及实验或高保真数值验证。