Movement and Heat Transfer of Particles in Inhomogeneous and Nonisothermal Rapidly Oscillating Fluid Flow

Based on the Krylov-Bogolyubov method of averaging the closed system of equations for particle motion and temperature in inhomogeneous rapidly oscillating velocity and temperature of fluid phase is derived. It is shown that the particle movement in a rapidly oscillating fluid velocity field occurs not only under the force of gravity and resistance, but also under force of migration. The migration force is the result of particle inertia and in homogeneity of oscillation of velocity field of the carrier phase. Effects of dynamic and thermal relaxation times of particle and gravity force have been studied. It is shown possibilities of accumulation of particles under the combined action of gravity and migration forces. For a linear dependence of the amplitude of velocity and temperature fluctuations of the fluid an analytical solution was presented. The analytical solutions have been found in good agreement with the results of numerical solution of system of equations of motion and heat transfer of particle.


Introduction
Fluid or gas turbulent flows with particles or droplets widely occur in natural phenomena and are intensively used in various technical applications.Owing to the slip condition for fluid phase intensity of velocity and tem-perature fluctuations has substantial heterogeneity.In this case, the additional force of migration occurs, and consequently the particle drifts to a channel wall [1]- [4].Temperature fluctuations and migration of particle to the wall lead to an additional heat flux.Magnitude of these effects is determined by dynamic and thermal relaxation times of the particle and integral time scale of turbulence.In theoretical studies of two-phase turbulent flows migration force is included on the basis of the probability density function [2]- [4].The main purpose of the present work is to investigate the physical mechanism of inertial particles drifting in the inhomogeneous turbulence.We use a simple model in which the fluctuations of velocity and temperature of the liquid phase are approximated by periodic oscillations.The amplitude of the oscillations is inhomogeneous in space.For separation of a particle motion into fast and slow components is involved the method of averaging developed by Krylov and Bogolyubov [5]- [7].Closed system of ordinary differential equations for particle trajectories and temperature averaged over the oscillation period was derived.It is shown that migration force is directed toward reducing the amplitude of the velocity fluctuations of liquid phase.Effect of particles drift reaches a maximum value for particles with dynamic relaxation time compared with the oscillation period.
For a linear dependence of amplitude of velocity and temperature oscillation on the distance from a wall, the analytical solution was obtained.Combined influence of gravity and migration forces on a particle trajectory was studied.The results of numerical integration of complete system of differential equations of dynamics and temperature of particles are in satisfactory agreement with the calculations obtained by analytical formulas.

Basic Equations
For the purpose of compactness of the presentation and elucidate the physical meaning of obtained results, we consider one-dimensional case.We study oscillating field of fluid velocity in the direction normal to the solid wall (see Figure 1).Coordinate axis x is directed along the normal to the wall.Fluid velocity ( ) , U x t and gravity vector g are parallel to the normal.Differential equations for a particle velocity ( ) Here ( ) Here , where Pr is Prandtl number of the fluid.It can be seen, that the distinction in the thermo-physical properties of liquid and material of particle leads to different values of thermal and dynamic relaxation times.
The initial conditions for the system of Equations ( 1) and ( 2) have the form ( ) ( ) ( ) Dependences of velocity and temperature of fluid on the distance from the wall (see Figure 1) are defined as The main purpose of this manuscript is to show that in an inhomogeneous, rapidly oscillating non-isothermal fluid flow appear migration force and an additional heat flux connected with inertia of the particles.

