Next Article in Journal
PV String-Level Isolated DC–DC Power Optimizer with Wide Voltage Range
Previous Article in Journal
Effect of LCA Data Sources on GBRS Reference Values: The Envelope of an Italian Passive House
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Heat Flux Distribution during Brush Seal Rubbing Using CFD with Porous Media Approach

Institute of Thermal Turbomachinery (ITS), Karlsruhe Institute of Technology, Kaiserstraße 12, 76131 Karlsruhe, Germany
*
Author to whom correspondence should be addressed.
Energies 2021, 14(7), 1888; https://doi.org/10.3390/en14071888
Submission received: 22 February 2021 / Revised: 22 March 2021 / Accepted: 25 March 2021 / Published: 29 March 2021
(This article belongs to the Special Issue Secondary Air Systems in Gas Turbines)

Abstract

:
This paper discusses the question of heat flux distribution between bristle package and rotor during a rubbing event. A three-dimensional Computational Fluid Dynamics (3D CFD) model of the brush seal test rig installed at the Institute of Thermal Turbomachinery (ITS) was created. The bristle package is modelled as a porous medium with local non-thermal equilibrium. The model is used to numerically recalculate experimentally conducted rub tests on the ITS test rig. The experimentally determined total frictional power loss serves as an input parameter to the numerical calculation. By means of statistical evaluation methods, the ma in influences on the heat flux distribution and the maximum temperature in the frictional contact are determined. The heat conductivity of the rotor material, the heat transfer coefficients at the bristles and the rubbing surface were identified as the dominant factors.

1. Introduction

Over the years, brush seals have been well established for use in stationary turbomachines and aircraft engines. The design processes are still based, to a large extent, on the manufacturers’ experience. The aim of the design process is to obtain a brush seal that is subject to minimal wear during operation, but it is still stiff enough to achieve the required pressure reduction in steady-state operation (see Ref. [1]). According to Aksit and Tichy [2], the wear of brush seals, similarly to other sliding contacts, depends on three ma in factors:
  • Material selection:
The choice of tribological partners determines, among other things, the coefficients of friction and wear properties. These are mainly defined by yield strength, modulus of elasticity, hardness, surface quality, and formation of oxide layers.
  • Contact mechanics:
The contact loads are of particular importance here. They are determined by the seal design, the operating parameters, and effects, such as blow-down, pressure stiffening, and hysteresis.
  • Heat transfer:
In particular, the temperature level in the friction contact has a direct effect on the mechanical and physical material properties of the sliding partners and the tribochemical reactions, such as the formation of oxide layers. The temperatures in the friction contact are determined by the level of the ambient air temperature and the heat input during the contact with the seal. These changes also indirectly affect the friction coefficient μ . In addition to the knowledge of the heat input in friction contact, the distribution of heat fluxes between the bristle package and the seal and, subsequently, the distribution within the bristle package are also important. In brush seals, these distributions are strongly influenced by cooling effects due to leakage. In the following, only the third factor, the heat transfer, will be discussed.
Already, in 1988, Gorelov et al. [3] recognized that a reduction of the leakage rate results in a significant heating of the brush seal. Accordingly, even small amounts of air are sufficient to cool the bristles adequately. The heat transfer within the package is comparable to that of a heat exchanger (see Ref. [4]). In this respect the diameter of the bristles and the cavities between the bristles play a decisive role.
With the help of a numerical flow simulation with modelling of the bristle package as a porous medium, Doğu and Aksit [5] calculated temperature profiles in the bristle package as a function of directly specified heat inputs in the friction contact. The maximum temperatures always occur at the bristle tips and decrease exponentially in radial direction to the level of the air inlet temperature. In the axial direction, a very uniform temperature distribution was determined. However, the authors assumed an isotropic heat conduction in the package. In reality, this is strongly anisotropic due to the orientation of the bristles and the bristle spaces. The highest radial temperature drop was found in the area of the back plate. The convective cooling effect increases due to higher leakage rates at high pressure differences. This, in turn, reduces the temperature level and increases the radial temperature gradient.
Demiroglu and Tichy [6] developed a semi-empirical equation in closed form to calculate the contact forces. This was used to calculate the friction power (see Ref. [7]):
P fric = µ · F N · u .
In order to validate the equation, rub tests were carried out and the temperatures at the rotor and at the downstream bristles were measured by infrared thermography. Because the part of the equation used is only valid without a pressure difference applied, no flow was applied in the validation measurements. The calculated frictional power served as an input variable for a finite element analysis. Therefore, the assumption was made that the heat flux between the rotor and seal is split in a ratio of 50:50. This was justified with similar heat conduction coefficients and an identical contact surface. However, further measurements have shown that this assumption is no longer valid under pressurisation. In this case, the heat flux distribution changes very much in favour of the bristles.
On the same test rig that was used by Demiroglu and Tichy [6], Ruggiero et al. [8,9] carried out steady-state, pressureless rub tests with brush seals with non-metallic bristles of aramid and carbon fibres. By means of a subsequent FE analysis, the heat input into the rotor was determined by iteratively adjusting this value until the calculated temperature gradients matched the measured ones. Following Demiroglu’s and Tichy’s hypothesis that, in the absence of flow through the seal, the heat flux distribution mainly depends on the thermal conductivity of the materials, Ruggiero et al. [8] assumed that, due to the low thermal conductivity of Kevlar, almost all of the heat should flow into the rotor.
Qiu and Li [10] used a numerical flow simulation to calculate the pressure forces. They did not resolve the individual bristles, but modelled the entire package as a porous medium. The contact forces were iteratively calculated with an FE model, taking into account the friction and deflection of the bristles. The frictional power was determined according to Equation (1). A constant friction coefficient of 0.24 was used. The resulting frictional power served as input for the Computational Fluid Dynamics (CFD) simulation. For the heat transfer in the bristle package, they assumed a local thermal equilibrium, i.e. the local air temperature corresponded to the bristle temperature. In a subsequent publication, the assumption of a local thermal equilibrium was replaced by an equation considering the convective heat transfer in the bristle package (see Ref. [11]). In order to determine the heat transfer coefficients at the bristles, they used correlations for banks of staggered tubes, as done by Doğu and Aksit [5]. The calculations showed an increase in the maximum temperatures in friction contact with increasing differential pressure, but a lower mean temperature in the package. With increasing speed, the maximum temperature and the mean temperature in the package also increased. As the temperature increased, a reduction in leakage was found. This is explained by a decreasing density of the streaming fluid.
Pfefferle et al. presented an instrumentation concept for the measurement of rotor temperatures during rubbing [12]. They placed thermocouples below a radially thin-walled rotor at six axial positions. For redundancy, the instrumentation was mirrored at a circumferential position offset by 180 . The thermocouple tips were welded to the rotor and the thermocouple lines were fixed with spot-welded sheets. In order to calculate the heat input into the rotor, an FE model was created and the heat input, as a boundary condition of the simulation, was iteratively varied until the best possible agreement between the experimental and numerical temperature curves was obtained. During the tests, the resulting total frictional power was calculated by subtracting the power measured before and during the rubbing. The difference between total frictional power and rotor heat input results in the heat input into the bristles and the air. Pfefferle [13] published the results of these tests. The tests were carried out with four seals of the clamped design, which are identical in design. Over all tests, it was calculated that 60–90% of the total frictional power goes into the rotor. Pfefferle [13] tried to model the transient behaviour of the experimental rubbing tests with a duration of 30 s in an additional simplified FE model in order to investigate these heat flux distributions in more detail. He was able to show that the transient test procedure leads to more heat being conducted into the rotor by up to 9% than would be the case in steady-state operation. A variation of the heat transfer coefficients showed that only for α b   <   100   W / ( m 2 K ) , more heat is conducted into the rotor. Pfefferle [13] explains this by an anisotropic arrangement of the bristles in axial and tangential direction, which results in little or no flow through parts of the bristle package (see also Ref. [14]). Pfefferle [13] believes that the worn bristle package itself provides a further explanation: agglomerated wear material at the bristle tips and severely deformed bristles reduce the leakage flow in the area in direct contact with the rotor. Both reasons lead to a reduction of convective heat transfer due to the reduced flow velocities. This is in contradiction to the correlations for banks of staggered tubes previously used for brush seals to calculate the heat transfer coefficients (see Refs. [11,15,16]).
The work by Pfefferle [13] was continued at the ITS in the subsequent years on a modified test rig (see Refs. [17,18,19,20]). The experimental investigations that were presented in these publications with subsequent calculation of the rotor heat input using an FE model demonstrated the influence of the geometry and operating parameters, as well as the degree of contamination of the bristle package on the heat flux distribution. When considering all of the experiments, it can be concluded that, in most cases, a large part of the frictional power generated is dissipated into the rotor.
The aim of this paper is to show the ma in factors influencing the heat flux distribution. The knowledge gained will be used to verify whether, as assumed by Pfefferle [13], the heat transfer coefficients at the bristles are significantly lower than that predicted by the correlations usually used for banks of staggered tubes. For this purpose, selected rub tests with seal 1 at different differential pressures are numerically simulated. The calculations include seal and rotor as well as the flow fields upstream and downstream of the seal. Hildebrandt previously published this numerical study (see Ref. [21]). The experimental results of the rub tests with this seal were published in [18].
Several approaches are available for modelling the flow in the bristle package. Besides the semi-empirical approaches (“bulk flow models”) and models with fully dissolved bristles, calculation approaches in which the bristle package is modelled as a porous medium have been established. When modelling the real geometry, the complete flow channels between the bristles within the package have to be meshed. Because of the small wire diameters in combination with high packing densities, this means an increased modelling and calculation effort. Furthermore, it is difficult to model the axially and radially variable bristle spacing according to a real brush seal. The approach of modelling the bristle packing by means of a porous medium overcomes the problems of this dynamic movement and the complex structure of the bristle package. For this reason, the following calculations are based on a modelling of the bristle package while using a porous medium.

