Prediction of Reformed Gas Composition for Diesel Engines with a Reformed EGR System Using an Artiﬁcial Neural Network

: Facing the reinforced emission regulations and moving toward a clean powertrain, hydrogen has become one of the alternative fuels for the internal combustion engine. In this study, the prediction methodology of hydrogen yield by on-board fuel reforming under a diesel engine is introduced. An engine dynamometer test was performed, resulting in reduced particulate matter (PM) and NOx emission with an on-board reformer. Based on test results, the reformed gas production rate from the on-board reformer was trained and predicted using an artiﬁcial neural network with a backpropagation process at various operating conditions. Additional test points were used to verify predicted results, and sensitivity analysis was performed to obtain dominant parameters. As a result, the temperature at the reformer outlet and oxygen concentration is the most dominant parameters to predict reformed gas owing to auto-thermal reforming driven by partial oxidation reforming process, dominantly.


Introduction
Recently, eco-friendly energy sources, such as electricity and fuel cells, which can replace internal combustion engines, have attracted significant attention due to extreme environmental pollution. However, eco-friendly energy sources involve several problems, such as insufficient driving range, lack of charging stations, and high vehicle prices. Researchers from MIT predicted that 60% of the world's vehicles will be powered by internal combustion engines until 2050 [1]. Therefore, research on internal combustion engines and renewable energy is continuously performed. Specifically, diesel engines, among internal combustion engines, are effective in reducing CO 2 [2]. They offer many advantages, such as high efficiency and torque. However, they also emit harmful exhaust gases, such as NOx and particulate matter (PM) [3]. NOx causes acid rain, and it is converted into nitrate when it comes in contact with moisture in air. This can cause respiratory diseases [4,5] in humans. Additionally, NOx is a precursor of fine dust, and thus, it reacts in air and forms extremely small particles. Diesel engines also emit significant amounts of PM [4]. To date, considerable research has been conducted to mitigate these emissions. Specifically, devices such as diesel particulate filters (DPFs) can be used for removing particulate matter. Furthermore, DPFs can be used for collecting most types of PM [6]. DPFs collect and burn the particles through post-injection. Exhaust gas recirculation (EGR) is a method in which a part of the exhaust gas is sent back to the intake air to return high-specific-heat pressure exhaust gas recirculation) valve was set to be inoperative, and the EGR rate was adjusted solely via the LP EGR valve and backpressure valve. The injection timing was adjusted such that the combustion rate was 50% between after top dead center (ATDC) 8 and 9 to set the combustion efficiency standard. The specifications of the engine are listed in Table 1.
Energies 2020, 13, x FOR PEER REVIEW 3 of 18 The HP EGR (High pressure exhaust gas recirculation) valve was set to be inoperative, and the EGR rate was adjusted solely via the LP EGR valve and backpressure valve. The injection timing was adjusted such that the combustion rate was 50% between after top dead center (ATDC) 8 and 9 to set the combustion efficiency standard. The specifications of the engine are listed in Table 1.  The sample injected from the sample inlet of the gas chromatograph was vaporized via heating and sent to the column filled with the carrier gas together with the gas phase in the liquid phase. The carrier gas moved to the outlet and came in contact with the filler while passing through the column and repeating the adsorption and distribution. At this time, a difference existed between the speeds at which each component in the sample moved through the column due to the difference in the partition coefficient from the moving phase to the fixed phase. Each component that was separated based on the speed of movement was detected by the detector, and a chromatogram based on the elution order was obtained. The calibration gas (consisting of 5% hydrogen, 5% CO, and 90% argon) was measured, as shown in Figure 2.  The sample injected from the sample inlet of the gas chromatograph was vaporized via heating and sent to the column filled with the carrier gas together with the gas phase in the liquid phase. The carrier gas moved to the outlet and came in contact with the filler while passing through the column and repeating the adsorption and distribution. At this time, a difference existed between the speeds at which each component in the sample moved through the column due to the difference in the partition coefficient from the moving phase to the fixed phase. Each component that was separated based on the speed of movement was detected by the detector, and a chromatogram based on the elution order was obtained. The calibration gas (consisting of 5% hydrogen, 5% CO, and 90% argon) was measured, as shown in Figure 2.  Figure 2a shows the area under the 5% hydrogen and 5% carbon monoxide peaks. The X-axis represents time, and the y-axis represents μV. The reformed gas obtained from downstream of the reforming device was measured, and the area under the peaks of the reformed gas was obtained, as shown in Figure 2b. Hence, it is possible to calculate the volume ratio of hydrogen and carbon monoxide contained in the exhaust gas.
A MEXA-9100DEGR gas analyzer (Horiba, Kyoto, Japan) was used. The analyzer can measure the concentrations of CO, CO2, total hydrocarbon (THC), O2, and NOx. Additionally, the EGR rate was calculated using the ECO2 value, which corresponds to the CO2 concentration in the intake air. The detailed specifications of the gas analyzer are listed in Table 2. The total hydrocarbon (THC)  Figure 2a shows the area under the 5% hydrogen and 5% carbon monoxide peaks. The X-axis represents time, and the y-axis represents µV. The reformed gas obtained from downstream of the reforming device was measured, and the area under the peaks of the reformed gas was obtained, as shown in Figure 2b. Hence, it is possible to calculate the volume ratio of hydrogen and carbon monoxide contained in the exhaust gas.
A MEXA-9100DEGR gas analyzer (Horiba, Kyoto, Japan) was used. The analyzer can measure the concentrations of CO, CO 2 , total hydrocarbon (THC), O 2 , and NOx. Additionally, the EGR rate was calculated using the ECO 2 value, which corresponds to the CO 2 concentration in the intake air. The detailed specifications of the gas analyzer are listed in Table 2. The total hydrocarbon (THC) concentration was measured using the number of C. Because we used C 3 H 8 gas, we set the span value by multiplying the value of the calibration gas by 3. The reforming performance improves when the fuel is atomized and uniformly distributed in the catalyst [22]. Therefore, the fuel in this study was atomized using an ultrasonic injector. An ultrasonic injector, S 40 (CERATORQ, Seoul, Korea), was used, and a controller was used to control the flow rate based on the ambient temperature and pressure conditions, as shown in Figure 3.
Energies 2020, 13, x FOR PEER REVIEW 5 of 18 concentration was measured using the number of C. Because we used C3H8 gas, we set the span value by multiplying the value of the calibration gas by 3. The reforming performance improves when the fuel is atomized and uniformly distributed in the catalyst [22]. Therefore, the fuel in this study was atomized using an ultrasonic injector. An ultrasonic injector, S 40 (CERATORQ, Seoul, Korea), was used, and a controller was used to control the flow rate based on the ambient temperature and pressure conditions, as shown in Figure 3. The in-cylinder pressure was measured using a Kistler's 6065A41(KISTLER Korea, Gyeonggido, Korea) product as opposed to the injector in the cylinder. The product is a glow plug type and can measure in-cylinder pressure in real-time. It can measure pressures in the range of 0-250 bar and has a sensitivity of approximately 20 pC/bar at 350 °C. The combustion pressure in the cylinder was monitored using a product termed Kibox (manufactured by KISTLER Korea, Gyeonggi-do, Korea), and the ignition timing was determined via combustion analysis. The engine specifications were input for precise measurement. The data were collected over 100 cycles, and the average pressure value was used.
A high-performance monolith catalyst with a high durability of 30% (5% Platinum (Pt)-Gadolinia doped Ceria (GDC) + 70% C11PR (commercial pre-reforming catalyst)) was used as the reformer. The hydrogen conversion efficiency of the catalyst was approximately 50%, and its durability was confirmed to be excellent based on a 100 h durability test performed after mounting the reformer on the diesel engine. The catalyst housing was designed and manufactured, as shown The in-cylinder pressure was measured using a Kistler's 6065A41(KISTLER Korea, Gyeonggi-do, Korea) product as opposed to the injector in the cylinder. The product is a glow plug type and can measure in-cylinder pressure in real-time. It can measure pressures in the range of 0-250 bar and has a sensitivity of approximately 20 pC/bar at 350 • C. The combustion pressure in the cylinder was monitored using a product termed Kibox (manufactured by KISTLER Korea, Gyeonggi-do, Korea), and the ignition timing was determined via combustion analysis. The engine specifications were input for precise measurement. The data were collected over 100 cycles, and the average pressure value was used.
A high-performance monolith catalyst with a high durability of 30% (5% Platinum (Pt)-Gadolinia doped Ceria (GDC) + 70% C11PR (commercial pre-reforming catalyst)) was used as the reformer. The hydrogen conversion efficiency of the catalyst was approximately 50%, and its durability was confirmed to be excellent based on a 100 h durability test performed after mounting the reformer on the diesel engine. The catalyst housing was designed and manufactured, as shown in Figure 4. A cartridge heater was used to heat the gas and fuel before they would reach the catalyst. In the lower part with the catalyst, a 4 kW ceramic band heater with ceramic insulation was used. The catalyst reacted via steam reforming, partial oxidation reforming, water shift reforming, and auto-thermal reforming. Thus, the reforming was performed through a highly complicated process, although the reforming efficiency was high. Well-known hydrocarbon fuel reforming reactions over Pt catalysts can be found in [23] referred to as auto-thermal reforming, which is combined reforming between steam reforming and partial oxidation.
in Figure 4. A cartridge heater was used to heat the gas and fuel before they would reach the catalyst. In the lower part with the catalyst, a 4 kW ceramic band heater with ceramic insulation was used. The catalyst reacted via steam reforming, partial oxidation reforming, water shift reforming, and auto-thermal reforming. Thus, the reforming was performed through a highly complicated process, although the reforming efficiency was high. Well-known hydrocarbon fuel reforming reactions over Pt catalysts can be found in [23] referred to as auto-thermal reforming, which is combined reforming between steam reforming and partial oxidation.
A 4 kW ceramic heater and four 500 W cartridge heaters were used to maintain the catalyst activation temperature of the reformer. The ceramic heater was controlled on/off to maintain the temperature, and a proportional-integral-derivative (PID) controller was used in the cartridge heater for precise temperature control. Specifically, a ceramic case was prepared and wrapped around the cartridge heater to prevent oxidation of the fuel. This ensured that the fuel did not come into contact with the cartridge heaters before reaching the catalyst. Additionally, the cartridge heater was controlled using a PID controller. In a further study, an additional heater used in this study, which can decrease system efficiency, can be substituted for a waste heat recovery strategy using hot temperature from the EGR cooler.

