A THEORETICAL STUDY OF FLOP-OVER OSCILLATIONS IN CONNECTED CHANNELS FOR NANOFLUIDS WITH CONCENTRATION DEPENDENT VISCOSITY Текст научной статьи по специальности «Физика»

ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА

2019

• ФИЗИКА •

Вып. 3

УДК 532.54; 534-14 PACS 47.55.pb; 47.57.J-

Теоретическое исследование перебросовых

каналах при учете зависимости вязкости среды от концентрации частиц

Д. З. Аллен, Р. Д. Бэнкс ь , В. А. Демин

a Новый колледж, Оксфордский университет, ул. Холивелл, Оксфорд, OX1 3BN email: dominic.allen@new.ox.ac.uk

b Санкт-Хью колледж, Оксфордский университет, ул. Св. Маргарет, Оксфорд, OX2 6LE email: robert.banks@st-hughs.ox.ac.uk

c Пермский государственный национальный исследовательский университет 614990, Пермь, ул. Букирева, 15 email: demin@psu.ru

В работе проведено теоретическое исследование тепловой конвекции наножидкости в подогреваемых снизу связанных каналах. Рассматривались связанные каналы конечной высоты с твердыми высоко теплопроводными границами. Прямое численное моделирование основывалось на уравнениях в частных производных, полученных с помощью галеркинского осреднения по продольной координате. Результирующие уравнения решались методом конечных разностей. По результатам расчетов показано, что нормальная термодиффузия в рассматриваемой конвективной среде оказывает решающее влияние на формирование специфических нелинейных перебросовых колебаний. В ходе численного моделирования были получены поля скорости, температуры и концентрации наночастиц в поперечном сечении на разной высоте каналов. Форма и период этих колебаний в зависимости от надкритичности были проанализированы для разных значений управляющих параметров. Квадратичная зависимость вязкости коллоидного раствора от концентрации наночастиц принималась во внимание в теоретической модели. Показано, что период перебросовых колебаний должен возрастать при уменьшении средней начальной концентрации наночастиц. Нормировка периода колебаний на критическое число Рэлея позволила получить универсальный закон, справедливый для наножидкостей c разными зависимостями вязкости от концентрации. Качественное сравнение показало согласие результатов расчета с известными экспериментальными данными.

Ключевые слова: наножидкости; связанные каналы; тепловая конвекция; термодиффузионное разделение; перебросовые колебания

Поступила в редакцию 13.08.2019; принята к опубликованию 30.08.2019

A theoretical study of flop-over oscillations in connected channels for nanofluids with concentration dependent viscosity

D. Z. Allena, R. J. Banksb, V. A. Deminc

наножидкости в связанных

© Аллен Д. З., Бэнкс Р. Д., Демин В. А., 2019

распространяется на условиях лицензии

Creative Commons Attribution 4.0 International (CC BY 4.0).

a New College, University of Oxford, Holywell Street, Oxford, OX1 3BN email: dominic.allen@new.ox.ac.uk

b St. Hughs College, University of Oxford, St. Margaret's Road, Oxford, OX2 6LE email: robert.banks@st-hughs.ox.ac.uk c Perm State University, Bukireva St. 15, 614990, Perm email: demin@psu.ru

Thermal convection of a nanofluid in connected channels heated from below was studied theoretically. Thin channels with boundaries of high thermal conductivity were considered. Direct numerical simulation was based on the set of partial differential equations derived with the help of Galerkin method. The final system was solved by the method of finite differences. It has been shown that normal thermodiffusion exerts principal influence on the formation of particular non-linear flop-over oscillations. Cross-sections of concentration, velocity and temperature were obtained at various heights in the channel. The form and period of these oscillations were analyzed with respect to su-percriticality for different values of governing parameters. A quadratic dependence of viscosity on concentration was taken into account. The period of oscillation was found to decrease with the growth of this effect. As a result of the normalization procedure the general law for the period of oscillation has been found. Qualitative agreement between calculation and experiment has been demonstrated.

Keywords: nanofluid; connected channels; thermal convection; thermodiffusive separation, flop-over oscillations

Received 13.08.2019; accepted 30.08.2019 doi: 10.17072/1994-3598-2019-3-12-19 1. Introduction

