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.

Nowadays quasi-3D approaches are included in many commercial and research 1D numerical codes, in order to increase their simulation accuracy in presence of complex shape 3D volumes, e.g. plenums and silencers. In particular, these are regarded as valuable approaches for application during the design phase of an engine, for their capability of predicting non-planar waves motion and, on the other hand, for their low requirements in terms of computational runtime. However, the generation of a high-quality quasi-3D computational grid is not always straightforward, especially in case of complex elements, and can be a time-consuming operation, making the quasi-3D tool a less attractive option. In this work, a quasi-3D module has been implemented on the basis of the open-source CFD code OpenFOAM and coupled with the 1D code GASDYN.

ELEVATE (European Low Emission V4 Automotive Two-stroke Engine) was a research project part funded by the European Commission to design and develop a compact and efficient gasoline two-stroke automotive engine. Five partners were involved in the project, IFP (Institut Français Du Pétrole) who were the project leaders, Lotus, Opcon (Autorotor and SEM), Politecnico di Milano and Queen's University Belfast. The general project targets were to achieve Euro 3 emissions compliance without DeNOx catalisation, and a power output of 120 kW at 5000 rev/min with maximum torque of 250 Nm at 2000 rev/min. Specific targets were a 15% reduction in fuel consumption compared to its four-stroke counterpart and a size and weight advantage over the four-stroke diesel with significant reduction in particulate and NOx emissions. This paper describes the design philosophy of the engine as well as the application of the various partner technologies used.

Partially premixed combustion (PPC) can be applied to decrease emissions and increase fuel efficiency in direct injection, compression ignition (DICI) combustion engines. PPC is strongly influenced by the mixing of fuel and oxidizer, which for a given fuel is controlled mainly by (a) the fuel injection, (b) the in-cylinder flow, and (c) the geometry and dynamics of the engine. As the injection timings can vary over a wide range in PPC combustion, detailed knowledge of the in-cylinder flow over the whole intake and compression strokes can improve our understanding of PPC combustion. In computational fluid dynamics (CFD) the in-cylinder flow is sometimes simplified and modeled as a solid-body rotation profile at some time prior to injection to produce a realistic flow field at the moment of injection. In real engines, the in-cylinder flow motion is governed by the intake manifold, the valve motion, and the engine geometry.

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.

This paper provides a review on state-of-art modern cooling systems employed for thermal cooling of electric motors for vehicle applications. In recent years, the pursue of a more sustainable and ecofriendly mobility has pushed the research towards the development of electric vehicle powertrain systems. Besides the evident advantages of the adoption of electric traction systems in terms of pollution and efficiency, the need of an effective cooling system for the electric machine components gained more and more importance in order to maintain high efficiency and ensure high durability. In fact, it is known that high temperatures can be harmful for the electric motor: besides the evident damages for mechanical parts, the influence on the permanent magnet properties is not negligible [1] [2]. In this fast-evolving environment, different solutions for the thermal problem have been researched and adopted, each one with its own pros and cons.

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.

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.

Multi-purpose agricultural tractors are vehicles that are usually used in rough paths and on off-road situations characterized by strong slope variations. The main feature of this kind of vehicles is the stability in working conditions to avoid overturning while it is on duty. This characteristic is given by the interaction between the suspension system and the vehicle frame. Due to the limited size of this kind of vehicle, the stability feature could be given by chassis deformation or using a two-piece frame connected by a spherical joint. This paper presents the validation of a numerical lumped-parameters model able to reproduce the vertical dynamics of a multi-purpose tractor featured by a yielding chassis. The unknown model parameters have been estimated firstly with static tests to study the vertical tire and suspension stiffnesses. The dynamic tests using a four-post-test rig have been performed to tune the unknown dynamic parameters.

The purpose of the ELEVATE (European Low Emission V4 Automotive Two-stroke Engine) industrial research project is to develop a small, compact, light weight, high torque and highly efficient clean gasoline 2-stroke engine of 120 kW which could industrially replace the relatively big existing automotive spark ignition or diesel 4-stroke engine used in the top of the mid size or in the large size vehicles, including the minivan vehicles used for multi people and family transportation. This new gasoline direct injection engine concept is based on the combined implementation on a 4-stroke bottom end of several 2-stroke engine innovative technologies such as the IAPAC compressed air assisted direct fuel injection, the CAI (Controlled Auto-Ignition) combustion process, the D2SC (Dual Delivery Screw SuperCharger) for both low pressure engine scavenging and higher pressure IAPAC air assisted DI and the ETV (Exhaust charge Trapping Valve).

