非线性常微分方程( 9 )和( 10 )根据边界条件(11),形成两点边值问题(BVP),并使用Cebeci和Bradshaw在书中描述的Keller box方法进行数值求解[ 29 ]. 在该方法中,通过以下四个步骤获得解决方案:
(i)
减少方程式( 9 )和( 10 )一级系统。
(ii)
用中心差分写出差分方程。
(iii)
用牛顿法将得到的代数方程线性化,并以矩阵-向量形式书写。
(iv)
利用块三元消元技术求解线性系统。
然后使用MATLAB R2010a软件对数值方法进行编程。 为了获得数值解,需要对步长进行适当的猜测 η , Δ η 和边界层的厚度
η
∞
(通常选择4到10之间的有限数)。 从的一些初始猜测值开始
η
∞
,方程式( 9 )和( 10 )根据边界条件(11)和一些特定的参数集进行求解,以获得速度剖面
(f)
′
( η ) 和感应磁场分布
克
′
( η ) 重复求解过程,直到中的进一步变化(增量)
η
∞
不会导致
(f)
″
( 0 ) 和
克
″
( 0 ) 或者,换句话说,结果独立于
η
∞
。采用的初始步长为 小时 = Δ η = 0.01 .表面摩擦系数
(f)
″
( 0 ) 、速度剖面
(f)
′
( η ) ,感应磁场分布
克
′
( η ) 以及感应磁场的变化率,我们以后称之为感应磁场梯度,
克
″
( 0 ) 得到了各种控制参数值,即运动参数 λ ,楔形参数 β 和磁性参数 S公司 为了评估所用数值方法的准确性,我们对非磁性情况下的一些结果进行了比较( S公司 = 0 )莱利和魏德曼获得的[ 14 ],拉贾戈帕尔 等。 [ 30 ],伊沙克 等。 [ 18 ]和郭[ 31 ]. 表 1 表示表面摩擦系数的值
(f)
″
( 0 ) 对于 λ = 0 以及楔形参数的各种值 β 对于非磁性外壳( S公司 = 0 ). 表 2 比较表面摩擦值
(f)
″
( 0 ) 对于三重解集,当 λ = 负极 0.4 , β = 0.05 莱利和魏德曼获得的[ 14 ]. 我们发现,本研究获得的结果与早期研究人员获得的结果非常一致。 因此,开发的代码可以放心地用于磁性外壳( 0 < S公司 < 1 ).
表1 的值
(f)
″
( 0 )
对于
λ = 0
,
S公司 = 0
和各种
β
表2 的值
(f)
″
( 0 )
对于
λ = 负极 0.4
,
β = 0.05
速度剖面的变化
(f)
′
( η ) 和感应磁场分布
克
′
( η ) 带有移动参数 λ ,楔形参数 β ,磁性参数 S公司 和倒数磁普朗特数 α 如图所示 1 到 4 所有样本剖面均渐近满足远场边界条件(11),从而支持所获得的数值结果。 从这些数字中,我们可以看出 λ 和 β 增加流体速度
(f)
′
( η ) 和感应磁场
克
′
( η ) 速度边界层厚度减小时也会增加。 相比之下,作为 S公司 和 α 随着速度边界层厚度的增加,流体速度和感应磁场减小。 我们还注意到倒数磁普朗特尔数的影响 α 在上更加明显
克
′
( η ) 与相比
(f)
′
( η ) .
图 5 表示表面摩擦系数的变化
(f)
″
( 0 ) 作为的函数 λ 对于各种值 S公司 当楔形参数和倒数磁普朗特数固定在 β = 0.5 和 α = 1 分别为。 发现对于所有磁性参数值 S公司 ,使用 0 ≤ S公司 < 1 ,该解决方案对于的所有值都是唯一的 λ ≥
λ
c(c)
,其中
λ
c(c)
是的最小值 λ 存在解决方案的。 临界值 |
λ
c(c)
| 随着 S公司 增加。 在我们对案件的计算中 β = 0.5 ,当感应磁梯度值达到时,溶液停止存在
克
″
( 0 ) = 0 如图所示 5 ,表面摩擦系数的值
(f)
″
( 0 ) 也会随着磁性参数的值而减小 S公司 增加。 此外,当 S公司 .
图 6 表示表面摩擦系数的变化
(f)
″
( 0 ) 和感应磁场梯度
克
″
( 0 ) 作为的函数 λ 对于各种磁参数值 S公司 当楔形参数固定在 β = 0.7 。该图表明,对于 S公司 ( 0 ≤ S公司 < 1 ),该解决方案对于的所有值都是唯一的 λ ≥ 0 ,虽然对于某些范围的值存在双重解决方案
λ
c(c)
≤ λ < 0 此外,皮肤摩擦系数
(f)
″
( 0 ) ,感应磁梯度
克
″
( 0 ) 和临界值 |
λ
c(c)
| 随着 S公司 增加。 如表所示 三 ,临界值
λ
c(c)
对于 β = 0 和 β = 1 在非磁性情况下( S公司 = 0 )我们在本研究中计算的结果与Klemp和Acrivos之前报告的结果非常一致[ 32 ]和侯赛尼 等。 [ 33 ]. 案例 β = 0 对应于平板,而 β = 1 指驻点流。
表3 的值
λ
c(c)
对于不同的值
米
什么时候
S公司 = 0
图 7 显示了速度剖面
(f)
′
( η ) 在临界值 λ ( =
λ
c(c)
)在分离之前,对于不同的值 S公司 什么时候 β = 0.7 我们观察到 S公司 增加,临界值 |
λ
c(c)
| ,感应磁梯度
克
″
( 0 ) 和表面摩擦
(f)
″
( 0 ) 减少,从而支持我们之前从图中观察到的结果 6 .我们还注意到感应磁场梯度
克
″
( 0 ) 随运动参数几乎线性变化 λ ,第二种解决方案的值大多很小 |
克
″
( 0 ) | <
10
负极
4
.图 8 显示速度
(f)
′
( η ) 和感应磁场
克
′
( η ) 当 β = 0.7 , S公司 = 0.3 和 λ = 负极 0.85 .
图 9 说明了表面摩擦系数的变化
(f)
″
( 0 ) 和感应磁场梯度
克
″
( 0 ) 作为的函数 λ 对于各种磁参数值 S公司 当楔形参数固定在 β = 0.03 。该图表明,对于所有值 S公司 ( 0 ≤ S公司 < 1 ),该解决方案对于的所有值都是唯一的 λ ≥ 0 ,虽然对于某些范围的值存在三重解决方案
λ
c(c)
≤ λ < 0 与之前考虑的情况类似,表面摩擦系数
(f)
″
( 0 ) ,感应磁梯度
克
″
( 0 ) 和临界值 |
λ
c(c)
| 也会随着 S公司 增加。 这里,我们还发现,对于第二和第三个解,感应磁梯度的值通常也很小 |
克
″
( 0 ) <
10
负极
4
| .
图 10 表示表面摩擦系数的变化
(f)
″
( 0 ) 作为的函数 λ 对于楔形参数的各种值 β 当磁性参数固定在 S公司 = 0.3 此处,倒数磁普朗特尔数的值也固定为 α = 1 结果显示出与非磁性情况下获得的结果相似的特性( S公司 = 0 )莱利和魏德曼报道[ 14 ]. 我们可以看到
(f)
″
( 0 ) 增加为 β 增加,并且有一个临界值
λ
c(c)
超过该移动参数不存在相似性解。 临界值的大小 |
λ
c(c)
| 也会随着楔形参数的增加而增加 β 增加。 图 7 还表示根据楔形参数的值有多种解决方案 β .跟随莱利和魏德曼[ 14 ]对于非磁性外壳( S公司 = 0 ),我们特别注意在存在磁场的情况下解集的以下有趣特征 S公司 = 0.3 。对于 0.6 < β ≤ 1 ,有一个独特的解决方案 λ ≥ 0 和一些范围的双重解决方案
λ
c(c)
< λ < 0 ; 对于 0.1 < β ≤ 0.6 ,解决方案对所有人来说都是独一无二的 λ ≥
λ
c(c)
; 对于 0 ≤ β ≤ 0.1 ,对于某些参数值范围,有三种解决方案 λ 更具体地说,我们的计算表明 β = 0.03 ,为该范围找到了唯一的解决方案 负极 0.8365 ≤ λ ≤ 负极 0.2502 和 λ ≥ 负极 0.336 ,虽然已经为该范围找到了三种解决方案 负极 0.335 ≤ λ ≤ 负极 0.2503 该结果与莱利和魏德曼报告的结果在质量上一致[ 13 ],其中找到了三种解决方案 0 < β < 0.14 ,为所有人提供独特的解决方案 λ 什么时候 0.14 < β < 0.5 和一些范围的双重解决方案 λ 什么时候 0.5 < β < 1 我们在这里提到莱利和魏德曼[ 14 ]据报道,对于非磁性外壳( S公司 = 0 ),的所有解决方案曲线 β > 0 有道理 ( 1 , 0 ) 作为限制点。 在本研究中,我们的计算表明,解曲线在点附近终止 ( 负极 0.8365 , 0 ) 如图所示 7 莱利和魏德曼[ 14 ]从“边缘”的角度解释了这个极限点的重要性
η
电子
边界层的厚度。 此外
η
电子
增加为
(f)
″
( 0 ) 减小,直到接近极限点,
η
电子
= ∞ 。如果参数的值为 S公司 增加。 数字 11 , 12 , 13 支持三重解存在的速度剖面和感应磁场剖面的现有样本 β = 0 , β = 0.03 和 β = 0.1 分别为。
图 14 显示了感应磁场梯度的变化
克
″
( 0 ) 作为运动参数的函数 λ 使用楔形参数 β .我们观察到感应磁场梯度
克
″
( 0 ) 随着楔形参数增加 β 较小值的增加 λ 但与 β 对于更大的值 λ .
按照早期研究人员采用的惯例,我们将解的前两个上层分支定义为
(f)
″
( 0 ) 对于给定值为 β ,而第三个分支的值最小
(f)
″
( 0 ) 我们注意到速度剖面
(f)
′
( η ) 对于解的前两个上分支,表现出相同的单调行为。 第一分支的边界层通常很薄,速度剖面
(f)
′
( η ) 迅速实现价值
(f)
′
( ∞ ) = 1 一般来说,与其他两个分支相比,第三分支的解通常涉及更大的边界层厚度。 它的特点通常是从一个相当小的值开始
(f)
″
( 0 ) > 0 速度剖面的发展具有非单调性
(f)
′
( η ) ,在假定其最终渐近值之前
(f)
′
( ∞ ) = 1 莱利和魏德曼报告了类似的非单调行为[ 14 ]当他们考虑速度剖面时
(f)
′
( η ) 的上分支解决方案 负极 1 ≤ β < 负极 0.5 .跟随Ishak 等。 [ 16 ],我们假设解的上分支具有最高值
(f)
″
( 0 ) (第一种解决方案)在物理上是稳定的,并且在实践中发生,因为它是 λ > 0 , 即。 当流体和固体表面沿同一方向移动时。
减少表面摩擦
(f)
″
( 0 ) 意味着阻力减小。 因此,磁场减少了阻力,加快了分离速度。 另一方面,增加楔块的夹角将增加阻力,从而延迟分离。 该结果与Ishak报告的结果一致 等。 [ 18 ]。
根据洛伦兹定律,感应磁场将反对原始磁场的变化,而不是磁场本身。 例如,如果原始场正在减小,那么感应磁场必须与原始场在同一方向上才能阻止减小。 从数字 6 和 9 ,我们看到感应磁场梯度
克
″
( 0 ) 随着值的增加而单调增加 λ 这种增加被认为是与原始磁场的减少相对立的。 此外,感应磁场梯度也随着 S公司 这符合洛伦兹定律。 我们还注意到 S公司 和 β 表面摩擦更明显
(f)
″
( 0 ) 与感应磁梯度相比
克
″
( 0 ) .