磁场作用下Oldroyd-B流体的非稳态驻点流动

Unsteady Stagnation-Point Flow of Oldroyd-B Fluid along a Stretching Sheet with Magnetic Field

摘要: 本文研究了拉伸板上Oldroyd-B流体的非稳态驻点流动，此外，还考虑了磁场对流动的影响。利用合适的相似变换将控制方程转化为常微分方程。结合双参数变形展开方法与基函数方法(DPTEM-BF)求得微分方程的近似解析解。通过图形展示了各物性参数对速度场的影响，值得注意的是：随着松弛时间参数的增加，流体的粘性力增大使得流体流动阻力增强，这导致了流体流速减小；延迟时间参数的变化使得速度场增大，并且动量边界层厚度变厚；较大的磁场参数对应于较大的洛伦兹力，这个力阻碍了流体流动，因而流体流动速度降低。

Abstract: This paper examines the influence of magnetic field on unsteady stagnation-point flow of Oldroyd-B fluid towards a stretching sheet. Suitable similarity transformations are put into use to yield the ordinary differential equation, which are dealt with double-parameter transformation expansion method with base function method (DPTEM-BF). Impacts of various physical parameters on the velocity field are explored via graphs. It is noteworthy that as relaxation time parameter enlarges, the viscous force of fluid increases, which causes larger resistance to fluid flow. The increase of retardation time parameter increases both the velocity and momentum boundary layer thickness. Larger magnetic parameter corresponds to the larger Lorentz force, which impedes fluid flow and thus slows it down.

1. 引言

近年来，由于非牛顿流体在工业、食品以及日常生活中的广泛应用而备受研究人员的关注，如石油、油漆、淀粉液、蛋清、果浆等都是非牛顿流体。其中，Oldroyd-B流体就是常见的一种非牛顿流体，它不仅能够反映流体的松弛特性，还能表征延迟特性对流体的影响。Khan等 [1] 研究了Oldroyd-B流体流动问题。Hayat等 [2] 对多孔拉伸板上Oldroyd-B流体的混合自然流动进行了讨论。Ullah等 [3] 得到了两同轴圆柱之间Oldroyd-B流体旋转流动的精确解。此外，Alshomrani等 [4] 还探索了在分层机制作用下，热辐射和化学反应对Oldroyd-B纳米流体流动和传热的影响。Hayat等 [5] 在Alshomrani等人的基础之上还考虑了热产生/吸收的影响并且推导了其关系式。Hafeez等 [6] 采用Cattaneo-Christov双扩散模型分析了多孔旋转圆盘上Oldroyd-B流体的流动、传热和传质，并且得到了数值解。

驻点流动在熔融纺丝工艺、拉丝、玻璃吹制、金属和聚合物挤出等应用中起着至关重要的作用，因此不少学者都对驻点流动问题进行了研究。Mustafa等 [7] 研究了幂律流体的驻点流动。Sharma等 [8] 对具有二阶速度滑移的拉伸板上微极流体的驻点流动进行了分析。在布朗运动和热泳扩散存在下，Bai等 [9] 对拉伸板上Maxwell纳米磁流体的驻点流动也进行了讨论。在此基础之上，Bai等 [10] 还采用Cattaneo-Christov双扩散模型探索了具有化学反应和热辐射效应的上随体Oldroyd-B纳米磁流体的驻点流动和传热。上述学者都只研究了稳态驻点流动，但实际上流体的流动与时间也有关，以下研究人员对不同流体的非稳态驻点流动也进行了分析。Hayat等 [11] 讨论了拉伸板上二级流体的非稳态驻点流动问题。Zaib等 [12] 研究了速度滑移和热泳效应对微极流体非稳态驻点流动的影响。Awais等 [13] 还分析了在焦耳加热下，径向拉伸圆盘上Maxwell磁流体的非稳态驻点流动。Mahdy [14] 对具有Biot数效应的加热拉伸板上Oldroyd-B纳米流体的非稳态驻点流动进行了建模，并且还考虑了纳米流体的动态粘度和导热率的影响。此外，Hayat等 [15] 也对Oldroyd-B纳米流体的非稳态驻点流动问题进行了探索，并且将热产生/吸收、非线性热辐射、布朗运动以及热泳扩散的影响考虑在内。

