Prognostic Assessment of the Performance Parameters for the Industrial Diesel Engines Operated with Microalgae Oil

: A study conducted on the high-speed diesel engine (bore/stroke: 79.5/95.5 mm; 66 kW) running with microalgae oil (MAO100) and diesel fuel (D100) showed that, based on Wibe parameters ( m and ϕ z ), the difference in numerical values of combustion characteristics was ~10% and, in turn, resulted in close energy efﬁciency indicators ( η i ) for both fuels and the possibility to enhance the NO x -smoke opacity trade-off. A comparative analysis by mathematical modeling of energy and traction characteristics for the universal multi-purpose diesel engine CAT 3512B HB-SC (1200 kW, 1800 min − 1 ) conﬁrmed the earlier assumption: at the regimes of external speed characteristics, the difference in P me and η i for MAO100 and D100 did not exceeded 0.7–2.0% and 2–4%, respectively. With the reﬁnement and development of the interim concept, the model led to the prognostic evaluation of the suitability of MAO100 as fuel for the FPT Industrial Cursor 13 engine (353 kW, 6-cylinders, common-rail) family. For the selected value of the indicated efﬁciency η i = 0.48–0.49, two different combinations of ϕ z and m parameters ( ϕ z = 60–70 degCA, m = 0.5 and ϕ z = 60 degCA, m = 1) may be practically realized to achieve the desirable level of maximum combustion pressure P max = 130–150 bar (at α ~2.0). When switching from diesel to MAO100, it is expected that the η i will drop by 2–3%, however, an existing reserve in P max that comprises 5–7% will open up room for further optimization of energy efﬁciency and emission indicators.


Introduction
Today, the world is challenged with the twin crises of fossil fuel reduction and environmental degradation. Unselective extraction and excessive consumption of fossil fuels have led to a decrease in underground-based carbon capitals. The hunt for alternative energy, which assures a positive correlation with sustainable growth, energy conservation and management, efficiency, and ecological protection, has become extremely marked over the last two decades.
The entire transport sector, including industries providing transportation, agree on the need to decarbonize traffic before 2050-2060 for the most developed and 2060-2080 for less developed economies [1,2]. Usually, this roadmap is primarily associated with the wider deployment of electric transport. It is likely that different energy vectors (CNG, LPG, synfuels, vegetable oils and biodiesel, GTL, H 2 , electricity-see Table 1) will play a role in transport decarbonization [3]. If properly allocated to hard-to-decarbonized modes of transport, synfuels and sustainable biofuels, coupled with the direct electricity consumption through either electrified railways or battery electric vehicles, will all be important in the process of reducing 'carbon intensity' in transport [3]. The European Union member countries alone consume approximately a fourth of the petroleum exploited globally per year. Global consumption of petroleum products has been growing as a Sustainability 2021, 13, 6482 2 of 23 result of the rapid development of Asian economies (China, India) as well. EU authorities have recently started referring to new pollution and climate change control measures more frequently. There is a unanimous consensus within the community on securing long-term clean energy supplies for Europe, in addition to the reduction of greenhouse gas (GHG) emissions from the transport and energy sectors. Lower environmental pollution and higher economic efficiency are probably the biggest advantages of fuel alternatives to gasoline and diesel. However, several studies [4,5] provide readers with the counter argument, that a massive replacement of combustion engine-powered vehicles by battery electric vehicles alone cannot deliver greenhouse gas reductions consistent with climate stabilization and, in the future, may lead to the depletion of key mineral deposits, such as magnesium and lithium. The producers of the high-power industrial diesel engines of low-and average-speeds see fuel flexibility and robustness (VLSFO/MGO, 20% H 2 in NG, biofuels, MeOH, NH3, H 2 as future alternatives) as the key advantage and offer a wide range of sector-specific scenarios, outlining the potential benefits of a particular fuel choice. Table 1. Alternative fuels currently being heavily studied for transport applications.

Types of Fuel Description
Oils and biodiesels (including microalgae oil and biodiesel), BTL (biomass-to-liquid), and alcohol fuel [6][7][8][9][10] Because they are produced in plants that chemically 'fix' or capture carbon dioxide, these types of fuel are characterized for their low environmental pollution. Nonetheless, the production of such fuels requires large amounts of energy if compared with gasoline or diesel fuel.
Gas-to-liquid (GTL) [11][12][13] Used as a substitute for diesel fuel, as GTL leads to a significant reduction in air pollution from internal combustion engines.
Synthetic fuels (or synfuels) [14,15] Production of synthetic diesel fuels obtained from biomass, household waste, and/or natural gas has begun approx. 20 years ago. Synfuels are not considered as alternative fuels since they do not require any modifications in the fueling infrastructure or engine design. However, synfuels expand the raw materials base as well as enhance biodiversity and restore the natural ecological balance due to their easy quality assurance during the production process.
Compressed natural gas (CNG) and liquefied natural gas (LNG) [16] CNG and LNG both are highly functional and efficient type of fuel gas. Theoretically, natural gas resources are vast, if not taking into account the global geopolitics. If we take into consideration the geopolitical situation in the world, specialists believe they may be depleted by 2060. Hence, natural gas is playing a large role in near-future energy prediction. The advantage of CNG/LNG compared with other types of alternative fuel is lower CO 2 emission and a higher heating value (48.7 MJ/kg) in comparison to diesel fuel (42.6 MJ/kg). Currently, nearly all European automobile manufacturers offer natural gas-powered vehicles to the market. Most of them can run both on gasoline and on natural gas, however, bi-fuel engines lead to higher environmental pollution if compared with those running only on natural gas.
Liquefied petroleum gas (LPG) [17] LPG (mainly propane and butane) is prepared by refining raw natural gas or crude oil and is a co-product of the refining process. This type of fuel is highly explosive. Moreover, LPG characteristics are different from those of the diesel fuel, which means that engines designed for both types of fuel are inefficient.
Hydrogen (H 2 ) [18] Hydrogen in gaseous or liquid form may be used in conventional internal combustion engines. This type of fuel carries three times more energy than gasoline; however, density of the former is significantly lower even when compressed. Moreover, a significant amount of electrical energy is required for H 2 generation.
There is a widespread consensus that diesel engine technology has not reached its full maturity and potential yet, in terms of efficiency or lower carbon impact [1]. This is an advanced technology through which synthetic, zero-emission fuels can be produced using only renewable energy and CO 2 [1]. Fuels such as H 2 , GTL, CNG, 3rd generation microalgae fuels, and synthetic hydrocarbons that are made using energy from renewables or other low-carbon energy sources could play a role in multiple hard-to-decarbonize subsectors of global transportation [19]. There is also another key, future-oriented reason that should motivate legislators and OEMs to keep diesel engines in the game: power-to-fuel or power-to-x [1,20]. Synthetic fuels made from carbon dioxide captured from the air or 3 rd generation microalgae fuels made from CO 2 captured from industrial power plants can be successfully used as transportation fuels in conventional engines.
Gaps in the literature, which we are trying to fill. The goal of the overall transport sector is to largely decarbonize and move from 7.7 metric gigatons of emissions per year to 3-2 metric gigatons by mid-century (2050), while ensuring climate resilience. Based on IEA data, predicted global demand for fuel and energy by the transport sector will increase by 140%, 75%, and 70% in aviation, freight transport, and passenger cars, respectively, between 2000 and 2050 [20]. As for the EU transport sector, the agreement was obtained in 2018 that a 14% RES target by 2030, including the gradual phase out of crop-based biofuels from 7% in 2020 to 3.8% in 2030 and a 3.5% share of advanced biofuels of 2nd and 3rd generation. As described previously, all kinds of alternative fuels being heavily studied for transport applications today will all be essential in the process of reducing 'carbon intensity' in transport. The main advantage of the 3rd generation microalgae oil-powered heavy-duty engine over other alternative fuels is that such a vehicle can be relatively called 'CO 2 neutral'-a feature which is characteristic to the very limited variety of fuels of the future (synfuels, power-to-x, etc.) [8]. Notwithstanding the large amount of research studies conducted so far, the majority of the dedicated works are limited to an assessment of diesel engines of a particular modification and presentation of the insights and recommendations obtained for the engine-specific scenario [6][7][8]. This situation usually leads to the disparity between the total quantity of recorded knowledge and the limited capacity of researchers, government bodies, and legislators to assimilate it as well as take action [6][7][8]. This was the main reason behind the prognostic assessment of the industrial diesel engine family for energy efficiency and CO 2 levels and to take a broader look at pure microalgae oil, a potentially carbon-free resource, as a candidate for future transportation energy mix. This study is a part of larger scale investigation conducted by co-authors in the field of transport decarbonization.

