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.

In the next years, the upcoming emission legislations are expected to introduce further restrictions on the admittable level of pollutants from vehicles measured on homologation cycles and real drive tests. In this context, the strict control of pollutant emissions at the cold start will become a crucial point to comply with the new regulation standards. This will necessarily require the implementation of novel strategies to speed-up the light-off of the reactions occurring in the after-treatment system, since the cold start conditions are the most critical one for cumulative emissions. Among the different possible technological solutions, this paper focuses on the evaluation of the potential of a burner system, which is activated before the engine start. The hypothetical burner exploits the lean combustion of an air-gasoline mixture to generate a high temperature gas stream which is directed to the catalyst section promoting a fast heating of the substrate.

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

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 OpenFOAM® CFD methodology is nowadays employed for simulation in internal combustion engines and a lot of work has been done for an appropriate description of all complex phenomena. At the moment in the RANS turbulence models available in the OpenFOAM® toolbox the turbulence modulation is not yet included, and the present work analyzes the predictive capabilities of the code in simulating high injection pressure fuel sprays after modeling the influence of the dispersed phase on the turbulence structure. Different experiments were employed for the validation. At first, non-evaporating diesel spray was considered in a constant volume and quiescent vessel. The validation was performed via the available experimental spray evolution in terms of penetrations and spatial/temporal fuel distributions. Then the Sandia combustion chamber was chosen for diesel spray simulation in non-reacting conditions.

The shift toward electrification and limitations in battery electric vehicle technology have led to high demand for hybrid vehicles (HEVs) that utilize a battery and an internal combustion engine (ICE) for propulsion. Although HEVs enable lower fuel consumption and emissions compared to conventional vehicles, they still require combustion of fuels for ICE operation. Thus, emissions from hybrid vehicles are still a major concern. Engine starts are a major source of emissions during any driving event, especially before the three-way catalyst (TWC) reaches its light-off temperature. Since the engine is subjected to multiple starts during most driving events, it is important to mitigate and better understand the impact of these emissions. In this study, experiments were conducted to analyze engine starts in a hybrid powertrain on different experimental setup.

The s.i. combustion process and its corresponding pollutant formation are investigated by means of a quasiD approach and a CFD model. This work has been motivated by the need to better understand the reliability of such models and to assess their accuracies with respect to the prediction of engine performances and emissions. An extended dissertation about the fundamental mechanisms governing the pollutant formation in the turbulent premixed combustion which characterizes the s.i. engines is given. The conclusion of such analysis is the definition of a new reduced chemical scheme, based on the application of partial-equilibrium and steady-state assumptions for the radicals and the solution of a transport equation for each specie which is kinetically controlled. For this purpose the CFD code OpenFOAM [1, 2, 3] and the thermo-fluid dynamic code GASDYN [4, 5] have been applied and enhanced.

Aerodynamic properties of a BMW car model taking over a truck model are studied computationally by applying the scale-resolving PANS (Partially-averaged Navier-Stokes) approach. Both vehicles represent down-scaled (1:2.5), geometrically-similar models of realistic vehicle configurations for which on-road measurements have been performed by Schrefl (2008). The operating conditions of the modelled ‘on-road’ overtaking maneuver are determined by applying the dynamic similarity concept in terms of Reynolds number consistency. The simulated vehicle configuration constitutes of a non-moving truck model and a car model moving against the air flow, the velocity of which corresponds to the car velocity.

Ducted Fuel Injection (DFI) has the potential to reduce soot emissions in Diesel engines thanks to the enhanced mixing rate resulting from the liquid fuel flow through a small cylindrical pipe located at a certain distance from the nozzle injector hole. A consolidated set of experiments in constant-volume vessel and engine allowed to understand the effects of ambient conditions, duct geometry and shape on fuel-air mixing, combustion and soot formation. However, implementation of this promising technology in compression-ignition engines requires predictive numerical models that can properly support the design of combustion systems in a wide range of operating conditions. This work presents a computational methodology to predict fuel-air mixing and combustion with ducted fuel injection. Attention is mainly focused on turbulence and combustion modelling.

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.

The combustion efficiency of large gas engines is limited by knocking combustion. Due to fact that the quality of the fuel gas has a high impact on the self-ignition of the mixture, it is the aim of this work to model the knocking combustion for fuel gases with different composition using detailed chemistry. A cycle-resolved knock simulation of the fast burning cycles was carried out in order to assume realistic temperatures and pressures in the unburned mixture Therefore, an empirical model that predicts the cyclic variations on the basis of turbulent and chemical time scales was derived from measured burn rates and implemented in a 1D simulation model. Based on the simulation of the fast burning engine cycles the self-ignition process of the unburned zone was calculated with a stochastic reactor model and correlated to measurements from the engines test bench. A good agreement of the knock onset could be achieved with this approach.

The aim of this paper is to present the adjoint equations for shape optimization derived from steady incompressible Navier-Stokes (N-S) equations and an objective functional. These adjoint Navier-Stokes equations have a similar form as the N-S equations, while the source terms and the boundary conditions depend on the chosen objective. Additionally, the gradient of the targeted objective with respect to the design variables is calculated. Based on this, a modification of the geometry is computed to arrive at an improved objective value. In order to find out, whether a more sophisticated approach is needed, the adjoint equations are derived by using two different approaches. The first approach is based on the frozen turbulence assumption and the second approach, which is advanced in this paper, is derived from the near wall k − ζ − f turbulence model.

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.

The definition of a robust methodology to perform a full-cycle CFD simulation of IC engines requires as first step the availability of a reliable grid generation tool, which does not only have to guarantee a high quality mesh but also has to prove to be efficient in terms of required time. In this work the authors discuss a novel approach entirely based on the OpenFOAM technology, in which the available 3D grid generator was employed to automatically create meshes containing hexahedra and split-hexahedra from triangulated surface geometries in Stereolithography (STL) format. The possibility to introduce local refinements and boundary layers makes this tool suitable for IC engine simulations. Grids are sequentially generated at target crank angles which are automatically determined depending on user specified settings such as maximum mesh validity interval and quality parameters like non-orthogonality, skewness and aspect ratio.