磁流体由于同时具有独特的磁性以及流动性在扬声器、阻尼器和选矿分离等方面有着独特的价值，这一特性使得不少学者对存在磁场情况下的流体的流动问题进行了深入思考。Asghar等 [16] 得到了在磁场作用下粘性流体的库埃特流动的解析解。Hayat等 [17] 研究了倾斜磁场对指数拉伸板上三级流体流动的影响。Mabood等 [18] 对在磁场和加热作用下，拉伸板上Jeffery流体的非稳态流动问题进行了讨论，并且还考虑了变导热系数的影响。Jakati等 [19] 分析了在斜向外磁场作用下，非均匀热源/散热器、布朗运动和热泳对线性拉伸板板上Maxwell纳米流体流动的影响。此外，Hashmi等 [20] 还研究了磁场对两个无限等温拉伸盘间Oldroyd-B纳米流体的混合对流以及轴对称流动的影响。然而，很少有学者研究了磁场对Oldroyd-B流体非稳态驻点流动的影响。

本文主要研究了在磁场作用下，拉伸板上Oldroyd-B流体的非稳态驻点流动。采用合适的相似变换将偏微分方程转化为常微分方程，并且使用双参数变形展开方法与基函数方法结合的方法(DPTEM-BF)进行求解。最后通过图形分析了松弛时间参数、延迟时间参数、磁场参数等对速度场的影响。

2. 数学模型

考虑拉伸板上Oldroyd-B流体的非稳态驻点流动问题。建立如图1所示的流动示意图，其中，拉伸板的方向沿着x轴的方向，y轴垂直于拉伸板的方向。假设拉伸板的速度为 $U_w=bx/\left(\text{1}-at\right)$ ，驻点流动速度为 $U_e=cx/\left(\text{1}-at\right)$ ，其中 $a,b,c>0$ 是常数。此外，沿着y轴方向施加了一个强度为B₀的均匀磁场。

Figure 1. The physical flow diagram

图1. 流动示意图

Oldroyd-B流体的非稳态驻点流动控制方程为：

$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$ (1)

$\begin{array}{c}\frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=-{\lambda }_1\left(\frac{{\partial }^{2}u}{\partial t^2}+2u\frac{{\partial }^{2}u}{\partial x\partial t}+2v\frac{{\partial }^{2}u}{\partial y\partial t}+u^2\frac{{\partial }^{2}u}{\partial x^2}+2uv\frac{{\partial }^{2}u}{\partial x\partial y}+v^2\frac{{\partial }^{2}u}{\partial y^2}\right)\\ +\nu \frac{{\partial }^{2}u}{\partial y^2}+{\lambda }_2\nu \left(\frac{{\partial }^{3}u}{\partial y^2\partial t}+v\frac{{\partial }^{3}u}{\partial y^3}+u\frac{{\partial }^{3}u}{\partial x\partial y^2}-\frac{\partial u}{\partial x}\frac{{\partial }^{2}u}{\partial y^2}-\frac{\partial u}{\partial y}\frac{{\partial }^{2}v}{\partial y^2}\right)\\ \text{+}\frac{\partial u_e}{\partial t}+u_e\frac{\partial u_e}{\partial x}+{\lambda }_1{u}_e^2\frac{{\partial }^2{u}_e}{\partial x^2}-\frac{\sigma B_0^2}{\rho }\left[u-u_e+{\lambda }_1\left(\frac{\partial u}{\partial t}+v\frac{\partial u}{\partial y}-\frac{\partial u_e}{\partial t}\right)\right]\end{array}$ (2)

边界条件为：

当 $y=0$ 时， $u=u_w\left(x,t\right)=\frac{bx}{1-at}$ ， $v=0$ (3)

当 $y\to \infty$ 时， $u=u_e\left(x,t\right)=\frac{cx}{1-at}$ (4)

其中，u和v分别为沿x和y轴方向上的速度分量， $\nu =\mu /\rho$ 指的是流体的运动粘度系数，λ₁和λ₂分别是松弛时间参数和延迟时间参数。

采用如下的无量纲函数以及相似变量：

$u=\frac{bx}{\left(1-at\right)}{f}^{\prime }\left(\eta \right),v=-\sqrt{\frac{\nu b}{\left(1-at\right)}}f\left(\eta \right),\eta =\sqrt{\frac{b}{\nu \left(1-at\right)}}y,\psi =\sqrt{\frac{\nu b}{\left(1-at\right)}}xf\left(\eta \right)$ (5)