2. Basic Equations for Calculating Flow Losses and Heat Transfer in Porous Media

The equations that are used to calculate the leakage flow of a brush seal are described below. Outside the bristle package, the usual Navier–Stokes equations are used. To account for the flow resistance that the bristles exert on the fluid, an additional resistance force F r per unit volume is defined within the bristle pack:
F r = A µ v B ρ | v | v .
This resistance is composed of a viscosity and an inertia term. Equation (2) is based on the results that were obtained by Darcy [22] and Forchheimer [23], and it is called the Forchheimer equation in its basic form. Because of the anisotropic composition of the bristle package, the flow losses are also directional. For example, the losses parallel to the longitudinal axis of the bristles are significantly lower than those transverse to the longitudinal axis of the bristles. However, they are largely responsible for the radial pressure gradient in the bristle package. As a result, the tensors A and B each consist of three elements:
A = a x 0 0 0 a s 0 0 0 a n , B = b x 0 0 0 b s 0 0 0 b n .
The coefficients are related to a local coordinate system, which is defined relative to the individual bristles. The ma in directions e x , e s , and e n are defined, as follows:
  • e x : parallel to the x-axis of the test rig,
  • e s : parallel to the bristle longitudinal axis, and
  • e n : perpendicular to x-axis and bristle longitudinal axis.
For the determination of the coefficients a i and b i , the findings from Ergun [24] can be applied, who investigated the flow through porous media in a Reynolds number range relevant to brush seals. According to Ergun [24], Equation (2) can be represented for the isotropic case as pressure loss over the height of the layer bed in the following way:
Δ p L = α ˜ · µ · ( 1 ε ) 2 d ˜ 2 · ε 3 · V + β ˜ · ρ · ( 1 ε ) d ˜ · ε 3 · V 2 .
The coefficients a i and b i are, thus, proportional to ( 1 ε ) 2 / ε 3 and ( 1 ε ) / ε 3 , respectively. The variable ε represents the porosity of the medium through which the medium flows. It is defined as the ratio of the void volume to its total volume:
ε = V f V tot = 1 V s V tot .
The porosity used for brush seals is described, in detail, in Section 2.2.
The diameter d ˜ shown in Equation (4) represents the equivalent sphere diameter:
d ˜ = 6 a v ,
a v = A V = d b · π · l 1 / 4 · d b 2 · π · l = 4 d b .
On the basis of measurements with packed beds of spherical particles, Ergun [24] suggests the factors α ˜ = 150 and β ˜ = 1.75 . With respect to the wire diameter d b of the brush seal, Equation (6) gives the factors α = α ˜ / d ˜ 2 · d b 2 = 66.7 and β = β ˜ / d ˜ 2 · d b 2 = 1.17 . In order to obtain a better match with the application under consideration of the brush seal, which resembles a cylinder packing, Chew and Hogg [25] and Proestler [26] suggest the factors α = 80 and β = 1.16 .
The coefficients a x and a n or b x and b n can be represented according to Equation (4), as follows:
a x = a n = α · ( 1 ε ) 2 d b 2 · ε 3 = 80 · ( 1 ε ) 2 d b 2 · ε 3 ,
b x = b n = β · ( 1 ε ) d b · ε 3 = 1 , 16 · ( 1 ε ) d b · ε 3 .
Chew et al. gave the loss coefficient a s along the bristle axis [27], based on a comparison of experimental and numerical results, with a ratio of 1:60 to the losses transverse to the bristle axis. In contrast, Proestler [26] calculated an analytical expression assuming the pressure reduction of a laminar flow in non-circular pipe sections. Proestler [26] calculated the factor a s , as:
a s = 32 · ε · ( 1 ε ) 2 d b 2 · ε 3 = 32 · ( 1 ε ) 2 d b 2 · ε 2 .
The dissipative losses along the bristle axis are neglected and the factor b s is set to zero. The numerical calculations performed in this paper use the loss coefficients according to Proestler [26].

2.1. Energy Conservation Equations

In general, two cases can be distinguished for heat transfer in porous media: The simple case of thermal equilibrium between the solid and the fluid and the case where there is a significant temperature difference between the two. In the case of a typical flow through the bristle package of a brush seal, strong forced convection occurs between the bristles and the leakage air, so that the second case applies. This is modelled using a double cell approach. Such an approach defines a fixed zone that spatially coincides with the porous fluid zone. This solid zone only interacts with the fluid in terms of heat transfer. The conservation equations for energy are separately solved for the fluid and solid zones. The conservation equation solved for the fluid zone is in accordance with Nield and Bejan [28]:
( 1 ε ) ( ρ c ) s T s t = ( 1 ε ) · ( λ s T s ) + α fs a fs ( T f T s )
and for the solid zone
ε ( ρ c p ) f T f t + ( ρ c p ) v · T f = ε · ( λ f T f ) + α fs a fs ( T s T f ) .
α fs is the heat transfer coefficient between the fluid and the solid and a fs is the specific surface area of the solid. Relating to the wire package of a brush seal, a fs equals to:
a fs = ( 1 ε ) d b · π · l 1 / 4 · d b 2 · π · l = ( 1 ε ) 4 d b = ( 1 ε ) · a v .
With regard to the heat conduction within the bristle package, it should be noted that, similar to the loss coefficients, it is also directional. Perpendicular to the longitudinal axis of the bristle, an effective heat conduction λ eff can be calculated, which results from a series connection of the heat conduction of the solid and fluid (see Ref. [29]):
1 λ eff = ε λ f + ( 1 ε ) λ s ,
λ eff = λ s , x = λ s , n = λ f , x = λ f , n = λ f · λ s ( 1 ε ) · λ f + ε · λ s .
Equations (11) and (12), in combination with Equation (15), lead to:
( 1 ε ) ( ρ c ) s T s t = ( 1 ε ) [ x ( k eff T s x ) + s ( k s T s s ) + n ( k eff T s n ) ] + α fs a fs ( T f T s ) ,
ε ( ρ c p ) f T f t + ( ρ c p ) v · T f = ε [ x ( k eff T f x ) + s ( k f T f s ) + n ( k eff T f n ) ] + α fs a fs ( T s T f ) .

2.2. Porosity of the Brush Seal

Referring to a cylinder packing, as in the case of a brush seal, the porosity can be expressed, as follows:
ε = 1 π · d b 2 · ρ p 4 · B · sin ( λ ) .
In Equation (18), it was neglected that the porosity is dependent on the radius. Because the circumferential area increases with increasing radius, but the number of bristles remains constant, the porosity increases in the outward radial direction. According to Qiu et al. [11], the porosity as a function of the radius can be calculated, as follows:
ε ( r ) = 1 ρ p π D R ( Δ r / s i n ( λ ) ) ( π d b 2 / 4 ) ( π ( r + Δ r ) 2 π r 2 ) B = 1 ρ p π D R d b 2 / s i n ( λ ) 4 B ( Δ r + 2 r ) .
For Δ r 0 follows from Equation (19):
ε ( r ) = 1 ρ p · π · D R · d b 2 8 · B · r · s i n ( λ ) .
In addition to the fixed geometric parameters, such as the wire diameter d b , the rotor outer diameter in relation to the package centre D R , and the laying angle λ , the porosity also depends on a size that changes depending on the operating pressure applied, namely the package width B (see Refs. [30,31]). Consequently, this size must be taken into special consideration when calibrating the model (see Section 4). Here, it is neglected that the laying angle λ also changes during rubbing, due to the radial deflection of the bristles. However, the changes are very small.

3. Description of the Model

In Figure 1, the computational doma in of the 0.5° section of the ITS brush seal test rig (see Ref. [18]).
Periodic boundary conditions are set in circumferential direction because of the rotational symmetry. The calculation doma in includes the rotor, the test seal, the axial brush seal, a part of the seal holder, and the flow areas upstream and downstream of the test seal. The test seal is modelled as a porous medium, assuming that there is no thermal equilibrium between the bristles and the fluid. The axial brush seal is also modelled as a porous medium, assuming that there is a thermal equilibrium. The calculation doma in starts at a radius of 50 mm. The recalculation of the experimental tests by means of the FE analyses (see Refs. [17,18,19]) has shown that the areas further inwards are only slightly thermally affected within the rubbing period of 30 s.
At the inlet, the total pressure and temperature are set according to the respective measured values from the experiments. At the outlets, the static pressures and temperatures are specified. At the outer surfaces ①–④ appropriate heat transfer boundary conditions are defined. Along the outer contour of the rotor, a circumferential speed corresponding to the tests is set. In the detailed section in Figure 1, it can be seen that the specific total frictional power loss q ˙ tot is applied to the contact area between the seal and rotor:
q ˙ tot ( t ) = P tot ( t ) A · ( 1 ε ) .
This quantity is based on the measurement of the total power loss during the experiments. Similar to the temperatures at the inlets and outlets, it is defined as time-dependent and specified as a boundary condition. The heat flux distribution will be self-adjusting due to the prevailing conditions. In the detailed section, it is also evident that the package width is not constant over the height. In the clamping area, the package width is defined to the minimum package width (see Section 4). The package width increases towards the inside and it reaches its maximum value at the contact point to the rotor. This specification is based on own observations and descriptions by other authors (e.g. Ref. [30]). A narrow gap was defined between the bristle package and the back plate, thus reducing the contact height of the bristles on the back plate h BP . In various publications (e.g. Refs. [32,33,34]), it has been described how the bristle package bends in axial downstream direction below the back plate when a pressure difference is applied. As a result, there is a bending in the upper area in the upstream direction and, thus, a reduction of the contact surface. This effect is to be replicated by reducing the contact surface in the model. The contact height at the back plate has an influence on the radial pressure drop and, thus, on the pressures in the pressure relief chamber, as well as on the heat transfer.
The axial position of the seal relative to the rotor is adjusted individually for all test cases, analogously to the experiment. In addition, the determined rotor elongation is taken into account to adjust the gap below the back plate as closely as possible to the gap during the rub tests.

