Search Results

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.

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.

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.

Today, Direct-Injection systems are widely used on Spark-Ignition engines in combination with turbo-charging to reduce the fuel-consumption and the knock risks. In particular, the spread of Gasoline Direct Injection (GDI) systems is mainly related to the use of new generations of multi-hole, high-pressure injectors whose characteristics are quite different with respect to the hollow-cone, low-pressure injectors adopted in the last decade. This paper presents the results of an experimental campaign conducted on the spray produced by a GDI six-holes injector into a constant volume vessel with optical access. The vessel was filled with air at atmospheric pressure. Different operating conditions were considered for an injection pressure ranging from 3 to 20 MPa. For each operating condition, spray images were acquired by a CCD camera and then post processed to evaluate the spray penetration and cone angles.

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.

Future human exploration roadmaps involve the development of temporary or permanent outposts on Moon and Mars. The capability of providing astronauts with proper conditions for living and working in extraterrestrial environments is therefore a key issue for the sustainability of those roadmaps, and closed-loop Environment Control and Life Support Systems (ECLSSs) and bio-regenerative plants represent the necessary evolution of current technologies for complying with the challenging requirements imposed. This paper presents the architectural design of a terrestrial plant to be exploited to test and validate air and water management technologies for a biological life support system in a closed environment. The plant includes a crew area and a plant growth area. These two spaces can be considered as either a unique volume or two separated environments with reduced contact, e.g. for plant harvesting or other up-keeping activities.

In this paper, an experimental and numerical analysis of combustion process and knock occurrence in a small displacement spark-ignition engine is presented. A wide experimental campaign is preliminarily carried out in order to fully characterize the engine behavior in different operating conditions. In particular, the acquisition of a large number of consecutive pressure cycle is realized to analyze the Cyclic Variability (CV) effects in terms of Indicated Mean Effective Pressure (IMEP) Coefficient of Variation (CoV). The spark advance is also changed up to incipient knocking conditions, basing on a proper definition of a knock index. The latter is estimated through the decomposition and the FFT analysis of the instantaneous pressure cycles. Contemporary, a quasi-dimensional combustion and knock model, included within a whole engine one-dimensional (1D) modeling framework, are developed. Combustion and knock models are extended to include the CV effects, too.

A swirl injector for GDI application was used to inject an iso-octane spray in a quiescent chamber, to study the effect of the air density on the spray behavior. Stroboscopic images are recorded at different delays from the injection trigger to study the spray shape and structure. The temporal evolution of different spray parameters, length, width, angle, volume, instantaneous global air-fuel ratio, is calculated from the images. The effect of the increasing air density is to shorten the time and length scale of the spray evolution.

The ride comfort of three Alfa Romeo cars, namely Giulietta (1955), Alfetta (1972) and 159 (2005) has been assessed both objectively and subjectively. The three cars belong to the same market segment. The aim is to let young engineers or graduate students understand how technology has evolved and eventually learn a lesson from the assessed trend. A number of cleats have been fixed at the ground and the three cars have traversed such uneven surface. The objective assessment of the ride comfort has been performed by means of accelerometers fixed at the seat rails, additionally a special dummy developed at Politecnico di Milano has been employed. The subjective assessment has been performed by a panel of passengers. The match between objective and subjective ratings is very good. Simple mathematical models have been employed to establish a (successful) comparison between experimental and computational results. The ride comfort differs substantially among the cars.

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.

Detailed chemistry represents a fundamental pre-requisite for a realistic simulation of combustion process in Diesel engines to properly reproduce ignition delay and flame structure (lift-off and soot precursors) in a wide range of operating conditions. In this work, the authors developed reduced mechanisms for n-dodecane starting from the comprehensive kinetic mechanism developed at Politecnico di Milano, well validated and tested in a wide range of operating conditions [1]. An algorithm combining Sensitivity and Flux Analysis was employed for the present skeletal reduction. The size of the mechanisms can be limited to less than 100 species and incorporates the most important details of low-temperature kinetics for a proper prediction of the ignition delay. Furthermore, the high-temperature chemistry is also properly described both in terms of reactivity and species formation, including unsaturated compounds such as acetylene, whose concentration controls soot formation.

