Search Results

To promote advanced combustion strategies complying with stringent emission regulations of CI engines, computational models have to accurately predict the injector inner flow and cavitation development in the nozzle. This paper describes a coupled 1D/2D/3D modeling technique for the simulation of fuel flow and nozzle cavitation in diesel injection systems. The new technique comprises 1D fuel flow, 2D multi-body dynamics and 3D modeling of nozzle inner flow using a multi-fluid method. The 1D/2D model of the common rail injector is created with AVL software Boost-Hydsim. The computational mesh including the nozzle sac with spray holes is generated with AVL meshing tool Fame. 3D multi-phase calculations are performed with AVL software FIRE. The co-simulation procedure is controlled by Boost-Hydsim. Initially Hydsim performs a standalone 1D simulation until the needle lift reaches a prescribed tolerance (typically 2 to 5 μm).

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 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.

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.

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.

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.

The development of thermally durable zeolite NH3/Urea-SCR formulations coupled with that of high porosity filters substrates has opened the way to integrate PM and NOx control into a single device, namely an SCR-coated Diesel Particulate Filter (SCRF). A few experimental works are already present in the literature regarding SCRF systems, mainly addressing the DeNOx performances of the system (in both presence and absence of soot) under both steady state and transient conditions. The purpose of the present work is to perform a simulation study focused on phenomena which are expected to play key roles in SCRF systems, such as coupling of reaction and diffusion phenomena, soot effect on DeNOx activity, SCR coating effect on soot regeneration and filtration efficiency and competition between soot oxidation and DeNOx processes involving NO2.

Wall flow particulate filters have been used as a standard exhaust aftertreatment device for many years. The interaction of particulate matter (PM) regeneration and catalytically supported reactions strongly depends on the given operating conditions. Temperature, species concentration and mass flow cause a change from advective to diffusive-controlled flow conditions and influence the rate controlling dominance of individual reactions. A transient 1D+1D model is presented considering advective and diffusive transport phenomena. The reaction scheme focuses on passive PM conversion and catalytic oxidation of NO. The model is validated with analytical references. The impact of back-diffusion is explored simulating pure advective and combined advective diffusive species transport. Rate approaches from literature are applied to investigate PM conversion at various operating conditions.

Engine simulation can be performed using model approaches of different depths in capturing physical effects. The present paper presents a comprehensive comparison study on seven different engine models. The models range from transient 1D cycle resolved approaches to steady-state non-dimensional maps. The models are discussed in the light of key features, amount and kind of required input data, model calibration effort and predictability and application areas. The computational performance of the different models and their capabilities to capture different transient effects is investigated together with a vehicle model under real-life driving conditions. In the trade-off field of model predictability and computational performance an innovative approach on crank-angle resolved cylinder modeling turned out to be most beneficial.

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.

The demand for high efficiency powertrains in automotive engineering is further increasing, with hybrid powertrains being a feasible option to cope with new legislations. So far hybridization has only played a minor role for L-category vehicles. Focusing on an exemplary high-power L-category on-road vehicle, this research aims to show a new development approach, which combines longitudinal dynamic simulation (LDS) with “Design of Experiments” (DoE) in course of hybrid electric powertrain development. Furthermore, addressing the technological aspect, this paper points out how such a vehicle can benefit from 48V-hybridization of its already existing internal combustion powertrain. A fully parametric LDS model is built in Matlab/Simulink, with exchangeable powertrain components and an adaptable hybrid operation strategy. Beforehand, characterizing decisions as to focus on 48V and on parallel hybrid architecture are made.

In 2017 the European authorities put into effect the first part of a new certification test procedure for Real Driving Emissions (RDE). Similar tests are planned in other regions of the world, such as the upcoming China 6a/6b standards, further tightening emission limits, and also the introduction of RDE tests. Both restrictions pose challenging engineering tasks for upcoming vehicles. RDE certification tests feature significantly more demanding engine operating conditions and thus, emit more pollutants by orders of magnitude compared to known cycles like NEDC. Here, especially the reduction of NOx is a specific technical challenge, as it needs to compromise also with reduction targets on carbon dioxide. The fulfilment of both emission limits requires a widening of the focus from an isolated engine or exhaust aftertreatment view to a system engineering view involving all hardware and software domains of the vehicle.

Turbulence modeling for fuel spray simulation plays a prominent role in the understanding of the flow behavior in Internal Combustion Engines (ICEs). Currently, a lot of research work is actively spent on Large Eddy Simulation (LES) turbulence modeling as a replacement option of standard Reynolds averaged approaches in the Eulerian-Lagrangian spray modeling framework, due to its capability to accurately describe flow-induced spray variability and to the lower dependence of the results on the specific turbulence model and/or modeling coefficients. The introduction of LES poses, however, additional questions related to the implementation/adaptation of spray-related turbulence sources and to the rise of conflicting numerics and grid requirements between the Lagrangian and Eulerian parts of the simulated flow.

This work describes the development of a computational model for the CFD simulation of compact heat exchangers applied for the oil cooling in internal combustion engines. Among the different cooler types, the present modeling effort will be focused on liquid-cooled solutions based on offset strip fins turbulators. The design of this type of coolers represents an issue of extreme concern, which requires a compromise between different objectives: high compactness, low pressure drop, high heat-transfer efficiency. In this work, a computational framework for the CFD simulation of compact oil-to-liquid heat exchangers, including offset-strip fins as heat transfer enhancer, has been developed. The main problem is represented by the need of considering different scales in the simulation, ranging from the characteristic size of the turbulator geometry (tipically μm - mm) to the full scale of the overall device (typically cm - dm).

Non-combustible particles, commonly summarized as ash, influence the lifetime performance of wall flow filters. This study aims to investigate this influence by means of simulation. An existing transient 1D+1D wall flow filter model is extended by dedicated transport balances for soot and ash (1), by a discrete cake model describing changing soot and ash compositions over the cake height (2), by a phenomenological cake filtration model (3), by dedicated cake property models (4) and by a phenomenological model capturing the radial mobility of solids within the cake (5). Results of three different types of simulations are shown. First, the various sub-models are assessed in isolated simulation configurations. The combination of these shall serve as theoretical model validation. Second, isolated loading and passive regeneration simulations are performed.