Investigation of the Process of Simple Distillation in Irrigated Pipe Elements

In modern chemical and oil refining complexes, separation processes are among the most popular and energy-intensive. Installations for their implementation should be equipped with nodes for creating vapor (evaporators) and liquid (deflegmators) irrigation. Evaporators of any type (film, thermosiphon, gas lift, cubic) belong to this class of devices. For example, in cubic evaporators, the gas flow is completely formed from flux bubbles that originate on the heat-conducting surface and float in the volume of the cubic liquid located in the apparatus. Due to the accompanying mass exchange, the bubbles are enriched with volatile components during ascent and noticeably increase in volume, and the growth of the bubble is determined, among other things, by the total flow. At the same time, in real bubbling-type equipment, the total mass transfer surface exceeds the cross-section of the device itself by more than two orders of magnitude. Thus, according to, the ratio of the internal cross-sectional area of the apparatus to the developed mass transfer surface is 0.0015–0.002. Based on the analysis of the integral equation of the diffusion boundary layer, it is shown that the presence of a resultant flow of substance through the phase interface (non-equimolarity of the process) in a two-phase gas (vapor)–liquid system leads to the transformation of the structure of the traditional mass transfer equation itself. The use of a new structure obtained for both binary and multicomponent mixtures makes it possible to significantly simplify the approach to the description and generalization of arbitrary mass transfer processes. The innovativeness of the proposed approach lies in its universality for non-equimolar processes. This simplifies the creation of models of any mass transfer devices and entire production lines. In addition, the proposed approach is a good auxiliary tool for various researchers and experimenters. It should also be noted that the separation processes of many products of organic origin during heating are characterized by the appearance of undesirable side reactions (thermal decomposition, condensation, polycondensation, formation of harmful impurities, etc.), which occur most intensively in the heating zones. At the same time, the evaporation and distillation units are subject to requirements for the minimum hydraulic resistance of the structure, the maximum separation capacity (efficiency), and the minimum residence time of the product in the apparatus (equivalent to the minimum holding capacity of the structure). It was noted that the specified requirements are most fully met by film-type devices.


Introduction
In view of the obvious influence of the velocity of the total phase flow flowing around the interface on the diffusive mass transfer between the core of the flow and the interface, the description of non-equimolar mass transfer becomes more complicated, starting from the very structure of the mass transfer equation [1][2][3][4]. This circumstance is explained by the transformation of velocity profiles and concentrations of the main flow under the influence of the transverse velocity component, which can occur under the influence of the diffusion stream, and under the influence of the transverse velocity component on the hydrodynamic Processes 2021, 9,2047 2 of 15 characteristics of the boundary layer. As a result, in turbulent streams, mass transfer is characterized not only by molecular diffusion, but also by turbulent diffusion, which in turn depends on the density of the transverse mass flow. Thus, bubbling columns and rectors were considered in [5], and the influence of the hydrodynamics of the process on the mass transfer coefficient, which is a key parameter for evaluating the operation of a gas-liquid contact device, was analyzed. P. Ham used a phenomenological relation that allows calculation in a more rigorous formulation of the problem.
To date, the kinetics of non-equimolar mass transfer is described by several basic methods:

1.
A method for determining the mass transfer coefficients in the vapor and liquid phases, during rectification based on the data of a physical experiment on the profiles of component concentrations along the height of the column, for the equimolar formulation of the problem [6][7][8].