In the most elementary treatment of plane-wave reflection at the open end of a duct system, it is often assumed that the ends are pressure nodes. This implies that pressure is assumed as a constant at the open end termination and that steady flow boundary condition is supposed as instantaneously established. While this simplifying assumption seems reasonable, it does not consider any radiation of acoustic energy from the duct into the surrounding free space; hence, an error in the estimation of the effects of the flow on the acoustical response of an open-end duct occurs. If radiation is accounted, a complicated three-dimensional wave pattern near the duct end is established, which tends to readjust the exit pressure to its steady-flow level. This adjustment process is continually modified by further incident waves, so that the effective instantaneous boundary conditions which determine the reflected waves depend on the flow history.

In the paper several application techniques of MonteCarlo (MC) method applied to thermal analysis of space vehicles are presented. Although these methods are widely used in other engineering domains, their introduction to the thermal one is quite recent and not fully developed in the industrial practice. This paper aims at showing that, even without demanding computation resources (all what presented has been obtained with a single processor PC) MonteCarlo analysis techniques, in a preliminary design phase, can support and integrate engineering judgment of the thermal designer. In particular, it is exploited the applicability of the method to reduced thermal models, with a clear advantage in terms of computation time. An original approach is proposed, and results are shown. The papers shows the applicability of the MC method to the case when experimental data of the uncertain parameters are available, using the bootstrap re-sampling techniques.

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.

Wind turbines in cold climates are likely to suffer from icing events, deteriorating the aerodynamic performances of the blades and decreasing their power output. Continuous ice accretion causes an increase in the ice mass and, consequently, in the centrifugal force to which the ice shape is subjected. This can result in the shedding of chunks of ice, which can jeopardize the aeroelastic properties of the blade and, most importantly, the safety of the surrounding people and of the wind turbine structure itself. In this work, ice shedding analysis is performed on a quasi-3D, multi-step ice geometry accreted on the NREL 5MW reference wind turbine. A preliminary investigation is performed by including the presence of an ice protection system to decrease the adhesion surface of the ice on the blade. A reference test case with a simple geometry is used as verification for the correct implementation of the procedure.

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.

This study examines the impact of slip in aero-thermal conditions of supercooled large droplets (SLD) produced in an Icing Wind Tunnel (IWT) on the ice accretion characteristics. The study identifies potential biases in the SLD model development based on IWT data and numerical predictions that assume the SLD to be in aerothermal equilibrium with the IWT airflow. To obtain realistic temperature and velocity data for each droplet size class in the test section of the Braunschweig Icing Wind Tunnel (BIWT), a Lagrangian droplet tracking solver was used within a Monte Carlo framework. Results showed that SLDs experience considerable slips in velocity and temperature due to their higher inertia and short residence time in the Braunschweig IWT. Large droplets were found to be warmer and slower than the flow in the test section, with larger droplets experiencing larger aerothermal slips.

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.

Aqueous Urea is a non-toxic and stable ammonia carrier and its injection and mixing represent the basis for the most common de-NOx technology for mobile applications. The reactant feed preparation process is defined by evaporation, thermolysis and hydrolysis of the liquid mixture upstream the Selective Catalytic Reduction reactor, and it is strongly dependent on the interaction between spray and gaseous flow. Low-pressure driven injectors are the common industrial standard for these applications, and their behavior in almost-ambient pressure cross flows is significantly different from any in-cylinder application. For this reason, two substantially different injectors in terms of geometry and design are experimentally studied, characterizing drop sizes and velocities through Phase Doppler Anemometry (PDA) and liquid mass spatial distribution through Shadow Imaging (SI).