Experimental Procedure
Experiments were performed in two ways. The operating conditions of the engine are listed in Table 3.  Table 3. In Figure 5, the box region corresponds to the desired A 4 kW ceramic heater and four 500 W cartridge heaters were used to maintain the catalyst activation temperature of the reformer. The ceramic heater was controlled on/off to maintain the temperature, and a proportional-integral-derivative (PID) controller was used in the cartridge heater for precise temperature control. Specifically, a ceramic case was prepared and wrapped around the cartridge heater to prevent oxidation of the fuel. This ensured that the fuel did not come into contact with the cartridge heaters before reaching the catalyst. Additionally, the cartridge heater was controlled using a PID controller. In a further study, an additional heater used in this study, which can decrease system efficiency, can be substituted for a waste heat recovery strategy using hot temperature from the EGR cooler.

Experimental Procedure
Experiments were performed in two ways. The operating conditions of the engine are listed in Table 3.  Table 3. In Figure 5, the box region corresponds to the desired region where the reforming catalyst can be used, as suggested by the OEM (Original equipment manufacturer). Oxygen is required for the partial oxidation reforming reaction; thus, the amount of catalyst that could be used in the high-load region was limited. Furthermore, the temperature was low in an exceedingly-low-load region, and thus, the activation temperature of the catalyst was also limited to a certain value. The hatched area corresponds to the experimental area that was used to train the ANN model, while the red dot denotes the point at which additional experiments were performed under other operating conditions that were not used for training. The yellow dot denotes the point used for the ANN model training and verification.
Energies 2020, 13, x FOR PEER REVIEW 7 of 18 region where the reforming catalyst can be used, as suggested by the OEM (Original equipment manufacturer). Oxygen is required for the partial oxidation reforming reaction; thus, the amount of catalyst that could be used in the high-load region was limited. Furthermore, the temperature was low in an exceedingly-low-load region, and thus, the activation temperature of the catalyst was also limited to a certain value. The hatched area corresponds to the experimental area that was used to train the ANN model, while the red dot denotes the point at which additional experiments were performed under other operating conditions that were not used for training. The yellow dot denotes the point used for the ANN model training and verification. The engine was controlled through a map of the ECU (Engine control unit or electronic control unit), and only the LP-EGR was controlled using a separate backpressure valve and an EGR valve. The experiment was performed with an EGR rate of 10%. First, the engine test was performed without heating the catalyst. Experimental data were obtained for each of the six BMEP values, and experiments were performed using the reformer to compare the results. The temperature of the reforming catalyst was controlled to between 500 °C and 600 °C via the waste heat of the exhaust gas and heater. Furthermore, the amount of fuel injected into the combustion chamber was adjusted to maintain the output because the output increased due to the addition of the reforming gas. The EGR rate and supplied-hydrogen ratio were calculated using the following equations: The EGR rate was calculated via Equation (1) using the ratio of the total CO2 concentration and CO2 concentration in the intake. The concentration of hydrogen was calculated as the ratio of the The engine was controlled through a map of the ECU (Engine control unit or electronic control unit), and only the LP-EGR was controlled using a separate backpressure valve and an EGR valve. The experiment was performed with an EGR rate of 10%. First, the engine test was performed without heating the catalyst. Experimental data were obtained for each of the six BMEP values, and experiments were performed using the reformer to compare the results. The temperature of the reforming catalyst was controlled to between 500 • C and 600 • C via the waste heat of the exhaust gas and heater. Furthermore, the amount of fuel injected into the combustion chamber was adjusted to maintain the output because the output increased due to the addition of the reforming gas. The EGR rate and supplied-hydrogen ratio were calculated using the following equations: EGR rate(%) = CO 2 (intake port) CO 2 (exhaust port) × 100 (1) The EGR rate was calculated via Equation (1) using the ratio of the total CO 2 concentration and CO 2 concentration in the intake. The concentration of hydrogen was calculated as the ratio of the lower heating value of hydrogen (LHV) of the injected total fuel to the LHV value of the injected hydrogen. The concentration of injected hydrogen was calculated using Equation (2).
Additionally, the reforming gas contains CO 2 , and thus, the ECO 2 value slightly increases when the reaction starts. However, this also occurs because of the reforming gas; hence, the experiment was performed without adjusting the valve.

