Next Article in Journal
Lanthanum Diffusion in Fluorapatite at 400 °C, 50 MPa and 4–16 wt %
Previous Article in Journal
New Concept for the Study of the Fluid Dynamics of Lithium Extraction Using Calix[4]arene Derivatives in T-Type Microreactor Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental and Theoretical Research on Pressure Drop Changes in a Two-Stage Air Filter Used in Tracked Vehicle Engine

Faculty of Mechanical Engineering, Military University of Technology, 2 gen. Sylwestra Kaliskiego St., 00-908 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Separations 2021, 8(6), 71; https://doi.org/10.3390/separations8060071
Submission received: 24 April 2021 / Revised: 12 May 2021 / Accepted: 17 May 2021 / Published: 21 May 2021

Abstract

:
The effect of mineral dust in the air sucked in by an engine on accelerated component wear and reduction in performance was presented. The necessity to use two-stage air filters (multicyclone-paper insert) for military vehicles was shown. The results showed that placing an air filter in the path of the air entering the engine causes an additional pressure drop (air filter resistance increase), which leads to engine power decrease and increased fuel consumption. An analysis of model filter beds’ pressure drop changes (depending on bed parameters, aerosol flow parameters, and dust content) was carried out. It was revealed that it is very difficult to model changes in pressure drop in filter beds for actual conditions that appear during vehicle operation. The air filter pressure drop measurement results of more than 20 tracked vehicles operating in variable air dust concentration conditions were presented. The forms of selected regression models of the “life curve” type, best suited to the actual changes in air filters pressure drop as a function of the vehicle mileage, were determined. Significant differences were found between the same model values for different units of the tested vehicles. The quality of forecasting pressure drop value by selected functions was assessed by extrapolating them to the value of the next measurement and comparing the forecast and actual value. It was found that for the performed experiment, sufficiently good results of experimental data approximation and forecasting were obtained for a simple linear model.

1. Introduction

Military motor vehicles and working machines are mainly operated off-road, in sandy and desert terrain, where mineral dust concentration in the air is particularly high and often exceeds 1 g/m3. This is very important for trucks and special vehicles, whose power units are high-power diesel engines with large (1000–6000 m3/h) volumetric air demand. A PT-91 tank engine travelling at a velocity of 20 km/h on military training area roads, where concentration of dust suspended in the air is 1 g/m3, sucks in more than 170 kg dust with the air within 1000 km (approximately 50 h of operation). Thus, the air filters are of large-size, heavy, time-consuming, and costly to maintain.
A common air pollutant sucked in by motor vehicle engines is mineral dust (road dust) carried from the area by moving vehicles or by wind. Basic road dust components constitute SiO2 silica and Al2O3 corundum, whose dust share reaches 95%; and Fe2O3, MgO, and CaO. Furthermore, the dust consists of K2O, Na2O, and SO3 [1,2]. The physico-chemical dust composition strictly depends on the composition and type of substrate, climatic factors (winds, rains, snow, frost, droughts, etc.), and sedimentation of industrial dust, and dust from forest fires and volcanic ash [3].
Road dust, like most naturally occurring dusts, is polydisperse. Dust grains with a diameter of dp < 2 µm remain in the air for a very long time (their mass share in the total amount of dust is usually small). Dust grains with dp = 2–10 µm stay in the air for a long time and are sucked in by the engines. Dust grains with dp = 10–50 µm constitute a significant share in the total air mass sucked in by the engine when the vehicle is operated in conditions of high dust concentration in the air. Dust grains larger than 50 µm are found in the air in training yards, military training areas, quarries, construction sites, and mines [4].
Dust concentration in the air around a moving vehicle is variable and depends on the type of area, traffic, vehicle movement conditions, meteorological conditions, aerodynamic properties of the vehicle, type of wheels and steering system, type of surrounding soils, and height.
Typical concentrations can range from 0.01 mg/m3 in clean rural environments and up to about 20 g/m3 in heavy desert vehicle convoys [5,6]. According to [7], dust concentration in dusty environments is in the range of 0.001–10 g/m3. The authors of [8] investigated that the maximum dustiness within a few meters from all-terrain vehicles driving on a sandy road is 0.05–10 g/m3. The authors of [9] state that the dust concentration on highways ranges from 0.0004–0.1 g/m3 to 0.03–8 g/m3 when a column of vehicles is moving on sandy terrain. Similar dust concentration values in the air occur during take-off or landing of a helicopter at an accidental landing site. The rotating propeller lifts large amounts of mineral dust from the ground, which the engine absorbs with the air. The highest dust concentration values are found low (0.5 m) above the ground. At the height of the propeller end of the CH-53 helicopter taking off, the dust concentration may have the value s = 3.33 g/m3, and at a distance of 30 m, s = 2.11 g/m3 [10]. However, the actual dust concentration at the inlet to the engine intake system does not often exceed the value of 2.5 g/m3 [5,7,8]. Dust air concentration within the limits of 0.6–0.7 g/m3 significantly reduces visibility, and at a concentration of about 1.5 g/m3, it is practically zero [11].
Dust suspended in the air is highly abrasive and is the most common cause of accelerated wear of two frictionally matching parts, such as the P-PR-CW engine components (piston-piston rings-cylinder walls). According to the hardness on the basis of the ten-level Mohs scale, where 10 corresponds to diamond hardness, silica has a hardness of 7, and corundum of 9. Dust sucked in with the air goes over the piston and the aforementioned constitutes a decisive factor that the upper part of the cylinder and piston and the upper piston rings wear the most. Abrasive engine component wear is caused mainly by particles of 1–40 µm, with the most harmful particles occurring within the range of 1–20 µm [2,7,10,12]. The authors of [12] state that approximately 30% of pollutants entering the engine may be released, essentially in the unchanged form, along with exhaust gases from the cylinders to the exhaust system, thus increasing the emission of particulate matter (PM) from the engine. Only 10–20% of the dust enters the engine with the air through the intake system sediments on cylinder liner walls. Together with the oil, this dust part forms a kind of the abrasive paste which, in contact with engine mating surfaces, for example P-PR-CW, causes abrasive wear. The remaining part of the dust is evaporated or oxidised during the fuel combustion process in the cylinder (under high pressure and high temperature conditions).
Some mineral dust particles, the melting point of which is much lower than the peak temperature in the combustion chamber (2000–2500 °C), melt and then, along with the exhaust gases, enter the exhaust system and sediment on the cooler surfaces with which they come into contact (for example, in a catalytic converter) forming a glassy surface. This results from the fact that each of the polymorphic forms of SiO2 quartz (tridymite and cristobalite) is stable in the appropriate temperature range. At 1600 °C, the quartz melts to give quartz glass [13].
The most dangerous for the two working components are dust particles featuring a diameter dp equal to the oil film thickness hmin between two surfaces at a concerned moment. In typical combustion engine combinations, oil film thicknesses vary within the range of hmin = 0.3–10 µm [14].
The excessive wear of P-PR-CW causes increased blow-by of fresh load into the crankcase, resulting in compression pressure and engine power drop [15], increased oil aging, and increased oil flowing consumption over the piston, which increases solid particles emission [16,17,18]. One way to reduce friction and wear of the P-PR-CW connection is to apply a chromium nitride coating to the cylinder liner surface [19].
Adequate air purity at the inlet to traction motors is ensured by air filters, which differ in their principle and mode of operation, design, separation partition type, and operating efficiency. Trucks, special vehicles (tanks, armoured personnel carriers, infantry fighting vehicles), and working machines engines operating in high dust concentrations conditions in air are equipped with two-stage filters.
The idea of using two-stage filters is to pre-separate in an inertial filter (multicyclone, monocyclone) dust grains of the largest size and mass (dp > 15–35 µm). Then, small-sized dust grains remain in the air stream, which are directed to the filter cartridge, usually made of pleated paper with a suitably selected surface, where they are retained with an accuracy of greater than dp = 2–5 µm. As a result, air is supplied to the engine cylinders with the required purity and filtration system operation time is extended; thus, the service interval is limited by reaching a certain value of the pressure drop determined for a given engine for the permissible pressure drop Δpfdop.
The first filtration stage is then an inertial multi-cyclone filter, and the second a porous barrier in the form of a cylindrical pleated paper insert, which is the dominant material that filters air sucked in by internal combustion engines. Filter papers are characterized by high efficiency (over 99.5%) and accuracy—over 5 µm, but low water absorption (in the range of 200–250 g/m2), which forces designers to construct filters with a large paper surface [12,20,21]. Modern air filters use filter cartridges made of filter paper or pleated nonwoven fabric (Figure 1), due to the fact that the largest possible surface of the paper is obtained in a limited space in the engine compartment.
Non-woven materials and metal mesh are also used to filter engine intake air. In many military vehicles, which are operated in high dust concentration conditions in air, irregular, oil-soaked filter beds made of thin metal wire are still used as second filtration [22]. Plastic filter materials soaked in oil are still used in many agricultural tractors. Since the early 1990s, it has been popular in most Asian countries to further improve filtration efficiency and to lubricate pleated disposable filter elements for intake air in motor vehicles [23]. It has been proven that such a procedure ensures greater efficiency and dust absorption as well as sufficiently long service intervals. In addition, the oil is both chemically stable and harmless to the engine itself, and suppresses solid particles re-entrainment. It was found that oil coating on flat fibrous filters resulted in increased dust holding capacity, reduced pressure drop over time and reduced efficiency at lower filtration speeds. The authors of [23,24] reported that flat fiber filters’ filtration efficiency is largely affected by the amount of oil. It was found that higher amounts of oil resulted in higher particle filtration efficiency, but at the cost of a higher pressure drop, which is not desired.
During engine operation, as a result of the deposition and accumulation of dust particles on the filter cartridge, filter pressure drop Δpf increases (Figure 2). The intensity of filter pressure drop increase depends on the conditions in which the vehicle is operated, mainly on dust concentration in air and vehicle (engine) operating time. Excessively contaminated air filter elements cause difficult air flow to the engine, which reduces engine cylinders from filling [25,26]. This directly affects conditions of fuel combustion process, which may deteriorate engine performance (power decrease and fuel consumption increase) or cause an unfavorable change in exhaust gas composition [27,28]. Therefore, manufacturers, depending on the engine type, set permissible pressure drop value Δpfdop, which should not be exceeded by the air filter while the vehicle is in use. For high power diesel engines, these values are in the range of 3.7–7.5 kPa [26]. For SI engines, Δpfdop values are within the range of 3.8–5 kPa [29], while in [30], it was stated that maximum pressure drop values limiting the air supply for diesel engines with a turbocharger usually assume values of 5–7, 6 kPa. For special vehicles engines, the admissible pressure drop values are in the range of 9–12 kPa [31]. From a technical point of view, air filter service life is commonly defined as the restriction level, which causes pressure drop in passenger car filter by about 2.5 kPa above the pressure drop of new (clean) filter [32]. For trucks and special vehicles, it is assumed that Δpfdop values are approximately 6.25–7.5 kPa above pressure drop of the clean air filter [33].
Truck and special vehicles driving units are self-ignition engines of high power and large (1000–6000 m3/h) volumetric air demand. Thus, air filters are large, heavy, time-consuming, and costly to maintain. These vehicles are operated in conditions of high and variable air pollutants concentrations, and thus pressure drop increase varies with intensity. Therefore, servicing the filter (performed after specified vehicle mileage standard) may be premature, which increases operating costs, or too late. In the second case, the increase in the filter pressure drop caused by the continuous accumulation of dust on the air filter element and resulting from the high dust concentration in the air, is more intense, which results in exceeding the set value of Δpfdop (Figure 1). Further operation of the filter is possible, but the rapidly increasing pressure drop causes an additional engine power decrease and a fuel consumption increase [33,34,35]. Moreover, the high negative pressure behind the filter and vibrations of the filter element caused by cyclically sucked air pulsations may also cause mechanical filter element destruction (Figure 3).
Pressure drop increase in the compression ignition engine inlet system by 1 kPa causes an average engine power decrease by 0.4–0.6% and fuel consumption increase by about 0.3–0.5% [4].
The challenge for scientists is not only to achieve a lower pressure drop and a longer filter life while removing particles with maximum efficiency, but also to know and then predict changes in the air filter pressure drop, depending on vehicle mileage operating in conditions of high and variable dust concentrations, which will allow servicing of the air filter to be performed once Δpfdop is reached.
The most reliable research method in this area are experimental tests in real operating conditions with the possibility of ongoing filter pressure drop registration. The high cost and difficulties of conducting such research are obvious. For this reason, and due to the large number of parameters determining air filter flow efficiency and resistance, as well as large computers computing capabilities, numerical tests in the form of CFD programs are widely used. Numerical testing eliminates the need for time-consuming and costly experiments during the cycle. For air filters, it is possible to forecast changes in flow resistance, depending on the vehicle mileage.

