Search Results

Although the application of cylinder pressure sensors to gain insight into the combustion process is not a novel topic itself, the recent availability of inexpensive in-cylinder pressure sensors has again prompted an upcoming interest for the utilization of the cylinder pressure signal within engine control and monitoring. Besides the use of the in-cylinder pressure signal for combustion analysis and control the information can also be used to determine related quantities in the exhaust or intake manifold. Within this work two different methods to estimate the pressure inside the exhaust manifold are proposed and compared. In contrary to first principle based approaches, which may require time extensive parameterization, alternative data driven approaches were pursued. In the first method a Principle Component Analysis (PCA) is applied to extract the cylinder pressure information and combined with a polynomial model approach.

In recent years, the increase of gasoline fuel injection pressure is a way to improve thermal efficiency and lower engine-out emissions in GDI homogenous combustion concept. The challenge of controlling particulate formation as well in mass and number concentrations imposed by emissions regulations can be pursued improving the mixture preparation process and avoiding mixture inhomogeneity with ultra-high injection pressure values up to 100 MPa. The increase of the fuel injection pressure in GDI homogeneous systems meets the demand for increased injector static flow, while simultaneously improves the spray atomization and mixing characteristics with consequent better combustion performance. Few studies quantify the effects of high injection pressure on transient gasoline spray evolution. The aim of this work was to simulate with OpenFOAM the spray morphology of a commercial gasoline injected in a constant volume vessel by a prototypal GDI injector.

The 4th Workshop of the Engine Combustion Network (ECN) was held September 5-6, 2015 in Kyoto, Japan. This manuscript presents a summary of the progress in experiments and modeling among ECN contributors leading to a better understanding of soot formation under the ECN “Spray A” configuration and some parametric variants. Relevant published and unpublished work from prior ECN workshops is reviewed. Experiments measuring soot particle size and morphology, soot volume fraction (fv), and transient soot mass have been conducted at various international institutions providing target data for improvements to computational models. Multiple modeling contributions using both the Reynolds Averaged Navier-Stokes (RANS) Equations approach and the Large-Eddy Simulation (LES) approach have been submitted. Among these, various chemical mechanisms, soot models, and turbulence-chemistry interaction (TCI) methodologies have been considered.

The modeling of fuel sprays under well-characterized conditions relevant for heavy-duty Diesel engine applications, allows for detailed analyses of individual phenomena aimed at improving emission formation and fuel consumption. However, the complexity of a reacting fuel spray under heavy-duty conditions currently prohibits direct simulation. Using a systematic approach, we extrapolate available spray models to the desired conditions without inclusion of chemical reactions. For validation, experimental techniques are utilized to characterize inert sprays of n-dodecane in a high-pressure, high-temperature (900 K) constant volume vessel with full optical access. The liquid fuel spray is studied using high-speed diffused back-illumination for conditions with different densities (22.8 and 40 kg/m3) and injection pressures (150, 80 and 160 MPa), using a 0.205-mm orifice diameter nozzle.

The implementation and the combination of advanced boundary conditions and subgrid scale models for Large Eddy Simulation (LES) in the multi-dimensional open-source CFD code OpenFOAM® are presented. The goal is to perform reliable cold flow LES simulations in complex geometries, such as in the cylinders of internal combustion engines. The implementation of a boundary condition for synthetic turbulence generation upstream of the valve port and of the compressible formulation of the Wall-Adapting Local Eddy-viscosity sgs model (WALE) is described. The WALE model is based on the square of the velocity gradient tensor and it accounts for the effects of both the strain and the rotation rate of the smallest resolved turbulent fluctuations and it recovers the proper y₃ near-wall scaling for the eddy viscosity without requiring dynamic procedure; hence, it is supposed to be a very reliable model for ICE simulation.

Modeling plume interaction and collapse for direct-injection gasoline sprays is important because of its impact on fuel-air mixing and engine performance. Nevertheless, the aerodynamic interaction between plumes and the complicated two-phase coupling of the evaporating spray has shown to be notoriously difficult to predict. With the availability of high-speed (100 kHz) Particle Image Velocimetry (PIV) experimental data, we compare velocity field predictions between plumes to observe the full temporal evolution leading up to plume merging and complete spray collapse. The target “Spray G” operating conditions of the Engine Combustion Network (ECN) is the focus of the work, including parametric variations in ambient gas temperature. We apply both LES and RANS spray models in different CFD platforms, outlining features of the spray that are most critical to model in order to predict the correct aerodynamics and fuel-air mixing.

