Optimalization of Design Parameters of Experimental Installation Concerning Preparation of Liquid Feed Mixtures

: The article describes the initial conditions for the development of universal mechanization means for the process of mixing dry and liquid components. The essence of the method is to study the motion of a particle with different constructive and physical properties of the medium. The mathematical model of particle motion is based on theoretical mechanics and hydraulics. In this case, the main purpose of the study is to ﬁnd the optimal design parameters for the installation. At the beginning, a theoretical analysis of the installation was carried out using the methods of classical mechanics and hydraulics. Experimental studies were carried out in several stages. At the beginning, one-factor experiments were conducted, followed by allocating the main factors and determining their interaction. Then, using the methods of planning the experiment, we obtained the regression equations and further optimized the parameters to summarize the main ﬁndings of the article. Modern installations should have versatility in any technological line; for example, an installation is presented that can not only mix, but further transport the mixture like a conventional pump, while providing a dosing device that is necessary for the feeding of dry components. Theoretical studies have been carried out in which the design of the impeller is substantiated at various speeds. Experimental studies to determine the design parameters of the installation are in continuous operation. The degree of homogeneity was Θ = 74%, with β 2 = 80 . . . 100 ◦ and β st = 65 . . . 102 ◦ , while the value of the consumption of electrical energy is equal to Eel = 0.265 . . . 0.28 kWh/t.


Introduction
For humans, as well as for society as a whole, nature is a living environment and necessary for the existence of the population. From the perspective of the mechanization of agricultural production, such machines and mechanisms should be created that satisfy the following requirement: economic, health-improving, hygienic, aesthetic, educational and scientific. With the development of farms in western countries-with well-established agricultural production, expectations have been increased in the solution of domestic food problems. It has been confirmed that very few small-scale energy resources for mechanization have been implemented in agriculture. However, the domestic industry can provide an almost complete mechanization of these processes and transform them into a complex one. It is also known that without solving this problem, that is, without a high level of mechanization of production processes, it is impossible to achieve high efficiency in the functioning of farms. The domestic industry has now developed and put into production various energy resources for small-scale mechanization. Their quantitative and qualitative composition can be chosen by the farmer, depending on natural conditions,

Materials and Methods
The installation is proposed as a universal device that can not only mix, but further transports the mixture like a conventional pump while providing a dosing device, which is necessary for supplying dry components. The installation diagram is shown in Figure 1. The basis of the creation was a conventional vane pump, but the radial arrangement of the blades was chosen as a universal impeller; that is, the blades were welded onto the disk. At the same time, the minimum number of blades Z = 6 pcs was observed, as in conventional turbine mixers, but a second row of blades was also installed, with the same number of blades. On the periphery of the working chamber, there is a disc on which blades are located, which help to mix the component. An opening was made in the impeller for the installation of a feeding device in the form of an auger that moves dry material from the loading chamber.
Processes 2021, 9, x FOR PEER REVIEW 2 of 14 various energy resources for small-scale mechanization. Their quantitative and qualitative composition can be chosen by the farmer, depending on natural conditions, specialization, structure of sown areas and the level of organization of production [1,2]. It is worth noting that the existing means of mechanization, namely, those used for the preparation of liquid feed mixtures, are not universal because using several devices to implement the technological process together is problematic.

