Engine Vibration Data Increases Prognosis Accuracy on Emission Loads: A Novel Statistical Regressions Algorithm Approach for Vibration Analysis in Time Domain

.B.); Abstract: Statistical regression models have rarely been used for engine exhaust emission parameters. This paper presents a three-step statistical analysis algorithm, which shows increased prediction accuracy when using vibration and sound pressure data as a covariate variable in the exhaust emission prediction model. The ﬁrst step evaluates the best time domain statistic and the point of collection of engine data. The univariate linear regression model revealed that non-negative time domain statistics are the best predictors. Also, only one statistic evaluated in this study was a statistically signiﬁcant predictor for all 11 exhaust parameters. The ecological and energy parameters of the engine were analyzed by statistical analysis. The symmetry of the methods was applied in the analysis both in terms of fuel type and in terms of adjustable engine parameters. A three-step statistical analysis algorithm with symmetric statistical regression analysis was used. Fixed engine parameters were evaluated in the second algorithm step. ANOVA revealed that engine power was a strong predictor for fuel mass ﬂow, CO, CO 2 , NO x , THC, CO Sick , O 2 , air mass ﬂow, t exhaust , whereas type of fuel was only a predictor of t air and t fuel . Injection timing did not allow predicting any exhaust parameters. In the third step, the best ﬁxed engine parameter and the best time domain statistic was used as a model covariate in ANCOVA model. ANCOVA model showed increased prediction accuracy in all 11 exhausted emission parameters. Moreover, vibration covariate was found to increase model accuracy under higher engine power (12 kW and 20 kW) and using several types of fuels (HVO30, HVO50, SME30, and SME50). Vibration characteristics of diesel engines running on alternative fuels show reliable relationships with engine performance characteristics, including amounts and characteristics of exhaust emissions. Thus, the results received can be used to develop a reliable and inexpensive method to evaluate the impact of various alternative fuel blends on important parameters of diesel engines. and A.K.; writing—original draft preparation, T.Ž.; writing—review and editing, T.Ž., J.M., A.R., Á .B. and A.K.; visualization, T.Ž., J.M., A.R., Á .B. and A.K.; supervision, A.K.; project administration, A.K.; funding acquisition, A.K. All authors have read and agreed to the published version of the manuscript.


Introduction
Vehicle emission indicators have become increasingly more stringent in pursuit of environmental benefits [1]. Decarbonization programs aimed at the use of clean, low-carbon fuels in vehicles have been used to this end. Ambitious targets set to reduce concentrations of hazardous compounds in exhaust gases have recently been supplemented with strict Papers with predictions of engine exhaust emissions using ANN are very sparse. Cay et al. [30] showed that ANN with input data such as fuel type, engine speed, torque, fuel flow, carbon monoxide, unburned hydrocarbon, break specific fuel consumption and airfuel ratio can be very useful and accurate in CO, THC, BSFC and AFR values predictions. Unfortunately, engine vibration was not included as input data.
One more comprehensive study for modeling of performance, emission, and vibration of a compressed ignition engine using ANN technique was performed by Hosseini et al. (2020) [31]. Even 12 parameters in input layer were used in ANN model for predicting another 12 parameters in output layer. However, as in previous study engine vibration data was not included in the input layer.
Literature review can imply that building a prediction models for exhausted emission of ignition engines through statistical point of view is still challenging. Although, there are developed several models which can be used for exhausted emission parameters predicting but these models require a lot of input variables which can cause data gathering issues.
In this paper, three-step statistical analysis algorithm is presented to develop optimal prognostic model for 11 exhausted emission parameters. Vibration and sound pressure data in combination with stationary engine parameters such as engine power, fuel type and injection timing were used as prediction model input parameters.
The objectives of this study are the following: (1) to investigate the prognostic impact of vibrational and sound pressure data on engine emission parameters; (2) to evaluate if vibration and sound pressure data together with fixed engine parameters, i.e., power, type of fuel and fuel injection timing can improve prediction accuracy; (3) to investigate the best prediction model for exhausted emission parameters dependent from the fixed engine parameters and vibration and sound pressure data.
The rest of the paper is organized as follows: methodology and data description are described in Section 2; in Section 3, the results of the investigated statistical analysis for emission prediction are presented, and the conclusions close the article in Section 4.
During statistical analysis, we have analyzed one of the eleven engine emission parameter. In such way, we have symmetrically repeated statistical analysis for the rest of 10 engine emission parameters. So, we have proposed three-step statistical analysis algorithm with symmetrical statistical regression analysis. Also, LRM, ANOVA and AN-COVA requires normally distributed (i.e., symmetrical) quantitative data sample. We have checked this assumption with Kolmogorov-Smirnov criteria and found that data satisfies this assumption. Data normality confirmation is not the key to our analysis, so authors decided to not define this analysis.