Artificial Neural Network
Neural networks can have structures with various shapes. However, the MLP is the most widely used type of neural network model for data analysis [20]. An MLP neural network model consists of an input layer, a hidden layer, and an output layer. Figure 6 shows the structure of a neural network and the training process used in this study.
Energies 2020, 13, x FOR PEER REVIEW 8 of 18 lower heating value of hydrogen (LHV) of the injected total fuel to the LHV value of the injected hydrogen. The concentration of injected hydrogen was calculated using Equation (2). Additionally, the reforming gas contains CO2, and thus, the ECO2 value slightly increases when the reaction starts. However, this also occurs because of the reforming gas; hence, the experiment was performed without adjusting the valve.

Artificial Neural Network
Neural networks can have structures with various shapes. However, the MLP is the most widely used type of neural network model for data analysis [20]. An MLP neural network model consists of an input layer, a hidden layer, and an output layer. Figure 6 shows the structure of a neural network and the training process used in this study. In this study, a backpropagation neural network was trained as follows: • Step 1: The input value was applied to the input node in the neural network, and multiplication or addition with the weight was repeated several times to calculate the output value based on the input. • Step 2: The output value differed from the target output value; thus, the error between the predicted output and target output was calculated in the learning process in the neural network. • Step 3: The weights were reduced or increased to reduce the error. • Step 4: The extent to which each weight should be changed was determined. • Step 5: The weight was updated to the value determined in Step 4. The processing of the neural network proceeded from the input layer to the hidden layer and then to the output layer. However, the weight update proceeded from the output layer to the hidden layer and then to the input layer. Hence, this process is termed backpropagation. • Step 6: Steps 1-5 were repeated until the error decreased to an appropriate level for all the learning data.
The learning model was constructed using MATLAB ® (V 2017, Natick, MA, USA)'s Neural Net fitting application. The parameters affecting the catalyst were selected as the input variables for the In this study, a backpropagation neural network was trained as follows: • Step 1: The input value was applied to the input node in the neural network, and multiplication or addition with the weight was repeated several times to calculate the output value based on the input.

