Droplet-Laden Flows in Multi-Stage Compressors an Overview of the Impact of Modeling Depth on Calculated Compressor Performance

Silvio Geist; Markus Schatz

doi:10.20944/preprints202506.0564.v1

Submitted:

06 June 2025

Posted:

09 June 2025

You are already at the latest version

Abstract

There are various mechanisms through which water droplets can be present in compressor flows, e.g. rain ingestion in aeroengines or overspray fogging used in heavy duty gas turbines to boost power output. For the latter, droplet evaporation within the compressor leads to a cooling of the flow as well as to a shift in the fluid properties which is beneficial to the overall process. However, due to their inertia, the majority of droplets is deposited in the first stages of a multistage compressor. While this phenomenon is generally considered in CFD-computations of droplet-laden flows, the subsequent re-entrainment of collected water, the formation of new droplets and the impact on the overall evaporation is mostly neglected because of the additional computational effort required, especially with regard to the modeling of films formed by the deposited water. The work presented here shows an approach which allows to integrate the process of droplet deposition and re-entrainment based on relatively simple correlations and experimental observations from literature. Thus, the two-phase flow in multistage compressors can be modelled and analyzed very efficiently. In this paper, the models and assumptions used are described first, then the results of a study performed based on a generic multistage compressor are presented, whereby the various models are integrated step by step to allow an assessment of their impact on the droplet evaporation throughout the compressor and overall performance. It can be shown that evaporation becomes largely independent of droplet size when deposition on both rotor and stator and subsequent re-entrainment of collected water is considered. In addition, open issues with regard to a future improvement of models and correlations of two-phase flow phenomena are highlighted based on the results of the current investigation.

Keywords:

multi-phase flow

;

compressors

;

CFD

Subject:

Engineering - Energy and Fuel Technology

Introduction

An effective and simple method to increase the power output of gas turbines is to spray a fine water mist into the compressor intake to cool the incoming airflow. In the course of this process, the inlet flow saturates, thus shifting evaporation towards the rear stages of the compressor which further boosts the process. Among the various 1D programs used to evaluate the impact of wet compression on compressor behavior and performance, only a few account for non-equilibrium momentum transfer [1,2]. White and Meacock [1] note that droplets larger than 5 µm in diameter experience a slip velocity when the gas flow is redirected within the blades, introducing additional two-phase flow dynamics. Droplets generated from current spray technology may reach sizes of 30 µm or above. As the evaporative time scale of 30 µm-droplets is relatively large, they may, in principle, travel eight or more stages deep into the compressor [3]. Due to the velocity slip, such droplets largely impact on the blades reducing the water mass available within the flow. The resulting films, their path through the machine and a possible re-entrainment of droplets into the gas path might then lead to a very different evaporation and intercooling behavior. However, CFD-studies often either do not model any of the inherently three-dimensional mechanisms described above, model only isolated phenomena or do consider only single stage or even single row setups so that the true impact of evaporating flows on compressor performance remains unclear.

This shortcoming was therefore addressed in [4] by the authors. The study aimed to provide a benchmark case for comparing different multiphase CFD-models, specifically focusing on droplet behavior, deposition and the effect of wet conditions on aerodynamic performance.

The main findings from the first paper were that:

The majority of droplets, typical for high fogging in terms of diameter hit the blading in the first stages. (at least 25 % and up to 98 % of the injected water mass for droplet diameters between 5 and 30 µm).
Droplet-wall interaction modeling is essential because deposition effectively decreases droplet water mass available for evaporation, while
Partial deposition (splashing) was shown to be particularly relevant for intermediate and large droplets, as the resulting secondary droplets strongly enhance evaporation and were identified to be the main source of phase change for those cases.

The paper at hand advances the previous approach by extending it with simple approaches to account for liquid film behavior and disintegration as illustrated in Figure 1. Liquid films in compressors can form on three key surfaces: the rotor, stator, and casing. These films can potentially reduce aerodynamic efficiency by increasing drag and flow resistance on the blades. Additionally, their presence can modify the heat transfer characteristics of the blades. Lastly and most important, these films transport deposited liquid towards the blade edges, thus determining the re-entrainment of the liquid back into the gas flow. White and Meacock also highlight a limitation of their 1D approach, noting that in reality vapor generated from casing film evaporation would remain confined within the boundary layers near the casing, rather than mixing completely across the blade height, as assumed in their model. This underscores the need for more detailed three-dimensional simulations to accurately capture the complex flow behavior and phase change dynamics within the compressor.

Multi-Phase Flow Modeling

The multiphase problem is modeled using an Eulerian framework for the gas phase and a Lagrangian approach for the droplet flow, allowing for the exchange of heat, momentum, and mass between the two phases. For further details about the CFD-settings and boundary conditions, the interested reader is referred to [4].

Gas Phase

The working fluid, humid air, is modeled as an ideal mixture (denoted by index g) of two perfect gases: dry air (index a) and water vapor (index v). It is presumed that thermodynamic properties at a certain location are solely dependent on the local mass fractions so that

Y^{*} = Y_{a} + w Y_{v}

(1)

holds for the mixture state variable

Y^{*}

e.g. heat capacity at constant pressure

c_{p}^{*}

, heat capacity ratio

γ^{*}

or gas constant

R^{*}

. The mass flow fraction

w = {\dot{m}}_{v} ⁄ {\dot{m}}_{a}

(2)

is the ratio of water vapor and dry air.

According to Dalton’s law, the total pressure of the mixture

p_{g} = p_{a} + p_{v}

(3)

is the sum of the component partial pressures

p_{a}

and

p_{v}

.

Droplet Phase

For simplicity, droplets are assumed to be spherical with a uniform temperature throughout their volume. They are grouped into parcels (denoted by index P), each representing the mass of several hundred individual droplets. All droplets within a parcel are assumed to have identical diameters and thermodynamic properties. The trajectory of each parcel through the gas field is calculated using a forward-stepping method that satisfies the force balance between parcel momentum and drag force, modeled by the Schiller and Naumann correlation [6]. In the rotor, additional forces, such as Coriolis and centrifugal forces, arising from the rotational domain, are also considered.

Parcels are either introduced at the compressor inlet with predefined diameter, temperature, and a velocity matching that of the gas, or they originate from splashing and re-entrainment events with respective properties. Each parcel is tracked until one of the following conditions is met:

the parcel exits the computational domain through the compressor outlet,
the parcel impacts a surface,
the droplet diameter decreases below $1 e - 09 m$ , or
the tracking time exceeds 20 seconds.

Heat Transfer

The amount of convective heat transferred is driven by the temperature difference

T_{g} - T_{P}

between the gas and the droplet while the thermal situation at the drop surface is accounted for by the Nusselt number. An empirical formulation by Ranz and Marshall [7] was used in this study which has proven to produce adequate results [8] for spherical droplets. Heat transfer resulting from the phase change is coupled directly to mass transfer via a constant latent heat of

L = 2256 k J / k g

equal to the difference of the enthalpies of saturated vapor and water at

1 b a r

.

Mass Transfer

The mass transfer is computed using the liquid evaporation model outlined in the CFX Solver Theory Guide [9]. A droplet parcel is considered to evaporate when its saturation pressure,

p_{P, s a t}

, is lower than the mixture pressure,

p_{g} (p_{P, s a t} < p_{g}

). Mass transfer occurs via diffusion, with the process influenced by the droplet surface conditions and the surrounding vapor layer. These effects are captured through the Sherwood number, with a correlation also proposed by Ranz and Marshall [10]. These heat and mass transfer models are already integrated into the latest version of the CFX software. The saturation pressure

p_{P, s a t}

primarily depends on the parcel temperature

T_{P}

