International Journal of Fluid Dynamics
Vol.04 No.04(2016), Article ID:19244,13 pages
10.12677/IJFD.2016.44009

Effects of Variable Thermal Conductivity and Viscosity on MHD Stagnation-Point Flow

Shengnan Li

College of Mathematics and Systems Science, Shenyang Normal University, Shenyang Liaoning

Received: Nov. 25th, 2016; accepted: Dec. 9th, 2016; published: Dec. 15th, 2016

Copyright © 2016 by author and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

ABSTRACT

The two-dimensional MHD stagnation-point flow with variable thermal conductivity and viscosity is studied. The thermal conductivity and viscosity are considered as functions of temperature. The governing Navier-Stokes equations are transformed into a set of ordinary differential equations (ODEs) by similarity transformation. The transformed ODEs are solved numerically using shooting method. Numerical calculations for various magnetic parameters, fluid viscosity parameters and Prandtl numbers are carried out, and the effects of magnetic strength and temperature changes to the skin friction coefficient, heat transfer near the wall and flow field characteristics are discussed in detail.

Keywords:Stagnation-Point Flow, MHD Flow, Similarity Solutions, Variable Thermal Conductivity, Variable Viscosity

考虑变热传导系数和粘性系数的MHD驻点流 问题研究

李胜男

沈阳师范大学数学与系统科学学院,辽宁 沈阳

收稿日期:2016年11月25日;录用日期:2016年12月9日;发布日期:2016年12月15日

摘 要

本文研究了考虑变热传导系数和粘性系数的二维定常MHD驻点流动问题。其中热传导系数和粘性系数是温度的函数。通过相似变换,将N-S控制方程组转化为常微分方程组;再利用打靶法计算此常微分方程组的数值解。分别对不同的磁场参数、流体粘性参数和普朗特数进行了数值计算,详细分析了磁场效应和温度变化对壁面摩擦系数、壁面传热及流场特征等的影响。

关键词 :驻点流,MHD流,相似解,变热传导系数,变粘性系数

1. 引言

二维驻点流问题是流体力学中的一个经典问题,并且被广泛应用于工业应用领域。例如,均匀磁场下导电流体的驻点流问题在现代冶金和金属加工过程中起着非常重要的作用,因为其最终产品的质量往往取决于金属冷却速率,而这个速率又与壁面边界层结构有关。

Hiemenz [1] 最先研究了经典二维平板驻点流问题,并且得到了N-S控制方程的精确解。从此,许多研究者对这个问题进行了不同方面的研究,例如三维驻点流问题、非定常驻点流问题、斜驻点流问题及带有热边界层的驻点流问题。Mahapatra [2] 等对平板二维定常不可压粘性导电流体的MHD边界层方程给出了精确的相似解,研究了常温表面上的流体温度分布情况。Grosan [3] 研究了固定平板MHD斜入射驻点流问题。在他的研究中,考虑了大小为常数的横向磁场。Lok [4] 等研究了拉伸/收缩平板上粘性导电流体的MHD斜入射驻点流问题,其中外加磁场是均匀的,方向平行于斜入射驻点流的分界流线。Mahapatra [5] 等考虑了收缩平板上斜入射驻点流的传热和辐射问题,研究发现流体内某一点的温度随着有效普朗特数的增加而减小。Chiam [6] 表示出了拉伸表面驻点流边界层传热问题的相似解,并讨论了其对流场的影响,其中热传导系数随着温度线性变化。Ali [7] 等扩展了此问题,考虑了感应磁场对流场的影响,研究了正交于拉伸平板的粘性不可压驻点流问题和常热传导系数的传热问题。Singh [8] 等研究了固定横向磁场下MHD斜入射驻点流问题。

本文研究变粘性系数二维粘性不可压驻点流问题,其中热传导系数和粘性系数都是温度的函数,同时考虑感应磁场对整个流场的影响。

2. 基本方程组

本文考虑二维定常粘性不可压流体在正交于固定平板y = 0平面上流动,如图1所示。导电流体放置在外加磁场H中,同时考虑感应磁场的影响。根据Davies [9] ,此问题的控制方程可以写为:

(1)

(2)

(3)

Figure 1. Physical model and coordinate system

图1. 物理模型及坐标系

(4)

