Next Article in Journal
Simulation of Optimal Driving for Minimization of Fuel Consumption or NOx Emissions in a Diesel Vehicle
Next Article in Special Issue
Dependence of Conjugate Heat Transfer in Ribbed Channel on Thermal Conductivity of Channel Wall: An LES Study
Previous Article in Journal
Perception of the Quality of Smart City Solutions as a Sense of Residents’ Safety
Previous Article in Special Issue
Experimental and Numerical Study of Heat Pipe Heat Exchanger with Individually Finned Heat Pipes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heat Transfer Analysis for Non-Contacting Mechanical Face Seals Using the Variable-Order Derivative Approach

Department of Manufacturing Engineering and Metrology, Faculty of Mechatronics and Mechanical Engineering, Kielce University of Technology, Aleja Tysiaclecia Panstwa Polskiego 7, 25-314 Kielce, Poland
Energies 2021, 14(17), 5512; https://doi.org/10.3390/en14175512
Submission received: 20 July 2021 / Revised: 24 August 2021 / Accepted: 31 August 2021 / Published: 3 September 2021
(This article belongs to the Special Issue Numerical Heat Transfer and Fluid Flow 2021)

Abstract

:
This article presents a variable-order derivative (VOD) time fractional model for describing heat transfer in the rotor or stator in non-contacting mechanical face seals. Most theoretical studies so far have been based on the classical equation of heat transfer. Recently, constant-order derivative (COD) time fractional models have also been used. The VOD time fractional model considered here is able to provide adequate information on the heat transfer phenomena occurring in non-contacting face seals, especially during the startup. The model was solved analytically, but the characteristic features of the model were determined through numerical simulations. The equation of heat transfer in this model was analyzed as a function of time. The phenomena observed in the seal include the conduction of heat from the fluid film in the gap to the rotor and the stator, followed by convection to the fluid surrounding them. In the calculations, it is assumed that the working medium is water. The major objective of the study was to compare the results of the classical equation of heat transfer with the results of the equations involving the use of the fractional-order derivative. The order of the derivative was assumed to be a function of time. The mathematical analysis based on the fractional differential equation is suitable to develop more detailed mathematical models describing physical phenomena.

1. Introduction

