-
In this paper we study the rotating electroosmotic flow of a power-law fluid with Navier slip boundary conditions under high zeta potential subjected to the action of a vertical magnetic field in a variable cross-section microchannel. Without using the Debye–Hückel linear approximation, the finite difference method is used to numerically calculate the potential distribution and velocity distribution of the rotating electroosmotic flow subjected to an external magnetic field. When the behavior index
$n = 1$ , the fluid obtained is a Newtonian fluid. The analysis results in this paper are compared with the analytical approximate solutions obtained in the Debye–Hückel linear approximation to prove the feasibility of the numerical method in this paper. In addition, the influence of behavior index n, Hartmann number Ha, rotation angular velocity$\Omega $ , electric width K and slip parameters$\beta $ on the velocity distribution are discussed in detail. It is obtained that when the Hartmann number Ha > 1, the velocity decreases with the increase of the Hartmann number Ha; but when the Hartmann number Ha < 1, the magnitude of the x-direction velocity u increases with the augment of Ha.-
Keywords:
- high zeta potential /
- Navier slip boundary condition /
- electro magneto-hydrodynamic /
- finite difference method
[1] Stone H A, Stroock A D, Ajdari A 2004 Annu. Rev. Fluid Mech. 36 381Google Scholar
[2] Patel M, Kruthiventi S S H, Kaushik P 2020 Colloids Surf. B 193 111058Google Scholar
[3] Srinivas, Bhadri 2016 Colloids Surf. A 492 144Google Scholar
[4] Nekoubin N 2018 J. Non-Newtonian Fluid Mech. 260 54Google Scholar
[5] Baños R D, Arcos J C, Bautista O, Méndez F, Merchán-Cruz E A 2021 J. Braz. Soc. Mech. Sci. 43 1Google Scholar
[6] Baños R, Arcos J, Bautista O, Méndez F 2020 Defect Diffus. Forum 399 92Google Scholar
[7] 姜玉婷, 齐海涛 2015 64 174702Google Scholar
Jiang Y T, Qi H T 2015 Acta Phys. Sin. 64 174702Google Scholar
[8] Ajdari A 2002 Phys. Rev. E 65 16301Google Scholar
[9] Chang C C, Wang C Y 2011 Phys. Rev. E 84 056320Google Scholar
[10] Song J, Wang S W, Zhao M L, Li N 2020 Z. Naturforsch. A: Phys. Sci. 75 649Google Scholar
[11] Shit G C, Mondal A, Sinha A, Kundu P K 2016 Colloids Surf. A 489 249Google Scholar
[12] 刘全生, 杨联贵, 苏洁 2013 62 144702Google Scholar
Liu Q S, Yang L G, Su J 2013 Acta Phys. Sin. 62 144702Google Scholar
[13] 段娟, 陈耀钦, 朱庆勇 2016 65 034702Google Scholar
Duan J, Chen Y Q, Zhu Q Y 2016 Acta Phys. Sin. 65 034702Google Scholar
[14] Weston M C, Gerner M D, Fritsch I 2010 Anal. Chem. 82 3411Google Scholar
[15] Jian Y J, Chang L 2015 AIP Adv. 5 057121Google Scholar
[16] Xie Z Y, Jian Y J 2017 Colloids Surf. A 529 334Google Scholar
[17] Habib U, Hayat T, Ahmad S, Alhodaly M S 2021 Int. Commun. Heat Mass Transfer 122 105111Google Scholar
[18] Sarkar S, Ganguly S 2017 J. Non-Newtonian Fluid Mech. 250 18Google Scholar
[19] Yang C H, Jian Y J, Xie Z Y, Li F Q 2020 Micromachines 11 418Google Scholar
[20] Xie Z Y, Jian Y J 2017 Energy 139 1080Google Scholar
[21] Wang S W, Li N, Zhao M L, Azese M N 2018 Z. Naturforsch. A: Phys. Sci. 73 825Google Scholar
[22] Xie Z Y, Jian Y J 2014 Colloids Surf. A 461 231Google Scholar
[23] Bird R B, Armstrong R C, Hassager O, Curtiss C F, Middleman S 1978 Phys. Today 31 54Google Scholar
-
图 2 目前数值解与Chang和Wang[9]解析解的比较, 其中
$ \beta = 0, $ $ K = 30, $ $\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 1{\text{ }}{\rm{V}}, {\text{ }}a = 0, {\text{ }}Ha = 0, {\text{ }}S = 0$ Figure 2. Comparison of the current numerical solution with the analytical solution of Chang and Wang [9],
$ \beta = 0, $ $ K = 30, $ $\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 1{\text{ }}{\rm{V}}, {\text{ }}a = 0, {\text{ }}Ha = 0, {\text{ }}S = 0$ 图 3 当无滑移边界条件时, 幂律流体行为指数n对外加磁场的旋转电渗流速度的影响, 其中
$\beta = 0, {\text{ }}K = 10, {\text{ }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, $ ${\text{ }}{\bar \psi _\omega } = 5 \;{\rm{V}}, {\text{ }}a = 0.05, {\text{ }}Ha = 1, {\text{ }}S = $ 1Figure 3. When there is a no-slip boundary condition, the influence of power-law fluid behavior index n on rotating electroosmotic flow velocity with the external magnetic field,
$\beta = 0, {\text{ }}K = 10, {\text{ }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5 \;{\rm{V}}, {\text{ }}a = 0.05, {\text{ }}Ha = 1, {\text{ }}S = 1$ 图 4 当存在滑移边界条件时, 幂律流体行为指数n对外加磁场的旋转电渗流速度的影响, 其中
$ \beta = 0.1, {\text{ }}K = 10, {\text{ }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, $ ${\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }}a = 0.05, {\text{ }}Ha = 1, {\text{ }}S =$ 1Figure 4. When there is a slip boundary condition, the influence of the power-law fluid behavior index n on the rotating electroosmotic flow velocity with an external magnetic field,
$\beta = 0.1, {\text{ }}K = 10, {\text{ }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }}a = 0.05, {\text{ }}Ha = 1, {\text{ }}S = 1$ 图 5 哈特曼数Ha对外加磁场的旋转电渗流速度的影响, 其中
$ n = 0.8, {\text{ }}K = 10, $ $ {\text{ }}\beta = 0.1, {\text{ }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, $ ${\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }}a = $ $ 0.05, {\text{ }}S =$ 1Figure 5. The influence of Hartmann number Ha on the velocity of rotating electroosmotic flow with external magnetic field,
$ n = 0.8, {\text{ }}K = 10, $ $ {\text{ }}\beta = 0.1, {\text{ }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }}a = 0.05, {\text{ }}S = 1 $ 图 6 哈特曼数Ha对外加磁场的旋转电渗流速度的影响, 其中
$ n = 1.2, {\text{ }}K = 10, $ ${\text{ }}\beta = 0.1, {\text{ }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, $ ${\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }}a = $ $ 0.05, {\text{ }}S =$ 1Figure 6. The influence of Hartmann number Ha on the velocity of rotating electroosmotic flow with external magnetic field,
$ n = 1.2, {\text{ }}K = 10, $ $ {\text{ }}\beta = 0.1, {\text{ }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }}a = 0.05, {\text{ }}S = 1 $ 图 7 旋转角速度
$\varOmega $ 对外加磁场的旋转电渗流速度的影响, 其中$ n = 0.8, $ $ {\bar \psi _\omega } = {\text{5}}\;{\rm{V}}, $ $ a = 0.05, {\text{ }}Ha = 1, {\text{ }}S = 1 $ (a)$K = 10, {\text{ }}\beta = $ $ 0.1;$ (b)$ K = 10, {\text{ }}\beta = 0.1; $ (c)$ K = 10, {\text{ }}\beta = 0; $ (d)$ K = 20, {\text{ }}\beta = 0.1. $ Figure 7. The influence of the rotational angular velocity
$\varOmega $ on the rotational electroosmotic flow velocity of the external magnetic field,$ n = 0.8, $ $ {\bar \psi _\omega } = {\text{5 }}V, $ $ a = 0.05, {\text{ }}Ha = 1, {\text{ }}S = 1 $ (a)$ K = 10, {\text{ }}\beta = 0.1; $ (b)$ K = 10, {\text{ }}\beta = 0.1; $ (c)$ K = 10, {\text{ }}\beta = 0; $ (d)$K = 20, $ $ {\text{ }}\beta = 0.1$ 图 8 旋转角速度
$\varOmega $ 对外加磁场的旋转电渗流速度的影响, 其中$ n = 1.2, $ $ {\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }} $ $ a = 0.05, {\text{ }}Ha = 1, {\text{ }}S = 1 $ (a)$\beta = 0.1, {\text{ }}K = $ $ 10.$ (b)$ \beta = 0.1, {\text{ }}K = 10. $ (c)$ \beta = 0, {\text{ }}K = 10. $ (d)$\beta = 0.1, {\text{ }}K = 30$ Figure 8. The influence of the rotational angular velocity
$\varOmega $ on the rotational electroosmotic flow velocity of the external magnetic field,$ n = 1.2, $ $ {\bar \psi _\omega } = 5{\text{ }}V, {\text{ }} $ $ a = 0.05, {\text{ }}Ha = 1, {\text{ }}S = 1 $ (a)$ \beta = 0.1, {\text{ }}K = 10. $ (b)$ \beta = 0.1, {\text{ }}K = 10. $ (c)$ \beta = 0, {\text{ }}K = 10. $ (d)$\beta = 0.1, $ $ {\text{ }}K = 30$ 图 9 电动宽度 K 对外加磁场的旋转电渗流速度分布的影响, 其中
$ n = 0.8, $ $ \beta {\text{ = 0}}{\text{.1, }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, $ ${\text{ }}a = 0.05, {\text{ }}Ha = $ $ 1, {\text{ }}S = $ 1Figure 9. The influence of the electric width K on the velocity distribution of rotating electroosmotic flow with external magnetic field,
$n = 0.8, {\text{ }}\beta {\text{ = 0}}{\text{.1, }}\varOmega = 100{\text{ rad/s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}V, {\text{ }}$ $ a = 0.05, {\text{ }}Ha = 1, {\text{ }}S = 1 $ 图 10 电动宽度 K 对外加磁场的旋转电渗流速度分布的影响, 其中
$ n = 1.2, $ $ \beta {\text{ = 0}}{\text{.1, }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }}a = 0.05, {\text{ }}Ha = $ $ 1, {\text{ }}S = $ 1Figure 10. The influence of the electric width K on the velocity distribution of rotating electroosmotic flow with external magnetic field,
$ n = 1.2, {\text{ }}\beta {\text{ = 0}}{\text{.1, }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }} $ $ a = 0.05, {\text{ }}Ha = 1, {\text{ }}S = 1 $ 图 11 电动宽度 K 对外加磁场的旋转电渗流速度分布的影响, 其中
$ n = 1.2, $ $\beta {\text{ = 0, }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }}a = 0.05, {\text{ }}Ha = $ $ 1, {\text{ }}S =$ 1Figure 11. The influence of the electric width K on the velocity distribution of rotating electroosmotic flow with external magnetic field,
$ n = 1.2, {\text{ }}\beta {\text{ = 0, }}\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }} $ $ a = 0.05, {\text{ }}Ha = 1, {\text{ }}S = 1 $ 图 12 滑移参数
$\beta $ 对外加磁场的旋转电渗流速度的影响, 其中$ n = 0.8, $ $ K = 10, $ $\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }}a = 0.05, {\text{ }}Ha = $ $ 1, {\text{ }}S =$ 1Figure 12. The influence of the slip parameter
$\beta $ on the rotating electroosmotic flow velocity with an external magnetic field,$ n = 0.8, $ $ K = 10, $ $ \varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }}a = 0.05, {\text{ }}Ha = 1, {\text{ }}S = 1$ 图 13 滑移参数
$\beta $ 对外加磁场的旋转电渗流速度的影响, 其中$ n = 1.2, $ $ K = 10, $ $ \varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }}a = 0.05, {\text{ }}Ha = $ $ 1, {\text{ }}S = 1 $ Figure 13. The influence of the slip parameter
$\beta $ on the rotating electroosmotic flow velocity with an external magnetic field,$ n = 1.2, $ $ K = 10, $ $\varOmega = 100{\text{ }}{\rm{rad}}/{\rm{s}}, {\text{ }}{\bar \psi _\omega } = 5{\text{ }}{\rm{V}}, {\text{ }}a = 0.05, {\text{ }}Ha = 1, {\text{ }}S = 1$ -
[1] Stone H A, Stroock A D, Ajdari A 2004 Annu. Rev. Fluid Mech. 36 381Google Scholar
[2] Patel M, Kruthiventi S S H, Kaushik P 2020 Colloids Surf. B 193 111058Google Scholar
[3] Srinivas, Bhadri 2016 Colloids Surf. A 492 144Google Scholar
[4] Nekoubin N 2018 J. Non-Newtonian Fluid Mech. 260 54Google Scholar
[5] Baños R D, Arcos J C, Bautista O, Méndez F, Merchán-Cruz E A 2021 J. Braz. Soc. Mech. Sci. 43 1Google Scholar
[6] Baños R, Arcos J, Bautista O, Méndez F 2020 Defect Diffus. Forum 399 92Google Scholar
[7] 姜玉婷, 齐海涛 2015 64 174702Google Scholar
Jiang Y T, Qi H T 2015 Acta Phys. Sin. 64 174702Google Scholar
[8] Ajdari A 2002 Phys. Rev. E 65 16301Google Scholar
[9] Chang C C, Wang C Y 2011 Phys. Rev. E 84 056320Google Scholar
[10] Song J, Wang S W, Zhao M L, Li N 2020 Z. Naturforsch. A: Phys. Sci. 75 649Google Scholar
[11] Shit G C, Mondal A, Sinha A, Kundu P K 2016 Colloids Surf. A 489 249Google Scholar
[12] 刘全生, 杨联贵, 苏洁 2013 62 144702Google Scholar
Liu Q S, Yang L G, Su J 2013 Acta Phys. Sin. 62 144702Google Scholar
[13] 段娟, 陈耀钦, 朱庆勇 2016 65 034702Google Scholar
Duan J, Chen Y Q, Zhu Q Y 2016 Acta Phys. Sin. 65 034702Google Scholar
[14] Weston M C, Gerner M D, Fritsch I 2010 Anal. Chem. 82 3411Google Scholar
[15] Jian Y J, Chang L 2015 AIP Adv. 5 057121Google Scholar
[16] Xie Z Y, Jian Y J 2017 Colloids Surf. A 529 334Google Scholar
[17] Habib U, Hayat T, Ahmad S, Alhodaly M S 2021 Int. Commun. Heat Mass Transfer 122 105111Google Scholar
[18] Sarkar S, Ganguly S 2017 J. Non-Newtonian Fluid Mech. 250 18Google Scholar
[19] Yang C H, Jian Y J, Xie Z Y, Li F Q 2020 Micromachines 11 418Google Scholar
[20] Xie Z Y, Jian Y J 2017 Energy 139 1080Google Scholar
[21] Wang S W, Li N, Zhao M L, Azese M N 2018 Z. Naturforsch. A: Phys. Sci. 73 825Google Scholar
[22] Xie Z Y, Jian Y J 2014 Colloids Surf. A 461 231Google Scholar
[23] Bird R B, Armstrong R C, Hassager O, Curtiss C F, Middleman S 1978 Phys. Today 31 54Google Scholar
Catalog
Metrics
- Abstract views: 3392
- PDF Downloads: 47
- Cited By: 0