Exhausted Emission Parameters
The investigations were elaborated on a naturally aspired, direct-injection 2.9 L IVECO AIFO Diesel engine provides the task of the machinery, which is driven by an M8B 160 generator. The engine set was made for generator purposes. Thus, it did not incorporate heat recovery. Therefore, a separate system had to be established for this purpose for the basic engine set. Exhaust gas analysis was performed using specialized measuring equipment, which was calibrated with a special gas before and after the tests. AVS 415 for FSN was used for particulate emission measurements. THC flame ionization detector, NOx chemiluminescence analyzer and CO (sick) and CO2 non-dispersive infrared detector and O2 para-magnetic were used to capture emissions of other gases. Table 1 summarizes instrumentation and their description, accuracy and measurement range used both for recording combustion characteristics and for emission tests, based on calculations the δ pi,t = 3.29%.

Type of Fuel
Four fuel blends and conventional diesel fuel was used in the research. Fuel blends consisted of conventional diesel fuel (D100) and hydrotreated vegetable oil (HVO100) or soybean oil methyl ester (SME100). Renewable biofuels accounted for 30% or 50% in the blends, and were marked as HVO30, HVO50, SME30 and SME50 in the article. All 4 biodiesel blends were volume-based. Physical and chemical properties of pure base fuels were analyzed in the laboratory and are presented in Table 2.

Experimental Engine
The fuel tests were carried out on an IVECO AIFO 2.9 liter direct injection, natu-rally aspirated compression ignition. The engine drives an M8B 160 generator, and the energy produced was used in a multi-power stage water heating boiler. The main pa-rameters of the engine are shown in Table 3. During the test series, using the pre-injection adjuster alone is used to set the pre-ignition angle, no other modifications were made to the engine. The engine operated in any case with the basic settings for fossil diesel.

The Engine Testing System
The test bench scheme is presented in Figure 1. The measurement system includes a number of static pressure and temperature sensors. Two piezo sensors are installed to measure the fast-changing pressures, one measuring the pressure in the combustion chamber and the other measuring the pressure in the high-pressure injection line. The sensors are connected to the data acquisition system after the charge amplifiers. The encoder on the crankshaft was used to trigger the measurements (1024 pulses/rev). 75 cycles were measured at each measurement point. The calibration of the system (pres-sure transducer, cables, charge amplifier and data acquisition system) before and after the tests (range 100 bar) and the accuracy was calculated.During the study. there were defined different starts of fuel injection timing, which were evaluated as crank angle degrees (CAD) before top dead center (BTDC) (5 CAD BTDC, 7.5 CAD BTDC, 10 CAD BTDC, 12.5 CAD BTDC, 15 CAD BTDC, and 17.5 CAD BTDC).  The combustion chamber pressures and injection tube pressures for D2 fuel and 20 kW load are shown in Figure 2. For both figures, the averages of 75 cycles at one operating point are shown. In the case of injection tube pressure, the pressure drop observed is associated with the first pressure drop at the start of injection and the pressure oscillations after injection are observed. The high ignition delay characteristic of this engine is clearly observed for the in-cylinder pressure. For pre-injections at 17.5, 15, and 12.5 degrees, the maximum pressure location is observed around 5 degrees (ATDC) for smaller pre-injections, the maximum location shifts significantly after TDC. The combustion chamber pressures and injection tube pressures for D2 fuel and 20 kW load are shown in Figure 2. For both figures, the averages of 75 cycles at one operating point are shown. In the case of injection tube pressure, the pressure drop observed is associated with the first pressure drop at the start of injection and the pressure oscillations after injection are observed. The high ignition delay characteristic of this engine is clearly observed for the in-cylinder pressure. For pre-injections at 17.5, 15, and 12.5 degrees, the maximum pressure location is observed around 5 degrees (ATDC) for smaller pre-injections, the maximum location shifts significantly after TDC.
3, x FOR PEER REVIEW