Materials and Methods
The installation is proposed as a universal device that can not only mix, but further transports the mixture like a conventional pump while providing a dosing device, which is necessary for supplying dry components. The installation diagram is shown in Figure  1. The basis of the creation was a conventional vane pump, but the radial arrangement of the blades was chosen as a universal impeller; that is, the blades were welded onto the disk. At the same time, the minimum number of blades Z = 6 pcs was observed, as in conventional turbine mixers, but a second row of blades was also installed, with the same number of blades. On the periphery of the working chamber, there is a disc on which blades are located, which help to mix the component. An opening was made in the impeller for the installation of a feeding device in the form of an auger that moves dry material from the loading chamber. The theoretical substantiation of such design should be considered, as the movement of a particle along the impeller depends on the different rotational speeds and viscosity of the medium, which changes with temperature. In this case, the temperature will be in the range 20-40 °C, and the maximum value here is dictated by zootechnical requirements for the preparation of liquid feed mixtures.
The calculation and modeling of the mixing process has a certain specificity, considering the properties of the medium, as well as the design of the mechanization device itself. The developed experimental setup [3,4] combines three functions-a batcher, a pump and a mixer-while its versatility allows it to be used in future, with applications beyond agriculture systems. The design of the impeller of the installation allows for a separate supply of dry and liquid components to the working chamber, where they meet at the periphery with fixed blades and the place where the mixing process is provided. When the material moves along the helical channel into the window on the cover disk, where it hits the impeller blade, the model of particle motion is considered a centrifugal field ( Figure 2) [5,6,7]. The theoretical substantiation of such design should be considered, as the movement of a particle along the impeller depends on the different rotational speeds and viscosity of the medium, which changes with temperature. In this case, the temperature will be in the range 20-40 • C, and the maximum value here is dictated by zootechnical requirements for the preparation of liquid feed mixtures.
The calculation and modeling of the mixing process has a certain specificity, considering the properties of the medium, as well as the design of the mechanization device itself. The developed experimental setup [3,4] combines three functions-a batcher, a pump and a mixer-while its versatility allows it to be used in future, with applications beyond agriculture systems. The design of the impeller of the installation allows for a separate supply of dry and liquid components to the working chamber, where they meet at the periphery with fixed blades and the place where the mixing process is provided. When the material moves along the helical channel into the window on the cover disk, where it hits the impeller blade, the model of particle motion is considered a centrifugal field ( Figure 2) [5][6][7]. The composition of the sum of projections of all external forces that act on the particle material, in vector form, is as follows: where: m-particle mass, kg;  -particle acceleration; mg -gravity, Н; where: ρl-fluid density, kg/m 3 ; V-body volume, m 3 . We transform the Archimedes force, accounting for the fact that there is a rigid body in the system. For this, we divide Equation (2) by the mass m = ρV and we obtain: where: ρ-body density, kg/m 3 .
In projections onto the Cartesian coordinate axes x and y, we obtain the following system of equations: The particle velocity can be represented in the form of system equations: The composition of the sum of projections of all external forces that act on the particle material, in vector form, is as follows: where: m-particle mass, kg; where: ρ l -fluid density, kg/m 3 ; V-body volume, m 3 . We transform the Archimedes force, accounting for the fact that there is a rigid body in the system. For this, we divide Equation (2) by the mass m = ρV and we obtain: where: ρ-body density, kg/m 3 . In projections onto the Cartesian coordinate axes x and y, we obtain the following system of equations: m ..
The particle velocity can be represented in the form of system equations: Based on the obtained mathematical model, we will construct the dependence of the particle motion at different viscosity coefficients (Figures 3-6). In this case, the dynamic viscosity takes the different temperatures of the liquid into account.
Based on the obtained mathematical model, we will construct the dependence of the particle motion at different viscosity coefficients (Figures 3-6). In this case, the dynamic viscosity takes the different temperatures of the liquid into account.      From the obtained dependences, we see that, with an increase in resistance, the trajectory of the particle along the x axis will increase; therefore, it will take longer to go in the radial direction. As the size of the working chamber is structurally limited by 75 mm, a rotation frequency of n = 1500 min −1 will be more effective since the particle will travel a shorter distance in the radial direction. If we consider all frequencies at one resistance ( Figure 6), then it can be seen that with a viscosity of k = 1.002 Pas (t = 20 °C) and a rotation frequency of n = 1500 min −1 , the particle fits into the design parameters of the working chamber. However, this only applies at the radial arrangement of the blades. In this case, the nature of the curve shows that, with such vortex formations, the mixing process will obtain a better result. To check the theoretical data, we will test the installation. In this case, it is necessary to restrict the range of rotation frequencies from n = 750 min −1 to n = 1750 min −1 , since it can be seen that, at large values, the particle will not interact with stationary blades. The setup was tested on a specially designed stand (Figure 7), which allows samples to be taken during the test. According to the flow through scheme, the installation works as follows: before starting work, valve 11 is closed and valves 9 and 10 are opened. After starting the electric motor, valve 6 is opened and dry components are poured into loading chamber 7. The resulting mixture enters tank 3 [6]. The quality of the mixture that is obtained by mixing the components is determined using the degree of homogeneity, which characterizes the completeness of the process as a whole [7]: From the obtained dependences, we see that, with an increase in resistance, the trajectory of the particle along the x axis will increase; therefore, it will take longer to go in the radial direction. As the size of the working chamber is structurally limited by 75 mm, a rotation frequency of n = 1500 min −1 will be more effective since the particle will travel a shorter distance in the radial direction. If we consider all frequencies at one resistance ( Figure 6), then it can be seen that with a viscosity of k = 1.002 Pas (t = 20 • C) and a rotation frequency of n = 1500 min −1 , the particle fits into the design parameters of the working chamber. However, this only applies at the radial arrangement of the blades. In this case, the nature of the curve shows that, with such vortex formations, the mixing process will obtain a better result. To check the theoretical data, we will test the installation. In this case, it is necessary to restrict the range of rotation frequencies from n = 750 min −1 to n = 1750 min −1 , since it can be seen that, at large values, the particle will not interact with stationary blades. The setup was tested on a specially designed stand (Figure 7), which allows samples to be taken during the test. According to the flow through scheme, the installation works as follows: before starting work, valve 11 is closed and valves 9 and 10 are opened. After starting the electric motor, valve 6 is opened and dry components are poured into loading chamber 7. The resulting mixture enters tank 3 [6]. The quality of the mixture that is obtained by mixing the components is determined using the degree of homogeneity, which characterizes the completeness of the process as a whole [7]: where: Φ 0 -normalized Laplace function. To assess the quality of the mixing plant, we use one of the indicators, such as specific energy consumption per unit mass of the finished product, determined by the formula (kWh/t): where: Р-power consumed by the installation, kW; Qh-hourly productivity of the installation, m 3 /h or t/h. Equation (8) is universal, but it does not consider the nature of the obtained product; therefore, it is better to introduce an indicator-such as the specific energy consumption of electrical energy-which considers the degree of homogeneity of the mixture, which is determined by the formula (kWh/t): The proportion of particles of the control component in the mixture is within the specified limits ± ∆: in this case, 0 < Θ < 1. The limiting case of complete mixing corresponds to the value Θ = 1. The recipe for the ratio of the components was adopted on the basis of zootechnical requirements in the range from 1:8 to 1:10; that is, for 8 L of water 1 kg of milk replacer, the deviation is permissible ∆ = ±20%. Thus, Equation (6) must be transformed, considering the assumption that the denominator will not contain the number 0.9973-which is close to mixing the components in a 1:1 ratio-but there will be an interval from 0.125 to 0.1, which corresponds to a ratio of 1:8 . . . 1:10.
Processes 2021, 9, 2104 7 of 12 To assess the quality of the mixing plant, we use one of the indicators, such as specific energy consumption per unit mass of the finished product, determined by the formula (kWh/t): where: P-power consumed by the installation, kW; Q h -hourly productivity of the installation, m 3 /h or t/h. Equation (8) is universal, but it does not consider the nature of the obtained product; therefore, it is better to introduce an indicator-such as the specific energy consumption of electrical energy-which considers the degree of homogeneity of the mixture, which is determined by the formula (kWh/t): The power consumed by the electric motor, which drives the installation, is determined by the equation [2,5]: where: U p -voltage phase of the network, V; I p -phase current, given that the "star" connection then; I p = I l -linear current, A; cos ϕ-motor power factor.