The flow field resulting from injecting a gas jet into a crossflow confined in a narrow square duct has been studied under steady regime using schlieren imaging and laser Doppler velocimetry (LDV). This transparent duct is intended to simulate the intake port of an internal combustion engine fueled by gaseous mixture, and the jet is issued from a round nozzle. The schlieren images show that the relative small size of the duct would confine the development of the transverse jet, and the interaction among jet and sidewalls strongly influences the mixing process between jet and crossflow. The mean velocity and turbulence fields have been studied in detail through LDV measurements, at both center plane and several cross sections. The well-known flow feature formed by a counter rotating vortex pair (CVP) has been observed, which starts to appear at the jet exit section and persists far downstream contributing to enhancing mixing process.

This paper is focused on the development and application of a CFD methodology that can be applied to predict the fuel-air mixing process in stratified charge, sparkignition engines. The Eulerian-Lagrangian approach was used to model the spray evolution together with a liquid film model that properly takes into account its effects on the fuel-air mixing process into account. However, numerical simulation of stratified combustion in SI engines is a very challenging task for CFD modeling, due to the complex interaction of different physical phenomena involving turbulent, reacting and multiphase flows evolving inside a moving geometry. Hence, for a proper assessment of the different sub-models involved a detailed set of experimental optical data is required. To this end, a large experimental database was built by the authors.

Dynamic tests have been performed on carbon fiber racing seats following the FIA regulations. The tests have shown, in rear impact tests, a relatively strong rebound leading to large forward bending of neck, and, in side impact tests, very large lateral displacement of the head, the latter protruding dangerously towards hard portions of the car structure. Stiffening the seat back by steel struts results in reducing strongly both the motion and the acceleration of the head. Simulations of the dynamics of the tests have been done with multi-body models, including the Hybrid III dummy and seat deflection, by means of the program VEDYAC. It has been found that computer simulation can predict very accurately the result of a test, provided the numerical models have been carefully calibrated to match the dummy tolerance bands. Once they have been calibrated and validated with a number of tests, the computer models can be very useful to extend the test results to different test conditions.

The scope of the work presented in this paper was to apply the latest open source CFD achievements to design a state of the art, direct-injection (DI), heavy-duty, natural gas-fueled engine. Within this context, an initial steady-state analysis of the in-cylinder flow was performed by simulating three different intake ducts geometries, each one with seven different valve lift values, chosen according to an estabilished methodology proposed by AVL. The discharge coefficient (Cd) and the Tumble Ratio (TR) were calculated in each case, and an optimal intake ports geometry configuration was assessed in terms of a compromise between the desired intensity of tumble in the chamber and the satisfaction of an adequate value of Cd. Subsequently, full-cycle, cold-flow simulations were performed for three different engine operating points, in order to evaluate the in-cylinder development of TR and turbulent kinetic energy (TKE) under transient conditions.

Extensive prior art within the Engine Combustion Network (ECN) using a Bosch single axial-hole injector called ‘Spray A’ in constant-volume vessels has provided a solid foundation from which to evaluate modeling tools relevant to spray combustion. In this paper, a new experiment using a Bosch three-hole nozzle called ‘Spray B’ mounted in a 2.34 L heavy-duty optical engine is compared to sector-mesh engine simulations. Two different approaches are employed to model combustion: the ‘well-mixed model’ considers every cell as a homogeneous reactor and employs multi-zone chemistry to reduce the computational time. The ‘flamelet’ approach represents combustion by an ensemble of laminar diffusion flames evolving in the mixture fraction space and can resolve the influence of mixing, or ‘turbulence-chemistry interactions,’ through the influence of the scalar dissipation rate on combustion.

Partially premixed compression ignition combustion is one of the low temperature combustion techniques which is being actively investigated. This approach provides a significant reduction of both soot and NOx emissions. Comparing to the homogeneous charge compression ignition mode, PPCI combustion provides better control on ignition timing and noise reduction through air-fuel mixture stratification which lowers heat release rate compared to other advanced combustion modes. In this work, CFD simulations were conducted for a low and a high air-fuel mixture stratification cases on a light-duty optical engine operating in PPCI mode. Such conditions for PRF70 as fuel were experimentally achieved by injection timing and spray targeting at similar thermodynamic conditions.

