1. Introduction
The evaluation of the lifetime cost of a seal, which, in theory, is extremely long [
1], must take into account not only the initial cost of acquisition, but also all expenses that occur throughout its operation [
2] (operating costs, replacement of parts, power consumption, fluid losses, cost of stocking spare parts, production losses, auxiliary aggregates, secondary defects, product contamination, and environmental pollution).
Identifying the pairs of materials from which the two main sealing components are made may have a positive impact on optimizing energy–economy–environment systems and ultimately contribute to global energy efficiency.
In particular, this happens by increasing the life expectancy of the seals, which automatically generates lower operating costs, reduces maintenance costs, reduces losses from equipment downtime during repair, avoids the risk of other equipment failures, and directly contributes to the reduction in CO2.
This research is related to the water industry in Romania, particularly to drinking water and sewage pumping systems. These systems, in addition to the initial costs of production, assembly, commissioning, and energy consumption for their operation, have annual expenses generated by maintenance. An important component of maintenance activity includes the replacement of mechanical seals, considered as wear parts.
According to a previous report, the market size of mechanical seals was valued at 3.64 billion EUR in 2022 [
3].
The costs of replacing the mechanical seals of the pumps and returning them to original operating conditions are summarized in
Table 1. We took into account the preventive maintenance only (which represents considerable costs every year).
Other costs that are essential for such an intervention were not included (bearings, painting, resistant abrasive solutions, other axle repairs, transport, etc.). We narrowed our focus to only consider the sewage pumps of a municipal sewage treatment plant serving over 300,000 inhabitants made by only one in four or five pump manufacturers (which can be up and running at any given moment); there are an even greater number of products and pumps than those listed below within the treatment plant.
In
Table 1, by adding all existing pumps at the city level, we can extrapolate that the derived costs and effects would be significantly larger.
SC Hidropump SRL, the local company in charge of the maintenance and repair of these pumps, provided initial data gathered from repair reports.
Table 1 shows the costs of preventive maintenance for a 1-year period. The costs for other manufacturers’ pumps and other equipment with mechanical seals incorporated in their design were not subject of this analysis. This also applies to the cost related to industrial consumption (which is massive in terms of equipment seals and energy consumption), pumps from residential and commercial housing, municipal buildings, water supply stations, fire pumping stations, or the problems generated by nonoperational equipment.
The CO
2 emissions due to the production of these seals are another important factor not considered in the current analysis [
4]. On a local level, the focus of public administrations on reducing CO
2 emissions can be seen in various projects on sustainable energy, such as the Action Plan for Sustainable Energy of the Municipality of Timișoara (Planul de Acţiune pentru Energia Durabilă a Municipiului Timişoara) [
5].
This research can make strong contributions over time in cutting down maintenance expenses, as well as in terms of a social, economic, and environmental impact with respect to sustainable energy.
The performance of front seals can influence the overall energy efficiency. In the future, a thorough evaluation of each mentioned sector can make an important contribution to reducing energy consumption [
3].
For small seals, 1–2 inches (25–50 mm) in diameter, we have a relatively low energy consumption of 1 hp (0.75 kW). However, larger seals can consume more than 3 hp (2.24 kW) of power [
6].
Analyzing the data within the table, we can find a total energy loss of 8.6 EUR/h for the conditions stated above, while the loss in inefficient systems can be as much as 4%.
There are other types of rotary shaft seals; however, according to
Table 2, mechanical face seals (MFSs) are clearly superior. A comparative analysis of MFSs and other rotary shaft seals based on technological, operating, and cost criteria is presented.
In the literature [
8,
9,
10,
11], problems of frictional couplings are analyzed, respectively, of MFSs under various aspects.
From a functional point of view, MFSs are considered by researchers and engineers zero-leakage seals. While being used to seal fluids that vary from water to propane propylene or diesel [
12], they are found to work successfully with both gas and liquids in extreme conditions such as those of aero-engines upstream pumps [
13] or other volatile media such as hydrocarbon mixtures in petrochemical plants if the thermophysical properties of the sealed medium are well expressed [
14]. Guidelines for estimating the wear of the MFS when functioning under friction instability conditions are found in the literature [
15].
From an analytical point of view regarding the thermal behavior, the importance of the thermal conductance factors is found in the literature [
16]; models for the heat fluxes between fluids and solids, enthalpy of the fluid, and the vapor mass fraction, using Newtons method [
17] and mathematical analysis for non-contacting (gas lubricated) MFSs [
18] develop more detailed mathematical models describing physical phenomena.
The tribological aspects of MFS are well explored in the literature, from the influence of the waviness of the rings to the viscosity–temperature effect in the main contact area of the MFS [
19,
20], the common tribological concerns of the influence of surface roughness on temperatures [
21], and research on reducing temperatures in dry contact through surface texturing with graphite [
22,
23,
24] to the potential use of superconducting magnetic force fields to improve carrying capacity and lubrication characteristics under heavy loading conditions [
25].
Thermodynamic analyses carried out with the finite element method (FEM) on dry friction couplings of general mechanical seals [
26] and general thermodynamic FEM simulations are found in the literature [
27], but there is a lack of research on MFSs with lubricated interfaces.
From an experimental perspective, we found data about various pairs of materials tested for the friction and wear behavior, in dry sliding contact seals or lubricated [
28,
29].
In conclusion, the thermodynamic aspects of the MFS analyzed with FEM for various pairs of materials of the primary seal of the MFS, and the connection between them and the wear, reliability and economic implications in the operation fields of the equipment in which MFS are found, are not addressed in the literature.
Active control of the MFS will overcome the limitation of reliability, allowing a particularly reliable seal and a priori estimation of wear over time. Increasing the life of MFSs offers great economic benefits in top fields, such as nuclear, aeronautical, chemical, petrochemical, pharmaceutical, food, domestic, and transportation fields [
30].
2. Fundamentals of Thermo-Hydrodynamics of an MFS in the Non-Isothermal Stationary Regime Analyzed by FEM
A realistic approach to the MFS phenomenology can only be made through a thermo-hydrodynamic (THD) or even thermo-elasto-hydrodynamic (TEHD) analysis [
27].
The decisive consequences of optimizing the operation of an MFS resulting from finding the compromise between high durability and low leakage rates are as follows:
A very small thickness of the fluid film that leads to the production of an appreciable amount of heat in the interface of the primary seal;
The thermoelastic deformations of the rings, although reduced in absolute value, are comparable to the thickness of the film, significantly influencing the geometric configuration of the interface.
The main source of heat production in an MFS is friction from the interface (
Figure 1). Heat removal from the seal is mainly achieved by the movement of the sealed fluid.
If the film in the interface is heated above the boiling point at the pressure in the film, then it will evaporate causing malfunction and eventually the seal is destroyed.
An elementary numerical analysis shows that the role of the leakage flow in cooling the interface is insignificant so that the entire amount of heat is practically taken over by the two rings of the primary seal.
In this context, an essential aspect of the heat transfer study is represented by the quantities produced in the interface between the pressure ring and the friction ring. If the temperature of the sealed fluid is known, the modeling of the temperature field can be restricted to the two rings with the consideration of convection with the fluid.
If the correlation of interface geometry with the thermal deformations of the rings (THD problem) is taken into account, then in the system of equations of MFS modeling with fluid film, the equations of thermo-elasticity, associated with the condition of resting the two rings, must be added.
Instabilities of any kind are the main cause of premature retirement of an MFS that theoretically works with fluid film [
7]. The effect of temperature manifests itself on three levels, all of which significantly influence the operation of the MFS:
By the phase change of the fluid in the interface (vaporization);
On the opening force through the variation of fluid viscosity;
On the elastic deformations of the rings.
All analytical THD models of fluid film MFSs resort to a simplification of the heat flux line trajectories in the pressure ring of the seal and neglect the heat conduction in the friction ring.
From this it follows that the thermal conductivity of the material of the friction ring would be unimportant.
Extending the use of FEM as an effective tool for analyzing the behavior of structures from a mechanical and thermal point of view has as a consequence the continuous enrichment of knowledge of the intimate phenomena that govern the behavior of MFS.
The MFS assembly indicates the same position for thermal and mechanical loading. Because of this, geometric symmetries will also be symmetries of thermal load.
When determining by calculation using FEM, there is the advantage of determining the temperature in a large number of internal points of the primary sealing rings, the results depending on the extent to which the contour conditions were adopted according to the actual operating conditions.
3. Theoretical Bases of the FEM in Two-Dimensional and Axially Symmetric Problems of Heat Transfer in the Stationary Regime
3.1. Two-Dimensional Problems
The steady-state heat exchange of a body with the environment, in the case of two-dimensional problems (
Figure 2), is expressed by the equation:
Heat exchange by convection:
with the boundary conditions:
In Equation (1):
ϴ—temperature;
λx, λy—thermal conductivity coefficients of the material from which the body is made;
M—unit thermal flow of internal heat sources;
q—unit thermal flow through the separation surface S2;
α—the coefficient of heat exchange by convection between the surface S2 and the surrounding environment;
ϴe—the temperature of the surrounding environment.
According to the theories of the calculus of variations, solving the differential Equation (1) is equivalent to the minimization of a functional quantity, to the determination of the values of the parameters
ϴ, respectively, for which the functional quantity takes a minimum value. The appropriate function in this case is
The first term of the functional corresponds to the heat exchange by conduction. In the differential element is the constant thickness of the body, which is usually taken as equal to 1. The second term corresponds to the heat exchange by convection through the S3 surface, and the third term corresponds to the heat exchange by the imposed thermal flux through the S2 surface.
The minimization of functional 1 is performed on finite elements, the overall effect being obtained by cumulating the partial effects corresponding to each finite element in part. Of the many types of finite elements used today, the triangular ones have the advantage of simplicity.
The temperature variation on the surface of the finite element is described by the shape or interpolation functions, which are characteristic of each type of finite element. These functions express the temperature on the element relative to the temperatures in the nodes. For the triangular finite element (
Figure 2):
in which the shape functions are linear relationships of the type:
Coefficients
a,
b,
c, are constant quantities, in the expression of which only the coordinates of the nodes of the finite element enter:
| | |
| | |
| | |
Minimizing the functional 1 means canceling the derivative of the functional in relation to the parameter
ϴ, that is
In the following, it is shown how the structural thermal balance equation was established starting from this condition.
Deriving 1 with respect to
ϴi and equating to zero, we obtain
From 2 and 3, it follows that
substituting in 5, we obtain
It was assumed that on the sides of the element that exchange heat by convection, α = constant and ϴe = constant, and on the sides on which heat is exchanged by imposed flow, q = constant.
Next, the terms of the two members of Equation (7) will be examined successively.
Left member. Term 1 represents the matrix of thermal conductivity coefficients. Assuming that the limits of a finite element
λx and
λz are constant and taking into account that
, the first term is written in the form:
Term 2 represents the contribution due to the heat exchange through convection.
In this expression, the integrals
I0…
I5 have the following values:
where:
is the length of the side that exchanges heat through convection, and h is the thickness of the body (
Figure 3).
The left member of Equation (7) is the first line of the matrix [k]e of the element.
The other lines are obtained analogously by making
and
, so that in the end the left member becomes
The right member. Term 1 represents the heat flux of the internal heat sources:
Assuming that on the surface of a finite element
M = const., and evaluating the shape function in the center of gravity of the coordinate element
, i.e., expressing the shape function as
, a sufficiently accurate value is obtained in the form:
which expresses the equal distribution of the flow
M on the three nodes of the element. If, instead of the unitary thermal flow of the internal sources, the flow distributed on nodes
M1,
M2,
M3 is known, then
so that the first term of the right member is written in the form:
Term 2 represents the convective exchange through the lateral surface:
Term 3 represents the heat exchange through the imposed heat flux:
Overall, for the entire element, the second member of Equation (7) is a column vector of the form:
Equating the left member with the right one, respectively (
a) with (
b), we obtain the thermal balance equation for a finite element:
Summing for all
m finite elements, the thermal balance equation for the entire section (body) is obtained:
or
where [
K] is the thermal conductivity matrix, {
ϴ} is the column-vector of the nodal temperatures, and {
F} the column-vector of the internal and external thermal fluxes in all nodes of the section. Solving the system of Equation (9) leads to the determination of the temperature in all nodes of the section.
3.2. Axisymmetric Problems
In the case of these problems, which include bodies with axial geometric and thermal symmetry, the functional to be minimized is
The finite element is a torus with a triangular cross-section (
Figure 4).
In the expression of the functional, the radius r appears variable on the surface of the finite element. Sufficiently accurate results are obtained using the radius at the center of gravity, . As the differential element is , from we can obtain .
More accurate results are obtained by using, instead of
r, the expression
R, given by
If the phenomenon of thermal radiation is also taken into account, the general form of the equation of the thermal balance in steady state becomes
where:
[K]—thermal conductivity matrix;
[R]—matrix of heat exchange by radiation;
{u}—vector of unknown temperatures;
Tabs—the absolute output data of the temperatures required for the calculations of heat exchange by radiation;
{P}—the vector of the applied constant heat flow;
{N}—the vector of temperatures dependent on the heat flows.
The equation is non-linear due to the presence of four terms of the power distribution law developed by radiation.
In addition, nonlinearities due to matrix coefficients and boundary conditions may appear.
These nonlinearities are introduced by specifying material properties and temperature-dependent boundary conditions.
The NASTRAN software applies the Newton–Raphson iteration scheme to solve these nonlinear equations.
The procedure leads to the following form of the thermal balance equation:
where:
[KT]i—tangential conductivity matrix dz/du;
[KT]i = [K]I + 4[R]i{ui + Tabs}3 − {dN/du}I;
{R}i—the residual vector;
{R}i = {P}I + {N}i − [K]{u} − [R]i{ui + Tabs}4.
At each iteration, the matrix in the left member and the vector in the right member are calculated on the basis of the temperature vector, solving for the unknown vector the new temperature values:
The program aims to achieve convergence solutions in the optimal sense, balancing various aspects of the solutions, including loading, prior residual data, prior tangent matrix data, etc.
4. Materials and Methods
From a constructive and functional point of view, the MFS seals by means of a class IV tribological coupling (coupling with plane surfaces) [
10].
The analyzed MFSs used in this paper are manufactured by ROSEAL S.A., namely EFS 103-S and EF 19 × 45.4 × 3.3. They are used in water pumps and are presented in
Figure 5.
The dimensions of MFS EFS-102-S used in the thermodynamical FEM analysis are presented in
Table 3.
The stationary ring commonly known as the stationary seat (6) (
Figure 5) is immobilized in relation to the housing and constitutes the friction element.
The floating ring (5), in pressed frontal contact with the seat because of the load applied by the spring (3), has the role of a pressure element and rotates together with the shaft.
For reasons of constructive proximity, the sealing function in this case is distributed on two levels. The first level is the primary sealing ensured by pressed and direct contact between the two rings with relative motion. The second level is the secondary sealing, represented by the rubber sleeve to the shaft (1) and to the rubber gasket to the housing (7); they prohibit leaks on the areas of immobile contact with the conjugate parts (bearing cover and shaft).
It is noticed that, in this way, the solution can achieve very high performances conveniently influencing each area.
There are a number of factors that lead to the wear and stress of a front seal, such as contact corrosion, overload, thermal overload, and wear.
Temperature. It appears in the interface of the primary seal due to the heat produced by friction.
Higher temperatures can cause severe wear, gasification of the lubricant film which results in a rapid increase in the coefficient of friction and wear, in some cases exceeding the sintering temperatures of materials, which has the effect of breaking, reciprocal welding or thermal cracking of the rings.
Thermal deformations. The temperature in the interface and its gradients also influence the geometry of the interface [
11].
In the case of elastic deformations, the size of the modulus of elasticity of the couple of materials used and the dimensions of the rings are determining factors. The quantifying values of the thermal properties of the materials, i.e., the thermal conductivity (λ), thermal expansion coefficient (α), and heat transfer coefficient (k), combined with the construction of the rings, in turn, influence the temperature gradient and, therefore, the shape of the gap.
Deformation of sealing surfaces. The effects of hydrostatic pressures and temperature differences cause deformations of the rings and of the sealing surfaces, respectively [
30].
The thermal effect is very important in the assessment of the total conicity of the gap because the increase in speed causes the increase in power which, in turn, determines the increase in conicity caused by the uneven heating of the interface.
Two of the factors that influence the level of friction of an annular coupler can be easily highlighted: the couple of materials in contact, namely, the friction rate () and the geometry of the interface.
In conclusion, the stable operation of an MFS can be modified when the thermal and mechanical distortions that can occur are combined with the angular misalignment.
The materials of the primary seal rings used by manufacturers and their physical, mechanical and thermal properties are presented in
Table 4 [
8,
32].
The materials are coded [
32] as follows:
G10—stainless steel AISI 316;
U2—silicon carbide (SiC);
Y2—polytetrafluoroethylene PTFE + 25% graphite powder;
B1—carbon graphitic material bonded with resin.
The modeling of frictional contact between the sealing components takes into account the relative movement that occurs between them during operation (n = 3000 rot/min).
The fluid to be sealed is water at a temperature of 80 °C under a pressure of 0.1 MPa.
The corresponding friction rates accepted in the modeling are those evaluated following the measurement of the roughness of the active surfaces of the friction ring coupling [
7].
The thermodynamical simulations were made on the two models of MFS earlier mentioned. They were modeled in CAD and imported into the finite element software MSC Nastran.
FEM is a numerical method that can be used to find the solutions to a wide range of engineering problems such as stress analysis, heat transfer analysis, electromagnetic field analysis, or fluid flow analysis.
The method is based on the idea of constructing complicated objects from simpler objects, or the division of complicated objects into simpler objects, called finite elements [FE], for which known calculation schemes can be applied. The operation of dividing the real model into parts of small dimensions is called meshing [
33].
In our study we focused on the thermodynamical phenomenon that occurs between the two main elements of the MFS.
The operating temperature of FE with boundary friction in case of natural cooling results from the thermal balance equation [
34]:
where:
k—heat transfer coefficient;
AT—the free surface (heat exchanger) of the primary seal.
Friction regime and thermal equilibrium are interdependent processes [
34]. As a result, the effective rate of friction (
μ) and the regime temperature of the friction ring coupling (
tf) are determined iteratively depending on the variation of the lubricants viscosity (
ν) with temperature (Δ
t) [
9].
In conclusion, the stable operation of a frontal seal can be modified when the thermal and mechanical distortions that can occur are combined with angular misalignment.
The model of the frontal seal made with axially symmetrical finite elements.
From a geometrical point of view, the rings of the primary seal of the MFS have axial symmetry, being composed of rotating bodies (
Figure 6). This symmetry means that the modeling problem can be reduced to the modeling of the half section of the seal.
The models analyzed with FEM correspond dimensionally and constructively to the component parts of the frontal seals (
Figure 5).
The results of the analyses are illustrated in
Figure 7a (EFS 103-S) and
Figure 7b (EF 19 × 45.4 × 13.3) where we can see the meshed parts.
The analysis of the model is repeated for pairs of materials from the 4 materials presented in
Table 4 and coded by the manufacturer.
The contact between the floating and stationary rings is described by a new “slide line” element type that includes the stiffness/friction data, the scale factor and the friction coefficient between the contacting surfaces. The geometric definition of the surfaces in contact includes a series of nodes called “master” and “slave” which belong to the finite sliding surfaces of the deformable bodies and which were defined together with the model discretization.
The “slide line” element is the interaction element between the bodies in contact; its representation in the analyzed model appears as the connecting line between the “master” and “slave” nodes. The coordinate system attached to the “slide line” element is defined by the local coordinate systems of the “master” and “slave” nodes. Each contact zone in the joint assembly is defined by a slide line element.
5. Results
Simulation results. The modeling reproduces the operating behavior of the MFS, respective to EFS 130-S and EF 19 × 45.4 × 13.3. Axial pressure applied to the floating ring will induce a contact pressure field in the interface of the friction ring coupling of the seals. Additionally, this determines the state of tension and the displacements that occur in the front sealing ring.
From a thermal point of view, modeling requires knowledge of:
The map of the temperatures and the heat flow in the primary sealing rings for the two types of MFS and for different pair of materials from the FEM analysis are illustrated in
Figure 8,
Figure 9,
Figure 10 and
Figure 11.
After the postprocessing, the distributions of temperatures and the thermal loading along the interface of the primary sealing are also obtained (
Figure 12,
Figure 13 and
Figure 14).
Based on the parameters determined by FEM analysis, the leakage rates can be evaluated. In the case of the non-isothermal flow of an incompressible fluid for an axial-symmetrical configuration [
35,
36], the Reynolds equation is
where:
r = d/2;
h—the dynamic viscosity of the working fluid.
Considering that at the limit, in the mentioned case, the interface leakage flow rate is given by the relation:
where:
dm = 0.5(d + D);
b = 0.5(D − d);
h = 0.5(Rz1 + Rz2);
Δp—the pressure drop in the interface;
Rz1/2—the contact roughness;
The dynamic viscosity varies depending on the temperature (T) of the fluid according to a polynomial dependence:
where:
ηf = the dynamic viscosity at reference temperature.
The power dissipated through friction in the interface is
where:
μ—the friction rate in the interface;
AT = π(D2 − d2)/4
ω—the angular velocity of the floating ring
And it transforms into heat:
where:
λ1/2—the thermal conductivity of the ring materials;
l1/2—the axial dimensions of the rings.
The debit relation becomes
or
From relations (20) and (21), we can see the direct dependence between the leakage flow rate Q, the product (pν) and the heat transfer coefficient kT, values that, in turn, characterize the operation of a seal, of a sliding bearing, respective to a heat exchanger.
Customizing relation (20) for the MFS analyzed and taking into account the FEM results, the flow rate-angular velocity dependences are obtained (
Figure 15 and
Figure 16).
The results of the analysis highlight the fact that the flow rate of leaks varies linearly with increasing the angular speed of the pressure ring and is independent of the materials used in the primary seal.
The values of leakage losses due to the leakiness of the MFS analyzed are far below the values allowed by Richtlinie VDI 2440, namely 1 cm
3/h (12 g/h) [
37].
6. Discussion
From a thermal point of view, the physical-mechanical and thermal properties of the working fluid and the materials of the primary seal of the MFS have a major influence on the operating behavior of the MFS. Thus, the materials of the primary seal must have a low coefficient of friction, to allow a quick evacuation of the heat produced, to have a low coefficient of expansion and not to damage each other.
From the FEM analysis, the results are as follows:
The maximum temperature is reached in the forecast area of the interface of the primary seal;
We can determine the limits of the rings’ ability to take over a certain power lost through friction without exceeding certain temperatures (targeting the choice of materials);
With the help of the model, it is possible to quickly and precisely determine the extreme temperatures of the rings;
The temperature distribution allows the analytical evaluation of the performance of the front seals (see the expression for calculating the loss flow rate);
The axial distribution of the temperature is not linear as approximated in the literature [
36,
38,
39];
The temperature gradient in the axial direction is a decisive factor in the thermal deformation of the rings with a direct effect on the geometry of the interface;
The temperature variation in the radial direction is insignificant;
There is a large temperature gradient in the axial direction near the interface;
Taking also into account the phenomenon of convection at the level of the pressure ring/medium, the relatively high values of the heat transfer coefficient kT resulting from the thermal analysis are observed.
The FEM results also allow the variation of the heat transfer coefficient kT in the axial direction to be determined and the value of the convection coefficient starting from the power is dissipated.
Until now, the focus has been on electric motor efficiency and this is seen in the trends in energy efficiency technologies widely used around the world (energy efficient motors, energy saving soft starters, variable speed drives, energy efficient transformers, etc.) This paper presents a new approach by taking a step further looking for energy efficiency methods resulting from extending the life of MFS [
40].
7. Conclusions
Increasing the life expectancy of the seals automatically generates lower operating costs, low maintenance costs, low losses from equipment downtime during repair, avoiding the risk of other equipment failures.
A good part of the maintenance activity also includes the replacement of mechanical seals, considered as wear parts. According to the report “Global Mechanical Seals Market Size by Design (Pusher Type Mechanical Seals, Non-Pusher Type Mechanical Seals), By End User Industry (Oil and Gas, Chemical, Mining), By Geographic Scope and Forecast,” the mechanical seals market size was valued at EUR 3.64 billion in 2022.
Active control of the MFS will overcome the limitation of reliability, allowing a particularly reliable seal and a priori estimation of wear over time. Increasing the life of MFS offers great economic benefits in the top fields, such as nuclear, aeronautical, chemical, petrochemical, pharmaceutical, food, domestic, and transportation fields.
There are a number of factors that lead to the wear and stress of a front seal, such as contact corrosion, overload, thermal overload, and wear. Temperature appears in the interface of the primary seal due to the heat produced by friction. The temperature in the interface and its gradients also influence the geometry of the interface.
A realistic approach to the MFS phenomenology can only be made through a thermo-hydrodynamic analysis.
From a thermal point of view, the physical-mechanical and thermal properties of the working fluid and the materials of the primary seal of the MFS have a major influence on the operating behavior of the MFS.
In the paper, the thermodynamic FEM analysis was made using the MSCNASTRAN software and was carried out for two types of MFS, namely EFS-103-S and EF 19 × 45.4 × 3.3, with an operation field of n = 3000 rot/min, water at 80 °C, pressure of 0.1 MPa and for various pairs of materials of the primary seal.
The results highlight the following aspects:
The maximum temperature is reached in the forecast area of the interface of the primary seal;
The limits of the rings’ ability to take over a certain power lost through friction without exceeding certain temperatures (targeting the choice of materials);
It is possible to quickly and precisely determine the extreme temperatures of the rings;
The temperature distribution allows the analytical evaluation of the performance of the front seals, namely the loss flow rate;
The axial distribution of the temperature is not linear as approximated in the literature;
The temperature gradient in the axial direction is a decisive factor in the thermal deformation of the rings with a direct effect on the geometry of the interface;
The temperature variation in the radial direction is insignificant; there is a large temperature gradient in the axial direction near the interface;
Additionally, taking into account the phenomenon of convection at the level of the pressure ring/medium, relatively high values of the heat transfer coefficient kT resulting from the thermal analysis are observed.
The FEM results also allow the determination of the variation of the heat transfer coefficient kT in the axial direction and the value of the convection coefficient starting from the power dissipated.
The conclusions are that energy saving solutions are available and there is potential to save many labor hours and waiting time caused by maintenance procedures, replace parts, and increase overall efficiency (factors that can all be quantified in energy according to the first law of thermodynamics) by simply choosing better pairs of materials, even if they have an initial higher cost.
Considering the global tendencies in energy demand and availability, researchers recommend using reliable design solutions whenever possible and FEM analysis in the detriment of experimental research because of the possibility of saving time and resources, not only in the case of MFSs, but also in any endeavor.
We recommend that energy efficiency techniques emphasize studying the link between the MFS (focused on increasing the lifetime of the seal) and the energy consumption of electric motors, and how significant savings can be achieved in this sector of energy, taking into account that, in our study, we have a total energy loss of 8 euros/hour (according to
Table 1).