Over the last few years, there have been many studies using fractional-order differential and integral operators to generalize classical differential and integral calculus with the aim of further understanding the nature of complex systems. Currently, attempts are being made to apply fractional calculus to solve various physical, mechanical, biological, or chemical problems. While classical integer-order operators are dependent only on the local behavior of the function, fractional-order operators accumulate all the information about the function. Another fundamental feature of fractional derivatives is that they are defined along a segment, not at a point, as is the case with classical derivatives. This feature ensures a more accurate and effective analysis of different phenomena. One of the shortcomings that most models have to overcome is that, mathematically, velocity is an instantaneous velocity defined at a point. Thus, it can be seen from the literature that differential calculus is used in the heat transfer theory.
It has been found, for example, that variable thermal conductivity is a key physical property of materials, especially when it is dependent on temperature. Variable thermal conductivity is of significance in a wide range of applications, including modern physics and mechanical engineering. It is taken into consideration primarily to determine the effect of temperature on the performance of machine elements. Special attention is paid to this property when rapid changes in temperature occur, as they may affect the operation of mechanical assemblies and subassemblies, and consequently the whole machine. The relationship between variable thermal conductivity and fractional differential calculus has been investigated extensively for problems related to fluid mechanics [1,2], viscous elasticity [3], relaxation [4], thermal elasticity and conductivity [5,6], control theory [7,8], and many others [9,10,11]. Over the years, researchers have extended the classical theory of elasticity and introduced the theory of thermoelasticity.
The mechanical properties, especially the relaxation, of additively manufactured materials were analyzed, for example, in [4]. A mathematical model was developed using fractional differential calculus to further fit the relaxation curves with the experimental data.
Warbhe et al. [5] used the quasi-static approach to the fractional-order theory of thermoelasticity to solve a two-dimensional problem for a thin circular plate with zero temperature across the lower surface and constant and evenly distributed temperature across the insulated upper surface. The integral transform technique was employed to solve the physical problem, while the displacement potential function was applied to calculate thermal stresses. Povstenko and Kyrylych [6] studied the interface between two solids. They considered generalized boundary conditions of a non-ideal thermal contact to solve the heat conduction equation of the fractional order as a function of time using the Caputo derivative.
The extensive monograph by Kaczorek [7] presented the theoretical and practical aspects of the application of fractional calculus to fractional discrete-time linear systems.
The research by Nowacki, described for instance in [12,13,14], was a breakthrough in this field. His findings are still valid and thus widely cited, as they can be applied to solve a large number of theoretical and practical engineering problems.
Knowledge of all theories available in this area is crucial to understand the deformations of elastic materials. Their thermal and mechanical behavior needs to be taken into account whenever thermoelastic deformations result in leaks, as is the case with non-contacting face seals, where the gap height may range from one to several micrometers.
Fractional differential calculus has proved well-suited to predict many physical phenomena associated with elastic media, e.g., thermal conductivity, heat transfer, and viscoelasticity. The classical differential equations describing physical phenomena are modified using the variable-order derivative (VOD) time fractional approach. Recently, there has been much research focusing on extending the applications of fractional differential calculus, e.g., [15,16,17,18]. These studies provide an insight into many important aspects associated with heat conduction. In [15], for example, Povstenko presented fundamental solutions to the time-fractional advection diffusion equation for two cases: the plane and the half-plane. He also considered Cauchy problems, inverse source problems, and Dirichlet problems. The solutions were expressed in terms of Bessel integral functions combined with Mittag–Leffler functions. In his other articles [16,17], Povstenko discussed solutions for heat transfer in composite media, in which he employed the time-fractional heat conduction equation with the Caputo derivative of fractional order to describe heat transfer in both constituent materials (0 < α ≤ 2 and 0 < β ≤ 2, respectively). The problem was solved under ideal contact conditions. This suggests that the temperature of the two materials and the heat fluxes were the same.
Liu et al. [19] proposed an approach to solving the time fractional nonlinear heat conduction equation. Similar solutions were proposed by Koca and Lotfy [20,21].
Another issue is an inverse problem in heat transfer. Liu and Feng [22], for instance, solved it using the fractional differential equation for a time-dependent derivative. The solution to the 2D time-fractional inverse problem of diffusion was based on an improved kernel technique. Maximum a posteriori estimation was required. Convergence was achieved using the regularization term and deriving a priori probability.
The primary aim of this article was to show how the results obtained by solving the classical equation of heat transfer differ from those obtained for the same physical phenomenon by employing the fractional differential equation. The study involved using the derivative order as a time-dependent function. The results from the two approaches were compared graphically. The analytical solutions of the main physical quantities were developed using the Marchi–Zgrablich transform, the finite Fourier cosine transform, and, finally, the Laplace transform.

2. Mathematical Model

The classical Fourier heat equation can be written as:
q = λ ϑ ,
where q—heat flux vector, ϑ = T T o and λ —thermal conductivity.
When combined with the energy conservation principle, the heat equation can be used to calculate the local deformation:
q t = ρ C p ϑ t t ,
with ρ being the density and Cp the specific heat capacity.
Many researchers have examined the application of differential calculus or integral calculus involving the use of derivatives of fractional orders. Fractional calculus is a natural extension of the notions of differentials and integrals used in classical differential calculus and integral calculus, respectively. Recently, there has been much interest in the use of fractional differential calculus to explain and model many physical phenomena. For example, fractional differential calculus has been used to design and/or model PIλDµ controllers, heat conduction [23,24,25], thermoelasticity [26], complex nonlinear systems [7], supercapacitors, electrical and mechanical systems [27,28,29], electrical filters, dielectric relaxation, diffusion, and viscoelasticity [30].
The time fractional heat conduction equation for the rotor can be written as [31]:
α ϑ t α = κ Δ ϑ 0 < α 2 a n d κ = K ρ C p .  
For the case considered here, it was assumed that α = α t . The function α t was defined as a function of time (see Figure 1).
The function describing the change in the coefficient α t can be written using the following formula:
α t = 1 + 1 1 + 10 3 t t 0 t 0 , 2 1 t > 2 ,
The transition to the time fractional heat conduction equation with the boundary conditions, was described, for example, in [32] as:
α ϑ t α = κ 1 r ϑ r + 2 ϑ r 2 + 2 ϑ z 2 ,   for   r i r r o ; 0 z L ; t > 0 ,
where L is the thickness of the sealing rings.
The Caputo derivative of the fractional order can be defined as described in [15]:
α f t t α = 1 Γ n α 0 t t τ n α 1 d n f τ d τ n d τ , n 1 < α < n , d n f t d t n , α = n , ,
where Γ α is the gamma function, with the initial conditions being: t = 0 ϑ = 0 , 0 < α 2 , t = 0 ϑ t = 0 , 1 < α 2 .
The heat conduction equation for the stator is thus rewritten as:
1 r θ r + 2 θ r 2 + 2 θ z 2 = 0 ,   for   r i r r o ; 0 z L ,
The heat transfer model analyzed here also needs to take into consideration the heat flux generated in the gap, which can be described with a relationship based on the simplified energy equation. A similar method was used in [32]:
μ ν ϕ z 2 + λ f 2 T f z 2 = 0
The velocity of the fluid particles ν φ is linearly variable. Ranging from zero on the stator surface to the value of ω r on the rotor surface, it can be described as:
ν ϕ z = ω r h
The energy Equation (8) was solved taking into account relationship Equation (9) to calculate the distribution of temperature in the gap.
Considering sealing rings with unmodified face surfaces and neglecting other factors affecting their geometry, e.g., mechanical deformations, we can assume that the height of the gap is constant: h = c o n s t . The gap height considered by other researchers is that discussed, for example, in reference [33].
As dynamic viscosity is a temperature-dependent parameter, here it can be defined using the formula proposed by Li et al. [1]:
μ = μ o exp b T m T o .
The average temperature of the medium in the gap can be determined from:
T m = 1 h r 0 h r T f d z .
Equation (10) is a function describing the distribution of dynamic viscosity μ r in the radial direction.
The seal performance is largely affected by the forces acting on the stator and the rotor. The closing force produced by the spring ensures the leak tightness of the device in the OFF mode and the contact of the stator and the rotor. However, when the device operates (is ON), the opening force is generated by the pressure of the medium in the gap. The distribution of pressure can be described using the one-dimensional Reynolds equation:
d d r r ρ h 3 μ d p d r = 0 ,
where the boundary conditions are:
p ( r ) r = r i = p i ,   p ( r ) r = r o = p o
Equation (12) is the simplified Reynolds equation describing the changes in pressure of the incompressible medium in the gap only in the radial direction.

3. Boundary Conditions

The schematic diagram in Figure 2 illustrates the non-contacting face seal with the flexibly mounted stator (1) and the rotor (2) connected to the shaft (6) of the turbo machine. The gap between the stator and the rotor, with a height ranging from one micrometer to several or even more than ten micrometers, is filled with a sealing medium, e.g., water. The distribution of temperature on the stator and rotor surfaces is dependent mainly on the heat flux in the gap.
The boundary conditions are represented graphically in Figure 3.
In this analysis, it is assumed that the inner and bottom surfaces of the rotor r i and the inner and top surfaces of the stator are not in contact with the surroundings; thus, heat transfer for these surfaces takes the general form: ϑ n = θ n = 0 . In a real system, these surfaces are in contact with other elements of the seal characterized by different physical properties. Under certain conditions, heat transfer taking place between these elements can be assumed to be negligible. Another reason to introduce the boundary conditions ( ϑ n = 0 ) for the surfaces is the considerable simplification of the calculations for the analytically solved model.
As the gap is limited by the ring faces, the boundary conditions for the rotor and the stator are as follows:
λ ϑ z z = L = λ f θ f z z = 0   and   ϑ = θ f ,  
λ f θ f z z = h = λ ϑ z z = 0   and   θ f = θ ,   respectively
For the outer surface of the rotor ( r o ), heat transfer is assumed to be by convection, and it can be expressed by:
λ ϑ r r = r o = α f ϑ r = r o ,   for   the   rotor
λ θ r r = r o = α f θ r = r o ,   for   the   stator .
where α f is the convection coefficient.
All the above boundary conditions are necessary to solve the system of three differential equations.

4. Problem Solution

For the case considered here, a cylindrical coordinate system corresponding to the geometry of the physical model was used, and thus Equation (3) can be written as:
α ϑ τ α = 1 r ϑ r + 2 ϑ r 2 + 2 ϑ z 2 ,   where   τ = κ t
Equation (18) was solved by adopting the Marchi–Zgrablich transform, the finite Fourier cosine transform, and the Laplace transform. The general form of the finite integral transform by Marchi-Zgrablich [34] is:
H f x = f ¯ n = r i r o r f x S p λ , α ,   k n r d r
The inverse integral transform proposed by Marchi-Zgrablich can be expressed as:
H 1 f ¯ x = f x = n = 1 a n S p λ , α ,   k n r ,  
where
a n = f ¯ p n C n ,
C n = r i r o r S p λ , α ,   k n r 2 d r .
Thus, Equation (18) is written as:
α ϑ ¯ τ α = 2 ϑ ¯ z 2 k 2 n ϑ ¯ .
The finite Fourier cosine transform is described in the general form as:
F c g x = g m = 0 L g x cos m π x L d x ,  
where m = 1 , 2 , 3 ,
F c d 2 g x d x 2 = 1 m d g x d x x = L d g x d x x = 0 m 2 π 2 L 2 g m g n = F c g x .  
The inverse transform is given in the general form as:
F c 1 g m = g x = g m = 0 L + 2 L m = 1 g m cos m π x L .  
Applying the finite Fourier cosine transform to write Equation (23) and assuming the boundary conditions given in Equation (14), we obtain:
α ϑ ¯ τ α = φ ¯ q λ m 2 π 2 L 2 ϑ ¯ m k 2 n ϑ ¯ m .
With the Laplace transform [25] being:
L D C α f t = s α f ^ s k = 0 n 1 f k 0 + s α 1 k , n 1 < α < n ,  
we can rewrite Equation (27) as:
s α ϑ ¯ ^ s = φ ¯ q λ m 2 π 2 L r 2 ϑ ¯ ^ s k 2 n ϑ ¯ ^ s .
Once the transformations are completed, Equation (18) has the following form:
ϑ ¯ ^ s = φ ¯ q λ 1 ω 2 1 s s α 1 s α + ω 2   where   m 2 π 2 + L 2 k 2 n L 2 = ω 2   and   ω 0 2 = k 2 n
After the three inverse integral transforms are employed, the equation describing the temperature distribution for the rotor takes the final form:
T r = T o + 1 L n = 1 r i r o r φ ¯ q S p λ s , α f ,   k n r d r λ r i r o r S p λ , α ,   k n r 2 d r S p λ , α ,   k n r 1 ω 0 2 1 ω 0 2 E α κ ω 0 2 t α + + 2 L n = 1 r i r o r φ ¯ q S p λ s , α f ,   k n r d r λ r i r o r S p λ , α ,   k n r 2 d r S p λ , α ,   k n r m = 1 cos m π x L 1 ω 2 1 ω 2 E α κ ω 2 t α ,
where E α , β z is the function of the Mittag–Leffler type [35,36]:
E α , β z = v = 1 z v Γ α v + β   α > 0 ,   β > 0 .
For α , β = 1 , E α , β z = e z , and thus Equation (31) is solved using the classical heat equation with predetermined initial and boundary conditions.
S p λ , α ,   k n r = k n Y 1 k n r i + α Y 0 k n r o λ k n Y 1 k n r o J 0 k n r + k n J 1 k n r i α J 0 k n r o + λ k n J 1 k n r o Y 0 k n r
Equation (33) is dependent on the boundary conditions assumed for the rotor or stator model. The distribution of temperature in the stator is defined as:
T s = T 0 + n = 1 cosh k n z r i r o r θ f S p λ s , α f ,   k n r d r cosh k n L r i r o r S p λ s , α f ,   k n r 2 d r S p λ s , α f ,   k n r .
where θ f is the excess temperature of the medium in the gap, calculated as θ f = T f T o .

5. Operating Conditions

The parameters characterizing the material used for the rotor and the stator are provided in Table 1.
The parameters concerning the seal geometry and performance are shown in Table 2. It is assumed that the surrounding fluid is water with a temperature of 20 °C.
The analytical solution to the mathematical model considered here shows the distribution of temperature in the cross-section of the seal when α in Equation (18) has a fractional order value of 1 α 2 . Moreover, the parameter α is time-dependent, which is illustrated graphically in Figure 1.

6. Results and Discussion

The simulations were performed to determine the influence of the factor α on the distribution of temperature in the seal (Figure 4).
Figure 4 illustrates the distributions of temperature in the rotor over time, ranging from 0.1 s to 1.5 s. When the classical heat equation is used (Figure 4a), the temperature of the rotor increases gradually because of the heat flux generated in the gap. The greatest difference between the diagrams is observed for the time range 0.1–1.25 s. With the assumption that α t , Figure 1 suggests that if values higher than t 0 = 1 s are used, the order of the derivative tends toward unity.
The heat flux provided by the sealing ring face at a given moment of time for t = 0.05 (s) causes a local increase in temperature of 2.1 °C, which may result in thermal deformations of the sealing ring in contact with the gap. In the case of the fractional differential equation (for t = 0.05 (s)), the equation is hyperbolic in nature, as illustrated in Figure 5b.
The solution based on the variable-order derivative (VOD) time fractional model provides a better correlation between the predicted data and the observed results than the classical approach.

7. Validation of Results

The results were validated using Ansys Workbench software (Figure 6).
Figure 6 depicts the temperature distributions calculated using the Ansys Workbench model (Figure 6a) and those obtained analytically (Figure 6b).
Figure 7 compares the distributions of temperature along the radius (r = 0.0525 m).
From Figure 7, it is apparent that the results obtained using both calculation methods are in good agreement.

8. Conclusions

This article has presented a mathematical model to determine the heat transfer phenomena occurring in the non-contacting face seal used in a turbo machine. The problem was solved analytically by applying three integral transforms. The resulting relationships described the changes in temperature in the face seal cross-section. The boundary conditions were taken into consideration. The classical Fourier law was extended by applying the time fractional derivative where α = α t . This model can be used to describe unstable heat transfer conditions or thermal shock.
Even small local changes in temperature of an order of 2–10 °C can cause thermoelastic deformations of the stator and the rotor. As a result, there is a change in the gap geometry. During the first moments of the startup, unstable conditions can be observed because the seal operates under dry friction conditions; thus, the heat flux generated between the rotor and the stator is much greater than during the non-contact operation of the seal assumed in this model. The deformations may contribute to a greater leak.
It should be mentioned that the mathematical analysis based on the fractional differential equation is suitable to develop more accurate models to analyze similar physical phenomena.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Li, C.; Zheng, L.; Zhang, X.; Chen, G. Flow and heat transfer of a generalized Maxwell fluid with modified fractional Fourier’s law and Darcy’s law. Comput. Fluids 2016, 125, 25–38. [Google Scholar] [CrossRef]
  2. Bartosik, A. Application of Rheological Models in Prediction of Turbulent Slurry Flow. Flow Turbul. Combust. 2010, 84, 277–293. [Google Scholar] [CrossRef]
  3. Kozior, T. The Influence of Selected Selective Laser Sintering Technology Process Parameters on Stress Relaxation, Mass of Models, and Their Surface Texture Quality. 3D Print. Addit. Manuf. 2020, 7, 126–138. [Google Scholar] [CrossRef]
  4. Bochnia, J.; Blasiak, S. Fractional relaxation model of materials obtained with selective laser sintering technology. Rapid Prototyp. J. 2019, 25, 76–86. [Google Scholar] [CrossRef]
  5. Warbhe, S.D.; Tripathi, J.J.; Deshmukh, K.C.; Verma, J. Fractional Heat Conduction in a Thin Circular Plate with Constant Temperature Distribution and Associated Thermal Stresses. J. Heat Transf. 2017, 139, 44502. [Google Scholar] [CrossRef]
  6. Povstenko, Y.; Kyrylych, T. Fractional heat conduction in solids connected by thin intermediate layer: Nonperfect thermal contact. Contin. Mech. 2019, 31, 1719–1731. [Google Scholar] [CrossRef] [Green Version]
  7. Kaczorek, T. Selected Problems of Fractional Systems Theory; Springer Berlin Heidelberg: Berlin/Heidelberg, Germany, 2011; ISBN 9783642205026. [Google Scholar]
  8. Laski, P.A. Fractional-order feedback control of a pneumatic servo-drive. Bull. Pol. Acad. Sci.-Tech. Sci. 2019, 67, 53–59. [Google Scholar] [CrossRef]
  9. Zmarzly, P.; Gogolewski, D.; Kozior, T. Design guidelines for plastic casting using 3D printing. J. Eng. Fibers Fabr. 2020, 15. [Google Scholar] [CrossRef]
  10. Adamczak, S.; Zmarzly, P.; Kozior, T.; Gogolewski, D. Assessment of Roundness and Waviness Deviations of Elements Produced by Selective Laser Sintering Technology. In Engineering Mechanics 2017; ACAD SCI Czech Republic, INST Thermomechanics: Prague, Czech Republic; ISBN 978-80-214-5497-2.
  11. Takosoglu, J.E.; Laski, P.A.; Blasiak, S. Innovative Modular Pneumatic Valve Terminal with Self-Diagnosis, Control and Network Communications. In Engineering Mechanics 2014; Fuis, V., Ed.; ACAD SCI Czech Republic, INST Thermomechanics: Prague, Czech Republic; pp. 644–647. ISBN 978-80-214-4871-1.
  12. Nowacki, W. Problems of thermoelasticity. Prog. Aerosp. Sci. 1970, 10, 1–63. [Google Scholar] [CrossRef]
  13. Nowacki, W. Dynamic Problems of Thermoelasticity; Noordhoof International: Leyden, CO, USA, 1976; ISBN 978-90-286-0045-4. [Google Scholar]
  14. Nowacki, W. Thermoelasticity, 2nd ed.; Elsevier Science: Kent, UK, 2014; ISBN 978-1-4831-6248-5. [Google Scholar]
  15. Povstenko, Y.Z. Fundamental Solutions to Time-Fractional Advection Diffusion Equation in a Case of Two Space Variables. Math. Probl. Eng. 2014, 2014, 705364. [Google Scholar] [CrossRef]
  16. Povstenko, Y. Fractional Heat Conduction in an Infinite Medium with a Spherical Inclusion. Entropy 2013, 15, 4122–4133. [Google Scholar] [CrossRef] [Green Version]
  17. Povstenko, Y. Fractional heat conduction in a semi-infinite composite body. Commun. Appl. Ind. Math. 2014, 6. [Google Scholar] [CrossRef]
  18. Raslan, W.E. Application of fractional order theory of thermoelasticity to a 1D problem for a spherical shell. J. Theor. Appl. Mech. 2016, 54, 295–304. [Google Scholar] [CrossRef]
  19. Liu, J.-G.; Yang, X.-J.; Feng, Y.-Y. On integrability of the time fractional nonlinear heat conduction equation. J. Geom. Phys. 2019, 144, 190–198. [Google Scholar] [CrossRef]
  20. Koca, I. Modeling the heat flow equation with fractional-fractal differentiation. Chaos Solitons Fractals 2019, 128, 83–91. [Google Scholar] [CrossRef]
  21. Lotfy, K. Analytical solution of fractional order heat equation under the effects of variable thermal conductivity during photothermal excitation of spherical cavity of semiconductor medium. Waves Random Complex Medium 2021, 31, 239–254. [Google Scholar] [CrossRef]
  22. Liu, S.; Feng, L. An Inverse Problem for a Two-Dimensional Time-Fractional Sideways Heat Equation. Math. Probl. Eng. 2020, 2020, 5865971. [Google Scholar] [CrossRef] [Green Version]
  23. Povstenko, Y. Time-fractional radial heat conduction in a cylinder and associated thermal stresses. Arch. Appl. Mech. 2012, 82, 345–362. [Google Scholar] [CrossRef] [Green Version]
  24. Povstenko, Y.Z. Axisymmetric Solutions to Time-Fractional Heat Conduction Equation in a Half-Space under Robin Boundary Conditions. Int. J. Differ. Equ. 2012, 2012, 154085. [Google Scholar] [CrossRef] [Green Version]
  25. Jiang, X.; Xu, M. The time fractional heat conduction equation in the general orthogonal curvilinear coordinate and the cylindrical coordinate systems. Phys. A 2010, 389, 3368–3374. [Google Scholar] [CrossRef]
  26. Povstenko, Y. Fractional Thermoelasticity; Springer International Publishing: Cham, Switzerland, 2015; ISBN 978-3-319-15334-6. [Google Scholar]
  27. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models; Imperial College Press: London, UK; Hackensack, NJ, USA, 2010; ISBN 978-1-84816-329-4. [Google Scholar]
  28. Blasiak, M.; Blasiak, S. Application of Fractional Calculus in Harmonic Oscilator. In Engineering Mechanics 2017; ACAD SCI Czech Republic, INST Thermomechanics: Prague, Czech Republic, 2017; pp. 146–149. ISBN 978-80-214-5497-2. [Google Scholar]
  29. Blasiak, M.; Blasiak, S. The Application of Integral Transforms to Solving Partial Differential Equations of the Fractional Order. In Engineering Mechanics 2017; ACAD SCI Czech Republic, INST Thermomechanics: Prague, Czech Republic, 2017; pp. 150–153. ISBN 978-80-214-5497-2. [Google Scholar]
  30. Li, Z.L.; Sun, D.G.; Han, B.H.; Sun, B.; Zhang, X.; Meng, J.; Liu, F.X. Response of viscoelastic damping system modeled by fractional viscoelastic oscillator. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2017, 231, 3169–3180. [Google Scholar] [CrossRef]
  31. Qi, H.-T.; Xu, H.-Y.; Guo, X.-W. The Cattaneo-type time fractional heat conduction equation for laser heating. Comput. Math. Appl. 2013, 66, 824–831. [Google Scholar] [CrossRef]
  32. Blasiak, S. Time-fractional heat transfer equations in modeling of the non-contacting face seals. Int. J. Heat Mass Transf. 2016, 100, 79–88. [Google Scholar] [CrossRef]
  33. Tournerie, B.; Danos, J.C.; Frene, J. Three-Dimensional Modeling of THD Lubrication in Face Seals. J. Tribol. 2001, 123, 196–204. [Google Scholar] [CrossRef]
  34. Ghonge, B.E.; Ghadle, K.P. Deflection of transient thermoelastic circular plate by Marchi-Zgrablich and Laplace integral transform technique. Appl. Mech. Lett. 2012, 2, 21004. [Google Scholar] [CrossRef] [Green Version]
  35. Haubold, H.J.; Mathai, A.M.; Saxena, R.K. Mittag-Leffler Functions and Their Applications. J. Appl. Math. 2011, 2011, 298628. [Google Scholar] [CrossRef] [Green Version]
  36. Rogosin, S. The Role of the Mittag-Leffler Function in Fractional Modeling. Mathematics 2015, 3, 368–381. [Google Scholar] [CrossRef]
Figure 1. Fractional derivative order α vs. time.
Figure 1. Fractional derivative order α vs. time.
Energies 14 05512 g001
Figure 2. Non-contacting face seal made up of the flexibly mounted stator (1), the rotor (2), the spring (3), the housing (4), the O-ring (5), and the shaft (6).
Figure 2. Non-contacting face seal made up of the flexibly mounted stator (1), the rotor (2), the spring (3), the housing (4), the O-ring (5), and the shaft (6).
Energies 14 05512 g002
Figure 3. Boundary conditions for the heat transfer in the non-contacting face seal. (a) Heat flux, (b) convection, and (c) insulated surfaces.
Figure 3. Boundary conditions for the heat transfer in the non-contacting face seal. (a) Heat flux, (b) convection, and (c) insulated surfaces.
Energies 14 05512 g003
Figure 4. Distribution of temperature in the cross-section of the non-contacting face seal at different values of α. (a) for α = 1 , (b) for α t = var .
Figure 4. Distribution of temperature in the cross-section of the non-contacting face seal at different values of α. (a) for α = 1 , (b) for α t = var .
Energies 14 05512 g004aEnergies 14 05512 g004b
Figure 5. (a) Distribution of temperature in the seal for t = 0.05 (s) and α = 1.998, (b) temperature vs. the rotor thickness for r = 0.0545 (m).
Figure 5. (a) Distribution of temperature in the seal for t = 0.05 (s) and α = 1.998, (b) temperature vs. the rotor thickness for r = 0.0545 (m).
Energies 14 05512 g005
Figure 6. Temperature distribution in the rotor calculated by (a) the Ansys Workbench model and (b) the proposed model for t = 300 s.
Figure 6. Temperature distribution in the rotor calculated by (a) the Ansys Workbench model and (b) the proposed model for t = 300 s.
Energies 14 05512 g006aEnergies 14 05512 g006b
Figure 7. Distributions of temperature along the rotor height (thickness) for r = 0.0525 m.
Figure 7. Distributions of temperature along the rotor height (thickness) for r = 0.0525 m.
Energies 14 05512 g007
Table 1. Properties of the rotor and stator materials.
Table 1. Properties of the rotor and stator materials.
PropertiesDensity
ρ (kg/m3)
Poisson’s Ratio νThermal Conductivity
λ (W/mK)
Thermal Expansion
τ   1 / ° C
Specific Heat
Cp (J/kgK)
Material
Silicon carbide31000.18130 5 · 10−6750
Resin-impregnated carbon18600.2015 4 · 10−6
Table 2. Parameters related to seal geometry and performance.
Table 2. Parameters related to seal geometry and performance.
Seal GeometrySeal Performance
Inner radius
ri (m)
0.050Convection coefficient (water)
αf (W/m2 K)
18,000
Outer radius
ro (m)
0.055Temperature of the surrounding fluid
To (°C)
20
Rotor/Stator thickness
L (m)
0.010Angular speed
ω (rad/s)
800
Gap height
h (m)
2 · 10−6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Blasiak, S. Heat Transfer Analysis for Non-Contacting Mechanical Face Seals Using the Variable-Order Derivative Approach. Energies 2021, 14, 5512. https://doi.org/10.3390/en14175512

AMA Style

Blasiak S. Heat Transfer Analysis for Non-Contacting Mechanical Face Seals Using the Variable-Order Derivative Approach. Energies. 2021; 14(17):5512. https://doi.org/10.3390/en14175512

Chicago/Turabian Style

Blasiak, Slawomir. 2021. "Heat Transfer Analysis for Non-Contacting Mechanical Face Seals Using the Variable-Order Derivative Approach" Energies 14, no. 17: 5512. https://doi.org/10.3390/en14175512

APA Style

Blasiak, S. (2021). Heat Transfer Analysis for Non-Contacting Mechanical Face Seals Using the Variable-Order Derivative Approach. Energies, 14(17), 5512. https://doi.org/10.3390/en14175512

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