•
Step 2: The output value differed from the target output value; thus, the error between the predicted output and target output was calculated in the learning process in the neural network.

•
Step 3: The weights were reduced or increased to reduce the error.

•
Step 4: The extent to which each weight should be changed was determined.

•
Step 5: The weight was updated to the value determined in Step 4. The processing of the neural network proceeded from the input layer to the hidden layer and then to the output layer. However, the weight update proceeded from the output layer to the hidden layer and then to the input layer. Hence, this process is termed backpropagation.

•
Step 6: Steps 1-5 were repeated until the error decreased to an appropriate level for all the learning data.
The learning model was constructed using MATLAB ® (V 2017, Natick, MA, USA)'s Neural Net fitting application. The parameters affecting the catalyst were selected as the input variables for the learning ANN model. Temperature significantly affects the activation of the catalyst. Therefore, the temperature of the front, middle, and rear ends of the catalyst housing were selected as the variables. Additionally, the amount of oxygen discharged from the engine was selected as a variable because the oxygen content significantly affects partial oxidation reforming. Specifically, THC, CO, and NOx emissions decreased when the reformed gas containing hydrogen was supplied to the cylinders. Moreover, the catalyst was significantly affected by the gas hourly space velocity (GHSV) value [24,25]. Therefore, the engine speed, load, and EGR rate (which are directly related to GHSV) constituted the variables. Hydrogen and CO were used as output variables and had the greatest influence on combustion among the components of the reformed gas. Hence, 11 parameters were selected as input variables for learning based on the test data from the operating conditions listed in Table 3.
Based on the measured data and the aforementioned input parameters, the ANN model was trained using 70% of the measured data as training data, 15% as test data (which was not used to train the ANN), and 15% as validation data to compensate for the error in the backpropagation process.
The accuracy and quality of the ANN model were assessed in terms of the root mean square error (RMSE), R 2 coefficient, and mean absolute prediction error (MAPE). RMSE denotes the average difference between predicted and experimental data. R 2 indicates the extent to which the regression line represents the actual dataset; it can vary from 0 to 1. In this study, R 2 was close to 1, and thus, it represents a perfect ANN prediction model. MAPE indicates the prediction error of the ANN model. The three parameters can be determined using the following equations [26]: The trained model was verified under the different operating conditions listed in Table 3. Simulink ® was used to construct the calculation control logic. If 11 input data were entered from signal 1 to signal 11 (engine speed, load, temperature, EGR rate, etc.), the neural network model predicted the values of hydrogen and carbon monoxide as the output values. Figure 7 shows the experimental results for cylinder pressure against the crank angle without and with reformed gas addition under 1500 rpm, for a BMEP of 6 bar. The hydrogen content was 0.5-1.2% by energy. It was confirmed that the ignition was delayed, and the maximum pressure decreased when the reformed gas was injected. The addition of the reformed gas led to a further reduction in NOx when compared with H2 addition alone. The results indicated that CO2, CO, and other gases with high heat capacity were generated and that the temperature of the cylinder decreased. The experimental results for NOx emission are shown in Figure 8.