and is determined using the semi-empirical Antoine equation

p_{P, s a t} = e x p (A - \frac{B}{(T_{P} + C)})

(4)

where

A, B

and

C

are substance-specific coefficients. For water, Rodgers and Hill [11] have determined the coefficients to be

A = 11.78, B = 3885.7 K

and

C = - 42.92 K

, applicable within the temperature range of

268.15 K < T < 408.15 K

. The error of this correlation, when compared to measured values, is

1 %

or less.

Droplet-Wall Interaction

The model specifically developed for compressor applications by Gomaa [12] is employed in this work. Similar to most other models, it incorporates the non-dimensional impact Weber number

{W e}_{I} = ρ_{P} D_{I} C_{I}^{2} ⁄ σ_{P}

(5)

as well as the impact Reynolds number

{R e}_{I} = ρ_{P} D_{I} C_{I} ⁄ μ_{P}

(6)

and the characteristic impingement parameter

K_{C}

[13].

As discussed in [4] and demonstrated by Cossali et al. [14] liquid water films, characterized by their physical thickness

h_{I}

, have no significant effect on the splashing process if the dimensionless film height in the impact zone is below or equal to unity

H_{I} = h_{I} ⁄ D_{I} \leq 1

(7)

Based on light interferometry measurements performed by Neupert et al. [15] in a compressor cascade, this condition is violated only when very small droplets impact near the leading edge. However, these droplets do not possess enough energy upon impact to cause splashing and are more likely to deposit. Therefore, it is sufficient to treat every impact as a dry splash, with the threshold set to its minimal value of

K_{C} = 85

according to [12], which also simplifies the modeling process.

When a droplet impacts the wall, as depicted in Figure 2, the model calculates a modified Weber number:

{\tilde{W e}}_{I} = {W e}_{I} {(K_{C}^{4} \sqrt{{W e}_{I}} / {R e}_{I})}^{- 2 / 5}

(8)

based on the impact parameters to determine the resulting regime. According to the method of Gomaa [12], the droplet fully deposits on the wall if

{\tilde{(W e)}}_{I} \leq 1

and the entire parcel mass is removed from the droplet phase and stored in a memory location corresponding to the blade row where the impact occurred. For

{\tilde{(W e)}}_{I} > 1

, depending on the impact strength, a fraction of the impacting liquid mass is ejected back as two classes of secondary droplets, with diameters:

D_{1} = 0.07 D_{I}

(9)

D_{2} = (0.07 + 0.04 / \sqrt{{\tilde{W e}}_{I}}) D_{I}

(10)

where

D_{I}

is the droplet diameter upon impact [12]. Table 1 shows examples of relative ejected mass according as a function of

{\tilde{W e}}_{I}

. For high-impact events with

{\tilde{W e}}_{I} > 50

, all secondary droplets are classified as

D_{1}

. The direction and velocity of the secondary droplets are not explicitly modeled; instead, it is assumed that these small droplets quickly adjust to the gas flow velocity and direction.

Films

Droplets primarily impact the blading on the leading edge and pressure side. While the surface structure of a water film definitely alters the blade surface, comparably to a change in surface roughness, this effect appears to have no significant overall impact on aerodynamic losses. Harbeck et al. [16] conducted probe measurements in a linear compressor cascade, comparing total pressure along a traverse behind the blade passage with and without water injection. No differences were observed, and the loss coefficient remained unchanged in both cases. This result aligns with the findings of Roumeliotis and Mathioudakis [17] in a compressor stage, and Day et al. [5] in a multi-stage compressor. In these studies, probe measurements were taken for a range of relative water loads. For high fogging operation, i.e. water loads below

4 %

of the dry air mass flow, no differences were found in flow deflection or total pressure loss coefficients compared to dry flow field measurements. As a result, aerodynamic losses are not considered in the modeling for this work. Instead, the focus is placed on the transport of liquid mass.

Numerical investigations of stator profiles, based on flat plate experiments [12,15,18], show that shear forces carry the liquid completely to the blade trailing edge. Similarly, results obtained for rotor film modeling studies [19,20,21] illustrate that water deposited on rotating blades is thrown almost fully radially towards the casing by centrifugal forces. Based on these observations, computationally intensive thin film modeling for each blade row to describe liquid movement on the blade surfaces is bypassed in the present work by simply transferring the deposited water mass to either the stator trailing edge or the casing above the rotor.

The liquid is then reintroduced into the flow in the form of droplets formed by disintegration of ligaments at a different location within the computational domain through injectors. Thus, the present work does not account for any transfer of momentum or heat between the film surface and the gas mixture, nor the mass transfer through evaporation or boiling of the film. However, these effects are of minor importance compared to droplet evaporation as will be elaborated in the next section.

Casing Films:

Quantitative measurements of film thickness on the casing are not available. Moreover, the film behavior on the casing is highly complex due to gas velocity gradients, potential droplet impacts and the change of temperature and pressure along the compressor stages. Accurate prediction requires spatial and temporal modeling approaches, especially for a correct description of mass and heat transfer to assess the casing film contribution to vapor formation and cooling through simulation.

In [22], photographs and videos of the rotor gap area of a compressor were taken through a transparent casing. When water droplets are introduced into the flow, a casing film over the rotor gap is clearly visible, moving almost entirely in the circumferential direction. As water loading increases, the film structure becomes irregular, spreading more on the casing and forming waves on its surface.

Matz et al. [2] estimate that

90 %

of compressor cooling occurs through droplet evaporation, emphasizing its dominant role compared to casing film. While the film surface is orders of magnitude smaller compared to the total surface area of all droplets, the residence time of the casing film in the compressor is significantly longer than the flight time of the droplets. Hill [23] argues that this relationship is crucial for high-fogging and that droplet evaporation plays a secondary role compared to film evaporation.

Zhluktov [24] and White and Meacock [25] modeled casing film in 1D compressor simulations, assuming liquid deposited on stators flows to the hub and casing, while liquid on rotors is fully allocated to the casing. Heat transfer depends on the surface-to-gas temperature gradient and boundary layer thickness, with different Reynolds number correlations for the casing and blades. Zhluktov’s model attributes

50 %

of the total vapor formation to the film, when assuming an adiabatic casing wall. White and Meacock compared cases with and without film heat and mass transfer modeling. They found that when the wall temperature matches the vapor saturation temperature, film cooling has little effect on the compressor exit temperature or vapor development, with results comparable to those without film evaporation modeling. However, if the film temperature exceeds the vapor saturation temperature, the vapor mass flow increases, reducing the compressor exit temperature significantly. Hence, the casing film's long residence time is not enough to cause a cooling effect, it also needs to have a corresponding temperature. Consequently, any accurate assessment of film evaporation requires the knowledge of the casing wall temperature distribution.

Due to the lack of knowledge on the behaviour of casing films as well as the uncertainty regarding its impact on overall evaporation, it will not be considered in the present work.

Film Re-Entrainment:

For the consideration of water re-entrainment at the stator trailing edge, studies based on flat plates ([12,27,28,29]) or airfoils ([30,31,32,33,34,35,36,37]) are particularly suitable, as experimental campaigns in turbomachines are significantly more complex and therefore less common [38]. In all studies, researchers observe a breakup process consisting of up to three phases, regardless of whether the liquid is present as a film or a rivulet. The liquid at the trailing edge grows into the surrounding gas flow in several finger-like water structures (ligaments) and is increasingly exposed to aerodynamic forces due to the velocity difference between the liquid and the gas. In the primary zone, these forces exceed the viscous forces within the ligament, causing it to break up into individual droplets (Figure 3). For blunt trailing edges, the liquid accumulates in the wake, forming larger primary droplets, which also undergo breakup in a secondary zone.