We present a framework for the robust optimization of the heat flux distribution for an anti-ice electro-thermal ice protection system (AI-ETIPS) and iced airfoil performance analysis under uncertain conditions. The considered uncertainty regards a lack of knowledge concerning the characteristics of the cloud i.e. the liquid water content and the median volume diameter of water droplets, and the accuracy of measuring devices i.e., the static temperature probe, uncertain parameters are modeled as uniform random variables. A forward uncertainty propagation analysis is carried out using a Monte Carlo approach. The optimization framework relies on a gradient-free algorithm (Mesh Adaptive Direct Search) and three different problem formulations are considered in this work. Two bi-objective deterministic optimizations aim to minimize power consumption and either minimize ice formations or the iced airfoil drag coefficient.

Very little is currently known of the aerodynamic interaction between neighboring cycloidal rotors. Such knowledge is, however, of crucial importance to tune the controller and rotor disposition of a cyclogiro aircraft. Thus, a three-dimensional computational fluid dynamics (CFD) model is developed, validated, and used to analyze the D-Dalus L1 four-rotor unmanned aircraft operating under several configurations. The model solves the Euler equations using the OpenFOAM toolbox in order to provide fast results on a desktop computer. Validation is performed against thrust forces and flow streamlines obtained during wind tunnel experiments at various flight velocities. Numerical results from CFD match the trends of the experimental data. Flow behavior matches the video footage of the wind tunnel tests. Although boundary layer effects are neglected, satisfactory results are obtained both qualitatively and quantitatively.

Prediction of in-cylinder flows and fuel-air mixing are two fundamental pre-requisites for a successful simulation of direct-injection engines. Over the years, many efforts were carried out in order to improve available turbulence and spray models. However, enhancements in physical modeling can be drastically affected by how the mesh is structured. Grid quality can negatively influence the prediction of organized charge motion structures, turbulence generation and interaction between in-cylinder flows and injected sprays. This is even more relevant for modern direct injection engines, where multiple injections and control of charge motions are employed in a large portion of the operating map. Currently, two different approaches for mesh generation exist: manual and automatic. The first makes generally possible to generate high-quality meshes but, at the same time, it is very time consuming and not completely free from user errors.

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.

In this study two different simulation approaches to large eddy simulation of spark-ignition engines are compared. Additionally, some of the simulation results are compared to experimentally obtained in-cylinder velocity measurements. The first approach applies unstructured grids with an automated meshing procedure, using OpenFoam and Lib-ICE with a mapping approach. The second approach applies the efficient in-house code PsiPhi on equidistant, Cartesian grids, representing walls by immersed boundaries, where the moving piston and valves are described as topologically connected groups of Lagrangian particles. In the experiments, two-dimensional two-component particle image velocimetry is applied in the central tumble plane of the cylinder of an optically accessible engine. Good agreement between numerical results and experiment are obtained by both approaches.

The increasing diffusion of gasoline direct injection (GDI) engines requires a more detailed and reliable description of the phenomena occurring during the fuel injection process. As well known the thermal and fluid-dynamic conditions present in the combustion chamber greatly influence the air-fuel mixture process deriving from GDI injectors. GDI fuel sprays typically evolve in wide range of ambient pressure and temperatures depending on the engine load. In some particular injection conditions, when in-cylinder pressure is relatively low, flash evaporation might occur significantly affecting the fuel-air mixing process. In some other particular injection conditions spray impingement on the piston wall might occur, causing high unburned hydrocarbons and soot emissions, so currently representing one of the main drawbacks of GDI engines.

Since vehicles are comprised of thousands of components, it is essential to reduce the Life Cycle Inventory (LCI) modelling workload. This study aims to compare different LCI modeling workload-reducing scenarios to provide a trade-off between the workload efforts and result accuracy. To achieve the optimal balance between computational effort and data specification requirements, the driver seat is used as a case study, instead of the entire vehicle. When all the components of a conventional light-duty commercial vehicle are sorted by mass descending order, seats are among the first five. In addition, unlike the other components, seats are comprised of metals as well as a wide range of plastics and textiles, making them a representative test case for a general problem formulation. In this way, methodology and outcomes can be reasonably extended to the entire vehicle.