Abstract

To increase the efficiency and accuracy of computing, an improved combined weight coefficient is used to develop an improved heat transfer model in AVL-BOOST environment. Similarly, a five-component biodiesel skeletal mechanism is employed to investigate the combustion process of biodiesel fuel. Then, the AVL-BOOST model is validated by the experimental results under different conditions. Finally, the improved heat transfer model is employed to investigate the propulsion and load characteristics of diesel engine fueled with biodiesel fuel in terms of power, BSFC, soot, and NOx emission. The result shows that the errors between experiment and simulation are less than 2% and the simulation model can predict the propulsion and load characteristics of the diesel engine. In addition, the comprehensive characteristic of case 5 is the best. Moreover, the big inject orifice is not beneficial to the fuel atomization and more soot is produced. Thus, it is very important to choose the appropriate injection rate reasonably.

1. Introduction

Due to the reliability and economy of diesel engine, the diesel engine has become the main power source of mechanical equipment [1], such as the ship, construction machinery [2], and heavy-duty truck [3]. The diesel engines have become the main power source for ships, construction machinery, heavy trucks, and other mechanical equipment, because of the reliability and economy of diesel engine. Particulate matter (PM) and NOx emissions from diesel engines pose a threat to the ecosystem and public health [4]. With the successive implementation of the “supplementary provisions” [5], the higher requirements have been put forward on how to reduce NOx and PM emissions and improve the performance of diesel engine [6]. To meet the requirements of the regulations, renewable energy sources should be developed to replace fossil fuels [7] and the system structure of existing diesel engines should be improved to adapt the needs of new fuels [8, 9].

As a new type of renewable energy, biodiesel fuel is a good alternative of petrochemical energy and has been highly concerned and expected by experts and scholars in recent years [10]. Biodiesel is not only a new type of renewable energy but also a good substitute of petrochemical energy, so it has attracted great attention and expectation in recent years. Compared with fossil fuels, biodiesel is nontoxic and biodegradable, and its emission is significantly reduced [11, 12]. Furthermore, unfortunately, some factors will result in quality degradation, such as light, metal contamination, and oxygen exposure, but these reductions can be mitigated by some technologies and additives [13]. Thus, lots of researchers had carried out the research on biodiesel fuel [14]. This is exactly what many researchers explore and research on biodiesel. Zhang et al. found that the high oxygen content of biodiesel improved combustion and greatly reduced HC and CO emissions [15]. The work carried out by Liu et al. showed that the addition of biodiesel caused the increase of BSFC and the reduction in torque. It is due to the shorter ignition delay and lower calorific of biodiesel [16]. In addition, the similar phenomena were found by Can [17] and Özener et al. [18]. Imdadul et al. had found that proper pentanol in biodiesel could reduce BSFC and improve the performance of biodiesel [19]. Further, the emission studies showed that the addition of biodiesel did significantly reduce CO and HC emissions [20], but NOx emissions increased significantly [21, 22]. Because the spray performance has a great impact on the combustion and emission characteristics of biodiesel, it is necessary to study the spray performance and improve the applicability and development of biodiesel [23].

Diesel engine combustion is a complex and changeable physical and chemical process [24]. Under the condition of limited resources in China, researchers often apply a numerical method to solve the scheme selection and design improvement, so as to improve the efficiency and reduce the design cost of diesel engine [25]. The designer only chooses the corresponding model and inputs the diesel engine parameters [26], then can make a preliminary demonstration of the model, because the modeling and simulation software provides a predesign platform [27]. For instance, Lino et al. had established a nonlinear injection system model for a common rail diesel engine [28]. They found that the model verified by experiment could predict accurately the spray process and the good spray performance improves the combustion process of diesel engine. The simulation model of an electronic control unit (ECU) pump injection system was established by Lähde et al., and the influence of various structural parameters on its economy and emission characteristics was studied [29]. The result showed that different fuel injection mass had a great impact on emission and performance characteristics of diesel engine. In addition, Wu et al. had developed a numerical model of diesel particulate filter and investigated on the effects of key factors on the performance of CeO2-based catalyzed diesel particulate filter. They found that the soot regeneration rate and the N2O concentration first increased and then decreased as the exhaust gas temperature increased and the higher exhaust gas temperature suppressed the production of N2O but raised the possibility of catalyst poisoning [30]. The research conclusions showed that the delay of the solenoid valve was an important factor. Frosina et al. had proposed a one-dimensional fuel injection system (FIS) model to explore the impact of biodiesel on the performance characteristic of common rail diesel engines [31]. The simulation results were consistent with the experimental results. There are many literatures about diesel engine emission and performance characteristics, which mainly focus on the fuel injection system model of diesel engine, but there are few researches on the mechanical injection system modification of electronic control injection system [32]. Therefore, the software provides a good idea for modeling diesel engine fuel injection system, and it can effectively analyze the fuel injection rate of diesel engine [33].