将(5)带入(2)~(4)，得到对应的常微分方程及边界条件如下：

$\begin{array}{l}S\left(f^{\prime }+\frac{1}{2}{f}^{″}\eta \right)+{\left(f^{\prime }\right)}^2-f{f}^{″}+{\beta }_1\{S^2\left(2f^{\prime }+\frac{7}{4}{f}^{″}\eta +\frac{1}{4}{f}^{‴}{\eta }^2\right)+2S\left[{\left(f^{\prime }\right)}^2+\frac{1}{2}{f}^{\prime }{f}^{″}\eta \right]\\ -S\left(3f{f}^{″}+f{f}^{‴}\eta \right)-2f{f}^{\prime }{f}^{″}+f^2{f}^{‴}\}-f^{‴}-{\beta }_2\left[S\left(2f^{‴}+\frac{1}{2}{f}^{iv}\eta \right)-f{f}^{iv}+{\left(f^{″}\right)}^2\right]\\ -AS-A^2+M\left\{\left(f^{\prime }-A\right)+{\beta }_1\left[S\left(f^{\prime }+\frac{1}{2}{f}^{″}\eta \right)-f{f}^{″}-SA\right]\right\}\end{array}$ (6)

对应的边界条件转化为：

当 $\eta =0$ 时， $f\left(0\right)=0$ ， $f^{\prime }\left(0\right)=1$ (7)

当 $\eta \to \infty$ 时， $f^{\prime }\left(\infty \right)=A$ (8)

其中， ${\beta }_{\text{1}}={\lambda }_{\text{1}}b/\left(\text{1}-at\right)$ 和 ${\beta }_{\text{2}}={\lambda }_{\text{2}}b/\left(\text{1}-at\right)$ 分别为松弛时间参数和延迟时间参数， $S=a/b$ 为不稳定参数， $A=c/b$ 是驻点参数， $M=\sigma B_0^2\left(\text{1}-at\right)/\rho b$ 指的是磁场参数。

3. DPTEM-BF方法求解

采用双参数变形展开方法与基函数结合的方法求解微分方程，根据双参数变形展开方法的原理，对自变量和因变量进行变形，引入小参数ε和新变量ξ，将函数f(η)转化为：

$f\left(\eta \right)={\epsilon }^{4}F\left(\xi \right)+f\left(0\right)+f^{\prime }\left(0\right)\eta +\frac{f^{″}\left(0\right)}{2!}{\eta }^2+\frac{f^{‴}\left(0\right)}{3!}{\eta }^3$ ， $\xi ={\epsilon }^{-1}\eta$ (9)

令 $f^{″}\left(0\right)=a_1$ ， $f^{‴}\left(0\right)=a_2$ ，其中a₁，a₂均为未知常数，F(ξ)为变量ξ的函数。

对应地，函数F(ξ)的初值条件为：

$F\left(0\right)=0$ ， $F^{\prime }\left(0\right)=0$ ， $F^{″}\left(0\right)=0$ ， $F^{‴}\left(0\right)=0$ (10)

将F(ξ)在 $\epsilon =0$ 附近展开成级数形式：

$F\left(\xi \right)=F_0\left(\xi \right)+{\displaystyle \underset{i=1}{\overset{\infty }{\sum }}{F}_i}\left(\xi \right){\epsilon }^i$ (11)

将(10)带入(11)中，得到：

$F_i\left(0\right)=0$ ， ${F^{\prime }}_i\left(0\right)=0$ ， $F_i{}^{\prime \text{}\prime }\left(0\right)=0$ ， $F_i{}^{\prime \text{}\prime \text{}\prime }\left(0\right)=0$ ， $i=0,\text{1},\text{2},\text{3,}\cdots$ (12)

将(9)和(12)带入式(6)中，把微分方程转化为关于ε的多项式，根据多项式恒为零的性质，使多项式中ε同次幂的系数为零，再利用初值条件，便可以计算出F_i(ξ)的表达式，紧接着便可以得到F(ξ)的表达式，消去小参数ε，得到含有a₁，a₂的关于f(η)的表达式。采用结合基函数来确定参数a₁，a₂的值，将函数f(η)表示成如下基函数的组合：

