Search Results

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.

Metallic foams or sponges are materials with a cell structure suitable for many industrial applications, such as reformers, heat catalytic converters, etc. The success of these materials is due to the combination of various characteristics such as mechanical strength, low density, high specific surface, good thermal exchange properties, low flow resistance and sound absorption. Different materials and manufacturing processes produce different type of structure and properties for various applications. In this work a genetic algorithm has been developed and applied to support the design of catalytic devices. In particular, two substrates were considered, namely the traditional honeycomb and an alternative open-cell foam type. CFD simulations of pressure losses and literature based correlations for the heat and mass transfer were used to support the genetic algorithm in finding the best compromise between flow resistance and pollutant abatement.

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.

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.

Objective of this work is the incorporation of the flame stretch effects in an Eulerian-Lagrangian model for premixed SI combustion in order to describe ignition and flame propagation under highly inhomogeneous flow conditions. To this end, effects of energy transfer from electrical circuit and turbulent flame propagation were fully decoupled. The first ones are taken into account by Lagrangian particles whose main purpose is to generate an initial burned field in the computational domain. Turbulent flame development is instead considered only in the Eulerian gas phase for a better description of the local flow effects. To improve the model predictive capabilities, flame stretch effects were introduced in the turbulent combustion model by using formulations coming from the asymptotic theory and recently verified by means of DNS studies. Experiments carried out at Michigan Tech University in a pressurized, constant-volume vessel were used to validate the proposed approach.

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.

Computational fluid dynamics represents a useful tool to support the design and development of Heavy Duty Engines, making possible to test the effects of injection strategies and combustion chamber design for a wide range of operating conditions. Predictive models are required to ensure accurate estimations of heat release and the main pollutant emissions within a limited amount of time. For this reason, both detailed chemistry and turbulence chemistry interaction need to be included. In this work, the authors intend to apply combustion models based on tabulated kinetics for the prediction of Diesel combustion in Heavy Duty Engines. Four different approaches were considered: well-mixed model, presumed PDF, representative interactive flamelets and flamelet progress variable. Tabulated kinetics was also used for the estimation of NOx emissions.

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.

A detailed understanding of Gasoline Direct Injection (GDI) techniques applied to spark-ignition (SI) engines is necessary as they allow for many technical advantages such as increased power output, higher fuel efficiency and better cold start performances. Within this context, the extensive validation of multi-dimensional models against experimental data is a fundamental task in order to achieve an accurate reproduction of the physical phenomena characterizing the injected fuel spray. In this work, simulations of different Engine Combustion Network (ECN) Spray G conditions were performed with the Lib-ICE code, which is based on the open source OpenFOAM technology, by using a RANS Eulerian-Lagrangian approach to model the ambient gas-fuel spray interaction.

The automotive industry is striving to adopt model-based engine design and optimization procedures to reduce development time and costs. In this scenario, first-principles gas dynamic models predicting the mass, energy and momentum transport in the engine air path system with high accuracy and low computation effort are extremely important today for performance prediction, optimization and cylinder charge estimation and control. This paper presents a comparative study of two different modeling approaches to predict the one-dimensional unsteady compressible flow in the engine air path system. The first approach is based on a quasi-3D finite volume method, which relies on a geometrical reconstruction of the calculation domain using networks of zero-dimensional elements. The second approach is based on a model-order reduction procedure that projects the nonlinear hyperbolic partial differential equations describing the 1D unsteady flow in engine manifolds onto a predefined basis.

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.

Increasing demands on the capabilities of engine simulation and the ability to accurately predict both performance and acoustics has lead to the development of several numerical tools to help engine manufacturers during the prototyping stage. The aid of CFD tools (3D and 1D) can remarkably reduce the duration and the costs of this stage. The need of achieving good accuracy, along with acceptable computational runtime, has given the spur to the development of a geometry based quasi-3D approach. This is designed to model the acoustics and the fluid dynamics of both intake and exhaust system components used in internal combustion engines. Models of components are built using a network of quasi-3D cells based primarily on the geometry of the system. The solution procedure is based on an explicitly time marching staggered grid approach making use of a flux limiter to prevent numerical instabilities.

