Analytical studies on deposition and entrainment present in the Venturi nozzle two-phase flow

Introduction

Bearing in mind the concern for dust and chemical pollutions removal from the atmosphere, searching of a so-called clean energy sources is obligatory. The Polish power system mostly relies on the carbon power plants thus Poland is struggling with the greenhouse gases (GHG) and dust particles high emission problem. Even though some works enabling carbon dioxide capture and flue gas purification are being undertaken, the problem is still not being sufficiently solved [ ¹ ]. In the light of this fact, the alternative heat carriers ought to be searched by researchers. The greenhouse gases emission reduction is a critical problem of the world which scientists should work on. Thus works on alternative fuels are of a highest importance [ ² ]. The synthesis gas (syngas), which can be produced through the biomass gasification process, seems to be a very attractive alternative to coal. However, syngas also needs to be purified from different chemicals and dust particles so it could cope the particular restrictions in the last step of its processing which is a key factor in its production technology. According to the literature [ ³ , ⁴ , ⁵ , ⁶ – ⁷ ], Venturi scrubbers are being widely used in the syngas purification process and biohazard materials cleanup, and others as these devices are found to be very attractive from both, efficiency and economical, points of view. Thus, authors have decided to propose an innovative analytical model for Venturi throat two-phase fluid flow pattern description which further can be used for Venturi designing methods refinement. Worth mentioning is the fact that in the current paper there are presented the preliminary studies in this field, and works are going to be continued. In the studied case, the two-phase two-component fluid flow is present in Venturi’s throat of a cylindrical shape. The gas (air) is a continuous phase, and the sorption liquid (water) is a dispersed phase. Thus, there occurs an annular flow with water as liquid wall film and a homogenous mixture in the core of the flow. The content of the solid particles has been omitted in the current considerations. It is obvious that the two-phase flow, its description and analysis is much more difficult, regarding the physics, than a single-phase flow. In the two-phase flow apart from inertia, viscosity and pressure forces present in a single-phase flow, there also occur surface tension forces. Additionally, wetting the internal lateral surface of the channel results in a liquid wall film forming. From one hand, the liquid wall film improves the droplet deposition but its existence also allows for drops entrainment. As a consequence of the mentioned interactions the equations of the movement of both phases are unstable. The two-phase flow pattern plays a significant role in the studied problem of a syngas cleaning process, thus the regimes of this flow cannot be neglected nor disregarded.

Literature review