For non-uniformly heated fluids, convective flows occur due to the dependence of density on temperature [1]. In the case of binary mixtures, the inhomogeneous concentration distribution also plays a role in convection, and this gives rise to more complex behavior [2]. These effects are dependent on the temperature and concentration gradients. In a mixture, the temperature gradients are responsible for the occurrence of concentration gradients (thermodiffusion). Thermodiffusion can be considered as a cross effect in the general theory of Onsager [3]. The case where a positive temperature gradient creates a negative concentration gradient of the heavier component is known as positive (normal) thermodiffusion.

In vertical connected channels heated from below, positive thermodiffusion and thermoconvection together give rise to observed "flop-over" oscillations in certain intervals of the Rayleigh number, where the velocity of the fluid in the direction of the channel has approximately a square-wave profile, Glukhov and Putin [4]. The calculation of flop-over oscillations was first carried out for binary mixtures in [5], where in accordance with experiment the boundaries were considered to have high thermal conductivity compared to the fluid.

In this paper, we investigate these flop-over oscillations in thin connected channels for nanofluids with heavy particles and viscosity dependent on concentration. Experiments have shown that concentration can

have a significant contribution to viscosity for many natural and artificial nanofluids [6, 7]. We suppose that this dependence may be important in providing a more accurate description of flop-over oscillations.

2. Statement of the problem

2.1. Geometry of set-up

Let us consider connected channels (convective loop) oriented vertically and heated from below, as in fig. 1. The channels have finite height h and half-width d, with the condition h >> d.

Fig.1. Geometry of the problem and coordinate system.

The boundaries are assumed to have high thermal conductivity, in accodance with experiments [4,5]. In the experiments, this was achieved by embedding the channels in a homogeneous metal block. The temperature profile along the centre of the opaque channels was determined using thermocouples. This allowed the flow intensity along the centre of the channel to be estimated quantatively. For small Ray-leigh number, it was found that the velocity of the fluid was one dimensional, directed along the channel. At the same time the fields of temperature and concentration were observed to be a function of all spatial coordinates. It was established that the complexity of these fields leads to a highly non-trivial time evolution of this hydrodynamic system.

In our system, the z-axis is oriented vertically upwards, the x and y axes lie in the channel's cross section. The temperature difference between the top and bottom of the channels is 0. Gravity is directed vertically downwards g = (0, 0,— g), see fig. 1.

2.2. Basic equations

Now let us take the well-known generalized equations of thermal convection in the Boussinesq approximation to describe the behavior of an incompressible fluid in connected channels. This system includes the generalized Navier-Stokes equation, mass conservation law, equations of heat and admixture transfer [1,8]:

dv _ 1 _

--h vV v =--Vp+vAv +2 V/yV v +

dt p

+Vz/x Vx® + g 3T-f3C 7, (2.1)

dT

-+ vV T — xAT .

dt

divu = 0 , (2.2)

dC ~dt

+ vV G = D AG + aAT .

(2.3)

Here, v , p, Tand C are dimensional fields of velocity, pressure, temperature and the concentration of the heavier component. Parameters /, pt, Pc are positive coefficients of thermal diffusivity, thermal expansion and dependence of density on concentration, respectively; D and a are the diffusion and thermodiffusion coefficients; p is the average density of the fluid; g is the gravitational acceleration; and 7 is the unit vector oriented vertically upward.

The kinematic viscosity, v, we take to be dependent on concentration quadratically [6], with the form:

(2.4)

where Ad, £d are positive dimensional parameters dependent on the properties of the nanofluid; vo is the kinematic viscosity of the pure liquid.

The sides of the channels are rigid therefore the velocity vanishes on the boundaries. As a result of the heating from below a linear temperature distribution takes place on the boundaries T0 = - z0/h + ©. The unknown fields of velocity v and temperature T = T0 + T' in equations (2.1), (2.2) must satisfy to the following boundary conditions:

v\ =0, lr

T' I =0.

Ir

These account for the existence of non-slip and nonpenetration conditions on all sides of the channels, and the high thermal conductivity of the boundaries.

