INTRODUCTIONA significant proportion of food products have viscous and quasi-viscous properties. Viscous media require new engineering knowledge on the development of their production technology and techniques, taking into account the rheological properties of processed raw materials and finished products.The data on the patterns of changes in rheological characteristics make it possible to affect the quality and structure of products both by adding supplements and by regulating the treating modes of a product. The method of intensifying heat exchange when treating viscous products by mixing the treated medium in apparatus with special blade-type mechanisms is used widely. It is necessary to study this method deeper, especially when heat-treating viscous media with non-Newtonian properties. Mixing viscous materials requires a lot of energy. Therefore, to reduce the power used for mixing, it is necessary to minimize the speed of rotation, taking into account the rheological characteristics of the medium treated and ensuring the required treating mode.In many cases, it may be reasonable to use mixing devices for such apparatus that come into contact with the heat exchange surface in order to clean a product layer. This method was studied by L.K. Nikolaev,B.L. Nikolaev, and a number of other authors [1, 2, 3, 4, 5]. Nevertheless, there are a lot of aspects of operation of apparatus with cleaning devices that require a deeper study, in particular, when heat-treating viscous media with non-Newtonian properties such as cream, animal and cooking fats, margarine emulsions, meat and fish mincemeat, and other viscous food products. There are also some ideas for changing the design of mixing mechanisms to activate the process that requires development.The statement of the problem of nonsteady heating or cooling of an infinite cylinder is based on the cases when:– axial heat fluxes are negligible as compared to radial ones;– the coefficient of convective transfer of heat from the outer surface α and the ambient temperature tc are constant;– there is a symmetry of the initial temperature distribution θg (r) round the radius of the cylinder;– the internal heat sources are assumed to be absent(q = 0).The assumptions made correspond to the following mathematical model of the process of nonsteady thermal conductivity of cylindrical bodies:  . (1)   However, this theoretical model does not take into account a number of significant practical circumstances that affect the process, such as the curvature of the surface, the frequency of scraper passages, the non-Newtonian nature of ryazhenka [fermented baked milk], viscosity, and a lot of other food products. The solution of the arising engineering problems can be implemented by experimental methods ensuring the extension of the obtained results using the similarity theory and dimensional analysis.The time-varying parameters of a number of technological processes are also important parameters of theoretical description of the processes carried out in a great number of food industry equipment. The starting modes of machines and apparatus can be considered a common example of such nonsteady processes. In many cases, it is necessary to take into account the effect of starting mode characteristics since the absence thereof results in inaccurate measurements, as well as in errors in the choice of operating modes of equipment. Thus, the estimation of the features of the starting mode is necessary [12, 13]. The studies of thermal and hydro-dynamic processes in apparatus with cleaning mixing devices are presented in the works by Begachev V.I., Braginskiy L.N., Glukhov V.P., Pavlushenko I.S., and Pavlov M.G. [7]; Konviser I.A. [8]; Konsetov V.V., Kudryavitskiy F.M., and Novichkov A.N. [9], Pavlushenko I.S. and Gluz M.D. [12], Abichandani H. [19]; de Goede R. and de Jong E. J. [20]; Dumont E., Fayolle F., and Legrand J. [21], et al.The studies of the mentioned above authors were carried out using the equipment that consisted of an apparatus with mixing devices that differed from the experimental plant used in our paper. Part of the considered results of the studies can be regarded only as a first approximation. STUDY OBJECTS AND METHODSTo successfully design and calculate the considered type of mixer, it is necessary to have the formulae for calculating the energy consumed in the flow apparatus with mixing devices obtained from the experimental data. Some data were obtained using directly an industrial mixer equipped with a set of sensors. To study the process of interaction of a rubber scraper with the cooling surface and the working body deeper, a special experimental laboratory plant based on a modified rotational rheometer was designed and tested [7, 8, 9]. A series of experimental studies were carried out using this modified rheometer. The manufactured scrapers were installed on the rotor cylinder of the rheometer. Alternative stator cylinders were used fo study for various volumes of operating area. When carrying out the studies, a set of methods was used that made it possible to determine the rotation speed of the rotor with a high accuracy, as well as to determine the dynamic pattern of the speed to estimate the starting modes.Experimental plant no. 1, a scraped-surface heat exchanger, is shown in Fig. 1, 2 and 3 where:1 – cooling jacket, 2 – working cylinder, 3 – inflator,4 – motor-reducer, 5 – product supply branch,6 – coolant outlet, 7 – elastic tape, 8 – metal perforated slat, 9 – hole in the blade, 10 – gap between the blade sections, 11 – crosspiece holder, 12 – shaft,13 – product outlet, 14 – coolant supply branch,15 – cleaning blade, 16 – inflator blade, 17 – elastic tape bolt, 18 – crosspiece hub bolt, 19 – crosspiece hub.  B B   Fig. 1. Experimental plant no. 1 for studying thermal and hydrodynamic processes.   Fig. 2. Inflator. B – B   Fig. 3. Cleaning devices of the experimental plant. The experimental plant is similar to a T1-OM-2T heat exchanger in design, but at the same time it is distinguished by the fact that in order to increase the productivity of the apparatus, the mixing device located in the working cylinder has some perforated cleaning-type blades installed at an angle of 90° to the surface of the working cylinder. Each of the blades is made of two metal plates, between which an elastic tape the upper edge of which rises above the plates by no more than 15–10 mm is clamped. The distance between the peripheral edges of diametrically located blades is greater than the internal diameter of the working cylinder Dc by 0.5–0.1 mm. The blades have holes of at least 0.05–0.1 Dc and are divided into sections with a length of no more than 3 Dc, with each successive blade displaced in relation to the previous blade in the longitudinal axial direction by 0.2 Dc, (the blade width is 0.2–0.4 Dc). There is also an inflator in the apparatus the diameter of which for the peripheral edges of the blades is equal to 0.9 Dc. This combination of features and their distinctive relationship make it possible to increase the productivity of the apparatus.The new design allowed us to conduct further experimental studies of heat exchange and consumed energy. As a model medium, 2.5% fat ryazhenka was used. The studies were performed within the range of the Reynolds criterion from 90.7 to 6380, while the effective viscosity values were found within the range from 0.03 to 0.14 Pa·s. The studies were carried out in the rheology laboratory of the Mega-department of biotechnologies and low-temperature systems of the ITMO University.The Table 1 presents the results of the studies of rheological processes. RESULTS AND DISCUSSIONStudies of the hydrodynamic processes in an apparatus with scraping mixing devices were conducted.The obtained data, when there is a change in the product temperature and there are different values of a shear rate gradient, show that the effective viscosity of the studied product depends significantly on the value of a product velocity and temperature gradient.When treating ryazhenka in a heat exchanger with cleaning devices at values of the Reynolds criterion below 2400, there is a laminar motion of the medium in the apparatus. With a further increase in the Reynolds criterion, there is a transitional flow range.Taking into account the previous large-scale studies carried out by a number of authors [1, 2, 5, 14, 15], to apply the obtained data to a wider range of equipment and treated media, a criterial heat exchange equation for Newtonian and non-Newtonian fluids was proposed taking into account the similarity theory: Table 1. Change in the rheological characteristics of 2.5% fat ryazhenka during heat treatment Producttemperaturetpr, °СTemperatureof the heat exchange surface tat, °СProductdensityρ ρ, kg/m3Scraper agitator rotation speedn, r/sConsumed energy N, WEffective viscosityμ μef, Pa∙sReynolds criterion ReHeating ryazhenka18.424.81031.10.7831.50.14090.722.427.21028.11.2502.90.105192.523.227.21027.91.8004.70.086338.324.527.81027.42.6178.00.070604.029.933.41024.83.3678.80.0481130.031.234.21024.34.30014.00.0421649.032.735.31023.75.20017.10.0362325.0Cooling ryazhenka31.926.61023.91.3672.30.071310.030.326.11024.71.9334.10.064486.729.625.51024.82.2175.10.062576.228.024.81025.82.5676.30.060690.126.723.91026.23.36710.80.057953.224.522.11027.34.36714.90.0491440.022.020.41028.48.50041.70.0393158.020.519.51028.911.83160.80.0306380.0      ,                 (2)   where:   Nu=α∙lλα∙lα∙lλλis the Nusselt criterion; Re=ω∙l∙ρμω∙l∙ρω∙l∙ρμμ is the Reynolds criterion; ω=π∙D∙n  is the speed of the cutting edge of the scraper, m/s; Pr=c∙μλc∙μc∙μλλ  is the Prandtl criterion;  is a characteristic dimension which is a distance between the outer edges of the cleaning devices, m;D is the internal diameter of the apparatus, m;z is the number of blades of the mixing device, pcs;μ μ and μстμμстстμat are the effective viscosity of the product at the average product temperature and at the wall temperature, respectively, Pa∙s;B is the empirical coefficient;a and b are the dimensionless coefficients determined in the course of the experiments.The use of l in the calculations, which is a distance between the edges of the blades of the mixing device, allows us to take into account the effect of the number of scrapers and the inner diameter of the working cylinder. When performing the theoretical description of the operation of a great number of apparatus and measuring equipment, an important problem is the technological processes the parameters of which can vary with time. Equipment starting modes are one of the common examples of such nonsteady processes. Variable parameters lead to deviations in the hydrodynamics of a medium, can affect the heat exchange in heat exchange equipment, they also lead to significant changes in a number of other processes carried out in food processing equipment. In many cases, ignoring the effect of starting mode parameters results in measurement errors and the selection of incorrect operating modes for machines and apparatus.To describe theoretically the operation process of this type of apparatus, it is possible to present the hydrodynamic pattern of a flow, which is a nonsteady flow in a flat gap, in a simplified form:  (3)    (4)     where νx, νy, νz are the projections of flow velocities; x, y, z are the coordinates; r is density; p is the pressure; gx is the projection of the acceleration of gravity; τ are the components of the stress tensor.Since there is a shear flow of a viscous fluid in the experimental apparatus, it is necessary to take into account an equation that relates the rate of shear deformation and stress.  .                              (5)   To simulate the motion of an incompressible viscous fluid between the plates of the apparatus, let us take the lower plate as a stable one and the upper plate is the one that constant voltage is applied to. TheEqs. 3–5 take the following form in the partial derivatives:  .                           (6)   The boundary and initial conditions are assumed to be equal to: (7)   . Let us transform the Eq. 6 in the partial derivatives to the equations represented in the ordinary derivatives:  .                       (8) Then: (9) .  Therefore:  (10) ,   (11) .  Hence: (12) ,   (13) .  In the ordinary derivatives, the equations take the form: (14) ,   (15) .  Let us integrate the Eqs. 14 and 15 with allowance for the boundary and initial conditions (7) that determine the constants C1, C2 and C3: ,               (16)   ,                     (17)  (18) .  The solution to the Eqs. 16–18 allows us to carry out the mathematical study of the development of the velocity profile of a flow in the space between the cylinders and to calculate the starting mode time of the apparatus:  (19)    . where τн is tangential stresses  constant voltage, η is the dynamic viscosity coefficient.The analysis of the data shows that when studying the initial phase of the apparatus with mixing devices, it is necessary to take into account the starting period.Considering a number of assumptions in deriving the formula (19) for the experimental studies of starting modes, Plant no. 2 shown in Fig. 4 and 5 was used. It is a modified Volarovich viscometer [10, 11] equipped with devices for measuring the rotor speed.In addition, the modification included the use of alternative external cylinders and the installation of scraping action simulating devices on the rotor cylinder. The experiments consisted in measuring the dependence of rotor speed on the time from the starting moment, while a constant torque was having an impact on the rotor. A series of experiments were carried out with external cylinders with a radius of 19 mm, 34 mm, and 72 mm with a radius of the inner cylinder of16 mm. As a model medium, 2.5% fat ryazhenka as a sample of a food product with non-Newtonian properties and 99.5% glycerin which is a medium with Newtonian properties were used. The use of glycerin is caused in particular by its well-known physical characteristics.As a result of processing the experimental data obtained using Experimental plant no. 1, the dimensionless coefficients a, b, c for the criterion equation (2) were obtained.The exponential factor of the Reynolds criterion was determined as a result of processing the experimental data obtained in the logarithmic coordinates:  ,                (20)   The study was carried out for the groupsof experiments the Prandtl criterion in which was:Pr1 = 388, Pr2 = 2462.The graphical-analytical processing of the obtained dependences presented in Fig. 6 showed that the both groups of experiments had the same slope to the abscissa axis. The exponent a of the formula (2) was determined as the slope of the graph to the axis of abscissae and was equal to a = 0.59. The value of the exponent indicated the significant effect of the Reynolds criterion on heat exchange.     Fig. 4. Experimental plant No. 2 for the study of starting modes in apparatus with cleaning devices:1 – digital non-contact tachometer; 2 – block;3 – weight; 4 – tachometer mark; 5 – pulley;6 – scraper; 7 – stator capacity; 8 – rotor cylinder.   1 5 6 2 9 3 4 8 7   Fig. 5. Diagram of the experimental plant basedon a modified rheometer: 1 – rotor cylinder; 2 – stator capacity; 3, 7 – weight; 4 – block; 5 – pulley;6 – electronic IT 5-ChMTermit contactless tachometer, 8 – scraper, 9 – holes in the blades of the working bodies.    Fig. 6. Dependence of the Nusselt criteria and the Reynolds criterion in an apparatus with mixingdevices – the determination of the exponent of degree of the Reynolds criterion. Dependence of the rotational speedon time when treating glycerin Dependence of the rotational speedon time when treating 2.5% ryazhenka Ω, c-1         Ω, c-1                                                          t, c                                                                                                       t, c Fig. 7. Dependence of the rotor speed on time from the starting moment for 2.5% fat ryazhenka and for 99.5% glycerin.  Similarly, the effect of the Prandtl criterion on heat exchange was determined. The obtained graphical dependence uses the Reynolds criterion coefficient obtained at the previous stage.  (21) .  As a result of graphical-analytical processing of the experimental data, the value b of the exponential factor for the Prandtl criterion for the formula (2) was found to be equal to 0.37.To obtain the coefficient B of the formula (5), the data are processed mathematically on the basis of the dependence (2):  .         (22)   As a result of the made experiments and mathematical processing of the experimental data, the criterial equation (2) takes the form:  (23) .  As a result of processing the experimental data obtained using Experimental plant no. 2, the starting mode parameters of the apparatus the character of which is close to an exponential dependence were obtained, they can be described in the form of a formula (24), the coefficients a and b in this formula have been determined by processing the data using CurveExpert (Fig. 7).  (24)   where Ω is the rotor speed, s-1;a and b are the empirical coefficients:a = 0.626∙10-1; b = 5.32 – the correlation coefficient for 2.5% fat ryazhenka is: R = 0.964; the rms deviationδ = 0.0607.For 99.5% glycerin, a = 0.814; b = 5.89; the correlation coefficient R = 0.957; the rms deviationδ = 0.0425.  CONCLUSIONThe effect of the Reynolds and Prandtl criteria on the heat exchange in an apparatus with modified cleaning devices has been determined. For the scraper devices with holes, there was an increase in the Reynolds criterion by 3–5% compared to the measurements made for the devices without holes.In order to apply the obtained data to a wider range of equipment and treated media, a criterial heat exchange equation for Newtonian and non-Newtonian fluids taking into account similarity simplexes was proposed that takes into consideration the similarity theory. The effect of the Reynolds and Prandtl criteria on the heat exchange in the apparatus were determined.The lawfulness of describing the characteristics of the starting mode were established using an exponential dependence of the form of Ω = a (1 – e-bt).The experimentally expected theoretical predominance of the effect of medium density on the starting mode compared to the effect of viscosity properties were established. The value of the coefficient of the exponent b differed by 8–15% for the measurements with ryazhenka and glycerin under other equal conditions. CONFLICT OF INTERESTThe authors declare that there is no conflict of interest related to this article.