For several decades, scientists have been conducting research into the effectiveness of Venturi scrubbers to remove pollutants form gas. In the construction of the Venturi scrubbers, there are known three different methods of liquid injection: as a film, as a jet and as a spray. The liquid is introduced as a spray in the ejector type Venturi scrubber which is discussed in works [ ⁸ , ⁹ – ¹⁰ ] whereas sorption liquid injected as a jet is typical for a Pease-Anthony Venturi scrubber which is studied in works [ ¹¹ , ¹² , ¹³ , ¹⁴ , ¹⁵ – ¹⁶ ]. Gonçalves et al. [ ¹¹ ] have presented a mathematical description of the trajectory of the injected liquid. The model takes into account the penetration of the liquid stream and its separation. This description was used by the authors in another article [ ¹² ] to determine the distribution of droplets density. The use of the Venturi scrubber in which a stream of a liquid is sprayed perpendicular to the throat, at the subsonic gas flow, finds many applications in various industrial processes such as combustion, evaporative cooling or in the medical and agricultural devices. However, the model described in the paper [ ¹¹ ], although it may find many applications in other types of scrubbers, according to the statement of the authors of this paper, is not valid for the scrubber used for the gas dedusting. The article [ ¹² ], which is a continuation of the paper [ ¹¹ ], focuses on the theoretical and experimental analysis of the droplets density distribution in the horizontal Venturi scrubber with rectangular cross-section. This issue is very important from the scrubber's efficiency point of view which depends on the liquid atomizing. By ensuring proper distribution of the liquid, greater gas dedusting efficiency can be achieved with a minimum water consumption. The researches were conducted for variable throat length. The liquid was injected through one hole, perpendicular to the throat which enabled the studies on the effect of penetration of the liquid stream on the droplets dispersion. In the article [ ¹² ], a mathematical model from the paper [ ¹¹ ] was used and developed. Unlike to article [ ¹¹ ], the subject of considerations in the paper [ ¹² ] is the Venturi Pease-Anthony scrubber, used to remove particles in the micron range. A very important aspect is the aforementioned liquid atomizing. In the scrubbers, where the droplets concentration is sometimes higher or lower, it has been observed that the purification process is not that effective as in the scrubbers in which the droplets density distribution in the throat area is approximately constant. Venturi scrubbers with liquid injected as a wall film through the confusor inlet, or as a jet, is the subject of investigations in the following works [ ⁵ , ¹⁵ , ¹⁷ , ¹⁸ , ¹⁹ – ²⁰ ]. The dynamics of mechanisms during the gas dedusting process in the Venturi scrubber, regardless of the method of liquid injection, is expressed by the pressure drop. One of the simplest pressure drop prediction model was proposed by Calvert [ ²¹ ]. The model assumes a one-dimensional, incompressible and adiabatic flow in which all the liquid is atomised into droplets of equal size. In addition, the Calvert model [ ²¹ ] assumes droplet acceleration to the gas velocity at the end of the throat while the effect of the gas acceleration in the throat and the wall friction are neglected. The preliminary Calvert model [ ²¹ ] was improved by Yung [ ²² ] (assuming that the droplet velocity at the end of the throat will not be equal to the gas velocity), and Leith [ ²³ ] (which takes into account the relative slowdown of the droplets in the diffuser). Other simple models of the pressure drop prediction are the Hesketh [ ²⁴ ] and the Boll [ ²⁵ ] models from which the first one assumes that droplets accelerate and finally reach the gas speed in the throat. The latter model [ ²⁵ ] includes the gas and droplets acceleration and the friction and assumes the fully atomised liquid. The model of Azzopardi and Govan [ ²⁶ ] was improved in works [ ²⁷ , ²⁸ – ²⁹ ] in which the growth and detachment of the gaseous boundary layer was determined. Obviously, regardless to the liquid injection method, a liquid wall film forms on the channel wall which more or less affects the overall pressure drop and the gas dedusting efficiency which was considered in papers [ ²⁴ ]– ²⁷ . The two-phase flow annular pattern was also considered in works of Viswanathan et al. [ ⁸ , ³⁰ ] in which the thickness of the liquid wall- film against the overall Venturi orifice length was determined. Additionally, the velocity of a liquid wall film and the pressure drop were estimated in a wide range of the Venturi operation conditions. The model of Azzopardi and Govan is described in papers [ ⁵ , ²⁶ ] and [ ⁹ ]. Azzopardi and Govan took into account the entrainment and deposition of the droplets and the growth of the gaseous boundary layer in the diffuser area to correctly predict the pressure drop. This model correctly explained that the excess pressure recovered from the previous models was a result of the increase of the gaseous boundary layer in the diffuser. Azzopardi and Govan [ ⁵ ] presented a one-dimensional model describing the flow of both liquid, dust and the gas. All previous Authors of models, like a Yung, Boll, Calvert, Miheisi, Chongvisal, Holland and Goel, Placek, Peters and Bayvel, as well as others, have considered the flow of the liquid only in a droplet form, and in the meantime, a liquid film is always formed on the walls of the scrubber, and it has a significant impact on the purification process [ ⁵ ], as it is discussed in this article. In addition, under certain conditions, when the dust exhibits particular adhesiveness or when there is a problem with a throat erosion, caused by a high velocity of the gas flow, the liquid can be intentionally injected as a film onto the channel walls and then atomized by the gas. Neglecting the existence of the liquid film on the walls of the scrubber, as it is usually done, some important factors are being ignored. First of all, the liquid that flows as a liquid film is not in the form of drops, so it does not take a part in the gas purification process. Secondly, the liquid film will behave like a rough surface in direct contact with the flowing gas which will increase the pressure drop due to friction. Thirdly, due to the presence of a liquid film along the entire length of the Venturi orifice, there will be a constant accumulation of droplets on the wall and their entrainment. Therefore, there will be new droplets of a velocity lower than the gas so they will be used effectively in the gas cleaning process. The model presented by Azzopardi and Govan [ ⁵ ] takes into account the flow of a part of the liquid in the form of a film on the wall, as well as the transition between the droplets and the liquid film. The paper describes a model of the annular flow which is based on the equation describing the variation of the liquid wall film mass flow rate, resulting from the droplets deposition and entrainment. An important discovery by Ekman, Johnstone and Muir and others were that the efficiency of the Venturi increases with the increase of the gas to liquid ratio at a constant pressure drop [ ⁵ ]. In the article [ ³¹ ], attention was drawn to the fact that the spatial distribution of the droplets and their size have a very large impact on the pollutants removal efficiency. The spatial distribution of the droplets began to be considered to improve the prediction model of the purification efficiency.

