Search Results

A numerical investigation on a series of Diesel spray experiments in constant-volume vessels is proposed. Non reacting conditions were used to assess the spray models and to determine the grid size required to correctly predict the fuel-air mixture formation process. To this end, not only computed liquid and vapor penetrations were compared with experimental data, but also a detailed comparison between computed and experimental mixture fraction distributions was performed at different distances from the injector. Grid dependency was reduced by introducing an Adaptive Local Mesh Refinement technique (ALMR) with an arbitrary level of refinement. Once the capabilities of the current implemented spray models have been assessed, reacting conditions at different ambient densities and temperatures were considered. A Perfectly Stirred Reactor (PSR) combustion model, based on a direct integration of complex chemistry mechanisms over a homogenous cell, was adopted.

Detailed chemistry and turbulence-chemistry interaction need to be properly taken into account for a realistic combustion simulation of IC engines where advanced combustion modes, multiple injections and stratified combustion involve a wide range of combustion regimes and require a proper description of several phenomena such as auto-ignition, flame stabilization, diffusive combustion and lean premixed flame propagation. To this end, different approaches are applied and the most used ones rely on the well-stirred reactor or flamelet assumption. However, well-mixed models do not describe correctly flame structure, while unsteady flamelet models cannot easily predict premixed flame propagation and triple flames. A possible alternative for them is represented by transported probability density functions (PDF) methods, which have been applied widely and effectively for modeling turbulent reacting flows under a wide range of combustion regimes.

A new multi-zone model for the simulation of HCCI engine is here presented. The model includes laminar and turbulent diffusion and conduction exchange between the zones and the last improvements on the numerical aspects. Furthermore, a new strategy for the zone discretization is presented, which allows a better description of the near-wall zones. The aim of the work is to provide a fast and reliable model for carrying out chemical analysis with detailed kinetic schemes. A preliminary sensitivity analysis allows to verify that 10 zones are a convenient number for a good compromise between the computational effort and the description accuracy. The multi-zone predictions are then compared with the CFD ones to find the effective turbulence parameters, with the aim to describe the near-wall phenomena, both in a reactive and non-reactive cases.

Swirling flows are very dominant in applied technical problems, especially in IC engines, and their prediction requires rather sophisticated modeling. An adaptive low-pass filtering procedure for the modeled turbulent length and time scales is derived and applied to Menter' original k - ω SST turbulence model. The modeled length and time scales are compared to what can potentially be resolved by the computational grid and time step. If the modeled scales are larger than the resolvable scales, the resolvable scales will replace the modeled scales in the formulation of the eddy viscosity; therefore, the filtering technique helps the turbulence model to adapt in accordance with the mesh resolution and the scales to capture.

Modeling combustion of transportation fuels remains a difficult task due to the extremely large number of species constituting commercial gasoline and diesel. However, for this purpose, multi-component surrogate fuel models with a reduced number of key species and dedicated reaction subsets can be used to reproduce the physical and chemical traits of diesel and gasoline, also allowing to perform CFD calculations. Recently, a detailed surrogate fuel kinetic model, named C3 mechanism, was developed by merging high-fidelity sub-mechanisms from different research groups, i.e. C0-C4 chemistry (NUI Galway), linear C6-C7 and iso-octane chemistry (Lawrence Livermore National Laboratory), and monocyclic aromatic hydrocarbons (MAHs) and polycyclic aromatic hydrocarbons (PAHs) (ITV-RWTH Aachen and CRECK modelling Lab-Politecnico di Milano).

This paper presents calculations of engine flows by using the Partially-Averaged Navier Stokes (PANS) method (Girimaji [1]; [2]). The PANS is a scale-resolving turbulence computational approach designed to resolve large scale fluctuations and model the remainder with appropriate closures. Depending upon the prescribed cut-off length (filter width) the method adjusts seamlessly from the Reynolds-Averaged Navier-Stokes (RANS) to the Direct Numerical Solution (DNS) of the Navier-Stokes equations. The PANS method was successfully used for many applications but mainly on static geometries, e.g. Basara et al. [3]; [4]. This is due to the calculation of the cut-off control parameter which requires that the resolved kinetic energy is known and this is usually obtained by suitably averaging of the resolved field. Such averaging process is expensive and impractical for engines as it would require averaging per cycles.