3.1. Investigated Test Cases

The numerical test cases that are examined are based on rub tests of seal 1 under variation of the differential pressure (see Ref. [18]). The bristle package consists of bristles of Haynes 25 with a bristle diameter d B of 0.07 mm, a package density ρ p of 200 Bpmm and a seal inner diameter of 300 mm. The laying angle of the bristles λ is 45°. The rotor is made of Inconel 718 and it has a diameter of 299.5 mm on the right front edge. The cone angle of the rubbing surface is 2.95°. The pressure difference was varied in the steps of 1.0, 2.5, 4.0, 5.5 and 7.0 bar. The choice was made for these test points, because the largest influences on the heat flux distribution have been observed in the experiment when varying the pressure difference and the heat transfer coefficients at the bristles α fs = α b must change over a wide range.
The heat transfer coefficients α b are determined according to correlations for banks of staggered tubes by Gnielinski [35], Žukauskas [36], and Mart in [37]. Appendix A describes the calculation of the coefficients. Figure 2 shows the calculated values. The mean values were used for the numerical simulations. The error bars, which span an interval of ±50%, show that the correlations yield very different results. The heat transfer coefficients were assumed to be constant in the entire bristle package.
Besides the heat transfer coefficients at the bristles, the porosity ε is a very important parameter of the modelling. The porosity was calculated according to Equation (20), i.e., the radial dependence of the porosity was taken into account. Furthermore, it was taken into account that the bristle package tapers with increasing radial height. In Figure 3a,b, it can be seen that the consideration of these effects has a very large influence on the porosity.
The first figure shows the values of porosity under sole consideration of the variable package width. Because of the package width B min defined as the minimum possible width in the area of the clamping, the porosity decreases strongly towards the inside. If it is also taken into account that the available space increases with increasing package height, the decrease in porosity is much less pronounced (see Figure 3b). To further take into account that the bristle package in the area of the back plate is compressed more than the bristles upstream when a pressure difference is applied, the porosity was modified according to Equation (22):
ε ( r , x ) = ε · 1 + B ( n m ) n + 1 ε x n B m , with n = 0.5 ; m = 0.45 .
Figure 3c shows the porosity when considering the package width and the weighting according to Equation (22). In the last contour plot, all of the effects are finally combined (see Figure 3d).

3.2. Mesh Independence Study

A mesh independence study was carried out in order to assess the influence of the created mesh on the quality of results. The analysis was performed for the medium pressure difference level of 4.0 bar. The mass flow, the heat flux distribution, and the maximum temperature in the friction contact serve as comparative values.
A total of four meshes with different resolutions were compared. The k ω -SST model was used to model the turbulence. Within the bristle package, laminar calculations were performed. This assumption is valid due to the low Reynolds numbers (see Refs. [38,39,40]). This is also confirmed in this study. Within the package, the Reynolds numbers based on bristle diameter were well below 1000. They only reach this order of magnitude in the region near the back plate gap where the highest velocities are present. At first, a steady-state solution was calculated without the imposition of a heat flux boundary condition. The fluid properties were assumed to be compressible. Subsequently, a transient calculation is performed for 30 s with time-dependent specification of the total frictional power loss and temperatures from the experiment. In Table 1, the results are summarized in relation to the mesh with the highest cell number:
For the following simulations, the computational mesh with 3.6 million cells was used.

4. Calibration of the Numerical Model

The adaptation or calibration to the experimental data is inherent in the modelling of the bristle package by means of a porous medium. Because the geometry of the bristle package changes during operation as a function of the differential pressure, the leakage characteristics must first be calibrated. Once the numerical model has been calibrated, it can be validated using further measurement data (see Section 5). In the present study, the leakage mass flow and pressures in the pressure relief chamber must be adjusted. For the calibration of the leakage mass flow, the porosity of the bristle pack ε is the decisive parameter. A variation of the porosity is achieved by adjusting the package width B. Because the package width in the clamping area was set to the minimum value B min , the calibration is done by varying the package width B max (see Figure 1). The minimum package width B min is calculated from:
B min = d + 3 2 · d b ρ p · d b s i n ( λ ) 1 .
Because the bristle package is compressed depending on the pressure difference applied, the porosity also changes in the same way. Therefore, calibration must be carried out separately for each pressure stage. Figure 4a shows the results of the calibration of the leakage mass flows as discharge coefficients. The comparison of the numerically calculated values with the values from the experiments shows a very good agreement and, thus, a good adaptation to the experiments. Only at a pressure difference of 2.5 bar, the calculated flow rate is below the experimental value. The experimental value is most likely an outlier. A comparison of the two experimental measurement runs carried out confirms this suspicion. Because the actual course of the leakage mass flow for this pressure is not known and a linear course between the values at differential pressures of 1.0 and 4.0 bar is not necessarily present, the average package width from the package widths at differential pressures of 1.0 and 4.0 bar was used. In Figure 4b, the respective underlying package widths B max are plotted over the pressure difference. It becomes apparent that, above a pressure difference of about 5.0 bar, the package is not compressed any further and the maximum sealing effect is achieved. Furthermore, it is clear that, even at the highest pressure differences, the package width is still significantly larger than the theoretical minimum package width B min .
When adjusting the pressures in the pressure relief chamber, the contact height of the package at the back plate h BP plays a decisive role in addition to the porosity of the package or the package width. The contact height has a significant influence on the radial pressure drop and, thus, the chamber pressures. In Figure 5a, the calculated and measured chamber pressures are plotted over the absolute pressure upstream of the seal. After adjusting the contact heights according to Figure 5b, the chamber pressures are within the measuring range. This is spanned by the deviations of the measured values at the circumferential positions of 0 and 180°. The decrease in the contact height with an increase in the differential pressure could, be explained by an increase in the axial deflection of the bristles and, therefore, also appears to be physically plausible, as described in Section 3.

5. Validation of the Numerical Model

In this section, the model will be validated. Several comparison variables are available in order to be able to check whether the model reflects reality accurately. One of these variables is the pressure distribution in axial direction. In the absence of own measurement data, the results by Schur et al. [31] and Bayley and Long [41] are used. In particular, the results shown by Schur et al. [31] are suitable for a comparison, since one of the tested bristle packages corresponds exactly to that of the seal 1. In Figure 6a–e, the normalized pressures below the bristle packs are plotted over the normalized package width.
In general, the measured values of Schur et al. [31] are well matched. The pressure gradient increases in the area of the back plate, whereby this effect increases with increasing pressure difference. In Figure 6f, this can be clearly seen when looking at the numerically calculated pressure profiles. The differences between the measured values and the calculations can be explained by different operating conditions. In the experiments conducted by Schur et al. [31], a circumferential speed of 47 m/s was present, while the measurements by Bayley and Long [41] were carried out purely statically. In the simulation, a circumferential speed almost twice as high as the value referred to by Schur et al. [31] is set. However, Schur et al. [31] could show that, with increasing speed, the pressure level also continues to rise. Furthermore, the measured values of the investigations by Schur et al. [31] are only available up to a pressure difference of 5.0 bar.
Altogether, the flow and pressure conditions can be simulated very well with the numerical model according to the experimental tests. The assumption of a rigid bristle package or the omission of the calculation of the bristle bending does not represent a significant impairment.
The measured temperatures within the rotor structure are used for validation because the flow variables are of less interest than the distribution of heat input in the rotor and seal. In addition, a comparison of the heat flux distributions are useful. The heat flux distribution is defined as the ratio of rotor heat input and the measured total frictional power loss:
h e a t f l u x d i s t r i b u t i o n = r o t o r h e a t i n p u t t o t a l f r i c t i o n a l p o w e r l o s s .
The values of the heat flow distributions of the current numerical simulation are compared with those of the corresponding experimental tests (see Ref. [18]). For the experiments, the total frictional powers were measured and the rotor heat inputs were calculated by means of FE analyses.
In Figure 7a–e, the experimental rotor temperatures are plotted together with the rotor temperatures that are determined from the numerical simulations for all five test cases. In addition, the temperature curves of the FE analysis are shown. In each case, the points in time within the contact duration of 30 s are shown, where the highest temperatures occurred.
While the temperature differences between the experiment and the FE analysis are very small as a matter of principle, there are stronger deviations between the experiment and numerical simulation, especially at low pressure differences. At higher pressure differences above 5.5 bar, the profiles aga in agree very well. In order to check at this point whether the hypothesis by Pfefferle [13], namely the overestimation of the heat transfer coefficients at the bristles by the correlations for banks of staggered tubes, is correct, calculations with heat transfer coefficients of α b = 100   W / ( m 2 K ) and α b = 1000   W / ( m 2 K ) have been performed for the first two pressure levels. The experimental temperature level is reached in the case of a heat transfer coefficient of α b = 100   W / ( m 2 K ) and a pressure difference of 2.5 bar. The calculated temperatures are too low for all other cases.
The findings from the comparison of the rotor temperatures are also well reflected by the differences in the heat flux distributions (see Figure 7f). The trend is quite well matched and the level of the heat flux distributions corresponds approximately to that of the experiment or FE analysis. Assuming that the slope of the experiment or FE analysis profile is correct, the heat flux distributions of the numerical simulation should follow the dotted line ( Energies 14 01888 i019). For this purpose, the values from the experiment or the FE analysis ( Energies 14 01888 i011) were transposed upwards until they corresponded to the value from the numerical simulation at the operating point with a differential pressure of 5.5 bar. This operating point was selected because the best agreement between numerical simulation and experiment was found there. If the dotted line ( Energies 14 01888 i019) is used as a comparison level, it becomes clear that too little heat is introduced into the rotor at low pressure differences in the case of numerical simulation. In the case of the additional calculations with heat transfer coefficients of α b = 100   W / ( m 2 K ) or 1000   W / ( m 2 K ) , this level ( Energies 14 01888 i019) is almost reached (not drawn in the diagram). However, other factors must also play a role, because, as shown above, the rotor temperatures are still calculated significantly too low.