Analytical model of the two-phase flow in Venturi’s throat

Authors of the paper have carried out a theoretical study on the two-phase air–water flow pattern in Venturi’s throat. In the analytical one-dimensional model presented in the paper, it was assumed that the liquid flows as droplets in the core of the flow and also as a film on the wall. In the case of the two-phase flow pattern presented in the paper, neglecting the liquid wall film flow results in lowering the correctness of the model. As it was mentioned by Azzopardi et al. [ ⁵ ], liquid film existing on the channel wall doesn’t take a part in the chemical pollutants removal and dust collection until there occurs the droplet entrainment from the wall film to the core of the flow. Additionally, the droplets which entrain from the liquid film flow to the two-phase core-flow have smaller velocity thus the pollutant removal efficiency increases [ ⁵ ]. Moreover, the liquid wall film affects the friction component of the momentum equation of the two-phase flow thus neglecting the liquid wall film results in incorrectness [ ⁵ ]. Current investigations based on the authors’ analytical model have led to the estimation of the pressure drop along Venturi’s throat and liquid fraction in the flow. The introduced authors’ model allows for the pressure drop calculation with the use of the core-flow parameters which are listed further in the paper. Additionally, the model allows to define the mass flux of the liquid film at the channel walls and the mass flow rate of the liquid in the core of the flow. Present model takes into consideration two important phenomena existing in the two-phase flow. Namely, the deposition of the liquid droplets within the liquid wall film and the entrainment of the drops from the wall film to the core of the two-phase fluid flow. The two-phase flow present in Venturi’s throat is a homogeneous air–liquid droplets mixture in the core of the flow and a liquid film at the channel wall. The total liquid (water) and gas (air) mass flow rates are given, and the superficial liquid velocities in the core of the flow and in the liquid wall film are the resultant ones. Obviously, the two-phase flow is being also characterized by the momentum exchange between both phases due to which the motion of the one phase may strongly affect the other’s flow. In the current preliminary studies, momentum exchange through a phase separation boundary was estimated, however, it was not included into calculus as the liquid content is low. Consequently, the calculation results have shown that momentum transfer can be neglected in the current analytical research studies. However, for highest values of the liquid to air ratio, it must be included as the shear stresses, resulting from the momentum exchange, also influence the value of the pressure drop.

As illustrated in the Fig.  ¹ , the following balance equation of the liquid film mass flow rate can be written:

Fig. 1

and hence

Equation ( ⁵ ) was first proposed by Whalley [ ³¹ ], however, he has included the heat flux term which is skipped in Eq. ( ⁵ ) as heating of the channel walls is not being considered in the current paper.

Introducing non-dimensional parameters:

The Eqs. ( ⁵ ) and ( ⁴ ) can be written in the following non-dimensional form