2. Literature Review in the Field of Research Changes in Filter Materials Pressure Drop

Forecasting is now recognized as a key feature of machinery and equipment operating strategies, giving operators an important decision-making tool by quantifying the amount of time remaining until functionality is lost. The remaining machinery service life can be estimated using the four main condition forecasting techniques. These are experience-driven, model-driven, data-driven, and hybrid activities. Experience-based (or reliability-based) actions use data and knowledge gained during the machine’s past service life. Forecasting is used comprehensively, for example to forecast pollutants in the atmosphere, and in particular PM 2.5 concentration, as an effective way to help in management of measures important for air quality [36,37]; forecasting monthly water streams flow for effective planning and water resources management [38]; forecasting rain-fed agricultural production in dry and semi-dry areas [39]; end point of soot filtration process in diesel particulate filters for diesel engines, using for this purpose changes in pressure drop caused by soot accumulation [40]; long-term electric vehicles charging load [41]; engine oil level using artificial neural networks, which can help classify the state of oil degradation with high accuracy and plan replacement time [42]; frequency of changing the tractor engine oil based on the properties and oil concentration of contaminants [43], which will help to prevent engine failures and time at the charging station for each individual charging process of electric vehicle batteries [44]; heat demand for individual buildings, which should enable demand management and lead to predictive optimization of heat demand at the city level [45].
In the available literature, it is hard to find materials describing research on forecasting air filter pressure drop changes, despite the importance of the task it has in protecting the engine against premature wear.
There is, however, a significant number of works presenting the results of air filters optimization in the direction of flow resistance minimizing and maintaining high filtration efficiency at the same time. This direction of research is caused by the continuous technological development of car engines, which are equipped with newer systems and devices, as a result of which the space available for the air filter is reduced. As a result, the filtration surface is reduced and the filtration speed increases. For filter paper used in car air filters, the maximum filtration speed should not exceed the permissible value (0.08–0.12 m/s) [20,46]. Exceeding this value may result in filtration efficiency decrease due to the phenomenon of reflection and particles re-entrainment [47]. The basic method of ensuring the required surface of the filter material is its pleating.
The quantity, height, and pleats width are the main filter cartridge dimensions that describe the pleating geometry and determine filtration paper surface area. They also have a significant influence on the flow field, both in respect of the air and dust particles, through and near the pleated surface of the filter material, which in turn affects the separation efficiency and pressure drop.
Pleated fibrous material separation efficiency and pressure drop of the intake air filters in vehicle engines have been researched by several researchers, both in numerical and experimental terms. It has been evidenced that each clean pleated filter made from the pleated material has been characterised by an optimal quantity of the pleats, with which it has shown the minimum pressure drop [48,49,50]. Within a course of the researches performed without dust presence, it has been determined that the pressure drop has initially decreased as the quantity of the pleats has grown and then it has increased.
Chen et al. [51] conducted a simultaneous pressure drop and separation efficiency optimisation by means of the pleat quality factor. They found that a value of the pleat quality factor qp initially increased as the quantity of the pleats has grown, and then decreased. The quality factor of the pleat qp was defined as the quotient of pleat height (vertical distance from the upper edge of the pleat to the base) and its spacing (the distance between pleat tops). Théron et al. [52] numerically and experimentally evaluated pleat geometrical characteristics influence on pleated fibrous filters performances in the depth filtration stage of sub-micronic aerosols. Wiegmann et al. [53] simulated the pleated air filter with different pleat shapes and filter media and their impact on the pressure drop. This paper [48] experimentally and numerically investigated the pleated air filter with different parameters leading to improve engine performance. A series of parameters were investigated including: pleat height, pleat spacing, pleat shape, filter medium thickness, air velocity, engine speed, and engine torque. The findings showed that lower pressure drop was obtained at the sine wave-shaped pleated air filter. Mahesh [54] numerically investigated the pressure drop characteristics of the air filter of a four-cylinder engine with a triangular and rectangular element.
Several researches have also been carried out to assess the separation efficiency of the pleated fibrous filters with dust presence. Rebai et al. [55] determined that the optimum quantity of the pleats resulting from the researches on filter pressure drop in clean state was significantly lower than the pressure drop based on dust mass captured in the filter. Fotovati et al. [56] researched pleated bed geometry influence on dust deposition in the pleats, as well as on the related pressure drop and separation efficiency. The latter conducted numerical researches maintaining a constant pleat height. Feng and Long [57] developed a CFD simulation on a macroscale to research an optimal design of the pleated air filter under dust loading conditions. It was found that optimum pleated density in the clean state can lead to a higher energy consumption during the dusting process. As has been observed, an optimal dust deposition resulting from pleats density depends on pleat height and is lower when the pleat height is greater.
Saleh and Tafreshi [58] numerically developed a pleated air filter with rectangular and triangular pleats in both depth and surface filtration regimes under dust loading. Pressure drop and filtration efficiency were examined too. To determine the flow field, numerical simulations were carried out on a pleated panel-type air filter with a pleat height of 26 mm, a pleat spacing of 4.5 mm, and a pleat angle of 2.5 degrees. A series of ISO 12103 A2 fine dust tests were carried out on flat cellulose filter material with different speeds in the range of 0.05–0.35 m/s [59,60], experimentally studied pleat geometry effects, and dust load on the pressure drop. The results indicated that the efficient filtration area was mainly influenced by the pleat ratio, so it should be kept below qp = 1.59.
Numerical models can also be used to forecast changes in filtration efficiency and filter pressure drop during operation. Available literature contains studies where models for predicting filter pressure drop were made for a specific project and relate to samples of filter material, filter cartridges sections, complete single-stage, and two-stage filters.
Paper [61] presents a numerical methodology related to pressure drop in an air filter with a pleated filter component prediction. Pleat geometry and inlet velocity have been determined to be the key optimisation parameters for the filter components. Optimal pleats pitch has been determined to obtain a minimum pressure drop and found to be highly dependent on pleat height. Nevertheless, the authors of [62] used a dimensionless semi-empirical model for the purposes of pleat geometry optimisation, based on vehicle filters experimental and numerical tests results. It was determined that for a given pleat height and a constant air stream, there is such a width between the pleats for which the pressure drop reaches the smallest value. Pressure drop [63] and separation efficiency in a filter featuring the triangular and rectangular pleats have been modelled depending on a number of parameters, including, inter alia, pleat geometry, flow properties, fibre diameter, filter porosity, and particle inertia. Advantage of the algebraic expressions developed in the paper is that they are easy to use when designing a pleated filter and can be used in the initial design stage.
The predictive model related to the characteristics of the non-woven filters [64] takes into account filter parameters (bed thickness, packing density, fibre diameter), operating conditions (separation velocity, temperature), and pollutants in nanoparticle form with sizes ranging between 0–700 nm. Numerical and laboratory tests at separation velocity υF = 0.038 m/s and, depending on dust mass deposited on the filter surface, within the range of 0–3 g/m2 have been conducted. Filter papers used in the vehicle filters achieve a dust load of 220–250 g/m2. There is no information in the paper on the possibility of using this model to predict fibrous bed pressure drop with coarse-grained dust, with a grain size considerably above 1 µm, as it is the case with air filters used in trucks and special vehicles. The possibility of considering dust when modelling the pressure drop for a PTFE polytetrafluoroethylene membrane and HEPA filters is presented in the paper [65]. Model test results have been verified experimentally. In an experiment in respect of dust loading on the PTFE, HEPA membranes were conducted with monodisperse SiO2 dust with regard to the various particle diameters ranging between 52–553 nm and polydisperse SiO2 dust. Experimental researches have evidenced that an intensity in pressure drop increase is greater with smaller particle diameter. The aforementioned is consistent with research results [66,67]. Furthermore, the smaller the particle diameter, the greater the efficiency obtained. Developed models can be applied to entire fibrous filter materials filtration process with polydisperse dust. The model presented in [68] was developed in order to predict efficiency and filter pressure drop, depending on the porosity and filter bed thickness as well as fibre diameter, and has included dust particles ranging between 0.1–1 µm. Model tests were carried out for various porosity filter (from 50% to 90%) and thickness: 8, 16, 32, 64, 128, and 512 µm. An equation that fits the pressure drop changes exponentially with the porosity and fibre diameter for the same filter bed thickness was developed. It has been determined that this equation can be used to predict the pressure drop at the filter.
Pressure drop [69] in filters with an additional layer of nanofibers applied to a standard filter bed was forecasted, depending on average fiber diameter (61–260 nm), thickness (0.33–2.17 μm), and filter bed packing density. It has been found that a pressure drop is proportional to the inlet velocity and filter thickness. Such inlet velocity values occur in passenger cars and trucks filter beds, which makes it possible to use results to predict the pressure drop in these filters. Regrettably, the model fails to predict the pressure drop depending on the dust load.
A two-layer filter [12] material made from microfibers and sub-microfibers was developed for use in vehicle engines air filters. Laboratory research of the composite in respect of the separation velocity υF = 0.11 m/s with fine-grained test dust (ISO 12103-1) with a concentration of 1 mg/m3 was conducted. The results were compared with the research results of other filter materials. When the pressure drop of the standard medium reached 2 kPa with a dust load of 64 g/m2, the pressure drop of the composite for the same dust mass was 38% lower than for the standard medium, featuring high strength. Test results of two trucks in a sandy terrain with high dustiness confirmed lower (about 45%) filter pressure drop with composite material than the filter with standard material.
Numerical and experimental researches [70] were conducted on two filter samples made from nanofibers (dw = 187 nm, dw = 283 nm) and microfibers (dw = 2.7 µm), respectively. Samples arranged individually as composites of two samples adjacent to each other and a kit of two samples positioned at a certain distance were subject to the research. Researches carried out with the use of sodium chloride particles with a size within the range of 40–500 nm evidenced a diversified, in terms of mileage and value, increase in samples pressure drop depending on their configuration.
As particle deposition on microfiber sample (M) increases, there is a slow exponential increase in the pressure drop until a linear trend is reached, the slope of which is the same as the total decrease in linear trend pressure in both samples (M) and (N).
Comparative tests [23] were carried out on standard dry filter papers and the same filter beds containing different oil amounts, corresponding to 0%, 40%, and 80% of dry filter weight. Pressure drop and flow resistance tests on filter paper with a thickness of about 0.6 mm and a grammage of 102 g/m2 were performed using fine dust at constant filtration speeds between 0.11 m/s, 0.17 m/s, 0.33 m/s. In the case of media without oil, pressure drop increases rapidly from the very beginning, while in the case of media with oil, the increase intensity is low over a long period of dust accumulation.
The presented analysis shows that filter materials pressure drop depends and can be modeled with high accuracy depending on filter material parameters, which significantly simplifies the process of designing air filters. In real operating conditions, pressure drop depends mainly on dust mass retained on the filter bed and the aerosol flow conditions. The conducted analysis showed that it is possible to model the flow resistance with dust particles participation at various air flow velocities, similar to conditions occurring in internal combustion engines filter beds—about 0.1 m/s.
Depending on filter configuration and separation conditions, changes in pressure drop are forecasted with a linear or exponential, concave or convex trend, corresponding to the observed changes in the air pressure drop.
However, dust grain sizes adopted for modeling (below 1 µm, and sometimes around 10 µm) and the dust concentration values (about 1 mg/m3) significantly differ from the conditions prevailing during special vehicles operation—concentration up to 2.5 g/m3, dust with polydisperse particle sizes below 100 µm. Depending on the filter configuration and filtration conditions, changes in the flow resistance are forecasted with a linear or exponential, concave or convex trend, according to the observed changes in the air flow resistance.
Understanding the changes in the air filter resistance depending on vehicle mileage operated under various operating conditions (high and variable dust concentrations in the air) is possible only during experimental tests in the field. Such studies are presented in this paper.
Using an appropriate methodology, pressure drops in air filters of tracked vehicles were recorded for specific mileages and under various operating conditions. Forms of selected regression models best suited to the actual changes in pressure drops (flow resistance) of tracked vehicles air filters during their operation were determined. The conducted research shows that sufficiently good results of experimental data approximation and forecasting were obtained for the linear model. The model can be used to forecast changes in air filter pressure drop and to determine the service date (vehicle mileage) when the permissible resistance is reached, which will avoid unnecessary servicing and reduce vehicle operating costs. In addition, tracked vehicle air filters are big and heavy, which makes handling difficult.