Calculations of the deformation Weber number [35], as well as high-speed imaging of this secondary breakup ([28,30,31,33]), demonstrate that this involves a bag breakup mechanism, where the ligament is inflated and torn apart by the gas flow. Melekidis [29] and Schlottke and Weigand [28] experimentally determined the spatial dispersion of the droplet spray downstream of the trailing edge. In both experimental setups, a single rivulet was placed on a simplified guide vane profile [29] or a plate [28], and its breakup process was captured in both top and side views using high-speed photography. From the envelope of the area covered by the droplets, Melekidis derived a cone angle to describe their spatial distribution.

In Figure 3, the cone angle is denoted by

θ

and represents the lateral spread of the spray in the context of axial compressors. The spray cone becomes larger on average at higher gas velocities because the oscillation of the ligaments during bag breakup increases. In [28], the droplet volumetric flow rate was evaluated at a plane downstream of the trailing edge. The distance from the trailing edge was chosen such that the secondary droplet breakup had fully developed. The results show that the spread of the droplet spray perpendicular to the trailing edge increases significantly with increasing gas velocity, consistent with the observations in [29]. Thus, on a fully wetted stator blade, disintegration can occur for the entire trailing edge, and the disintegrations overlap, resulting in a basically uniform distribution of secondary droplets along the full blade span.

The trailing edge disintegration model must describe the characteristics of the spray following ligament breakup and possibly secondary droplet breakup. Ansys CFX provides "Particle Injection Regions" (see Chapter 8.8 of the Ansys Modeling Guide [39]), hereafter referred to as "injectors", to achieve this. These injectors can be freely positioned within the computational domain and introduce a defined spray into the gas flow and will be elaborated in section 4.2. To calculate the final droplet diameters downstream of the trailing edge, a correlation from Esquivias et al. [31] is used:

D_{T E} = 3.3 {(1 / {W e}_{T E})}^{0.72} h_{T E}

(11)

where

h_{T E}

is the trailing edge thickness. This correlation provides the Sauter mean diameter (

D_{32}

) of the spray generated for both sharp and blunt trailing edges, and aligns well with experimental results for velocities at the trailing edge (

C_{T E}

) between

50

and

175 m / s

, which is well within the range encountered in axial compressors. The trailing edge Weber number based on the trailing edge thickness

{W e}_{T E} = ρ_{g} C_{T E}^{2} h_{T E} / σ_{P}

(12)

is the key parameter. Additionally, a correlation for the ligament length

L_{b} = (1.7 - 0.0021 {W e}_{T E}) h_{T E}

(13)

also based on the trailing edge Weber number was proposed.

The ligament detachment frequency is neglected, as only steady-state simulations were performed.

In the CFD simulations, the procedure is such that the liquid mass flow absorbed by the wall-droplet interaction model is calculated for each blade, casing, and hub using multiple junction box routines at the beginning and end of each iteration step. A user-defined function links the deposited mass flow on a surface to the mass flow of a specified injector. As schematically illustrated in Figure 4, an injector is placed in the flow field behind the stationary blade, and the deposited mass on the blade is linked to the injector at the following time step.

Benchmark Compressor

The compressor geometry utilized for the numerical experiments is the LSM Benchmark Compressor published in [4], which features ten stages and an Inlet Guide Vane (IGV). Data on blade curves, grid files, and case boundary conditions for the LSM Benchmark Compressor are available upon request. The aim is to establish a comprehensive database of results obtained with various multi-phase models and methodologies, all using the same geometry. Key parameters from a simulation at the main operating point (OP) of the compressor – using saturated air with a vapor content of

w_{i n} = 6.11 ‰

– are summarized in Table 2.

This operating point is also marked on the compressor performance map shown in Figure 5. The simulation results indicate that the compressor achieves an overall total pressure ratio

Π = 7.02

with an inlet mass flow rate of

{\dot{m}}_{i n} = 157.9 k g / s

.

When calculating the specific compression work

a_{c}

, it is important to consider not only the contribution from the air but also from the water vapor and liquid droplets to accurately assess the High-Fogging process. White and Meacock [19] suggest the following approach:

a_{c} = (c_{p, a} + w c_{p, v}) ∆ T_{t, g} + f_{P} c_{p, l} ∆ T_{P} + L ∆ w

(14)

where the first term is the total enthalpy change of the moist air, the second term is the sensible heat change of the liquid fraction and third term is the latent heat contribution from evaporation. Note that this is a purely thermodynamic approach. Any mechanical loss due to droplet impact or due to centrifuging of water is not considered here.

The polytropic efficiency

η_{P o l}

is generally defined as the ratio of reversible enthalpy change

v d p

to the overall specific enthalpy change

d h

of the flow. Through

v d p = d h - T d s

(15)

the specific entropy change ds can be determined. This is approximated by:

Δ s_{12} = Δ s_{a, 12} + {w_{2} s}_{v, 2} - w_{1} s_{v, 1} + {f_{2} s}_{l, 2} - f_{1} s_{l, 1} + Δ s_{m i x}

(16)

between stations 1 and 2 where

s_{a}, s_{v}, s_{l}

are the specific entropies of the air, vapor and liquid phase. Further details on these terms can be found in the appendix. The last term

Δ s_{m i x}

describes the change of entropy due to the mixing of the phases. To calculate this fraction, Baehr and Kabelac [40] suggest the following approach:

Δ s_{m i x} = R_{v} [(\frac{R_{a}}{R_{v}} + w) l n (\frac{R_{a}}{R_{v}} + w)] - R_{v} [w l n (w) + \frac{R_{a}}{R_{v}} l n (\frac{R_{a}}{R_{v}})]

(17)

Note that

Δ s_{m i x}

does neither depend on pressure nor temperature but increases with the vapor fraction

w

.

R_{a}

and

R_{v}

are the specific gas constants of air and water vapor.

Next, in Equation (15) the specific enthalpy change

d h

is equated to the specific compressor work

a_{c}

from Equation (14). Together with Equation (16) this leads to the following expression for polytropic efficiency:

η_{P o l} = v d p / a_{c} = 1 - T_{o u t} Δ s / a_{c}

(18)

This formulation is valid for mixtures of arbitrary composition, including air, vapor, and droplets. The polytropic efficiency for the main operating point OP without inlet fogging calculated from Equation (18) is

η_{P o l} = 0.904

, as listed in Table 2.

Modeling Approach

Assumptions and Simplifications

To reduce the number of free parameters and computational effort, a number of assumptions and simplifications have been made. The working fluid is assumed to be an ideal and homogeneous mixture of air and water vapor, both of which have been modelled as perfect gases. Vapor diffusivity has been assumed to be independent of pressure and temperature. For liquid water, all properties including latent heat and surface tension were set to be constant. The impact of these assumptions on the accuracy of the results is discussed at the end of this paper.

Cases Studied

In order to assess the impact of the various two-phase flow phenomena discussed in the previous sections, several cases with increasing level of complexity have been simulated:

Case 1: Dry Flow:

The converged solution of the gas flow field without inlet fogging, i.e. without droplet injection at operating point OP was obtained.

Case 2: Full Deposition:

A combined Eulerian/Lagrangian simulation was performed for which the flow field computed in the dry case was used as the initial solution. Monomodal sprays in the form of 50,000 droplet parcels, with diameters of

5, 15

or

30 µ m

, were released at the compressor inlet and tracked through the flow field, representing the fine, intermediate, and the large mass-carrying, portions of the spray, respectively.

All parcels at the inlet are initialized without slip relative to the gas flow and with a temperature corresponding to the gas total temperature at the compressor inlet. The total mass of all packages corresponds to a water loading of

1.07 %