A three dimensional numerical simulation of the ECN “Spray A” is presented. Both primary and secondary breakup of the spray are included. The fuel is n-Dodecane. The n-Dodecane kinetic mechanism is modeled using a skeletal mechanism that consists of 103 species and 370 reactions [9]. The kinetic mechanism is computationally heavy when coupled with three dimensional numerical simulations. Multidimensional chemistry coordinate mapping (CCM) approach is used to speedup the simulation. CCM involves two-way mapping between CFD cells and a discretized multidimensional thermodynamic space, the so called multidimensional chemistry coordinate space. In the text, the cells in the discretized multidimensional thermodynamic space are called zone to discriminate them from the CFD cells. In this way, the CFD cells which are at the similar thermodynamic state are identified and grouped into a unique zone. The stiff ODEs operates only on the zones containing at least one CFD cell.

Highway and roadway surface measurement is a practice that has been ongoing for decades now. This sort of measurement is intended to ensure a safe level of road perturbances. The measurement may be conducted by a slow moving apparatus directly measuring the elevation of the road, at varying distance intervals, to obtain a road profile, with varying degrees of resolution. An alternate means is to measure the surface roughness at highway speeds using accelerometers coupled with high speed distance measurements, such as laser sensors. Vehicles out rigged with such a system are termed inertial profilers. This type of inertial measurement provides a sort of filtered roadway profile. Much research has been conducted on the analysis of highway roughness, and the associated metrics involved. In many instances, it is desirable to maintain an off-road course such that the course will provide sufficient challenges to a vehicle during durability testing.

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.

The air entrainment in a transient diesel spray was studied using laser Doppler anemometry to provide information on the effect of gas density and temperature. The spray was injected vertically into a confined quiescent atmosphere and the entrained mass flow rate was evaluated by measuring the air velocity component normal to a cylindrical geometric surface surrounding the spray, and extending to about 200 nozzle diameters (50 mm). The experimental results, relative to a density range from 0.84 to 7.02 kg/m3 and a temperature range from 293 to 473 K, indicate that the non dimensional entrainment rate, averaged in time over the main injection period, depends on the distance from the nozzle and both gas density and temperature. A first analysis, based on the available data, allowed to quantify the dependence and provided a correlation with such variables.

This work presents a comprehensive numerical model for ice accretion and Ice Protection System (IPS) simulation over a 2D component, such as an airfoil. The model is based on the Myers model for ice accretion and extended to include the possibility of a heated substratum. Six different icing conditions that can occur during in-flight ice accretion with an Electro-Thermal Ice Protection System (ETIPS) activated are identified. Each condition presents one or more layers with a different water phase. Depending on the heat fluxes, there could be only liquid water, ice, or a combination of both on the substratum. The possible layers are the ice layer on the substratum, the running liquid film over ice or substratum, and the static liquid film between ice and substratum caused by ice melting. The last layer, which is always present, is the substratum. The physical model that describes the evolution of these layers is based on the Stefan problem. For each layer, one heat equation is solved.

This paper presents a novel fully-automatic remeshing procedure, based on the level-set method and Delaunay triangulation, to model three-dimensional boundary problems and generate a new conformal body-fitted mesh. The proposed methodology is applied to long-term in-flight ice accretion, which is characterized by the formation of extremely irregular ice shapes. Since ice accretion is coupled with the aerodynamic flow field, a multi-step procedure is implemented. The total icing exposure time is subdivided into smaller time steps, and at each time step a three-dimensional body-fitted mesh, suitable for the computation of the aerodynamic flow field around the updated geometry, is generated automatically. The methodology proposed can effectively deal with front intersections, as shown with a manufactured example.

Active prechamber combustion systems for SI engines represent a feasible and effective solution in reducing fuel consumption and pollutant emissions for both marine and ground heavy-duty engines. However, reliable and low-cost numerical approaches need to be developed to support and speed-up their industrial design considering their geometry complexity and the involved multiple flow length scales. This work presents a CFD methodology based on the RANS approach for the simulation of active prechamber spark-ignition engines. To reduce the computational time, the gas exchange process is computed only in the prechamber region to correctly describe the flow and mixture distributions, while the whole cylinder geometry is considered only for the power-cycle (compression, combustion and expansion). Outside the prechamber the in-cylinder flow field at IVC is estimated from the measured swirl ratio.