3. Operational Air Filters Research

3.1. Object and Research Subject

The objective of the operational researches was to determine changes in air filter pressure drop as a function of vehicle mileage.
Pressure drop measure Δpf of the air filter is a difference in the static pressure before (p1) and behind (p2) filter at a set Q value:
Δpf = p1p2
The pressure drop presented as a function of the stream of the air volume Q is called the aerodynamic characteristic Δpf = f(Q).
The flow resistance presented as a function of the air volume flow Q is called the aerodynamic characteristic Δpf = f(Q) and increases parabolically with increasing Q. The air flow through the filter Q is equal to the air demand of the Qsil engine. For a given type of diesel engine and under certain conditions of its operation (for example, at a stabilized thermal state), the stream value Q flowing through the filter depends mainly on the rotational speed of the engine n.
Since the aerodynamic air filter characteristics depend on the stream Q, which is equal to the air demand QSil of the engine, air filter resistance being part of the engine intake system can be represented as a function of the engine rotational speed Δpf = f(n). With engine rotational speed increase, the air stream Q increases its value and at the rotational speed of the engine’s maximum power nN it assumes the nominal value QN equal to engine air demand QSilmax. This value is taken as the starting parameter for air filter laboratory tests’ characteristics and diagnostic measurements.
During vehicles use, due to filter insert contamination with dust, the air filter flow resistance Δpf increases its value. Diagnostic assessment of partition air filter technical condition consists of measuring its flow resistance Δpf at the rotational speed of the maximum engine power nN, which corresponds to the nominal air stream value QN = QSilmax.
Measured, after certain vehicle runs si, at the same engine speed n = nN, filter flow resistance Δpf assumes greater values:
Δpf0 < Δpf1 < Δpf2 < ... < Δpfi < ... ≤ Δpfdop
where Δpf0 is initial air filter flow resistance, Δpfi air filter flow resistance after a specified run si, Δpfdop allowable resistance.
Performed measurements are the basis for identifying filter technical condition and making a decision about the need to service it, at the moment of reaching or exceeding the permissible resistance value specified for a given filter and the engine cooperating with it. The research was conducted in summer (July–August) on military training ground with a rolling sandy terrain. The vehicles were operated in variable weather conditions (sunny days alternating with rainfall), which determined dust concentration in air sucked in by the engines. Unfortunately, the conditions for conducting the experiment, i.e., testing real vehicles, without equipping them with the additional, non-factory fitted equipment and moving over a large area, excluded the possibility of continuous monitoring environmental parameters; for example, temporary dustiness in the vicinity of the air intake vents in these vehicles.
Two-stage air filters (Figure 4) working in the “multi-cyclone-porous partition” system and ensuring clean air intake to the engine, which is the power unit of the military caterpillar vehicles, constituted the subject of the research.
A multicyclone is made of several dozen return cyclones with a tangential inlet, while the second stage of filtration is a moistened filter bed in the form of three (bottom—B, middle—M, top—T) cassettes arranged in a series filled with an irregular, compressed bed of steel wire with a diameter of ds = 0.25 mm. Filter bed thickness of the individual cassettes had the following values: 70 mm, 55 mm, 55 mm, respectively. Filtration speed in the cassettes, with the maximum air stream flowing through the filter, Qmax = 1800 m3/h, had a value of υF = 6 m/s. Such a large value is due to the small filtering cassettes area from the limited space available for the air filter inside the closed drive compartment.
The dust separated by cyclones falls to the bottom of the common settling tank, from where it is removed by ejection to the engine exhaust pipe, and then outside. The remaining dust part is transferred to filter cassettes, which must be serviced on a periodical basis. The servicing primarily consists of washing the filter cassettes (pouring a washing agent under pressure through them), drying them, and then submerging them (in order to moisten the wire surface) with engine or diesel oil heated to 60–100 °C. The presented technology of servicing air filters “multicyclone—moistened mesh bed” is labour-intensive, dangerous, and expensive.
The servicing should only be performed when filter cartridges loading with dust generates a pressure drop that results in excessive engine power loss. Due to the variable vehicle operating conditions, and mainly the variable value of dust concentration depending on terrain type, weather conditions, or the nature of vehicle use, it is difficult to predict a vehicle mileage at which such a moment will occur. For this reason, the air supply systems for many vehicles engines are equipped with sensors signalling that the filter has reached the admissible pressure drop Δpf = Δpfdop. This is the signal to perform the servicing. A lack of such signalling forces the filters to be operated in the planned-preventive system, i.e., after the mileage specified by the manufacturer. Due to different vehicle operating conditions, air filters’ maintenance is subject to maintenance, the current resistance of which slightly exceeds the initial filter resistance Δpf0 and is significantly lower than the admissible resistance Δpfdop, and the pressure drop of which has increased several times in relation to the initial one and is comparable and in extreme cases even higher than the admissible one.
Air filter pressure drop measurements were carried out on an everyday basis upon vehicle operation end, ensuring a stable engine thermal condition and measurements repeatability. Filter pressure drop measurements Δpf were determined at the maximum engine vehicle rotational velocity nN = 2000 rpm with an accuracy of 50 rpm, by repeating the measurement three times. These measurements’ average was used to analyse the results. The engine velocity was determined based on the readings of the on-board speedometer of the vehicle.
Pressure drop measurements were made with the use of an instrument, the basic component of which is a vacuum gauge with a measuring range of 0–20 kPa and an indication accuracy of 0.1 kPa.

3.2. Air Filters Research Results Analysis