Vibrations and Sound Measurment System
Vibrations generated by internal combustion (IC) engine depend on disbalanced return motion and rotating parts, cyclic variation in gas pressure, dynamic excitation forces from rotating parts of the engine, and structural properties of the engine mounting mechanism. Stiffness and damping of the engine mounting structure must be high in the low frequency range, and low-in the high frequency range. Proper engine mounting must be used to reduce engine vibrations. Sometimes using mounting elements of appropriate characteristics at the engine-frame point of contact is necessary. Various types of vibration insulation materials are often used to reduce forces transmitted from the engine to the mounting structure.
When conducting experimental tests of the engine, sound of the surrounding environment and vibrations did not affect measurements in the room. The walls of the test room had an acoustic lining made of sound-insulating material, an acoustic door, and a soundproof double-glazed window for viewing inside the engine test chamber from the operator's room. A GRAS 46AE microphone (Frequency range: 3.15 Hz to 20 kHz; Dynamic range: 17 dB (A) to 138 dB; Sensitivity: 50 mV/Pa) was used to measure sound pressure (Figure 3a

Vibrations and Sound Measurment System
Vibrations generated by internal combustion (IC) engine depend on disbalanced return motion and rotating parts, cyclic variation in gas pressure, dynamic excitation forces from rotating parts of the engine, and structural properties of the engine mounting mechanism. Stiffness and damping of the engine mounting structure must be high in the low frequency range, and low-in the high frequency range. Proper engine mounting must be used to reduce engine vibrations. Sometimes using mounting elements of appropriate characteristics at the engine-frame point of contact is necessary. Various types of vibration insulation materials are often used to reduce forces transmitted from the engine to the mounting structure.
When conducting experimental tests of the engine, sound of the surrounding environment and vibrations did not affect measurements in the room. The walls of the test room had an acoustic lining made of sound-insulating material, an acoustic door, and a soundproof double-glazed window for viewing inside the engine test chamber from the operator's room. A GRAS 46AE microphone (Frequency range: 3.15 Hz to 20 kHz; Dynamic range: 17 dB (A) to 138 dB; Sensitivity: 50 mV/Pa) was used to measure sound pressure ( Figure 3a

Methodology of Statistical Analysis
Sixteen statistical parameters of time domain were evaluated before performing statistical regression analysis (Table 4). Descriptive analysis for engine exhaust parameters and time domain statistics were also done. Three-step statistical analysis was performed to develop optimal prognostic model for exhausted emissions.
In the first step the best time domain parameter and outcome data acquisition engine point was identified using univariate linear regression model (LRM):

Methodology of Statistical Analysis
Sixteen statistical parameters of time domain were evaluated before performing statistical regression analysis (Table 4). Descriptive analysis for engine exhaust parameters and time domain statistics were also done. Three-step statistical analysis was performed to develop optimal prognostic model for exhausted emissions.

Parameter Expression
Parameter Expression In the first step the best time domain parameter and outcome data acquisition engine point was identified using univariate linear regression model (LRM): where Y E -exhausted emission counts (dependent variable), E-type of exhausted emission, T k i ith time domain parameter estimate (independent variable) for kth outcome data acquisition engine point, α-intercept value, β-regression parameter for independent variable, ε-random error.
In the second step, one way analysis of variance (ANOVA) model was used to evaluate fuel, engine power and injection timing impact for exhausted emissions. Factor with the highest R 2 was determined as the strongest predictor for exhausted emissions and further was included into analysis of covariance (ANCOVA) model together with the best time domain parameter: where Y E -exhausted emission counts (dependent variable), E-type of exhausted emission, T k i ith time domain parameter estimate (covariate) for kth outcome data acquisition engine point, Z j -categorical variable with j levels (independent variable), α-regression intercept value, β-regression parameter for covariate, γ-regression parameter for independent variable, µ-regression parameter for covariate and independent variable interaction, ε-random error.
In the third step, accuracy between the best prognostic model and real data was evaluated using mean absolute percentage error (MAPE): where R E i -real exhausted emission counts for ith experiment,Ŷ E i -prognostic exhausted emission counts for i th experiment, E-type of exhausted emission.

