Search Results

The interaction of turbulent premixed methane combustion with the surrounding flow field can be studied using optically accessible test rigs such as a rapid compression expansion machine (RCEM). The high flexibility offered by such a test rig allows its operation at various thermochemical conditions at ignition. However, limitations inherent to such test rigs due to the absence of an intake stroke do not allow turbulence production as found in IC-engines. Hence, means to introduce turbulence need to be implemented and the relevant turbulence quantities have to be identified in order to enable comparability with engine relevant conditions. A dedicated high-pressure direct injection of air at the beginning of the compression phase is considered as a measure to generate adjustable turbulence intensities at spark timing and during the early flame propagation.

Partially premixed combustion (PPC) can be applied to decrease emissions and increase fuel efficiency in direct injection, compression ignition (DICI) combustion engines. PPC is strongly influenced by the mixing of fuel and oxidizer, which for a given fuel is controlled mainly by (a) the fuel injection, (b) the in-cylinder flow, and (c) the geometry and dynamics of the engine. As the injection timings can vary over a wide range in PPC combustion, detailed knowledge of the in-cylinder flow over the whole intake and compression strokes can improve our understanding of PPC combustion. In computational fluid dynamics (CFD) the in-cylinder flow is sometimes simplified and modeled as a solid-body rotation profile at some time prior to injection to produce a realistic flow field at the moment of injection. In real engines, the in-cylinder flow motion is governed by the intake manifold, the valve motion, and the engine geometry.

The current development to set up an automatic procedure for automatic mesh generation and automatic mesh motion for internal combustion engine simulation in OpenFOAM®-2.2.x is here described. In order to automatically generate high-quality meshes of cylinder geometries, some technical issues need to be addressed: 1) automatic mesh generation should be able to control anisotropy and directionality of the grid; 2) during piston and valve motion, cells and faces must be introduced and removed without varying the overall area and volume of the cells, to avoid conservation errors. In particular, interpolation between discrete fields is frequent in computational physics: the use of adaptive and non-conformal meshes necessitates the interpolation of fields between different mesh regions. Interpolation problems also arise in areas such as model coupling, model initialization and visualisation.

The influence of spray-wall interaction on air entrainment in an unsteady non-evaporating diesel spray was studied using laser Doppler anemometry. The spray was injected into confined quiescent air at ambient pressure and temperature and made to impact on a flat wall. The air velocity component normal to a cylindrical surface surrounding the spray was measured during the entire injection period, allowing to evaluate the time history of the entrained air mass flow rate. The influence of wall distance and spray impingement angle on air entrainment characteristics has been investigated and the results indicate that the presence of a wall increases the entrained mass flow rate in the region close to the surface, during the main injection period. Normal impingement appears to produce stronger effects than oblique incidence at 30 and 45 deg. A qualitative explanation of the results is also proposed, based on the drop-gas momentum exchange mechanism.

Prediction of in-cylinder flows and fuel-air mixing are two fundamental pre-requisites for a successful simulation of direct-injection engines. Over the years, many efforts were carried out in order to improve available turbulence and spray models. However, enhancements in physical modeling can be drastically affected by how the mesh is structured. Grid quality can negatively influence the prediction of organized charge motion structures, turbulence generation and interaction between in-cylinder flows and injected sprays. This is even more relevant for modern direct injection engines, where multiple injections and control of charge motions are employed in a large portion of the operating map. Currently, two different approaches for mesh generation exist: manual and automatic. The first makes generally possible to generate high-quality meshes but, at the same time, it is very time consuming and not completely free from user errors.

The paper focuses on the development of a mesh moving method based on non-conformal topologically changing grids applied to the simulation of IC engines, where the prescribed motion of piston and valves is accomplished by rigidly translating the sub-domain representing the moving component. With respect to authors previous work, a more robust and efficient algorithm to handle the connectivity of non-conformal interfaces and a mesh-motion solver supporting multiple layer addition/removal of cells, to decouple the time-step constraints of the mesh motion and of the fluid dynamics, has been implemented as a C++ library to extend the already existing classes for dynamic mesh handling of the finite-volume, open-source CFD code OpenFOAM®. Other new features include automatic decomposition of large multiple region domains to preserve processors load balance with topological changes for parallel computations and additional tools for automatic preprocessing and case setup.

The dynamics and evolution of turbulent structures inside an engine-like geometry are investigated by means of Large Eddy Simulation. A simplified configuration consisting of a flat-top cylinder head with a fixed, axis-centered valve and low-speed piston has been simulated by the finite volume CFD code OpenFOAM®; the standard version of the software has been extended to include the compressible WALE subgrid-scale model, models for the generation of synthetic turbulence, some improvements to the mesh motion strategy and algorithms for LES data post-processing. In order to study both the initial transient and the quasi- steady operating conditions, ten complete engine cycles have been simulated. Phase and spatial averages have been performed over cycles three to ten in order to extract first and second moment of velocity; these quantities have then been used to validate the numerical procedure by comparison against experimental data.

In this study two different simulation approaches to large eddy simulation of spark-ignition engines are compared. Additionally, some of the simulation results are compared to experimentally obtained in-cylinder velocity measurements. The first approach applies unstructured grids with an automated meshing procedure, using OpenFoam and Lib-ICE with a mapping approach. The second approach applies the efficient in-house code PsiPhi on equidistant, Cartesian grids, representing walls by immersed boundaries, where the moving piston and valves are described as topologically connected groups of Lagrangian particles. In the experiments, two-dimensional two-component particle image velocimetry is applied in the central tumble plane of the cylinder of an optically accessible engine. Good agreement between numerical results and experiment are obtained by both approaches.

The increasing diffusion of gasoline direct injection (GDI) engines requires a more detailed and reliable description of the phenomena occurring during the fuel injection process. As well known the thermal and fluid-dynamic conditions present in the combustion chamber greatly influence the air-fuel mixture process deriving from GDI injectors. GDI fuel sprays typically evolve in wide range of ambient pressure and temperatures depending on the engine load. In some particular injection conditions, when in-cylinder pressure is relatively low, flash evaporation might occur significantly affecting the fuel-air mixing process. In some other particular injection conditions spray impingement on the piston wall might occur, causing high unburned hydrocarbons and soot emissions, so currently representing one of the main drawbacks of GDI engines.

Previously published research [1] covering the role of piston material properties in brake torque variation sensitivity and roughness concluded that phenolic pistons have significantly higher low-pressure range compliance than steel pistons, which promotes lower roughness propensity. It also determined that this property could be successfully characterized using a modern generation of direct-acting servo hydraulically actuated brake component compression test stands. This paper covers a subsequent block of research into the role of the caliper piston in brake torque variation sensitivity (BTV sensitivity) and thermal roughness of a brake corner. It includes measurements of hydraulic stiffness of pistons in a “wet” fixture, both with and without a brake pad and multi-layer bonded noise shim.