Effect of Reformed Gas Addition
PM emission also decreased significantly. The number, as well as mass of the smaller particulates, decreased. The size of the larger particulates decreased when the oxidation was accelerated, and the smaller particulates quickly disappeared. Specifically, the PM emission decreased at high loads, and this was considered to promote PM oxidation because the temperature of the cylinder under high loads was higher than that under low loads. The PM emission results are shown in Figures 9 and 10.   The addition of the reformed gas led to a further reduction in NOx when compared with H2 addition alone. The results indicated that CO2, CO, and other gases with high heat capacity were generated and that the temperature of the cylinder decreased. The experimental results for NOx emission are shown in Figure 8.
PM emission also decreased significantly. The number, as well as mass of the smaller particulates, decreased. The size of the larger particulates decreased when the oxidation was accelerated, and the smaller particulates quickly disappeared. Specifically, the PM emission decreased at high loads, and this was considered to promote PM oxidation because the temperature of the cylinder under high loads was higher than that under low loads. The PM emission results are shown in Figures 9 and 10.

ANN Training Results
A greater number of neurons in the learning model increases the probability of over-fitting, which consequently reduces data prediction accuracy because the training model focuses only on training data. To prevent over-fitting, we constructed a learning model after comparing the number of neurons from 1 to 15. A Bayesian algorithm was used because the number of data points was low. The Bayesian algorithm is used when the number of data points is low. The most accurate prediction was obtained when the number of neurons was 13. Therefore, the number of neurons for the learning model used in this study was set as 13. Figure 11 shows the regression results with 13