of the air flow without fogging (Case 1, see Table 2). Droplets hitting any surface were assumed to be deposited completely and hence removed from the flow, while the impact information was written to an output list for later evaluation. Source terms between the phases were calculated every 5th gas phase iteration and were applied using a relaxation factor up until the next coupling of the phases and until 400 gas phase iterations were reached.

Case 3: Partial Deposition:

Case 2 was repeated with the splashing model of Gomaa [12]. Thus, if a droplet hits a surface, the splashing model determines the deposited mass, and secondary droplet parcels are tracked in the next iteration. Deposited water mass is removed from the flow.

Case 4: Partial Disintegration:

Five injectors were positioned behind the IGV and the stator of each stage to re-entrain the water mass deposited on the latter with a monomodal spray of given diameter DTE, using the disintegration model presented in the next section. Figure 6 shows an illustration of that approach for the IGV and corresponding droplet trajectories originating from the injectors until they enter the rotor passage.

The simulation of the disintegration process within the rotor gap is omitted. Instead, the spray angle of the injectors in each stage was chosen such that uniform coverage of the blade height with droplet parcels was achieved at the downstream mixing plane with all five injectors activated. The droplet velocity was set to the average axial gas velocity at the downstream mixing plane.

For Case 4, any water mass deposited on the rotor is still thought to be centrifuged out towards the casing and subsequently removed from the control volume.

With regard to the droplets introduced at the compressor inlet, in addition to the monodisperse sprays with diameters of

5, 15

or

30

µm, a polydisperse spray was also analyzed during this stage. Figure 7 (top) presents the number distribution of a polydisperse spray used to augment stationary gas turbine compressors, as described in [41]. It was created by approximating the spray mass distribution (Figure 7 bottom) using the Rosin-Rammler function:

{R R}_{D} = f (D) = {e x p (- D / D_{R})}^{n_{R}}

(19)

with parameters

D_{R} = 26 µ m

and

n_{R} = 3.0

. This approach allows the discrete number distribution to be represented as a continuous function to the flow solver, as indicated by the black line in Figure 7. Consequently, a variety of droplet classes is injected that closely resemble the illustrated spray.

Case 5. Full Disintegration:

In this case the water deposited on the rotors is finally accounted for and it is fed back to the flow by an additional sixth injector placed behind the adjacent downstream stator (except the IGV) on the outermost normalized radius used for the injectors (corresponding to position 5 in Figure 6), because this strategy assumes that the casing water mass is completely drawn off by the next stator and disintegrated at the outer span of its trailing edge, as observed by Day et al. in [5].

To assess simulation quality, pressure ratio, temperature ratio and inlet mass airflow were monitored for every simulation, with final variations below

1 % .

Simulating Case 0 on a 16.4-million-cell grid required 14 GB of memory and 40 hours of computation. Including droplet trajectory calculations (Case 2 and 3) doubled the computation time, while accounting for disintegration (Case 4 and 5) increased both time and memory needs three- to fourfold.

Properties of the Injector

This section explains the adaptation of the injector model approach used to model re-entrainment of deposited water.

Representative Diameter

The Sauter mean diameter DTE of the secondary spray after re-entrainment is determined using Equation (11) by inserting the mass flow averaged mixture density and axial velocity in the absolute system at the mixing planes behind each stator. In order to assess the size range of the droplets formed by disintegration for a wide operating range of the compressor as shown in Figure 5 the results of simulations of operating points D, F, and H (see map in Figure 5) for inlet guide vane (IGV) positions −10°, +10°, and N were used.

The mean values and the corresponding standard deviations are plotted across the corresponding stages in Figure 8. It can be observed that

D_{T E}

is largest, approximately

55 µ m

, at the IGV. The diameter steadily decreases downstream and reaches a value of approximately

12 µ m

behind the last stator. Although the operating points used for the calculation cover the entire operating range of the compressor shown in Figure 5, the standard deviation for the droplet diameters formed due to disintegration at each stator is only between

2

and

4 µ m

. Therefore, in a first approximation, it is sufficient to use the average values obtained for each stage for the Sauter mean diameter of the droplets re-introduced into the flow by the injector. The relative uncertainty is greatest in the rear section, where the resulting diameters are small. A comparison with the experiments of Doerr et al. [38] shows that the diameters calculated from Equation (11) are valid in the context of axial compressors. In the publication, droplets in the range of

25 µ m

to

35 µ m

behind the second rotor of a 4.5-stage test compressor were measured using phase Doppler anemometry, with their origin explicitly attributed to the disintegration at the stator trailing edge.

Ligament Length

To ensure a proper positioning of the injectors, the ligament length

l_{b}

at each stator was estimated using Equation (13). If this length is larger than the distance to the subsequent rotor, the ligament would collide with the rotor itself, rendering the disintegration model invalid. Again, several operating points were considered to assess the range of ligament lengths across the whole performance map. The results are plotted in Figure 9. The trend shows a decreasing pattern, with the ligaments being estimated to be extremely short, approximately

0.6 m m

at the IGV and around

0.2 m m

at Stator 10. Between stages 2 and 5, no physically reasonable result was obtained because the validity range of the correlation (i.e.

{W e}_{T E} < 800

) was exceeded. Due to the overall downward trend, it is not expected that the ligaments for these stages are significantly longer than in the adjacent stages.

In comparison, the axial gaps between the cascades in each stage are significantly larger than the ligaments. Based on these results, the injectors were placed at a distance of

2 %

of the stator axial chord length in the direction of the stagger angle.

Droplet Temperature

This passage discusses the assumptions and limitations of modeling droplet temperature in an axial compressor simulation. In reality, the temperature of the disintegrated water droplets corresponds to the temperature of the thin water film present on the trailing edge of the stator blades. However, the approach used in the present work does not calculate this directly because simulating the heat transfer between the water film, the gas flow, and the blade surface would require a complex film model. Instead, the model assumes that the droplet temperature is equal to the average droplet temperature at the mixing plane in front of each stator, and that the heat transfer from and to the film is neglected. It is further assumed that the droplet temperature remains constant from the mixing plane to the point where they impact the stator. This approach simplifies the simulation but may underestimate the actual droplet temperature, because droplets could absorb additional heat between the mixing plane and the stator. As discussed before, the heat transfer from the water film might also affect the overall cooling of the gas.

Injector Position

In the current work, it is assumed that liquid re-entrainment at the trailing edge occurs at a similar blade span as the droplet impact along the blade. This has to be considered for the “Partial Disintegration” model (Case 4). In order to account for the change of the range of droplet impact depending on droplet size, an average radius of droplet impact was calculated for all droplet sizes by weighing the individual impact radii with the impacting droplet volume at each stator. The left-hand side of Figure 10 shows the droplet volume-weighted radius at the mixing planes in front of the stators for each droplet size considered. An averaged radius of 0.5 means that the droplets impact along the full span of the blade; larger averaged radii indicate that deposition is shifted towards the casing, thus the hub part of the blade is not impacted by droplets. From Figure 10 it is evident that the center of mass shifts for droplets of different sizes.

While for small droplets (

D_{P} = 5 µ m

) the center of mass largely remains at the mid-radius (

r_{r e d} = 0.5

), larger droplets are centrifuged outwards, shifting the center of mass towards the upper third (e.g. to

r_{r e d} ~ 0.6 - 0.7

for

D_{P} = 15 µ m

) or even further towards the casing after passing only a few stages. To account for this effect in the disintegration modeling, five virtual injector positions downstream of each stator are implemented, spaced evenly across the span. Based on this approach, five configurations are introduced.