其中:x和y是沿平板表面及与其垂直平面建立的笛卡尔坐标;是沿x和y轴的速度分量;是沿x和y轴的磁场强度分量;是流体密度;是磁导率;是电导率;是常压下的比热;是边界层边界的x轴方向速度分量和磁场强度分量。

本文中,粘性系数为T的函数,定义为

(5)

(6)

其中,都是常数,大小取决于的参比状态及小参数,反映流体的热性能。

已经发现,对于0˚F至400˚F的液体金属,热传导系数是关于T的线性函数。我们假设热传导系数有以下形式 [10]

(7)

其中,是一个小量级参数,分别是壁面上及远离壁面处的热传导系数,是壁面上及远离壁面处的流体温度。

边界条件如下

, c (8)

(9)

其中,是大于零的常数,是外流场均匀磁场的大小。

引入下面的相似变换

(10)

(11)

(12)

非线性偏微分方程组(1)~(4)可转化为如下常微分方程组

(13)

(14)

(15)

边界条件(8)~(9)转化为

(16)

(17)

(18)

其中:参数分别表示磁普朗特数的倒数,普朗特数,磁参数和流体粘性参数。其定义为

(19)

(20)

(21)

(22)

当流体粘性系数为常数时,方程(13)退化为

(23)

该结果与Ali et al. [7] 方程的结果一致。

当热传导系数为常数时,方程(15)退化为

(24)

该结果与Mahapatra et al. [5] 方程的结果一致。

我们关注的量为表面摩擦系数和努赛尔特数,分别定义为

(25)

(26)

其中,壁面剪切应力,壁面热流量为当地雷诺数。

3. 结果与讨论

我们采用龙格–库塔方法和打靶法对方程(13)~(15)及其边界条件(16)~(18)进行数值求解。为了检验所使用数值求解方法的正确性和准确性,我们计算了粘性系数取常数且不考虑磁场,即的值。数值结果与Howarth [11] 和Li et al. [12] 的结果进行了对比,如表1所示。可以发现,计算结果与之前发表过的结果吻合良好,说明我们的数值方法是准确的。

首先计算了表面摩擦系数随磁系数和流体粘性参数的变化。如图2所示,表面摩擦系数随着的增加而减小,也随着的增加而减小。图3给出了当地努赛尔特数-数变化情况,可以发现,数的增加而增加,随着的增加而减小。

图4~6表示,当时,磁系数对水平速度剖面、感应磁场剖面、温度剖面的影响。可以看出,速度和感应磁场随着的增加而减小,而温度随着的增加而增加。

图7~9表示,当时,粘性系数对水平速度剖面、感应磁场剖面和温度剖面的影响。可以看出,速度和感应磁场随着的增加而减小,而温度随着的增加而增加。

图10~12表示,当时,普朗特数对水平速度剖面、感应磁场剖面和温度剖面的影响,可以看出,速度、感应磁场和温度都随着的增加而减小。从图12

Table 1. Value of when = const and

表1.的数值(为常数,)

Figure 2. Variation of with for different when and

图2.时,不同值对应的的变化

Figure 3. Variation of with for different when and

图3.时,不同值对应的的变化

Figure 4. Velocity profiles for different values of when, and

图4.时,不同值对应的剖面

Figure 5. profiles for different values of when, and

图5.时,不同值对应的剖面

Figure 6. Temperature profiles for different values of when, and

图6.时,不同值对应的剖面

Figure 7. Velocity profiles for different values of when, and

图7.时,不同值对应的剖面

Figure 8. profiles for different values of when, and

图8.时,不同值对应的剖面

Figure 9. Temperature profiles for different values of when, and

图9.时,不同值对应的剖面

Figure 10. Velocity profiles for different values when, and

图10.时,不同值对应的剖面

Figure 11. profiles for different values of when, and

图11.时,不同值对应的剖面

Figure 12. Temperature profiles for different values of when, and

图12.时,不同值对应的剖面

Figure 13. Velocity profiles for different values of when, and

图13.时,不同值对应的剖面

Figure 14. profiles for different values of when, and

图14.时,不同值对应的剖面

Figure 15. Temperature profiles for different values of when, and

图15.时,不同值对应的剖面

可以发现,热边界层厚度随着的增加而减少。这个结果是我们所期望得到的,因为传热能力会随着的增加而减小,即热边界层厚度将会减少。

