The occurrence of knocking combustion is limiting the efficiency of modern spark ignition engine operation. Thus, an understanding of the processes at the knock limit is required for further optimization of the combustion process. In this work, the combustion of a multicomponent Toluene Reference Fuel (TRF) in a single-cylinder research engine is investigated under knocking conditions. The fuel exhibits a negative temperature coefficient (NTC) regime for thermodynamic conditions relevant to the engine operation. A precursor model is used to capture the auto-ignition process. Under homogeneous conditions, a two-stage auto-ignition is observed. Inside the NTC regime, the temperature affects both first-stage and second-stage auto-ignition delay times. With a subsequently conducted multi-cycle engine LES, the effects of temperature stratification and turbulent flame propagation on the local auto-ignition process are investigated. It is observed, that the NTC behavior leads to a widespread two-stage auto-ignition. The knock intensity observed in the experiments is directly related to the mass consumed by auto-ignition. This is due to the fast consumption of the auto-ignited mass by the flame front. With that, the NTC behavior affects the local auto-ignition process in the unburned mixture while the flame propagation determines the knock intensity for the operating conditions at the knock limit.

In this study, a single-cylinder direct-injection spark-ignition research engine with full optical access was used to investigate the influence of the flow field and fuel/air mixing on cyclic variability, in particular in the early flame propagation. The engine was operated under lean-burn conditions at 1500 rpm. Two different injection strategies were compared, port-fuel injection (PFI) and direct injection (DI), the latter with early and late injection split about 2:1 in fuel mass. High-speed particle image velocimetry captured the flow in the tumble plane in the compression stroke. The velocity fields and the movement of the tumble vortex are analyzed. Simultaneously, a second camera detected the chemiluminescence of the flame, and the projected area of the line-of-sight-integrated flame luminosity was extracted through morphological image processing. By combining pressure-based combustion analysis and high-speed optical diagnostics, the early flame propagation and the flow field are correlated. In separate experiments the equivalence ratio was imaged for the DI at selected crank angles and correlated with CA10 to learn about the influence of mixture inhomogeneity on early flame propagation. With PFI, the flow near the spark plug just before ignition is closely related to the subsequent speed of combustion. The combustion-relevant flow features can be traced back in time to about –90 °CA. In contrast, the chosen DI scheme results in a decorrelation of flow and flame, and the equivalence ratio distribution at ignition becomes more important. For both flow and mixture fields, regions of high correlation with early-combustion metrics are typically associated with gradients in the multi-cycle average fields.

The cross-scalar dissipation rate of reaction progress variable and mixture fraction widetilde{{varepsilon_{cxi } }} plays an important role in the modelling of stratified combustion. The evolution and statistical behaviour of widetilde{{varepsilon_{cxi } }} have been analysed using a direct numerical simulation (DNS) database of statistically planar turbulent stratified flames with a globally stochiometric mixture. A parametric analysis has been conducted by considering a number of DNS cases with a varying initial root-mean-square velocity fluctuation u^{prime } and initial scalar integral length scale ell_{phi }. The explicitly Reynolds averaged DNS data suggests that the linear relaxation model for widetilde{{varepsilon_{cxi } }} is inadequate for most cases, but its performance appears to improve with increasing initial ell_{phi } and u^{prime } values. An exact transport equation for widetilde{{varepsilon_{cxi } }} has been derived from the first principle, and the budget of the unclosed terms of the widetilde{{varepsilon_{cxi } }} transport equation has been analysed in detail. It has been found that the terms arising from the density variation, scalar-turbulence interaction, chemical reaction rate and molecular dissipation rate play leading order roles in the widetilde{{varepsilon_{cxi } }} transport. These observations have been justified by a scaling analysis, which has been utilised to identify the dominant components of the leading order terms to aid model development for the unclosed terms of the widetilde{{varepsilon_{cxi } }} transport equation. The performances of newly proposed models for the unclosed terms have been assessed with respect to the corresponding terms extracted from DNS data, and the newly proposed closures yield satisfactory predictions of the unclosed terms in the widetilde{{varepsilon_{cxi } }} transport equation.

Regulations to reduce greenhouse gas emissions of passenger vehicles are becoming increasingly stringent. The aerodynamic drag is a major contributor to the vehicle’s total energy consumption where a large portion is attributed to the base wake. This paper optimises the angles of small trailing edge flaps on a base cavity of a full-scale sports utility vehicle placed in a wind tunnel. The trailing edge flaps are controlled using servos mounted inside the cavity. The flap angles are optimised using a surrogate model based optimisation algorithm with the objective of reducing the aerodynamic drag at different yaw angles and to create a yaw-insensitive geometry by considering several weighted yaw angles to form the driving cycle averaged drag. Low drag designs are further investigated using base pressures and wake measurements. The results show that the base pressures are symmetrised by reducing the crossflow in the wake. As the model is yawed the wake becomes increasingly downwash dominated by a large rotating windward structure which is reduced by the optimised flaps. The cycle averaged drag optimised design has a smaller increase in drag when yawed compared to a design optimised without considering yaw.

The effects of varying turbulence intensity and turbulence length scale on premixed turbulent flame propagation are investigated using Direct Numerical Simulation (DNS). The DNS dataset contains the results of a set of turbulent flame simulations based on separate and systematic changes in either turbulence intensity or turbulence integral length scale while keeping all other parameters constant. All flames considered are in the thin reaction zones regime. Several aspects of flame behaviour are analysed and compared, either by varying the turbulence intensity at constant integral length scale, or by varying the integral length scale at constant turbulence intensity. The turbulent flame speed is found to increase with increasing turbulence intensity and also with increasing integral length scale. Changes in the turbulent flame speed are generally accounted for by changes in the flame surface area, but some deviation is observed at high values of turbulence intensity. The probability density functions (pdfs) of tangential strain rate and mean flame curvature are found to broaden with increasing turbulence intensity and also with decreasing integral length scale. The response of the correlation between tangential strain rate and mean flame curvature is also investigated. The statistics of displacement speed and its components are analysed, and the findings indicate that changes in response to decreasing integral length scale are broadly similar to those observed for increasing turbulence intensity, although there are some interesting differences. These findings serve to improve current understanding of the role of turbulence length scales in flame propagation.

The generation of initial or inflow synthetic turbulent velocity or scalar fields reproducing statistical characteristics of realistic turbulence is still a challenge. The synthetic eddy method, previously introduced in the context of inflow conditions for large eddy simulations, is based on the assumption that turbulence can be regarded as a superposition of coherent structures. In this paper, a new type of synthetic eddy method is proposed, where the fundamental eddy is constructed by superposing three Hill's vortices, with their axes orthogonal to each other. A distribution of Hill's vortices is used to synthesize an anisotropic turbulent velocity field that satisfies the incompressibility condition and match a given Reynolds stress tensor. The amplitudes of the three vortices that form the fundamental eddy are calculated from known Reynolds stress profiles through a transformation from the physical reference frame to the principal-axis reference frame. In this way, divergence-free anisotropic turbulent velocity fields are obtained that can reproduce a given Reynolds stress tensor. The model was tested on both isotropic and anisotropic turbulent velocity fields, in the framework of grid turbulence decay and turbulent channel flow, respectively. The transition from artificial to realistic turbulence in the proximity to the inflow boundary was found to be small in all test cases that were considered.