Configuration 1 covers deposition and subsequent re-entrainment along the whole blade span; while in the following configurations injector locations are “switched off”, beginning with the hub location, until only the injector close to the casing is left in configuration 5. The resulting average radii of re-entrainment for this approximation is depicted in Figure 10 using gray symbols for intermediate and large droplets (square:

15 µ m

; triangle:

30 µ m

). It can be seen that the average radius of re-entrainment matches that of droplet deposition with good accuracy. For small droplets, only configuration 1 is used.

For the full disintegration case (Case 5), which considers the re-entrainment of droplets which deposit on the rotor blades and are then centrifuged towards the casing, an additional injector is placed at

83 %

span (position 5 in Figure 6).

Results

Overall Impact of Re-Entrainment Modeling

The following section analyzes how the different cases considered affect the simulation results of the flow with inlet fogging. Figure 11 shows the vapor fraction increase from compressor inlet to outlet,

Δ w = w_{o u t} - w_{i n}

, relative to the injected droplet liquid mass,

f_{i n}

. This ratio can be regarded as a measure of evaporation efficiency as it states which percentage of the injected water mass evaporates.

The upper plot shows that for a spray of the smallest droplet size,

D_{P} = 5 µ m

, evaporation efficiency is similar for all cases, i.e. re-entrainment and disintegration of deposited water is of minor importance for this droplet size. Starting from stage 7, the value of

Δ w ⁄ f_{i n}

remains constant for both the splashing and partial disintegration cases. For the full disintegration case, this value slightly increases.

For sprays of medium and large droplets or that featuring the Rosin-Rammler distribution, the full disintegration model (Case 5), re-injecting any deposited water back into the flow, significantly increases the evaporation efficiency from stage 4 onwards and vapor fraction increases linearly. Yet, the overall amount of evaporated water for these droplet classes with

Δ w ⁄ f_{i n} \approx 0.44

is smaller than for fine droplets with

Δ w ⁄ f_{i n} \geq 0.75

. For the partial disintegration case (Case 4), the effect is less pronounced compared to the simulation which only accounts for splashing (Case 3). Overall, the impact of the droplet size is minimal at the same modeling depth except for the smallest diameters, demonstrating that evaporation efficiency is actually more or less independent of droplet size but rather depends strongly on re-entrainment and disintegration for typical spray droplets and polydisperse sprays. Thus, the modeling approach plays a crucial role in the assessment of compressor inlet fogging.

Figure 12 depicts the ratio of droplet water mass fraction

f_{P}

to injected water mass

f_{i n}

to assess how much liquid in the form of droplets is present at a certain compressor stage.

When this parameter approaches zero, there are no droplets left in the gas flow, either due to evaporation or due to deposition. This absence of droplets is visible for the fine droplets in the splashing and partial disintegration case from stage 7 onwards. This is the reason why the vapor content (visible in the upper part of Figure 11) remains constant.

When examining the full disintegration case, there is a droplet population of about

10 %

of the injected mass flow which reaches the compressor outlet. Because these droplets originate from re-entrainment, which actually yields droplets with a larger diameter than

5 µ m

(see Figure 8) this has only a minor effect on evaporation efficiency as they are less agile to evaporate.

For the remaining droplet classes, it becomes evident once again that droplet size has little influence on the overall course of evaporation, yet the modeling approach significantly impacts the droplet mass fraction. While the droplet mass fraction increases only slightly for the partial disintegration model (Case 4), the behavior changes fundamentally for the full disintegration approach. Sawtooth-like patterns can be observed, with the minima marking the planes behind the rotors at which deposited droplets are centrifuged to the casing. The peaks correspond to the planes behind the stators, where the water mass deposited in the entire preceding stage is re-introduced into the gas flow. For Case 5, the fraction of water mass exiting the compressor is approximately

20 %

, while about

45 %

evaporates. The remainder is still removed from the flow. This is due to the fact that some of the relatively large droplets generated by the injectors near the casing are immediately re-deposited on the casing, especially in the rear stages.

Deposition can be assessed based on Figure 13. Again, modeling partial disintegration at the stator only (Case 4) yields only small changes compared to the case without re-entrainment (Case 3) because the majority of the relatively large droplets formed due to disintegration at the trailing edge is deposited on the following rotor blades and centrifuged towards the casing. The full disintegration model (Case 5), however, has an enormous impact on the amount of deposited water already in the front stages. Nevertheless, part of the liquid is always deposited at the casing, leading to a deposited mass fraction of

10 % - 35 %

at the compressor outlet, depending on the size of the droplets injected at the inlet.

Next, the key parameters specific compressor work

a_{C}

(Equation (14)) and polytropic efficiency

η_{p o l}

(Equation (18)) are utilized for further comparison. These quantities are plotted in Figure 14 as a percentage deviation from the values given for the case without fogging (Case 1) in Table 2. The size of the data point symbols is proportional to the diameter of the droplet class being analyzed. Simulations with the Rosin-Rammler spray are represented by squares. The colorbar illustrates the simulation case used.

Compression work reduces most, about

3 %

, when small droplets are injected. This also leads to a significant decrease in polytropic efficiency by at least

3 %

. A reduction of compression work and a simultaneous decrease of efficiency is not intuitive, but it has to be kept in mind that this is the polytropic efficiency obtained for the mixture. Figure 15 illustrates the contributions of each phase to the total specific entropy change, as calculated using Equation (16), for the smallest droplet size class. Evaporative cooling reduces the specific entropy of the air compared to the dry flow case (

Δ s_{c a s e 1} \approx 60 J ⁄ k g / K

). However, the processes of compression, heating, and vapor generation significantly increase the specific entropy of the vapor, resulting in a comparable overall entropy change. Additionally, the entropy change due to mixing (Equation (17)) further increases the total entropy, ultimately leading to lower mixture polytropic efficiencies.

When injecting larger droplets or a typical spray, data points cluster depending on modeling depth and lead to a similar decrease of specific work and efficiency. The full disintegration approach (Case 5) causes the largest deviations while the difference between the splashing and partial disintegration approach (Case 3 and Case 4) is marginal. The highest reduction of compression work is obtained for the smallest droplets, however, these constitute only a very small fraction of a typical spray.

The most realistic assessment of the potential benefits of high fogging for the LSM Benchmark compressor is obtained based on the two simulations using the Rosin-Rammler droplet distribution. The first simulation, based on a partial disintegration model, (Case 4) assumes disintegration at stators only, while the second, full disintegration model (Case 5) accounts for the re-entrainment of all but one-third of the injected water mass.

Spatial Distribution of Flow Quantities

The radial distributions of vapor fraction (Figure 16), temperature (Figure 17) and mass flow density (Figure 18) are analyzed as a function of model depth. The figures show the radial profiles behind the rotors of stages 1, 5, and 9 for all model depths and for the droplet size

D_{P} = 30 µ m

as these droplets are most likely to traverse towards the casing during their journey through the compressor. The plotted values are normalized to the results of the dry flow case (Case 1) with

w_{c a s e 1} = w_{i n} = 0.00611

.

Behind the first rotor, minor differences can be observed at

r_{r e d} > 0.9

, where higher vapor content due to increased evaporation and correspondingly lower temperatures, and higher mass flow density are noted especially for the full disintegration case (Case 5).

In later stages, more vapor is generated as re-entrainment is considered, and the profiles behind rotors R5 and R9 are strongly affected by the extent of re-entrainment modelled. In general, the shift of droplets towards the casing, or rather to the upper third of the annulus, leads to higher evaporation there, causing the mass flow rate to increase in this region. This also affects the flow temperature distribution due to uneven cooling. The strongest impact can be seen for the full disintegration case (Case 5) with the highest mass flow density in the upper annulus (

r_{r e d} > 0.6

) behind rotor R9, while mass flow near the hub (

r_{r e d} < 0.2

) decreases.