Let us continue the numerical analysis of the equations (2.1)-(2.3) in terms of non-dimensional variables. We shall use following set of units during the simulation: length [Ax, Ay, Az] - d; time [i] - d2/vo; velocity [v] - vo/d; temperature [T] - &, concentration [C] - Pt&/Pc; pressure [p] - pvo2/d2. The equations of thermal convection in non-dimensional form can be written as

dv ~dt

+ v\7 v = -Vp + vAv + 2 V/A7 v +

ünU

x Vxv +- T-C 7, (2.5)

Pr

dT 1

-+ iJv T = —AT.

dt Pr

divv = 0 , (2.6)

dC 1

-+ vV C = — AG + eAT

dt Sc

v = 1 + AC + £C2.

(2.7)

(2.8)

The following non-dimensional governing parameters are used:

Ra :

£/?<9d3 v v aß

^-, Pr = —, Sc = —, £ — —L

D ßt

X

ß©

ßt0

ß

1 r.

H = —. d ■

where Ra, Pr and Sc are the Rayleigh, Prandtl and Schmidt numbers respectively; s is the non-dimensional thermodiffusion parameter; A and £ describe the dependence of viscosity on concentration; H is the non-dimensional height.

In the calculations, we use the straight trajectory approximation for the velocity v — (0, 0, u(.r. 1/. / )).

The form of u(x, y, t) automatically satisfies the in-compressibility requirement. As h >> d, we limit our considerations to the vertical channels of the convec-tive loop.

2

d

The non-dimensional flux of matter is given by J = - V C' + eT . The non-penetration boundary condition implies

J =0. «Ii

(2.9)

We impose the condition that the net density flux through the two channels is zero:

If

u(2) dxdy = 0 ,

VF

0.

In the following calculations we omit the primes for greater clarity. In the terms of the deviation fields the equation system has the form

du dp

— =--- + vA, и +

dt dz L

where the superscripts denote the left and right channels correspondingly.

3. Mechanical equilibrium

For small Rayleigh number, a mechanical equilibrium state exists, where v 0 and djdt = 0. The

non-dimensional temperature distribution at equilibrium is

T = - — + 1. 0 H

Accounting for the condition on flux and equation (2.7), the concentration field is described by

The integration constant C can be expressed in terms of the mean concentration C :

Let us emphasize the net flux of matter is equal to zero in mechanical equilibrium.

4. Method of solution

In non-equilibrium states, we have 3D fields of concentration, temperature and velocity, which in general depend on time. The full concentration and temperature fields can be written as C = Co + C and T = To + T', where C and T' are the deviations from mechanical equilibrium. We reformulate the equations (2.5)-(2.7) in terms of these deviations and solve numerically.

To simplify the application of boundary condition (2.9), we introduce the new function F = C + sT. The associated boundary condition is

11„U

+ Vz/V u +- (1 + s)T — F , (4.1)

Pr

dT dT 1 A ^ и „

--h и-= —ДГ + —, (4.2)

dt dz Pr H

dF dF 1 e

— + u— = —AF + —AT. (4.3) dt dz Sc Pr

The plane Laplace operator is used in equation (4.1)

д2 Э2

A = ■

dx2 dy2

Experimental results [4] indicate the vertical distribution of temperature can be written in terms of two harmonic functions

T = Tt (x, y, t) sin

7TZ

H

+ T2(x, y,t) cos

7TZ

H

Here Ti and T2 are amplitudes dependent on the cross section coordinates and time only. Considering the symmetries of equations (4.i)-(4.3) we suggest the following expansion for F:

27TZ

F = F1(x,y,t) + F2 (x, y, t) cos

+F (x, y, t) cos

H

We integrate (4.1) over both channels to remove from this equation the pressure gradient as it does not affect buoyancy forces. For equations (4.2) and (4.3) we apply the Galerkin procedure to find the equations for amplitudes. The resultant system of partial differential equations for the amplitudes u, T1, T2, Fl, F2 and F3 can be found in the appendix. The equation for u involves only terms anti-symmetric for left and right channels.

