#### 1. Introduction

Colloidal suspensions or, as they are now called, nanofluids, are widely used in various areas of modern technology. For example, magnetic fluids are used to polish optical components, and suspensions of silica particles in liquid crystals significantly improve the characteristics of optical storage devices [1-2]. We also note their use in chemical processes (catalysis), in the creation of new drugs, lubricants and so on [3-7]. With the growth of the productivity of electronic devices and the development of high-energy technologies, it has become necessary to create efficient cooling systems and control large heat flows [8]. Promising developments have emerged related to molecular computers based on switchable bistable molecules or their aggregates. One way to intensify heat transfer is to increase the thermal conductivity of a liquid by adding solid particles with high thermal conductivity [9-12].

Particularly interesting for the creation of such suspensions are nanoparticles. Recent studies have shown that liquid-phase media, in which nanoparticles from wide-gap semiconductors or dielectrics are used as a dispersive component, are very effective for the realization of a number of nonlinear optical effects [13]. In these media, unlike homogeneous media, the nonlinear optical response arises from the change in the refractive index induced by the light wave and the absorption coefficient due to the phenomena of thermal diffusion and electrostriction of particles [14-18]. At the same time, in our opinion, physical mechanisms connected in particular with the processes of heat and mass transfer in such media require additional investigation. It should be added that a number of mathematical problems associated with analytical methods for solving the corresponding equations are also far from being resolved.

#### 2. Theoretical model

The paper considers a liquid-phase medium with nanoparticles irradiated by a light beam with a Gaussian intensity profile. As a result of the action of the light field in the medium, gradients of temperature and concentration arise, which cause the processes of heat and mass transfer. These phenomena are described by a system of balance equations for temperature and particles, taking into account the concentration convection [19], written in one-dimensional form (without using the Boussinesq approximation):

In the heat equation, the term responsible for the Dufour effect is omitted, in view of its smallness. The processes of sedimentation are also neglected. Here, the following designations are used: *T* is the temperature of the medium, *C = C(r,t) = m*
*m* is the mass of the whole medium), *C
*

*C*) is the coefficient of thermal conductivity

*, I*

*v*= η/h is the rate of concentration convection, which we assume to be constant. In this case, η is the kinematic viscosity, and

*h*is the characteristic distance (in our case, we take

*h*≈
10

*x*

*k*B
is Boltzmann's constant, c is the speed of light in vacuum,

*n*e f
is the effective refractive index of the medium. The last term expresses the contribution from the particle flux and is related to the effect of the gradient force from the side of the light field (electrostriction of particles).

Let us take into account the fact that the temperature establishment processes go faster than the diffusion processes. This makes it possible to study the latter against a background of a stationary temperature. We will consider the case of small concentrations: C

Further, the concentration dependence of the thermal conductivity can be represented in the form:

where **
**. This kind of dependence was theoretically found in [20] and was experimentally confirmed in publications [21,22]. In the equations of heat conduction and diffusion, divergent terms are taken in the form:

(4)
The validity of this approximation for the first formula is based on the use of equality (3), which leads to a reasonable inequality,

(5)
Assuming that C

(6)
where

(7)

(8) We rewrite the problem (7)–(8) in dimensionless variables:

(9)

(10)
where **
**Peclet number,

b = λ 0 S T α 0 I ˜ 0 ( 1 - p C 0 ) x 0 2 ,

It is not possible to find the exact analytic solution of equation (7) with initial-boundary conditions (8). The Peclet number estimate for typical values of the quantities in our case gives

This value indicates the predominance of convection processes over molecular diffusion. In view of this fact, we rewrite problem (7)–(8) in a truncated form:

(11) Its exact solution can be written in the form:

(12)
where
*A* can be determined from the condition of conservation of the number of particles in the system.

Figure 2 shows the graphs of the function u (ρ) at increasing times τ.

Thus, accounting for the concentration convective flow, even in the one-dimensional approximation, leads to amplitude-modulated travelling waves.

Let us now consider the state of mechanical equilibrium of a nanofluid, that is, when its macroscopic flow is absent and the equilibrium distributions of temperature and concentration are described by the equations:

(13)

(14) The solution of this system is the expression:

(15)
where

The constant *B* in formula (15) can be determined from the condition of constancy of the average concentration.

#### 3. Discussion

Note that the convective processes in a colloidal suspension of nanoparticles taking into account sedimentation (in the Benard scheme) were studied in [23-25]. The mode of running concentration waves was detected. In our approach, the effect of radiation due to the action of gradient forces on the side of the light field, being an additional source of convection, leads to a completely different spatial distribution of concentration. We note that the amplitudes of the modulated waves in solution (12) are damped and have a spatially complex structure. In particular, an analysis of the exponent shows that with positive thermal diffusion (b

Solution (15) shows that in our case the convection of a nanofluid has two processes: thermal diffusion and electrostriction, which can be either unidirectional or multidirectional. It is easy to see from (15) that if the contribution from electrostriction is small (m

We understand that a number of questions require more detailed study, for example, analysis of convection stability, numerical modelling of these processes and their study in a 2D measurement. We will devote our attention to these problems in future research.

This work was supported by Ministry of Education and Science of Khabarovsk krai (Russia).