When a malfunction occurs in a marine main engine system, the impact of the anomaly will propagate through the system, affecting the performance of all relevant components in the system. The phenomenon of fault propagation in the system caused by induced factors can interfere with fault localization, making the latter a difficult task to solve. This paper aims at showing how the "characteristic curves method" is able to properly locate malfunctions also when more malfunctions appear simultaneously. To this end, starting from the working principle of each component of a real marine diesel engine system, comprehensive and reasonable thermal performance parameters are chosen to describe their characteristic curves and include them in a one-dimensional thermodynamic model. In particular, the model of a low-speed two stroke MAN 6S50 MC-C8.1 diesel engine is built using the AVL Boost software and obtaining errors lower than 5% between simulated values and test bench data. The behavior of the engine is simulated considering eight multi-fault concomitant phenomena. On this basis, the fault diagnosis method proposed in this paper is verified. The results show that this diagnosis method can effectively isolate the fault propagation phenomenon in the system and quantify the additional irreversibility caused by the Induced factors. The fault diagnosis index proposed in this paper can quickly locate the abnormal components.
Citation: Nan Xu, Longbin Yang, Andrea Lazzaretto, Massimo Masi, Zhenyu Shen, YunPeng Fu, JiaMeng Wang. Fault location in a marine low speed two stroke diesel engine using the characteristic curves method[J]. Electronic Research Archive, 2023, 31(7): 3915-3942. doi: 10.3934/era.2023199
Related Papers:
[1]
Xiaoyan Wu, Guowen Ye, Yongming Liu, Zhuanzhe Zhao, Zhibo Liu, Yu Chen .
Application of Improved Jellyfish Search algorithm in Rotate Vector reducer fault diagnosis. Electronic Research Archive, 2023, 31(8): 4882-4906.
doi: 10.3934/era.2023250
[2]
Xu Chen, Wenbing Chang, Yongxiang Li, Zhao He, Xiang Ma, Shenghan Zhou .
Resnet-1DCNN-REA bearing fault diagnosis method based on multi-source and multi-modal information fusion. Electronic Research Archive, 2024, 32(11): 6276-6300.
doi: 10.3934/era.2024292
[3]
Jichen Hu, Ming Zhu, Tian Chen .
The nonlinear observer-based fault diagnosis method for the high altitude airship. Electronic Research Archive, 2025, 33(2): 907-930.
doi: 10.3934/era.2025041
[4]
Zhenzhong Xu, Xu Chen, Linchao Yang, Jiangtao Xu, Shenghan Zhou .
Multi-modal adaptive feature extraction for early-stage weak fault diagnosis in bearings. Electronic Research Archive, 2024, 32(6): 4074-4095.
doi: 10.3934/era.2024183
[5]
Jian Yuan, Hao Liu, Yang Zhang .
Automatic detection method of abnormal vibration of engineering electric drive construction machinery. Electronic Research Archive, 2023, 31(10): 6327-6346.
doi: 10.3934/era.2023320
[6]
Zhuanzhe Zhao, Mengxian Wang, Yongming Liu, Zhibo Liu, Yuelin Lu, Yu Chen, Zhijian Tu .
Adaptive clustering algorithm based on improved marine predation algorithm and its application in bearing fault diagnosis. Electronic Research Archive, 2023, 31(11): 7078-7103.
doi: 10.3934/era.2023359
[7]
Yu Chen, Qingyang Meng, Zhibo Liu, Zhuanzhe Zhao, Yongming Liu, Zhijian Tu, Haoran Zhu .
Research on filtering method of rolling bearing vibration signal based on improved Morlet wavelet. Electronic Research Archive, 2024, 32(1): 241-262.
doi: 10.3934/era.2024012
[8]
Chunyan An, Wei Bai, Donglei Zhang .
Meet-in-the-middle differential fault analysis on Midori. Electronic Research Archive, 2023, 31(11): 6820-6832.
doi: 10.3934/era.2023344
[9]
Zhizhou Zhang, Yueliang Pan, Weilong Zhao, Jinchu Zhang, Zheng Zi, Yuan Xie, Hehong Zhang .
Frequency analysis of a discrete-time fast nonlinear tracking differentiator algorithm based on isochronic region method. Electronic Research Archive, 2024, 32(9): 5157-5175.
doi: 10.3934/era.2024238
[10]
Xingyue Liu, Kaibo Shi, Yiqian Tang, Lin Tang, Youhua Wei, Yingjun Han .
A novel adaptive event-triggered reliable H∞ control approach for networked control systems with actuator faults. Electronic Research Archive, 2023, 31(4): 1840-1862.
doi: 10.3934/era.2023095
Abstract
When a malfunction occurs in a marine main engine system, the impact of the anomaly will propagate through the system, affecting the performance of all relevant components in the system. The phenomenon of fault propagation in the system caused by induced factors can interfere with fault localization, making the latter a difficult task to solve. This paper aims at showing how the "characteristic curves method" is able to properly locate malfunctions also when more malfunctions appear simultaneously. To this end, starting from the working principle of each component of a real marine diesel engine system, comprehensive and reasonable thermal performance parameters are chosen to describe their characteristic curves and include them in a one-dimensional thermodynamic model. In particular, the model of a low-speed two stroke MAN 6S50 MC-C8.1 diesel engine is built using the AVL Boost software and obtaining errors lower than 5% between simulated values and test bench data. The behavior of the engine is simulated considering eight multi-fault concomitant phenomena. On this basis, the fault diagnosis method proposed in this paper is verified. The results show that this diagnosis method can effectively isolate the fault propagation phenomenon in the system and quantify the additional irreversibility caused by the Induced factors. The fault diagnosis index proposed in this paper can quickly locate the abnormal components.
1.
Introduction
Due to the advantages of fuel economy and reliability, diesel engines have attained a leading position in marine applications, and the large two-stroke engine is widely used as the main power plant of civil ships [1,2].
To meet the requirements of Ocean transportation, Marine diesel engine has been developing towards the direction of intelligence and high efficiency [3,4], which also increases the complexity of diesel engine structure. This change increases the possibility of Marine diesel engine system failure. When a failure occurs in a marine main engine system, it will not only produce economic losses but also risk personnel safety
In the past decades, the maintenance of marine diesel engines has evolved from corrective actions to current trends. A lot of studies focus on predictive measures [5] that may improve the reliability of these engines. The purpose of these diagnostic systems is to detect and diagnose diesel engine anomalies before the anomalies cause undesired consequences.
A suitable diagnostic system requires a complete and reliable database that helps identify and diagnose anomalies when symptoms characterizing the abnormality are activated. In order to build a complete and reliable database, diesel engine simulation models have been used to reproduce the fault phenomena [6,7,8]. With the development of technology, real-time monitoring parameters are becoming more and more comprehensive [9].
A marine diesel engine has a complex structure. When a diesel engine breakdowns, a fault may unstably lead to various symptoms and a detected symptom may be caused by many different faults, which brings considerable difficulties to locate the real root causes of the diesel engine.
In recent years, machine learning and artificial intelligence technology have been widely used in the field of diesel engine fault diagnosis, which provides an effective tool to study the relationship between diesel engine faults and abnormal symptoms [17,18]. The deep learning method can use small data to develop fault diagnosis systems [19,20].
The fault tree method is an effective method to analyze the causal relationship between diesel engine faults and abnormal symptoms, where the diagnosis principle of this method is based on the expert's knowledge of the abnormalities [21,22].
The Marine diesel engine system needs to work in wet, vibratory and other harsh conditions for a long time, which leads to various faults of the Marine diesel engine system during operation. In extreme cases, multiple faults may occur at the same time. When multiple faults occur simultaneously in several different components, the coupling phenomenon between faults will affect the relationship between diesel engine faults and abnormal symptoms.
Therefore, there is a technical problem in the research field of diesel engine fault diagnosis to solve the problem of multi-fault diagnosis [23,24].
With the development of fault diagnosis based on the structural theory of thermoeconomics [25], the useful concepts of "intrinsic" and "induced" malfunctions have emerged [26,27], which may help better understand the causes and effects of failures. Compared with the theory of perturbations [28], the fault diagnosis method based on the thermoeconomic structure theory can obtain more accurate fuel impact values. Valero et al. [29] applied this method to the fault diagnosis of a 160MW coal-fired power plant and analyzed the influence of each component on the fuel consumption of the system. The authors used thermodynamic steady-state simulations to analyze the induced effects (induced malfunction and dysfunction) on other components caused by the intrinsic malfunction of a certain component, and they found that the highest additional fuel plant consumption was attributable to the component where the highest inefficiency occurs.
Verda et al. [30,31,32] studied the influence of the control system on thermoeconomic diagnosis using simulations to establish the corresponding "free condition" for different operating conditions. This "free condition" was mathematically determined using a specific model of the system. This particular condition is characterized by the same position of the governing parameters as that of the reference condition, but contains the anomalies occurring in the actual operating condition. On this basis, the influence of extra fuel consumption and irreversiblele losses caused by control system intervention is studied. It has been found that the control system intervention may generate higher induced malfunction in some components.
With the deepening of thermal-economic fault diagnosis, researchers [33] found that the complex interaction between components and the intervention of the control system may be the origin of induced malfunctions. These malfunctions are not actual faults within the corresponding components, but they appear because of the mutual interaction between the components' behavior. Thus, the reliability of this fault diagnosis technique may be strongly compromised if the induced malfunction cannot be identified and separated. However, the identification and separation of induced malfunctions are not easy.
Lazzaretto and Stoppato et al. [34] applied this method to the fault diagnosis of multiple complex energy systems to verify the reliability of the method. They found that when the degree of interaction between components in the system is low, this method is effectively able to locate the fault. On the other hand, when the system configuration is complex, this method cannot strictly guarantee the identification of the true source of the fault. The main reason is that it cannot isolate the propagation phenomenon of anomalies, and in turn cannot effectively identify the intrinsic malfunction and induced malfunction, because the variation of an operation parameter generated by the malfunction affects the exergy variables of both the components in which the malfunction occurs and those that suffer an induced malfunction.
On this basis, Toffolo and Lazzaretto et al [35] proposed a new method, named "Characteristic curve method" [36], to detect the malfunctions, based on the idea that the intrinsic malfunction of a component results in a change of its characteristic curve. Therefore, they proposed a new index to check whether the actual operating point of the component moves away from its reference operating condition along the characteristic curve or outside of it. Thus, it can be judged whether the anomaly of the component operation is due to anomalies in other system components (induced malfunctions) or to anomalies within the component (intrinsic malfunction).
In this paper, we propose a diagnosis method based on the characteristic curve of each component for isolating the propagation phenomena of faults in the system. The goal of this method is to realize the fault location of the diesel engine system When multiple faults occur simultaneously in several different components.
This paper describes the path and mode of fault propagation in the system. Then, based on the characteristics of each component in the diesel engine system, an accurate model of a marine diesel engine system is built after the selection of a suitable set of thermal performance parameters for the accurate description of the characteristic curve of each engine component. The method proposed in this paper is then verified by eight case studies for identifying the failure in the marine diesel engine system.
According to the manufacturer specifications, the diesel fuel used in the engine certification tests is Marine Gas Oil.
The model simulates the indicated cycle. The gas composition and thermodynamic properties in each element of the intake and exhaust systems piping discretization at any simulation time step (i.e., crank angle) are calculated by solving the one-dimensional conservation equation set commonly used in engine one-dimensional gas-dynamic simulations [37], i.e., mass Eq (1), momentum Eq (2) and energy Eq (3) conservation Eq (4) coupled with the ideal-gas constitutive equation, which writes:
∂ρ∂t+∂ρu∂x=0
(1)
∂ρu∂t+μ∂ρu∂x+∂p∂x=0
(2)
∂(ρe)∂t+[∂(ρuh)∂x]=0
(3)
pV=mRT
(4)
The average cell size is set to 100 (Target Average Cell Size for Spatial Pipe Discretization), and the convergence control is set in several components. The simulation calculation finishes when at least three engine cycles meet the convergence criteria. The convergence criterion is that the variation of the cycle averaged values ("transients") of the IMEP of each cylinder in BOOST™ elements over the last three consecutive cycles is less than a prescribed threshold (500 Pa). In addition, the fuel used in the test bench has the same chemical composition as the vessel.
The components of the model are shown in Figure 1. The system consists of the following components connected by pipes: air filter, compressor, air cooler, intake manifold, cylinders, exhaust manifold and turbine. In the simulation model, measuring points are located in the most relevant positions.
Figure 1.
The marine diesel engine model AVL BOOST 2020.
The components that make up the model are connected by unidimensional pipes. The pipe parameters include length, equivalent diameter, friction coefficientand wall temperature. The local pressure drop of the air filter is set to a fixed value in the range of 0~0.02 bar. The type of turbocharger is TCA66-21. The complete turbine and compressor maps are input into the marine diesel engine model through tools provided by the software to make sure that the turbocharging process can respond to the actual situation when the boundary conditions change or failure conditions are introduced in the model. The cooling capacity of the air cooler is determined by the inlet temperature of the cooling water and the efficiency of the air cooler. The simulation model considers the actual volume of intake and exhaust manifolds. The heat transfer coefficient of both intake and exhaust manifold models is 0.
Cylinders are defined in terms of their design dimensions, geometric compression ratio, combustion, heat transfer and scavenging port and exhaust valve data. The heat release rate of the combustion process is simulated by a Wiebe law [38,39,40,41]. Blow-by is not considered in this article. The lift curve of the scavenging port and exhaust valve is used to determine the filling and emptying of a cylinder. The Woschni 1990 heat transfer model is used to simulate convective heat losses.
2.2. Model validation
The model was verified by comparing the simulation data with the engine shop test data. Table 2 shows the comparison between the mean values of the simulation and experimental main parameters [42]. The deviation between simulation and experimental values is generally less than 5%. In addition, the verification of the relationship Speed-Power, Speed-SFOC, Speed-Boost Pressure and Speed-Turbocharger Outlet Temperature was also performed, as shown in Figures 2−5.
Table 2.
Comparison between simulated and experimental.
Therefore, it is considered that the model has good accuracy in simulating engine behavior, and can be used as an abnormal simulation platform without the need to generate failures in a real engine and with a consequent remarkable saving of fuel and time.
3.
Fault location method
3.1. Propagation phenomenon of failure
The mass and energy flow of the diesel engine system is shown in Figure 6. When a component of the system has a failure, the thermodynamic quantities of the mass and energy flow associated with the abnormal component are therefore altered. Due to the interaction among components, this modification affects the operations of other components in turn. Therefore, the impact of the anomaly will propagate through the system, affecting the performance of all relevant components in the system.
For better and easy understanding, the compressor failure (5% reduction in compressor isentropic efficiency) case is taken as an example. As shown in Figure 7, in the case of a compressor failure, the performance of other components in the system where the anomaly does not exist also changed. This phenomenon occurs because compressor failure changes the operating conditions of other components [36,43]. The failure of the compressor will change the mass flow, pressure and temperature of the air at the outlet of the compressor, resulting in changes in the operating conditions of the air cooler. Although there is no fouling or blockage in the air cooler, the operating state of the air cooler will change accordingly. Based on similar principles, the working performance of components such as cylinders and turbines will also change.
In other words, due to the interaction among components, the changes caused by the intrinsic malfunction can spread throughout the system, creating induced malfunction in the components where the anomaly does not exist. This phenomenon of fault propagation in the system can interfere with fault localization.
3.2. Definition of characteristic curve and diagnostic index
The failure or anomaly of a component (such as changes in the compressor blade geometry, blockage of an air cooler, deposits on heat exchange surfaces, etc.) will affect the functional relationships among the thermodynamic variables (temperature, pressure, mass flow rates, etc.) on which the mass and energy streams involved in the operation of that component depend. These functional relationships are often referred to as component characteristic curves (performance maps in machines, heat transfer models for heat exchangers).
As shown in Figure 8, the operating state of the component will deviate from its original reference state characteristic curve (broken line) and move from the reference state (point A) to a new state (point C). Due to the interaction among components, the alterations caused by the failure of the component will spread through the whole system, affecting the operating states of other components, which react to the changes imposed by the faults in another component according to their non-modified characteristic curves. As shown in Figure 8, the operating state of a component in perfect order will not deviate from the original characteristic curve and move from the reference state (point A) to a new state (point B). In other words, external factors or induced factors do not affect the characteristic curve of the component that can be affected only by internal factors. External, induced and internal factors are defined as follows [33]: External factors-modification of external conditions such as environment variables. Induced factors-variations caused by the change of other system variables, e.g., component interaction or the intervention of the control system. Internal factors-degradation or failures.
The location of malfunctioning components depends on the knowledge of the characteristic curve of each component in mathematical form. In particular, Eq (5) formalizes the characteristic curve of the component ith by a set of relationships f that define a performance parameter or a thermodynamic variable π characterizing the component behavior as a function of a subset δk of the independent variables involved in the component operation.
πi,ref=fi,ref(δi,refk)
(5)
On the other hand, it is always possible to quantify with πi,real the actual value of the performance parameter πi when the component operates at a specific real condition, and with Iiindex the difference between πi,real and the value of the performance parameter πi as expected from Eq (5) for the corresponding operating condition. The difference Iiindex writes as reported in Eq (6).
Iiindex=πi,real−fi,ref(δi,realk)
(6)
Therefore, during the real operation of the component under normal behavior of the system that embeds the component, the result of Eq (6) is zero, since the working point of the component ith, corresponds to the operating condition predicted in accordance with the reference characteristic curve. In contrast, when a performance degradation or fault changes the characteristic curve of the component, the result of Eq (6) will be different from zero. Thus, Iiindex can be used as a diagnostic index.
The characteristic curve of the component can be linearly approximated by using its derivative as calculated in the reference state (as shown in Figure 8, the tangent AB1 in the reference state, approximates the point B with the point B1, the slope of the tangent AB1 is the derivative ∂fi,ref/∂δik, at point A, the quantity ϵ stands for the "residual effects"). In mathematical terms, such approximation can be formalized as reported in Eq (7).
As shown in Figure 9, if the real operating point of the component is on the original characteristic curve, moving from reference state A to the new state B, the malfunction of the component is only an induced effect. In this case, neglecting the approximation introduced by the linearization of the characteristic curve, Iiindex = 0. In contrast, if the component has an intrinsic malfunction, the real operating point of the component will depart from the characteristic curve, moving from reference state A to the new state C, and the value of Iiindex is expected to be non-zero.
The main components in the diesel engine system can be divided into two categories. The first category includes all the "work components", i.e., the components whose product can be expressed as useful exergy. For example, the turbocharger compressor consumes the mechanical work provided by the turbine to increase the internal energy of the air. The second category includes all the "dissipative components", whose product cannot be expressed as useful exergy. For example, the air coolers dissipate heat and reduce the temperature of the cylinder intake air, and the dissipated internal energy is the target product of the air cooler.
In this paper, the irreversibility is selected as parameter π for the diagnostic index Iiindex of the components embedded into the diesel engine system. This because:
The irreversibility Ii includes the knowledge related to the mass and energy flows and directly reflects the efficiency of the energy conversion process.
For the "work component", the existence of anomalies will always increase the diagnostic index Iiindex.
For the "dissipative component", the existence of anomalies will always reduce the diagnostic index Iiindex.
In summary, if the component has an intrinsic malfunction (degradation or failure of the component), the indicator Iiindex of the component will change.
3.4. The characteristic parameters of component
The basic rule for the selection of variables δik is that they must be a set of independent variables characterizing the behavior of the component. Therefore, the component characteristic curves fi,ref, and so the derivatives ∂fi,ref/∂δik, can be defined as functions of δik[36].
For an energy system, all performance variables can be expressed as a function of thermodynamic variables, so thermodynamic variables are the actual independent variables of the energy system.
In addition, exergy represents a synthesis of thermodynamic information that is useful in describing the outcomes of component behavior.
As shown in Section 3.2 of this article, a component's performance is determined by its own physical constraints (characteristic curves) and operational parameters.
Therefore, the natural choice for these variables is a set of independent component thermodynamic variables (including mass flow, pressure and temperature). The number of independent variables is equal to the number of component degrees of freedom.
However, for components with complex production structures or energy conversion, such as cylinders of diesel engines, the required number of independent variables available cannot be monitored. Therefore, specific output parameters should be used to replace the available arguments that cannot be monitored.
As shown in Figure 6, exergy flowing into the compressor includes air exergy flow at compressor inlet Exa1 and mechanical exergy WC.
Exa1=f(mxa1,Txa1,Pxa1)
(10)
WC=f(nC,TC)
(11)
Therefore, the independent variable group representing the compressor characteristic curve includes mxa1, Txa1, Pxa1, nC and TC.
However, the signal monitoring of compressor speed and torque is difficult and costly. The relationship between compressor outlet thermal parameters and WC is shown in Eqs (12)–(15).
ηC=(Exa2−Exa1)WC
(12)
π=Pxa2Pxa1
(13)
Txa2=Txa1[1+1ηC(πk−1k−1)]
(14)
mxa1=mxa2
(15)
So, the independent variable group representing the compressor characteristic curve includes mxa1, Txa1, Pxa1, Txa2 and Pxa2.
The selection principle of the independent variable group of the air cooler is similar to the superheater [44]. Considering that the flow rate and pressure of the cooling water are constant. The independent variable group representing the air cooler characteristic curve includes mxa2, Txa2, Pxa2 and Tcoolant.
For the diesel engine system, the function of the cylinder is to use part of the high-pressure air sent by the compressor to burn the fuel and to mix the combustion products with the rest of the high-pressure air to form high-temperature and high-pressure gas, which drove the piston to produce mechanical exergy as shown in Figure 6. Therefore, the performance of the cylinder is also related to its mechanical structure (design parameters).
Take cylinder No. 1 as an example.
Exf1 indicates the chemical exergy of the fuel.
Exf1=mf1LHV
(16)
Exa4 represents air exergy flow at the inlet.
Exa4=f(mxa4,Txa4,Pxa4)
(17)
Exg1 indicates exhaust gas exergy flow.
Exg1=f(mxg1,Txg1,Pxg1)
(18)
Thermodynamic processes, chemical reactions, mechanical work, heat transfer, mass diffusion and friction processes occur simultaneously within this component. In other words, the irreversible losses in the working process of the cylinder must include irreversibility losses due to fuel combustion, cylinder heat transfer, exhaust gas, mechanical friction, etc. The description of each irreversibility loss requires monitoring a large number of characteristic parameters, such as compression end temperature, atomization of the fuel in the cylinder, etc., which are very demanding to be managed and some of the parameters are not easy to measure with the standard marine engine technology. On the other hand, the overall irreversibility loss in the cylinder can be calculated by the cylinder input energy and the cylinder output work.
The indicated mean effective pressure (IMEP) is the work output of one cycle for unit swept volume and relates to the cylinder power in accordance with Eq (19).
Wi=n60IMEPi(πSD24)
(19)
Wi——Power of the ith cylinder, [kW]
n——diesel engine speed, [r/min]
IMEPi——Mean Effective Pressure of the ith cylinder, [kPa]
D——bore, [mm]
S——stroke, [mm]
Therefore, for the diesel engine system, the independent variable group representing the cylinder No. 1 characteristic curve includes mxa4, Txa4, Pxa4, and mf1.
The selection principle of the independent variable group of the turbine is similar to the compressor. However, according to the differences between the compressor map and the turbine map, the independent variable group representing the turbine characteristic curve includes mxa7, Txa7, Pxa7 and Pxa8.
The characteristic parameters selected for the cylinder and all of the other engine components considered in this study are shown in Tables 3 to 6. It should be noted that the fuel quantity injected in each cylinder is assumed here as a constant value related to the diesel engine load because, in current diesel engine systems, it is difficult to monitor the flow rate in each fuel injector.
Table 3.
Characteristic parameters of the compressor.
The ΔIi, calc is the expected variation of component irreversibility due to a change Δδk of the component independent variables, according to the reference component characteristic curve.
Since the analytic form of each characteristic curve is generally unknown, the derivatives ∂Ii,ref/∂δik of every component can be calculated by generating several virtual operating conditions near the reference operation using a simulator.
In order to make the reference state close to the actual state, it is necessary to change the diesel engine speed and fuel injection quantity.
n=ne(WWe)1/3
(22)
We——Nominal engine power, [kW]
ne——Nominal engine speed, [r/min]
W——Diesel engine power, [kW]
n——Diesel engine speed, [r/min]
Δh=100geW60γLn
(23)
Δh——The amount of oil supplied to one cylinder in 100 engine cycles, [ml];
ge——Fuel consumption rate, [g/(kW·h)];
W——Diesel engine power, [kW];
γ——Density of fuel oil, 0.825[g/(cm3)] for diesel;
L——Number of cylinders;
n——Diesel engine speed, [rpm].
In the reference state, Eq (22) allows calculating the speed of the diesel engine, whereas Eq (23) allows estimating the amount of fuel injected in the cylinder.
Taking the turbine as an example, four virtual operating conditions (ref1, ref2, ref3 and ref4) are required to calculate the derivatives of the turbine, since the number of component independent variables is four (PMP10, TMP10, PMP11 and mMP10). as shown in Eq (24). Note that the turbine outlet pressure is an environmental variable (atmospheric pressure). To avoid the rank of the equations being less than the number of independent variables of the equation, resulting in countless solutions to the equation, at least two different turbine outlet pressures should be set among the four reference states generated by the simulation simulator.
Once the model has been adjusted and validated, failures are introduced one by one.
The compressor failure F1 is usually caused by dust accumulation in the impeller or diffuser as well as damages that produce changes in geometry. The compressor failure is simulated by reducing the isentropic efficiency of the compressor.
The air cooler failure F2 is usually caused by the increase of fouling on the inner wall of the air cooler, which will produce a reduction of cooling capacity. The air cooler failure is simulated by reducing the isentropic efficiency.
This kind of cylinder failure F3 (Blocking of the injector hole of the cylinder) is usually caused by the fault of the fuel injection system or the carbon accumulation of the nozzle, resulting in a reduction in the corresponding fuel mass flow rate.
This kind of cylinder failure F4 (Excessive blow-by) is usually caused by abnormal wear of the piston ring, which is simulated by increasing the clearance between the piston ring and the sleeve.
The turbine failure F5 is usually caused by dust accumulation in the impeller or diffuser as well as damages that produce changes in geometry. The turbine failure is simulated by reducing the isentropic efficiency of the turbine.
In this paper, the method has been verified by cases real 1~8 as shown in Table 7.
Table 7.
Cases real 1−8.
F1
F2
F3
F4
F5
Real 1
a 5% reduction in compressor efficiency
a 10% decrease in air cooler efficiency
Real 2
a 5% reduction in compressor efficiency
Cylinder No 1 a 10% reduction in fuel mass flow rate
Real 3
a 5% reduction in compressor efficiency
a 10% reduction in turbine efficiency
Real 4
a 10% decrease in air cooler efficiency
Cylinder No 1 a 10% reduction in fuel mass flow rate
Real 5
a 10% decrease in air cooler efficiency
a 10% reduction in turbine efficiency
Real 6
Cylinder No 1 a 10% reduction in fuel mass flow rate
a 10% reduction in turbine efficiency
Real 7
Cylinder No 2 a 6% reduction in fuel mass flow rate
Cylinder No 1 a 0.03mm increase in the gap
Real 8
a 5% reduction in compressor efficiency
a 5% decrease in air cooler efficiency
Cylinder No 1 a 4% reduction in fuel mass flow rate
The diagnosis results of cases real 1–8 operating conditions are listed in Tables 10−17, in which each row referring to the malfunctioning components is highlighted in bold.
For example, looking at Table 10, first, it is worth noting that all the failures are accurately located and identified. In particular, the expected irreversibility variation of all components is higher than 0, which means that the propagation of the induced effects caused by failures involves all components in the system.
In addition, it can be found in Table 16 that all the failures are accurately located and identified. However, this method cannot identify the type of failure (piston ring abnormal wear or fuel pump wear).
Focusing on the malfunctioning components (i.e., the compressors and air coolers), the value of the index Iiindex (see the last but one column in the table) is high. This means that the actual irreversibility of such components is noticeably different from the expected irreversibility, indicating that their characteristic curve has changed because of their intrinsic operation anomalies. In contrast, the corresponding values of index Iiindex for the other components are close to 0 (note that the non-zero value is mostly due to the error ε in the approximation of the derivative), which means that the actual irreversibility of these components is basically equal to the expected irreversibility, i.e., the components are affected by malfunction induced effects and their characteristic curves do not deviate from the original one.
Finally, since the exergy flow rate of different components of the diesel engine system varies greatly, the ratios Iiindex/|ΔIi, calc| that are reported in the last column of the table, are the more proper diagnostic indexes to take into consideration. Tables 9 to 16 show that this index amplifies the anomaly of the components in which the intrinsic malfunction occurs increasing the gap existing between these components and the others, which are affected only by the approximation error in the evaluation of the derivatives of the characteristic curves. Thus, this index greatly simplifies the failure location and identification process.
5.
Conclusions
This paper proposes a diagnosis method for diesel engines based on the characteristic curve of each component. When multiple faults occur simultaneously in several different components, the advantage of this method is to effectively locate the malfunctions, so as to maintain the reliability of the system. The author establishes and verifies the reliable model of the 6S50 engine and each of its components, and then describes the characteristic curve of each component. On this basis, eight multiple malfunctions cases have been diagnosed and discussed to show the ability of the proposed method.
The following conclusions can be drawn from this application:
1) The validity and reliability of the method were demonstrated by eight test cases. The method is able to effectively isolate the propagation phenomena of faults in the system, and accurately identify the location of faults in the marine diesel engine system.
2) The author selected comprehensive and practical thermal performance parameters (that can be monitored) to describe the characteristic curves of the main components (compressor, turbine, cylinder and air cooler). Therefore, the proposed fault location method has high practical application value.
3) Both the irreversibility in the ith component, Iiindex and Iiindex/|ΔIi, calc| and can be used as fault location indicators for the marine diesel engine system. Considering the exergy flow rate of different components of diesel engine systems varies greatly, it is more effective to choose Iiindex/|ΔIi, calc| as the fault location indicator.
4) The fault diagnosis method proposed in this paper can effectively locate malfunctions, but this method cannot identify the type of failure (piston ring abnormal wear or fuel pump wear). This method still needs to be combined with other fault identification methods to identify the type of failure.
Nomenclature
Acronyms
IMEP
Indicated Mean Effective Pressure (bar)
LHV
Lower Heating Value (J/kg)
Greek symbols
Δ
Increment
π
dependent variable
δ
independent variable
ϵ
residual effects
γ
Density of fuel oil [g/(cm3)]
ρ
density
u
velocity of fluid
ƞ
isentropic efficiency
Symbols
I
Irreversibility
Iindex
Diagnostic index
E
Exergy flow
n
Diesel engine speed (rpm)
P
Diesel engine power (kW)
m
Mass flow rate (kg/s)
p
Pressure (bar)
T
Temperature (K)
T
Torque (N·m)
f
characteristic curve function
W
Power (kW)
We
Nominal engine power (kW)
ge
Fuel consumption rate g/(kW·h)
Δh
The amount of oil supplied per 100 cycles (ml)
ne
Nominal engine speed (rpm)
D
Bore (mm)
S
Stroke (mm)
F
Failure type
C1
Air compressor
C2
Air cooler
C3~C8
Cylinder 1~6
C9
Turbine
L
Number of cylinders
t
time
x
means flow longitudinal dimension
e
means total internal energy
h
total enthalpy
R
gas constant
k
adiabatic exponent air
Subscripts
MP
Measure point location
ref
reference operating condition
real
real operating condition
C
Compressor
T
Turbine
xa
air exergy flow
xg
gas exergy flow
xf
fuel exergy flow
xq
cooling water exergy flow
calc
calculated
i
Index for numerating of components
Conflict of interest
The authors declare there is no conflict of interest.
References
[1]
V. T. Lamaris, D. T. Hountalas, A general purpose diagnostic technique for marine diesel engines–Application on the main propulsion and auxiliary diesel units of a marine vessel, Energy Convers. Manage., 51 (2010), 740–753. https://doi.org/10.1016/j.enconman.2009.10.031 doi: 10.1016/j.enconman.2009.10.031
[2]
D. W. Wang, L. Shi, S. P. Zhu, B. Liu, Y. H. Qian, K. Y. Deng, Numerical and thermodynamic study on effects of high and low pressure exhaust gas recirculation on turbocharged marine low-speed engine, Appl. Energy, 261 (2020), 114346. https://doi.org/10.1016/j.apenergy.2019.114346 doi: 10.1016/j.apenergy.2019.114346
[3]
J. Carlton, J. Aldwinkle, J. Anderson, Future Ship Powering Options: Exploring Alternative Methods of Ship Propulsion, London: Royal Academy of Engineering, 2013.
[4]
Brent, Haight, 2011 marine propulsion order survey: a review of mechanical drive, auxiliary and diesel-electric marine propulsion orders in 2010, Diesel & Gas Turbine Worldwide, 43 (2011), 30–30.
[5]
A. Jardine, D. Lin, D. Banjevic, A review on machinery diagnostics and prognostics implementing condition-based maintenance, Mech. Syst. Signal Process., 20 (2006), 1483–1510. https://doi.org/10.1016/j.ymssp.2005.09.012 doi: 10.1016/j.ymssp.2005.09.012
[6]
J. Galindo, J. R. Serrano, F. Vera, C. Cervello, M. Lejeune, Relevance of valve overlap for meeting Euro 5 soot emissions requirements during load transient process in heavy duty diesel engines, Int. J. Veh. Des., 41 (2006), 343–367, https://doi.org/10.1504/IJVD.2006.009675 doi: 10.1504/IJVD.2006.009675
[7]
G. Vera, P. Rubio, H. Grau, Improvements of a failure database for marine diesel engines using the RCM and simulations, Energies, 13 (2019), 104. https://doi.org/10.3390/en13010104 doi: 10.3390/en13010104
[8]
A. Jose, V. G. Francisco, H. G. Jose, M. C. Jose, A. H. Daniel, Marine diesel engine failure simulator based on thermodynamic model, Appl. Therm. Eng., 144 (2018), 982–995. https://doi.org/10.1016/j.applthermaleng.2018.08.096 doi: 10.1016/j.applthermaleng.2018.08.096
[9]
R. Pawletko, S. Polanowski, Evaluation of current developments and trends in the diagnosis of marine diesel engines based on the indicator diagrams analysis, J. KONES,21(2014), 389–396. https://doi.org/10.5604/12314005.1130492 doi: 10.5604/12314005.1130492
[10]
A. J. Bayba, D. N. Siegel, K. Tom, Application of Autoassociative Neural Networks to Health Monitoring of the CAT 7 Diesel Engine, Army Research Laboratory, 2012.
[11]
D. T. Hountalas, Prediction of marine diesel engine performance under fault conditions, Appl. Therm. Eng., 20 (2000), 1753–1783. https://doi.org/10.1016/S1359-4311(00)00006-5 doi: 10.1016/S1359-4311(00)00006-5
[12]
J. Rubio, F. Vera-García, Sistema de Diagnóstico de Motor Diesel Marino Basado en Modelo Termodinámico y de Inteligencia Artificial, Ph.D thesis, Universidad Politécnica de Cartagena, 2017. https://doi.org/10.13140/RG.2.2.32226.63688
[13]
X. W. Li, Study on fault diagnosis of narine diesel engine, Comput. Simul., 12 (2012), 21–32. https://doi.org/10.1108/eb052507 doi: 10.1108/eb052507
[14]
Z. Jian, Simulation research on EGR reducing Nox emission of diesel engine, Int. J. Energy Power Eng., 4 (2015), 275–279. https://doi.org/10.11648/j.ijepe.20150405.16 doi: 10.11648/j.ijepe.20150405.16
[15]
C. Iclodean, N. Burnete, Computer simulation of ci engines fuelled with biofuels by modelling injection iRate law, Res. J. Agric. Sci., 44 (2012), 249–257.
[16]
T. Firsa, AVL Boost simulation of engine performance and emission for compressed natural gas direct injection engine, J. Energy Environ., 6 (2014).
[17]
O. C. Basurko, Z. Uriondo, Condition-Based maintenance for medium speed diesel engines used in vessels in operation, Appl. Therm. Eng., 80 (2015), 404–412. https://doi.org/10.1016/j.applthermaleng.2015.01.075 doi: 10.1016/j.applthermaleng.2015.01.075
[18]
G. Qi, Z. Zhu, K. Erqinhu, Y. Chen, Y. Chai, J. Sun, Fault-diagnosis for reciprocating compressors using big data and machine learning, Simul. Modell. Pract. Theory, 80 (2018), 104–127. https://doi.org/10.1016/j.simpat.2017.10.005 doi: 10.1016/j.simpat.2017.10.005
[19]
Z. Zhu, Y. Lei, G. Qi, Y. Chen, Y. Chai, Y. An, et al., A review of the application of deep learning in intelligent fault diagnosis of rotating machinery, Measurement, 206 (2022), 112346. https://doi.org/10.1016/j.measurement.2022.112346 doi: 10.1016/j.measurement.2022.112346
[20]
X. Huang, G. Qi, N. Mazur, Y. Chai, Deep residual networks-based intelligent fault diagnosis method of planetary gearboxes in cloud environments, Simul. Modell. Pract. Theory, 116 (2022), 102469. https://doi.org/10.1016/j.simpat.2021.102469 doi: 10.1016/j.simpat.2021.102469
[21]
V. Kneevi, J. Orovi, L. Stazi, J. Ulin, Fault tree analysis and failure diagnosis of marine diesel engine turbocharger system, J. Mar. Sci. Eng., 8 (2020), 1004. https://doi.org/10.3390/jmse8121004 doi: 10.3390/jmse8121004
[22]
Z. Jiang, D. Wei, L. Wang, Z. Zhao, J. Zhang, Fault diagnosis of diesel engines based on a classification and regression tree (CART) decision tree in Chinese, J. Beijing Univ. Chem. Technol., 45 (2018), 71–75. https://doi.org/10.13543/j.bhxbzr.2018.04.013 doi: 10.13543/j.bhxbzr.2018.04.013
[23]
F. Elamin, Y. Fan, F. Gu, A. Ball, Diesel engine valve clearance detection using acoustic emission, Advances in Mechanical Engineering, 2 (2010), 495741. https://doi.org/10.1155/2010/495741 doi: 10.1155/2010/495741
[24]
K. Chen, Z. Mao, H. Zhao, J. Zhang, A variational stacked autoencoder with harmony search optimizer for valve train fault diagnosis of diesel engine, Sensors, 20 (2019), 223. https://doi.org/10.3390/s20010223 doi: 10.3390/s20010223
[25]
A. Valero, C. Torres, L. Serra, A general theory of thermoeconomics. Part Ⅰ: Structural analysis, in 1992 International Symposium ECOS, ASME, New York, USA, (1992), 137–145.
[26]
A. Valero, C. Torres, F. Lerch, Structural theory and thermoeconomic diagnosis. Part Ⅲ: intrinsic and induced malfunctions, in 1999 International Symposium ECOS, ASME, New York, USA, (1999), 8–10.
[27]
A. Valero, C. Torres, L. Serra, Structural theory and thermoeconomic diagnosis. Part Ⅰ: On malfunction and dysfunction analysis, Energy Convers. Manage., 43 (2002), 1503–1518. https://doi.org/10.1016/S0196-8904(02)00032-8 doi: 10.1016/S0196-8904(02)00032-8
[28]
A. Valero, M. A. Lozano, J. L. Bartolomé, On-line monitoring of power-plant performance, using exergetic cost techniques, Appl. Therm. Eng., 16 (1996), 933–948. https://doi.org/10.1016/1359-4311(95)00092-5 doi: 10.1016/1359-4311(95)00092-5
[29]
A. Valero, F. Lerch, L. Serra, Structural theory and thermoeconomic diagnosis. Part Ⅱ: Application to an actual power plant, Energy Convers. Manage., 43 (2002), 1519–1535. https://doi.org/10.1016/S0196-8904(02)00033-X doi: 10.1016/S0196-8904(02)00033-X
[30]
V. Vittorio, L. M. Serra, V. Antonio, Effects of the productive structure on the results of the thermoeconomic diagnosis of energy systems, Int. J. Thermodyn., 5 (2002), 127–137. https://doi.org/10.5541/ijot.95 doi: 10.5541/ijot.95
[31]
V. Vittorio, L. M. Serra, A. Valero, Zooming procedure for the thermoeconomic diagnosis of highly complex energy systems, Int. J. Thermodyn., 5 (2002), 75–83. https://doi.org/10.5541/ijot.90 doi: 10.5541/ijot.90
[32]
V. Vittorio, L. M. Serra, V. Antonio, The effects of the control system on the thermoeconomic diagnosis of a power plant, Energy, 29 (2004), 331–359. https://doi.org/10.1016/j.energy.2003.10.003 doi: 10.1016/j.energy.2003.10.003
[33]
A. Stoppato, A. Lazzaretto, Exergetic analysis for energy system diagnosis, in 1996 Biennial JointConference on Engineering Systems Design and Analysis, ASME, New York, USA, (1996), 191–198.
[34]
A. Lazzaretto, A. Toffolo, A critical review of the thermoeconomic diagnosis methodologies for the location of causes of malfunctions in energy systems, J. Energy Res. Technol., 128 (2003), 345–354. https://doi.org/10.1115/1.2358148 doi: 10.1115/1.2358148
[35]
A. Stoppato, C. Carraretto, A. Mirandola, A diagnosis procedure for energy conversion plants: Part Ⅰ—Description of the method, in ASME International Mechanical Engineering Congress and Exposition, American Society of Mechanical Engineers, New York, USA, (2001), 493–500. https://doi.org/10.1115/IMECE2001/AES-23658
[36]
A. Stoppato, C. Carraretto, A. Mirandola, A diagnosis procedure for energy conversion plants: Part Ⅱ—Application and results, in ASME International Mechanical Engineering Congress and Exposition, American Society of Mechanical Engineers, New York, USA, (2001), 501–508. https://doi.org/10.1115/IMECE2001/AES-23659
[37]
A. Toffolo, A. Lazzaretto, A new thermoeconomic method for the location of causes of malfunctions in energy systems, J. Energy Res. Technol., 129 (2007), 1–9. https://doi.org/10.1115/1.2424960 doi: 10.1115/1.2424960
[38]
A. Lazzaretto, A. Toffolo, R. Passuello, The characteristic curve method in energy systems diagnosis: analysis of uncertainties in a real plant, in ASME International Mechanical Engineering Congress and Exposition, ASME, New York, USA, (2005), 379–390.
H. M. Nahim, R. Younes, H. Shraim, M. Ouladsine, Oriented review to potential simulator for faults modeling in diesel engine, J. Mar. Sci. Technol., 21 (2016), 533–551. https://doi.org/10.1007/s00773-015-0358-6 doi: 10.1007/s00773-015-0358-6
[41]
G. Radica, Expert system for diagnosis and optimisation of marine diesel engines, Strojarstvo, 50 (2008), 105–116.
[42]
P. Karpiński, K. Pietrykowski, L. Grabowski, Turbocharging the aircraft two-stroke diesel engine, Combust. Engines, 178 (2019). https://doi.org/10.19206/CE-2019-319 doi: 10.19206/CE-2019-319
[43]
G. Cong, T. Gerasimos, C. Hui, Analysis of two stroke marine diesel engine operation including turbocharger cut-out by using a zero-dimensional model, Energies, 8 (2015), 5738–5764. https://doi.org/10.3390/en8065738 doi: 10.3390/en8065738
[44]
L. Wang, P. Fu, N. Wang, T. Morosuk, Y. Yang, G. Tsatsaronis, Malfunction diagnosis of thermal power plants based on advanced exergy analysis: The case with multiple malfunctions occurring simultaneously, Energy Convers. & Manage., 148 (2017), 1453–1467. https://doi.org/10.1016/j.enconman.2017.06.086 doi: 10.1016/j.enconman.2017.06.086
This article has been cited by:
1.
Hong Je-Gal, Young-Seo Park, Seong-Ho Park, Ji-Uk Kim, Jung-Hee Yang, Sewon Kim, Hyun-Suk Lee,
Time-Series Explanatory Fault Prediction Framework for Marine Main Engine Using Explainable Artificial Intelligence,
2024,
12,
2077-1312,
1296,
10.3390/jmse12081296
2.
Guolong Li, Yanjun Li, Site Li, Shengdi Sun, Haotong Wang, Jian Su, Jianxin Shi, Xin Zhou,
Research on anomaly detection of steam power system based on the coupling of thermoeconomics and autoencoder,
2025,
318,
03605442,
134819,
10.1016/j.energy.2025.134819
Nan Xu, Longbin Yang, Andrea Lazzaretto, Massimo Masi, Zhenyu Shen, YunPeng Fu, JiaMeng Wang. Fault location in a marine low speed two stroke diesel engine using the characteristic curves method[J]. Electronic Research Archive, 2023, 31(7): 3915-3942. doi: 10.3934/era.2023199
Nan Xu, Longbin Yang, Andrea Lazzaretto, Massimo Masi, Zhenyu Shen, YunPeng Fu, JiaMeng Wang. Fault location in a marine low speed two stroke diesel engine using the characteristic curves method[J]. Electronic Research Archive, 2023, 31(7): 3915-3942. doi: 10.3934/era.2023199