The research included 22 filters (in 22 caterpillar vehicles) that had been previously serviced. For such prepared air filters, their initial pressure drops Δpf0 were determined. Considerable differences in the values of these initial pressure drops were found within 5–8 kPa (Figure 5).
The higher values of Δpf0 (after servicing—washing the filter cassettes) significantly exceeding the initial pressure drop of the brand-new filters (which is within a narrow range of Δpf0 = 4.1–4.2 kPa) were recorded for the filters operated in the vehicles for several years. This is due to the difficulty of completely washing away the dust from the densely packed and thick filter bed cassettes and the formation of hard-to-dissolve deposits on the bed wires of dust and oil, hardly removable by the washing technology applied (assumed as the instructions). With air filter operation time expiration, the mass of the unremoved sediments increases, and thus the initial pressure drop of the filter becomes higher and higher values frequently exceed twice the value of the pressure drop of the new filter; for example, in the case of filter No. 12—Δpf0 = 8.23 kPa.
The high initial pressure drop virtually reduces the mileage reserve of the vehicle conditional upon the technical efficiency of the air filter, i.e., difference between the admissible pressure drop Δpf0 and the initial pressure drop Δpf0. For this type of filter and the mating engine of the caterpillar vehicle, the admissible pressure drop conditioned by a 5% decrease in engine power had value of Δpf0 = 12 kPa [35].
The number of kilometres travelled each day was different and resulted from the tasks performed. The vehicles in the observed research period achieved various mileages, not exceeding 1000 km. The research results are presented in Figure 6.
Among the analysed mileages, there are the ones whose shape for the entire curve line is visually the same—the nature of the increase in the value of the pressure drop is similar, whether in relation to the same values or to the changes in these values. Regrettably, for most of the mileages subject to the analysis, the value ranges with the clearly different natures are visible. Both ranges, with linear (with a different slope) and non-linear (sometimes concave, sometimes convex, with the inflection points) shapes of these mileages, can be distinguished. This may indicate radically different conditions of vehicle use between the measurements, i.e., driving with the different levels of air pollution: on grass, on sandy or rocky off-road, on dry or wet dirt road, platooning as the first vehicle (low air polluted) or as a subsequent vehicle (high dustiness in the air). The possibility of using the researched vehicles in such diametrically different conditions basically excludes drawing general conclusions about the intensity of pollution of the filtration barrier of the air filters, even for a single specimen. It is possible to assess this intensity and build a given model of the changes in pressure drop as a function of the vehicle mileage for the different filters used only under the same (possibly similar, but not significantly different) conditions. Each significant change in these conditions should result in the abandonment of the current model and the development of a new one. Information on the weather and road conditions between the measurements was unfortunately not recorded (which has resulted from research conditions), so visually noticeable changes in the value of the pressure drop of the filter cannot be related to the change in conditions of use of the vehicle. When modelling the course of these changes, a constancy of these conditions was arbitrarily assumed, and the models were built on the basis of all available experimental values, without dividing them into sets, giving visually the same nature of changes in the pressure drop.
It is easy to notice (Figure 6) that the mileages for which the measurements were made are very diverse, from less than 200 to almost 1000 km. The mileage values between measurements (from just a few to nearly 200 km) and measurements quantity for a single vehicle (from 9 for a total mileage of 257 km) to 17 (for a total mileage of 609 and 702 km) were also very irregular. This results in a highly differentiated density of performance, for example: for the 176 km mileage, 9 measurements were made, and for the approximately 5.5 times longer distance, 985 km—14 measurements, an increase by only 5.
A difference in increase intensity in filter pressure drop value is also seen. The obtained mileages Δp = f(s) are not parallel (although it is noticeable for a large part of the mileages), but they intersect. These statements lead to the conclusion that these filters cannot be considered the same, which in turn makes it necessary to build individual models for each of them and excludes the possibility of building a generalised model for the filters of a concerned type. The conclusion about the impossibility of building a general model for the filter type is also confirmed by the curve line of the pressure drop value of the researched filters, shown in Figure 7 and Figure 8. In both of these figures, the mileages create the bunches of the curves arranged in a fan pattern. The observed changes in the value of the pressure drop of the researched filters are different, both in absolute and relative terms.

4. General Comments on the Usefulness of Modelling