The equation system for amplitudes was solved numerically by the method of finite differences in combination with the explicit scheme [9]. One sided and central differences were used to approximate the time and spatial coordinate derivatives respectively, including the convective and diffusive terms [10]. Numerical stability was achieved by choosing an appropriate time step. The initial transitional flows in the results were ignored and stable regimes were investigated.

5. Results and discussion

In the course of direct numerical simulation, we focused our attention on the regime of flop-over oscil-

s

n I г

lations. The characteristic form of these oscillations is illustrated in figure 2.

10 и

-10

4000

t 8000

Fig. 2. The velocity at the center of the channel for A = 0.2, £ = 0.5, mean concentration C0 = 0.1, Ra= 53.

We observe that an approximately steady flow is established with periodic sharp changes in direction. The mechanism responsible for the observed effects is mainly attributed to the horizontal thermodiffusive separation of the mixture. From experiment [4], the horizontal temperature gradients d/d = 3 K/cm are much larger than the vertical gradients 0/h = 0.3 K/cm, which has a characteristic diffusive separation time of h2/D ~ 300 days. The horizontal gradients have characteristic separation time d2/D ~ 3 hours, which coincides in order of magnitude with the oscillation period of the ferrofluid from experiment [4]. Considering first the channel with a flow moving up, the fluid in the center has greater velocity and temperature than the surrounding fluid (fig.3 (a)). Therefore, as a result of thermodiffusion, the heavier admixture diffuses towards the slower fluid layer on the boundary. Over time, the channel with upwards flow accumulates the heavier admixture, whilst the opposing channel loses it. When the channel becomes sufficiently heavy, the flow stops abruptly and changes direction.

For our simulation the following values for the governing set of parameters were used: Pr = 10, Sc = 75, s = 0.3 and H = 30.5. Between simulations the parameters A, £, C and Ra were varied in order to

analyze the effect of concentration dependence of viscosity.

In figure 3 the channel with upward velocity is presented. Figure 3(a) shows the temperature deviation profile at various heights along the channel. In our theoretical model we neglect the horizontal parts of the convective loop therefore the full temperature at the ends of the channels must be continuous. At the bottom of the channel the temperature profile has a smaller amplitude than in the upper parts of the channel. This behavior is expected as the inflow at the bottom of the channel has just exited the colder channel.

The concentration deviation field at various heights is illustrated in fig. 3(b). For z = h/4 and lower the concentration is greater in the center of the channel than near the boundaries. We can see that the concentration of nanoparticles diffuses towards the boundaries as fluid moves up the channel. Towards the top of the channel (z = 3 h/4) the concentration in the center is smaller. These calculations support the thermodiffusion mechanism outlined earlier in the text and in [5].

The aim of our work is to take into account the effect of concentration dependence of viscosity (2.8). It was not necessary to consider this effect in molecular solutions [5], but experimental data suggest that this may be significant for nanofluids [6].

In our model, the parameters A and £ determine the strength of this dependence. The mean concentration Ca also affects the convective flows of nanofluids. We

are interested in the period of flop-over oscillations for different values of these parameters, as this is accurately measured by experiment [4, 5]. We can obtain from fig. 2 the dimensional value of oscillation period т taking the viscosity v = 2.7-10-6 m2/s and channel width d = 2 mm, we have т~ 1.6 hours, which is qualitatively verified by the result from experiment [4]. The quantitative data of period dependence on Ray-leigh number from simulation are presented in fig. 4. It is seen that the period of oscillations increases with Rayleigh number (which is proportional to the temperature difference) approximately linearly. Properties of the nanofluid A and £ are considered to be constant. With the decrease of mean concentration the period becomes greater in agreement with experiment [4]. The threshold of flop-over oscillations depends on the mean concentration, A and £. Flop-over oscillations start for each mean concentration at a different critical Rayleigh number Rac, which are the lowest values of Ra in fig. 4. Let us attempt to find a universal law for these dependences.

To visualize final data we introduce supercriticality U = (Ra - Rac)/Rac. By plotting the period against su-percriticality for each mean concentration value we produce a generalized law (fig. 5).

3h 4

h 2

h 4

3h 4

(b)

(a)

Fig. 3. The fields of temperature (a) and concentration (b) in a cross section of the channel on different heights for X = 0.2, c = 0.5, mean concentration C = 0.1, Ra = 56.