$f\left(\eta \right)\approx f_{N_1,N_2}\left(\eta \right)=f_{0,0}\left(\eta \right)+{\displaystyle \underset{j=3}{\overset{N_1}{\sum }}{\displaystyle \underset{i=1}{\overset{N_2}{\sum }}{a}_{ij}{f}_{i,j}\left(\eta \right)}}=f_{0,0}\left(\eta \right)+{\displaystyle \underset{j=3}{\overset{N_1}{\sum }}{\displaystyle \underset{i=1}{\overset{N_2}{\sum }}{a}_{ij}{\eta }^j{\text{e}}^{i{a}_0\eta }}}$ (13)

其中，f_0,0(η)为：

$f_{0,0}\left(\eta \right)=m_1+s\eta -m_1{\text{e}}^{a_0\eta }+\left(1-s+a_0{m}_1\right)\eta {\text{e}}^{a_0\eta }+m_2{\eta }^2{\text{e}}^{a_0\eta }$ (14)

令 $N_{\text{1}}=N_{\text{2}}=\text{3}$ ，其中m₁，m₂，a₀，a₁₃，a₂₃，a₃₃为假设的未知数，最后使用牛顿迭代法便可以计算出每个未知参数的值，将其带入函数f(η)的表达式，就可以得到原微分方程的近似解析解。

4. 结果与讨论

采用双参数变形展开方法与基函数结合的方法求得了微分方程的近似解析解，由表1可知，将所求结果与已有的研究结果进行比较，具有很好的一致性。此外，通过图形详细分析讨论了各物性参数对速度场的影响。

Table 1. Comparison of − f ″ ( 0 ) with different values of S while β 1 = β 2 = A = M = 0

表1. 当 ${\beta }_1={\beta }_2=A=M=0$ 时，S取不同值， $-f^{″}\left(0\right)$ 值的比较结果

图2展示了松弛时间参数β₁对速度场的影响，由图2可知，当β₁增加时，流体的粘性力增大，这使得流体的流动阻力增大，因此流体的流动速度减小。

图3绘制了不同延迟参数值β₂下的速度分布情况，从图中可以发现，延迟时间参数对速度场有着相反的影响，随着β₂的增大，流体动量穿透深度加深，拉伸板上流体的速度分量增加，这意味着流体流速增加。此外，β₂的增大也会使得流体动量边界层厚度变厚。

图4反映了不稳定参数S变化时速度场的变化趋势，由于S的增加，流体的流动速度以及边界层厚度都均匀地减小。

图5给出了驻点参数A对速度的影响，当A增加时，流体的动量损失减小，边界层内流体速度越快趋于周围流体速度，因此速度边界层厚度会变薄。

图6显示了磁场参数M变化时对速度的影响，随着M值的增大，流体的流动速度减慢，并且边界层厚度变薄。从物理上来讲，M的增加使得磁场作用下的洛伦兹力增大，这个力阻碍了流体流动，因此流体流动速度减小。

Figure 2. Influence of β₁ on the velocity field $f^{\prime }(\; \eta \; )$

图2. β₁对速度的影响

Figure 3. Influence of β₂ on the velocity field $f^{\prime }(\; \eta \; )$

图3. β₂对速度的影响

Figure 4. Influence of S on the velocity field $f^{\prime }(\; \eta \; )$

图4. S对速度的影响

Figure 5. Influence of A on the velocity field $f^{\prime }(\; \eta \; )$

图5. A对速度的影响

Figure 6. Influence of M on the velocity field $f^{\prime }(\; \eta \; )$

图6. M对速度的影响

5. 总结

本文研究了磁场对拉伸板上Oldroyd-B流体非稳态驻点流动的影响。通过合适的相似变换将控制方程转化为常微分方程，并且使用双参数变形展开方法与基函数方法结合的方法(DPTEM-BF)进行求解。采用图形详细分析了各物性参数对速度的影响，值得注意的是：随着驻点参数的增大，流体的动量损失减小，边界层内流体速度越快趋于周围速度，所以速度边界层厚度变薄；而磁场参数的增加会使得磁场作用下的洛伦兹力增大，这个力阻碍了流体流动，因此流体流动变缓。

基金项目

国家重点研究计划(2016YFC0700601)，国家自然科学基金(21878018)，北京市自然科学基金和北京市教育委员会联合资助项目(KZ201810016018)，北京建筑大学市属高校基本科研业务费项目(X20142)。