In accordance with [ ³² ], it was assumed that

The deposition rate can be determined by the following formula [ ³² , ³⁵ , ³⁶ ]:

The model described by Whalley was incorporated into the Harwell annular flow modelling code HANA in which the simplified correlations for deposition Eq. ( ¹⁰ ) and entrainment were first used [ ³⁷ ]. The non-dimensional form of Eq. ( ¹⁰ ) is given by:

For droplets of a diameter smaller than one micron, the deposition mass transfer coefficient is defined as follows [ ³⁸ ]:

For the droplets of a diameter greater than one micron, the deposition mass transfer coefficient is estimated from the following equation [ ³⁸ ]:

The Darcy-Weisbach drag coefficient is presented in a form of:

Equation ( ¹⁸ ) is valid for 3·103>Re>105 {\text{Re}} > 10^{5}$$\end{document}]]> . For the Reynolds Number values in the range of 105>Re>108 {\text{Re}} > 10^{8}$$\end{document}]]> , the following Nikuradse correlation should be used:

In the calculus having been carried out for the purpose of the studied issue, Reynolds Number was estimated for both, the liquid and the vapour, for varying liquid to gas ratios. The Reynolds Number value for the air, for the cases studied in the paper, was in the following range Reair∈2.7·105,1.12·106 . For these conditions, drag coefficient computed using ( ¹⁸ ) and ( ²⁰ ) varies only by few percent which is negligible as the current studies are of a quality character instead of quantity.

After some mathematical transformations, formula describing the deposition mass transfer coefficient for droplets larger than 0.1 micron is obtained

Bearing in mind Eqs. ( ¹¹ ), ( ¹³ ) and ( ²³ ), the following equation describing the non-dimensional deposition flux may be obtained

The quality of the core-flow, xc , is defined using the equation presented below

The liquid concentration in the core of the flow is exposed through the following formula:

Therefore, substituting ( ²⁶ ) and ( ²⁷ ) to ( ²⁴ ), with respect of ( ²⁵ ), leads to the final expression of D+ which is presented below:

Through the following equation, with respect to Eqs. ( ⁹ ) and ( ²⁸ ), the non-dimensional entrainment rate is being introduced

Regarding Eq. ( ⁸ ), with introducing Eqs. ( ²⁸ ) and ( ²⁹ ), the final form of the first order differential equation being considered is exposed and then solved analytically:

Differential Eq. ( ³⁰ ) was solved analytically, and its solution for z ⁺ is following

As it can be seen, the solution is implicit function of the dependent variable gf . The proposed model enables prediction of the throat length along which there occur both the deposition and entrainment of the droplets. Moreover, the wall film mass flux variation in the channel is being estimated. The distance at which there occurs the droplet entrainment is an individual parameter which depends on the liquid and gas inlet mass flow rates and the channel diameter. It is commonly known that the droplet entrainment may occur only if the following condition determining the critical wall film mass flux is being fulfilled [ ⁴⁰ ]:

The scheme of estimating the momentum exchange at the phase separation boundary, through which there occur liquid deposition and entrainment, is introduced with the use of the following equations and was calculated regarding deposition and entrainment fluxes [ ⁴¹ ]. This analysis is being presented to prove that in this particular case the momentum transfer at the phase separation boundary can be neglected:

τw [kg/(m s ² )] is the shear stresses at the channel wall in the case of a mass flux transfer through the phase separation boundary, τ [kg/(m s ² )] is the shear stresses at the equivalent case with mass flux transfer through the phase separation boundary disabled, VDep [m/s] is the velocity of the droplets which deposit at the liquid wall film, VEntr [m/s] is the velocity of the droplets which entrain from the liquid film. The value of the velocity, VDep/Entr , was estimated based on the results obtained through the proposed analytical model as it is the velocity of the droplets deposition/entrainment. Therefore, VDep and VEntr , were estimated for the purpose of the current calculus and were computed with the use of Eqs. ( ²⁸ ) and ( ²⁹ ), defining deposition and entrainment flux, respectively. In the following set of equations the parameters τDep/Entr+ , τw,Dep/Entr and VDep/Entr for both, deposition and entrainment fluxes are presented:

Through Eqs. ( ³⁷ , ³⁸ , ³⁹ , ⁴⁰ , ⁴¹ , ⁴² , ⁴³ ), a scheme of Δτ estimation is introduced, where Δτ=1.08·10-4 Pa is the shear stress increase, resulting from the momentum transfer through the phase separation boundary. As it can be seen, it is a negligible value. Equation ( ³⁴ ) is valid for low b values. The coefficient ξ was computed in a way presented earlier in the text with the use of the Reynolds Number, calculated regarding the boundary layer thickness, δ , presented in the form of the formula:

δ [m] is the boundary layer thickness [ ⁴² ], νν [m ² /s] is the air kinematic viscosity, uv,∞ [m/s] is the undisturbed continuous phase velocity.

The model assumes two-phase air–water flow in Venturi’s throat. The air is a continuous phase, and the water is a dispersed one, and it is introduced to the flow as droplets of a uniform diameter of 10 ^–5  m. The atomised liquid is being injected at Venturi’s throat inlet, and the two-phase flow is thus being formed. The drops which flow in the core of the flow deposit from the very beginning of the channel, and so, the liquid film is formed at the wall. Obviously, the thickness of the film varies with the number of the droplets being deposited on or entrained from the liquid wall film. The wall film thickness increases until it reaches a critical value at which the rolling waves are being formed, and the droplets entrainment phenomenon starts to occur. Regarding the analytical model validation case in which there were applied water and air mass flow rates, m˙air  = 0.483 kg/s and m˙w  = 0.013 kg/s, coming from the literature, the droplet entrainment starts to occur at the distance of 0.084 m of Venturi’s throat. Figure  ² shows how do the non-dimensional parameters Gf+ and Gc+ defined by ( ⁷ ) vary with the non-dimensional throat length. In Fig.  ² , there can be seen a small curvature of the graphs. This curvature increases with the liquid to gas ratio, with the throat length and with the throat diameter decrease which is discussed further in the paper.

Fig. 2

Fig. 3

Fig. 4

Analytical model validation

The model presented here concerns the liquid superficial velocity variation within the wall film and within the core of the two-phase flow along Venturi’s throat. According to the authors’ knowledge, there are no available data in the literature which could be used to verify the accuracy of the model predictions. Nevertheless, the works of Viswanathan [ ⁸ ] and Azzopardi et al. [ ⁵ ] confirm the correctness of the deposition-entrainment two-phase fluid flow modelling approach. Thus, the analytical model validation has been carried out indirectly, through the pressure gradient analysis and comparison to the experimental data published by Silva et al. [ ¹⁷ ] for a large-scale Venturi. Accordingly, in the analytical model calculations it was assumed that the two-phase air–water mixture flows at the rate of m˙air  = 0.483 kg/s and m˙w  = 0.013 kg/s on the inlet, respectively. Also, the diameter and the length of Venturi’s throat are equal d  = 0.1225 m and L  = 0.3 m, respectively (Fig. ³ ). Liquid to gas ratio, L/G , in this case is equal 0.027. The experimental pressure drop in Venturi’s throat was equal to 114 Pa. The pressure drop was calculated using the following equation [ ⁴⁰ ] (Fig. ⁴ ):

In the first, frictional term,

and X is the Lockhart-Martinelli parameter, defined as follows:

and α1 and α0 were computed using Eq. ( ⁵² ). The parameters BB , n , n2 and n3 present in Eq. ( ⁵² ) can be found in [ ⁴⁰ ], in Table 10.2. In the current studies, the following values of the parameters BB=1 , n=1 , n2=1 , n3=0 , corresponding to the homogenous two-phase flow regime, were used. The analytical two-phase fluid flow pressure drop was calculated using Eq. ( ⁴⁶ ) which is common in the literature, e.g. [ ⁴⁰ ]. In the current analytical studies Venturi’s throat is being considered. However, experiments, on the basis of which the model is being validated, have been carried out for the whole Venturi orifice. Thus, in the pressure drop calculus there must be included the confusor local pressure drop, resulting from the variation of the cross-section, to ensure the accordance of the theoretical model with the experiment. In the pressure drop calculus all needed parameters such as the friction factor for vapour flowing alone in the tube, the quality, the two-phase multiplier and the void fraction come from the analytical model results. The pressure drop was calculated in two ways. In the first one, the two-phase fluid flow parameters were used for the total flow in Venturi’s throat (Eq. ( ⁴⁶ )). In the second one, only the parameters of the core of the flow were used. The total pressure drop and the pressure drop in the flow core should be equal. Calculating the pressure drop based on the core-flow parameters requires the knowledge about the flow pattern in the core which results directly from the analytical model results presented in the paper. The core pressure drop was calculated with the use of the following equation (Fig. ⁴ ):

As a part of the research on the Venturi gas dynamics, the CFD analysis has been carried out for the single-phase air flow, m˙air  = 0.483 kg/s, [ ³ ] and compared to the literature data [ ¹⁷ ].

Analytical model calculation results

The calculations have been carried out also for m˙air  = 2 kg/s and m˙w  = 0.7 kg/s and other values, listed in the Table ² . The diameter of Venturi’s throat is the same as used in the previous chapter, that is d  = 0.1225 m. The calculations results presented in this section have been carried out for two Venturi throat lengths, which are 0.3 m and 1.35 m. It can be seen in Fig.  ⁵ how that the pressure drop varies with liquid (water) to gas (air) ratio in the throat of 0.3 m length. Additionally, it is shown there how the core void fraction varies with the L/G ratio. The character of the two graphs in Fig.  ⁵ is rational as the void fraction in the core of the flow decreases with L/G ratio increase, and the overall pressure drop increases with L/G increase.

Fig. 5

Droplet entrainment activated

Figures from Figs.  ⁵ , ⁶ , ⁷ present the results for the case of the entrainment enabled. Obviously, as it was mentioned earlier in the text, the droplet deposition is present from the very beginning of the flow channel while the droplet entrainment starts to occur only if the condition defined by Eq. ( ³³ ) is fulfilled.

Fig. 6

Fig. 7

Figure  ⁶ presents the variation of the water mass flux in the liquid wall film, Gf+ , and in the core of the flow, Gc+ , through 1.35 m Venturi throat. These characteristics have been obtained for L/G  = 0.35. As it can be clearly seen the curvature of the lines is greater than in the case of the L/G  = 0.027 shown in Fig.  ² .

The homogenous model of a two-phase flow, applied for the purpose of the current analytical research studies yields reasonably accurate results for limited range of the flow conditions (bubbly flow or mist flow). As it is shown in Fig.  ¹ , it is assumed that the mist flow exists in the core of the investigated two-phase flow domain (Table ³ ). The void fraction, α, was calculated using Eq. ( ⁵² ). In Figs.  ⁶ , ⁷ there are presented the results for the two-phase air–water flow, of the L/G ratio equal 0.35, with droplet entrainment and deposition activated as the critical mass flux was reached from the very beginning of Venturi’s throat.

Table 3

d [m] →  L/G↓	0.05	0.06	0.12	0.18
0.027	3.181•10–5	3.219•10–5	3.299•10–5	3.302•10–5
0.35	4.094•10–4	4.097•10–4	4.22•10–4	4.28•10–4
0.45	5.3075•10–4	5.3081•10–4	5.36•10–4	5.39•10–4

Droplet entrainment deactivated