ANN Training Results
A greater number of neurons in the learning model increases the probability of over-fitting, which consequently reduces data prediction accuracy because the training model focuses only on training data. To prevent over-fitting, we constructed a learning model after comparing the number of neurons from 1 to 15. A Bayesian algorithm was used because the number of data points was low. The Bayesian algorithm is used when the number of data points is low. The most accurate prediction was obtained when the number of neurons was 13. Therefore, the number of neurons for the learning model used in this study was set as 13. Figure 11 shows the regression results with 13 PM emission also decreased significantly. The number, as well as mass of the smaller particulates, decreased. The size of the larger particulates decreased when the oxidation was accelerated, and the smaller particulates quickly disappeared. Specifically, the PM emission decreased at high loads, and this was considered to promote PM oxidation because the temperature of the cylinder under high loads was higher than that under low loads. The PM emission results are shown in Figures 9 and 10.

ANN Training Results
A greater number of neurons in the learning model increases the probability of over-fitting, which consequently reduces data prediction accuracy because the training model focuses only on training data. To prevent over-fitting, we constructed a learning model after comparing the number of neurons from 1 to 15. A Bayesian algorithm was used because the number of data points was low. The Bayesian algorithm is used when the number of data points is low. The most accurate prediction was obtained when the number of neurons was 13. Therefore, the number of neurons for the learning model used Energies 2020, 13, 5886 12 of 17 in this study was set as 13. Figure 11 shows the regression results with 13 neurons with respect to training data, test data, and integrated data between training and test data. As shown in the figure, suitable regression values were obtained with 13 neurons.
Energies 2020, 13, x FOR PEER REVIEW 12 of 18 neurons with respect to training data, test data, and integrated data between training and test data. As shown in the figure, suitable regression values were obtained with 13 neurons. The actual measured amounts of H2 and CO in the reformed gas and the amounts predicted using the neural network model are shown in Figure 12. The sample number on the x-axis represents the number of tests based on the operating conditions listed in Table 3. The actual measured amounts of H 2 and CO in the reformed gas and the amounts predicted using the neural network model are shown in Figure 12. The sample number on the x-axis represents the number of tests based on the operating conditions listed in Table 3.   Figure 13 shows a comparison between the predicted and target values for the ANN model verification under the operating conditions listed in Table 3. As shown in the figure, the H2 reforming tendencies were well-matched at 1500 rpm, which is within the operating range used for training; however, the tendencies were underestimated below 2200 and 2400 rpm, which was beyond the operating range used in training. In contrast, the CO reforming tendencies were wellmatched under every operating condition. However, the values were underestimated below 2200 and 2400 rpm, which was beyond the operating range used in training. This indicates that ANN  Figure 13 shows a comparison between the predicted and target values for the ANN model verification under the operating conditions listed in Table 3. As shown in the figure, the H 2 reforming tendencies were well-matched at 1500 rpm, which is within the operating range used for training; however, the tendencies were underestimated below 2200 and 2400 rpm, which was beyond the operating range used in training. In contrast, the CO reforming tendencies were well-matched under every operating condition. However, the values were underestimated below 2200 and 2400 rpm, which was beyond the operating range used in training. This indicates that ANN training with extremely few operating conditions (six in this study) leads to a reduction in prediction accuracy despite the optimization of the number of neurons. Furthermore, it is necessary to consider more factors that affect the reforming reaction under high-rpm and high-load operating conditions. Energies 2020, 13, x FOR PEER REVIEW 14 of 18 training with extremely few operating conditions (six in this study) leads to a reduction in prediction accuracy despite the optimization of the number of neurons. Furthermore, it is necessary to consider more factors that affect the reforming reaction under high-rpm and high-load operating conditions. A sensitivity analysis was performed to determine the most sensitive of the selected 11 variables. The sensitivity of each parameter was compared to determine the variables that affect the catalyst performance most significantly. The variables were excepted one by one during the A sensitivity analysis was performed to determine the most sensitive of the selected 11 variables. The sensitivity of each parameter was compared to determine the variables that affect the catalyst performance most significantly. The variables were excepted one by one during the training process, Energies 2020, 13, 5886 15 of 17 and the training results were compared with the test data. How far the trained data was from the test data can be deduced as the R-value shown in Table 4. As shown in Table 4, the largest error between the training and test data was observed when the T_housing_out referring temperature at the downstream of the catalyst was excepted. The sensitivity ranking results indicate that T_housing_out and O 2 are the most effective variables for improving reforming performance. After all, results from sensitivity analysis suggest that the main reforming reaction in this study is mainly driven by partial oxidation, which is reasonable.

