Mathematical Model of Thermodynamic Processes in the Intake Manifold of a Diesel Engine with Fuel and Water Injection

: The article presents the results of a study aimed at creating a mathematical model of thermodynamic processes in the intake manifold of a forced diesel engine, taking into account the features of simultaneous injection of fuel and water into the collector. In the course of the study, the tasks of developing a mathematical model were solved, it was implemented in the existing software for component simulation “Internal combustion engine research and development” (ICE RnD), created using the Modelica language, and veriﬁcation was undertaken using the results of bench tests of diesel engines with injection fuel and water into the intake manifold. The mathematical model is based on a system of equations for the energy and mass balances of gases and includes detailed mathematical submodels of the processes of simultaneous evaporation of fuel and water in the intake manifold; it takes into account the e ﬀ ect of the evaporation of fuel and water on the parameters of the gas state in the intake manifold; it takes into account the inﬂuence of the state parameters of the working ﬂuid in the intake manifold on the physical characteristics of fuel and water; it meets the principles of component modeling, since it does not contain parameters that are not related to the simulated component; it describes the process of simultaneous transfer of vapors and non-evaporated liquids between components; and it does not include empirical relationships requiring data on the dynamics of fuel evaporation under reference conditions. According to the results of a full-scale experiment, the adequacy of the mathematical model developed was conﬁrmed. This model can be used to determine the rational design parameters of the fuel and water injection system, the adjusting parameters of the forced diesel engine that provide the required power, and economic indicators, taking into account the limitations on the magnitude of the mechanical and thermal loads of its parts.


Introduction
The modern economy requires the maximum reduction of time and money spent on the development and preparation of production, the analysis of the object's technical condition at all stages of the life cycle, and the fast introduction of changes to the products design parameters and operation modes [1]. The development and implementation of digital twin technology can significantly increase the effectiveness of research and development by maximizing the replacement of costly and time-consuming real product tests with virtual ones without sacrificing their quality. This fully applies to reciprocating internal combustion engines-the main sources of mechanical energy for operating land vehicles [2]. Most researchers believe that the model should implement simulation methods [3], In theory it is possible to create custom components in which the processes of simultaneous evaporation of fuel and water are adequately simulated; however, it is impossible to simulate the processes of transfer of fuel droplets between components, since this requires changes to the basic foundations implemented in programs, which is impossible because of their proprietary code.
Thus, solving the problem of synthesizing simulation models of engines boosted by simultaneous injection of water and fuel into the intake manifold is difficult due to the lack of reliable mathematical models describing the processes of simultaneous evaporation of water and fuel, transfer of their vapors and non-evaporated liquid between components. The aim of the study was to solve the scientific problem of developing such a model. This model can be implemented into well-known mathematical models describing in a one-dimensional formulation the processes in the gas-air duct, which are developed, validated and implemented in the form of software (that is, the task of developing a complete mathematical model of the engine is not undertaken).
The model developed includes detailed mathematical submodels of the processes of simultaneous evaporation of fuel and water in the intake manifold; it takes into account the effect of the evaporation of fuel and water on the parameters of the gas state in the intake manifold; it takes into account the influence of the parameters of the state of the working fluid in the intake manifold on the physical characteristics of fuel and water; it meets the principles of component modeling, since it does not contain parameters that are not related to the simulated component; it describes the process of simultaneous transfer of vapors and non-evaporated liquids between components; and it is based on basic equations of physics, and not on empirical dependencies that require data on the dynamics of fuel evaporation under reference conditions. Together, these features of the mathematical model distinguish it from the known ones.