Materials and Methods
A simplified scheme, known as an 'experiment planning chart', for this investigation is outlined in Figure 1. Below is a description of the basic steps in this research process. A reliable set of CI engine parameters that could be transferable to the entire family of industrial engines for a possible switch to pure microalgae oil was obtained to examine factors that may contribute to the acceptance of or resistance to this carbon-neutral type of fuel.

Fuels and Test Engines
The methodological part discusses first the 'pilot' study as the first step of the entire research trial, and then the 'main study' [21]. The purpose of the pilot study was to extract a number of parameters (factors) from an 1Z engine map of microalgae oil use coupled with the assessment of smoke opacity and CO 2 levels obtained at various loads and an in depth analysis of the thermodynamic engine cycle simulation. Pure microalgae oil (MAO100), as a less investigated type of fuel, was subjected to CAT 3512B HB-SC and FPT Industrial Cursor 13 engines performance modelling. Conventional diesel fuel (D100) not containing a required 5% biodiesel additive was used as a reference fuel to compare and to contrast the engine indicators obtained for both fuels (Table 2). reliable set of CI engine parameters that could be transferable to the entire family of industrial engines for a possible switch to pure microalgae oil was obtained to examine factors that may contribute to the acceptance of or resistance to this carbon-neutral type of fuel.

Fuels and Test Engines
The methodological part discusses first the 'pilot' study as the first step of the entire research trial, and then the 'main study' [21]. The purpose of the pilot study was to extract a number of parameters (factors) from an 1Z engine map of microalgae oil use coupled with the assessment of smoke opacity and CO2 levels obtained at various loads and an in depth analysis of the thermodynamic engine cycle simulation. Pure microalgae oil (MAO100), as a less investigated type of fuel, was subjected to CAT 3512B HB-SC and FPT Industrial Cursor 13 engines performance modelling. Conventional diesel fuel (D100) not containing a required 5% biodiesel additive was used as a reference fuel to compare and to contrast the engine indicators obtained for both fuels (Table 2). During the 'pilot' study, the sensitivity of microalgae oil to commercial diesel fuel, using the primary and derivative parameters of the bench test engine (type-1Z) that were published by one of the co-authors in 2019 [21], were assessed to build a data matrix, whose rows represent different repetitions of an experiment, and whose columns represent different kinds of data taken for each repetition. A more detailed description of the Obtaining complex input parameters set (primary and derrivative) from 1Z engine performance on microalgae oil + AVL boost modelling of in-cylinder processes from a study published by one of co-authors in 2019 [21].
Model calibration and evaluation of its ability to generate a reliable set of outputs from a supplied set of obtained inputs.
The use of a calibrated model for the assessment of parameters for CAT 3512B HB-SC industrial engine running with microalgae oil. Energy efficiency, process parameters and CO 2 emission levels were simulated and compared with the experimental data (from a study published by one of the co-authors in 2021 [23]) for the same engine that was exploited on diesel fuel (as an interim concept).
With refinement and development of the interim concept, the model led to prognostic evaluation of suitability of microalgae oil as fuel for FPT Industrial Cursor 13 engine family  During the 'pilot' study, the sensitivity of microalgae oil to commercial diesel fuel, using the primary and derivative parameters of the bench test engine (type-1Z) that were published by one of the co-authors in 2019 [21], were assessed to build a data matrix, whose rows represent different repetitions of an experiment, and whose columns represent different kinds of data taken for each repetition. A more detailed description of the specific test procedures, methodologies applied, technical characteristics of the measurement equipment, and uncertainties of the measured parameters is presented in Ref. [21]. The derivative parameters necessary to build a calculation matrix are as follows: P e , p me , BSFC, net efficiency, mass air flow rate, exhaust gas temperature, air excess ratio, etc. They were established at three different loads and with the advancement or retardation of the start of fuel injection timing every 2 degrees of the crankshaft rotation angle from −16 • to +4 • relative to TDC. To obtain a full dataset for the mathematical modelling of thermochemical processes, the following derivative parameters were also assessed: incylinder pressure, pressure rise rate, in-cylinder temperature, temperature rise rate, and heat release rate. Technical specifications of the engines subjected to the 'pilot' and 'main' studies are presented in Table 3, while Figure 2 depicts the schemes of the discussed above type of engine. start of fuel injection timing every 2 degrees of the crankshaft rotation angle from −16° to +4° relative to TDC. To obtain a full dataset for the mathematical modelling of thermochemical processes, the following derivative parameters were also assessed: in-cylinder pressure, pressure rise rate, in-cylinder temperature, temperature rise rate, and heat release rate. Technical specifications of the engines subjected to the 'pilot' and 'main' studies are presented in Table 3, while Figure 2 depicts the schemes of the discussed above type of engine.

IMPULS Modelling
Mathematical model (MM) calibration by standard techniques required a primary (historical) database (an outcome of a pilot study) for key variables and the primary techniques of statistical process control. MM calibration was mainly carried out manually using many adjustments of the engine parameters obtained during its operation on microalgae oil to improve the model (both natural state and history matching).
A set of inputs used for adjusting the parameters of the 'Wiebe' [22] combustion function included in 1-D thermodynamic engine cycle simulation tool IMPULS, developed at the Central Diesel Research Institute (CNIDI) [23] was analyzed. The structure of this software is constantly being improved with the supplementation of sub-models intended to assess the formation and combustion of the fuel-air mixture, the dynamics of fuel injec-