6. Results of the Sensitivity Study

So far, the comparison of the calculated and experimental results has shown that, at low pressure differences, the amount of heat into the rotor is underestimated by the numerical simulation. This seems to be a first indication that the heat transfer coefficients at the bristles were chosen too high. However, the fact that the temperature differences persist if the heat transfer coefficients are significantly reduced suggests that other effects or factors must also have a significant influence on the heat transfer. A possible explanation would be that the contact areas have been chosen too large and, thus, the specified specific total frictional power loss becomes too small (see Equation (21)). The package width or porosity are coupled with the leakage flow and they have been calibrated. However, a closer look at the bristle packages shows that in reality, not all of the bristles have exactly the same length. Therefore, it is possible that individual bristles do not come into contact with the rotor and consequently do not contribute to the conversion of the total frictional power loss. Photographs of the undersides of the bristle packages confirm this assumption (see Figure 8). This becomes particularly clear in the photograph shown in Figure 8a). In the area marked with the letter A, the bristles show clear signs of rubbing, whereas large parts (B) show no traces. Additionally, in the other pictures, a clear gradation of the bristles can be seen. It seems quite plausible that this effect especially occurs at low pressure differences.
A screening test was carried out for the test case at a medium pressure difference of 4.0 bar in order to investigate and quantify the influences. In addition to the influencing factors already identified, namely the heat transfer coefficients at the bristles α b and the contact surface A, the thermal conductivities of the materials of the bristles, the rotor and the seal housing were varied. The variation of the contact area was achieved by adjusting the specific total frictional power loss. Furthermore, it was examined whether the adjustment of the contact height of the bristles on the back plate h BP has an influence. The latter factor was varied in three steps from −10% over 0% to +10%, while the variation of the thermal conductivities covers the range of ±5% supplemented by an additional step at 0%. This corresponds approximately to the scattering of the available material data of the materials used. The heat transfer coefficients were varied in steps from −50% over 0% to +50%. This is roughly in accordance with the fluctuations in the predictions of the correlations under consideration. The maximum reduction of the contact area was estimated with −30% and the maximum increase in the contact area with +5%. For the central point, the contact area was −12.5%. The underlying test design complies with a definitive screening design, therefore comprising 13 tests.
The influences of the factors on the heat flux distribution and the maximum temperature in friction contact were investigated. The results of the investigation are summarized in a ma in effect diagram in Figure 9. Regarding the heat flux distribution, it can be concluded that the heat conductivities of the two direct friction partners ( λ Bristle , λ Rotor ) as well as the heat transfer coefficients at the bristles α b have a clear influence. The influences of the contact surface A and the thermal conductivity of the seal housing λ Seal play a minor role. Adjusting the contact height of the bristles on the support ring h BP does not result in a change in the heat flux distribution. This factor also has no influence on the maximum temperature in the friction contact. The contact temperature is mainly determined by the thermal conductivity of the rotor material λ Rotor and, to a greater extent, by the contact surface A.
A Pareto diagram of the standardized effects can be considered in order to be able to check which factors are significant (see Figure 10a,b). In this diagram, the absolute values of the standardized effects are shown in the order of their influence. The standardized effect is the observed effect divided by its standard deviation. The result corresponds to the t-value and is thus a measure of the significance of the effect. The selected significance level is α = 0.05 . This results in a t-value of the reference line of 2.45. If the bar of an effect extends beyond the reference line, it is significant. In the case of the heat flux distribution (see Figure 10a), this means that the heat transfer coefficients at the bristles are the decisive factor. The thermal conductivities of the rotor and bristle materials have a similarly strong effects and are also statistically significant. From the Pareto diagram for the maximum rubbing temperature in the friction contact as response variable (see Figure 10b), it is clear that the rubbing surface is the decisive factor by a large margin. The thermal conductivity of the rotor material is slightly significant.
The results of the screening test, in which the ma in influencing factors have been identified, are now to be transferred to all five test cases in order to obtain the influence of the pressure difference. The factors heat transfer coefficients at the bristles α b , thermal conductivity of the rotor material λ Rotor and the contact surface A are to be examined more closely. The thermal conductivity of the bristle material λ bristles was excluded, because the number of simulations to be performed had to be reduced. This limitation seems to be justified, since a similar effect on heat flux distribution is expected, as the thermal conductivity of the rotor material and the magnitudes of the effects are nearly the same.
The simulations were varied in accordance with a response surface design in the form of a central composite test plan. The design of a central composite design consists of a cube with two levels per factor. In addition, test points in the form of a star are also investigated. The star, starting from the central point, is created by varying the individual factors. Because the experimental points of the star exceed those of the cube, each factor is varied on five levels. With this type of experimental design, first- and second-order terms can be estimated. The range of the experimental design, in relation to the star points, corresponds to that of the screening test.
Figure 11 shows the effects of the three factors on the response variables heat flux distribution and maximum temperature in friction contact for all five test cases. As expected, the amount of heat that is introduced into the rotor increases as the thermal conductivity of the rotor material increases and decreases accordingly as the convective heat transfer at the bristles increases. The influence of the contact surface is only significant at low differential pressures. As a tendency, the changes of the varied parameters have less effect on the target value with increasing pressure difference. Thus, the effect of the pressure difference outweighs the three factors examined. This is not the case when evaluating the effects on the maximum temperature in the rubbing contact. Still apparent is the disproportionately strong influence of the contact surface on the temperature in the rubbing contact.
In Figure 12a,b, the results of the simulation are plotted as standardized effects over the pressure difference. Instead of a Pareto chart, a line chart was chosen. The reference line at t = 2.23 still represents the limit above which an effect is considered significant. Only those factors are shown in the diagrams whose effects are significant at least one pressure difference. In the case of heat flux distribution (see Figure 12a), it is clear, as outlined above, that the effects decrease significantly with increasing pressure difference. Except for the contact area A ( Energies 14 01888 i005), all applied factors rema in significant. The heat flux distribution is, however, increasingly determined by the increased pressure level and the convective heat transfer in the bristle package, which has increased in absolute numbers. Interaction effects between the individual factors are not significant. With regard to the maximum temperature in the contact area, the effects are not dependent on the pressure difference (see Figure 12b).
If the knowledge gained is transferred to the results that are discussed in Section 5, this means that a reduction of the heat transfer coefficients at the bristles, as well as a reduction of the contact area, should mainly affect the temperatures and heat flux distributions at low pressure differences. The corresponding values at high pressure differentials should show significantly smaller changes when varying the factors.
The parameter values listed in Table 2 were transferred to the five test cases and the numerical calculations were repeated in order to reduce the deviations in the temperature profiles between experiment and numerical simulation. In a first step, the thermal conductivity of the rotor was assumed to be constant. In any case, there is no change between the test cases, since the experiments always involved the same seal and the same rotor. However, if the thermal conductivity varies, the assumed values from Table 2 would probably have to be adjusted again, in order to achieve a match of the temperature profiles.
In Figure 13a–e, the experimental rotor temperatures are plotted with the temperatures determined from the numerical simulations for all five test cases over the axial measuring position. The time points within the rubbing period of 30 s at which the highest temperatures occurred are shown. It is clear that the simulated temperatures at low pressure differences now correspond to the experimental values. At high pressure differences, the numerically calculated temperatures increase as compared to the original calculations, but they are still very close to the measured temperature values.
Additionally, the curve of the heat flux distribution over the pressure difference (see Figure 13f) now follows the slope resulting from the evaluation of the experimental data by means of the FE analysis. The values of the heat flux distribution that result from the results of the numerical simulation are now always higher than those of the FE analysis. There can be many reasons for this. Finally, it has to be considered that two numerical calculations are compared with each other, which can never exactly represent reality, despite careful verification and validation of the models. In Figure 13f, error bars are shown for a range of ±10%.