7000

6000

5000

4000

3000

X A = 0, £ = 0 OA = 0.2, Ç = 0.5, = 0.1

*A = 0.2, £ = 0.5,C0 = 0.2 * 1

r

x * O

x o o * LQ-it-i-'- Ra

50

55

60

65

70

75

Fig. 4. Non-dimensional period of oscillations against Rayleigh number for different values of parameters.

7000

T 6000

5000

4000

300C

2000

x A = 0, £= 0 O A = 0.2, £= 0.5, C = 0.1 *A = 0.2,i = 0.5,C0 = 0.2

ß

0.1

0.2

0.3

0.4

0.5

Fig. 5. Normalized period of oscillations against Rayleigh number for some values ofparameters.

The linear trend line in fig.5 is given by:

T(fi) = 7.5/./ + 3.0 ■ 103.

The law suggests that, for a given nanofluid, there is a minimum period for flop-over oscillations regardless of mean concentration. For our parameters, this is ~ 3-103. The period of flop- over oscillation in a given nanofluid (constant Sc, Pr, and s) for small values of supercriticality can to a good approximation be described as having a linear dependence on 0, a convenient parameter to measure in experiments. The flopover regime's dependence on concentration must be contained in the critical value of temperature 0c = 0c(l, £ C'o,...). Our numerical results so far suggest 0c increases withCo. Many further simulations would

have to be performed to determine a more explicit form of 0c; however its value is easily attainable from experiment.

Though the model provided by this paper allows for a description of flop-over oscillations and dependence on Rayleigh number and mean concentration it is by no means complete. In this simulation we increased the Rayleigh number (experimentally realized by increasing the temperature difference) without changing

z

other parameters dependent on temperature difference

(e.g. A, C ). Further research may wish to take into

account the temperature dependence of the other governing parameters. Possible future work could even include varying the Schmidt number with potential concentration dependence.

6. Conclusion

Non-steady thermal convection of a nanofluid in connected channels with boundaries of high heat conductivity was investigated theoretically guided by experiment. In our model we assume the liquid moves predominantly vertically, adjacent to heat-conducting boundaries. We suggest that the separation of components in the horizontal plane by thermodiffusion is responsible for the complex "flop-over" oscillatory regime. Our direct numerical simulation based on hydrodynamic equations for the convective loop is in qualitative agreement with the results of experiments

[4, 5]. The simulations in this paper have demonstrated that the introduction of concentration dependent viscosity clearly affects the period of flop-over oscillations in a convective loop. This brings numeric simulations qualitatively in line with the results of experiments on nanofluids.

We are grateful for the financial support of the scholarship provided by British Petroleum. We thank the University of Oxford's Career Service and the International Department at Perm State University for providing the research opportunity.

Appendix

The following final partial differential equation system gives the time evolution for amplitudes Ti(x, y, t), 72(x, y, t), u(x, y, t), Fi(x, y, t), F2(x, y, t), F=(x, y, t). After solving for these amplitudes, we return to the full values of temperature, velocity and concentration.

du

~dt

l + xc +ÇC +

£Ü 12

А и

2ЯЯа

Pr-,

1 + ét T.--F. 1 2 1

Щ2 + F* + F* - 2eF2T2 + e2 T2 + т22 -

4e

Зтг

3F1+F3 Tl Alu +

2F^Fî+F^F2+F^F3 e T2VF2+F2VT2 + S2 T^Ti+T2VT2

4e

Зтг

3 TJVFJ+FJVTJ - TJVFJ+FJVTJ V±u\

dt h 2 Pr

A±Ti -

t:

4 и дТ2 ижт _ J_ TV H * dt. H 1 ~ Pr

AlT2-

z

dt, H

Sc

тгРг

dF3 4v_F =l_ dt m 2 ~ se

A±Tl

ALF3- 4

7Г

H J

2

F„

dF2 16 и F dt 3H 3 ~~ Pr

A T

_L 2

4e

Зтг Pr

AlTi -

Тг

To--

2 Se

Fn

V -— e +—e

± dx x dy y

2

2

Список литературы