Typically, the researchers often used GT Power and ALV-BOOST software to investigate the effect of zero-dimensional or dual-zone combustion model on diesel engine performance [34] and emission characteristics and input different parameters to study its impact on overall parameters of diesel engine [35, 36]. For example, Ma et al. established a simulation model by using AVL-BOOST and neural network software to study the relative contribution of operating parameters on the emission and performance characteristics of common rail diesel engine in terms of torque, smoke, NOx, and braking specific fuel consumption (BSFC) [37]. Similarly, some scholars used the AVL-BOOST software [38] to establish the simulation model of the diesel engine and analyzed the influence of different fuel ratios on the comprehensive performance of the diesel engine [39, 40]. In addition, other scholars had carried out the modeling and made the analysis of four main factors which could impact the performance and emission characteristics of diesel engine, including the first and the second injection mass ratio [40], interval time [41], main injection start time [42], and exhaust gas recirculation (EGR) ratio [43]. However, the one-dimensional models are simpler, faster and included more detailed physical models [44]. Therefore, more advanced models should be adopted, including simulations of combustion, heat transfer and pollutant formation processes [45].

In this paper, an improved heat transfer model is established by using the combined weight coefficients to simulate the combustion and heat transfer process of diesel engine fueled with biodiesel in AVL-BOOST environment. In addition, the five-component biodiesel skeletal mechanism is used to predict the combustion process of biodiesel. The AVL-BOOST model is validated by the experimental results under different conditions. Moreover, the new model is employed to analyze the engine characteristic and investigate the effect of fuel injection rate on the engine characteristic of diesel engine fueled with biodiesel fuel.

2. Materials and Methods

2.1. Mathematical Model
2.1.1. AVL MCC Combustion Model

To improve the calculation accuracy, the paper uses the quasidimensional MCC-AVL combustion model to predict the combustion in the cylinder. The model considers the influence of diffusion combustion and premixing, and it can predict the emission characteristic of diesel engine accurately. The expression of heat release rate is expressed as follows: where is the combustion heat release rate, is the total heat release rate of diffusion combustion, is the total heat release rate of premixed combustion, and is the crank angle.

The total heat release rate of diffusion combustion is expressed as where is the total heat release rate of diffusion combustion, is the combustion constant, F is the evaporative fuel mass, is the low calorific value of fuel oil, is the turbulent energy density, is the effective air mass fraction, is the EGR influence constant, and is the mixing ratio constant.

The actual heat release rate of premixed combustion is expressed as where is the total heat release rate of premixed combustion, is the crank angle, is the premixed combustion duration, is the start of combustion angle, and is the shape parameter.

2.1.2. Heat Transfer Model

Heat transfer in a diesel engine cylinder is a very complex process. Woschni 1978 heat transfer model is generally used for calculation of high-pressure circulating heat transfer [32]. In the process of heat transfer, the influence of convective heat transfer in the cylinder is mainly considered [33]. The Woschni model published in 1978 for the high-pressure cycle is summarized as follows: where is the cylinder temperature, is the cylinder pressure, is the cylinder diameter, is the gas velocity coefficient, is the mean velocity of piston, is the model constant, is the inlet valve closing cylinder volume, is the actual cylinder volume, is the inlet valve closing cylinder pressure, and is the inverted cylinder pressure.

The heat transfer coefficient of the Hohenberg heat transfer model is summarized as follows: where is the cylinder pressure, is the mean velocity of piston, and is the actual cylinder volume.