Reynolds-averaged Navier-Stokes (RANS) computations of heat transfer involving wall bounded flows at elevated Prandtl numbers typically suffer from a lack of accuracy and/or increased mesh dependency. This can be often attributed to an improper near-wall turbulence modeling and the deficiency of the wall heat transfer models (based on the so called P-functions) that do not properly account for the variation of the turbulent Prandtl number in the wall proximity (y+< 5). As the conductive sub-layer gets significantly thinner than the viscous velocity sub-layer (for Pr >1), treatment of the thermal buffer layer gains importance as well. Various hybrid strategies utilize blending functions dependent on the molecular Prandtl number, which do not necessarily provide a smooth transition from the viscous/conductive sub-layer to the logarithmic region.

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 present work introduces an extended particulate filter model focusing on capabilities to cover catalytic and surface storage reactions and to serve as a virtual multi-functional reactor/separator. The model can be classified as a transient, non-isothermal 1D+1D two-channel model. The applied modeling framework offers the required modeling depth to investigate arbitrary catalytic reaction schemes and it follows the computational requirement of running in real-time. The trade-off between model complexity and computational speed is scalable. The model is validated with the help of an analytically solved reference and the model parametrization is demonstrated by simulating experimentally given temperatures of a heat-up measurement. The detailed 1D+1D model is demonstrated in a concept study comparing the impact of different spatial washcoat distributions.

The present work describes an existing transient, non-isothermal 1D+1D particulate filter model to capture the impact of different types of particulate matter (PM) on filtration and regeneration. PM classes of arbitrary characteristics (size, composition etc.) are transported and filtered following standard mechanisms. PM deposit populations of arbitrary composition and contact states are used to describe regeneration on a micro-kinetical level. The transport class and deposit population are linked by introducing a splitting deposit matrix. Filtration and regeneration modes are compared to experimental data from literature and a brief numerical assessment on the filtration model is performed. The filter model as part of an exhaust line is used in a concept study on different coating variants. The same exhaust line model is connected to an engine thermodynamic and vehicle model. This system model is run through a random drive cycle in office simulation.

A computational, three-dimensional approach to investigate the behavior of diesel soot particles in the micro-channels of wall-flow Diesel Particulate Filters is presented. The KIVA3V CFD code, already extended to solve the 2D conservation equations for porous media materials [1], has been enhanced to solve in 2-D and 3-D the governing equations for reacting and compressible flows through porous media in non axes-symmetric geometries. With respect to previous work [1], a different mathematical approach has been followed in the implementation of the numerical solver for porous media, in order to achieve a faster convergency as source terms were added to the governing equations. The Darcy pressure drop has been included in the Navier-Stokes equations and the energy equation has been extended to account for the thermal exchange between the gas flow and the porous wall.

In the last decades the continuously tightening limitations on pollutant emissions has led to an extensive adoption of after-treatment devices on the exhaust systems of modern internal combustion engines. While these devices are primarily introduced for reducing and controlling the emissions, they also play an important role influencing the wave motion inside the exhaust system and so affecting the acoustics and the performances of the engine. In this paper a novel approach is proposed for the modeling of two after-treatment devices: the catalyst and the Diesel Particulate Filter. The models are based on a fast quasi-3D approach, named 3Dcell, originally developed by the authors for the acoustic modeling of silencers. This approach allows to model the wave motion by solving the momentum equation along the three directions.