1. Гершуни Г. З., Жуховицкий Е. М. Конвективная устойчивость несжимаемой жидкости. М.: Наука, 1972. 392 с.

2. Nield D. The thermohaline Rayleigh-Jeffreys problem. Journal of Fluid Mechanics. 1967. Vol. 29. N. 3. P. 545-558. DOI.10.1017/S0022112067001028

3. Onsager L. Theories and problems of liquid diffusion // Annals of the New York Academy of Sciences. 1945. Vol. 46. N. 5. P. 241-265.

4. Glukhov A. F., Putin G. F. Convection of magnetic fluids in connected channels heated from below // Fluid Dynamics. 2010. Vol. 45. N. 5. P. 713-718.

5. Glukhov A. F., Demin V. A., Putin G. F. Binary-mixture convection in connected channels heated from below // Fluid Dynamics. 2007. Vol. 42. N. 2. P. 160-169.

6. Rudyak V. Ya. Viscosity of nanofluids - Why it is not described by the classical theories // Advances in nanoparticles. 2013. N. 2. P. 266-279.

7. Sprenger L., Lange A., Odenbach S. Thermodiffusion in concentrated ferrofluids: A review and cur-

rent experimental and numerical results on nonmagnetic thermodiffusion // Physics of Fluids. 2013. Vol. 25, 122002.

8. Ландау Л. Д., Лифшиц Е. М. Гидродинамика. М.: Наука, 1986. 736 с.

9. Роуч П. Вычислительная гидродинамика. М.: Мир, 1980. 616 с.

10. Fletcher C.A.J. Computational techniques for fluid dynamics. Springer series in computational physics. Vol. 1, 2005, 406 p.

References

1. Gershuni G. Z., Zhukhovitskii E.M. Convective stability of incompressible fluids. Jerusalem: Keter Publishing House, 1976, 330 p.

2. Nield D. The thermohaline Rayleigh-Jeffreys problem. Journal of Fluid Mechanics, 1967, vol. 29, no. 3, pp. 545-558. D0I:10.1017/S0022112067001028

3. Onsager L. Theories and problems of liquid diffusion. Annals of the New York Academy of Sciences, 1945, vol. 46, no. 5, pp. 241-265.

4. Glukhov A. F., Putin G. F. Convection of magnetic fluids in connected channels heated from below. Fluid Dynamics, 2010, vol. 45, no. 5, pp. 713-718.

5. Glukhov A. F., Demin V. A., Putin G. F. Binary-mixture convection in connected channels heated from below. Fluid Dynamics, 2007, vol. 42, no. 2, pp. 160-169.

6. Rudyak V. Ya. Viscosity of nanofluids - Why it is not described by the classical theories. Advances in nanoparticles, 2013, no. 2, pp. 266-279.

7. Sprenger L., Lange A., Odenbach S. Thermodiffusion in concentrated ferrofluids: A review and current experimental and numerical results on nonmagnetic thermodiffusion. Physics of Fluids, 2013, vol. 25, 122002.

8. Landau L. D., Lifshitz E. M. Course of Theoretical Physics, vol. 6. Fluid Mechanics. Oxford: Butterworth-Heinemann, 1987. 554 p.

9. Roache P. Computational fluid dynamics. Albuquerque, New Mexico, Hermosa Pub., 1976. 446 p.

10. Fletcher C. A. J. Computational techniques for fluid dynamics. Springer series in computational physics. Vol. 1, 2005, 406 p.

Просьба ссылаться на эту статью в русскоязычных источниках следующим образом:

Аллен Д. З., Бэнкс Р. Д., Демин В. А. Теоретическое исследование перебросовых колебаний наножидко-сти в связанных каналах при учете зависимости вязкости среды от концентрации частиц // Вестник Пермского университета. Физика. 2019. № 3. С. 12-19. doi: 10.17072/1994-3598-2019-3-12-19

Please cite this article in English as:

Allen D. Z., Banks R. J., Demin V. A. A theoretical study of flop-over oscillations in connected channels for nanofluids with concentration dependent viscosity. Bulletin of Perm University. Physics, 2019, no. 3, pp. 12-19. doi: 10.17072/1994-3598-2019-3-12-19