No abstract

Cycle-to-cycle variability is numerically simulated for high-speed, full-load operation of a turbocharged gasoline direct injection engine. Large-Eddy simulation is adopted to replicate the fluctuations of the flow field affecting the turbulent combustion. Experimental data were provided at knock onset, and large-Eddy simulation was validated for the same condition. In the original engine configuration, the spark plug is displaced toward the exhaust side, while the electrodes orientation is arbitrary. A 90 rotation is imposed to evaluate the effects of the aerodynamic obstruction caused by the electrode with respect to the flow field and the flame kernel growth. A second speculative analysis is performed modifying the position of the spark plug. The electrodes are shifted 2mm toward the intake side since this variation is compatible with the cylinder head layout. For both variations in orientation and position, the effects on the flow field around the spark plug are investigated. Statistical analysis is carried out on early flame kernel formation and knock tendency. The results highlight that the orientation of the electrodes affects the flow field for each cycle but plays a negligible role on the statistical cyclic variability, indirectly justifying the lack of an imposed orientation. As for the spark plug position, the numerical analysis indicate that the shifting of the electrodes toward the intake side slightly improves the knock limit mainly because of a reduction in in-cylinder peak pressure. In general, it is inferred that improvements may be achieved only through a simultaneous modification of the fuel jet orientation and phasing

No abstract

No abstract

A zonal hybridization of the RNG [Formula: see text]-[Formula: see text] URANS model is proposed for the simulation of turbulent flows in internal combustion engines. The hybrid formulation is able to act as URANS, DES or LES in different zones of the computational domain, which are explicitly set by the user. The resulting model has been implemented in a commercial computational fluid dynamics code and the LES branch of the modified RNG [Formula: see text]-[Formula: see text] closure has been initially calibrated on a standard homogeneous turbulence box case. Subsequently, the full zonal formulation has been tested on a fixed intake valve geometry, including comparisons with third-party experimental data. The core of the work is represented by a multi-cycle analysis of the TCC-III experimental engine configuration, which has been compared with the experiments and with prior full-LES computational studies. The applicability of the hybrid turbulence model to internal combustion engine flows is demonstrated, and PIV-like flow statistics quantitatively validate the model performance. This study shows a pioneering application of zonal hybrid models in engine-relevant simulation campaigns, emphasizing the relevance of hybrid models for turbulent engine flows.

As known, reliable information about underlying turbulence intensity is a mandatory prerequisite to predict the burning rate in quasidimensional combustion models. Based on 3D results reported in the companion part I paper, a quasi-dimensional turbulence model, embedded under the form of "user routine" in the GT-Power™ software, is here presented in detail. A deep discussion on the model concept is reported, compared to the alternative approaches available in the current literature. The model has the potential to estimate the impact of some geometrical parameters, such as the intake runner orientation, the compression ratio, or the bore-to-stroke ratio, thus opening the possibility to relate the burning rate to the engine architecture. Preliminarily, a well-assessed approach, embedded in GT-Power commercial software v.2016, is utilized to reproduce turbulence characteristics of a VVA engine. This test showed that the model fails to predict tumble intensity for particular valve strategies, such LIVC, thus justifying the need for additional refinements. The model proposed in this work is conceived to solve 3 balance equations, for mean flow kinetic energy, tumble vortex momentum, and turbulent kinetic energy (3-eq. concept). An extended formulation is also proposed, which includes a fourth equation for the dissipation rate, allowing to forecast the integral length scale (4-eq. concept). The impact of the model constants is parametrically analyzed in a first step, and a tuning procedure is advised. Then, a comparison between the 3-and the 4-eq. concepts is performed, highlighting the advantages of the 3-eq. version, in terms of prediction accuracy of turbulence speed-up at the end of the compression stroke. An extensive 3-eq. model validation is then realized according to different valve strategies and engine speeds. The user-model is then utilized to foresee the effects of main geometrical parameters analyzed in part I, namely the intake runner orientation, the compression ratio, and the bore-to-stroke ratio. A two-valve per cylinder engine is also considered. Temporal evolutions of 0D-and 3D-derived mean flow velocity, turbulent intensity, and tumble velocity present very good agreements for each investigated engine geometry and operating condition. The model, particularly, exhibits the capability to accurately predict the tumble trends by varying some geometrical parameter of the engine, which is helpful to estimate the related impact on the burning rate. Summarizing, the developed 0D model well estimates the in-cylinder turbulence characteristics, without requiring any tuning constants adjustment with engine speed and valve strategy. In addition, it demonstrates the capability to properly take into account the intake duct orientation and the compression ratio without tuning adjustments. Some minor tuning variation allows predicting the effects of bore-to-stroke ratio, as well. Finally, the model is verified to furnish good agreements also for a two-valve per cylinder engine, and with reference to two different h...

The paper reports an activity aiming at characterizing cycle-to-cycle variability (CCV) of the spark-ignition (SI) process in a high performance engine. The numerical simulation of spark-ignition and of early flame kernel evolution are major challenges, mainly due to the time scales of the spark discharge process and to the reduced spatial scales of flame kernel. Typical mesh resolutions are insufficient to resolve the process and a dedicated treatment has to be provided at a subgrid level if the ignition process is to be properly modelled. The focus of this work is on the recent ISSIM-LES (Imposed Stretch Spark-Ignition Model) ignition model, which is based on an extension of the flame surface density (FSD) transport equation for a dedicated flame kernel treatment at subgrid scales. The FSD equation is solved immediately after spark discharge. The interaction of the flame kernel with the flow field is fully accounted for since spark formation and a transition is provided from ignition to propagation phase. The comparison is carried out with the AKTIM-Euler ignition model in terms of flame interaction with the flow field (e.g. arc convection, flame blow-off, flame holder effect). A multiple cycle LES activity provided a set of cycle-resolved conditions for spark-ignition comparisons, and the flame kernel development is carefully analyzed for the two ignition models on a wide range of thermo-physical conditions. Spark-ignition cyclic\ud variability and combustion traces are compared with experiments. Results confirm that the simulated cycle-to-cycle variability increases through the adoption of the ISSIM-LES ignition model

No abstract