Developed Mathematical Model
Analysis of the existing methods of gas and thermodynamics, known mathematical models for one-dimensional gas-dynamic calculation of the gas-air path of piston engines, indicates the appropriateness of using phenomenological models such as Lyshevsky's [14], Razleitseva's [15], Kuleshov's [16], Hiroyasu's [17] for modeling processes in the intake manifold of a diesel engine with fuel and/or water injection, developed for the conditions of fuel injection into the combustion chamber. According to these models, the camera volume is divided into separately modeled zones in which interconnected processes take place. For example, according to Kuleshov, the relative rate of evaporation of the liquid for the i-th zone: where τ is time; τ s is the current time; τ s0i is the time of fuel getting into i-th zone; τ ui is current time from the start of injection; σ zi is the fraction of fuel in the i-th zone; K i is the evaporation constant in the i-th zone; D 32 is the Sauter average droplet diameter; b ui is relative constant of evaporation in the i-th zone; Y is a function that corrects unaccounted factors and errors in assumptions made. The total evaporation rate is defined as the sum of the evaporation rates in all zones. The constants of fuel evaporation in different zones are determined as follows: where: Nu D is the Nusselt criterion for the diffusion process; D p is the vapor diffusion coefficient; D p0 is vapor diffusion coefficient under the conditions p 0 = 0.1 MPa, T 0 = 293 K; T k is the temperature of the drop; p is the current pressure; ρ l is the fluid density; p s is the vapor pressure. Most researchers assume that the initial diameter of a liquid droplet is equal to the diameter of the spray hole [18]. To describe the frequency curve of the droplet size distribution, the P. Rosin and E. Rammler equation [19] is well suited: where: z is the characteristic droplet size; i is the number of drops of size z; n and a m are empirical indicators.
In models similar to (1)(2)(3)(4)(5)(6), periodic injection of liquid is considered and the diameter of all drops is assumed to be the same at every moment of time and changing simultaneously. According to Razleitsev [20], under diesel conditions, the dispersion of spraying is characterized by an average droplet diameter, which makes it possible to significantly simplify the task of modeling the dynamics of liquid evaporation in the intake manifold. However, the model of Razleitsev does not explicitly show the temperature of the medium, in addition, it is designed for cyclical rather than continuous fuel supply (i.e., τ evap → ∞). It can be taken as the basis for a model of fuel and water injection into the intake manifold, since, unlike other known models (AVL Boost, GT-Suite, etc.), and it takes into account the characteristics of the atomizer, the fuel density and the parameters of the gas state. In the simulated volume, fully complies with the principles of the component approach in the synthesis of simulation models, does not require the determination of parameters that depend on the characteristics of the gas flow in the reservoir.
The basis of the mathematical model of thermodynamic processes in the intake manifold of a diesel engine is a system of equations of energy and mass balances [21] of a mixture of air, exhaust gases (entering the collector when being thrown from the combustion chamber), water vapor and fuel (hereafter gas): where m is the manifold gas mass; U is the internal gas energy; u is specific internal energy of gas; dH dτ is the enthalpy flow due to mass transfer; dQ int dτ is the drain of energy due to heating and evaporation of water and fuel; dQ wall dτ is the energy flow through the walls; dU dτ is the change in the internal energy of the gas.
where dm in dτ is the mass flow of gas flowing into the reservoir; dm out dτ is the gas mass flow from the manifold; dm int dτ is the change in mass of gas due to the evaporation of water and fuel. The system of equations of energy and mass gas balances (7) and (8)  The amount of apparent heat spent on heating the fuel droplets to a boiling point: where m l f is the mass of liquid fuel; c l f is the specific heat of liquid fuel; T e f is the boiling point of the fuel; T l f is the temperature of the liquid fuel. The amount of latent heat spent on the evaporation of fuel: where m e f is the mass of evaporated fuel; q b f is the specific heat of vaporization of the fuel. The equations for water can be defined similar to (10) and (11): where m lw is the mass of liquid water; c lw is the specific heat of liquid water; T ew is the boiling point of water; T lw is the temperature of liquid water; m ew is the mass of evaporated water; q bw is the specific heat of vaporization of water. Enthalpy flow due to mass transfer: where h in and h out are the specific enthalpy of inflowing and outflowing gas, respectively. The physical properties of a mixture of air, water vapor and fuel depend on its temperature. For a mixture of air, water vapor and fuel in the intake manifold, the average values of the specific physical indicators are determined using the weighted average.
The boiling point of water, depending on pressure (MPa): Since diesel fuel is a mixture of hydrocarbons having a different boiling point (180-360 • C [22]), when performing calculations it is recommended to use a temperature at which 50% of the fractions are distilled off (for diesel summer and winter fuels T ef = 280 • C). For water in the liquid phase, the specific heat capacity in the temperature range 0-90 • C varies insignificantly (from 4213.8 J/(kg·K) to 4205.8 J/(kg·K) [23]), therefore it can be assumed constant and equal to c v = 4210 J/(kg·K). For diesel fuel in the liquid phase, the specific heat is: The fuel temperature in Equation (17) can be taken equal to the arithmetic mean between the initial and T ew . The specific heat of water vaporization in the pressure range 0.1-0.5 MPa is approximated by the linear dependence [24]: q bw = (1.54 − 0.1241 · p) · 10 6 J/kg.
For diesel fuel in the liquid phase, the specific heat of vaporization at the boiling point is: Energies 2020, 13, 4315 6 of 18 To simulate the process of evaporation of jets of liquids in the intake manifold, we used the assumption that the dynamics of the process of evaporation of a liquid repeats the dynamics of its injection, but with a time shift equal to the time of complete evaporation of the droplet: The time of complete evaporation of a drop is inversely proportional to the relative constant of evaporation: (20) which, in turn, is determined by the Razleitsev method adapted for the conditions of continuous injection [15]. The average rate of fluid flow from the nozzle atomizer with continuous injection: where ρ l is the density of the liquid; D c is the diameter of the nozzle holes; i c is the number of coop holes. M-criteria: where µ l is the dynamic viscosity coefficient of the liquid; ρ g is the gas density in the combustion chamber; σ l is the surface tension coefficient of the liquid. Weber criteria: The average surface diameter of liquid droplets: where ρ g is the gas density in the combustion chamber; E 32 = 1.7 is the empirical coefficient in the original Razleytsev's methodology [15]. Relative constant of fuel evaporation: where k evap is the constant depending on evaporation conditions (assumed k evap = 1). In Equations (21)-(25), the physical characteristics of a specific liquid (fuel or water) should be used, including as a function of temperature. The dependencies obtained on the basis of the IAPWS (International Association for the Properties of Water and Steam) formulations are the most accurate, however, according to Lyshevsky [14], and in the range of conditions of the combustion chamber and especially the intake manifold, these properties can be considered constant.
When calculating, it must be taken into account that not all injected fuel or water can evaporate in the intake manifold. To control this situation, the time is calculated during which the gas flow passes the distance L inj from the place of installation of the nozzle to the exit of the intake manifold: where U g is flow rate at the nozzle installation site; A inj is cross-sectional area of the intake manifold at the nozzle installation site.
If τ e > τ L , the non-evaporated fuel or gas flows into the neighboring element of the gas-air path, where the evaporation process continues.
The temperature of the gases in the manifold is determined from the equation: where → C-component concentration vector; In general terms, the flow of mass of gas flowing in or out of the intake manifold is described by the well-known formula: where µA is effective cross-sectional area; p in and T in are gas pressure and temperature before and after the cross-section, respectively; k is adiabatic exponent; p out -gas pressure over the cross section (must not exceed critical for a given flow regime): For several interconnected elements of the gas-air path, the mathematical model is supplemented by the equations of the balance of gas flow and pressure at the joints: where c is the number of connections; dm c dt -gas consumption through each connection. Through the connection, the values of the temperature variables, the concentration vector of the components and the specific enthalpy of the gas are transmitted, taking into account the direction of flow.