7. Conclusions

The aim of the numerical analysis of the heat flux distribution was to identify the crucial factors for the heat transfer at the seal. The chosen approach was to model the bristle package as a porous medium. For the calculation of the porosity, in addition to the package width and radial height, an axial weighting was used for the first time to take that the bristle package is more strongly compressed under pressure in the vicinity of the back plate into account. Another new feature is the quantitative validation of such a model with experimental data from brush seal rub tests. Among other things, the calculations should serve to verify the applicability of the often used correlations for calculating the heat transfer coefficients at the bristles on the basis of banks of staggered tubes.
After the evaluation of the results, it has been shown that it is quite possible that the correlations for banks of staggered tubes used for brush seals significantly overestimate the heat transfer coefficients α b . Besides the reasons already mentioned by Pfefferle [13] (see Section 1), the oblique flow of the bristles also has an influence here. When applying the correlations, a straight inflow is always assumed. Publications by Moreno and Sparrow [44] and Žukauskas [36] show that the Nusselt numbers decrease with increasing inclination. According to Žukauskas [36], the Nusselt numbers can decrease by approximately 20% with an oblique flow of 45°. If the oblique flow is increased to 80°, the values decrease to just above 40%. As a result of the swirl that is caused by the rotor rotation and especially by the laying angle of the bristles, a straight inflow or flow through the bristle pack is not ensured to a large extent. Especially in the area close to the back plate, the flow follows very strongly the laying direction of the bristles (see Ref. [45]). Furthermore, the assumption of isotropic heat transfer coefficients throughout the package is a rough assumption. In addition, the correlations for calculating the heat transfer coefficients assume an average flow velocity in the package. In all literature sources applying these correlations to brush seals, this mean velocity is calculated by relating the total leakage mass flow to the gap area between the back plate and rotor. This results in average velocities that are too high. When compared to a bank of staggered tubes, the bristle package of a brush seal exhibits locally strong pressure and velocity differences. As a result, although very high convective heat transfer coefficients are conceivable on a localized basis, on average, the correlations overestimate the values.

Author Contributions

Conceptualization, M.H.; Methodology, M.H.; Formal analysis, M.H.; Investigation, M.H.; Resources, C.S. and H.-J.B.; Writing–original draft preparation, M.H.; Writing–review and editing, C.S. and H.-J.B.; Supervision, C.S. and H.-J.B.; Project administration, M.H., C.S. and H.-J.B.; Funding acquisition, M.H., C.S. and H.-J.B. All authors have read and agreed to the published version of the manuscript.

Funding

This experimental work was funded by FVV (“Forschungsvereinigung Verbrennungskraftmaschinen e.V.”), project number 1156. The APC was funded by the KIT-Publication Fund of the Karlsruhe Institute of Technology.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Acknowledgments

The experimental investigations that served as a basis for this scientific publication are the result of a research project which has been initiated by the FVV (“Forschungsvereinigung Verbrennungskraftmaschinen e.V.”) and carried out at the Institute of Thermal Turbomachinery, Karlsruhe Institute of Technology. The work has been supported financially by the FVV and the participating companies. This support is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
aNormalized transverse bristle spacing, -
ARubbing area, m 2
A Viscosity tensor, 1/m 2
a i Loss coefficient relating to viscosity, 1/m 2
a v , a fs Specific surface, -
A 1 Transversal bristle spacing, m
A 2 Diagonal bristle spacing, m
bNormalized longitudinal bristle spacing, -
BPackage width, m
B Inertia tensor, 1/m
b i Loss coefficient relating to inertia, 1/m
cNormalized diagonal bristle spacing, -
cSpecific heat capacity, kJ/(kg K)
c p Specific heat capacity at constant pressure, kJ/(kg K)
C 1 Constant, -
C 2 Constant, -
d ˜ Equivalent sphere diameter, m
d b Bristle diameter, m
D R Rotor diameter referred to the centre of the package, m
f A Empirical arrangement factor, -
F N Total contact force of the bristles in radial direction, N
F r Resistance in porous medium per volume unit, kg/(m 2 s 2 )
h fh Gap height under back plate (fence height), m
h BP Contact height of the bristles on the back plate, m
lBristle length, m
L char Characteristic length, m
mConstant, -
m ˙ BS Leakage mass flow through bristle package, kg/s
nConstant, -
N x Number of bristles in axial direction, -
N Θ Number of bristles in tangential direction, -
P fric Friction power, W
P tot Total frictional power loss, W
q ˙ tot Area-specific heat flux, W/m 2
rRadial coordinate, m
SBristle centre distance, m
tTime, s
tt-value, -
TTemperature, K
uRotor surface velocity, m/s
v , V Flow velocity, m/s
VVolume, m 3
wAverage flow velocity in the porous medium related to the back plate gap, m/s
Greek Symbols
α Level of significance, -
α b Heat transfer coefficient at bristles, W/(m 2 K)
α fs Heat transfer coefficient between fluid and solid, W/(m 2 K)
α Rotor Heat transfer coefficient at the rotor, W/(m 2 K)
β Bristle spacing, m
Δ p Pressure difference, N/m 2
ε Porosity of the bristle package, -
λ Bristle laying angle, deg
λ Thermal conductivity, W/(m K)
λ eff Effective thermal conductivity of bristle material, W/(m K)
λ Bristle Thermal conductivity of bristle material, W/(m K)
λ Seal Thermal conductivity of seal casing material, W/(m K)
λ Rotor Thermal conductivity of rotor material, W/(m K)
µ Dynamic viscosity, kg/(m s)
µ Coefficient of friction, -
ν Kinematic viscosity, m 2 /s
ρ Density, kg/m 3
ρ p Package density (bristles per millimetre of circumference, unit in text: Bpmm), 1/m
Dimensionless quantities
H g Hagen number
L q Lévêque number
N u Nusselt number
N u 1 Nusselt number of a single tube in cross-flow
N u Bundle Nusselt number of a plain tube bundle in cross-flow
N u D Nusselt number of a tube bundle in cross-flow
P r Prandtl number
R e Reynolds number of a tube bundle according to Gnielinksi [35]
R e D , max Reynolds number of a tube bundle according to Žukauskas [36]
Subscripts
DDiagonal
fFluid
LLongitudinal
l a m Laminar
m a x Maximal value
m i n Minimal value
nDirection perpendicular to the x-axis and longitudinal bristle axis
sSolid
sDirection parallel to the bristle longitudinal axis
s t a t Static quantity
TTransversal
t o t Total
t , t o t a l Total quantity
t u r b Turbulent
xDirection parallel to the x-axis of the test rig
Abbreviations
BPBack Plate
CFDComputational Fluid Dynamics
FEFinite Element
HFDHeat Flux Distribution

Appendix A. Calculation of the Heat Transfer Coefficients at the Bristles

The heat transfer coefficients α b are determined according to correlations for banks of staggered tubes by Gnielinski [35], Žukauskas [36], and Mart in [37]. Knowledge of the bristle spacing is required to apply the correlations. The characteristic quantities are shown in Figure A1. The assumptions are made that the bristles are staggered and the bristle spacing is identical at the same radial height in all directions ( A 1 = A 2 = δ ).
Figure A1. Arrangement of bristles in bristle package.
Figure A1. Arrangement of bristles in bristle package.
Energies 14 01888 g0a1
To calculate the bristle spacing δ , four equations with four unknowns are available. The bristle width B is known in relation to the radial bristle height.
1.
Number of rows in flow direction N x :
N x = B d b S L + 1
2.
Bristle spacing in flow direction S L :
S L = ( d b + δ ) 2 1 / 4 ( d b + δ ) 2 = 3 2 ( d b + δ )
3.
Bristle spacing δ :
δ = π D R N Θ d b sin ( λ ) sin ( λ )
4.
Number of bristles in circumferential direction N Θ :
N Θ = N N x

Appendix A.1. Calculation of Heat Transfer Coefficients According to Gnielinski

The average Nusselt number of a cross-flow plain tube bundle is calculated according to Gnielinski [35] from the average Nusselt number of a cross-flow single tube multiplied by an arrangement factor f A :
N u Bundle = f A N u 1 .
For a single tube, Equation (A6) applies:
N u 1 = 0.3 + N u 1 , lam 2 + N u 1 , turb 2 ,
with
N u 1 , lam = 0.664 R e Ψ , 1 P r 3
and
N u 1 , turb = 0.037 R e Ψ , 1 0.8 P r 1 + 2.443 R e Ψ , 1 0.1 ( P r ( 2 / 3 ) 1 ) .
The Reynolds number is calculated according to:
R e Ψ , 1 = w L char Ψ ν , 10 < R e Ψ , 1 < 10 6 .
The characteristic length is L char = π 2 d b . As void fraction Ψ the average porosity of the bristle pack is used. The velocity w related to the fence height is calculated as follows:
w = m ˙ BS ρ π ( D R 2 + h fh ) 2 D R 2 2 .
In P r a is the thermal conductivity of the flow medium.
P r = ν a , 0.6 < P r < 10 3
In the case of staggered arrangement, the arrangement factor f A is:
f A = 1 + 2 d b 3 S L .
The material properties of the flow medium are related to the average temperature of the inlet and outlet of the bristle package.

Appendix A.2. Calculation of Heat Transfer Coefficients According to Žukauskas