2.
The method of coupled physical and mathematical modeling based on the variational formulation of conservation laws [9]. 3.
The first direction under consideration is based on the determination of mass transfer coefficients in the vapor and liquid phases from experimental data on mass transfer in two-phase systems and the determination of component concentration profiles by column height. The authors [2][3][4] note the need for a more rigorous formulation of the problem, especially when calculating the mass transfer coefficients in the vapor and liquid phases from experimental data. At the same time, difficulties arise when extracting mass transfer coefficients from a physical experiment in a distillation column [6][7][8].
The second direction-"the method of coupled physical and mathematical modeling based on the variational formulation of conservation laws" is based on the relationship of hydrodynamic and mass transfer characteristics of the process. The method consists in representing the processes occurring in the equipment in the form of elementary zones with a certain structure formed by the design of the apparatus and depending on the space-time frame. The mathematical description of each zone is carefully verified during the physical experiment on laboratory models. "The invariance of the structure of the description of the elementary zone to the scale of the apparatus as a whole is postulated" [9]. The third direction is based on an algorithm for generating a vapor potential based on the chemical characteristics of molecules to close the design of devices. For this purpose, the most important parameters of the system under consideration were identified: parameters such as the depth of the potential well and the effective diameter of the potential molecule, as the most approved in thermophysical practice, are determined [10,11]. The interaction within molecules is determined by quantum mechanical parameters: the length and the binding energy [12].
In the work of L. A. Serafimov [13], using the theory of thermodynamics of irreversible processes, it is shown that even for an equimolar process in a strict formulation of the problem, the structure of the equations of multicomponent mass transfer differs very significantly from the phenomenological relations. The calculation was based on multiple regression of logarithmic dimensionless variables and was coupled with experimental data; the deviation was less than 10%.
In the article [14], the authors conducted a numerical analysis of evaporation from the surface of a liquid film of binary composition (a mixture of water and ethylene glycol). The results of the analysis reflect the influence of the input parameters of the gas and liquid film on the temperature and concentration profiles. Thus, it was found that the concentration of ethylene glycol in the liquid phase has a great influence on heat and mass transfer. Moreover, when its concentration in the mixture is less than 40%, more water evaporates than if the film consisted only of pure water under the same conditions.
In [15], A. Stankevich and Ya. Mulein considered the existing methods used to intensify the processes of mass transfer. The authors considered new designs of mass transfer equipment for the chemical and petrochemical industries. New engineering solutions were Processes 2021, 9,2047 3 of 15 also offered for the construction of contact devices for absorbers and rectification columns. The proposals considered "promise impressive improvements in the productivity and efficiency of equipment, as well as noticeably decreasing their linear dimensions" [15].
In [16], mass transfer is described by criterion equations using Nusselt and Sherwood numbers. The authors investigated the features of mass transfer during condensation of vapor-gas mixtures. The experiments were carried out in circular channels filled with moist air. There was a temperature gradient on the channel walls, thereby causing changes in saturation conditions, in addition to condensation or evaporation from the walls. Moreover, the paper clearly demonstrates how the film layer of condensate changes over time, increasing during condensation and thinning during evaporation.
The article [17] provides a solution to the problem of heat and mass transfer under conditions of phase transformations or chemical reaction, which gives an approximate method for solving a system of kinetic equations. Using the example of calcining limestone particles, it is shown that continuous solutions can be obtained for all stages of the process [17]. The proposed method is applicable to solving specific problems of heat and mass transfer under conditions of chemical reactions or structural transformations, which can be useful for theoretical analysis of the kinetics of technological processes.
When calculating mass transfer coefficients, researchers quite often use criterion numbers, and, consequently, similarity theory, which in general, as mentioned above, certainly affects the accuracy of the problem statement during its study. For example, in [18,19], mass transfer was studied during the reaction of a coke solution (C coke + CO 2 = 2CO), which is essentially a special case of a non-equimolar process. However, non-equimolar diffusion in the process of external and internal diffusion was often neglected when studying the reaction of loss of coke solution. Junchen Huang mentions that the solution losses were 25% by mass, which indicates the tangible contribution of non-equimolarity in modeling mass transfer.
A similar approach is used in the work [20], in which a new process for the production of ultrapure gasoline was presented, which solves the problem of loss by loss of octane number (ON) during desulfurization. In the production of pure gasoline, using naphtha obtained by catalytic cracking (FCC), this problem is well known. The authors propose an approach to using a unique approach to separation: regression calculation of LLE results [19].
The analysis of works devoted to non-equimolar transfer of matter shows that relatively reliable calculation methods have been developed only for cases of evaporation and condensation of pure components in an inert medium. It should be noted that in these works, studies were carried out either of the processes of heat and mass exchange or evaporation of a liquid at a temperature below its boiling point, and the influence of the Stefan flow on mass transfer was also studied.
The scheme of interaction of phase flows in a single tube element of the film apparatus is shown in Figure 1. Mass transfer between the phases occurs on the contact surface of the phases between the liquid film flowing down the inner wall of the pipe and the steam flow rising up the pipe. This scheme corresponds to the countercurrent mode of phase interaction. For the direct flow scheme, the picture does not change fundamentally. It should also be taken into account that, for a significant class of devices (film evaporators, distillers), the entire vapor irrigation flow is formed due to partial evaporation of the liquid film, and vice versa-for partial condensers-the entire liquid irrigation flow is formed due to partial condensation of the flux flow. Evaporation or condensation is achieved at the same time due to the external supply (removal) of heat through the pipe wall. At the interface itself, in each section, it is assumed that a thermodynamic equilibrium is achieved between the interacting phase flows [P x = P y ; T x = T y ; y f = y * (x f )]. yf-molar concentrations of gas in the interface of the phases, mol. fraction; xf -molar conce tions of liquid in the interface of the phases, mol. fraction; y∝-molar concentrations of gas in phase (core) center, mol. fraction; x∝ -molar concentrations of liquid in phase (core) center, Fraction.