One of the defects of the single heat transfer model is that the heat transfer in the diesel cylinder cannot be predicted accurately. In order to improve the calculation accuracy and make up for the shortage of a single model and based on the theory of optimized weighted array prediction, an optimized weighted array model of heat transfer in the diesel engine cylinder is established in this paper. The optimal weighting coefficients are derived by the minimum variance method. Thus, the advantages of the two models are combined. The combined model makes the model more accurate. The author uses the minimum variance method to find the optimal weighting coefficient in this paper. Therefore, the advantages of these two modes are combined. The combined model makes the model more accurate.

Assume that is the expected predictive value, and the experimental observation value is . The predicted value of the Woschni 1978 heat transfer model is ; the error between the model value and the test value is . The predicted value of the second model is ; the error between the second model value and test value is . Respectively, the expected predictive value of and error are as follows: where is the weight coefficient of the Woschni 1978 heat transfer model and is the weight coefficient of the Hohenberg heat transfer model. where , , , and are defined as , , and .

The minimum quadratic sum of the combined expected forecasting error is as follows: with and .

It can be proved that  ≤ min() and min() ≤ min(). It shows that the combined heat transfer model is better than the single model prediction.

2.2. Fuel Injection System

The instantaneous injection rate (IFR) and the fuel injection pressure can be obtained by the oil pump testing device (6PSDW300). The IFR system is composed of the electric control part and mechanical part. Moreover, the IFR was measured by the EFS8427 measurement units of French EFS Company. The main experiment steps can be expressed as follows:

Step 1. Preparation work. Check the oil pump testing device, fill the low-pressure oil tank with fuel, and install the electronic unit pump in the mechanical part. The electric control part is supplied by a single-phase current and connected to the computer through the Ethernet by a specific software as a user interface.

Step 2. Set the relevant operating parameters of the oil pump test bed. Start motor drive oil pump operation and transport oil from the low-pressure oil tank to the high-pressure common rail pipe. Fuel which is injected from the fuel injector will be collected in a constant static pressure fuel collection chamber.

The main parameters of FIS of the ECU pump is shown in Table 1. In the experimental process, the injection pressure and fuel supply pressure were measured by the pressure sensor installed in the mechanical part, and the original parameters such as the fuel injection rate and fuel injection mass were provided to the FIR control system.

2.3. Simulation Model of Entire Diesel Engine

In order to investigate the steady-state and transient-state performance characteristic of diesel engine, the entire model of diesel engine was developed in the AVL-BOOST environment. Figure 1 shows the layout of the entire diesel engine. In addition, the main parameters of diesel engine are shown in Table 2. The air cooled by the air cooler enters the cylinder through the air intake passage. After combustion process, the exhaust gases are discharged and flow into the turbine inlet. The turbine drives the turbocompressor, which delivers the compressed air to the cooler.

2.4. Fuel Sample Preparation

Rapeseed is widely grown in southern China, where nearly half the population uses the rapeseed oil as cooking oil. Therefore, the diesel engine fueled with the rapeseed methyl ester (RME) was studied experimentally. The RME can be obtained by the transesterification method, and the biodiesel can be widely used in diesel engines [16]. In addition, in the experimental reactor, the biodiesel transesterification was carried out by alkali catalysis about 1.2 h. Actually, the transesterification conducted 1 : 6 oil molar to methanol ratio with 1%wt/wt potassium hydroxide as the alkaline catalyst. The physical properties of biodiesel are shown in Table 3. In the early research work, our team has carried out the related research on biodiesel. Therefore, the detailed biodiesel information can be obtained from references [16, 21].

2.5. Model and Validation

In order to validate the simulation result, the experiment had been carried out on the experimental bench of diesel engine. Horiba MEXA-1600 was used for measuring generated NOx with 1% error, while AVL Dismoke-4000 was employed for measuring generated soot. FCMM-2 was used for measuring the BFSC. DEWE-2010CA was employed for monitoring the combustion of diesel engine. The list of measurements, the measuring range, and accuracy is shown in Table 4. In addition, an ECU control system was employed for controlling the electronically controlled diesel engine. EFS-IFR600 was employed for measuring fuel injection rate with 0.5% measuring error. A hydraulic dynamometer was used for measuring the diesel engine load. Moreover, the temperature, flow, and pressure were measured by suitable sensors. The schematics of experimental device are shown in Figure 2.