Accordingly, the aerodynamic load of the compressor stages is impaired based on that change of the radial mass flow, which also affects operating behavior.

Discussion of Modeling Assumptions

Finally, a critical assessment shall be made of the assumptions and simplifications used in the various approaches, as well as their impact on the results.

The approach chosen is absolutely sufficient for comparing various evaporation, droplet, and droplet-wall interaction models. However, the results show that a significant amount of liquid reaches the middle and rear stages of the compressor, where pressure and temperature have increased noticeably compared to the inlet conditions. Therefore, an estimate is made on how the rise in pressure and temperature affects splashing and trailing edge disintegration for these stages. To this end, the lower part of Figure 19 shows the evolution of static pressure and droplet temperature. The upper part shows the evolution of surface tension, latent heat and diffusivity based on correlations presented within this section and compared to the respective constant value set for the numerical simulations.

Surface tension

σ_{l}

directly influences the calculation of the impact Weber number (Equation (8)) and the trailing edge Weber number (Equation (12)). According to the Eötvös equation

σ_{l} = k ρ_{l} (T_{l}^{c} - T_{l}) / M_{l}

(20)

the surface tension of a liquid is linearly dependent on the liquid temperature

T_{l}

. This equation applies to all pure liquids and allows for the prediction of surface tension if the molar mass

M_{l}

, the density

ρ_{l}

, and the critical temperature

T_{l}^{c}

of the liquid are known. The universal Eötvös constant is

k = 2.10 {* 10}^{- 07} J / K / {m o l}^{2 / 3}

. For water droplets considered here, the surface tension is calculated as follows, with

ρ_{l} = 997.05 k g / m^{3}

,

M_{l} = 18.02 k g / k m o l

and

T_{l}^{C} = 647.1 K

as a function of the droplet temperature

T_{P}

. Droplet temperature and the surface tension relative to the compressor inlet are shown in Figure 19 along the stages. As the droplet temperature rises above

T_{P} = 360 K

, the surface tension at the compressor exit decreases to

85 %

of its initial value. This would intensify splashing in the rear section, where the droplet temperature is higher. The ligament length at the stator trailing edges would be shorter, and the resulting secondary droplets would be finer. Compared to the simulations conducted here, this would likely improve the cooling effect due to the increase in droplet surface.

Figure 19 further shows the latent heat L computed from the pressure data of the reference case without droplets, using the steam tables of XSteam [42]. Referenced to the assumed constant value of

L = 2265 k J ⁄ k g

, it becomes evident that less energy is required to evaporate the droplets in compressor stages further downstream of the inlet. This means that the same amount of water can induce a greater cooling effect in the front section of the compressor than in the rear. Therefore, without accounting for pressure-dependent phase change enthalpy, the cooling effect is overestimated in the cases discussed here.

Finally, the dependence of diffusivity on pressure and temperature will be analysed. Diffusivity has a direct impact on mass transfer and, consequently, on the vapor content in the mixture. Fuller et al. [43] provide an empirical correlation that determines the diffusivity of substance

v

in substance

a

as:

D_{v - a} = \frac{T^{1.75} {(1 / M_{v} + 1 / M_{a})}^{0.5}}{1000 p {(v_{v}^{1 / 3} + v_{a}^{1 / 3})}^{2}}

(21)

based on pressure and temperature. The parameters

v_{v}

and

v_{a}

are referred to as diffusion volumes. These values are determined using a nonlinear least squares method for common substance pairs. For the water vapor-air pair, the parameters are

v_{v} = 12.7

and

v_{a} = 20.1

. As droplet temperature increases, the water vapor generated at the droplet surface can diffuse more quickly into the surrounding gas mixture, enhancing mass transfer. However, the diffusion coefficient decreases as ambient pressure rises, because the vapor particles must move against this increased pressure. The effect of pressure is much more pronounced than that of temperature, leading to reduced evaporation in the middle and rear sections of the compressor.

Finally, the impact of large droplets on the rotor blades of a compressor will eventually lead to significant material erosion and hence to a drastic change in aerodynamic efficiency. This aspect has been omitted in the current work as it does not give insight to the actual influence of inlet fogging on the compression itself, however, it must be considered to evaluate the benefits and drawbacks of inlet fogging as a whole.

Conclusions

This paper has introduced a tool based on Ansys CFX for simulating high-fogging in compressors based on a generic, open-source geometry which is intended to serve as a benchmark for different models.

By assuming an adiabatic film flow and using a semi-empirical approach to model the water mass deposited on blades, suitable model parameters for trailing edge disintegration were derived. Extensive simulations, both with and without droplets, provided a comprehensive dataset for model validation. It is recommended that droplets larger than

5 µ m

be modelled with a splashing and trailing edge disintegration approach, as the fine secondary droplets generated during impact are crucial to capture the cooling effect. The trailing edge disintegration model effectively links the stators and rotors and intense droplet impact near the blade tips enhances evaporation, leading to a significant decrease in gas temperature and an increase in water vapor content along the compressor, especially near the casing. As a result, specific compression work reduces while the mixture polytropic efficiency of the process decreases. Based on the results presented, the benefit of droplet cooling alone lies within 1 and

1.5 %

for the LSM benchmark compressor examined.

However, as outlined during discussion of the results, elaborate modeling of the substance properties may lead to a different evaporation behavior especially in the later stages. Moreover, droplet impact erosion and its impact on efficiency has not been considered here.

Lastly, it would be feasible to integrate a model for heat and mass transfer due to evaporation or boiling of the film on the casing and blades into the proposed concept with manageable effort for further research. The flexibility in defining droplet injectors for re-entry modeling also enables detailed investigation of alternative injection methods, such as interstage fogging.

Author Contributions

Conceptualization, S.G. and M.S.; methodology, software, validation and formal analysis, S.G.; writing—original draft preparation, S.G.; writing—review and editing, M.S.; supervision, M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Geometry data and CFD-meshes for the LSM benchmark compressor are available upon request. Please contact M.S.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

Quantities
$C$	absolute velocity (m/s)	$T$	temperature (K)
$\dot{m}$	mass flow rate (kg/s)	$L$	latent heat (J/kg)
$p$	pressure (Pa)	R	specific gas constant (J/kg/K)
$ρ$	density (kg/m^3)	D	droplet diameter (m)
$σ$	surface tension (N/m)		diffusivity (m^2/s)
h	thickness (m)	$μ$	dynamic viscosity (Pa s)
$a_{c}$	specific compression work (J/kg)	$l_{b}$	ligament length (m)
r	radius (m)	n	rev per minute (1/min)
M	molar mass (g/mol)	$T_{l}^{c}$	critical temperature (K)
s	specific entropy (J/kg/K)	$k$	Eötvös constant ( $J / K / {m o l}^{2 / 3}$ )
Numbers
$H$	relative film thickness	$R e$	Reynolds number
$K_{C}$	Splashing threshold	$δ$	Deviation
$W e$	Weber number	$η_{P o l}$	polytropic efficiency
w	vapor mass per dry air mass	f	liquid mass per dry air mass
$Π$	total Pressure ratio	$τ$	total temperature ratio
${R R}_{D}$	Rosin-Rammler function
Sub-/Superscripts
$g$	ideal mixture	$l$	liquid phase
$a$	pure air	$P$	parcel
$v$	pure vapor	$32$	Sauter mean
t	total	TE	trailing Edge
sat	saturation	I	Impact
~	modified	in	inlet
out	outlet	red	normalized
ax	axial direction	depos	deposited mass
tr	triple point value	s	static

Appendix

Specific entropies for the phases are calculated using:

s_{a} = c_{p, a} l n (T / T_{t r}) - R_{a} l n (p / p_{t r})