The Averaging Method
Equation ( 1) is reduced to the single equation for displacement of a particle To solve the Equation ( 5) we involved the Krylov-Bogolyubov method of averaging [4]- [6].Displacement and temperature of particles composed of slow averaged components ( ) and fast oscillations of displacement ( ) By averaging over the oscillation period 2π T ω = we obtain the following expressions Fluctuations of particle coordinate is substantially lesser than the averaged displacement Functions representing the amplitude of velocity and temperature of carrier fluid (4) along the particle trajectory are written in the form of an expansion in the sense of small oscillations of the particle coordinate ( With these formulas and the equation for temperature (2) and coordinate (5) of the particle take the form Fluctuations of velocity and temperature of a particle represented as Correctness of the above presentations are illustrates in Appendix.Period of oscillation T is substantially less than the time scale of variation of averaged parameters for the particle.Time derivatives of velocity and temperature fluctuations of the particle have the form Upon substituting these expressions into the Equation ( 6) we write down For particle temperature (7) following equation is rewritten  ( )

( )
cos t ω and integrated over the period of oscillation T .As a result, from Equation (9) we obtain expressions for the coefficients ( ) ( ) , 1 1 From Equation (10) we obtain expressions for the coefficients ( ) ( ) ( ) .
After averaging the Equations ( 9) and (10) over the period of oscillation we obtain the following equations for averaged velocity and temperature of the particle ( ) From the Equations ( 11) and ( 12) is seen that in homogeneity in amplitude of velocity and temperature fluctuations of carrier fluid leads to migration force and an additional heat transfer for particles.The drift velocity in Equation ( 11) is proportional to derivative of the square of amplitude of fluid velocity fluctuations and directed toward reducing intensity of fluctuations ( ) 12) one can notice that at inhomogeneous temperature of fluid on the particle acts additional cooling heat flux for ( ) It follows from the Equation (11) that, for particles of low inertia V T τ  additional migration force and additional heat flux disappears.Value of migration force and additional heat flux is also reduced for particles with high inertia V T τ  .Migration force is maximal for particles with a relaxation time 1 Actual coordinate and velocity of particles are calculated based on the following algebraic equations Actual temperature of particle is From the formulas (3) it is follow the initial conditions for the averaged velocity, temperature and particle coordinate

Linear Dependence of Velocity and Temperature Amplitudes
Consider the analytical solution for the case of linear dependence of amplitude of fluctuations of velocity and temperature of liquid with the distance from the wall (see Figure 1)

Solution of the Equations for a Particle Velocity and Coordinate
In considered case the equation for averaged particle coordinate has the following form Solution of the last equation can be written as Expression for actual coordinate of particle follows from (13) In this considered special case equation for averaged particle velocity has the form Solution of this equation is e e e e 1 e 1 1 Initial values for the averaged coordinate and velocity of a particle are Actual velocity of a particle follows from expression (

Solution of Equation for a Particle Temperature
Equation for averaged temperature of a particle follows from Equation (12) and has the form Initial value of average temperature of a particle is equal Solution of Equation ( 16) represented By substituting in the above integral expression for the averaged particle displacement ( ) X s , we find out the following formula ( ) ( ) ( ) 1 e e e e 0 e 1 e 2 1 1 1 Formula for actual temperature of a particle follows from expression (15) and takes the form

Calculation Results
We give comparison of the results of analytical research with the data of direct numerical calculation of the system of Equations ( 1) and (2).As dimensionless variables we used Stokes number St V T τ = , dimensionless temporary variable t t T * = and dimensionless coordinates . Dimensionless gravity acceleration is g Tg U * =  and dimensionless thermal relaxation time is 2 illustrates a satisfactory agreement between results of calculations using analytical formulas from the previous paragraph and numerical integration the ordinary differential Equations ( 1) and (2).Migration of particles towards the wall is a result of its inertia.One can see that particle drifts to the wall.This effect can be explained as follows.On the one hand, a particle leaves the region with high level of velocity fluctuations and enters in a region with low level.On the other hand, viscous resistance leads to loss of kinetic energy of the particle.Particle due to its inertia cannot return to the region with a higher level of velocity fluctuations.
The Figure 3 illustrates the effect of preferential concentration of particles at some distance from the wall when the gravity force is directed along the normal to the wall (see Figure 1).Unlike the case ( 0 g * = ), the particle periodically fluctuates at a certain distance from the wall.
The Figure 4 presents a satisfactory agreement between the results of numerical integration of differential Equations ( 1) and (2) and calculations using analytical formulas for the particle temperature.

Conclusions
Simple model of a particle motion and heat transfer in inhomogeneous field of fluctuating velocity and temperature of fluid was suggested.The proposed model reflects the main features of the turbulent flow near the wall.For separation of the fast and slow components of velocity and temperature fluctuations of a particle averaging method of Krylov-Bogolyubov was used.
Physical interpretation of migration force and additional feat flux in inhomogeneous rapidly oscillating velocity and temperature of fluid phase is established.It is shown that the force of migration and additional heat flux disappear for low inertia particles and decrease with increasing particles Stokes number.
For a linear dependence of the amplitude of velocity and temperature fluctuations on the distance to the wall analytical solutions was found.The results of numerical solution of the system of equations of motion and heat transfer of a particle are in satisfactory agreement with the calculations by analytical formulas.
The possibility of formation regions with preferential concentration of particles in a gravitational field and in inhomogeneous velocity fluctuations is illustrated.

Figure 1 .
Figure 1.Sketch of the problem.
Figure represent the alteration of temperature with time (a) and change the particle temperature along its trajectory (b).The particle drift toward the wall and its temperature decreases.