The average Nusselt number for heat transfer in a cross-flow tube bundle is calculated according to Žukauskas [36] as:
N u D = C 1 C 2 R e D , max m P r 0.36 P r P r s 1 / 4 .
The correlation is valid for: 0.7 P r 500 and 10 R e D , max 2 × 10 6 and a number of tubes in axial direction of N x 20 . All material properties of the flow medium, with the exception of P r s , are related to the average temperature of the inlet and outlet of the bristle package.
For a staggered arrangement with S T / S L < 2 and 10 3 < R e D , max < 2 × 10 5
C 1 = 0.35 S T S L 1 / 5 and m = 0.60 .
applies.
C 2 = 0.99 since N x < 20 . The Reynolds number is calculated as:
R e D , max = w max d b ν ,
with
w max = S T S T d b w .
For w, Equation (A10) applies.

Appendix A.3. Calculation of Heat Transfer Coefficients According to Martin

An alternative to the approaches by Gnielinski [35] and Žukauskas [36] is the calculation approach by Mart in [37]. The method used is based on the so-called Lévêque analogy. This links the pressure loss and the heat transfer coefficients.
The following applies to staggered tube bundles:
N u = 0.404 L q 1 / 3 ,
with the Lévêque number for b = S L / d b < 1 :
L q = 0.92 H g P r ( 4 a b / π 1 ) ( b c ) , with a = S T d b = S D d b = c .
The Hagen number is:
H g = H g lam + H g turb 1 e x p ( 1 ( R e D , max + 200 ) / 1000 ) ,
with
H g lam = 140 R e D , max ( b 0.5 0.6 ) 2 + 0.75 / a 1.6 ( 4 a b / π 1 )
and
H g turb = f t , s R e D , max 1.75 + f t , n R e D , max 2 ,
where f t , n for N x > 10 becomes zero
For f t , s
f t , s = 1.25 + 0.6 / ( a 0.85 ) 1.08 + 0.2 ( b / a 1 ) 3 0.005 ( a / b 1 ) 3 .
applies.
For R e D , max , Equation (A15) is valid.

References

  1. Dinc, O.S.; Reluzco, G.; Turnquist, N.A.; Lawen, J.; Kerber, O.; Brunner, F.; Crum, G.; Stuck, A.; Cromer, R.; Marks, P. Brush Seals in Industrial Gas Turbines—Turbine Section Interstage Sealing. In Proceedings of the 34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Cleveland, OH, USA, 13–15 July 1998; Volume AIAA98-3175. [Google Scholar] [CrossRef]
  2. Aksit, M.F.; Tichy, J.A. Wear of Brush Seals: Background and New Modeling Approach. Tribol. Trans. 1998, 41, 368–374. [Google Scholar] [CrossRef]
  3. Gorelov, G.M.; Reznik, V.E.; Tsibizov, V.I. Experimental Study of Brush Seal Flow Characteristics and Comparison with a Labyrinth Seal. Izv. Vuz. Aviatsionnaia Tekhnika 1988, 31, 43–46. [Google Scholar]
  4. Stephen, D.; Hogg, S.I. Development of Brush Seal Technology for Steam Turbine Retrofit Applications. In Proceedings of the 2003 International Joint Power Generation Conference Power, Atlanta, GA, USA, 16–19 June 2003; pp. 505–512. [Google Scholar] [CrossRef]
  5. Doğu, Y.; Aksit, M.F. Brush Seal Temperature Distribution Analysis. J. Eng. Gas Turbines Power 2006, 128, 599. [Google Scholar] [CrossRef]
  6. Demiroglu, M.; Gursoy, M.; Tichy, J.A. An Investigation of Tip Force Characteristics of Brush Seals. In Proceedings of the ASME Turbo Expo 2007, Montreal, QC, Canada, 14–17 May 2007; Volume 4, Parts A and B. pp. 1249–1260. [Google Scholar] [CrossRef]
  7. Demiroglu, M.; Tichy, J.A. An Investigation of Heat Generation Characteristics of Brush Seals. In Proceedings of the ASME Turbo Expo 2007, Montreal, QC, Canada, 14–17 May 2007; Volume 4, Parts A and B. pp. 1261–1270. [Google Scholar] [CrossRef]
  8. Ruggiero, E.J.; Allen, J.; Demiroglu, M.; Lusted, R.M. Heat Generation Characteristics of a Kevlar Fiber Brush Seal. In Proceedings of the 43rd AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Cincinnati, OH, USA, 8–11 July 2007; Volume AIAA2007-5738. [Google Scholar] [CrossRef]
  9. Ruggiero, E.J.; Allen, J.; Lusted, R.M. Heat Generation Characteristics of a Carbon Fiber Brush Seal. In Proceedings of the 44th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Hartford, CT, USA, 21–23 July 2008; Volume AIAA2008-4508. [Google Scholar] [CrossRef]
  10. Qiu, B.; Li, J. Numerical Investigations on the Heat Transfer Behavior of Brush Seals Using Combined Computational Fluid Dynamics and Finite Element Method. J. Heat Transf. 2013, 135, 122601. [Google Scholar] [CrossRef]
  11. Qiu, B.; Li, J.; Feng, Z. Investigation of Conjugate Heat Transfer in Brush Seals Using Porous Media Approach Under Local Thermal Non-Equilibrium Conditions. In Proceedings of the ASME Turbo Expo: Turbine Technical Conference and Exposition—2014, Montreal, QC, Canada, 15–19 June 2015; Volume 5C: Heat Transfer, p. V05CT15A009. [Google Scholar] [CrossRef]
  12. Pfefferle, D.; Dullenkopf, K.; Bauer, H.J. Design and Validation of a New Test Rig for Brush Seal Testing Under Engine Relevant Conditions. In Proceedings of the ASME Turbo Expo 2011, Vancouver, BC, Canada, 6–10 June 2011; Volume 5: Heat Transfer, Parts A and B. pp. 679–689. [Google Scholar] [CrossRef]
  13. Pfefferle, D. Untersuchung des Wärmeeintrags und der Reibleistung Anstreifender Bürstendichtungen: Dissertation; Forschungsberichte aus dem Institut für Thermsiche Strömungsmaschinen, Logos Verlag: Berlin, Germany, 2017; Volume 65. [Google Scholar]
  14. Braun, M.J.; Hendricks, R.C.; Canacci, V. Flow Visualization in a Simulated Brush Seal. In Proceedings of the ASME 1990 International Gas Turbine and Aeroengine Congress and Exposition, Brussels, Belgium, 11–14 June 1990; p. V005T16A008. [Google Scholar] [CrossRef] [Green Version]
  15. Owen, A.K.; Jones, T.V.; Guo, S.M.; Hogg, S.I. An Experimental and Theoretical Study of Brush Seal and Shaft Thermal Interaction. In Proceedings of the ASME Turbo Expo 2003, Atlanta, GA, USA, 16–19 June 2003; Volume 5: Turbo Expo 2003, Parts A and B. pp. 977–984. [Google Scholar] [CrossRef]
  16. Doğu, Y. Investigation of Brush Seal Flow Characteristics Using Bulk Porous Medium Approach. In Proceedings of the ASME Turbo Expo 2003, Atlanta, GA, USA, 16–19 June 2003; Volume 5: Turbo Expo 2003, Parts A and B. pp. 1091–1101. [Google Scholar] [CrossRef]
  17. Hildebrandt, M.R.; Schwarz, H.; Schwitzke, C.; Bauer, H.J.; Friedrichs, J. Effects of the Back Plate Inner Diameter on the Frictional Heat Input and General Performance of Brush Seals. Aerospace 2018, 5, 58. [Google Scholar] [CrossRef] [Green Version]
  18. Hildebrandt, M.; Schwitzke, C.; Bauer, H.J. Experimental Investigation on the Influence of Geometrical Parameters on the Frictional Heat Input and Leakage Performance of Brush Seals. J. Eng. Gas Turbines Power 2019, 141, 042504. [Google Scholar] [CrossRef]
  19. Hildebrandt, M.; Schwitzke, C.; Bauer, H.J. Contamination of Brush Seals by Oil and Salt and Its Impact on Rubbing and Hysteresis Behaviour. Int. J. Turbomach. Propuls. Power 2019, 4, 40. [Google Scholar] [CrossRef] [Green Version]
  20. Hildebrandt, M.; Munz, O.; Schwitzke, C.; Bauer, H.J. Experimental Study on the Waer- in Behaviour of Brush Seals. In Proceedings of the ASME Turbo Expo: Turbine Technical Conference and Exposition—2020, Online, 21–25 September 2020; p. 999. [Google Scholar]
  21. Hildebrandt, M. E in Beitrag zum Verständnis des Anstreif- und Leckageverhaltens von Bürstendichtungen bei Variation geometrischer Parameter, Dissertation; Forschungsberichte aus dem Institut für Thermsiche Strömungsmaschinen, Logos Verlag: Berlin, Germany, 2021; Volume 82. [Google Scholar]
  22. Darcy, H. Les Fontaines Publiques de la Ville de Dijon: Exposition et Application des Principes a Suivre et des Formules a Employer dans les Questions de Distribution D’Eau; Victor Dalmont: Paris, France, 1856. [Google Scholar]
  23. Forchheimer, P.H. Wasserbewegung durch Boden. Z. Acker Pflanzenbau 1901, 49, 1736–1749. [Google Scholar]
  24. Ergun, S. Fluid Flow Through Packed Columns. Chem. Eng. Prog. 1952, 48, 89–94. [Google Scholar]
  25. Chew, J.W.; Hogg, S.I. Porosity Modeling of Brush Seals. J. Tribol. 1997, 119, 769. [Google Scholar] [CrossRef]
  26. Pröstler, S. Modellierung und Numerische Berechnung von Wellenabdichtungen in Bürstenbauart: Zugl. Dissertation an der Ruhr-Universität Bochum 2005; Dr. Hut: München, Germany, 2005. [Google Scholar]
  27. Chew, J.W.; Lapworth, B.L.; Millener, P.J. Mathematical Modeling of Brush Seals. Int. J. Heat Fluid Flow 1995, 16, 493–500. [Google Scholar] [CrossRef] [Green Version]
  28. Nield, D.A.; Bejan, A. Convection in Porous Media, 5th ed.; Springer: Cham, Switzerland, 2017. [Google Scholar]
  29. Kaviany, M. Principles of Heat Transfer in Porous Media, 2nd ed.; Mechanical Engineering Series; Springer: New York, NY, USA, 1995. [Google Scholar]
  30. Chen, L.H.; Wood, P.E.; Jones, T.V.; Chew, J.W. Detailed Experimental Studies of Flow in Large Scale Brush Seal Model and a Comparison With CFD Predictions. J. Eng. Gas Turbines Power 2000, 122, 672. [Google Scholar] [CrossRef]
  31. Schur, F.; Friedrichs, J.; Flegler, J.; Georgakis, C.; Polklas, T. Pressure Distributions Below Brush Seals at Varying Operating Conditions. In Proceedings of the ASME Turbo Expo: Turbine Technical Conference and Exposition—2018, Oslo, Norway, 11–15 June 2018; Volume GT2018-76194. [Google Scholar] [CrossRef]
  32. Guardino, C.; Chew, J.W. Numerical Simulation of 3D Bristle Bending in Brush Seals. In Proceedings of the ASME Turbo Expo 2004, Vienna, Austria, 14–17 June 2004; Volume 4: Turbo Expo 2004, pp. 277–288. [Google Scholar] [CrossRef]
  33. Aksit, M.F. Analysis of Brush Seal Bristle Stresses With Pressure-Friction Coupling. In Proceedings of the ASME Turbo Expo 2003, Atlanta, GA, USA, 16–19 June 2003; Volume 5: Turbo Expo 2003, Parts A and B. pp. 1061–1069. [Google Scholar] [CrossRef]
  34. Aksit, M.F.; Doğu, Y.; Kandemir, I.; Kizil, H. An Investigation of Pressure Stiffness Coupling in Brush Seals. In Proceedings of the Eleventh International Symposium on Transport Phenomena and Dynamics of Rotating Machinery 2006, Honolulu, HI, USA, 26 February–2 March 2006; Volume ISROMAC-11-4. [Google Scholar]
  35. Gnielinski, V. Querumströmte Einzelne Rohreihen und Rohrbündel. In VDI-Wärmeatlas; Springer Reference, Springer Vieweg: Berlin, Germany, 2013; pp. 819–823. [Google Scholar]
  36. Žukauskas, A. Heat Transfer from Tubes in Crossflow. In Advances in Heat Transfer; Hartnett, J.P., Irvine, T., Eds.; Academic Press: New York, NY, USA, 1972; Volume 8, pp. 93–160. [Google Scholar] [CrossRef]
  37. Martin, H. The Generalized Lévêque Equation and its Practical use for the Prediction of Heat and Mass Transfer Rates from Pressure Drop. Chem. Eng. Sci. 2002, 57, 3217–3223. [Google Scholar] [CrossRef]
  38. Pekris, M.J.; Franceschini, G.; Gillespie, D. Effect of Geometric Changes in an Idealised Contacting Brush Seal Bristle Pack on Typical Key Performance Measures. In Proceedings of the ASME Turbo Expo 2011, Vancouver, BC, Canada, 6–10 June 2011; Volume 5:Heat Transfer, Parts A and B. pp. 999–1010. [Google Scholar] [CrossRef]
  39. Kudriavtsev, V.V.; Braun, M.J. Model Developments for the Brush Seal Numerical Simulation. J. Propuls. Power 1996, 12, 193–201. [Google Scholar] [CrossRef]
  40. Hendricks, R.C.; Kudriavtsev, V.V.; Braun, M.J.; Athavale, M.M. Flows in Pinned Arrays Simulating Brush Seals. In Proceedings of the Sixth International Congress on Fluid Dynamics and Propulsion, Cairo, Egypt, 29–31 December 1996; pp. 258–266. [Google Scholar]
  41. Bayley, F.J.; Long, C.A. A Combined Experimental and Theoretical Study of Flow and Pressure Distributions in a Brush Seal. J. Eng. Gas Turbines Power 1993, 115, 404. [Google Scholar] [CrossRef]
  42. Proctor, M.P.; Walker, J.F.; Perkins, H.D.; Hoopes, J.F.; Williamson, G.S. Brush Seals for Cryogenic Applications: Performance, Stage Effects, and Preliminary Wear Results in LN2 and LH2: NASA Technical Paper 3536; NASA: Washington, DC, USA, 1996. [Google Scholar]
  43. Büscher, S. Untersuchung von mehrstufigen Bürstendichtugnen für Dampfturbinen: Dissertation; Mitteilungen des Pfleiderer-Instituts für Strömungsmaschinen; Verlag und Bildarchiv W.H. Faragallah: Sulzbach, Germany, 2010; Volume 15. [Google Scholar]
  44. Moreno, A.; Sparrow, E.M. Heat Transfer, Pressure Drop, and Fluid Flow Patterns in Yawed Tube Banks. Int. J. Heat Mass Transf. 1987, 30, 1979–1995. [Google Scholar] [CrossRef]
  45. Franceschini, G.; Jones, T.V.; Gillespie, D. Improved Understanding of Blow-Down in Filament Seals. J. Turbomach. 2010, 132, 41004. [Google Scholar] [CrossRef]