Each technical product, including a machine, is manufactured for a specific purpose. This means that during the operation period, it should fulfil certain functions and meet the requirements defined by its designer or manufacturer. The scope and type of these requirements are specified in machine technical documentation in the form of basic operating parameters. In vehicles, these can be, for instance, engine power, maximum velocity, and fuel consumption. A set of values of such economic and machine technical parameters, determined during its construction and manufacturing for a given measure value of operation or within its specified range, determines machine technical condition. In the simplest case, two-state is assumed and states of fitness and unfitness for use are distinguished. They are unambiguously defined by the requirements contained in the technical documentation, and fitness is generally defined as the state in which the machine is capable of performing its functions according to the parameters specified in the technical documentation. The transitions between fitness and unfitness for use are called failures and repairs.
Technical diagnostics involves a problem of recognising the current, past, and future machine technical condition, while a diagnostic process itself is based on the results of the measurement values of the specific operating parameters, called the diagnostic parameters. It is worth recalling here that one of the properties required for a diagnostic parameter is its uniqueness. A threshold value of a diagnostic parameter dividing a range of the possible values of this parameter into ranges appropriate for the state of fitness and unfitness for use is called a limit value of this parameter. An output of the recorded values of the diagnostic parameter beyond the area appropriate for the condition of fitness, i.e., exceeding its limit value, corresponds to the machine entering the state of unfitness, i.e., its damage.
While repeating the measurements of the diagnostic parameter y of the machine for the successively larger values of the measure s of its operation (calendar time, working time, mileage, work cycles, fuel used, work done, etc.) and archiving their results, you have a certain set of data—the pairs of the values (si, yi), for which a regression model can be built with the aj parameters:
y = f(s, aj)
The forms of function f are usually elementary functions known from mathematical analysis (different functions are adequate to the experimental course to a different degree); due to the requirement of unambiguousness of the diagnostic parameter, monotonic functions are used here. The values of the aj parameters are determined using data approximation methods, the most popular of which seems to be the method of least squares (the given form of the approximating function, i.e., the values of the aj coefficients, is maximally adequate to the experimental course in terms of the approximation method used). The practical use of the obtained regression model is its inversion and determination of operation measure value, for which this model reaches any value within the range of its values; for example, the aforementioned limit value ygr. The knowledge of operation measure value for which, according to the adopted regression model, the diagnostic parameter will reach its limit value
sgr = f−1(ygr, aj)
is extremely important for operation system management, when deciding on the allocation of a given machine to perform a specific work; if the sgr is greater than machine value operation measure expected for assigned task completion (current sum value of the operation measure and operation measure value corresponding to the assigned task), it is highly probable that the machine will be subject to failure in the course of this task, resulting in the necessity to replace it with another machine and the necessity of its emergency transport to a repair service. Use of this probability concept is related to the adopted limit value of the diagnostic parameter, usually determined by some probabilistic method.
Assessment of modelling changes possibility in selected parameter value of the partition condition (metal, entangled wire) in vehicle air filter (i.e., the pressure drop of the filtered air reflecting a degree of its polluted with the captured, deposited matter) as a function of vehicle mileage (i.e., building the so-called regression models) was carried out based on data obtained for the 22 off-road vehicles designed for operation in various (difficult) operating conditions, including use in the areas with high dustiness. The vehicles were equipped with the same design filters. Due to the researched filters’ impossibility of self-cleaning, non-periodic increasing functions were selected for the modelling.
The method of least squares was applied to build the regression models, describing the obtained experimental mileage. In this sense, the obtained models are maximally adequate. The calculation algorithm contained in the MS Excel spreadsheet was used for a linear function, i.e., with two parameters—a slope factor and a shift factor. Therefore, it was assumed that all models will have functions with two parameters, while the data modelled with the non-linear functions required preliminary recalculation that enables functions’ use. When modelling the normalised mileages, whether absolutely or relatively, the initial value Δp0 is a natural third parameter.
Due to the obtained experimental mileages (visually: linear, concave and convex), for the assumed limitation of the function parameters to two, the following functions were used as modelling functions.
1. Linear function:
y = as + b
The values yi used for modelling purposes were: Δpi, Δpi − Δp0 and Δpip0. Parameter values a and b were obtained with the use of MS Excel methods. The model values are determined from:
Δpi = a1si + b1
Δpi = Δp0 + a2si + b2
Δpi = Δp0(a3si + b3)
according to the modelled values.
2. Power function:
y = bsa
For s = 0, function value y = 0. The values yi used for modelling were Δpi − Δp0. It helped easily avoid the problem resulting from a required zero value of the model for s = 0. Then, the model described by the relationship
Δp = Δp0 + bsa
assumes the initial value Δp0. This model use requires its linearization to the form:
ln(Δp − Δp0) = ln(bsa) = lnb + alns = B + aS = Y
where a—slope of the obtained straight line, B—shift factor of this straight line, Y and S—values obtained from:
Y = ln (Δp − Δp0)
S = lns
It is a double logarithmic data representation. Model linearisation excludes the use of the starting point for modelling, because ln0 does not exist. A model shift factor after its linearisation B is related to the corresponding factor of the non-linear model b by an obvious relationship:
B = lnb.
Parameter values a and B are obtained according to MS Excel procedures. The value b is determined from:
b = eB
The values of the model are determined from:
Δpi = Δp0 + bsia
3. Exponential function 1a:
y = bas = beslna = beAs = elnbeAs = eBeAs = eB + As
For s = 0, function value y = b. The values yi used for modelling were Δpi or Δpip0. Model use requires its linearisation to the form:
lny = ln(bas) = lnb + slna = B + As = Y
where A—slope of the obtained straight-line, B—shift factor of this straight line, Y—values obtained from the relationship:
Y = lny
This is a logarithmic data representation. Model linearisation excludes the use of the starting point for modelling, because ln0 does not exist. Model slope and shift factor after linearisation of A and B are related to the corresponding factors of the non-linear model a and b with obvious relationships:
A = lna
B = lnb
The values of the parameters A and B are obtained according to MS Excel procedures. The values a and b are determined from the relationships:
a = eA
b = eB
The values of the model are determined from the relationship:
Δpi = b1a1si
Δpi = Δp0b2a2si
according to the modelled values.
4. Exponential function 1b:
y + 1 = bas = beAs = eB + As
Values yi used for modelling purposes had values of Δpi − Δp0. Since the initial data value is zero, it is required to add any positive number to all data (number 1 has been assumed in the research). Such data shift operation value has enabled to obtain a non-zero model value for s = 0. This model use requires its linearisation to the form:
ln(y + 1) = ln(bas) = lnb + slna = B + As = Y
where A—slope of the obtained straight-line B—shift factor of this straight line, Y—values obtained from:
Y = ln(y + 1)
It is a logarithmic shifted values representation. Model linearisation does not exclude starting point use for modelling. The slope and model shift factor after its linearisation A and B are related to the corresponding factors of the non-linear model a and b by the relationships (17). Parameter values A and B are obtained with the use of MS Excel procedures. Values a and b are appropriately determined from the relationship (18). Model values are determined from:
Δpi = Δp0 + basi − 1
5. Exponential function 2a:
y = ebs^a
For s = 0, function value is y = 1. The values yi used for the purposes of modelling were Δpip0, because only they equal 1 for s = 0. This model use requires its linearisation to the form:
ln(lny) = ln(bsa) = lnb + alns = B + aS = Y
where: a—slope of the obtained straight line, B—shift factor of this straight line, Y and S—values obtained from the relationship:
Y = ln(lny)
S = lns
It is a double-logarithmic logarithm representation of relatively corrected value. Model linearisation excludes starting point use for modelling, because ln (ln1) = ln0 does not exist. The shift factor of the model after its linearisation B is related to the corresponding factor of the non-linear model b by the relationship (11). Parameter values a and B are obtained by MS Excel procedures. The value b is determined from the relationship (12). Model values are determined from:
Δpi = Δp0ebsi^a
6. Exponential function 2b:
y + 1 = ebs^a
The values yi used for modelling purposes were Δpi − Δp0. Since the initial data value is zero, it is required to add a positive number to all data (number 1 has been assumed in the research). Such data shift operation value has enabled to obtain a non-zero model value for s = 0. This model use requires its linearisation to the form:
ln(ln(y + 1)) = ln(bsa) = lnb + alns = B + aS = Y
where a—slope of the obtained straight line, B—shift factor of this straight line, Y and S—values obtained from the relationship:
Y = ln(ln(y + 1)
S = lns
It is a double-logarithmic logarithm value representation corrected relatively and shifted by 1. Model linearisation excludes starting point use for modelling, because ln(ln(0 + 1)) = ln0 does not exist. Model shift factor after its linearisation B is related to the corresponding factor of the non-linear model b by the relationship (11). Parameter values a and B were obtained by MS Excel procedures. The value b is determined from the relationship (12). Model values are determined:
Δpi = Δp0 + ebsi^a − 1
7. Logarithmic function:
y = ln(bs)a = aln(bs) = a(lnb + lns) = alnb + alns = B + aS
The values yi used for the modelling purposes were Δpi, Δpi − Δp0 or az Δpip0. This model use requires its linearisation to the form:
y = alnb + alns = B + aS
where A—slope of the obtained straight-line B—shift factor of this straight line, S—values obtained from the relationship:
S = lns
This is a logarithmic argument representation. Model linearisation excludes starting point use for modelling, because ln0 does not exist. The shift factor of the model after its linearisation B is related to the factors of the non-linear model a and b by an obvious relationship:
B = alnb
Parameter values a and B are obtained by MS Excel procedures. The value b is determined from the relationship (12). Model values are determined from:
b = eB/a
Model values are determined from:
Δpi = a1ln(b1si)
Δpi = Δp0 + a2ln(b2si)
Δpi = Δp0a3ln(b3si),
according to the modelled values.
It is obvious that increasing exponential models are concave, except for those with a power index (called exponential 2 in this paper) for 0 < a <1 (it is taken a root from the argument) and b > 0 and for a, b < 0. For these exceptions, the models used feature an inflection point. In these two exceptions, the smaller the a is, the more explicit the initial convexity and the less explicit the end concavity. The logarithmic models are convex. The power models are concave, linear or convex, depending on index value (convex for less than 1, linear for 1 and concave for greater than 1).
Changes of examples (for filter No. 14 according to Figure 5, Figure 6, Figure 7 and Figure 8) in the factors occurring in these models for an increasing number of the modelled data (according to an occurrence of the subsequent recorded filter pressure drop) are given in Table 1.
In Figure 9 and Figure 10, changes in selected models shape (i.e., exponential 2a: y = ebs^a and linear y = as + b) an increasing measurement results number used for modelling. It is clearly visible that as a range of the modelled experimental data increases (obtaining the newer and newer measurement results), model forms are corrected. This correction value is not large (changes in the first significant digit only by 1, sometimes changes only in the fourth significant digit). The exception here are logarithmic models for which changes by several orders of the value were obtained. However, logarithmic function nature is too non-linear to be used to model changes in air filters pressure drop, which will be pointed out. Model forms obtained for the full analysed mileage in case of this example are as follows:
  • Δp = 0.004138s + 5.602588 ≈ 0.0041s + 5.60,
  • Δp = 5.65 + 0.002461s1.079358 ≈ 5.65 + 0.0025s1.08,
  • Δp = 5.636023 × 1.000634s ≈ 5.64 × 1.000634s,
  • Δp = 5.65 + 1.049925 × 1.002305s − 1 ≈ 4.65 + 1.05 × 1.0023s,
  • Δp = 5.65e0.000511s^1.033640 ≈ 5.65e0.000511s^1.03,
  • Δp = 5.65 + e0.004501s^0.900222 − 1 ≈ 4.65 + e0.0045s^0.9,
  • Δp = 0.580821ln (479.015s) ≈ 0.58ln(479s).
Model graphs against measurement result values are presented in Figure 11. On visual assessment basis of the presented curves, it can be concluded that in spite of the differences in the form of the functions describing them, all, except for the logarithmic model, model the experimental mileage quite well.
In numerical terms, a matching degree of the considered models to the experimental data was stated by value determination of the correlation coefficient between the experimental values and the values obtained on the basis of the designated models. Table 2 presents this coefficient values for the examples of experimental mileages characterised by a visually similar (constant) shape for the entire mileage curve, sequentially for an increasing modelled data number, starting from the datasets ending with measurement 5—out of the 22 analysed experimental mileages, only 7 were identified (Table 2 shows values for two of them). Table 3 presents coefficient values for two selected mileages with different (visually variable) mileage curves (with the intervals of the linear (with the different gradients) and non-linear (sometimes concave, sometimes convex, with the inflection points) mileage curves shapes.
The data contained in these tables confirms good modelling visual assessment of the experimental mileages and its similar quality—correlation coefficient obtained values are high (except for those specified for logarithmic models, they are not less than 0.9 and in a significant number of the cases they are close to one), the differences between them occur at the 2nd and 3rd, and even at the 4th and 5th fractional decimal places. For experimental mileages with visually constant shapes, model orders from the point of view of the degree of their matching to the experimental mileages is rather constant as modelled data scope increases, and for mileages with the variable shapes, the order changes depending on the changes in experimental mileage curves shapes.

5. Regression Models Use to Forecast Changes in Filter Pressure Drop and Technical Condition

As stated in the introduction, a practical sense of building regression models is diagnostic parameter values’ predicting possibility, and thus inferring about the future technical machine condition, and in case of having diagnostic parameter limit value, also about operation margin until the machine reaches failure state. For the purposes of assessing, predicting quality can be used the simplest, apparently obvious measure, which is the difference between the predicted value Δpi by a given model for the value of the operation measure si, for which the next research has been performed, and the recorded value Δpie in this research, which can be called the absolute prediction error Bi:
Bi = Δpi − Δpie.
The relative prediction error bi can also be determined:
bi = (Δpi − Δpie)/Δpie.
With such elements defined order in these differences, the obtained error sign value is significant. If this value is positive, it means that the concerned model overstates the predicted value in relation to the real one, and therefore it is too close to the fact that the failure occurs. If the value is negative, it is otherwise. Since the error values are determined on an ongoing basis, according to the occurrence of the successive measurement results (for a set of n measurements, the model form is determined, after performing another, n + 1, measurement for the operation measure sn + 1, this model value is determined, which is compared with the value measured, and a new model form is determined, from which the value for the operation measure sn + 2 is determined, etc.), then for the obtained set of error values, their average level can be determined as their modules’ arithmetic mean (where the determined errors feature the different signs, this average will not be zero). Table 4 shows absolute and relative error values (multiplied by 100%) for predicting the value of the selected air filter with a visually constant shape (see Table 1 and Figure 9, Figure 10 and Figure 11) based on two selected models: linear and exponential 2a, updated on an ongoing basis, with first values row corresponding to the models determined on the basis of only two recorded values and the prediction with regard to the third value.
As can be noticed, for the presented example, the prediction is quite good—the prediction error does not exceed 1% of the real value, except for predicting with a large prognosis horizon (nearly nine times longer than others). While comparing the obtained prediction error values to correlation coefficient value between the experimental and modelled values for this filter (the fourth in Table 2—Filter No. 14), it can be perceived that the better model correlation with the experimental data will not necessarily result in better predicting—a correlation concerns known data, and predictive data refers to the data not known yet, with the values depending on the unpredictable changes in the machine operating conditions. In the example presented, correlation coefficient values are greater for the exponential model 2a than for the linear one, but prediction average relative error based on the linear model is less than this level for the exponential model 2a.

6. Conclusions

Non-decreasing change course in air pressure drop value in vehicle filters, as a vehicle mileage function, imposes increasing monotonic functions’ selection for development of regression models in respect to these changes. Linear functions, power, exponential and logarithmic with two parameters were applied in the research, which enabled to use the least squares method, implemented in a popular and broadly available spreadsheet, in order to determine these parameters’ value.
All monotonic functions with parameters obtained by the use of least squares method as pressure drop regression models for vehicle filters resulted in large correlation coefficient values between the modelled (experimental) and model (calculated) values. The lowest correlation for the modelled experimental data was usually obtained for the logarithmic models.
Pressure drop diversification value of the researched filters forces the development of each model on an individual basis. It is obvious that these models must be modified for constantly obtained new modelled data values. Observed variability in the shape of experimental air pressure drop curve lines of some researched filters probably results from vehicles usage in considerably changing road conditions (i.e., in significantly different level conditions of polluted filtered air, with vehicle usage periods in clean air and vehicle usage periods in very dusty air). Conditional mood usage here results from a lack of unequivocal information about road conditions by vehicles between filter pressure drop measuring moments. Obtaining such information would require continuous road monitoring conditions (e.g., air dust in air intake vicinity), which was impossible during the experiment. In subsequent papers, developing regression models possibility on the basis of a certain number (necessary to determine) of the last measurements should be considered. Certainly, this approach is also not entirely satisfactory and will not provide models that perfectly reflect road conditions variability.
The correlation degree of linear and non-linear models with experimental data is similar. Similarly, error rates related to predicting the actual air filter pressure drop values are similar. Therefore, it seems reasonable to state that a linear function can be used to model this pressure drop (pressure drop change, as vehicle mileage function is linear). This conclusion simplifies modelling and predicting operations as it eliminates the need to linearise non-linear models (the number of mathematical operations necessary to define the model is lower). Resignation from the use of non-linear models is also substantiated by vehicle use in considerably changing road conditions. This, in turn, causes increased rate differences in air pressure drop value at the filter, related to its clogging by the filtered pollutants—for vehicle use in air that is gradually cleaner, pressure drop convex shape changes are obtained, and in air that is gradually more polluted, pressure drop concave shape changes are obtained. Therefore, for cyclical change in the degree of air purity over a longer usage period, linear model adoption as a certain averaging of the cyclically adequate concave and convex models seems to be substantiated.
According to the authors, dynamically variables’ determination, individual regression models in respect of individual vehicles—linear or non-linear ones—is of great importance, as it gives the opportunity to move from a planned and preventive strategy for air filters’ use in special vehicles used in variable dustiness conditions (often high, for example, in deserts) to strategies according to their actual technical condition. This will make the vehicle maintenance system more flexible. The presented mathematical apparatus itself is so simple that it can be used by any engineer supervising the vehicle operation process, employing commonly available computing environments, for example, spreadsheets that are part of office suites. It is also possible to develop simple applications intended for the mobile devices that carry out appropriate calculations in order to obtain information about the estimated vehicle mileage up to the value in which it is necessary to perform servicing (for example, cleaning) of a concerned air filter.

Author Contributions

Conceptualization, T.D.; methodology, T.D.; software, G.B.; validation, T.D. and G.B.; formal analysis, T.D. and G.B.; data curation, T.D. and G.B.; writing—original draft preparation, T.D.; writing—review and editing, T.D.; visualization, T.D. and G.B.; funding acquisition, T.D. All authors have read and agreed to the published version of the manuscript.

Funding

The article was written as a part of university research grant implementation No. UGB 881/2021.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained with in the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bojdo, N.; Filippone, A. Effect of desert particulate composition on helicopter engine degradation rate. In Proceedings of the 40th European Rotorcraft Forum, Southampton, UK, 2–5 September 2014. Conference Paper. [Google Scholar] [CrossRef]
  2. Smialek, J.L.; Archer, F.A.; Garlick, R.G. Turbine Airfoil Degradation in the Persian Gulf War. J. Miner. Met. Mater. Soc. (TMS) 1994, 46, 39–41. [Google Scholar] [CrossRef]
  3. Vogel, A.; Durant, A.J.; Cassiani, M.; Clarkson, R.J.; Slaby, M.; Diplas, S.; Krüger, K.; Stohl, A. Simulation of Volcanic Ash Ingestion Into a Large Aero Engine: Particle–Fan Interactions. ASME J. Turbomach. 2019, 141, 011010. [Google Scholar] [CrossRef] [Green Version]
  4. Dziubak, T. Intake Air Filtration for Internal Combustion Engines of Motor Vehicles; Military University of Technology: Warsaw, Poland, 2012. [Google Scholar]
  5. Barris, M.A. Total FiltrationTM: The Influence of Filter Selection on Engine Wear, Emissions, and Performance; SAE International Fall Fuels and Lubricants Meeting and Exhibition. SAE Technical Paper Series: 952557; SAE: Warrendale, PA, USA, 1995. [Google Scholar] [CrossRef]
  6. Schaeffer, J.W.; Olson, L.M. Air Filtration Media for Transportation Applications. Filtr. Sep. 1998, 35, 124–129. [Google Scholar]
  7. Jaroszczyk, T.; Pardue, B.A.; Heckel, S.P.; Kallsen, K.J. Engine air cleaner filtration performance—Theoretical and experimental background of testing. In Proceedings of the AFS Fourteenth Annual Technical Conference and Exposition, Tampa, FL, USA, 1 May 2001. [Google Scholar]
  8. Pinnick, R.G.; Fernandez, G.; Hinds, B.D.; Bruce, C.W.; Schaefer, K.W.; Pendelton, J.D. Dust Generated by Vehicular Traffic on Unpaved Roadways: Sizes and Infrared Extinction Characteristics. Aerosol. Sci. Technol. 1985, 4, 99–121. [Google Scholar] [CrossRef]
  9. Barbolini, M.; Di Pauli, F.; Traina, M. Simulation der luftfiltration zur auslegung von filterelementen. MTZ Mot. Z. 2014, 75, 52–57. [Google Scholar] [CrossRef]
  10. Bojdo, N. Rotorcraft Engine Air Particle Separation. A Thesis Submitted to the University of Manchester for the Degree of Doctor of Philosophy in the Faculty of Engineering and Physical Sciences. 2012. Available online: https://www.escholar.manchester.ac.uk/uk-ac-man-scw:183545 (accessed on 19 April 2021).
  11. Szczepankowski, A.; Szymczak, J.; Przysowa, R. The Effect of a Dusty Environment upon Performance and Operating Parameters of Aircraft Gas Turbine Engines. In Proceedings of the Specialists’ Meeting—Impact of Volcanic Ash Clouds on Military Operations NATO AVT-272-RSM-047, Vilnius, Lithuania, 17 May 2017. [Google Scholar] [CrossRef]
  12. Long, J.; Tang, M.; Sun, Z.; Liang, Y.; Hu, J. Dust Loading Performance of a Novel Submicro-Fiber Composite Filter Medium for Engine. Materials 2018, 11, 2038. [Google Scholar] [CrossRef] [Green Version]
  13. Bojdo, N.; Filippone, A. A Simple Model to Assess the Role of Dust Composition and Size on Deposition in Rotorcraft Engines. Aerospace 2019, 6, 44. [Google Scholar] [CrossRef] [Green Version]
  14. Wróblewski, P. Technology for Obtaining Asymmetries of Stereometric Shapes of the Sealing Rings Sliding Surfaces for Selected Anti-Wear Coatings; SAE Powertrains, Fuels & Lubricants Meeting, SAE Technical Paper 2020-01-2229; SAE International: Warrendale, PA, USA, 2020. [Google Scholar] [CrossRef]
  15. Koszałka, G.; Suchecki, A. Changes in blow-by and compression pressure of a diesel engine during a bench durability test. Combust. Engines 2013, 154, 34–39. [Google Scholar]
  16. Bastuck, T.; Böhnke, F.; Hoppe, S.; Mittler, R. Systemische Kolbenringauslegung zur Reduzierung von Partikelrohemissionen. MTZ Mot. Z. 2020, 81, 50–55. [Google Scholar] [CrossRef]
  17. Lensch-Franzen, C.; Gohl, M.; Scholl, P.; Paoloni, F. Einfluss der Flüchtigkeit von Schmierölen auf die Öl- und Partikelemissionen. MTZ Mot. Z. 2019, 80, 46–55. [Google Scholar] [CrossRef]
  18. Gunkel, M.; Frensch, M.; Robota, A.R. Gelhausen, Innermotorische Emissionsreduzierung Zusammenhang zwischen Partikelemissionen und Ölverbrauch. MTZ Mot. Z. 2018, 79, 46–51. [Google Scholar]
  19. Singh, S.K.; Chattopadhyaya, S.; Pramanik, A.; Kumar, S.; Basak, A.K.; Pandey, S.M.; Murtaza, Q.; Legutko, S.; Litak, G. Tribological properties of chromium nitride on the cylinder liner under the influence of high temperature. Materials 2020, 13, 4497. [Google Scholar] [CrossRef] [PubMed]
  20. Dziubak, T.; Dziubak, S.D. Experimental Study of Filtration Materials Used in the Car Air Intake. Materials 2020, 13, 3498. [Google Scholar] [CrossRef]
  21. Maddineni, A.K.; Das, D.; Damodaran, R.M. Numerical investigation of pressure and flow characteristics of pleated air filter system for automotive engine intake application. Sep. Purif. Technol. 2019, 212, 126–134. [Google Scholar] [CrossRef]
  22. Rybak, P. Operating loads of impulse nature acting on the special equipment of the combat vehicles. Eksploat. Niezawodn. Maint. Reliab. 2014, 16, 347–353. [Google Scholar]
  23. Müller, T.K.; Meyer, J.; Thébault, E.; Kasper, G. Impact of an oil coating on particle deposition and dust holding capacity of fibrous filters. Powder Technol. 2014, 253, 247–255. [Google Scholar] [CrossRef]
  24. Maddineni, A.K.; Das, D.; Damodaran, R.M. Inhibition of particle bounce and re-entrainment using oil-treated filter media for automotive engine intake air filtration. Powder Technol. 2017, 322, 369–377. [Google Scholar] [CrossRef]
  25. Toma, A.; Fileru, C.S. The restriction produced by the air filtration system versus the restriction produced by the air filter. In Proceedings of the MATEC Web Conference, Chisinau, Moldova, 31 May–2 June 2018; Volume 178. [Google Scholar]
  26. Fileru, I.F.; Toma, A. Research on the Air Filters’ Maintenance for Diesel Engines. Procedia Technol. 2016, 22, 969–975. [Google Scholar]
  27. Synák, F.; Kalašová, A.; Synák, J. Air Filter and Selected Vehicle Characteristics. Sustainability 2020, 12, 9326. [Google Scholar] [CrossRef]
  28. Komorska, I.; Wolczynski, Z. Fault diagnostics of air intake system of the internal combustion engine. Appl. Cond. Monit. 2018, 10, 91–100. [Google Scholar]
  29. Jaroszczyk, T.; Pardue, B.A.; Holm, C.E. Recent Advances in Engine Air Cleaners Design and Evaluation. J. KONES Intern. Combust. Engines 2004, 11, 259–275. [Google Scholar]
  30. Bugli, N.-J. Service Life Expectations and Filtration Performance of Engine Air Cleaners; SAE Technical Paper Series 2000-01-3317; SAE International: Warrendale, PA, USA, 2000. [Google Scholar] [CrossRef]
  31. Bugli, N.J.; Green, G.S. Performance and Benefits of Zero Maintenance Air Induction Systems. In Proceedings of the SAE World Congress, Detroit, MI, USA, 11–14 April 2005. [Google Scholar]
  32. Norman, K.; Huff, S.; West, B. Effect of intake air filter condition on vehicle fuel economy. U.S. Dep. Energy (DOE) Inf. Bridge. 2009. [Google Scholar] [CrossRef] [Green Version]
  33. Bugli, N.J. Automotive Engine Air Cleaners—Performance Trends; SAE Technical Paper Series 2001-01-1365; SAE International: Warrendale, PA, USA, 2001. [Google Scholar] [CrossRef]
  34. Boruta, G.; Dziubak, T. Off-road vehicle air filter flow resistance forecast model. Zagadnienia Eksploat. Maszyn. 2002, 37, 125–146. [Google Scholar]
  35. Dziubak, T.; Szwedkowicz, S. Operating properties of non-woven fabric panel filters for internal combustion engine inlet air in single and two-stage filtration systems. Eksploat. i Niezawodn. Maint. Reliab. 2015, 17, 519–527. [Google Scholar] [CrossRef]
  36. Liu, H.; Xu, Y.; Chen, C. Improved pollution forecasting hybrid algorithms based on the ensemble method. Appl. Math. Model. 2019, 73, 473–486. [Google Scholar] [CrossRef]
  37. Choi, S.W.; Kim, B.H.S. Applying PCA to Deep Learning Forecasting Models for Predicting PM2.5. Sustainability 2021, 13, 3726. [Google Scholar] [CrossRef]
  38. Li, H.; Huang, G.; Li, Y.; Sun, J.; Gao, P. A C-Vine Copula-Based Quantile Regression Method for Streamflow Forecasting in Xiangxi River Basin, China. Sustainability 2021, 13, 4627. [Google Scholar] [CrossRef]
  39. Rezapour, S.; Jooyandeh, E.; Ramezanzade, M.; Mostafaeipour, A.; Jahangiri, M.; Issakhov, A.; Chowdhury, S.; Techato, K. Forecasting Rainfed Agricultural Production in Arid and Semi-Arid Lands Using Learning Machine Methods: A Case Study. Sustainability 2021, 13, 4607. [Google Scholar] [CrossRef]
  40. Zhang, B.; Zuo, H.; Huang, Z.; Tan, J.; Zuo, Q. Endpoint forecast of different diesel-biodiesel soot filtration process in diesel particulate filters considering ash deposition. Fuel 2020, 272, 117678. [Google Scholar] [CrossRef]
  41. Zheng, Y.; Shaoa, Z.; Zhang, Y.; Jianb, L. A systematic methodology for mid-and-long term electric vehicle charging load forecasting: The case study of Shenzhen, China. Sustain. Cities Soc. 2020, 56, 102084. [Google Scholar] [CrossRef]
  42. Rodrigues, J.; Costa, I.; Farinha, J.T.; Mendes, M.; Margalho, L. Predicting motor oil condition using artificial neural networks and principal component analysis. Eksploat. i Niezawodn. Maint. Reliab. 2020, 22, 440–448. [Google Scholar] [CrossRef]
  43. Hönig, V.; Procházka, P.; Obergruber, M.; Kučerová, V.; Mejstřík, P.; Macků, J.; Bouček, J. Determination of Tractor Engine Oil Change Interval Based on Material Properties. Materials 2020, 13, 5403. [Google Scholar] [CrossRef] [PubMed]
  44. Huber, J.; Dann, D.; Weinhardt, C. Probabilistic forecasts of time and energy flexibility in battery electric vehicle charging. Appl. Energy 2020, 262, 114525. [Google Scholar] [CrossRef]
  45. Hietaharju, P.; Ruusunen, M.; Leiviskä, K. Enabling demand side management: Heat demand forecasting at city level. Materials 2019, 12, 202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Gervais, P.C.; Poussier, S.; Bardin-Monnier, N.; Karcher, G.; Thomas, D. Combination of Single-Photon Emission and X-Ray Computed Tomography to visualize aerosol deposition in pleated filter. Sep. Purif. Technol. 2014, 126, 52–61. [Google Scholar] [CrossRef]
  47. Maddineni, A.K.; Das, D.; Damodaran, R.M. Air-borne particle by fibrous filter media under collision effect: A CFD-based approach. Sep. Purif. Technol. 2018, 193, 1–10. [Google Scholar] [CrossRef]
  48. Allam, S.; Elsaid, A.M. Parametric study on vehicle fuel economy and optimization criteria of the pleated air filter designs to improve the performance of an I.C diesel engine: Experimental and CFD approaches. Sep. Purif. Technol. 2020, 241, 116680. [Google Scholar] [CrossRef]
  49. Feng, Z.; Pan, W.; Wang, Y.; Long, Z. Modeling filtration performance of pleated fibrous filters by Eulerian-Markov method. Powder Technol. 2018, 340, 502–510. [Google Scholar] [CrossRef]
  50. Rebai, M.; Prat, M.; Meireles, M.; Schmitz, P.; Baclet, R. A semi-analytical model for gas flow in pleated filters. Chem. Eng. Sci. 2010, 65, 2835–2846. [Google Scholar] [CrossRef]
  51. Chen, C.W.; Huang, S.H.; Chiang, C.M.; Hsiao, T.C.; Chen, C.C. Filter quality of pleated filter cartridges. Ann. Occup. Hyg. 2008, 52, 207–212. [Google Scholar]
  52. Théron, F.; Joubert, A.; Coq, L. Numerical and experimental investigations of the influence of the pleat geometry on the pressure drop and velocity field of a pleated fibrous filter. Sep. Purif. Technol. 2017, 182, 69–77. [Google Scholar] [CrossRef] [Green Version]
  53. Wiegmann, A.; Rief, S.; Kehrwald, D. Computational study of pressure drop dependence on pleat shape and filter media. In Proceedings of the International Conference & Exhibition for Filtration and Separation Technology (Filtech), Rhein-Main-Hallen Wiesbaden, Germany, 27 February–1 March 2007; Volume 1, pp. 79–86. [Google Scholar]
  54. Mahesh, J. Parametric study and CFD analysis of air filter. Asian J. Converg. Technol. (AJCT) 2019, 5, 1–9. [Google Scholar]
  55. Rebai, M.; Prat, M.; Meireles, M.; Schmitz, P.; Baclet, R. Clogging modeling in pleated filters for gas filtration. Chem. Eng. Res. Des. 2010, 88, 476–486. [Google Scholar] [CrossRef]
  56. Fotovati, S.; Tafreshi, H.V.; Pourdeyhimi, B. A macroscopic model for simulating pressure drop and collection efficiency of pleated filters over time. Sep. Purif. Technol. 2012, 98, 344–355. [Google Scholar] [CrossRef]
  57. Feng, Z.; Long, Z. Modeling unsteady filtration performance of pleated filter. Aerosol Sci. Technol. 2016, 50, 626–637. [Google Scholar] [CrossRef] [Green Version]
  58. Saleh, A.; Tafreshi, H. A simple semi-numerical model for designing pleated air filters under dust loading. Sep. Purif. Technol. 2014, 137, 94–108. [Google Scholar] [CrossRef]
  59. Maddineni, A.K.; Das, D.; Damodaran, R.M. Experimental and Numerical Study on Automotive Pleated Air Filters; SAE Technical Paper Series 2016-28-0100; SAE International: Warrendale, PA, USA, 2016. [Google Scholar] [CrossRef]
  60. Li, S.; Hu, S.; Xie, B.; Xin, J.; Wang, F.; Zhou, F. Influence of pleat geometry on the filtration and cleaning characteristics of filter media. Sep. Purif. Technol. 2019, 210, 38–47. [Google Scholar] [CrossRef]
  61. Maddineni, A.K.; Das, D.; Damodaran, R.M. Oil-treated pleated fibrous air filters for motor vehicle engine intake application. J. Automob. Eng. 2019, 234, 702–713. [Google Scholar] [CrossRef]
  62. Del Fabbro, L.; Laborde, J.C.; Merlin, P.; Ricciardi, L. Air Flows and Pressure Drop Modelling for Different Pleated Industrial Filters. Filtr. Sep. 2002, 39, 35–40. [Google Scholar] [CrossRef]
  63. Saleh, A.M.; Tafreshi, H.V.; Pourdeyhimi, B. An analytical approach to predict pressure drop and collection efficiency of dust-load pleated filters. Sep. Purif. Technol. 2016, 161, 80–87. [Google Scholar] [CrossRef]
  64. Thomas, D.; Pacault, S.; Charvet, A.; Bardin-Monnier, N.; Appert-Collin, J.C. Composite fibrous filters for nano-aerosol filtration: Pressure drop and efficiency model. Sep. Purif. Technol. 2019, 215, 557–564. [Google Scholar] [CrossRef]
  65. Zhang, W.; Deng, S.; Wang, Y.; Lin, Z. Modeling the surface filtration pressure drop of PTFE HEPA filter media for low load applications. Build. Environ. 2020, 177, 106905. [Google Scholar] [CrossRef]
  66. Dziubak, T. Performance characteristics of air intake pleated panel filters for internal combustion engines in a two-stage configuration. Aerosol Sci. Technol. 2018, 52, 1293–1307. [Google Scholar] [CrossRef]
  67. Thomas, D.; Penicot, P.; Contal, P.; Leclerc, D.; Vendel, J. Clogging of fibrous filters by solid aerosol particles Experimental and modelling study. Chem. Eng. Sci. 2001, 56, 3549–3561. [Google Scholar] [CrossRef]
  68. Wang, X.; Kim, K.; Lee, C.; Kim, J. Prediction of air filter efficiency and pressure drop in air filtration media using a stochastic simulation. Fibers Polym. 2008, 9, 34–38. [Google Scholar] [CrossRef]
  69. Bian, Y.; Zhang, L.; Chen, C. Experimental and modeling study of pressure drop across electrospun nanofiber air filters. Build. Environ. 2018, 142, 244–251. [Google Scholar] [CrossRef]
  70. Leung, W.W.F.; Hau, C.W.Y.; Choy, H.F. Microfiber-nanofiber composite filter for high-efficiency and low pressure drop under nano-aerosol loading. Sep. Purif. Technol. 2018, 206, 26–38. [Google Scholar] [CrossRef]
Figure 1. Passenger car air filter element: (a), pleated filter paper, (b) non-woven fabric, (c) filter element made of pleated filter paper used in trucks.
Figure 1. Passenger car air filter element: (a), pleated filter paper, (b) non-woven fabric, (c) filter element made of pleated filter paper used in trucks.
Separations 08 00071 g001
Figure 2. Changes in filtration efficiency and pressure drop of the partition air filter during operation: 1—low dust in the air, 2—heavy dust in the air.
Figure 2. Changes in filtration efficiency and pressure drop of the partition air filter during operation: 1—low dust in the air, 2—heavy dust in the air.
Separations 08 00071 g002
Figure 3. Effects of using air filter after exceeding permissible resistance: (a) excessively contaminated filter element. (b) insert destruction.
Figure 3. Effects of using air filter after exceeding permissible resistance: (a) excessively contaminated filter element. (b) insert destruction.
Separations 08 00071 g003
Figure 4. “Multicyclone—moistened filter bed” filtration system in tracked vehicles used for the research: (a) complete air filter (multicyclone indicated), (b) cyclone, (c) filter cassettes, (d) filter cassette cutout, QG—polluted air stream, QG—clean air stream, QG—air stream sucked in with dust.
Figure 4. “Multicyclone—moistened filter bed” filtration system in tracked vehicles used for the research: (a) complete air filter (multicyclone indicated), (b) cyclone, (c) filter cassettes, (d) filter cassette cutout, QG—polluted air stream, QG—clean air stream, QG—air stream sucked in with dust.
Separations 08 00071 g004
Figure 5. Pressure drop initial values Δpf0 of the researched air filters.
Figure 5. Pressure drop initial values Δpf0 of the researched air filters.
Separations 08 00071 g005
Figure 6. Pressure drop values Δpf of the researched air filters depending on the mileage s of the caterpillar vehicles equipped with these filters.
Figure 6. Pressure drop values Δpf of the researched air filters depending on the mileage s of the caterpillar vehicles equipped with these filters.
Separations 08 00071 g006
Figure 7. The obtained pressure drop values Δp normalised absolutely (differently) to their initial value Δp0 for researched vehicles mileage value s.
Figure 7. The obtained pressure drop values Δp normalised absolutely (differently) to their initial value Δp0 for researched vehicles mileage value s.
Separations 08 00071 g007
Figure 8. The obtained pressure drop values Δp u normalised relatively (quotient) to their initial value Δp0 for mileage value s of the researched vehicles.
Figure 8. The obtained pressure drop values Δp u normalised relatively (quotient) to their initial value Δp0 for mileage value s of the researched vehicles.
Separations 08 00071 g008
Figure 9. Obtained exponential model changes 2a in pressure drop value Δp of selected vehicle air filter (with a visually constant shape of these changes) for an increasing number of modelled data.
Figure 9. Obtained exponential model changes 2a in pressure drop value Δp of selected vehicle air filter (with a visually constant shape of these changes) for an increasing number of modelled data.
Separations 08 00071 g009
Figure 10. Obtained linear changes models in pressure drop value Δp of selected vehicle air filter (with a visually constant shape of these changes) for an increasing number of modelled data.
Figure 10. Obtained linear changes models in pressure drop value Δp of selected vehicle air filter (with a visually constant shape of these changes) for an increasing number of modelled data.
Separations 08 00071 g010
Figure 11. Obtained changes models in pressure drop value Δp of the air filter for the entire analysed selected vehicle mileage s.
Figure 11. Obtained changes models in pressure drop value Δp of the air filter for the entire analysed selected vehicle mileage s.
Separations 08 00071 g011
Table 1. Obtained change models in pressure drop value Δp in selected vehicle air filters as mileage value function s.
Table 1. Obtained change models in pressure drop value Δp in selected vehicle air filters as mileage value function s.
Measurement Values ResultsModel(1)Δp= as + b(3)Δp= bas(4)Δp=Δp0 +bas− 1
nos [km]Δp [kPa]
105.65Measurementsabbaba
2175.711–20.0035295.6500005.6500001.0006221.0000001.003433
3435.771–30.0027515.6549895.6550631.0004821.0053641.002596
4755.911–40.0033775.6460185.6467611.0005840.9995811.003000
51026.011–50.0035315.6426385.6440641.0006060.9991711.003019
61406.161–60.0036575.6385435.6414091.0006211.0007771.002969
71606.251–70.0037565.6347195.6387991.0006331.0015191.002950
81856.291–80.0036405.6402075.6448211.0006101.0087281.002798
93607.111–90.0040345.6106285.6330811.0006381.0304051.002515
103807.21–100.0041175.604035.6310101.0006421.0362201.002443
114107.331–110.0041585.6004885.6308121.0006431.0411051.002389
124607.491–120.0041385.6025885.6360231.0006341.0499251.002305
Model(2)Δp=Δp0+ bsa(5)Δp=Δp0ebs^a(6)Δp=Δp0 +ebs^a− 1(7)Δp= aln(bs)
MeasurementsbababAab
2–30.0072290.7469360.0012930.7412900.0076450.7168460.0646561.33 × 1037
2–40.0037080.9632770.0006740.9522400.0042560.9066690.1275801.42 × 1018
2–50.0031701.0113770.0005820.9974090.0038010.9413270.1624299.01 × 1013
2–60.0029041.0368410.0005401.0193690.0036660.9518560.2049445.61 × 1010
2–70.0027461.0526520.0005151.0326090.0035990.9571050.2365031.29 × 109
2–80.0028311.0443240.0005341.0228470.0037850.9432500.2543962.28 × 108
2–90.0026091.0647990.0005141.0321240.0040450.9265520.40306243,424.3
2–100.0025121.0742320.0005071.0358590.0041950.9174700.4794544185.44
2–110.0024601.0794600.0005041.0370610.0043220.9101040.5343541172.35
2–120.0024611.0793580.0005111.0336400.0045010.9002220.580821479.0153
Table 2. Correlation values coefficient between the experimental and modeled pressure drop values Δp of selected air filters with a visually constant shape (for increasingly lower correlation coefficient values, bold black, black, blue, green, purple, brown and red colors were used).
Table 2. Correlation values coefficient between the experimental and modeled pressure drop values Δp of selected air filters with a visually constant shape (for increasingly lower correlation coefficient values, bold black, black, blue, green, purple, brown and red colors were used).
Experimental Values in The Coordinate System
s [km] − Δp [kPa]
Data RangeModel Correlation Coefficient Values (×0.001)
1234567
Separations 08 00071 i001to 5992998991988998998972
to 6995995995992996997944
to 7994997993986997998933
to 8991998989977999999918
to 9994998993984999999918
to 10991999989977999998923
Separations 08 00071 i002to 5996996996998996997942
to 6997998998999998999935
to 7998999998999999999935
to 8997997997995997996943
to 9997998999996999998864
to 10998999999997999999881
to 119989559999971000999894
to 12999999999995999997904
information contained in Table 1 and Table 4, and shown in
in Figure 6, Figure 7 and Figure 8 relates to this vehicle
Table 3. Coefficient values between the experimental and modeled pressure drop values Δp of selected air filters with a visually variable shape (for increasingly smaller correlation coefficient values, bold black, black, blue, green, purple, brown and red were used).
Table 3. Coefficient values between the experimental and modeled pressure drop values Δp of selected air filters with a visually variable shape (for increasingly smaller correlation coefficient values, bold black, black, blue, green, purple, brown and red were used).
Experimental Values in the Coordinate System
s [km] − Δp [kPa]
Data RangeModel Correlation Coefficient Values (×0.001)
1234567
Separations 08 00071 i003to 5975984974970983983949
to 6983979984986979981892
to 7991988992994988990893
to 8995994995994994995905
to 9991994988973992987914
to 10981987976958985976933
to 11981948976961985978944
to 12979987974959985977952
to 13974985968948982972957
to 14980969975958986979951
Separations 08 00071 i004to 5792993790783993993920
to 6863976862857976976860
to 7973979871866979980896
to 8915900916917900901761
to 9927877929936877879747
to 10947904949955904907782
to 11958910959966920923800
to 12965932967974933937816
to 13971943973978944947830
Table 4. Pressure drop values Δp of the selected air filter with a visually constant shape registered and predicted on the basis of two selected models.
Table 4. Pressure drop values Δp of the selected air filter with a visually constant shape registered and predicted on the basis of two selected models.
Measurement Results
Δpie [kPa]
Linear Model: Δp = as + bExponential Model 2a: Δp = Δp0ebs^a
Prediction
Δpi [kPa]
Error
Bi [kPa]
Error
bi [%]
Prediction
Δpi [kPa]
Error
Bi [kPa]
Error
bi [%]
5.775.8020.0320.551---
5.915.861−0.049−0.8245.832−0.078−1.316
6.015.99−0.020−0.3255.97−0.040−0.659
6.166.137−0.023−0.3746.123−0.037−0.594
6.256.224−0.026−0.4216.215−0.035−0.563
6.296.330.040.6296.3260.0360.568
7.116.951−0.159−2.2407.039−0.071−1.000
7.27.143−0.057−0.7867.158−0.042−0.585
7.337.292−0.038−0.5167.311−0.019−0.253
7.497.5130.0230.3127.5590.0690.927
average error level 0.04660.698 0.04740.718
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dziubak, T.; Boruta, G. Experimental and Theoretical Research on Pressure Drop Changes in a Two-Stage Air Filter Used in Tracked Vehicle Engine. Separations 2021, 8, 71. https://doi.org/10.3390/separations8060071

AMA Style

Dziubak T, Boruta G. Experimental and Theoretical Research on Pressure Drop Changes in a Two-Stage Air Filter Used in Tracked Vehicle Engine. Separations. 2021; 8(6):71. https://doi.org/10.3390/separations8060071

Chicago/Turabian Style

Dziubak, Tadeusz, and Grzegorz Boruta. 2021. "Experimental and Theoretical Research on Pressure Drop Changes in a Two-Stage Air Filter Used in Tracked Vehicle Engine" Separations 8, no. 6: 71. https://doi.org/10.3390/separations8060071

APA Style

Dziubak, T., & Boruta, G. (2021). Experimental and Theoretical Research on Pressure Drop Changes in a Two-Stage Air Filter Used in Tracked Vehicle Engine. Separations, 8(6), 71. https://doi.org/10.3390/separations8060071

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