Verification of the Mathematical Model
Developed mathematical models were implemented into developed software. Since the complete model of a diesel engine, in addition to the mathematical description of the gas-air path, should include submodels of processes in the combustion chamber, the crank and gas distribution mechanisms, heat transfer in structural elements and with the environment, a code analyzer-optimizer, a mathematical solver, etc., the "Combustion Engine Research and Development" (ICE RnD) software, previously developed by A.A. Malozyomov and G.A. Malozyomov were taken as the basis. A detailed description of the program and the mathematical models implemented in it is given in [25,26]. The program is written in the language Modelica [27], which makes it easy enough to make changes to the code, and it is a class library from which simulation models of a diesel engine of any configuration are compiled using the graphical interface. The program classes that describe the processes in the elements of the gas-air path (volumes, pipes, connecting channels, combustion chamber) were modernized and supplemented to provide the possibility of modeling a diesel engine with fuel and water injection into the intake manifold [28], namely: Energies 2020, 13, 4315 8 of 18 • the base class "volume" was supplemented by a description of the process of evaporation of fuel and water and its effect on the parameters of the state of gases in the combustion chamber; • in the remaining base classes, including component relationships, additional variables were included-the concentration of fuel and water vapor; • the base class "cylinder" was supplemented by descriptions of the combustion process of fuel vapor coming from the intake manifold, and the effect of water vapor and combustion products of fuel vapor on the combustion performance the Vibe equation of the main portion of fuel injected into the combustion chamber [28].
When calculating a diesel engine with fuel and/or water injection into the intake manifold, the problem of choosing the main design parameters of the fuel and water injection system and rational values of their mass supply is usually solved in order to ensure the specified power values with the minimum possible increase in the maximum values of temperature, pressure and rate of pressure rise gases in the combustion chamber. This article presents the calculation results confirming the adequacy of the developed mathematical model.
The initial data, boundary and initial conditions for the calculation corresponded to the design and operational characteristics of the 12CHN15/18 diesel engine (Uraltrak, Chelyabinsk, Russia). The dimensions of the collector and the gas and thermodynamic conditions in it by default corresponded to the averaged parameters of a real diesel engine operating at maximum power at a nominal crankshaft speed. The diameter of the spraying hole of water is 0.5 mm, of fuel is 0.1 mm, the continuous supply is reduced to one operating cycle of a diesel engine, water is 0.15 g, and fuel is 0.05 g. The temperature of the injected water and fuel is 50 • C. Both nozzles are installed at a distance of 0.6 m in front of the inlet valve seat of the first cylinder. The number of spray holes for both fuel and water is one at a time.
The engine design scheme created in the modified ICE RnD program is shown in Figure 1. The parameters of the gas turbine boost were specified by the boundary conditions; for the simulation of fast-flowing processes lasting less than 1 s, and this is acceptable, taking into account the inertia of the turbocharger. The injection of fuel and water was carried out after stabilization of the diesel parameters, after 2 s. after the start of the calculation. The parameters and indicators of the components of the left and right row of cylinders were assumed to be identical, therefore, the calculation was carried out only for one row, shown in Figure 1. Since the calculation task did not include determining the loads on the elements of the crank mechanism, this simplification is acceptable.
To verify the mathematical model developed, we used the data obtained by M. Grankin [5,28] during experimental studies of a 12ChN15/18 diesel engine equipped with a fuel and water injection system into the intake manifold. Since the developers of the engine are mainly interested in the main parameters of the purpose (power and fuel efficiency) and the thermal-mechanical loading of the parts of the crank mechanism and the cylinder-piston group, during the experiment, the parameters of the purpose and indicator diagrams of the gas pressure in the combustion chamber were determined.
The appearance of the nozzle for injection of water and fuel is shown in Figure 2. The experimental conditions corresponded to the settings of the test simulation model. Fuel and water were supplied to the nozzle separately and mixed after leaving the spray hole. The diameter of the spray hole of water is 0.5 mm, fuel is 0.1 mm, continuous supply reduced to one operating cycle of diesel, water is 0.15 g, and fuel is 0.05 g. Both nozzles were installed at a distance of 0.6 m in front of the inlet valve seat of the first cylinder. Figure 3 shows a comparison of the calculated and experimental indicator diagrams of the gas pressure in the combustion chamber for the maximum diesel power mode at a speed of 2100 min −1 .
To fulfill the restrictions on the maximum gas pressure in the combustion chamber, the fuel injection advance angle was reduced by 4 degrees of crankshaft rotation compared to the base diesel engine.
In addition, the figure shows the calculated and experimental indicator diagrams of the gas pressure of a diesel engine without the injection of fuel and water into the intake manifold. The appearance of the nozzle for injection of water and fuel is shown in Figure 2. The experimental conditions corresponded to the settings of the test simulation model. Fuel and water were supplied to the nozzle separately and mixed after leaving the spray hole. The diameter of the spray hole of water is 0.5 mm, fuel is 0.1 mm, continuous supply reduced to one operating cycle of diesel, water is 0.15 g, and fuel is 0.05 g. Both nozzles were installed at a distance of 0.6 m in front of the inlet valve seat of the first cylinder.   To verify the mathematical model developed, we used the data obtained by M. Grankin [5,28] during experimental studies of a 12ChN15/18 diesel engine equipped with a fuel and water injection system into the intake manifold. Since the developers of the engine are mainly interested in the main parameters of the purpose (power and fuel efficiency) and the thermal-mechanical loading of the parts of the crank mechanism and the cylinder-piston group, during the experiment, the parameters of the purpose and indicator diagrams of the gas pressure in the combustion chamber were determined.
The appearance of the nozzle for injection of water and fuel is shown in Figure 2. The experimental conditions corresponded to the settings of the test simulation model. Fuel and water were supplied to the nozzle separately and mixed after leaving the spray hole. The diameter of the spray hole of water is 0.5 mm, fuel is 0.1 mm, continuous supply reduced to one operating cycle of diesel, water is 0.15 g, and fuel is 0.05 g. Both nozzles were installed at a distance of 0.6 m in front of the inlet valve seat of the first cylinder.   To fulfill the restrictions on the maximum gas pressure in the combustion chamber, the fuel injection advance angle was reduced by 4 degrees of crankshaft rotation compared to the base diesel engine.
In addition, the figure shows the calculated and experimental indicator diagrams of the gas pressure of a diesel engine without the injection of fuel and water into the intake manifold. From the results of the analysis, it was concluded that at the level of significance α = 0.01, there is no significant difference between the experimentally determined values of the gas pressure in the combustion chamber and their simulations; the mathematical model adequately describes the experimental data. The calculated value of the pair correlation coefficient is 0.992, which indicates the degree of reliability of the developed mathematical model sufficient for engineering calculations. Figure 4 shows a comparison of the experimental and calculated values of brake-specific fuel consumption of the engine at speed characteristic modes with and without fuel and water injection into the intake manifold.    From the results of the analysis, it was concluded that at the level of significance α = 0.01, there is no significant difference between the experimentally determined values of the gas pressure in the combustion chamber and their simulations; the mathematical model adequately describes the experimental data. The calculated value of the pair correlation coefficient is 0.992, which indicates the degree of reliability of the developed mathematical model sufficient for engineering calculations. Figure 4 shows a comparison of the experimental and calculated values of brake-specific fuel consumption of the engine at speed characteristic modes with and without fuel and water injection into the intake manifold.  From the data shown in Figures 3 and 4, it follows that at a given level of significance, the pair correlation coefficient of the experimental and calculated values of the quantities is 0.994 and 0.985, respectively. The absolute error in calculating the effective engine power does not exceed 13 kW (1.6%). Figure 5 shows the dynamics of changes in the concentration of water vapor obtained using the developed model and the CFD model created by M. Grankin [28] using the OpenFOAM software [29]. In the course of the calculation, water injection was simulated with a temperature of 50 • C through a hole with a diameter of 0.5 mm and a flow rate of 0.015 kg/s (the moment of injection start-0.1 s from the beginning of calculation) into a pipe with a diameter of 0.14 m and a volume of 0.017 m 3 , into which, under a pressure of 375 kPa air was supplied with a temperature of 530 K, the outlet counter pressure was 365 kPa (which approximately corresponds to the conditions of the intake manifold of the tested diesel engine operating at maximum power). through a hole with a diameter of 0.5 mm and a flow rate of 0.015 kg/s (the moment of injection start-0.1 s from the beginning of calculation) into a pipe with a diameter of 0.14 m and a volume of 0.017 m3, into which, under a pressure of 375 kPa air was supplied with a temperature of 530 K, the outlet counter pressure was 365 kPa (which approximately corresponds to the conditions of the intake manifold of the tested diesel engine operating at maximum power). As can be seen from Figure 5, the developed mathematical model makes it possible to obtain results that are quite close to the data of CFD calculations. The differences are primarily due to the fact that the developed model does not take into account the statistics of the distribution of drops by diameter, that is, the presence of small drops at the initial moment of injection, which begin to evaporate almost immediately.
Thus, the developed mathematical model adequately reflects the features of the processes of an engine equipped with a fuel and water injection system into the intake manifold and can be used in engineering calculations to create such systems.