This work has the objective to present the extension of a novel quasi-dimensional model, developed to simulate the combustion process in diesel Compression Ignition (CI) engines, to describe this process when Dimethyl ether (DME) is used as fuel. DME is a promising fuel in heavy-duty CI engines application thanks to its high Cetane Number (CN), volatility, high reactivity, almost smokeless combustion, lower CO2 emission and the possibility to be produced with renewable energy sources. In this paper, a brief description of the thermodynamic model will be presented, with particular attention to the implementation of the Tabulated Kinetic Ignition (TKI) model, and how the various models interact to simulate the combustion process. The model has been validated against experimental data derived from constant-volume DME combustion, in this case the most important parameters analyzed and compared were the Ignition Delay (ID) and Flame Lift Off Length (FLOL).

Turbocharging is playing today a fundamental role not only to improve automotive engine performance, but also to reduce fuel consumption and exhaust emissions for both Spark Ignition and Diesel engines. Dedicated experimental investigations on turbochargers are therefore necessary to assess a better understanding of its performance. The availability of experimental information on turbocharger steady flow performance is an essential requirement to optimize the engine-turbocharger matching, which is usually achieved by means of simulation models. This aspect is even more important when referred to the turbine efficiency, since its swallowing capacity can be accurately evaluated through the measurement of mass flow rate, inlet temperature and pressure ratio across the machine.

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.

Today, multi-hole Diesel injectors can be mainly characterized by three different nozzle hole shapes: cylindrical, k-hole, and ks-hole. The nozzle hole layout plays a direct influence on the injector internal flow field characteristics and, in particular, on the cavitation and turbulence evolution over the hole length. In turn, the changes on the injector internal flow correlated to the nozzle shape produce immediate effects on the emerging spray. In the present paper, the fluid dynamic performance of three different Diesel nozzle hole shapes are evaluated: cylindrical, k-hole, and ks-hole. The ks-hole geometry was experimentally characterized in order to find out its real internal shape. First, the three nozzle shapes were studied by a fully transient CFD multiphase simulation to understand their differences in the internal flow field evolutions. In detail, the attention was focused on the turbulence and cavitation levels at hole exit.

The accurate prediction of pollutant emissions generated by IC engines is a key aspect to guarantee the respect of the emission regulation legislation. This paper describes the approach followed by the authors to achieve a strict numerical coupling of two different 1D modeling tools in a co-simulation environment, aiming at a reliable calculation of engine-out and tailpipe emissions. The main idea is to allow an accurate 1D simulation of the unsteady flows and wave motion inside the intake and exhaust systems, without resorting to an over-simplified geometrical discretization, and to rely on advanced thermodynamic combustion models and kinetic sub-models for the calculation of cylinder-out emissions. A specific fluid dynamic approach is then used to track the chemical composition along the exhaust duct-system, in order to evaluate the conversion efficiency of after-treatment devices, such as TWC, GPF, DPF, DOC, SCR and so on.

In this work the 3Dcell method, a quasi3D approach developed by the Internal Combustion Engine Group at Politecnico di Milano, has been extended and applied to the fluid dynamic simulation of turbocharging devices for internal combustion engines, focusing on the compressor side. The 3Dcell is based on a pseudo-staggered leapfrog method applied to the governing equation of a 1D problem arbitrarily oriented in space. The system of equations is solved referring to the relative system in the rotating zone, whereas the absolute reference system has been used elsewhere. The vaneless diffuser has been modelled resorting to the conservation of the angular momentum of the flow stream in the tangential direction, combined with the solution of the momentum equation in the radial direction.

Gasoline, direct injection engines represent one of the most widely adopted powertrain for passenger cars. However, further development efforts are necessary to meet the future fuel consumption and emission standards imposing an efficiency increase and a reduction of particulate matter emissions. Within this context, computational fluid dynamics is nowadays a consolidated tool to support engine design; this work is focused on the development of a set of CFD models for the prediction of combustion in modern GDI engines. The one-equation Weller model coupled with a zero-dimensional approach to handle initial flame kernel growth was applied to predict flame propagation. To account for mixture fraction fluctuations which might lead to the presence of soot precursor species, burned gas chemical composition is computed using tabulated kinetics with a presumed probability density function.