In the AVL-BOOST simulation environment, a five-component biodiesel skeletal mechanism was employed to investigate the combustion process of biodiesel fuel consisting of 134 species and 475 reactions [46], the friction is calculated with the Chen-Flynn model, and the combustion in the cylinder is fitted with the AVL MCC combustion model [47, 48]. The combustion model is effective to consider the influence of the jet kinetic energy and fuel injection rate on the instantaneous heat release and can accurately predict the formation of nitrogen oxides and soot. In addition, the combined heat transfer model, Woschni1978 heat transfer model, and Hohenberg heat transfer model are employed to simulate the heat transfer process in the cylinder, respectively.

In order to validate the combined model, the experiment was carried out and the comparisons of heat transfer coefficient are shown in Figure 3 at 100% load. It can be found that predicted result of the combined model is in agreement with the experimental results. Through contrast experiment and correction, the predicted result of the combined model has the maximum error of 2.1%, but the predicted result of the Woschni 1978 heat transfer model has the maximum error of 3.2%. Thus, the combined model can better predict the performance characteristic of diesel engine.

3. Results and Discussion

The performance and emission characteristics of diesel engine and the effect of fuel injection rate on diesel engine fueled with RME were investigated.

3.1. Load and Propulsion Characteristic of Diesel Engine

In the work, the diesel engine was performed at the E3 and D2 cycle modes in accordance with international ISO8178 standards. More specifically, the diesel engine was investigated at five different speeds (628 rpm, 710 rpm, 799 rpm, 911 rpm, and 1000 rpm), corresponding to 25%, 50%, 75%, and 100% engine load. In addition, the diesel engine was operated at 25%, 50%, 75%, and 100% engine loads when the engine speed is 1000 rpm. The E3 and D2 cycle modes for different test points are shown in Table 5. The improved heat transfer model and Woschni 1978 model are used to investigate the load and propulsion characteristics of diesel engine fueled with RME. The load and propulsion characteristics are discussed in the part in terms of torque, BSFC, temperature of behind the intercooler, NOx emission, and exhaust temperature.

3.1.1. Torque

The variation of torque with different models is shown in Figure 4. Figure 4(a) shows that the torque increases with the increase of power, torque is proportional with power, and all the points are on the same line. In addition, the result also shows that the torque can be calculated accurately by the numerical modes. More specifically, the maximum difference between the experimental result and predicted result calculated by the Woschni 1978 model is 3.8% and the maximum difference between the experimental result and predicted result calculated by the improved model reached 2.6%. In addition, Figure 4(b) shows that the torque increases with the increase of engine speed. The improved model has higher accuracy. More specifically, the maximum difference between the experimental result and predicted result with the Woschni 1978 model is 3.5% and the maximum difference between the experimental result and predicted result with the improved model reached 2.7%. Thus, the improved combination model can better predict the torque of diesel engine.

3.1.2. Brake Specific Fuel Consumption

The variation of BSFC with different models is shown in Figure 5. Figure 5(a) shows that the BSFC firstly decreases and then increases with the increase of the power. In addition, the BSFC reached the minimum value at 75% engine load. Because the operating point of 75% engine load was the design operating point, the BSFC reached the minimum value at 75% load. Actually, the predicted results of the improved model have better accuracy compared with those of the Woschni 1978 model. More specifically, the maximum difference of BSFC between the experiment and the improved model is 2.8%. In addition, the maximum difference between the experiment and Woschni 1978 model is 3.5%. Figure 5(b) shows that the BSFC firstly decreases and then increases with the increase of the engine speed. Between 799 rpm and 910 rpm, the BFSC is lower than the other operating points. The BSFC reaches the minimum value at 75% engine load. It is due to the fact that the operating point is the design operating point. Therefore, the efficiency is relatively large. In addition, it can also be found that the improved model can better predict the BSFC and the maximum difference between improved model and experiment is slightly lower than the maximum difference between the experiment and Woschni 1978 model.

3.1.3. NOx Emission

The variation of BSFC with different models is shown in Figure 6. Figure 6(a) shows that the NOx emission firstly decreases, then increases and finally decreases with the increased engine power. The NOx emission increases with the increasing power, but the increase of power is slower than the increase of NOx emission. Thus, the NOx emission per unit power increases. When the load is higher than 75%, the NOx emission deceases. It is due to the fact that combustion is deteriorated by the increased fuel. In addition, Figure 6(b) shows that the NOx emission firstly increases and then decreases with the increased engine speed. It is due to the higher cylinder temperature caused by fuel injection mass. When the speed is higher than 799 rpm, the effects of high cylinder temperature would be dominated by the impact of the decreased oxygen content. Thus, the NOx emission per unit power decreases. It also can be found that the improved combination model can better predict the NOx emission.