IMPULS Modelling
Mathematical model (MM) calibration by standard techniques required a primary (historical) database (an outcome of a pilot study) for key variables and the primary techniques of statistical process control. MM calibration was mainly carried out manually using many adjustments of the engine parameters obtained during its operation on microalgae oil to improve the model (both natural state and history matching).
A set of inputs used for adjusting the parameters of the 'Wiebe' [22] combustion function included in 1-D thermodynamic engine cycle simulation tool IMPULS, developed at the Central Diesel Research Institute (CNIDI) [23] was analyzed. The structure of this software is constantly being improved with the supplementation of sub-models intended to assess the formation and combustion of the fuel-air mixture, the dynamics of fuel injection, evaporation, flame spread, and the possibility to evaluate fuels with different chemical composition, etc. Most of the phenomenological sub-models implemented in this program are similar in their nature with the other widely used and adopted software-AVL BOOST [24]. The first law of thermodynamics for closed systems in the form of differential equation is formulated for the CI engine as follows: Therefore, dU considers the change in internal energy in the system (J); dQ re quantifies the combustion heat released (J); dQ e calculates the energy exchange (wall heat transfer from the cylinder gas) (J); pdV is the volumetric work (p-pressure (Pa); V-volume (m 3 )); h·dm sum up all enthalpy flows (index ex over the exhaust valves, s over intake valves); m s is the supply (intake) air mass (kg); m ex is the mass of exhaust gas (kg); and τ is the time (s). As described in Ref. [22], depending on the flow direction, a negative sign indicates that the enthalpy leaves the system, and a positive sign indicates an entering of enthalpy flow. The mass balance equation (Equation (2)) and state equation (Equation (3)) are as follows [25]: where: m inj is the mass of injected fuel (kg); R is the gas constant (J/kg·K); and T is the temperature (K). The rate of heat release according to the 'Wiebe' model [22,26] was determined by the following equation: Thereby, Q describes the total heat input; τ is the angle between initial and current time; C is a function parameter that is equal to 6.9 for the case of complete combustion; and τ z is the relative time of combustion.
The fraction of heat released since the start of combustion can be assessed by using Equation (5) [26]: The heat release was determined by the 'Wiebe' model, and parameters were then converted to partial load modes by using the empirical correlations by Woschni [22,27] in the form of a simplified relationship between the heat release parameters and operational parameters of an engine (fuel/air ratio, injection timing, crankshaft rotational speed, charge in air pressure, charge in air temperature, etc.). A simplified 'Wiebe' heat release diagram is characterized by two parameters: m-form factor and ϕ z -conditional duration of combustion [23,27].
This set of equations is solved by using the selected method, which is programmed in the computer code together with the boundary conditions [26,28]. When the iterative process was finished, the validation of the measured and predicted results from the IMPULS model is performed. Experimental results in the form of datasets of a 1Z engine running with (i) microalgae oil and (ii) diesel fuel that were obtained during the pilot study were compared and contrasted with the ones derived by IMPULS software. The following parameters were compared at p me = 12.6, 18.9, and 25.1 bar: air pressure after compression (p k ), air temperature after compression (T K ), pressure of compression in the cylinder (p c ), maximum cycle pressure (combustion pressure) (p max ), exhaust gas temperature (T g ), hourly air supply (G air ), excess air coefficient (α), in-cylinder pressure increase rate (λ), indicated thermal efficiency (η i ), and effective efficiency (η e ). A full description of each interim step of MM calibration to assess the 1Z engine running on pure diesel fuel is described in detail in Ref. [23]. Later on, MM suitability for accurate assessment of the complete cycle of combustion for a 1Z engine powered by MAO100 was investigated assuming that ϕ inj = 2 CAD BTDC, n = 2000 min −1 , BMEP = 4, 6 and 8 bar. After checking the overlap between the experimental and modeling outcomes which show good agreement between the results, we made a prerequisite for the next step in this study-to transfer and to adopt the successfully established heat release characteristics for the 1Z engine of a passenger car to a heavy-duty industrial diesel engine equipped with a similar fuel injection system. Hence, the CAT 3512B HB-SC engine was assessed for six different speeds (n = 1000, 1200, 1300, 1400, 1500, 1600, and 1800 min −1 ). Variation in 10 parameters was analyzed for high load mode: brake effective mean pressure, indicated thermal efficiency, effective efficiency, brake-specific fuel consumption, exhaust gas temperature, excess air Sustainability 2021, 13, 6482 7 of 23 coefficient, maximum cycle pressure (combustion pressure), pressure after compression, mechanical efficiency, and the change in CO 2 emissions.
The potential of this work extends beyond the creation of engine maps to allow investigations into the transferability of heat release characteristics from a passenger car diesel engine to the industrial one; this study offers a guideline as a modeling outcome for prognostic assessment of the engine parameters obtained for the entire family of engines [29,30]. Finally, modeling of the combustion process for the CI engine (see Table 3), representing the entire FPT Industrial Cursor 13 engine family has been performed to evaluate the compatibility of microalgae oil with a wide range of industrial CI engines. This part of the research was dedicated to further interpretation of engine operational parameters (combustion pressure, indicated efficiency, hourly fuel consumption, maximum combustion pressure, exhaust gas temperature, and CO 2 emissions) in order to describe the character of heat release according to various combinations of the form factor (m) and combustion duration (ϕ z ). A relative change in the aforementioned variables when switching an engine from diesel fuel to microalgae oil was established, as well as the zones of rational combination of m and ϕ z for each specific parameter were identified. The fitness of this model for the entire family of industrial engines was demonstrated by the successful adoption of data sets obtained during the pilot study and their successful transfer creating candidate engine maps and identifying the zones of rational combination of m and ϕ z parameters. On the basis of the tests and modelling conducted, the IMPULS model, due to the described above peculiarities of the algorithm realization, seems a potentially useful tool for compatibility analysis of various size diesel engines to run on microalgae oil with the described above physico-chemical properties (see Table 2).

