1. Introduction
The bearing chamber is an important part in the aero-engine lubrication system. The lubricating oil is shed into the bearing chamber in the form of oil droplet via under-race lubrication method [
1]. After the high-speed moving oil droplets collide with the wall, some of them are deposited to form the oil film, and others are splashed to form a large number of secondary oil droplets. Then, the complex air–oil two phase flow and heat transfer state of air, oil droplet, and oil film coexist in the bearing chamber. Because the oil droplet size is very small, many secondary oil droplets are suspended in the bearing chamber under the action of the air phase flow field. Furthermore, the interaction between oil droplets induces many dynamic events, such as collision, coalescence, deposition, evaporation, breakup, etc. The coalescence and breakup of oil droplets will affect the size distribution of secondary oil droplets and further restrict the research of flow and heat transfer in the bearing chamber. Moreover, the oil droplet size distribution directly affects the mixture state of oil and gas. Previously the Rosin–Rammler (R–R) distribution is widely used for the calculation of oil droplet size distribution in bearing chamber. However, experimental research indicates that R–R distribution is not consistent with the actual condition. Therefore, the study of droplet diameter distribution of the secondary oil droplets, which considers the coalescence and breakup of oil droplets, is of great significance in the aero-engine lubrication design.
At present, the research on coalescence and breakup of oil droplets mostly adopts the population balance model, which is earliest used in chemical process by Hulbert et al. [
2]. Since then it has become a general method to analyze the particle diameter distribution of dispersed phase in multiphase system. The population balance model describes the time and space distribution of discrete phase by tracing the change of number density function of discrete phase with convection, diffusion and the source term of coalescence and breakup in the flow field. This method has been widely used because it remains part of the size information of discrete phase, and uses the averaging method to describe the motion of continuous phase, which reduces computation greatly. The core of population balance model lies in the profound understanding of coalescence and breakup mechanism of oil droplets. The coalescence and breakup of oil droplets is affected by the flow field around the droplets and the physical parameters (such as surface tension) of the oil droplet, which involves a very complicated multi-level and multiscale physical mechanism. The mechanism description of oil droplet breakup can be divided into three categories: (1) the droplet breakup caused by turbulence, (2) the droplet breakup caused by viscous shear force, and (3) the droplet breakup caused by the instability of the liquid–liquid interface. The coalescence mechanism of oil droplet is more complex. Moreover, the probability model is used more often, as the coalescence rate of oil droplet is equal to the product of collision frequency and coalescence efficiency between oil droplet. The collision mechanism of oil droplet can be divided into four categories: (1) collision caused by random turbulence fluctuation, (2) collision caused by buoyancy, (3) collision caused by shear force, and (4) collision induced by wave vortex. The coalescence efficiency model can be divided into three categories: energy model, velocity model, and liquid membrane drainage model.
During the past decades, many scholars have done many researches on the coalescence and breakup of bubbles/droplet and the size distribution of bubbles/droplet. Liao and Lucas et al. [
3,
4] reviewed the coalescence and breakup model of liquid, and introduced the research process of oil droplet coalescence and breakup mechanism. Ramkrishna et al. [
5] investigated the effect of coalescence and breakup of oil droplet on the size distribution by using the population balance model, and further studied the flow behavior of liquid–liquid system from mechanism. Coulaloglou and Tavlarides et al. [
6] assumed that the oscillating droplet will break up when the turbulent kinetic energy transmitted by the turbulent vortex to the droplet exceed the surface energy of the droplet, but the two droplets will coalesce when the contact time exceed the film discharge time of the two droplets. Tsouris and Tavlarides et al. [
7] studied the droplet coalescence behavior in the turbulent using high-speed imaging technology and population balance model, proposed a novel collision frequency model, and compared the predicted droplet diameter with the experimental data. The results show that the model is applicable under different discrete phase fraction. Luo and Svendsen et al. [
8] established a bubble breakup model, in which the bubble breakup rate is equal the product of the collision frequency with the turbulent vortex and the bubble breakup efficiency.
Regarding the solution of the population balance equation, Li et al. [
9] used the quadrature method of moments (QMOM) to study the droplet breakup and coalescence behavior of the liquid–liquid dispersion system in the stirred reactor. The results show that the droplet coalescence behavior occurs with some difficultly in the low discrete phase fraction, due to the small probability of droplet collision. Li et al. [
10] reviewed the CFD-PBM simulation in extraction column, and the principles and advantages and disadvantages of different solving methods of population balance model were compared. Xing et al. [
11] used the experimental method and CFD-PBM coupled model to analyze the influence of viscosity in the bubble bed on the hydrodynamic behavior in the bubble area, such as the total gas holdup, the big and small bubble holdup, and the bubble size distribution. The results show that the CFD-PBM coupled model has a good prediction ability for the hydrodynamic behavior of the bubble bed in a wide range of gas velocities and viscosities. Wang et al. [
12,
13] used the population balance model to describe the bubble size distribution in the gas–liquid system, considered the different mechanism of bubble coalescence and breakup, established a relatively mature bubble coalescence and breakup model, and calculated the bubble size distribution in the gas–liquid system by numerical solution of PBM. Zhao et al. [
14] used a multi-Monte Carlo algorithm to simulate the collision and coalescence process of particles in nanoparticle flow, and the simulation results were in good agreement with the direct numerical simulation results. Chen et al. [
15] analyzed the parameter relationship between the structural condition of bearing chamber and the oil droplet size distribution, clarified the oil droplet mass distribution based on the continuous oil drop size distribution, and introduced the concept of continuous oil droplet diameter into the analysis of oil droplet movement. Glahn et al. [
16] measure the size distribution of oil droplet using the phase Doppler particle analyzer (PDPA) technology under various operating condition in the rotating disk chamber. In conclusion, the population balance theory has been widely applied in the analysis of coalescence and breakup between bubbles and droplets. However, no relevant research work has been done to consider the influence of coalescence and breakup on the oil droplet size distribution in the bearing chamber.
Therefore, based on the population balance model, the coalescence and breakup behaviors of oil droplets are described through the coalescence and breakup model of the oil droplet. Then, the air–oil two-phase flow model and the population balance model are solved by Fluent software. The numerical simulation of the coalescence and breakup processes of the oil droplet in a bearing chamber is performed. The air phase velocity distribution and the volume fraction of the oil droplets at different times are discussed. Moreover, the variation rules of oil droplet size distribution are obtained in different initial oil droplet diameter, air inflow, and oil inflow. The coalescence and breakup of oil droplets can change the size distribution of secondary oil droplets. Moreover, the oil droplet size distribution directly affects the mixture state of oil and gas. The research work in this paper can improve the accuracy of air–oil two-phase flow and heat transfer analysis, and provide reference for further research on oil mist concentration test in the aero-engine bearing chamber.
2. Numerical Model
2.1. Geometric Model and Physical Parameters
The geometric structure of the bearing chamber is shown in
Figure 1. The structural parameters and values of the bearing chamber involved in the analysis mainly include the rotor radius of the bearing,
rs = 60 mm; the height of the bearing chamber,
hb = 30 mm; the width of the bearing chamber,
wb = 30 mm; and the diameter of oil scavenge
ds (scavenge is shown in
Figure 2) and the air vent
dv are both 16 mm. Seal air flows into the bearing chamber from air inlet and discharges through vent. The gap of oil inlet is 2 mm. Moreover, the lubricating oil is shed into the bearing chamber from the clearance between the bearing outer ring and the cage (oil inlet in
Figure 1) and discharges through scavenge. The equivalent size of oil inlet is 4 mm (oil inlet in
Figure 2).
In the analysis, the oil brand [
17] is aviation lubricating oil 4109. The oil density
ρl is 926 kg/m
3, the oil dynamic viscosity
μl is 0.007 Pa·s, and the oil surface tension coefficient
σl is 0.035 N/m. The gas density
ρg is 1.225 kg/m
3 and the air dynamic viscosity
μg is 1.789 × 10
−5 Pa·s.
2.2. Mesh Model and Grid Independence Verification
According to the structural parameters of the bearing chamber, Gambit software is used to build up the geometric model of the flow region, hexahedral grids are used to discrete the flow region, and the grid of air inlet and oil inlet is refined. The grid model is shown in
Figure 2.
In order to eliminate the influence of mesh on the calculation, three groups mesh with different size are generated on the gambit platform, and the number of meshes is 487,585, 609,719, and 1,530,258, respectively. Under the condition of air inflow 12 g/s, the air phase velocities at four position are extracted on the section z = 0, and the four coordinate points are as follows; (−0.11,0,0), (−0.09,0,0), (0.09,0,0), (0.11,0,0). The calculation results are shown in
Figure 3. From the figure, the velocities at four position obtained by three groups mesh are approximate. In order to save computing cost, the grid model of 609,719 is chosen in this paper.
2.3. Mathematical Model
In the air–oil two-phase flow, the dynamic events of oil droplet include collision, coalescence, breakup, deposition, etc. Under the combined effect of these dynamic events, the oil droplet diameter distribution changes with time. Under a series of physical actions, not only do the conservation of mass, momentum, and energy between oil droplet and air need to be considered, but also the population balance equation needs to describe the oil droplet diameter distribution. Especially, the characteristics of oil droplet collision, coalescence, and breakup are significant in the bearing chamber. It is necessary to focus on the change of oil diameter before and after coalescence and breakup. Therefore, the parameters, such as the volume fraction and velocity of each phase, are obtained by solving the Eulerian–Eulerian two-fluid model, and the dynamic events of oil droplet are described by using the population balance model. The oil droplet diameter distribution with the change of time and space is obtained by combining the two models. The Eulerian–Eulerian two-fluid model and the population balance model in this paper are as follows.
2.3.1. Two-Fluid Model
In two-fluid model [
18], the air and oil phase are both treated as interconnected continuous media. As the volume occupied by one phase cannot be occupied by the other, the phase volume fraction is introduced. The volume fraction is a continuous function of space and time, and the sum of the volume fractions of each phase is 1. The volume fraction equation is as follows,
where
V is the total volume of air oil two-phase,
Vq is the volume of the
q phase,
αq is the volume fraction of the
q phase, and
n has two phases: gas is labeled as phase 1 and oil droplet is labeled as phase 2. It should be noted that a two-fluid model can be considered a valid approximation only up to Stokes number ~0.1–0.2. A similar set of continuity equations and momentum equations can be derived for each phase.
The continuity equation is
where
is the density of the
q phase;
is the velocity vector of the
q phase.
The momentum equation is
where
is the viscous-stress tensors of the
q phase,
p is the pressure share by gas phase and the liquid phase,
is the Reynolds stress tensors,
is the gravitational acceleration vector, and
is the interfacial force term between the gas phase and liquid phase. When q is oil droplet,
is positive; when q is air,
is negative. The term
is usually composed by three forces: drag, virtual mass, and lift. In present work, only the drag force is taken into account, whereas the others are neglected. The drag force is calculated as follows,
where
is the slip velocity and
,
are the average vector of the air and oil, respectively.
CD is the drag coefficient;
d32 is the mean Sauter diameter, calculated as the ratio between the moments of order three and two of the droplet size distribution.
In the two-fluid model, the turbulent Reynolds stress should be closed. The Reynolds stress arises from the operation of Reynolds averaging on the momentum equation. In this paper, the
k-ε turbulence model [
19,
20] is selected to close the turbulence term in two fluid model. The turbulence model
k-ε can be represented by the following equations,
where
k is turbulent kinetic energy;
ε is turbulent energy dissipation term;
is turbulent viscosity, which can be expressed as
,
,
,
,
, and
are model constants; and
is the generating terms of turbulent kinetic energy.
2.3.2. Population Balance Model
The application of the population balance model [
21] is used to calculate the oil droplet diameter distribution in the bearing chamber. The transport equation of the oil droplet number density function can be expressed as
In the right-hand side of the equation in turn are the source term (I) generated by coalescence, the sink term (II) generated by coalescence, the source term (III) generated by breakup, and the sink term (IV) generated by breakup. Where is the number distribution function (1/m6), which means the number of oil droplets in the range of V and in unit volume is ; is the coalescence rate function (m3/s) of oil droplets in volume V and V′ due to collision. is the breakup rate function (1/s) of oil droplets in volume V′. is the probability distribution density function (1/m3) of oil droplets in volume , which breaks up from oil droplets in volume V′.
The population balance equation (PBE) is a hyperbolic integral-partial differential equation; only a few simple cases have an analytical solution, so it is necessary to use a numerical method to solve the PBE. At present, the usually used methods include the discrete method, moment method, and Monte Carlo method. Based on the calculation cost and accuracy, the discrete method is used to solve the population balance equation in this paper.
The basic idea of the discrete method is dividing the continuous diameter distribution of oil droplets into
N discrete subintervals, and the size of all the droplets in the subinterval is equal to the node value
. Meanwhile, the distribution function of the oil droplet number can be approximately as follows,
where
Ni is the number of droplets in unit volume.
In most cases, the oil droplet size calculated by the coalescence and breakup model is not consistent with the node value of the subinterval, so it is necessary to allocate such droplets to the node value of the interval in a certain proportion, and the same time guarantee the conservation of the total oil droplet mass and the total number of oil droplet. Based on the above principles, Ramkishna proposes that the ratio of the newly generated droplets in the interval
to the node
and
is
and
, and satisfies that
By substituting Equations (9), (11), and (12) into Equation (8) and through a series of algebraic transformations, the discrete droplet transport equation is given by Equation (13).
where
After the droplet is broken up, the distribution proportion of each sub interval is expressed as
Luo coalescence model [
8] is adopted in this paper. The coalescence rate in model is defined as the frequency of new droplets generated after the collision of two droplets, and the volume of the two oil droplets is
and
. The coalescence rate is given by
where
is the collision frequency;
is the probability of coalescence. The collision frequency is obtained by
where
is the characteristic velocity of collision between two oil droplets with diameter
and
, and the number density is
and
, which is defined as
where
is constant,
is the density of each phase.
Laakkonen breakup model [
22] is adopted in this paper, which considers that the breakup of oil droplet is affected by both surface tension and viscous force; the degree of influence depends on the magnitude of two forces; and the collision model can be expressed as
where
.
The probability distribution function of sub droplets is given by
The final breakup frequency is
2.4. Initial and Boundary Conditions
In the transient analysis, the boundary conditions include the inlet and outlet conditions of gas phase and lubricating oil, as well as the relevant wall conditions. The boundary conditions are treated as follows.
- (1).
Inlet boundary conditions: mass flow is adopted at the air inlet and oil inlet. At the air inlet, the initial air mass flow is 12 g/s. At the oil inlet, the initial oil mass flow is 20 g/s, and initial oil droplet diameter is 28 μm.
- (2).
Outlet boundary conditions: outlet boundary conditions adopt pressure-outlet and set the liquid reflux rate as 0 to ensure the stability of calculation.
- (3).
Wall boundary conditions: all walls adopt a no-slip velocity boundary condition. The shaft wall adopts a moving wall and the rotor speed is 16,000 r/min.
The calculation in this paper does not consider the effect of temperature. The oil droplets collision, coalescence, breakup, evaporation and deposition, nucleation and other dynamic events, as well as the geometry of the bearing chamber will affect the oil droplet size distribution in the bearing chamber. However, taking all the factors into account at the same time can make the analysis more complicated. Therefore, the calculation in this paper assumes that the temperature is constant and does not take into account the evaporation of oil droplets and the change of fluid properties.
2.5. Calculation Method
The coalescence and breakup model of the oil droplet in the bearing chamber is solved by using ANSYS Fluent 17.0. The 3D unsteady model is adopted, the first-order upwind scheme is used, and the phase Coupled SIMPLE algorithm modified by pressure and velocity is employed for solving the transient equation. The coalescence and breakup model are as follows. The coalescence model uses the Luo model; the breakup model uses Laakkonen breakup model, in which both surface tension and viscous force are considered; and the oil droplet breakup is combine constrained by two forces, which is closer to the reality. The computation is performed on a Microsoft Windows 10 enterprise 64-bit operation system, and the simulations are run using 40 cores with 192 GB RAM machine.