Conclusions
In this study, a reformed gas addition method using a reforming catalyst was employed to reduce harmful emissions in the exhaust gas of a diesel engine. Experiments were performed using a 1.6-L engine and a reforming catalyst on the LP-EGR line. A neural network was used to predict the composition of the reforming gas, which changes continuously in real-time. The results are summarized as follows: 1.
The addition of a reforming gas to a diesel engine can reduce NOx emissions. A longer ignition delay reduces the maximum pressure and temperature in the combustion chamber; this, in turn, mitigates the production of NOx, which occurs mainly at high temperatures following the Zeldovich mechanism. Furthermore, the addition of reforming gas increases the premixing rate because the gas enters the intake line through the LP EGR line. Therefore, the rich region decreases, and the lean region temperature rises, which promotes PM oxidation.

2.
The EGR rate, catalyst temperature, exhaust gas composition, and engine operating conditions were selected as input data for the ANN learning model. The amount of hydrogen and carbon monoxide production was selected as the output data. The training model was composed of 134 datasets, and the prediction accuracy is 92.2%. When the number of neurons was 13, the model exhibited the highest prediction accuracy. 3.
The ANN model exhibited a low prediction accuracy for experimental conditions that were not used for training. Therefore, a prediction model for a wider range of operating conditions can be constructed by acquiring and learning more data. 4.
The sensitivities of the variables affecting the accuracy of the catalyst performance prediction model were compared. Partial oxidation reforming is an exothermic reaction, and thus, a high temperature at the rear end of the catalyst indicates effective reforming. Therefore, the temperature and amount of oxygen downstream of the catalyst affect the reforming performance most significantly.