Fuel Properties
A study of physical and chemical properties of fuel revealed that HVO has~2% higher and SME~13% lower heating value compared to conventional diesel fuel, thus the use of SME30 and SME50 increases fuel consumption. A significantly higher SME fuel kinematic and dynamic viscosity are likely to degrade injection quality of these fuel blends, prolong the time of evaporation and mixing with air, as well as the combustion time. These changes in SME fuel injection and combustion may lead to a longer rate of heat release as well as a pressure rise during the kinetic (premixed) combustion phase [32,33]. This increases the impact induced by rapidly increasing pressure on the piston and other surfaces of the combustion chamber. However, SMEs contain~11% oxygen which improves combustion quality, reduces CO, THC emissions and smoke content. Increased oxygen concentration in SME fuels and a more intense rate of heat release increases NOx emissions. First-generation biodiesel (including SME) concentration in fuel is still limited to 7% in Europe according to EN 590 diesel standard.
HVO is another type of fuel that may be used to prepare fuel blends. It has a lower density, but its kinematic and dynamic viscosity is close to that of diesel fuel. High HVO cetane number shortens the ignition delay phase of fuel blends, which in turn reduces the rate of heat release and pressure increase in the kinetic (premixed) combustion phase, reducing the mechanical load on the engine and NO x emissions upon a decrease in temperatures [34,35]. HVO fuels have a lower C/H ratio, which reduces smoke content. SME has a high cold filter plugging point (CFPP) and pour point temperature; therefore, this fuel cannot be used in winter. HVO CFPP and pour point temperatures are low, and these fuels can be mixed with diesel under winter conditions as well. HVO also has a high oxidative stability, good lubricity, thus can be used both pure and in blends with diesel in various proportions.

Vibration and Sound Pressure of the Engine
Determining vibrations and noise emitted by components of a diesel engine is one of the most difficult environmental tasks because each engine mechanism affects vibrations and noise separately.
For each experiment, vibration and sound pressure data were collected from the engine unit with a 3.2 kHz sampling frequency for 2 s. These results are presented in Appendix A.

Descriptive Statistics of Exhausted Emission Parameters
Descriptive statistical analysis revealed that engine power was strongly associated with almost all exhausted emission parameters. O 2 and mass flow of the air were inversely proportional to increased engine power. The average O 2 decreased up to 42.0% from 17.4 V% at 4 kW engine power to 10.1 V% at 20 kW engine power, while mass flow of the air had more slightly decrease (1.7%). Other parameters were directly proportional. Huge value jump regarding increased power was observed in NO x (289.5%), CO (221.5%) and CO 2 (200.0%) ( Table 5).
Any type of fuel had a major impact for CO and THC: The mean (SD) variated from 361 (233.4) ppm to 514 (292.7) ppm and from 102 (29.3) ppm to 139 (30.9) ppm for CO and THC, respectively. Fuels based on SME blend showed higher levels of fuel mass flow, CO, CO 2 , NO x , THC, t air , t fuel , t exhaust than HVO blends, while CO Sick and air mass flow were slightly higher in HVO based blends. O 2 value remained independent from fuel blends ( Table 5).
Increased injection timing was extremely associated with increased NO x . Mean (SD) of NO x value increased from 352 (194.5) ppm at 5 CAD BTDC to 1000 (495.0) ppm at 17.5 CAD BTDC for injection timing at 12 kW. Interestingly, THC achieved minimal values at 10 CAD BTDC with mean (SD) 111 (28.4) ppm, while moving towards to lower and higher injection timing category THC loads started to grow. The same tendency was observed in CO loads where minimal value with mean (SD) 423 (311.3) ppm was achieved at 12.5 CAD BTDC. Lower mean t exaust parameter values were associated with increased injection timing. Mean (SD) of t exaust value decreased from 320 (104.1) • C at 5 CAD BTDC to 301 (107.0) • C at 17.5 CAD BTDC. Other exhausted emission parameters remained stable at various injection timing at 12 kW (Table 5).

Step 1: Significant Vibro-Acoustic Parameters in Time Domain
A total of 16 time domain parameters for each 5 vibro-acoustic engine outcomes parameters (longitudinal 1 (LONG 1 ), longitudinal 2 (LONG 2 ), traverse 1 (TRAV 1 ), traverse 2 (TRAV 2 ) and sound pressure (SP) (Figure 3)) were evaluated as independent prognostic variables using univariate LRM. Regression analysis revealed that even 8-11 exhausted emission parameters could be prognosed using parameters T 2 -T 11 from all 5 vibro-acoustic engine outcomes. Also, 10 exhausted emission parameters could be prognosed using T 15 in LONG 1 and TRAV 2 point and T 16 in TRAV 1 and TRAV 2 point, however T 15 was weak prognostic parameter at LONG 2 and TRAV 1 engine points. Only one parameter T 11 calculated from TRAV 1 engine data acquisition point had statistically significant impact for all 11 engine exhausted emission parameters and was defined as the best prognostic parameter in univariate LRM (Table 6). Due to strong prognostic value only T 11 gathered from TRAV 1 was included in the further step 2 statistical ANCOVA analysis as model covariate variable. Descriptive analysis showed that engine power was strongly and directly associated with increased T 11 at TRAV 1 : the median (Q1-Q3) T 11 at TRAV 1