3.1.4. Exhaust Temperature

The variation of exhaust temperature with different models is shown in Figure 7. Figure 7(a) shows that the exhaust temperature increases with the increase of the engine power. It is due to the fact that the increased fuel injection mass results in the increase of engine load. Similarly, the maximum difference between the experiment and the improved model is slightly lower than the maximum difference between the experiment and Woschni 1978 model. In addition, Figure 7(b) shows that the exhaust temperature increases with the increase of the engine speed. The engine load increases with the increase of speed, then results in the increase of exhaust temperature. It can also be found that the improved model can better compute the temperature. Thus, the improved combination model can better predict the exhaust temperature.

3.1.5. Temperature of behind the Intercooler

The variation of temperature of behind the intercooler with different models is shown in Figure 8. Figure 8(a) shows that the temperature of behind the intercooler increases with the increase of the engine power. In addition, Figure 8(b) shows that the temperature of behind the intercooler increases with the increase of the engine speed. The main reason is that the increase of load leads to the increase of combustion temperature, which increases the cooling load. It can also be found that the improved model can better compute the temperature. More specifically, the maximum difference between the experiment and the improved model is slightly lower than the maximum difference between the experiment and Woschni 1978 model.

As mentioned, the Woschni 1978 model and improved model also can better predict the propulsion and load characteristics of diesel engine. However, the improved model is more accurate in prediction. It is due to the fact that the advantages of the two models are combined.

3.2. Fuel Inject Rate

In order to optimize and match the structural parameters of the inject system, the engine is carried out at 1000 rpm and the fuel injection mass remains constant in the experiment. Four parameters were considered mainly in this work. To further investigate deeply the detailed optimized planning, the experiments were carried out by changing the plunger diameter from 13 mm to 14.5 mm (4 levels), the injection orifice diameter from 0.22 mm to 0.28 mm (4 levels), the length of high-pressure oil pipe from 850 to 900 (2 levels), and the cam profile velocity from 0.43 mm/° to 0.46 mm/° (2 levels).

The measurement precision and calibration of fuel inject rate are very essential to analyze the formation of mixed gas and combustion processes in the engine cylinder. In order to improve the combustion, the high-pressure oil pipe, injector orifice, plunger diameter, injection duration angle, and injection pressure should be considered in the paper. The injection duration is less than 35°CA, and the injection pressure is greater than 100 MPa. According to the advanced plan, five cases were selected and were employed to investigate the effects of variable injection rates on the engine characteristic in the paper.

The five cases selected are shown in Table 6, and the fuel injection rates are shown in Figure 9. It can be found that the effect of the length of high-pressure oil pipe on injection rate is insignificant and the injection pressure decreases with the increase of length of the high-pressure oil pipe in Table 5. It is due to the increased resistances caused by the increased length of the high-pressure pipe. In addition, Figure 9 shows that the fuel injection rate and the injection pressure increase with the increased of the diameter of the plunger when the fuel injection mass keeps constant. The injection duration angle decreases with the increase of the diameter of the plunger. However, the injection pressure and injection duration angle decrease with the increase of the injector orifice diameter. Injector orifice diameter and diameter of plunger are significant to fuel injection rate.

The improved heat transfer model is used to investigate the effect of fuel injection rate on the engine characteristic of diesel engine fueled with RME in terms of cylinder pressure, heat release rate, soot emission, NOx emission, and cylinder temperature. In the experiment, the engine speed is 1000 rpm and the fuel injection mass remains constant.

As shown in Figure 10(a), the cylinder pressure of case 2 is the highest. It is due to the fact that the fuel atomization and cylinder combustion are improved by the high fuel injection pressure. In addition, the maximum peak pressure corresponding to the lower due to the delayed injection pressure and the time to reach the peak pressure also have a certain delay. As shown in Figure 10(b), as the fuel injection rate is larger, the heat release rate is larger at an early stage. It is due to the larger injector orifice diameter. The smaller the fuel injection pressure, the greater the fuel injection map, the longer the delay of heat release rate.