Mass Transfer in the Film
In this case, in the film column, the intensity of one of the flows will increase the initial (zero) value to the final (maximum) value, and the intensity of the other will progressively decrease. The effect of the convective (total) mass flow on the transfer is maximized in this case. Therefore, when calculating the mass transfer pr it becomes necessary to specify the distribution of the total flow along the interface phases. In relation to evaporators, condensers, and distillers, this condition must b on the basis of appropriate thermal calculations. Another formulation of the problem also be used: an initial task of the distribution of the total flow along the interface, foll by a thermal calculation to determine the appropriate profile of the heat supply (heat along the interface to ensure a given technological condition.
As already noted above, the interface itself is characterized by a complex wave ture, the state of which is determined by the flow regime of the liquid film. Howev most studies, this circumstance is not taken into account, and the surface itself is tak a smooth cylindrical film: F = π − 2 ̅ . When calculating this surface, the av thickness of the film is used. According to [2], the error in determining the interfac resulting from this assumption can reach 15%. At the same time, because in the k generalizations the mass transfer coefficients were also calculated taking into accou same assumption, errors in determining the contact area of the phases are, to a c extent, compensated by a corresponding increase in the calculated mass transfer c cients included in the generalizing equations. Obviously, the mass transfer coefficie this case should be considered as some conditional coefficients related to some condi surface, which is rigidly related to the height of any arbitrary section of the pipe.

Mass Transfer in the Film
In this case, in the film column, the intensity of one of the flows will increase from the initial (zero) value to the final (maximum) value, and the intensity of the other flow will progressively decrease. The effect of the convective (total) mass flow on the mass transfer is maximized in this case. Therefore, when calculating the mass transfer process, it becomes necessary to specify the distribution of the total flow along the interface of the phases. In relation to evaporators, condensers, and distillers, this condition must be met on the basis of appropriate thermal calculations. Another formulation of the problem can also be used: an initial task of the distribution of the total flow along the interface, followed by a thermal calculation to determine the appropriate profile of the heat supply (heat sink) along the interface to ensure a given technological condition.
As already noted above, the interface itself is characterized by a complex wave structure, the state of which is determined by the flow regime of the liquid film. However, in most studies, this circumstance is not taken into account, and the surface itself is taken as a smooth cylindrical film: F = π d − 2δ H. When calculating this surface, the average thickness of the film is used. According to [2], the error in determining the interface area resulting from this assumption can reach 15%. At the same time, because in the known generalizations the mass transfer coefficients were also calculated taking into account the same assumption, errors in determining the contact area of the phases are, to a certain extent, compensated by a corresponding increase in the calculated mass transfer coefficients included in the generalizing equations. Obviously, the mass transfer coefficients in this case should be considered as some conditional coefficients related to some conditional surface, which is rigidly related to the height of any arbitrary section of the pipe.
Most researchers [2] write the mass transfer processes as: Or: As is evident here, the mass transfer equations are special cases of the equations for the transition from a multicomponent mixture to a binary one and for neglecting the total flow, i.e., for an equimolar process [21][22][23]: When taking a linear dependence for the equilibrium line (y f = mx f + b), a transition to the mass transfer equations can be implemented: (6) or where 1/K oy = 1/β y + m/β x (8) and In Equation (6), the transfer coefficients and driving forces are expressed in terms of the corresponding characteristics of the vapor phase, and, in Equation (7), in terms of the characteristics of the liquid phase. The phase resistances to the mass transfer process are determined by the values of the mass transfer coefficients, so Equations (8) and (9) are called the phase resistance additivity equations.
When conducting experiments on the study of mass transfer in devices of any type, the researcher can obtain information only on the mass transfer coefficients. The transition to information about the mass transfer coefficients can then be obtained only on the basis of model representations, for example, from the analysis of Equations (8) and (9). To date, a large number of methods have been developed for decomposing the total mass transfer coefficients into partial mass transfer coefficients [2,24].
Generalization of the experimentally obtained information about the mass transfer coefficients in the vast majority of cases was carried out by researchers using the apparatus of similarity theory. Generalizations were carried out for narrow groups of processes: rectification, absorption, evaporation, etc. In turn, the calculation procedure assumes the use of mass transfer coefficients of only one type-equimolar. In the future, only generalizations for the rectification process will be mainly used. When calculating the processes occurring in multicomponent mixtures, binary mass transfer coefficients are also used for all binary pairs that make up the multicomponent mixture. It is assumed that the binary mixtures are located in a hydrodynamic environment (Re = idem, µ = idem, σ = idem, etc.), corresponding to the conditions of separation of the multicomponent mixture, but differing from each other in their diffusion characteristics.
Certain difficulties arise in this case, especially for devices operating in conditions of pronounced non-equimolarity. Thus, in film evaporators of the countercurrent type, the flow modes of each of the phases, as already noted, vary very widely. Moreover, in some cases, the physically realizable conditions for the existence of film flows may be violated. In these cases, it is necessary to use pipes of variable cross-section according to the height of the device in order to avoid the occurrence of choking or tearing of the film. Accordingly, the calculation procedure should provide for checking the flow modes of the phases and the selection of criteria equations corresponding to specific conditions of phase interaction.
In Equations (11) and (12), E k is the total flow expressed in mass units (kg/(m 2 ·s)), and these expressions themselves are analogs of the Reynolds criterion for the transverse flow of a substance in the pipe. The authors note that the expression (10) is valid only for local pipe sections with specific values of the variables included in it. In [2], an attempt is made to obtain the flow using these criteria and mass transfer equations, which are already characterized by significant complexity for binary mixtures. Therefore, the practical use of this approach is problematic. The authors draw an interesting conclusion that "partial condensation in a wide range of thermal loads has the opposite effect on the kinetics of mass transfer in the vapor and liquid phases." It can be noted that this conclusion directly follows from the structure of the Equations of non-equimolar mass transfer for the vapor and liquid phases (3)-(5) proposed in this paper. Moreover, in this case, the effect of continuous redistribution of phase resistances in the process of non-equimolar mass transfer is easily explained and taken into account for any processes for both binary and multicomponent mixtures.
Most reliable generalizations in the field of the kinetics of mass transfer processes for both vapor and liquid phases have been made for turbulent flow regimes [2][3][4][21][22][23][24][25][26][27]. This is natural, because it is in these areas that real mass transfer devices are operated. At the same time, within the framework of the proposed approach, the field of actual operating modes of the equipment has significantly expanded and covers low-intensity flow modes of one of the phases, up to laminar modes. However, the presence of a convective flow of matter through the interface greatly affects the diffusion transfer. Moreover, under conditions of significant non-equimolarity in low-intensity flow regimes, the diffusion process in the corresponding phase can be completely blocked, and the mass transfer can be completely controlled by the diffusion resistance of another interacting phase, which in this case must be in a turbulent flow mode. These circumstances reduce the loss of calculation accuracy caused by possible errors in the calculation of mass transfer coefficients. Thus, when calculating the mass transfer processes, the following criteria equations were used to calculate the mass transfer coefficients.

•
Laminar flow mode (Re ≤ 2300): which was obtained theoretically and confirmed experimentally [27]. • Turbulent mode (Re > 2300) equation: which was obtained in [26] on the basis of a numerical experiment on a mathematical model, and correlates well with the data of a number of researchers [21][22][23].

For the Liquid Phase
• Laminar wave-free mode (Re G < 12). In this mode, the stability of the film flow and the very possibility of its existence is easily violated. Therefore, in this mode, the mass flow of any component was calculated without taking into account the diffusion flow (only by the convective component of the total flow). • Laminar wave modes (12 < Re G < 320)-V. A. Malyusov's equation was used [28]: which was obtained by the generalization of experimental data on rectification in a film column [29].

Development of an Algorithm for Calculating Mass Transfer in Irrigated Pipe Elements
Based on the developed mathematical description of the processes of non-equimolar mass transfer, a method was developed for calculating the mass transfer process in a film apparatus by integrating the mass transfer equations over the interface. For an infinitesimal surface element, the scheme can be shown as in Figure 2: Processes 2021, 9, x FOR PEER REVIEW

•
The turbulent mode ( < 320) is represented by the equation: which was obtained by the generalization of experimental data on rectificati column [29].

Development of an Algorithm for Calculating Mass Transfer in Irrigated Pipe E
Based on the developed mathematical description of the processes of no mass transfer, a method was developed for calculating the mass transfer proc apparatus by integrating the mass transfer equations over the interface. For a mal surface element, the scheme can be shown as in Figure 2: The system of Equations (17) and (18) allows us to calculate the characteristics of the output flows for the known values of the characteristics of the input phase flows (G n , y i,n , L n , x i,n ), the availability of information about the transfer coefficients B * ,y i,j , B * ,x i,j , and the regularities of the formation of vapor-liquid equilibrium. For numerical integration, the apparatus (pipe) is divided into a number of sections with a length within which all the kinetic coefficients, and hence the mass flows, can be assumed constant. Trial calculations have shown that the length of a separate section can be taken in the range of 100-200 mm. Further reduction in the length does not lead to an increase in the accuracy of the result.
Because the total flow occurs in the apparatus as a result of heat supply (heat removal) through the pipe wall, it is necessary to supplement the system (17) and (18) with heat balance equations or use any additional conditions to determine the temperature of divergent flows. It is most convenient to use the assumption that the contacting phases, the composition of which is determined by the average discharge concentrations for the corresponding cross-section, are under saturation conditions, i.e., their temperature corresponds to the dew point.
The algorithm is designed to solve the problem in a verification statement: the flow rate and composition of the power flow and the geometry of the device (diameter, length, and number of contact tubes) are set as initial data. The separation calculation is carried out for a single tube, i.e., the assumption is made about the uniformity of the distribution of power across all tubes. The calculation is carried out in the following sequence: Based on technological considerations, the distribution of the total flow along the mass transfer surface is set (for example, uniform).
For all sections of the pipe, an initial approximation is built for the costs and compositions of liquid irrigation. When forming an initial approximation of the flow rate, the characteristics of the power flow are used.
Further calculation of the separation process is carried out from the bottom up (from the lower section of the pipe to the upper one).
For each calculated section of the pipe, according to the equations, the thickness of the liquid film is calculated, the conditions for ensuring the integrity of the liquid film are checked, and its flow mode (Re L ) is calculated. For the gas phase, the condition for the presence of a reserve in terms of the choking rate is checked and the mode of its flow (Re G ) is calculated. All the necessary physical and chemical properties (temperature, density, viscosity, diffusion coefficients, etc.) are calculated for the average section of the site. The criterion equations are selected and the mass transfer coefficients are calculated for both phases of the system β * ,y i,j , β * ,x i,j . The joint solution of the mass transfer Equations (23) and (24) provides the mass flows of components formed in the calculated section N i, f . The calculation is reduced to finding the concentrations of phases on the interface that ensure the fulfillment of conditions (25) and (26). When searching for the roots of a solution, for example, the Newton method can be used.
According to the known values of mass flows, the characteristics of the output phase flows are calculated: The previously accepted values of expenses and flow compositions are replaced by the newly calculated ones. Based on the assumption that the average discharge compositions of the interacting phases correspond to the saturation conditions, the temperatures of divergent flows are calculated. The newly calculated characteristics of the phase flows are placed in all the calculated cells instead of the previously accepted ones.
After completing the calculation procedure for all sections of the tube (cells), the condition for achieving a stationary state of the system is checked, which is characterized by the fulfillment of the condition: where ζ i = 10 −4 is the specified accuracy of the solution. In this case, y D,i and x W,I are defined, respectively, as the composition of the vapor leaving the last (upper) cell and the composition of the liquid leaving the first (lower) cell. If condition (31) is not met, the calculation is performed again. When the condition is met, the calculation is completed.
The method used to find a solution for a multi-element apparatus is a kind of method of simple iterations and has received the name "relaxation method" in the literature [24]. It is proved that its application provides reliable and fairly fast convergence. It is also very economical from the point of view of programming and in machine implementation. In this paper, the calculation program for the distillation apparatus was developed using the ANSYS Fluent software and computing complex. Because this is compatible with the universal modeling program HoneyWell UNISIM, all information on physical and chemical properties and vapor-liquid equilibrium was extracted directly from it.

Experimental Study of the Simple Distillation Process
To verify the developed method of kinetic calculation of film mass transfer devices, data from an experimental study of the simple distillation process performed at the Tallinn Polytechnic Institute by Yu. I. Kallas was processed [3,4].
The experiment was carried out in a single-tube film column with a heated outer wall; the scheme of the experimental installation is shown in Figure 3.
The previously accepted values of expenses and flow compositions are replaced by the newly calculated ones.
Based on the assumption that the average discharge compositions of the interacting phases correspond to the saturation conditions, the temperatures of divergent flows are calculated. The newly calculated characteristics of the phase flows are placed in all the calculated cells instead of the previously accepted ones.
After completing the calculation procedure for all sections of the tube (cells), the condition for achieving a stationary state of the system is checked, which is characterized by the fulfillment of the condition: where = 10 is the specified accuracy of the solution. In this case, yD,i and xW,I are defined, respectively, as the composition of the vapor leaving the last (upper) cell and the composition of the liquid leaving the first (lower) cell. If condition (31) is not met, the calculation is performed again. When the condition is met, the calculation is completed.
The method used to find a solution for a multi-element apparatus is a kind of method of simple iterations and has received the name "relaxation method" in the literature [24]. It is proved that its application provides reliable and fairly fast convergence. It is also very economical from the point of view of programming and in machine implementation. In this paper, the calculation program for the distillation apparatus was developed using the ANSYS Fluent software and computing complex. Because this is compatible with the universal modeling program HoneyWell UNISIM, all information on physical and chemical properties and vapor-liquid equilibrium was extracted directly from it.

Experimental Study of the Simple Distillation Process
To verify the developed method of kinetic calculation of film mass transfer devices, data from an experimental study of the simple distillation process performed at the Tallinn Polytechnic Institute by Yu. I. Kallas was processed [3,4].
The experiment was carried out in a single-tube film column with a heated outer wall; the scheme of the experimental installation is shown in Figure 3. In the experiment, the liquid heated to the boiling point flowed down the inner surface of the distillation pipe (d in = 0.016 m, L = 0.85 m, L/d in = 53) and evaporated due to the heat flow supplied to the wall.
The heating conditions of the wall provided a relative constancy of the transverse mass flow along the height of the apparatus. In experiments in the steady-state mode, for a certain period of time, the distillate vapors leaving the top of the column were selected, the relative selection of the distillate was measured, and its composition was analyzed. The study was conducted on four mixtures (benzene-toluene, n. hexane-n. heptane, ethanolwater, and n. hexane-benzene), which significantly differed in their physico-chemical properties. Several series of experiments were carried out on each mixture with different compositions of the initial mixture. Within each series of experiments (25-35 points), the relative selection of distillate (D/F) and the duration of the experiment changed, which made it possible to change the evaporation intensity of the mixture (N c ) and the modes of movement of both the vapor and liquid phases (Re y , Re x ) within a wide range. The total number of conducted experiments exceeded 220.
Yu. I. Kallas also carried out a joint solution of the equations of convective diffusion and continuity using the following assumptions [3,4]:

•
The laminar flow regime of the gas phase (ε D = 0) is considered; • The transverse mass flow (N c ) is constant over the entire height of the pipe; • The equilibrium dependence is linear (y * = mx + b); • A linear profile of concentrations along the pipe height is adopted.
The accepted assumptions allowed the author to obtain an expression that determines the concentration profile in an arbitrary pipe section: Equation (32) allows us to find the value of the concentration gradient on the interface of the phases (r = R), and hence the value of the mass flow of the component through the interface of the phases: The assumptions made by Yu. I. Kallas very much coarsen the decision. In particular, this applies to the neglect of the turbulent diffusion coefficient. Similarly, there is little justification for the assumption of the linearity of the concentration profile along the pipe height, especially for long pipes that are used in real evaporators and distillers.
In addition, as was proved in [21][22][23], the total flow also transforms the diffusion component of the mass flow through the phase interface, both for laminar and turbulent modes, which is not taken into account in the work of Yu. I. Kallas. Nevertheless, with the help of the assumptions made, the author was able to obtain an expression for calculating the boundary concentration (y f ), to relate the kinetic regularities of the concentration of the oncoming vapor and liquid flows in the film column, and using balance ratios to obtain an expression for calculating the composition of the distillate leaving the pipe: where A, B, C, E, I are expressions that take into account the balance relations, the diffusion properties of the vapor and liquid phases of the separated system, and the geometric dimensions (d, L) of the apparatus; m and b are the coefficients of the equilibrium dependence equation. Experimental verification of Equation (34) confirmed its acceptable accuracy [25].
It is of interest to analyze the last assumption used by Yu. I. Kallas [3,4,25], because it allowed the author to obtain an analytical solution to the problem. For a film machine: It is evident here that the concentration gradient in the film apparatus depends both on the current value of the mass flow in the section under consideration and on the longitudinal coordinate. Therefore, the assumption that the concentration gradient is constant over the pipe height is not obvious. It is known that for monotone functions, the second derivative of which preserves the sign on the considered section (0 < z < H), the best linear approximation for the derivative is its value determined in the middle of the section (z = H/2). Then, by solving the system of Equations (35)-(38), we obtain: Equation (39) should be considered as an approximate solution of the film distillation problem in the binary formulation of the problem. The calculation of the flows is carried out according to the mass transfer equations for only two sections of the pipe (z = 0 and z = H). When using the condition on the linearity of the equilibrium dependence, analytical dependences similar to (35) can also be obtained from (39). In the multicomponent formulation, Equation (35), of course, is also not applicable. The use of (39) in this case, although it is possible in principle, is not justified enough. In a more rigorous formulation, the calculation should be carried out directly by numerical integration of the complete system of Equations (23)-(26) according to the algorithm described above. To implement a specific calculation procedure for a simple distillation process, this algorithm can be reduced to the form shown in Figure 4, in which the condition of constancy of N c in the height of the pipe is established. The algorithm can also be easily transformed for other conditions.
To compare the adequacy of the considered models, experimental data were processed [3,4] For comparability of the results, all the physico-chemical properties of the mixtures were calculated using the formulas used in [3,4,30]. The equilibrium state of the systems at the interface was described using the database of the HoneyWell UNISIM program (license 17d3-8598-20a7-b2ec-WEB) [30,31]. A comparison of the calculated values of distillate concentrations with experimental data is presented in Table 1. As a criterion for the adequacy of each of the models, the value of the standard deviation of the calculated values of distillate concentrations from the experimental data for each of the series of experiments was taken:

Results and Discussion
It is clear that, in all cases, when using model No. 2, the standard deviation was found be significantly smaller (by 2-3 times) than when using special models developed only for the simple distillation process and only in the binary formulation of the problem. This fact confirms a certain universality of the proposed mathematical description and the expediency of its use for describing arbitrary non-equimolar mass transfer processes.
The effects of non-equimolarity are determined both by the direction of convective mass flows (N c ) and by the relations between N c and the coefficients of equimolar mass transfer, which would be formed in the considered apparatus in a similar hydrodynamic situation (Re = idem, Sc = idem). This relation for binary systems was previously called by us the non-equimolarity factor of the process, which can be considered as a similarity criterion that does not depend on either the type of the modeled process or its hardware design.
In the process of simple distillation, the effect of redistribution of phase resistances to the mass transfer process is very clearly visible, depending on the degree of nonequimolarity of the process. Indeed, under non-equimolar conditions, the diffusion flow is only a part of the total mass flow. In accordance with the structure of the mass transfer Equations (23) and (24), with a positive direction Nc (distillation process), the relative fraction of the resistance of the liquid phase increases, which is manifested in an increase in the ratio of the concentration difference x − x f / y f − y in the diffusion boundary layers of the contacting phases in comparison with the equimolar process. Our calculations showed that the proportion of the resistance of the liquid phase in the total resistance to mass transfer for the mixtures under consideration was 92-100%. Interestingly, the same effect was noted by Kallas Yu. I. [24]. At the same time, for equimolar conditions for the same mixtures, more than 80% of the mass transfer resistance is controlled by the resistance of the gas phase. It is obvious that neglecting the effects of non-equimolarity in this case can lead to significant errors in the calculations of mass transfer processes. For processes with the opposite direction of N c (condensation, desorption processes), the effects will obviously be reversed.

Conclusions
Thus, based on the generalization and analysis of the relatively extensive experimental material accumulated in previous studies, the adequacy of the proposed structure of the mass transfer equations can be considered proven, in addition to the independence of this structure from the conditions of flow around the interface of the phases, i.e., from the hardware design of the process.
The adequacy of the proposed structure of the mass transfer equations and its independence from the flow conditions of the interface were confirmed by generalization of the experimental material and a numerical experiment on a mathematical model. The program for calculating the mass transfer process in film apparatuses, which is suitable for calculating arbitrary mass transfer processes (absorption, rectification, desorption, evaporation, etc.) performed in a multicomponent formulation, is used in the practice of design work.

Conventional Designations
H-the height of the device, or the pipe section, m; B*-elements of the square matrix of coefficients of multicomponent equimolar mass transfer, m 3 /(m 2 s); kg-mol/(m 2 s); D-binary diffusion coefficients, m 2 /s; d-the diameter of the device, or the pipe section, m; F, G, L-power consumption, steam, liquid phases, respectively, m 3 /s; kg-mol/s; g-acceleration due to gravity, m/s 2 ; idem-the property of an object or operation, when reapplying the operation to the object, give the same result as when the first k-the number of components in a multicomponent mixture;