(A1)

s_{v} = c_{p, a} l n (T / T_{t r}) - R_{a} l n (p / p_{t r}) + L (T_{t r}) / T_{t r}

(A2)

s_{l} = c_{p, l} l n (T_{P} / T_{t r})

(A3)

For this purpose, a reference state with the temperature

T_{t r}

and the pressure

p_{t r}

must be set in which the entropy of air and water vapor should be zero. For the liquid phase, the evaporation entropy for this state is also considered. The reference state used is the triple point of water, which has

T_{t r} = 273.16 K

and

p_{t r} = 611.657 P a

.

References

A.J. White and A. J. Meacock. An evaluation of the effects of water injection on compressor performance. J. Eng. Gas Turb. Power 126, 4 (2004), 145–158. [CrossRef]
C. Matz, G. Cataldi, W. Kappis, G. Mundinger, S. Bischoff, E. Helland and M. Ripken. Prediction of Evaporative Effects Within the Blading of an Industrial Axial Compressor. Journal of Turbomachinery 132, 4 (2010). [CrossRef]
J.R. Khan and T. Wang. Implementation of a non-equilibrium heat transfer model in stage-stacking scheme to investigate overspray fog cooling in compressors. International journal of thermal sciences 68 (2013), 63–78. [CrossRef]
S. Geist and M. Schatz. A generic multistage compressor design for multiphase code and model benchmark. Proc. ETC 15, Budapest, Hungary (2023). [CrossRef]
I. Day, J. Williams and C. Freeman. Rain Ingestion in Axial Flow Compressors at Part Speed. Journal of Turbomachinery 130 (1) (2008). [CrossRef]
Z. Naumann and L. Schiller. A drag coefficient correlation. Z. Ver. Deutsch. Ing 77, 318 (1935), 323.
W. Ranz and W. Marshall. Evaporation from droplets. Chem. Eng. Prog 48, 3 (1952), 141–146.
A. Seck, S. Geist, J. Harbeck, B. Weigand and F. Joos. Evaporation modeling of water droplets in a transonic compressor cascade under fogging conditions. International Journal of Turbomachinery, Propulsion and Power 5,1 (2020). [CrossRef]
Ansys Inc.. Ansys CFX Solver Theory Guide, R2021 R1. 2021.
W. Ranz and W. Marshall. Evaporation from droplets—Part 2. Chem. Eng. Prog., 48 (1952) 173–180.
R. Rodgers and G. Hill. Equations for vapour pressure versus temperature: derivation and use of the Antoine equation on a hand-held programmable calculator. British journal of anaesthesia 50 (5) (1978), 415–424. [CrossRef]
H. Gomaa. Modeling of liquid dynamics in spray laden compressor flows. Doctoral Thesis, University of Stuttgart (2014).
C. Mundo, M. Sommerfeld and C. Tropea. On the modeling of liquid sprays impinging on surfaces. Atomization and sprays 8, 6 (1998). [CrossRef]
G. E. Cossali, M. Marengo and M. Santini. Splashing characteristics of multiple and single drop impacts onto a thin liquid film. Proc. Int. Conf. Multiphase Flows (2007).
N. Neupert, J. Harbeck and F. Joos. An experimentally derived model to predict the water film in a compressor cascade with droplet laden flow. Journal of Engineering for Gas Turbines and Power 141, 9 (2019). [CrossRef]
J. Harbeck, S. Geist and M. Schatz. An Approach to Measure Total-Head in Wakes and Near End Walls at High Fogging Conditions”. Proc. ASME Turbo Expo (2021). [CrossRef]
I. Roumeliotis and K. Mathioudakis. Water Injection Effects on Compressor Stage Operation. Journal of Engineering for Gas Turbines and Power 129 (3) (2006), 778–784. [CrossRef]
H. J. Dohmen, F.-K. Benra, C. Schepers and S. Schuster. Analysis of the Size of Secondary Droplets due to Film Separation at the Blade Trailing Edge. Proc. ASME Fluids Engineering Division Summer Meeting, Chicago, IL, USA (2014).
J. Williams and J.B. Young. Movement of Deposited Water on Turbomachinery Rotor Blade Surfaces. Journal of Turbomachinery 129, 2 (2006), 394–403. [CrossRef]
T. Nikolaidis, P. Pilidis, J. A. Teixeira and V. Pachidis. Water Film Formation on an Axial Flow Compressor Rotor Blade. Proc. ASME Turbo Expo, Berlin, Germany (2008), 79–87.
T. Nikolaidis and P. Pilidis. The Effect of Water Ingestion on an Axial Flow Compressor Performance. Proc. Inst. Mech. Eng., Part G: Journal of Aerospace Engineering 228 (3) (2014), 411– 423. [CrossRef]
L. Liu, H. Zhang, J. Li, C. Yu, F. Lin and C. Nie. Measurements and Visualization of Process from Steady State to Stall in an Axial Compressor with Water Ingestion. Proc. ASME Turbo Expo San Antonio, TX, USA (2013).
P. G. Hill. Aerodynamic and thermodynamic effects of coolant injection on axial compressors. The Aeronautical Quarterly 15, 4 (1963), 331–348. [CrossRef]
S. Zhluktov, S. Bram and J. de Ruyck. Injection of Water Droplets in an Axial Compressor. Proc. 5th World Conference on Experimental Heat Transfer, Fluid Mechanics and Thermodynamics. Thessaloniki, Greece (2001), 1415–1420.
A. J. White and A.J. Meacock. Wet compression analysis including velocity slip effects. Journal of engineering for gas turbines and power 133, 8 (2011). [CrossRef]
J. Loebig, B. Vittal and M. Booher. Numerical Simulation of Water/Methanol Evaporation in an Axial Flow Gas Turbine Compressor. 34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit. Cleveland, OH, USA (1998). [CrossRef]
T. Sattelmayer and S. Wittig. Internal Flow Effects in Prefilming Airblast Atomizers: Mechanisms of Atomization and Droplet Spectra. Journal of Engineering for Gas Turbines and Power 108 (3) (1986), 465–472. [CrossRef]
A. Schlottke and B. Weigand. Two-Phase Flow Phenomena in Gas Turbine Compressors with a Focus on Experimental Investigation of Trailing Edge Disintegration. Aerospace 8 (4) (2021), 91. [CrossRef]
S. Melekidis. Untersuchung der Tropfenentstehung an Schaufelhinterkanten. Master's Thesis, University of Stuttgart, (2017).
G. Chaussonnet, S. Gepperth, S. Holz, R. Koch and H.-J. Bauer. Influence of the Ambient Pressure on the Liquid Accumulation and on the Primary Spray in Prefilming Airblast Atomization. International Journal of Multiphase Flow 125 (2020). [CrossRef]
B. Esquivias, B. Hickey and V. McDonell. Investigation of Water Films Shed from an Airfoil in a High-Speed Flow. Proc. ASME Turbo Expo. Rotterdam, Netherlands (2022), GT2022-80350. [CrossRef]
S. Gepperth, D. Guildenbecher, R. Koch and H.-J. Bauer. Pre-filming Primary Atomization: Experiments and Modeling. 23rd European Conference on Liquid Atomization and Spray Systems (ILASS-Europe). Brno, Cech Republik (2010).
S. Gepperth, A. Müller, R. Koch and H. Bauer. Ligament and Droplet Characteristics in Prefilming Airblast Atomization. ICLASS, 12th Triennial International Annual Conference on Liquid Atomization and Spray Systems. Heidelberg, Germany (2012).
S. Gepperth, R. Koch and H.-J. Bauer. Analysis and Comparison of Primary Droplet Characteristics in the Near Field of a Prefilming Airblast Atomizer. Proc. ASME Turbo Expo, San Antonio, TX, USA (2013).
T. Inamura, N. Katagata, H. Nishikawa, T. Okabe and K. Fumoto. Effects of Prefilmer Edge Thickness on Spray Characteristics in Prefilming Airblast Atomization. International Journal of Multiphase Flow 121 (2019). [CrossRef]
B. Javed, T. Watanabe, T. Himeno and S. Uzawa. Effect of Trailing Edge Size on the Droplets Size Distribution Downstream of the Blade. Journal of Thermal Science and Technology 12 (2017). [CrossRef]
B. Javed, T. Watanabe, T. Himeno and S. Uzawa. Experimental Investigation of Droplets Characteristics after the Trailing Edge at Different Angle of Attack. International Journal of Gas Turbine, Propulsion and Power Systems 9 (3) (2017), 32–42. [CrossRef]
T. Doerr, S. Schuster and D. Brillert. PDA Laser Measurements of Droplet Laden Flows in a Four-Stage Axial Compressor. Journal of Turbomachinery 145 (2) (2022). [CrossRef]
Ansys inc. Ansys CFX Solver Modeling Guide, R2021 R1. (2021).
H. D. Baehr and S. Kabelac. Gas-Dampf-Gemische. Feuchte Luft. In: Thermodynamik: Grundlagen und technische Anwendungen. Springer, Berlin Heidelberg, (2012).
C. B. Meher-Homji and T. R. Mee. Gas Turbine Power Aug¬men¬tation by Fogging Of Inlet Air. Proc. 28th Turbomachinery Symposium. Houston, TX, USA (1999).
M. Holmgren. X Steam, Thermodynamic properties of water and steam. https: //www.mathworks.com/matlabcentral/fileexchange/ 9817-x-steamthermodynamic-properties-of-water-and-steam. MATLAB Central File Exchange. Last visited: 21.07.2024.
G. N. Fuller, P. D. Schettler und J. C. Giddings. New Method for Prediction of Binary Gas-phase Diffusion Coefficients. Industrial & Engineering Chemistry 58 (5) (1966), 18–27. [CrossRef]