Selection of Strategies to Improve NO x -Smoke Opacity Trade-Off
The findings of a pilot study showed that, although the compression ignition engines are designed for optimum operation with fossil fuel, further advancements in lowering exhaust gas emissions are possible through the use of CO 2 -neutral fuel with favorable fuel characteristics and the proper adjustment of high reaction fuel injection time (ϕ inj ). Figure 3 shows the direction of possible NO x -smoke opacity improvement for a 1Z engine. The abscissa denotes the specific emission of nitrogen compounds and the ordinate denotes the specific smoke opacity levels. The NO x -smoke opacity trade-offs were assessed for the engine running at average load conditions, with the advancement or retardation of the start of fuel injection timing every 2 degrees of the crankshaft rotation angle (CA) from −2 to +16 degCA relative to the top dead center (TDC). For the examined 1Z engine, the optimum approach was to retard the main injection timing (by 2 degCA) [31]. In that case, in addition to the reduction in smoke opacity levels predetermined by the simple replacement of diesel fuel with microalgae oil, further emission reductions were achieved leading to 41.2% (4.3 g/kWh vs. 2.5 g/kWh) lower engineout NOx emissions for both unary fuels. Moreover, the interval of −2… 0 degCA can be described as the best setting of an engine for smoke and NOx stabilization and reduction: the first part of the interval related to 0 degCA demonstrated the same (and the second lowest) levels of NOx emissions for both tested fuels, while the second one showed a 22.2% and 2.0% decrease in smoke opacity levels for diesel fuel and MAO100, respectively. It is worth noting that, initially, the difference between smoke emission levels at 0 degCA for D100 and MAO100 was 49%, while the retardation of the SOI by two degrees relative to TDC reduced this difference to 29% (see Figure 3). The findings presented in Figure 3 coincide well with the published works [32][33][34], which reported that one of the most significant CI engine parameters affecting the NOx emissions is retarding injection timing. As far as eight different advanced injection timings at 2000 min −1 are concerned (from +2 to +16 degCA), the earlier was the injection timing, the lower were the NOx emissions due to lower in-cylinder temperature and heat release rate peak of the injection. The difference between the NOx levels obtained at +2 degCA  As shown in Figure 3, the original high values of smoke opacity in the case of diesel fuel operation are reduced by two-fold by replacing it with microalgae oil. As microalgae oil produces lower smoke opacity levels compared with D100 at the same engine control setting, there is a big potential for nitrogen oxide emission reduction through a proper engine readjustment aiming to significantly improve the NO x -smoke opacity trade-off. For the examined 1Z engine, the optimum approach was to retard the main injection timing (by 2 degCA) [31]. In that case, in addition to the reduction in smoke opacity levels predetermined by the simple replacement of diesel fuel with microalgae oil, further emission reductions were achieved leading to 41.2% (4.3 g/kWh vs. 2.5 g/kWh) lower engine-out NO x emissions for both unary fuels. Moreover, the interval of −2 . . . 0 degCA can be described as the best setting of an engine for smoke and NO x stabilization and reduction: the first part of the interval related to 0 degCA demonstrated the same (and the second lowest) levels of NO x emissions for both tested fuels, while the second one showed a 22.2% and 2.0% decrease in smoke opacity levels for diesel fuel and MAO100, respectively. It is worth noting that, initially, the difference between smoke emission levels at 0 degCA for D100 and MAO100 was 49%, while the retardation of the SOI by two degrees relative to TDC reduced this difference to 29% (see Figure 3).
The findings presented in Figure 3 coincide well with the published works [32][33][34], which reported that one of the most significant CI engine parameters affecting the NO x emissions is retarding injection timing. As far as eight different advanced injection timings at 2000 min −1 are concerned (from +2 to +16 degCA), the earlier was the injection timing, the lower were the NO x emissions due to lower in-cylinder temperature and heat release rate peak of the injection. The difference between the NO x levels obtained at +2 degCA and +16 degCA was 14 g/kWh. The results of a pilot study revealed that for the particular type of test engine with gradually retarded degrees of SOI (start of injection), MAO100 is in many cases comparable to most of the critical engine parameters and indicators, compared with D100 [21]. For some of the parameters (air mass flow, pressure provided by the turbocharger, air-fuel equivalence ratio, the level of nitrogen oxides, and slightly reduced oxygen content in the exhaust gases), depending on the load applied, the tests indicated a very small difference in measured values, that fell within the measurement error. In the form of increased levels of emission, negative trends were observed only for HC and CO 2 [21]. In general, the decreasing trends of soot concentration in the exhaust gases of microalgae oil were mostly associated with better carbon-oxygen balance, leading to an improvement in the combustion reaction and better promotion of soot oxidation, especially at the end of the cycle. The increasing trend of CO 2 gases for microalgae oil indicates the higher fuel consumption rates for MAO100 being predetermined by its lower heating value compared with diesel fuel. Another strong argument for a wider adoption of microalgae oil in diesel engines lies in the fact that when the test engine was powered by pure P.moriformis oil, the indicated thermal efficiency (η i ) was approximately equal to D100-this trend was inversely proportional to the BSFC character. The use of MAO100 led to the slightly higher, but falling within the margin of error, indicated thermal efficiency (0.355) in comparison to diesel fuel (0.350) at high load (0.8 MPa), followed by the similarly close values of η i at 0.6 MPa (0.350 vs. 0.345) and 0.4 MPa (0.325 vs. 0.320) loads [35]. This allows a wider window for energy and emission enhancement for industrial diesel engines using microalgae oil as fuel.

Model Calibration Outcomes
1-D predictive engine model was developed and calibrated. Developing MM for better accuracy needs a lot of data related to component characteristics and computational power. Data requirements vary with fuel characteristics, engine architecture, and depending on how detailed the model being developed was, based on the requirements and information available [36]. According to the engine IMPULS simulation model and the acquired operating points of 1Z engine running with pure diesel fuel and pure microalgae oil, the engine input and output performance database was obtained, including 16 variables juxtaposed for three different loads (see Table 4). The calibration of MM for the 1Z engine running on pure diesel fuel is described in detail in Ref. [23]. Measured in-cylinder pressure is used to calculate the rate of heat release (see Figure 4). Variables such as G air , P k , T k , ϕ inj , and others were used to run IMPULS. 1-D predictive engine model also calculated an in-cylinder pressure; however, since that was a tool that also needs calibration, only the measured in-cylinder pressure was used. Hence, a calibration of constants for various form factor (m) and conditional duration of combustion (ϕ z ) values was performed to match the measurement data as accurately as possible. As described in Ref. [37], an automated calibration script was used for a set of 1Z operation points together with a range for the model constants. The regiment evaluated the output from IMPULS and compared it with measurement datasets. The result of this routine is provided as a new dataset of constants for further iterations. The whole procedure was enabled to stop automatically when the regiment meets one of its stopping criteria. The main output of MM calibration is a derivation of the optimum constants with the smallest difference between the measured and predicted results [37]. According to Stenberg (2008) [37], a set of constants are applied to the validation points to see the performance of MM for all feasible regions of an optimization problem rather than those used for the training of the MM.     Whether the established model is accurate or not could be evaluated by the fitting degree between the modeling data and the experimental dataset. An obtained model accuracy for various parameters is as follows: p me , 0-4.3%; p K , 2.5-4.5%; T K , 0.8-1.7%; p c , 0.6-1.7%; p max , 1.6-3.9%; T g , 1.9-3.2%; α, 5.1-10.0%; λ, 0-3.9%; η e , 0.3-3.4%; and g cycl , 0-1.7%. In general, the biggest difference between the measured and predicted parameters did not exceed 3.5% on average. The modelling results of the same parameters for average and high loads have shown a close to 100% percent coincidence with the experimentally obtained ones.
As described earlier, we simultaneously collected in-cylinder pressure data for three different loads coupled with the datasets of pressure provided by the turbocharger. Firstly, the variation in in-cylinder pressure was determined experimentally (see Figure 2) when the engine was running with both unary fuels, averaged over 10.000-13.300 cycles per each load setting. Then, the MM was calibrated for the reference point (75% load) of an engine running on diesel fuel [23] and microalgae oil using also some data from the manufacturer to select the appropriate maps for the AC compressor AFN TH 614 and turbine CHRA 454082 of the turbocharger due to the lack of their data. They are rescaled from existing maps to adjust the corrected mass flow rate and the corrected speed with the given pressure ratio exported from the in-cylinder pressure diagram to fit the current 1-D model [38]. This interim undertaking is important in calibrating the outlet pressure of the compressor by fitting the computed pressure diagram with the experimental one along the compression process. It also helps to verify the experimentally established amount of air and thus the volumetric efficiency against the modelling outcomes. This process is repeated for the other two cases presented in Figure 4 (0.6 MPa (50% load) and 0.4 MPa (25% load) for better calibration of the maps of the turbocharger and for improving the accuracy of the computed results along the different cases. The initial conditions of atmospheric pressure and temperature are adapted for each of the three loads, according to the pilot study data available for better accuracy [38].  As previously described in Ref. [38], to minimize the mean absolute percentage deviation, we used a derivative version of the goodness-of-fit method that describes the correlation between the in-cylinder pressure diagram of a 1Z engine and the corresponding data from the pilot study along the combustion process, by finding the optimal values of injection system parameters due to the lack of their availability. Regarding the validation of MM, Figure 4 displays the comparison of experiments and numerical results for three load points of a 1Z engine running on pure microalgae oil. At the high load (0.8 MPa) operating point where the magnitudes of the brake mean effective pressure are greater, the values obtained during the pilot study and MM are within 1.3% of each other, and at the lower load points where the magnitudes of the BMEP are much smaller, the values are within 2.1% of each other for P me = 0.6 MPa and 1.6% for P me = 0.4 MPa (see Figure 4). The plots indicate a very good agreement between the measured and predicted values for brake mean effective pressure. It is worth nothing that the range of calculation errors (1.3-2.1%) was mainly characteristic to these parts of the curves, which show the pressure peaks exceptionally. Again, this is due to the magnitude of the BMEP values. In terms of peak pressure, the error remains lower than +/−1.5 bar and can be considered as acceptable. Overall, this work shows a methodology able to establish some parameters for optimum engine operation in high-pressure direct injection combustion mode. This was achieved through the analysis of combustion chamber pressure (Figure 4a-c) and heat release rate (Figure 4d,e). The resulting plots of heat release have high precision, leading to a reliable estimation of the process, that is critical to a good modeling of the selected types of diesel engines. The insignificant differences were obtained only around the piston position corresponding to the top dead center, which in turn was predetermined by the application of the single-phase vibe combustion model for diesel engines. In other words, microalgae oil performance was accurately approximated by a single linear segment, which means a single-phase form. The differences between the modelling and experimental outcomes of the heat release characteristics on average did not exceeded 2%, with the exception of dx/dϕ peaks that show the 5-7% discrepancy in values. Figure 4 illustrates that the simulation data of the cylinder pressure and the heat release rate were in good agreement with the pilot study results; thus, the model was reasonable.
The changes in heat release characteristics HRR (dQ/dϕ) when switching the engine from diesel fuel to microalgae oil is shown in Figure 5. Commonly in compression ignition engines, the four stages in the combustion process are as follows: ignition delay period, premixed combustion phase (all fuel has been injected, and the pressure increases rapidly), mixing controlled combustion phase (the fuel is burned and produces power), and late combustion phase (the pressure is going down) [35].
In this section, the impact of microalgae oil and diesel fuel on engine combustion characteristics is discussed in terms of HRR. Depending on the research aims, the characterization of the ongoing processes within engine cylinders are calculated using the experimentally obtained in-cylinder process. These are the heat release characteristics dx dϕ = f (ϕ), expressed by differential equations, and the parameters of the in-cylinder processes, such as ignition delay, factor cycle dynamics (σ, λ, dp/dϕ max ), characteristic temperature and pressure (T max , P max ), and others [39]. The dx/dϕ is considered as the basic parameter describing heat release during the process of fuel combustion within the cylinders. The reliability of the results of mathematical modelling of the diesel engine parameters depends on the appropriate setting of dx/dϕ [39]. HRR is affected by the amount of air-fuel ratio, SOI, and the thermodynamic properties inside the cylinder. The premixed phase of combustion ranges between the start of combustion (362 degCA for MAO100 and 363 degCA for D100) and the CA corresponding to the maximum value of heat flux (367 degCA for D100 and 374 degCA for MAO100) [21]. The early stage of the premixed combustion phase is associated with the rapid increase in in-cylinder pressure showing the highest peaks of HRR for diesel fuel at P me = 0.  (Figure 5c,d). During the assessment of heat release characteristics for the selected range of loads ( Figure 5), the m and ϕ z parameters demonstrated, on average, 10% higher values in the case of microalgae oil if compared with D100 as the reference fuel. The explanation for this lies in a fact that, during the early stage of mixing controlled combustion phase, a smaller portion of MAO was prepared for combustion during the ignition delay period due to the lower LHV, which in turn reduced the HRR. The values of m and ϕ z will be applied for the modelling of CAT 3512B HB-SC engine as described in Section 3.3. The combustion cycle of a diesel engine: juxtaposition of pilot study outcomes vs mathematical modelling (ϕ inj = 2 degCA BTDC, n = 2000 min −1 ). the exception of dx/dφ peaks that show the 5-7% discrepancy in values. Figure 4 illustrates that the simulation data of the cylinder pressure and the heat release rate were in good agreement with the pilot study results; thus, the model was reasonable.
The progress of the engine differential experimental values HRR ( / ) when switching the engine from diesel fuel to microalgae oil is shown in Figure 5. Commonly in compression ignition engines, the four stages in the combustion process are as follows: ignition delay period, premixed combustion phase (all fuel has been injected, and the pressure increases rapidly), mixing controlled combustion phase (the fuel is burned and produces power), and late combustion phase (the pressure is going down) [35]. In this section, the impact of microalgae oil and diesel fuel on engine combustion characteristics is discussed in terms of HRR. Depending on the research aims, the characterization of the ongoing processes within engine cylinders are calculated using the experimentally obtained in-cylinder process. These are the heat release characteristics = ( ), expressed by differential equations, and the parameters of the in-cylinder processes, such as ignition delay, factor cycle dynamics (σ, λ, dp/dφmax), characteristic temperature and pressure (Tmax, Pmax), and others [39]. The dx/dφ is considered as the basic parameter describing heat release during the process of fuel combustion within the cylinders. The

The Use of a Calibrated MM for the Assessment of CAT 3512B HB-SC Industrial Engine Running on Microalgae Oil
One of the most efficient methods to examine the critical challenges in the internal combustion engine research, when diesel engines of various types and modifications are subject of successful adoption of alternative and sustainable fuels, is the application of modeling software, including IMPULS, in which important indicators of the HRR (m and ϕ z ) were determined according to logarithmic anamorphosis method [40,41] and modified methodology by Bulaty and Glanzman [42]. Meanwhile, the successfully established heat release characteristics for the 1Z engine of a passenger car were successfully adopted to the industrial diesel engine CAT 3512B HB-SC, equipped with a similar fuel injection system. A list of variables experimentally obtained for the 1Z load close to the nominal regime, served as a prerequisite for MM to assess the brake effective mean pressure, indicated thermal efficiency, effective efficiency, brake-specific fuel consumption, exhaust gas temperature, excess air coefficient, maximum cycle pressure (combustion pressure), pressure after compression, mechanical efficiency, and the change in CO 2 emissions at six different speeds (n = 1000, 1200, 1300, 1400, 1500, 1600, and 1800 min −1 ) of a CAT 3512B HB-SC engine (see Table 5). For the simulation of CAT 3512B HB-SC by using the statistical data of an internal combustion engine, the datasets of p max (maximum combustion pressure) were assumed according to the nominal regime of an engine (n = 1800 min −1 ) running with diesel fuel as well as the threshold value of the exhaust gas temperature (before entering the turbine) was set to T g ≤ 973 K (700 • C) and served as a limiting parameter. In other words, MM calibration for all obtained parameters was performed at the nominal regime of a CAT engine, while, for the other seven speeds, the calibrating parameters were not further adjusted. Instead of this, the obtained parameters were juxtaposed with the calculated threshold values of P max and T g to obtain a clear view when it is necessary to limit the volume of g cikl as a way to prevent a further increase in process parameters if their threshold value was already achieved. The following indicative parameters were obtained for diesel fuel (I) and microalgae oil (II) at two different injection settings of an engine (Table 5). For engine operation with diesel fuel, we employed m and ϕ parameters that were established through the use of the modified methodology by Bulaty and Glanzman [42], while, for MAO100, these parameters were assumed to be 10% higher compared with D100. To prevent the rise of the exhaust gas temperature of an engine running with microalgae oil above the 973 K level (see Table 5), the third (III) mode of injection was a subject of IMPULS modeling, where the cyclic portion (g cycl ) of microalgae oil was reduced while other parameters (m and ϕ z ) remained unchanged. For the in-cylinder pressure profiles (P me ) obtained at each selected speed, we apply the torque reserve factor of 1.1-a value that is typical to industrial diesel engines. A precondition of g cycl × LHV = idem D100 has been met, which is based on the already discussed above impermissible level of the exhaust gas temperature (T g ) caused by the prolonged fuel burning process and, especially, its higher intensity during the second phase of combustion at n = 1000-1300 min −1 (see Table 5). The juxtaposition of (I) and (III) outcomes shows that, notwithstanding the fact that T g reached the desired level at n = 1000-1300 min −1 in the case of MAO100 (see Table 5), a 4.0-4.5% reduction in engine torque/brake mean effective pressure (p me ) was observed in the range of crankshaft speeds starting from its nominal value which corresponds to the OEM settings up to the speed that generates the maximum torque, followed by 7.1-15.3% decrease starting from n = 1300 min −1 (M t max ) and further on until n = 1000 min −1 (see Table 5). At the same time, brake-specific fuel consumption of MAO100 increased on average by 19-21% compared with diesel fuel. This phenomenon was predetermined by the difference between LHV values for pure diesel fuel (42.8 MJ/kg) and MAO100 (36.8 MJ/kg) that comprises 13.8% and making it necessary to increase the volume of g cikl to compensate the lower energy content of microalgae oil. These numbers suggest that practically it is essential to ensure that 15.96% larger amounts of fuel are delivered to each engine cylinder when microalgae oil is burnt [35]. Another reason for the exceeded fuel consumption lies in the decline of the effective pressure (η e ) in the speed range of n ≤ 1300 min −1 (due to restriction measures for exhaust gas temperature). Nevertheless, the precondition for the induction period ϕ i = idem was applied for both D100 and MAO100, microalgae oil was responsible for a longer duration of the combustion process that led to the appearance of the reserve next limiting factor-P max (maximum cycle pressure)-showing a decline on average by 9.0-10.5 bar along the range of higher speeds (1800-1400 min −1 ), and by 14.5-20.0 bar for the range of lower speeds (1300-1200 min −1 ). Two influencing factors had a direct impact on the obtained results: (i) adjustment of a combustion process (the uniform duration of the induction phase ϕ i led to a longer conditional combustion duration ϕ z in the case of MAO100), and (ii) forced reduction in g cikl to prevent further increase in T g above the threshold value of 973 K. If compared with diesel fuel, CO 2 emissions of the engine running with microalgae oil (full life-cycle carbon dioxide emissions are absent from analysis) have shown an increasing trend (by 4.0-4.5%) for the higher speed range of an engine (1800-1400 min −1 ) and, on the contrary, for the lower speed range of 1300-1200 min −1 , the 4.5-8.5% reduction was observed.
From the simulation outcomes of the CAT 3512B HB-SC engine running with MAO100, which clearly shows the untapped margin of maximum permissible values of the restrictive criterion P max (see Table 5), we found reasonable grounds for an improvement of the dynamic characteristics of ICE via further advancement of SOI as depicted in Figure 3 (Section 3.1). Hence, the MM outcome using the (IV) injection setting for CAT 3512B HB-SC engine running with pure microalgae oil (T g ≤ 973 K, ϕ inj = 2 degCA BTDC) was juxtaposed against (I) diesel fuel and presented in Table 6. Figure 6 presents the variation in break mean effective pressure, maximum cycle pressure (combustion pressure), indicated thermal efficiency, and exhaust gas temperature for the entire range of engine speeds (1800-1000 min −1 ) for the I (D100), II (MAO100), and IV (MAO100) injection strategies. Results revealed that the traction characteristics of a CAT 3512B HB-SC engine were significantly improved after the last correction of the injection settings. In comparing diesel fuel (I) and microalgae oil (IV), we found that the difference in indicated thermal efficiency now comprises only 0.7-2.0% (see Table 6) and the threshold value of 973 K for exhaust gas temperature was successfully reduced below the threshold level or, in the case of n = 1300 min −1 , an obtained temperature level felt within the margin of experiment error (see Table 6 and Figure 6).
The mathematical modeling results of P me , P max , η i , and T ex for the universal multipurpose diesel engine CAT 3512B HB-SC, which are reported in Figure 6, generally confirm the benefit of using the IV-th injection setting for pure microalgae oil (T g ≤ 973 K, ϕ inj = 2 degCA BTDC) over the II-nd one. A careful look at variation in P me and P max under the range of higher speeds (1800-1400 min −1 ) have shown the increase in these parameters for the case IV on average by 1.6-2.2% and 1.4-6.7%, respectively (Figure 6a,c). The range of lower speeds (1300-1200 min −1 ) was responsible for in 0.5-1.8% lower P me rates and in 6.0-10.8% higher P max values. It was found that at the regimes of external speed characteristics, the difference in P max and P me for MAO100 (IV) and D100 (I) did not exceed 0.7-2.0% and 2-4%, respectively. The exhaust gas temperature was reduced by 8-12 K within the range of engine speeds starting from n = 1400 m −1 to n = 1800 m −1 . A drastic reduction in T g by 18-40 K was observed for the IV-th setting of injection at n = 1300-1000 m −1 (Figure 6d), which in turn enabled us to suppress the further rise of temperature above the threshold value of T g ≤ 973 K and to diminish its negative effects. Figure 6d shows the variation in indicated thermal efficiency against the entire range of engine speeds for the I, II, and IV injection strategies. If comparing the modelling outcomes for the II and IV injection strategies, the later demonstrated 0.4-2.3% higher η i values, thus being an indicator of better fuel conversion efficiency achieved through the advancement of SOI.
The modelling outcomes ( Figure 6, Table 6) revealed that the IV-th injection strategy applied for microalgae oil is more or less equally sensitive to certain engine parameters (P me , P max , η i , η e , α, p k , η m , and T ex ) compared with diesel fuel (I) as well as prevents temperature rise and BMEP losses-negative trends emerged from the modelling of the II-nd and III-rd scenarios. This makes a serious prerequisite for a more rigorous evaluation of the suitability of P. moriformis microalgae oil as fuel for FPT Industrial Cursor 13 engine family, thus bringing a potentially carbon-free resource into the mix, in the form of renewable fuel for the broader range of CI engines. The mathematical modeling results of Pme, Pmax, ηi, and Tex for the universal multipurpose diesel engine CAT 3512B HB-SC, which are reported in Figure 6, generally confirm the benefit of using the IV-th injection setting for pure microalgae oil (Tg ≤ 973 K, φinj = 2 degCA BTDC) over the II-nd one. A careful look at variation in Pme and Pmax under the range of higher speeds (1800-1400 min −1 ) have shown the increase in these parameters for the case IV on average by 1.6-2.2% and 1.4-6.7%, respectively (Figure 6a,c). The range of lower speeds (1300-1200 min −1 ) was responsible for in 0.5-1.8% lower Pme rates and in 6.0-10.8% higher Pmax values. It was found that at the regimes of external speed characteristics, the difference in Pmax and Pme for MAO100 (IV) and D100 (I) did not exceed 0.7-2.0% and 2-4%, respectively. The exhaust gas temperature was reduced by 8-12 K within the range of engine speeds starting from n = 1400 m −1 to n = 1800 m −1 . A drastic reduction in Tg by 18-40 K was observed for the IV-th setting of injection at n = 1300-1000 m −1 (Figure 6d), which in turn enabled us to suppress the further rise of temperature above the threshold value of Tg ≤ 973 K and to diminish its negative effects. Figure 6d shows the variation in indicated thermal efficiency against the entire range of engine speeds for the I, II, and IV injection strategies. If comparing the modelling outcomes for the II and IV injection strategies, the later demonstrated 0.4-2.3% higher ηi values, thus being an indicator of better fuel conversion efficiency achieved through the advancement of SOI. The modelling outcomes ( Figure 6, Table 6) revealed that the IV-th injection strategy applied for microalgae oil is more or less equally sensitive to certain engine parameters (Pme, Pmax, ηi, ηe, α, pk, ηm, and Tex) compared with diesel fuel (I) as well as prevents temperature rise and BMEP losses-negative trends emerged from the modelling of the II-nd and III-rd scenarios. This makes a serious prerequisite for a more rigorous evaluation of the

Prognostic Assessment of Industrial Diesel Engine Family for Energy Efficiency and CO 2 Levels
This work sought to further understand the engine efficiency and prospects for improvement in heavy-duty diesel engines having a comparatively high compression ratio (and brake mean effective pressure, P me ) and broad field of application-the necessary prerequisites for facilitating an uptake and ensuring wider deployment of CO 2 -neutral fuel, pure microalgae oil P. moriformis. The research approach involved the modelling and analysis of a heavy-duty diesel engine, representing the entire FPT Industrial Cursor 13 engine family, to evaluate the compatibility of microalgae oil with a wide range of industrial CI engines. The FPT Industrial Cursor 13 family engines comply with statutory emission standards since they have an optimized combustion process with quick and complete combustion that allows high levels of efficiency even with a lean gas mixture. The performance of the motors ranges from 180 to 420 kW and these power units are broadly applied for mobile (tractors, forestry machinery) and stationary applications (gen-sets). Engines of this family can be operated on diesel fuel and on natural gas. For the tuning of parameters of the Wibe combustion model included in 1-D thermodynamic engine cycle simulation tool IMPULS, we selected a water-cooled six-cylinder direct injection diesel engine of 353 kW capacity at rated speed n = 1900 min −1 . The two primary reference engines tested were a 1Z model 1.9-L passenger car diesel engine, a representative engine for Volkswagen Group, and a CAT 3512B HD-SC model 51.8-L heavy-duty engine, representative of Caterpillar. A simulation study of engine performance taking into account different ratios of heat release parameters m and ϕ z was performed (m = 0-1.5, ϕ z = 50-80 degCA). To determine the energy efficiency indicator η i of a CAT 3512B HD-SC running on diesel fuel and microalgae oil (p me = idem) in a given load mode, it is sufficient to evaluate the air excess coefficient (α) that has a direct influence on the duration of ϕ z . Air-access coefficient values (α = 1.8-2.2) were taken from a range, which corresponds to the values provided in a technical engine specification. Model calibration and validation during the initial stages of this study enabled us to reach the difference threshold of 2-3% between the simulation outcomes of the main energy indicators (η i , η e ) and the same parameters provided in the technical documentation from OEM. The outcomes of the prognostic assessment, in a form of relative change in the main FPT Industrial Cursor 13 engine family parameters (changes (δ) of: indicated thermal efficiency (η i ), hourly fuel consumption (G f ), maximum cycle pressure (P max ), exhaust gas temperature (T g ), and CO 2 emissions) were presented in Table 7 and Figure 7. Table 7 contains the calculation outcomes of η i , G f , P max , and T g for different ϕ z and m combinations (48 combinations, 16 points of each parameter per single chart) as well as the relative changes of these parameters when diesel fuel and microalgae oil are compared. The dataset obtained for an engine's indication of the thermal efficiency when diesel fuel was used (see Table 7) served as a mean to build a reference chart (Figure 7a) to compare and to contrast the changes of other parameters when the engine was switched to run with MAO100 (Figure 7b-f).  Combination of these parameters in similar proportions is characteristic to industrial heavy duty diesel engines equipped with an accumulator fuel-injection system common rail [23]. If comparing the changes in indicated efficiency (δ ) for the reference fuel and alternative fuel, the differences obtained at m = 0.5 did not exceed 2-3% (Figure 7b).  Juxtaposition of D100 (ϕ inj = 0 degCA TDC) vs. MAO100 (ϕ inj = 0 degCA TDC, ϕ z is in 10% ↑ compared to D100) The obtained results revealed that, if considering the smallest changes in the indicated thermal efficiency values as an outcome of the best compatibility of m and ϕ z parameters, this indicator can be characterized by a relatively short period of heat release (50-60 degCA) and moderate dynamics (m = 0-0.5).
Combination of these parameters in similar proportions is characteristic to industrial heavy duty diesel engines equipped with an accumulator fuel-injection system common rail [23]. If comparing the changes in indicated efficiency (δη i ) for the reference fuel and alternative fuel, the differences obtained at m = 0.5 did not exceed 2-3% (Figure 7b).
Further increasing the time of conditional combustion duration (ϕ z ) up to 70-80 degCA, a range that is characteristic to many sorts of conventional injection systems in industrial CI engines, also gives a qualitatively positive non-changing impact of m factor on η i . However, quantitatively, the change in the indicated efficiency values reaches 6-7% ( Figure 7b); this forewarns about the necessity of advancing the engine's timing (ϕ inj ) in order for the injection process to occur earlier. A curved up and elevated area of the plot indicates the 4-6% change in P max values (Figure 7d) if compared with diesel fuel, and the expected drop in exhaust gas temperature T g (Figure 7e) predetermines the increase in energy efficiency indicators without compromising the cylinder-piston group parts and turbine parts, that might be affected by thermal and mechanical overloads.
The differences in CO 2 emissions (full life-cycle carbon dioxide emissions are absent from analysis) within the zone of rational combination of m and ϕ z parameters did not exceed 4-5% (Figure 7f) if juxtaposed to diesel fuel. Furthermore, an increase in ϕ z up to 70-80 degCA, that indeed led to the deterioration of the indicated efficiency of the engine, resulted in higher CO 2 emission levels by about 5.5-6.0%.
In general, the nature of the combustion process of conventional and alternative fuels is determined by the quality of the fuel spray and its distribution and mixing within the combustion chamber [43,44]. Breakup and distribution of the spray are largely determined by the physical and chemical properties of the fuel, air motion, and temperature that occur within the cylinders, coupled with injection pressure, nozzle design and geometry, and spray angle [42]. Heavy duty, industrial and passenger car diesel engines still require fuel injection equipment to be adopted according ever-renewing stringent requirements, especially for efficient operation at lower load and speed operating points. Since emission legislation drive cycles and EU and US federal exhaust emission standards for stationary engines require operation within these regimes, research into the qualitative variation in the working process parameters within diesel engines of various types and modifications is of great importance. At the same time, the possibility to achieve a desirable characteristic of fuel combustion for the selected compression ratio of an engine is limited by the magnitude of the thermal and mechanical stress acting on the internal parts of ICE. This conclusion is based on the detailed monitoring of the obtained P me and T g values describing the different cases of the heat release process. For the selected value of the indicated efficiency η i = 0.48 − 0.49 of a Cursor 13 engine, two different combinations of ϕ z and m parameters derived from a function X = f(ϕ) (ϕ z = 60-70 degCA, m = 0.5 and ϕ z = 60 degCA, m = 1) may be practically realized to achieve the desirable level of maximum combustion pressure P max = 130-150 bar (at α~2.0). When switching from diesel to MAO100, it is expected that the value of the indicated engine efficiency will drop by 2-3%, however, an existing reserve in P max that comprises 5-7% (7-10 bar for the particular level of P max = 130-150 bar) will open up room for further optimization of injection timing. These examples were embedded in Figure 7b-f, where the RCoP zones were identified for each operational parameter of the engine to facilitate the smooth transition to microalgae oil. Realization of the interim steps described in this study lays down the foundation for the accurate prognostic assessment of the expected operational parameters and CO 2 emissions for the engine family of similar design and a comparable range of P e values for widespread adoption of MAO100 usage.

Conclusions
The authors of this study raise a concern regarding the mass usage of fossil fuels mined from ancient deposits and consumed by industrial engines which are extensively used in various industries and sectors of the economy by offering to replace conventional fuel as a constituting part of the final energy mix with a novel type of biofuel produced from the less investigated microalgae specie P. moriformis. The following conclusions can be drawn as a result of the research:

•
The interval of −2 . . . 0 degCA was found to be the best setting of an engine for smoke and NO x stabilization and reduction, nevertheless D100 or MAO100 were used. That leaves many opportunities for the wider deployment of their binary blends of various ratios to be consumed in diesel engines. Moreover, the pilot study showed that the use of microalgae oil in passenger car engine positively affected the indicated thermal efficiency (η i ) of the prime mover, finding it very similar to that of diesel fuel: 0.355 and 0.350 (P me = 0.8 MPa), 0.350 and 0.345 (P me = 0.6 MPa), 0.325 and 0.320 (P me = 0.4 MPa). • Following accuracy of the 1-D predictive engine model was obtained for various parameters: p me : 0-4.3%, p K : 2.5-4.5%, α: 5.1-10.0%, λ: 0-3.9%, η e : 0.3-3.4% and g cycl : 0-1.7%, T K : 0.8-1.7%, p c : 0.6-1.7%, p max : 1.6-3.9%, T g : 1.9-3.2%.

•
For the CAT 3512B HB-SC engine running with microalgae oil, we proposed a boundary condition for the injection modelling settings (T g ≤ 973 K, ϕ inj = 2 degCA BTDC) that led to improvement of the overall traction characteristics: the difference in η i was almost eliminated and comprised only 0.7-2.0% without any compromise in exceeding the threshold value of 973 K for exhaust gas temperature. • An extensive simulation of the FPT family engine, type Cursor 13 was performed by taking into account different strategies of a combustion process duration and its dynamics through the adjustment of m and ϕ z parameters within the broad range of variation: m = 0-1.5, ϕ z = 50-80 degCA. The obtained results revealed that, if considering the smallest changes in the indicated thermal efficiency values as an outcome of the best compatibility of m and ϕ z parameters, this indicator can be characterized by a relatively short period of heat release (50-60 degCA) and moderate dynamics (m = 0-0.5).

•
The zones of rational combination of m and ϕ z were identified for each operational parameter of the engine to facilitate the smooth transition to microalgae oil. It was found that the differences in carbon dioxide emissions within the zone of rational combination of m and ϕ z parameters did not exceed 4-5% if compared with D100.

•
The study found that microalgae oil is more or less equally sensitive to key engine parameters, compared with diesel fuel, and can be successfully adopted to the entire families of industrial diesel engines.