The presented paper deals with a methodology to model cycle-to-cycle variations (CCV) in 0-D/1-D simulation tools. This is achieved by introducing perturbations of combustion model parameters. To enable that, crank angle resolved data of individual cycles (pressure traces) have to be available for a reasonable number of engine cycles. Either experimental data or 3-D CFD results can be applied. In the presented work, experimental data of a single-cylinder research engine were considered while predicted LES 3-D CFD results will be tested in the future. Different engine operating points were selected - both stable ones (low CCV) and unstable ones (high CCV). The proposed methodology consists of two major steps. First, individual cycle data have to be matched with the 0-D/1-D model, i.e., combustion model parameters are varied to achieve the best possible match of pressure traces - an automated optimization approach is applied to achieve that.

The Particle Measurement Programme (PMP) developed an exhaust particle number measurement protocol that has been adopted by current light duty vehicle emission regulations in Europe. This includes thermal treatment of the exhaust aerosol to isolate solid particles only and a number counting device with a lower cutpoint of 23 nm to avoid measurement of smaller particles that may affect the repeatability of the measurement. In this paper, we examine a potential alternative to the PMP system, where the thermal treatment is replaced by a catalytic stripper (CS). This offers oxidation and not just evaporation of the volatile components. Alternative sampling systems, either fulfilling the PMP recommendations or utilizing a CS, have been explored in terms of their volatile particle removal efficiency. Tests have been conducted on diesel exhaust, diesel equipped with DPF and gasoline direct injection emissions.

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 describes the findings of a design, simulation and test study into how to reduce particulate number (Pn) emissions in order to meet EU6c legislative limits. The objective of the study was to evaluate the Pn potential of a modern 6-cylinder engine with respect to hardware and calibration when fitted to a full size SUV. Having understood this capability, to redesign the combustion system and optimise the calibration in order to meet an engineering target value of 3×1011 Pn #/km using the NEDC drive cycle. The design and simulation tasks were conducted by JLR with support from AVL. The calibration and all of the vehicle testing was conducted by AVL, in Graz. Extensive design and CFD work was conducted to refine the inlet port, piston crown and injector spray pattern in order to reduce surface wetting and improve air to fuel mixing homogeneity. The design and CFD steps are detailed along with the results compared to target.

The European Emissions Stage 5b standard for diesel passenger cars regulates particulate matter to 0.0045 g/km and non-volatile part/km greater than 23 nm size to 6.0x10₁₁ as determined by the PMP procedure that uses a heated evaporation tube to remove semi-volatile material. Measurement artifacts associated with the evaporation tube technique prevents reliable extension of the method to a lower size range. Catalytic stripper (CS) technology removes possible sources of these artifacts by effectively removing all hydrocarbons and sulfuric acid in the gas phase in order to avoid any chemical reactions or re-nucleation that may cause measurement complications. The performance of a miniature CS was evaluated and experimental results showed solid particle penetration was 50% at 10.5 nm. The sulfate storage capacity integrated into the CS enabled it to chemically remove sulfuric acid vapor rather than rely on dilution to prevent nucleation.

Particulate emission reduction has long been a challenge for diesel engines as the diesel diffusion combustion process can generate high levels of soot which is one of the main constituents of particulate matter. Gasoline engines use a pre-mixed combustion process which produces negligible levels of soot, so particulate emissions have not been an issue for gasoline engines, particularly with modern port fuel injected (PFI) engines which provide excellent mixture quality. Future European and US emissions standards will include more stringent particulate limits for gasoline engines to protect against increases in airborne particulate levels due to the more widespread use of gasoline direct injection (GDI). While GDI engines are typically more efficient than PFI engines, they emit higher particulate levels, but still meet the current particulate standards.

A new fast and efficient parallel numerical solver for reacting and compressible flows through porous media has been developed in the OpenFOAM® (Open Field Operation and Manipulation) CFD Toolbox. With respect to the macroscopic model for porous media originally available in OpenFOAM®, a different mathematical approach has been followed: the new implemented solver makes use of the physical normal components resulting from the velocity expansion in the unit orthogonal vector basis to compute the Darcy pressure drop across the porous medium. Also, an additional sink term to account for the increased flow friction over the porous wall has been included into the momentum equation. In the new solver, the pressure correction equation is still able to achieve a faster convergency at very low permeability of the medium, also when it is associated with grid non-orthogonality.