Bearing on mind the validity of revealing the real influence of the droplet entrainment on both superficial water velocities, Gf+ and Gc+ , the case of the droplet entrainment neglecting is introduced in this paragraph (Figs.  ⁸ , ⁹ ). Therefore, figure exposed below (Fig.  ⁸ ) shows the similar characteristics as in the previous figure (Fig.  ⁶ ), however, the droplets entrainment is deactivated.

Fig. 8

Fig. 9

Analysis of Venturi’s throat geometrical parameters influence on the fluid flow pattern

The extreme channel length of 12.25 m was chosen intentionally to clearly visualise the overall Gf+ and Gc+ variation along such along channel. With the help of Fig.  ¹⁰ a, remembering that z+=zd , one can determine the distance from the throat inlet where the steady state with respect to Gf+ and Gc+ is reached. Additionally, Fig.  ¹¹ can be used to find the critical length of Venturi’s throat at which a highest liquid fraction exists. Based on the current analytical model results it can be stated that, for the specific fluid flow parameters, assumed in the paper, the large-scale Venturi throat length could be less than 1 m. Large-scale Venturi nozzle can even reach the length of 0.65 m so, the results presented in the current paper are in a good accordance with the available literature data [ ⁴⁴ ]. From Figs.  ¹⁰ and ¹¹ , it can be seen that the differential Eq. ( ³⁰ ) is of a nonlinear character.

Fig. 10

Fig. 11

The volumetric liquid fraction variation in Venturi’s throat is presented in Fig.  ¹² . As it can be seen, lowering the throat’s diameter strongly affects the liquid content, and it decreases with the diameter reduction. Thus, for the purpose of the syngas purification process large Venturi orifices should be used.

Fig. 12

Based on Eq. ( ⁵⁴ ), the mean core liquid fraction was estimated from the relationship

βc=1-φc

It was found that βc varies significantly with L/G ratio and slightly with the throat diameter d . This parameter plays a significant role in the syngas purification process as it gives the information about a droplet content in the in the core of the flow. Each droplet may collect particular amount of solid particles. Estimation of a droplet content in the core of the flow is therefore necessary in the estimation of dust particles collection efficiency. In general, the greater βc , the higher efficiency is being reached.

Conclusions

Current investigations based on a novel analytical model, presented in the paper, have led to the estimation of the pressure drop and liquid fraction in the flow, along the Venturi throat. The introduced model allows for the pressure drop calculation with the use of the core-flow parameters. Additionally, the model allows for determining the mass flux of the liquid film at the channel walls and the mass flow rate of the liquid in the core of the flow. Present model takes into consideration two important phenomena existing in the two-phase flow. Namely, the deposition of the liquid droplets within the liquid wall film and the entrainment of the drops from the wall film to the core of the two-phase fluid flow. The innovative concept presented in the paper leads to set of equations which is solved analytically in the paper. The proposed model enables prediction of the throat length at which there occur both the deposition and entrainment of the droplets. Moreover, the wall film mass flux variation in the channel is being estimated. The distance at which there occur the droplet entrainment is an individual parameter which depends on the liquid and gas inlet mass flow rates and the channel diameter. The analytical studies on the two-phase fluid flow with droplets deposition and entrainment, presented in the paper, have led to the conclusion that the liquid mass flux variation in the film and in the core of the flow is nonlinear. The results also show that neglecting the droplets entrainment leads to the linearization of the mentioned characteristics (Figs. ⁸ , ⁹ .) In the paper, the core-flow parameters were determined and can be used in the pressure drop estimation with a satisfactory accuracy. Investigations on the influence of the throat length and diameter variation on the studied two-phase flow have led to the statement that large diameter (Fig.  ¹² ) and throat length less than one metre (Fig.  ¹¹ ) may lead to higher cleaning efficiency. Additionally, the liquid fraction in the core-flow was estimated which is essential when it comes to the cleaning efficiency estimation. The analytical model gives reasonable results, and it is planned to be further expanded by adding confusor and diffuser to the Venturi cylindrical throat which was analysed in the current paper. Moreover, based on the presented analytical considerations, dust collection efficiency is planned to be estimated in the future work.