Figure 1. Computational doma in of the numerical model; Fluid ( Energies 14 01888 i001), Solid ( Energies 14 01888 i002), Porous media ( Energies 14 01888 i003).
Figure 1. Computational doma in of the numerical model; Fluid ( Energies 14 01888 i001), Solid ( Energies 14 01888 i002), Porous media ( Energies 14 01888 i003).
Energies 14 01888 g001
Figure 2. Heat transfer coefficients at the bristles α b ; Gnielinski [35] ( Energies 14 01888 i004), Žukauskas [36] ( Energies 14 01888 i005), Mart in [37] ( Energies 14 01888 i006), Average values ( Energies 14 01888 i007).
Figure 2. Heat transfer coefficients at the bristles α b ; Gnielinski [35] ( Energies 14 01888 i004), Žukauskas [36] ( Energies 14 01888 i005), Mart in [37] ( Energies 14 01888 i006), Average values ( Energies 14 01888 i007).
Energies 14 01888 g002
Figure 3. Distribution of porosity ε in the bristle package with consideration of: (a): width, (b): radius, (c): width + ax. compression, and (d): width + radius + ax. compression
Figure 3. Distribution of porosity ε in the bristle package with consideration of: (a): width, (b): radius, (c): width + ax. compression, and (d): width + radius + ax. compression
Energies 14 01888 g003
Figure 4. (a): Discharge coefficients from experiment ( Energies 14 01888 i008) and numerical simulation (CFD) ( Energies 14 01888 i009) under variation of the pressure difference; (b): Package width depending on the pressure difference ( Energies 14 01888 i009), B min ( Energies 14 01888 i010).
Figure 4. (a): Discharge coefficients from experiment ( Energies 14 01888 i008) and numerical simulation (CFD) ( Energies 14 01888 i009) under variation of the pressure difference; (b): Package width depending on the pressure difference ( Energies 14 01888 i009), B min ( Energies 14 01888 i010).
Energies 14 01888 g004
Figure 5. (a): Pressures in the pressure relief chamber as a function of the inlet pressure, experiment ( Energies 14 01888 i011), CFD ( Energies 14 01888 i009); (b): Contact height at the back plate h BP depending on the pressure difference.
Figure 5. (a): Pressures in the pressure relief chamber as a function of the inlet pressure, experiment ( Energies 14 01888 i011), CFD ( Energies 14 01888 i009); (b): Contact height at the back plate h BP depending on the pressure difference.
Energies 14 01888 g005
Figure 6. Norm. pressure under the package plotted against the norm. package width, (a): CFD ( Energies 14 01888 i012), Schur et al. [31] ( Energies 14 01888 i011), both Δ p = 1.0 bar; (b): CFD Δ p = 2.5 bar ( Energies 14 01888 i012), Schur (2018) Δ p = 2.0 bar ( Energies 14 01888 i011), Schur et al. [31] Δ p = 3.0 bar ( Energies 14 01888 i013), Bayley and Long [41] R p < 3.8 bar ( Energies 14 01888 i014); (c): CFD Δ p = 4.0 bar ( Energies 14 01888 i012), Schur (2018) Δ p = 3.0 bar ( Energies 14 01888 i011), Schur et al. [31] Δ p = 5.0 bar ( Energies 14 01888 i013); (d): CFD Δ p = 5.5 bar ( Energies 14 01888 i012), Schur et al. [31] Δ p = 5.0 bar ( Energies 14 01888 i011); (e): CFD Δ p = 7.0 bar ( Energies 14 01888 i012), Schur et al. [31] Δ p = 5.0 bar ( Energies 14 01888 i011); (f): CFD Δ p = 1.0 bar ( Energies 14 01888 i015) – Δ p = 7.0 bar ( Energies 14 01888 i016).
Figure 6. Norm. pressure under the package plotted against the norm. package width, (a): CFD ( Energies 14 01888 i012), Schur et al. [31] ( Energies 14 01888 i011), both Δ p = 1.0 bar; (b): CFD Δ p = 2.5 bar ( Energies 14 01888 i012), Schur (2018) Δ p = 2.0 bar ( Energies 14 01888 i011), Schur et al. [31] Δ p = 3.0 bar ( Energies 14 01888 i013), Bayley and Long [41] R p < 3.8 bar ( Energies 14 01888 i014); (c): CFD Δ p = 4.0 bar ( Energies 14 01888 i012), Schur (2018) Δ p = 3.0 bar ( Energies 14 01888 i011), Schur et al. [31] Δ p = 5.0 bar ( Energies 14 01888 i013); (d): CFD Δ p = 5.5 bar ( Energies 14 01888 i012), Schur et al. [31] Δ p = 5.0 bar ( Energies 14 01888 i011); (e): CFD Δ p = 7.0 bar ( Energies 14 01888 i012), Schur et al. [31] Δ p = 5.0 bar ( Energies 14 01888 i011); (f): CFD Δ p = 1.0 bar ( Energies 14 01888 i015) – Δ p = 7.0 bar ( Energies 14 01888 i016).
Energies 14 01888 g006
Figure 7. (ae): Measuring points from experiment ( Energies 14 01888 i017), Fit of the measured values from experiment ( Energies 14 01888 i016), FE-Analysis ( Energies 14 01888 i010), CFD with α b = adjusted ( Energies 14 01888 i018), CFD with α b = 100   W / ( m 2 K ) ( Energies 14 01888 i019), CFD with α b = 1000   W / ( m 2 K ) ( Energies 14 01888 i020); (f): Heat flux distribution under variation of the pressure difference, Experiment ( Energies 14 01888 i011), CFD ( Energies 14 01888 i012), CFD (target) ( Energies 14 01888 i019), (a) Δ p = 1.0 bar, (b) Δ p = 2.5 bar, (c) Δ p = 4.0 bar, (d) Δ p = 5.5 bar, (e) Δ p = 7.0 bar.
Figure 7. (ae): Measuring points from experiment ( Energies 14 01888 i017), Fit of the measured values from experiment ( Energies 14 01888 i016), FE-Analysis ( Energies 14 01888 i010), CFD with α b = adjusted ( Energies 14 01888 i018), CFD with α b = 100   W / ( m 2 K ) ( Energies 14 01888 i019), CFD with α b = 1000   W / ( m 2 K ) ( Energies 14 01888 i020); (f): Heat flux distribution under variation of the pressure difference, Experiment ( Energies 14 01888 i011), CFD ( Energies 14 01888 i012), CFD (target) ( Energies 14 01888 i019), (a) Δ p = 1.0 bar, (b) Δ p = 2.5 bar, (c) Δ p = 4.0 bar, (d) Δ p = 5.5 bar, (e) Δ p = 7.0 bar.
Energies 14 01888 g007
Figure 8. Zones with different contacts on the running surface; (a): Proctor et al. [42], (b): @ ITS, (c): Buscher [43].
Figure 8. Zones with different contacts on the running surface; (a): Proctor et al. [42], (b): @ ITS, (c): Buscher [43].
Energies 14 01888 g008
Figure 9. Main effect diagrams for heat flux distribution ( Energies 14 01888 i016) and maximum temperature in rubbing contact ( Energies 14 01888 i010).
Figure 9. Main effect diagrams for heat flux distribution ( Energies 14 01888 i016) and maximum temperature in rubbing contact ( Energies 14 01888 i010).
Energies 14 01888 g009
Figure 10. Pareto diagrams of standardized effects, significance level α = 0 , 05 , A: h BP , B: λ Bristle , C: λ Rotor , D: λ Seal , E: α b , F: A, (a): Response variable: Heat flux distribution, (b): Response variable: Max. temperature in rubbing contact.
Figure 10. Pareto diagrams of standardized effects, significance level α = 0 , 05 , A: h BP , B: λ Bristle , C: λ Rotor , D: λ Seal , E: α b , F: A, (a): Response variable: Heat flux distribution, (b): Response variable: Max. temperature in rubbing contact.
Energies 14 01888 g010
Figure 11. Main effect diagrams for heat flux distribution (top) and maximum temperature in rubbing contact (bottom); Δ p = 1.0 bar ( Energies 14 01888 i015) − Δ p = 7.0 bar ( Energies 14 01888 i016).
Figure 11. Main effect diagrams for heat flux distribution (top) and maximum temperature in rubbing contact (bottom); Δ p = 1.0 bar ( Energies 14 01888 i015) − Δ p = 7.0 bar ( Energies 14 01888 i016).
Energies 14 01888 g011
Figure 12. Standardized effects over variation of the pressure difference, significance level α = 0.05 , λ Rotor ( Energies 14 01888 i009), α b ( Energies 14 01888 i004), A ( Energies 14 01888 i005), λ Rotor × λ Rotor ( Energies 14 01888 i021), α b × α b ( Energies 14 01888 i006), A × A ( Energies 14 01888 i010), λ Rotor × A ( Energies 14 01888 i022), Reference line ( Energies 14 01888 i015), (a): Response variable: Heat flux distribution, (b): Response variable: Max. temperature in rubbing contact.
Figure 12. Standardized effects over variation of the pressure difference, significance level α = 0.05 , λ Rotor ( Energies 14 01888 i009), α b ( Energies 14 01888 i004), A ( Energies 14 01888 i005), λ Rotor × λ Rotor ( Energies 14 01888 i021), α b × α b ( Energies 14 01888 i006), A × A ( Energies 14 01888 i010), λ Rotor × A ( Energies 14 01888 i022), Reference line ( Energies 14 01888 i015), (a): Response variable: Heat flux distribution, (b): Response variable: Max. temperature in rubbing contact.
Energies 14 01888 g012
Figure 13. (ae): Measuring points from experiment ( Energies 14 01888 i017), Fit of the measured values from experiment ( Energies 14 01888 i016), FE-Analysis ( Energies 14 01888 i010), CFD ( Energies 14 01888 i018); (f): Heat flux distribution under variation of the pressure difference, Experiment ( Energies 14 01888 i011), CFD ( Energies 14 01888 i012), (a) Δ p = 1.0 bar, (b) Δ p = 2.5 bar, (c) Δ p = 4.0 bar, (d) Δ p = 5.5 bar, (e) Δ p = 7.0 bar.
Figure 13. (ae): Measuring points from experiment ( Energies 14 01888 i017), Fit of the measured values from experiment ( Energies 14 01888 i016), FE-Analysis ( Energies 14 01888 i010), CFD ( Energies 14 01888 i018); (f): Heat flux distribution under variation of the pressure difference, Experiment ( Energies 14 01888 i011), CFD ( Energies 14 01888 i012), (a) Δ p = 1.0 bar, (b) Δ p = 2.5 bar, (c) Δ p = 4.0 bar, (d) Δ p = 5.5 bar, (e) Δ p = 7.0 bar.
Energies 14 01888 g013
Table 1. The results of mesh independence study.
Table 1. The results of mesh independence study.
# CellsMass Flow RateHFDMax. Temperature at Rubbing Contact
890,000 4.43 % 1.96 % + 5.7 K
1,760,000 3.18 % 0.84 % + 4.9 K
3,600,000 0.94 % + 1.26 % 0.4 K
7,960,000---
Table 2. Parameters of the numerical simulation
Table 2. Parameters of the numerical simulation
Test CaseHeat Transfer CoefficientRubbing AreaThermal Conductivity
Δ p α b A λ Rotor
1.0 bar −96.5% 102 W / ( m 2 K ) −17.5%±0%
2.5 bar −85% 585 W / ( m 2 / K ) −10%±0%
4.0 bar −77.5% 1148 W / ( m 2 / K ) −5%±0%
5.5 bar −65% 2170 W / ( m 2 / K ) ±0%±0%
7.0 bar −50% 3650 W / ( m 2 / K ) ±0%±0%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hildebrandt, M.; Schwitzke, C.; Bauer, H.-J. Analysis of Heat Flux Distribution during Brush Seal Rubbing Using CFD with Porous Media Approach. Energies 2021, 14, 1888. https://doi.org/10.3390/en14071888

AMA Style

Hildebrandt M, Schwitzke C, Bauer H-J. Analysis of Heat Flux Distribution during Brush Seal Rubbing Using CFD with Porous Media Approach. Energies. 2021; 14(7):1888. https://doi.org/10.3390/en14071888

Chicago/Turabian Style

Hildebrandt, Manuel, Corina Schwitzke, and Hans-Jörg Bauer. 2021. "Analysis of Heat Flux Distribution during Brush Seal Rubbing Using CFD with Porous Media Approach" Energies 14, no. 7: 1888. https://doi.org/10.3390/en14071888

APA Style

Hildebrandt, M., Schwitzke, C., & Bauer, H. -J. (2021). Analysis of Heat Flux Distribution during Brush Seal Rubbing Using CFD with Porous Media Approach. Energies, 14(7), 1888. https://doi.org/10.3390/en14071888

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop