1. Introduction
According to [
1], today’s world population is more than three times larger than it was in the mid-twentieth century, and it is still growing. The continuous growth of population, together with the expansion of urbanization, rapid industrialization and economic development, is associated with constantly increasing energy consumption, consumer demands, and an increased need for industrial goods. These factors necessitate an increase in the future demand for services provided by the road freight sector. Indeed, it has been reported in [
2] that while the oil consumption and energy use of the road passenger vehicle fleet have begun to plateau and decline, oil use by the road freight sector has continued to increase. Although the electrification of vehicles is mandated by economic, social and geopolitical trends, no electrification strategy for heavy-duty vehicles is foreseen. One of the favored polices to reduce significant harmful climate and air quality impacts from road freight transport is the utilization of the dual-fuel internal combustion (IC) engine. Indeed, as stated in [
3], dual-fuel IC engines can be utilized in applications where electrification is not deemed to be a viable solution to tackle emission issues. These applications include cargo ships, heavy-duty trucks and marine engines. The concept of dual-fuel combustion has attracted growing interest, primarily due to its fuel flexibility and its ability to achieve lower NO
x and soot emissions compared to conventional diesel engines [
4]. The focus of this paper is the diesel ignited gas engine. This type of engine is operated by burning a premixed natural gas/air mixture ignited with a small amount of directly injected diesel fuel. This concept has several disadvantages, which are discussed in [
5]. Nevertheless, dual-fuel engines are the subject of active research, and establishing an effective design process for dual-fuel engines and its fuel injection equipment represents a key priority.
By leveraging advanced simulation technologies, the cost and time required for the development of dual-fuel IC engine and its fuel injection equipment can be significantly reduced. In particular, 3D computational fluid dynamics (CFD) simulation tools have long been utilized for the design of IC engines due to their potential to offer accurate predictions of the combustion process with a reduced number of prototypes in the development and test phases. Nevertheless, the combustion process involves strong coupling between chemistry, transport and fluid dynamics [
6]. Hence, simulating turbulent combustion is, to date, a challenging task. In addition, the dual-fuel engine combines the characteristics of both the compression and spark ignition engine operational modes. Hence, all combustion regimes (i.e., autoignition, diffusion combustion and premixed flame front propagation) have to be modelled simultaneously [
7]. Furthermore, the introduction of gaseous fuel in the combustion chamber strongly modifies mixture formation and the combustion properties of the pilot diesel spray, which is utilized as an ignition agent [
8]. The work of several authors, both experimental [
9,
10] and numerical [
4,
7] has shown that the ignition delay time of pilot diesel fuel strongly depends on the fuel mixture fraction of natural gas and diesel. These challenges strongly rely on the chemical kinetics of the gaseous fuel–air mixture, which, in turn, varies based on the type and concentration of used gaseous fuel, as well as the quantity of the employed pilot fuel. Thus, numerical modeling of dual-fuel combustion is a complex task. The process can be accurately depicted by utilizing detailed chemical reaction schemes that involve necessary reactions to the estimated consumption rates of both fuels. A detailed understanding of chemistry is essential for the correct prediction of the ignition, stabilization or extinction of reaction zones, as well as for understanding the formation of harmful pollutants. Since detailed chemical mechanisms can involve thousands of chemical species and elementary reactions, a direct resolution of the combustion process’s chemistry requires significant computational power. For the numerical simulation of industrial interest, this computational cost is not affordable. This has led to the development of reduced chemical kinetic mechanisms and combustion models capable of providing faster results. In the work of several authors [
11,
12], reduced chemical kinetic mechanisms for CFD simulation of natural gas/diesel dual-fuel engines have been developed. However, even with the reduced chemistry, simulation of the combustion process in IC engines of industrial interest is limited to the Reynolds-averaged Navier–Stokes (RANS) turbulence modeling framework due to the long computational time [
4]. On the other hand, combustion models can provide valuable results under significantly lower computational effort compared to the detailed chemistry simulation. A successful framework for the simulation of the dual-fuel combustion process based upon the three-zone extended coherent flame model is presented in [
4,
13]. Nevertheless, simplified combustion models sacrifice some level of accuracy to achieve computational efficiency. Indeed, to gain a deeper understanding of the dual-fuel combustion process, including pollutant emissions and flame characteristics, under various operating conditions, detailed chemical kinetics should be employed. With the aim of combining detailed chemical kinetics with computational efficiency, in this study, the tabulated chemistry combustion modeling approach has been employed to predict the combustion process in a large, single-cylinder, dual-fuel IC engine. More specifically, the focus of this paper is the flamelet-generated manifold (FGM) tabulated chemistry approach. The FGM combustion model makes use of the lookup tables generated using the chemical kinetics mechanism with any degree of detail. The pre-tabulation approach ensures affordable computational times without loss of accuracy due to combustion modeling. Thus, the FGM model is ideally suited for 3D CFD simulation of IC engines. The model has proven its capability to predict the combustion process in both diesel [
14,
15] and gasoline [
16] IC engines. To adequately capture turbulence chemistry interactions (TCI), the model relies on a priori assumptions regarding the occurring combustion regime, either premixed or non-premixed. As already stated, dual-fuel combustion involves the combination of combustion regimes. However, simultaneous consideration of the premixed and non-premixed combustion regimes in the same lookup table would significantly increase computational efforts with respect to the memory and table generation time demands. Therefore, the work presented here makes use of two lookup tables (one for the premixed and one for the non-premixed regime) and interpolates the solution from the tables based on the value of the combustion regime indicator calculated during CFD simulation. The methodology for handling multiple fuels within the FGM tabulated chemistry approach is described in detail in
Section 2.2.
In addition, resolving the entire range of spatial and temporal scales that characterize turbulent flames in realistic geometries is not computationally tractable even with state-of-the-art computing resources. Consequently, a common practice for turbulent simulation of industrial interest is to adopt the RANS turbulence modeling approach. However, RANS models are single-point closures, relying on the assumption of self-similarity of the turbulence spectrum [
17]. This assumption suggests that the complete turbulence spectrum can be defined by a single characteristic turbulent length scale. Hence, the physics of the flow dominated by the organized, large-scale coherent structures cannot be captured satisfactorily [
18]. On the other hand, large eddy simulation (LES) can offer a significant advantage over RANS modeling approaches, such as studying cycle-to-cycle variations, which provide more design sensitivity for investigating geometrical and operational changes and produce more detailed results [
19]. Nevertheless, the penalty for the increased range of resolved flow structures is the substantial increase in computational demand. Therefore, the work presented here aims to employ a hybrid RANS/LES modeling approach to combine the benefits of a computationally affordable RANS approach and high-fidelity and detailed LES approach for practical engineering applications. An emerging hybrid turbulence modeling approach that has demonstrated its potential in recent years is the partially averaged Navier–Stokes (PANS) method. PANS is designed to capture essential, large-scale fluctuations while modeling the remaining flow scales. As a result, it offers enhanced outcomes compared to the RANS method, while demanding significantly less computational resources than LES calculations. The successful implementation of the PANS technique has been demonstrated across a wide range of applications (see [
17,
20,
21,
22]). Various PANS variants have been developed to date; for this study, a PANS model based on the
model proposed by Basara et al. [
23] is coupled with the FGM tabulated chemistry combustion model.
The aim of this study is to propose a cost-effective 3D CFD simulation workflow for predicting a dual-fuel combustion process suitable for everyday industrial use. The proposed simulation workflow enables the consideration of detailed chemistry effects during dual-fuel combustion and provides more insight into the turbulent combustion process due to resolving the portion of flow fluctuations. Potential benefits of the proposed FGM methodology for dual-fuel combustion based on two lookup tables have been tested on a wide range of chemical mechanisms with different levels of detail. Additionally, the FGM PANS 3D CFD simulation workflow for numerical simulation of dual-fuel combustion has been validated against available experimental data from a single-cylinder, large diesel ignited gas engine.
4. Results and Discussion
In order to highlight the advantages of the FGM combustion modeling approach in terms of computational efficiency, comparisons were made between FGM simulations and detailed chemistry simulations employing a detailed diesel surrogate chemical kinetics mechanism developed by the Lawrence Livermore National Laboratory (LLNL). This mechanism describes the oxidation of an n-dodecane/m-xylene mixture and incorporates 2885 species and 11,754 elementary reactions [
35]. The numerical results are validated against experimental in-cylinder pressure traces and the rate of heat release curves for three engine operating points.
As a first step, in
Figure 5, a comparison is made between the temperature curves obtained from a 0D calculation using the LLNL chemical mechanism and those from a 3D CFD simulation employing the FGM dual-table combustion modeling approach, utilizing tables generated from the same chemical mechanism. The CFD simulation utilized a straightforward cubic domain comprising 50 hexahedral cells. The results reveal that the ignition delay time from the CFD simulation closely aligns with the ignition delay time determined through the 0D calculation. This congruence suggests that the accuracy of predicting the ignition delay time is maintained even in the presence of chemistry tabulation.
Figure 6 illustrates the results of the 0D perfectly stirred reactor (PSR) simulation using the LLNL chemical kinetics mechanism. The simulations were conducted for a fuel ratio of zero, representing pure diesel, and a fuel ratio of 0.8, corresponding to a dual-fuel mixture with 80% methane content. The conditions for both cases were a temperature of 850 K and a pressure of 55 bar. The progress variable source term is mapped onto a predefined progress variable and mixture fraction grid.
The results indicate that the distribution of the progress variable source is similar for both cases, but the peak values attained for the dual-fuel mixture are lower. Additionally, in both cases, the progress variable source term has its maximum at the stoichiometric mixture fraction. As previously mentioned, the presence of premixed gaseous fuel alters the mixture formation and, consequently, the combustion characteristics of the pilot diesel fuel. Specifically, in this case, methane significantly extends the ignition delay time of diesel fuel. This phenomenon is well-known, having been confirmed both experimentally and numerically. This behavior is further evident in
Figure 7, which depicts the progress variable source plotted against the fuel ratio. As the methane content increases in the dual-fuel mixture, the progress variable source term decreases and consequently prolongs the ignition delay time. To provide an accurate description of the autoignition of the dual-fuel mixture, it is crucial to represent both fuels—methane and diesel—in the reaction chemistry. It is evident that the thermochemical effects of the second fuel are at the very least qualitatively correctly captured within the FGM lookup tables.
In addition,
Figure 8 and
Figure 9 visualize the output of 0D PSR simulations for a table boundary temperature of 1250 K for the pure diesel case and dual-fuel configuration with an 80% methane content. In these figures, the CO
2 and soot mass fraction are mapped on a predefined progress variable and mixture fraction grid, respectively. The results demonstrate an already proven trend [
40]: the introduction of gaseous fuel leads to a reduction in both CO
2 and soot mass fractions. Notably, the peak values occur at a mixture fraction approximately corresponding to the stoichiometric mixture value. Over the years, there has been significant research into the nature and extent of exhaust emissions from dual-fuel engines. Dual-fuel-type compression ignition engines, fueled with various gaseous fuel resources, produce less exhaust emissions than conventional diesel engines without any substantial operating and capital costs [
40]. The tabulated chemistry approach has demonstrated its ability to preserve the inherent characteristics of species formation from the chemical kinetics mechanism.
In
Figure 10,
Figure 11 and
Figure 12, the in-cylinder pressure traces and rate of heat release curves are presented, which were obtained from 3D CFD simulations, utilizing the FGM combustion model and employing distinct approaches to consider turbulence chemistry interactions. The results are validated against experimental results from the single-cylinder, dual-fuel research engine described in
Section 3. The same trend is consistently observed across all three engine operating points. When only the progress variable variance is taken into consideration, the combustion process starts too early, and the peak from the premixed combustion and peak pressure values are overpredicted. On the other hand, when only mixture fraction variance is considered, all of the values are significantly underpredicted, resulting in incomplete fuel combustion. Simulations where the reactor solution is interpolated from both tables align well with experimental data in terms of peak values and the onset of combustion. This concordance indicates that, for an accurate depiction of the dual-fuel combustion process, it is essential to account for both premixed and non-premixed combustion regimes. Furthermore, the presented results indicate that the dual-table FGM combustion modeling approach, which accounts for different combustion regimes through PDF averaging of control variables, holds potential for accurately predicting the dual-fuel combustion process. This approach also demonstrates the advantage of lower computational costs in comparison to the single-table approach that attempts to simultaneously consider both regimes.
As stated previously, the principal aim of this study is to propose an optimal numerical simulation workflow for designing dual-fuel IC engines and their associated injection systems, tailored for everyday industrial applications. Therefore, the dual-table FGM combustion modeling approach has been coupled with the PANS
turbulence model. The work presented in [
22] has already demonstrated the potential of the PANS
turbulence model to accurately predict the dual-fuel combustion process with reduced computational demands, as compared to the approach in which the LES turbulence modeling approach is employed. Moreover, the same study highlighted the PANS
model’s capability to predict cycle-to-cycle variations. Following this, a similar study has been performed in this work. For the late injection engine operating point, the dual-fuel combustion process is simulated using three distinct turbulence modeling approaches coupled with the dual-table FGM combustion model: the RANS
model, PANS
model, and the LES-with-subgrid model based on the coherent structure function as developed by Kobayashi [
41,
42].
Figure 13 shows the comparison of temperature and velocity fields as obtained with these different turbulence modeling methods, all on the same computational mesh. Observing the results, it becomes evident that RANS simulation results are very smooth, with no visible fluctuations or small-scale structures. In contrast, PANS results show more pronounced wrinkle-like patterns due to resolving the portion of the turbulent flow scales. When comparing the PANS results and LES results on the same mesh, it becomes apparent that PANS yields an equivalent level of detail in terms of resolving a portion of fluctuating flow scales. However, a noteworthy distinction arises: the PANS simulation shows more accurate prediction than the LES simulation in the near-wall region. This can be attributed to the fact that the PANS model can be effectively coupled with the universal wall approach, which combines the integration up to the wall with the wall functions. In contrast, for wall-resolved LES calculation, a very fine numerical grid is needed, and most industrial computations cannot meet this stringent requirement [
27]. Based on the results, it can be concluded that the employed numerical mesh resolution is inadequate to accurately capture the involved physical phenomena with the LES CSM model.
Figure 14 illustrates the PANS model’s resolution parameter along with the total, resolved or SSV turbulence, as well as the unresolved turbulent kinetic energy. Notably, the highest value of the resolution parameter is found near the nozzle exit, signifying that this portion of the flow remains beyond the resolution capability of the employed mesh. Furthermore, it is worth mentioning that the resolution parameter values are notably below unity in the central part of the cylinder and equal to unity near the wall. Consequently, within regions where the resolution parameter equals unity, the PANS model transitions to its RANS parent model, leading to a completely modeling turbulent flow. This ensures an adequate description of in-cylinder processes, even in the cases where the mesh resolution is insufficient for resolved simulation in those areas. Correspondingly, the maximum unresolved turbulent kinetic energy is situated in regions where the resolution parameter exhibits higher values. In contrast, the maximum resolved energy is found in areas where the resolution parameter has lower values. The total turbulent kinetic energy is obtained by summing the resolved and unresolved turbulent kinetic energy components.
Finally, the following instantaneous pictures in
Figure 15 illustrate the relationship between the combustion regime indicator and the mass fraction of evaporated pilot fuel, as obtained using a dual-table FGM PANS simulation approach. Evidently, this demonstrates that the implemented combustion regime indicator is capable of distinguishing the premixed, non-premixed, and partially premixed regions within the domain. It is worth noting that this scalar accounts solely for the gaseous phase of the mixture fraction, transitioning to the partially premixed regime only when the spray starts to evaporate.
Up to this point, the dual-table FGM combustion modeling approach has effectively demonstrated its capacity to accurately capture the thermochemical impacts of the chemical kinetics mechanism. Notably, the dual-table approach has proven to be computationally less demanding compared to the single-table approach. Moreover, the PANS turbulence model has exhibited its potential to yield a more detailed solution compared to RANS, on coarser meshes that are typically not sufficient for LES requirements. Hence, this simulation workflow has been employed to simulate another important aspect of the IC engine, namely cycle-to-cycle variations. In order to perform this numerical study on cyclic variations, initial and boundary conditions were perturbed based on the work of [
43,
44,
45]. This test is mainly performed to investigate the feasibility of the dual-table FGM PANS modeling approach to predict cycle-to-cycle variations compared to detailed chemistry simulation in terms of the computational cost and accuracy. In this study, thirty individual cycles were simulated in parallel on 90-degree sector mesh utilizing 60 CPUs. The outcome of this study is shown in
Figure 16. It is evident from the presented results that numerically obtained pressure traces lie in the range of the experimentally measured minimum and maximum in-cylinder pressure. It is evident from
Figure 17 that the numerically obtained pressure trace and rate of heat release averaged over 30 cycles are in good agreement with the average pressure trace and rate of heat release from the measurement. In the same figure, a pressure trace and rate of heat release as obtained using detailed chemistry simulation coupled with a PANS turbulence model are shown. Detailed chemistry simulation is performed utilizing the same chemical kinetics mechanism, namely the previously mentioned LLNL, on the same numerical mesh utilizing 60 CPUs. The solver turnaround time for one high-pressure phase combustion cycle on 90-degree sector mesh for detailed chemistry simulation is 77 days. Performing the multicycle simulation utilizing a detailed chemistry approach was unfeasible; therefore, the required computational time can only be estimated. The computational time required for FGM and detailed chemistry simulation are compared in
Table 8. The FGM tables have to be generated prior to the CFD simulation, and then the same table is used for every cycle. It is evident from the presented results that the FGM combustion modeling approach enables a substantial reduction in the computational time compared to detailed chemistry simulation.