Calculation Results and Discussion
Analysis of the calculation results, using the developed mathematical model, implemented into a complex simulation model of a diesel engine, makes it possible to reveal the general patterns of the influence of various factors on the characteristics of engine processes with water and fuel injection into the intake manifold, which can be illustrated with specific examples.
A preliminary calculation showed that the farther from the nozzle along the length of the gas-air duct its element is, the greater the delay in the change and subsequent stabilization of the vapor concentration. For the inlet channels of the head of the first (counting from the injector installation site) cylinder, the delay is 0.2-0.3 s. (Figure 6). As can be seen from Figure 5, the developed mathematical model makes it possible to obtain results that are quite close to the data of CFD calculations. The differences are primarily due to the fact that the developed model does not take into account the statistics of the distribution of drops by diameter, that is, the presence of small drops at the initial moment of injection, which begin to evaporate almost immediately.
Thus, the developed mathematical model adequately reflects the features of the processes of an engine equipped with a fuel and water injection system into the intake manifold and can be used in engineering calculations to create such systems.

Calculation Results and Discussion
Analysis of the calculation results, using the developed mathematical model, implemented into a complex simulation model of a diesel engine, makes it possible to reveal the general patterns of the influence of various factors on the characteristics of engine processes with water and fuel injection into the intake manifold, which can be illustrated with specific examples.
A preliminary calculation showed that the farther from the nozzle along the length of the gas-air duct its element is, the greater the delay in the change and subsequent stabilization of the vapor concentration. For the inlet channels of the head of the first (counting from the injector installation site) cylinder, the delay is 0.2-0.3 s. (Figure 6).
The amplitude of the change in vapor concentration depends on the volume of the element, gas and thermodynamic conditions, the largest amplitude (±0.00016 kg/kg or ±2.4%)-at the place where the nozzle is installed. In the channels of the cylinder head, the concentration of vapors is practically constant and decreases by up to 0.8% at the moment the intake valve opens. The slight difference in the rate of concentration change is caused by the longer duration of the water evaporation delay, which is clearly expressed for the location of the nozzles and is practically zero for the inlet of the head. The stepped nature of the change in the concentration of vapors in the inlet channel of the cylinder head is explained by the cyclic opening and closing of valves and, accordingly, the dynamics of gas flow in its volume.
The effect of the evaporation of fuel and water on the temperature of the gases in the elements of the gas-air duct is shown in Figure 7. It should be noted that the average temperature of the gases in front of the inlet bodies is lower than at the place of installation of the nozzles. This is caused by heat loss through the walls of the elements of the gas-air path, which are taken into account in the calculation. The gas pressure does not change as noticeably as the temperature, although some reduction occurs. The amplitude of the change in vapor concentration depends on the volume of the element, gas and thermodynamic conditions, the largest amplitude (±0.00016 kg/kg or ±2.4%)-at the place where the nozzle is installed. In the channels of the cylinder head, the concentration of vapors is practically constant and decreases by up to 0.8% at the moment the intake valve opens. The slight difference in the rate of concentration change is caused by the longer duration of the water evaporation delay, which is clearly expressed for the location of the nozzles and is practically zero for the inlet of the head. The stepped nature of the change in the concentration of vapors in the inlet channel of the cylinder head is explained by the cyclic opening and closing of valves and, accordingly, the dynamics of gas flow in its volume.
The effect of the evaporation of fuel and water on the temperature of the gases in the elements of the gas-air duct is shown in Figure 7. It should be noted that the average temperature of the gases in front of the inlet bodies is lower than at the place of installation of the nozzles. This is caused by heat loss through the walls of the elements of the gas-air path, which are taken into account in the calculation. The gas pressure does not change as noticeably as the temperature, although some reduction occurs.  Figure 8 shows that the duration of the transition process from the moment the fuel and water supply to the intake manifold starts to stabilize the torque value is 0.2-0.3 s. This must be taken into account if the results of a one-dimensional calculation of processes in the gas-air duct are used to set the boundary conditions when simulating processes in a combustion chamber (for example, in CFD calculations).  The amplitude of the change in vapor concentration depends on the volume of the element, gas and thermodynamic conditions, the largest amplitude (±0.00016 kg/kg or ±2.4%)-at the place where the nozzle is installed. In the channels of the cylinder head, the concentration of vapors is practically constant and decreases by up to 0.8% at the moment the intake valve opens. The slight difference in the rate of concentration change is caused by the longer duration of the water evaporation delay, which is clearly expressed for the location of the nozzles and is practically zero for the inlet of the head. The stepped nature of the change in the concentration of vapors in the inlet channel of the cylinder head is explained by the cyclic opening and closing of valves and, accordingly, the dynamics of gas flow in its volume.
The effect of the evaporation of fuel and water on the temperature of the gases in the elements of the gas-air duct is shown in Figure 7. It should be noted that the average temperature of the gases in front of the inlet bodies is lower than at the place of installation of the nozzles. This is caused by heat loss through the walls of the elements of the gas-air path, which are taken into account in the calculation. The gas pressure does not change as noticeably as the temperature, although some reduction occurs.  Figure 8 shows that the duration of the transition process from the moment the fuel and water supply to the intake manifold starts to stabilize the torque value is 0.2-0.3 s. This must be taken into account if the results of a one-dimensional calculation of processes in the gas-air duct are used to set the boundary conditions when simulating processes in a combustion chamber (for example, in CFD calculations).  Figure 8 shows that the duration of the transition process from the moment the fuel and water supply to the intake manifold starts to stabilize the torque value is 0.2-0.3 s. This must be taken into account if the results of a one-dimensional calculation of processes in the gas-air duct are used to set the boundary conditions when simulating processes in a combustion chamber (for example, in CFD calculations).
The period of delay in the evaporation of fuel and water should not exceed the time for drops of fuel to reach the outlet from the manifold. Droplets entering the combustion chamber can adversely affect engine performance due to additional cooling of the cylinder walls due to their evaporation. Based on this condition, as well as depending on the value of the maximum cyclic fluid supply, the maximum permissible diameter of the spray holes and the distance from the nozzles to the outlet from the manifold should be determined. Figure 9 shows the influence of the temperature of gases at the inlet to the intake manifold on the concentration and the period of delay in the evaporation of fuel and water.
For the engine under consideration, the time for fuel droplets to reach the exit from the manifold is 0.024-0.026 s, and the maximum permissible diameter of the spray holes is 0.6 mm, the time for water droplets to exit the manifold is 0.024-0.030 s, the maximum permissible diameter spray holes 0.7 mm. As can be seen from Figure 9, the time of evaporation of liquid droplets in the considered range of temperatures of gases at the inlet to the intake manifold does not exceed the time for the non-evaporated droplets to reach the outlet from the manifold.
When developing an engine, it is important to comply with the limitations on the mechanical loading of its parts, which can be estimated by the value of the maximum gas pressure in the combustion chamber. Figure 10 shows the dependence of the maximum gas pressure in the combustion chamber on the cyclic supply of fuel and water to the intake manifold. Restrictions can be fulfilled by changing The period of delay in the evaporation of fuel and water should not exceed the time for drops of fuel to reach the outlet from the manifold. Droplets entering the combustion chamber can adversely affect engine performance due to additional cooling of the cylinder walls due to their evaporation. Based on this condition, as well as depending on the value of the maximum cyclic fluid supply, the maximum permissible diameter of the spray holes and the distance from the nozzles to the outlet from the manifold should be determined. Figure 9 shows the influence of the temperature of gases at the inlet to the intake manifold on the concentration and the period of delay in the evaporation of fuel and water. For the engine under consideration, the time for fuel droplets to reach the exit from the manifold is 0.024-0.026 s, and the maximum permissible diameter of the spray holes is 0.6 mm, the time for water droplets to exit the manifold is 0.024-0.030 s, the maximum permissible diameter spray holes 0.7 mm. As can be seen from Figure 9, the time of evaporation of liquid droplets in the considered range of temperatures of gases at the inlet to the intake manifold does not exceed the time for the non-evaporated droplets to reach the outlet from the manifold.
When developing an engine, it is important to comply with the limitations on the mechanical loading of its parts, which can be estimated by the value of the maximum gas pressure in the combustion chamber. Figure 10 shows the dependence of the maximum gas pressure in the combustion chamber on the cyclic supply of fuel and water to the intake manifold. Restrictions can be fulfilled by changing the angle of the start of injection of the main portion of fuel into the combustion chamber. The rational value of the angle (4 degrees of crankshaft rotation less than that The period of delay in the evaporation of fuel and water should not exceed the time for drops of fuel to reach the outlet from the manifold. Droplets entering the combustion chamber can adversely affect engine performance due to additional cooling of the cylinder walls due to their evaporation. Based on this condition, as well as depending on the value of the maximum cyclic fluid supply, the maximum permissible diameter of the spray holes and the distance from the nozzles to the outlet from the manifold should be determined. Figure 9 shows the influence of the temperature of gases at the inlet to the intake manifold on the concentration and the period of delay in the evaporation of fuel and water. For the engine under consideration, the time for fuel droplets to reach the exit from the manifold is 0.024-0.026 s, and the maximum permissible diameter of the spray holes is 0.6 mm, the time for water droplets to exit the manifold is 0.024-0.030 s, the maximum permissible diameter spray holes 0.7 mm. As can be seen from Figure 9, the time of evaporation of liquid droplets in the considered range of temperatures of gases at the inlet to the intake manifold does not exceed the time for the non-evaporated droplets to reach the outlet from the manifold. When developing an engine, it is important to comply with the limitations on the mechanical loading of its parts, which can be estimated by the value of the maximum gas pressure in the combustion chamber. Figure 10 shows the dependence of the maximum gas pressure in the combustion chamber on the cyclic supply of fuel and water to the intake manifold. Restrictions can be fulfilled by changing the angle of the start of injection of the main portion of fuel into the combustion chamber. The rational value of the angle (4 degrees of crankshaft rotation less than that Using a simulation model, the power dependences (relative to the nominal value for a serial diesel engine) of the engine on the supply of fuel and water to the intake manifold were obtained, shown in Figure 11. of the engine without fuel and water injection into the intake manifold) was justified during the calculation. Using a simulation model, the power dependences (relative to the nominal value for a serial diesel engine) of the engine on the supply of fuel and water to the intake manifold were obtained, shown in Figure 11.  Using a simulation model, the power dependences (relative to the nominal value for a serial diesel engine) of the engine on the supply of fuel and water to the intake manifold were obtained, shown in Figure 11. It should be noted that the dependences were obtained for specific gas and thermodynamic conditions (pressure, temperature and gas velocity in the intake manifold) corresponding to the maximum power mode of a particular diesel engine. For other engine models and modes of operation, the rational values of the design parameters may be different.
In general, the nature of the dependences obtained is expected, consistent with practical experience and available experimental data. The developed mathematical model, described in this article and implemented in a complex engine simulation model, makes it possible to evaluate the quantitative characteristics of these dependencies and, on their basis, to determine the rational design parameters of the fuel and water injection system into the intake manifold. For example, in the framework of this study, it was possible to increase the maximum engine power by 13%, while reducing the specific effective fuel consumption by 2% and ensuring that the maximum pressure of the gases in the combustion chamber is met. It should be noted that the dependences were obtained for specific gas and thermodynamic conditions (pressure, temperature and gas velocity in the intake manifold) corresponding to the maximum power mode of a particular diesel engine. For other engine models and modes of operation, the rational values of the design parameters may be different.

Conclusions
In general, the nature of the dependences obtained is expected, consistent with practical experience and available experimental data. The developed mathematical model, described in this article and implemented in a complex engine simulation model, makes it possible to evaluate the quantitative characteristics of these dependencies and, on their basis, to determine the rational design parameters of the fuel and water injection system into the intake manifold. For example, in the framework of this study, it was possible to increase the maximum engine power by 13%, while reducing the specific effective fuel consumption by 2% and ensuring that the maximum pressure of the gases in the combustion chamber is met.

Conclusions
In the course of a theoretical study, a mathematical model of working processes in the elements of the gas-air duct of a diesel engine with forcing by simultaneously supplying fuel and water to the intake manifold was developed, which: • includes detailed mathematical submodels of the processes of simultaneous evaporation of fuel and water in the intake manifold; • takes into account the influence of the process of evaporation of fuel and water on the parameters of the gas state in the intake manifold; • takes into account the influence of the parameters of the state of the working fluid in the intake manifold on the physical characteristics of fuel and water; • complies with the principles of component modeling, since it does not contain "unnecessary" parameters that are not related to the component being modeled; • describes the process of simultaneous transfer of vapors and non-evaporated liquids between components; • is based on basic equations of physics, and not on empirical dependencies that require data on the dynamics of fuel evaporation under reference conditions.
Together, these features of the mathematical model distinguish it from the well-known ones implemented in the form of software widely used in the engine building industry (AVL Boost, GT-Suite, Ricardo Wave, Siemens Amesim, etc.).
The practical value of the mathematical model created for the engineering and scientific community lies in the fact that it can be used in conjunction with known submodels of processes in systems and mechanisms of a reciprocating internal combustion engine for the synthesis and subsequent use of complex simulation models and digital twins of diesel modifications with forcing by changeable injection of fuel and water into the intake manifold. By contrast with the aforementioned known mathematical models, it allows us with a sufficient degree of reliability to determine the rational design parameters of the fuel and water injection systems, the adjusting parameters of the forced diesel engine, while providing the required power and economic indicators and taking into account the limitations on the magnitude of the mechanical and thermal loads of its parts.
As part of the research, the mathematical model created was introduced into the existing software "Internal combustion engine research and development", designed for the simulation of engines and power plants based on them, previously developed at South Ural State University. Comparison of the results of test calculations and field experiments confirmed the adequacy of the developed mathematical model and software.
In the course of test calculations, general regularities of the influence of various factors on the characteristics of the 12ChN15/18 engine processes with water and fuel injection into the intake manifold were revealed, and rational parameters of the injection system were determined to provide the required power and economic indicators, taking into account the permissible mechanical loads on its parts.