As shown in Figure 10(c), the soot generation rate of case 3 is the first to reach the peak. The injector orifice diameter is larger so that the injection pressure becomes lower. Finally, the fuel atomization becomes worse. Therefore, case 3 generates the most soot emission, followed by case 1 and case 4. As shown in Figure 10(d), the NOx generation rate of case 4 is the minimum and the NOx generation rate of case 3 is the maximum. The high cylinder temperature is the important parameter, resulting in the NOx generation. In the early stage, a large amount of fuel is injected into the cylinder in case 3. In addition, the fuel injection pressure is low so that the fuel atomization is relatively bad. It is easy to form the local high temperature area in case 3. Thus, the NOx generation rate of case 3 is the maximum. It can be found that the curve of NOx generation is delayed with the decrease of fuel injection rate and NOx generation value is in turn to reach the maximum value.

As shown in Figure 10(e), due to the high fuel injection pressure and fine atomization, the peak in-cylinder temperature of case 2 is the highest, followed by case 5. The faster the fuel pressure is built, the faster the temperature will reach the maximum value. In other words, with the delay of building pressure time, the peak in-cylinder temperature is also delayed. The fuel injection pressure and atomization effect of fuel have the most efficient influence on the cylinder temperature.

3.3. Experimental Verification

In the experimental process, the engine speed is 1000 rpm and the fuel injection mass remains constant. The engine characteristic of diesel engine fueled with RME is investigated in terms of power, BSFC, NOx emission, and soot emission. Figure 11 shows the comparisons of the simulation and experiment results. It can be found that the error between the experiment result and simulation result is less than 2%. Thus, the simulation model can accurately forecast the trend of the engine characteristic.

The BFSC and soot are the minimum, and the effective power of case 5 is the maximum. More specifically, compared with case 4, NOx emission was more than 0.94%, but the soot was reduced by 5%, and BFSC decreased by 0.42% in case 5. In addition, it can be found that the economy is worst in case 3. The NOx emission is the least in case 4. The efficiency power and torque are relatively large, and the formations of BFSC and soot are the least in case 5. Compared with case 3, the effective power and the effective torque are increased by 1.02% and 1.023%, respectively; the soot formation decreased by 7.4%, but the nitrogen oxides only increased by 1.1% in case 5. Based on an overall consideration of various factors, case 5 is the most reasonable. Therefore, it is very important to choose the injection rate reasonably. The big inject orifice is not beneficial to the fuel atomization, and more soot is produced.

4. Conclusion

Nowadays, the energy crisis [4956] and environmental pollution [5561] are more and more serious accompanied by global economic development and manufacturing production. How to improve the combustion efficiency of diesel engine and reduce harmful gas emission of diesel engine is a very important work. In the work, an improved heat transfer model in the cylinder was developed in the AVL-BOOST environment. The corresponding entirety model was performed by AVL-BOOST coupled CHEMKINII code which consists of 475 reactions and 134 species, then validated by the experiment. Finally, the propulsion and load characteristics of diesel engine fueled with biodiesel fuel were simulated by the improved heat transfer model. In addition, the effect of fuel injection rate on the engine characteristic of diesel engine fueled with RME was investigated in terms of exhaust temperature, BSFC, soot, and NOx emission. The main conclusions are as follows: (1)Compared with the Woschni 1978 model, the improved model is more accurate in prediction. The maximum difference between the experiment and Woschni 1978 model is reduced. It is due to the fact that the improved model combines the advantages of the two models(2)The higher engine power will increase the effective torque, exhaust temperature, and temperature of behind the intercooler. More specifically, the calculated results are very similar and the maximum difference between the experiment and the improved model is lower than the maximum difference between the experiment and the Woschni 1978 model(3)Compared with case 3, the effective power and the effective torque are increased by 1.02% and 1.023%, respectively; the soot emission decreased by 7.4%, but the nitrogen oxides only increased by 1.1% in case 5. Based on an overall consideration of various factors, case 5 is the most reasonable

Data Availability

All data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This work is supported by the Natural Science Foundation of Guangxi under the research grant 2018GXNSFAA281267, the Guangxi Young and middle-aged college Teachers Basic Research Ability Promotion Project under the research grant 2020KY39008, and the Key Scientific Research Projects of Guangxi Vocational & Technical Institute of Industry under the research grant 2019015KY025.