Figure 1. Schematic representation of the water structures occurring in a compressor stage and the associated phenomena of drop impact, liquid films and liquid disintegration. The figure is adapted from [5].

Figure 2. Schematic representation of the drop-wall interaction, its regime and important parameters for the impact on a dry wall.

Figure 3. Mechanism of liquid re-entrainment on a thick trailing edge according to Inamura et al. [35].

Figure 4. Program sequence based on two consecutive iterations of the droplet phase.

Figure 5. Compressor Map without inlet fogging (from [4])

Figure 6. Calculated parcel trajectories originating from the injectors behind the IGV.

Figure 7. Number distribution (top) and cumulate mass (bottom) of a droplet spray according to [41] and the corresponding Rosin-Rammler approximation.

Figure 8. Secondary drop diameters DTE formed by disintegration downstream of each compressor stage according to Equation (11).

Figure 9. Mean and standard deviation of ligament length downstream of each stator row according to Equation (13).

Figure 10. Droplet volume-weighted average dimensionless radius at the entrance of each stator passage (solid line, black) and corresponding average radius of re-entrainment (gray symbols). The illustrations on the right show the corresponding injector configuration.

Figure 11. Change of vapor fraction relative to injected water mass for different model cases. Operating point: OP as marked in Figure 5, water load:

1.07 %

Figure 11. Change of vapor fraction relative to injected water mass for different model cases. Operating point: OP as marked in Figure 5, water load:

1.07 %

Figure 12. Change in droplet mass relative to injected water mass for different model depths. Operating point: OP as marked in Figure 5, water load:

1.07 %

Figure 12. Change in droplet mass relative to injected water mass for different model depths. Operating point: OP as marked in Figure 5, water load:

1.07 %

Figure 13. change of deposited liquid mass to injected water mass for different model depths. Operating point: OP as marked in Figure 5, water load:

1.07 %

Figure 13. change of deposited liquid mass to injected water mass for different model depths. Operating point: OP as marked in Figure 5, water load:

1.07 %

Figure 14. Changes in specific work and polytropic efficiency for different modeling depths and droplet classes. Operating point: OP (see Figure 5), water load:

1.07 %

Figure 14. Changes in specific work and polytropic efficiency for different modeling depths and droplet classes. Operating point: OP (see Figure 5), water load:

1.07 %

Figure 15. Specific entropy change by phase for different model depths and droplet class

D_{P} = 5 µ m

. Operating point: OP (see Figure 5), water load:

1.07 %

Figure 15. Specific entropy change by phase for different model depths and droplet class

D_{P} = 5 µ m

. Operating point: OP (see Figure 5), water load:

1.07 %

Figure 16. Radial distribution of vapor fraction for different modeling approaches relative to Case 1. Operating point: OP (see Figure 5), water load:

1.07 %, D_{P} = 30 µ m

Figure 16. Radial distribution of vapor fraction for different modeling approaches relative to Case 1. Operating point: OP (see Figure 5), water load:

1.07 %, D_{P} = 30 µ m

Figure 17. Radial distribution of total temperature for different modeling approaches relative to Case 1. Operating point: OP (see Figure 5), water load:

1.07 %, D_{P} = 30 µ m

Figure 17. Radial distribution of total temperature for different modeling approaches relative to Case 1. Operating point: OP (see Figure 5), water load:

1.07 %, D_{P} = 30 µ m

Figure 18. Radial distribution of mass flow density for different modeling approaches relative to Case 1. Operating point: OP (see Figure 5), water load:

1.07 %, D_{P} = 30 µ m

Figure 18. Radial distribution of mass flow density for different modeling approaches relative to Case 1. Operating point: OP (see Figure 5), water load:

1.07 %, D_{P} = 30 µ m

Figure 19. Dependence of important model quantities of the multiphase modeling (top) on pressure and temperature (bottom) across the stages.

Table 1. overview of ejected droplet mass after impact.

${\tilde{W e}}_{I}$	1	5	10	25	50
rel. ejected mass	0,00	0,18	0,36	0,70	0,91

Table 2. Compressor performance at design point OP without inlet fogging.

${\dot{m}}_{i n}$ $[k g / s]$	τ [-]	n [1/min]	$η_{P o l}$ $[-]$	$w_{i n}$ $[‰]$
$157.9$	1.85	5000	0.904	6.11
$a_{c}$ $[k J / k g]$	Π [-]	$p_{t, i n}$ $[m b a r]$	$T_{t, i n}$ $[K]$	$f_{i n}$ $[%]$
247.6	7.02	1013.25	288.15	0.0

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Droplet-Laden Flows in Multi-Stage Compressors an Overview of the Impact of Modeling Depth on Calculated Compressor Performance

Abstract

Keywords:

Subject:

Introduction

Multi-Phase Flow Modeling

Gas Phase

Droplet Phase

Heat Transfer

Mass Transfer

Droplet-Wall Interaction

Films

Casing Films:

Film Re-Entrainment:

Benchmark Compressor

Modeling Approach

Assumptions and Simplifications

Cases Studied

Case 1: Dry Flow:

Case 2: Full Deposition:

Case 3: Partial Deposition:

Case 4: Partial Disintegration:

Case 5. Full Disintegration:

Properties of the Injector

Representative Diameter

Ligament Length

Droplet Temperature

Injector Position

Results

Overall Impact of Re-Entrainment Modeling

Spatial Distribution of Flow Quantities

Discussion of Modeling Assumptions

Conclusions

Author Contributions

Funding

Data Availability Statement

Conflicts of Interest

Nomenclature

Appendix

References

MDPI Initiatives

Important Links

Subscribe