Step 2: ANCOVA Model for Engine Exhausted Emission Parameters
In this step one-way ANOVA analysis was performed to evaluate prognostic impact of engine power, type of fuel and injection timing for all 11 exhausted emission parameters. Analysis revealed that engine power had the highest R 2 values for fuel mass flow, CO, CO2, NOx, THC, COSick, O2, air mass flow, texhaust. Type of fuel had the highest R 2 value for tair and tfuel while fuel injection timing had very low R 2 values and was not included in It should be noted that variation was very low in power groups. Distribution of T 11 at TRAV 1 remained to be independent from various type of fuels. Also, the variation was increased in each of fuels compared with variation in engine power groups. It could be explained that engine power is very strong factor for the shape of T 11 at TRAV 1 distribution and have significant impact in fuel blends groups (Figure 4b).
Distributions T 11 at TRAV 1 in various injection timing groups were slightly different. Interestingly, there was no any linear tendency injection timing and T 11 at TRAV 1 values. The minimum T 11 at TRAV 1 was observed in 7.5 CAD BTDC group with median (Q1-Q3) 1.38 (1.37-1.43), while maximum was observed in 12.5 CAD BTDC group 1.48 (1.45-1.50) (Figure 4c). Figure 4 shows parameter T 11 distribution in engine power, type of fuel and injection timing categories. T 11 had moderate (Spearman's ρ = 0.688, p < 0.001) and weak (Spearman's ρ = 0.194, p = 0.036) correlation with engine power and injection angle, respectively, and there was no correlation with type of fuel.

Step 2: ANCOVA Model for Engine Exhausted Emission Parameters
In this step one-way ANOVA analysis was performed to evaluate prognostic impact of engine power, type of fuel and injection timing for all 11 exhausted emission parameters. Analysis revealed that engine power had the highest R 2 values for fuel mass flow, CO, CO 2 , NO x , THC, CO Sick , O 2 , air mass flow, t exhaust . Type of fuel had the highest R 2 value for t air and t fuel while fuel injection timing had very low R 2 values and was not included in further analysis (Table 7). Engine power as independent factor has very high predictive model fit for fuel mass flow, CO, CO 2 , O 2 and t exhaust with R 2 > 90 (Table 7). R 2 for T 11 at TRAV 1 was somewhat lower than engine power or type of fuel but remained still statistically significant for all exhausted emission parameters (Table 7).
Finally, from Table 7 the best predictive models for exhausted emission were developed of T 11 at TRAV 1 , engine power and type of fuel. Regarding ANCOVA Equation (2) (3) and (4). All R 2 values were higher in ANCOVA model than in ANOVA model and much higher than in univariate LRM. The major R 2 improvement between ANOVA and ANCOVA models was reached in NO x and t air with 12.4% and 9.9% increment, respectively. Exhaust parameters with R 2 > 90% in ANOVA model had very low (<1%) R 2 increment in ANCOVA model.

Step 3: Prediction Model Accuracy
Regression model with independent predictor T 11 at TRAV 1 alone remained the worst prediction model regarding relatively high MAPE values. On average, ANOVA and ANCOVA models were 3.3 times more accurate that univariate LRM and remained as strong prediction models for exhausted emissions. Tremendous prediction accuracy improvement switching from univariate LRM to in ANOVA and ANCOVA models was fixed in fuel mass flow, CO, CO 2 , O 2 and t exhaust parameters. However, THC, t air , and t fuel had lower accuracy improvement than other parameters but still ANOVA and ANCOVA models had showed significant accuracy improvement. Only prediction for air mass flow parameter had not been improved through all prediction models. (Figure 5).
Step 2 ANCOVA model with two factors (T 11 at TRAV 1 and engine power or type of fuel) remained the best predictive model for all 11 exhausted emission parameters with the lowest MAPE values for all exhausted emission parameters ( Figure 5). The most MAPE reduction was defined in NO x , t air , and t fuel comparing ANOVA and ANCOVA models. Figure 6a shows that increased prediction accuracy was observed in higher engine power regardless type of fuel, whereas vibration covariate has not had any influence in low (4 kW and 8 kW) engine power. Furthermore, vibration covariate has also not had any influence in pure diesel engines. However, T 11 at TRAV 1 influence greatly improved while diesel mixtures were used (Figure 6b,c).         Table 8 parameter estimates. One-way ANOVA prognosis was made form one-way ANOVA with only power or type of fuel predictor.