Results
First, we will investigate the effectiveness of the preparation of the mixture at a fixed speed in the ranges of 750, 1500 and 1750 min −1 , at a temperature of 40 • C in the ratio of components 1:8 . . . 1:10 with continuous mixing of the components. Therefore, the preparation time will be tpr = 0 min. In Figure 8, we present the graphs of the dependence at a speed of rotation of n = 750, 1500 and 1750 min −1 , with a different number of fixed blades Z = 12, 16, 20, 24 and 28 pcs. In this case, the basic installation angle of the impeller blades was β2 = 90 • , and the angle of fixed blades was βst = 90 • . Subsequently, the angles were changed from 30 • to 150 • (Figures 9 and 10). The power consumed by the electric motor, which drives the installation, is determined by the Equation [ where: Up-voltage phase of the network, V; Ip-phase current, given that the "star" connection then; Ip = Il-linear current, А; сosφ-motor power factor.

Results
First, we will investigate the effectiveness of the preparation of the mixture at a fixed speed in the ranges of 750, 1500 and 1750 min −1 , at a temperature of 40 °C in the ratio of components 1:8 ... 1:10 with continuous mixing of the components. Therefore, the preparation time will be tpr = 0 min. In Figure 8, we present the graphs of the dependence at a speed of rotation of n = 750, 1500 and 1750 min −1 , with a different number of fixed blades Z = 12, 16, 20, 24 and 28 pcs. In this case, the basic installation angle of the impeller blades was β2 = 90°, and the angle of fixed blades was βst = 90°. Subsequently, the angles were changed from 30° to 150° (Figures 9 and 10).   As can be seen from the obtained dependences, the degree of homogeneity will have the greatest value at Ѳ = 78 ... 82% with the number of blades at Z = 20 pcs and at speeds from 750 to 1750 min −1 , while the angle of installation of the blades is 30°; that is, the blades are bent back against the rotation of the impeller.
To optimize such parameters, it is necessary to implement the experimental design and find the optimal parameters through the analysis of two-dimensional sections.
We will implement the experiment plan in several stages. The first plan will assess the implementation of the frequency of rotation of the impeller -x -1 , and the number of stationary blades -x -2 , with the basic values of the angles at β2 = 90° and βst = 90°.
The experiment used to determine the optimal parameters 2 is shown in Table 1.   As can be seen from the obtained dependences, the degree of homogeneity will have the greatest value at Ѳ = 78 ... 82% with the number of blades at Z = 20 pcs and at speeds from 750 to 1750 min −1 , while the angle of installation of the blades is 30°; that is, the blades are bent back against the rotation of the impeller.
To optimize such parameters, it is necessary to implement the experimental design and find the optimal parameters through the analysis of two-dimensional sections.
We will implement the experiment plan in several stages. The first plan will assess the implementation of the frequency of rotation of the impeller -x -1 , and the number of stationary blades -x -2 , with the basic values of the angles at β2 = 90° and βst = 90°.
The experiment used to determine the optimal parameters 2 is shown in Table 1.  As can be seen from the obtained dependences, the degree of homogeneity will have the greatest value at Θ = 78 . . . 82% with the number of blades at Z = 20 pcs and at speeds from 750 to 1750 min −1 , while the angle of installation of the blades is 30 • ; that is, the blades are bent back against the rotation of the impeller.
To optimize such parameters, it is necessary to implement the experimental design and find the optimal parameters through the analysis of two-dimensional sections.
We will implement the experiment plan in several stages. The first plan will assess the implementation of the frequency of rotation of the impeller -x −1 , and the number of stationary blades -x −2 , with the basic values of the angles at β2 = 90 • and βst = 90 • .
The experiment used to determine the optimal parameters 2 is shown in Table 1. The results of the experiment are presented in Table 2. The regression equation was obtained, which describes the influence of the rotational speed and the number of blades: Using the software applications Microsoft Office Excel 2007 and Statgraphis Plus 5.0, the estimates of the regression coefficients were calculated, their significance was assessed and the adequacy of the obtained models was checked. These results were used to construct two-dimensional sections of the response surfaces. The calculations carried out to determine the average value of the response and calculated value of the optimization criterion were determined in Microsoft Office Excel 2007. In Equations (11) and (12), was calculated the number of blades (b2 = +0.36) and frequent rotation (b2 = +0.04) for specific energy consumption of electrical energy.
The estimates of the regression coefficients were considered significant with a 95% confidence level when the p-Value [8,9], given in the ANOVA table, does not exceed 0.05.
Analyzing the section (Figure 11), it can be concluded that, with the number of blades set at Z = 20 . . . 24 pcs and the rotation frequency at n = 1250 . . . 1500 min −1 , the maximum value of the degree of homogeneity Θ = 73.7% was reached. The two-dimensional section shows that the smallest consumption of electrical energy was achieved with a combination of factors Z = 26 . . . 38 • C, rotation frequency n = 1250 . . . 1500 min −1 and the value was Eel = 0.25 . . . 0.265 kWh/t.
shows that the smallest consumption of electrical energy was achieved with a combination of factors Z = 26 ... 38 °C, rotation frequency n = 1250 ... 1500 min −1 and the value was Eel = 0.25 ... 0.265 kWh/t. For future tests, we will choose a rotational speed of n = 1500 min −1 , and the number of blades will be fixed at Z = 24 pcs.
The second plan will evaluate the implementation of the values of the angles β2о and βst.
The experiment used to determine the optimal parameters 2 is shown in Table 3 and Table 4.  For future tests, we will choose a rotational speed of n = 1500 min −1 , and the number of blades will be fixed at Z = 24 pcs.
The second plan will evaluate the implementation of the values of the angles β2o and βst.
Using the software applications Microsoft Office Excel 2007 and Statgraphis Plus 5.0, the regression coefficient estimates were calculated, their significance was assessed and the adequacy of the obtained models was checked. These results were used to construct twodimensional sections of the response surfaces. The calculations carried out to determine the average value of the response and the calculated value of the optimization criterion were determined in Microsoft Office Excel 2007. In Equations (13) and (14), we calculated the angle of the impeller blades (b2 = +1. 41) and the same angle (b2 = +0.04) for specific energy consumption. The estimates of the regression coefficients were considered significant with a 95% confidence level when the p-Value [8,9], given in the ANOVA table, did not exceed 0.05. Table 4. Results of research according to plan to determine the optimal parameters 2.

Factors
Optimization Criterion

Discussion
Analyzing the section (Figure 12), we can conclude that, at impeller angles in the range of 80 . . . 120 degrees and fixed blade angles ranging from 60 to 102 degrees, the maximum value of the degree of homogeneity Θ = 74.0% was reached.

Conclusions
1. Modern installations should have versatility in any technological line; for example, the installation, which can not only mix, but further transports the mixture like a conventional pump while providing a dosing device that is necessary to supply dry components. 2. When conducting theoretical studies concerning design of the impeller [12,13] we can The two-dimensional section shows that the smallest consumption of electrical energy was achieved with a combination of factors β2 = 50 . . . 80 • , βst = 48 . . . 110 • and the value was Eel = 0.25 kW h/t. Considering that the homogeneity was 74%, we will then choose