图13~15表示,当时,对水平速度剖面、感应磁场剖面和温度剖面的影响,可以看出,速度、感应磁场和温度都随着的增加而增加。从图13~14中也可以发现对水平速度剖面和感应磁场剖面的影响很微弱,这是因为没有出现在方程(13)~(14)中,导致没有直接通过温度影响

4. 结论

本文研究了考虑感应磁场的二维粘性导电流体驻点流问题,详细分析了变热传导系数和变粘性系数下流场随磁系数、粘性参数及普朗特数的变化情况。主要结果总结如下:速度和感应磁场随着的增加而减小,而温度随着的增加而增加。速度、感应磁场和温度都随着的增加而减小,随着的增加而增加。对速度剖面和感应磁场剖面的影响较对温度剖面的影响弱。

基金项目

辽宁省自然科学基金资助项目(2015020650)。

文章引用

李胜男. 考虑变热传导系数和粘性系数的MHD驻点流问题研究
Effects of Variable Thermal Conductivity and Viscosity on MHD Stagnation-Point Flow[J]. 流体动力学, 2016, 04(04): 69-81. http://dx.doi.org/10.12677/IJFD.2016.44009

参考文献 (References)

  1. 1. Hiemenz, K. (1991) Die Grenzschicht an einem in den gleichförmigen Flüssigkeitsstrom eingetauchten geraden Kreiszylinder. Dingler’s Polytech., 3326, 321-324.

  2. 2. Mahapatra, T.R. and Gupta, A.S. (2001) Magnetohydrodynamic Stagnation-Point Flow towards a Stretching Sheet. Acta Mechanica, 152, 191-196. https://doi.org/10.1007/BF01176953

  3. 3. Grosan, T., Pop, I., Revnic, C. and Ingham, D.B. (2009) Magnetohydrodynamic Oblique Stagnation-Point Flow. Meccanica, 44, 565-572. https://doi.org/10.1007/s11012-009-9196-0

  4. 4. Lok, Y.Y., Merkin, J.H. and Pop, I. (2015) MHD Oblique Stagnation-Point Flow towards a Stretching/Shrinking Surface. Meccanica, 50, 1-13. https://doi.org/10.1007/s11012-015-0188-y

  5. 5. Mahapatra, T.R., Nandy, S.K. and Gupta, A.S. (2012) Oblique Stagnation-Point Flow and Heat Transfer towards a Shrinking Sheet with Thermal Radiation. Meccanica, 47, 1325-1335. https://doi.org/10.1007/s11012-011-9516-z

  6. 6. Chiam, T.C. (1996) Heat Transfer with Variable Conductivity in a Stagnation-Point Flow towards a Stretching Sheet. International Communications in Heat and Mass Transfer, 23, 239-248.

  7. 7. Ali, F.M., Nazar, R., Arifin, N.M. and Pop, I. (2011) MHD Stagnation-Point Flow and Heat Transfer towards Stretching Sheet with Induced Magnetic Field. Applied Mathematics and Mechanics (English Edition), 4, 409-418. https://doi.org/10.1007/s10483-011-1426-6

  8. 8. Singh, P., Tomer, N.S., Kumar, S. and Sinha, D. (2010) MHD Oblique Stagnation-Point Flow towards a Stretching Sheet with Heat Transfer. Journal of Applied Mathematics and Mechanics, 6, 94-111.

  9. 9. Davies, T.V. (1963) The Magneto-Hydrodynamic Boundary Layer in the Two-Dimensional Steady Flow Past a Semi- Infinite Flat Plate I, Uniform Conditions at Infinity. Proceedings of the Royal Society A, 273, 496-508. https://doi.org/10.1098/rspa.1963.0105

  10. 10. Prasad, K.V., Vajravelu, K. and Datti, P.S. (2010) The Effects of Variable Fluid Properties on the Hydro-Magnetic Flow and Heat Transfer over a Non-Linearly Stretching Sheet. International Journal of Thermal Sciences, 49, 603-610. https://doi.org/10.1016/j.ijthermalsci.2009.08.005

  11. 11. Howarth, L. (1934) On Calculation of the Steady Flow in the Boundary Layer near the Surface of a Cylinder in a Stream. Aeronautical Research Committee Reports and Memoranda No. 1632.

  12. 12. Li, G. and Gao, Z. (2008) Exact Solution for Two Dimensional Unsteady Oblique Stagnation Point Flow—An Application of Interacting Shear Flows Theory. Acta Aerodynamica Sinica, 26, 83-86.

期刊菜单