Conclusions
Experimental data showed that engine power has the biggest influence on generated vibrational and sound pressure. Root mean square (RMS) values increased 20-30% while engine power was changed from 4 kW to 20 kW (data not shown). Type of fuel and injection timing was not associated with significant RMS differences (data not shown).
Vibration and sound pressure data have shown a high predictive power for exhausted emissions in univariate LRM analysis. All non-negative time domain statistics were associated with larger number of exhausted emissions which can be predicted using only vibration and sound pressure data. Furthermore, best predictor was defined RMS related function T11. Vibration data gathered from TRAV1 and aggregated with T11 had statistically significant impact for all 11 exhausted emission parameters.
ANOVA analysis revealed that engine power and type of fuel can be used in exhausted emissions prognostic equation development. In our case engine power was strong predictor for fuel mass flow, CO, CO2, NOx, THC, COSick, O2, air mass flow, texhaust, whereas type of fuel was only a predictor of tair and tfuel. 17 Table 8 parameter estimates. One-way ANOVA prognosis was made form one-way ANOVA with only power or type of fuel predictor.

Conclusions
Experimental data showed that engine power has the biggest influence on generated vibrational and sound pressure. Root mean square (RMS) values increased 20-30% while engine power was changed from 4 kW to 20 kW (data not shown). Type of fuel and injection timing was not associated with significant RMS differences (data not shown).
Vibration and sound pressure data have shown a high predictive power for exhausted emissions in univariate LRM analysis. All non-negative time domain statistics were associated with larger number of exhausted emissions which can be predicted using only vibration and sound pressure data. Furthermore, best predictor was defined RMS related function T 11 . Vibration data gathered from TRAV 1 and aggregated with T 11 had statistically significant impact for all 11 exhausted emission parameters.
ANOVA analysis revealed that engine power and type of fuel can be used in exhausted emissions prognostic equation development. In our case engine power was strong predictor for fuel mass flow, CO, CO 2 , NO x , THC, CO Sick , O 2 , air mass flow, t exhaust , whereas type of fuel was only a predictor of t air and t fuel .
Consolidation of univariate LRM and ANOVA analyses shown increased prediction power for all 11 exhausted emission parameters. MAPE value reduced approximately three times from univariate LRM MAPE value. Prediction accuracy has grown in higher engine power modes by adding vibration data as a covariate variable in ANCOVA model while prediction remained independent from vibration data in lower engine power (4 kW and 8 kW). The same tendency was seen in different type of fuel. Prediction of exhausted emission parameters has grown in fuels HVO30, HVO50, SME30 and SME50 by adding vibration data as a covariate variable in ANCOVA model while prediction remained independent from vibration data in fuel D100.
Study findings conclude that adding vibration parameter to prognostic model helps to achieve higher prediction accuracy rate for various biodiesel fuels exhausted emissions in a regression model. The findings of the study allow concluding that adding vibration parameter to prognostic model helps to achieve higher prediction accuracy rate for various biodiesel fuel exhaust emissions in the regression model. Characteristics of noise and vibrations of diesel engines running on alternative fuels show reliable relationships with performance characteristics of the engine, volumes, and characteristics of emissions. Thus, the results received allow creating a reliable and inexpensive method for assessing the impact of various alternative fuel blends on important parameters of diesel engines.

Appendix A
Typical results of measurement of vibrations and sound pressure (D100, HVO50 and SME50) under 20 Nm load and injection timing 10 CAD BTDC, respectively, marking: red-1 point transverse direction; blue-1 point longitudinal direction; green-2 points transverse direction; orange-2 points longitudinal direction; orange-sound pressure, please see Figure A1.

Appendix A
Typical results of measurement of vibrations and sound pressure (D100, HVO50 and SME50) under 20 Nm load and injection timing 10 CAD BTDC, respectively, marking: red-1 point transverse direction; blue-1 point longitudinal direction; green-2 points transverse direction; orange-2 points longitudinal direction; orange-sound pressure, please see Figure A1.