Search Results

Transient engine operations are modeled and simulated with a 1-D code (GT Power) using heat release and emission data computed by a 3-D CFD code (Kiva3). During each iteration step of a transient engine simulation, the 1-D code utilizes the 3-D data to interpolate the values for heat release and emissions. The 3-D CFD computations were performed for the compression and combustion stroke of strategically chosen engine operating points considering engine speed, torque and excess air. The 3-D inlet conditions were obtained from the 1-D code, which utilized 3-D heat release data from the previous 1-D unsteady computations. In most cases, only two different sets of 3-D input data are needed to interpolate the transient phase between two engine operating points. This keeps the computation time at a reasonable level. The results are demonstrated on the load response of a generator which is driven by a medium-speed diesel engine.

Hydrogen engines are currently considered as a viable solution to preserve the internal combustion engine as a power unit for vehicle propulsion. In particular, lean-burn gasoline Spark-Ignition (SI) engines have been a major subject of investigations due to the reduced emission levels and high thermodynamic efficiency. This strategy is suitable for the purpose of passenger car applications and cannot be tailored in the field of high performance engine, where the air mass delivered would require oversized turbocharging systems or more complex charging solutions. For this reason, the range of stoichiometric feeding condition is explored in the high performance engine, leading to the consequent issue of abatement of pollutant emissions. In this work a 1D model will be applied to the modeling of a V8 engine fueled with DI of hydrogen. The engine has been derived by a gasoline configuration and adapted to hydrogen in such a way to keep the same performance.

The continuous pursuit of higher combustion efficiencies, as well as the possible usage of synthetic fuels with different properties than fossil-ones, require reliable and low-cost numerical approaches to support and speed-up engines industrial design. In this context, SI engines operated with homogeneous ultra-lean mixtures both characterized by a classical ignition configuration or equipped with an active prechamber represent the most promising solutions. In this work, for the classical ignition arrangement, a 3DCFD strategy to model the impact of the ignition system type on the CCV is developed using the RANS approach for turbulence modelling. The spark-discharge is modelled through a set of Lagrangian particles, whose velocity is modified with a zero-divergence perturbation at each discharge event, then evolved according to the Simplified Langevin Model (SLM) to simulate stochastic interactions with the surrounding gas flow.

The demand of game-changing technologies to improve efficiency and abate emissions of heavy-duty trucks and off-road vehicles promoted the development of novel engine concepts. The Recuperated Split-Cycle (R-SC) engine allows to recover the exhaust gases energy into the air intake by separating the compression and combustion stages into two different but connected cylinders: the compressor and expander, respectively. The result is a potential increase of the engine thermal efficiency. Accordingly, the 3D-computational fluid dynamics (CFD) modelling of the gas exchange process and the combustion evolution inside the expander becomes essential to control and optimize the R-SC engine concept. This work aims to address the most challenging numerical aspects encountered in a 3D numerical simulation of an R-SC engine.

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.

A 2-D CFD model was developed to describe the heat transfer, and reaction kinetics in a honeycomb structured ceramic diesel particulate trap. This model describes the steady state as well as the transient behavior of the flow and heat transfer during the trap regeneration processes. The trap temperature profile was determined by numerically solving the 2-D unsteady energy equation including the convective, heat conduction and viscous dissipation terms. The convective terms were based on a 2-D analytical flow field solution derived from the conservation of mass and momentum equations (Opris, 1997). The reaction kinetics were described using a discretized first order Arrhenius function. The 2-D term describing the reaction kinetics and particulate matter conservation of mass was added to the energy equation as a source term in order to represent the particulate matter oxidation. The filtration model describes the particulate matter accumulation in the trap.

Active prechamber combustion systems for SI engines represent a feasible and effective solution in reducing fuel consumption and pollutant emissions for both marine and ground heavy-duty engines. However, reliable and low-cost numerical approaches need to be developed to support and speed-up their industrial design considering their geometry complexity and the involved multiple flow length scales. This work presents a CFD methodology based on the RANS approach for the simulation of active prechamber spark-ignition engines. To reduce the computational time, the gas exchange process is computed only in the prechamber region to correctly describe the flow and mixture distributions, while the whole cylinder geometry is considered only for the power-cycle (compression, combustion and expansion). Outside the prechamber the in-cylinder flow field at IVC is estimated from the measured swirl ratio.

A further development of the ETAB atomization and drop breakup model for high pressure-driven liquid fuel jets, has been developed, tuned and validated. As in the ETAB model, this breakup model reflects a cascade of drop breakups, where the breakup criterion is determined by the Taylor drop oscillator and each breakup event resembles experimentally observed breakup mechanisms. A fragmented liquid core due to inner-nozzle disturbances is achieved by injecting large droplets subject to this breakup cascade. These large droplets are equipped with appropriate initial deformation velocities in order to obtain experimentally observed breakup lengths. In contrast to the ETAB model which consideres only the bag breakup or the stripping breakup mechanism, the new model has been extended to include the catastrophic breakup regime. In addition, a continuity condition on the breakup parameters has lead to the reduction of one model constant.

High fidelity engine simulation requires realistic fuel models. Although typical automotive fuels consist of more than few hundreds of hydrocarbon species, researches show that the physical and chemical properties of the real fuels could be represented by appropriate surrogate fuel models. It is desirable to represent the fuel using the same set of physical and chemical surrogate components. However, when the reaction mechanisms for a certain physical surrogate component is not available, the chemistry of the unmatched physical component is described using that of a similar chemical surrogate component at the expense of accuracy. In order to reduce the prediction error while maintaining the computational efficiency, a method of on-the-fly reactivity adjustment (ReAd) of chemical reaction mechanism along with fuel re-distribution based on reactivity is presented and tested in this study.

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.

Time-averaged temperatures at critical locations on the piston of a direct-fuel injected, two-stroke, 388 cm3, research engine were measured using an infrared telemetry device. The piston temperatures were compared to data [7] of a carbureted version of the two-stroke engine, that was operated at comparable conditions. All temperatures were obtained at wide open throttle, and varying engine speeds (2000-4500 rpm, at 500 rpm intervals). The temperatures were measured in a configuration that allowed for axial heat flux to be determined through the piston. The heat flux was compared to carbureted data [8] obtained using measured piston temperatures as boundary conditions for a computer model, and solving for the heat flux. The direct-fuel-injected piston temperatures and heat fluxes were significantly higher than the carbureted piston. On the exhaust side of the piston, the direct-fuel injected piston temperatures ranged from 33-73 °C higher than the conventional carbureted piston.

A correct prediction of the initial stages of the combustion process in SI engines is of great importance to understand how local flow conditions, fuel properties, mixture stratification and ignition affect the in-cylinder pressure development and pollutant formation. However, flame kernel growth is governed by many interacting processes including energy transfer from the electrical circuit to the gas phase, interaction between the plasma channel and the flow field, transition between different combustion regimes and gas expansion at very high temperatures. In this work, the authors intend to present a comprehensive, multi-dimensional model that can be used to predict the initial combustion stages in SI engines. In particular, the spark channel is represented by a set of Lagrangian particles where each one of them acts as a single flame kernel.

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.

The paper illustrates and validates a novel predictive combustion model for the estimation of performances and pollutant production in CI engines. The numerical methodology was developed by the authors for near real-time applications, while aiming at an accurate description of the air mixing process by means of a multi-zone approach of the air-fuel mass. Charge stratification is estimated via a 2D representation of the fuel spray distribution that is numerically derived by an axial one-dimensional control-volume description of the direct injection. The radial coordinate of each control volume is reconstructed a posteriori by means of a local distribution function. Fuel mass clustered in each zone is further split in ‘liquid’, ‘unburnt’ and ‘burnt’ sub-zones, given the local properties of the fuel spray control volumes with respect to space-time location of modelled ignition delay, liquid length, and flame lift-off.

In order to comply with 2002 EPA emissions regulations, cooled exhaust gas recirculation (EGR) will be used by heavy duty (HD) diesel engine manufacturers as the primary means to reduce emissions of nitrogen oxides (NOx). A feedforward controlled EGR cooling system with a secondary electric water pump and proportional-integral-derivative (PID) feedback has been designed to cool the recirculated exhaust gas in order to better realize the benefits of EGR without overcooling the exhaust gas since overcooling leads to the fouling of the EGR cooler with acidic residues. A system without a variable controlled coolant flow rate is not able to achieve these goals because the exhaust temperature and the EGR schedule vary significantly, especially under transient and warm-up operating conditions. Simulation results presented in this paper have been determined using the Vehicle Engine Cooling System Simulation (VECSS) software, which has been developed and validated using actual engine data.

Nowadays a great attention is paid to the level and quality of noise radiated from the tailpipe end of intake and exhaust systems, to control the gas dynamic noise emitted by the engine as well as the characteristics of the cabin interior sound. The muffler geometry can be optimized consequently, to attenuate or remark certain spectral components of the engine noise, according to the result expected. Evidently the design of complex silencing systems is a time-consuming operation, which must be carried out by means of concurrent experimental measurements and numerical simulations. In particular, 1D and multiD linear/non-linear simulation codes can be applied to predict the silencer behavior in the time and frequency domain. This paper describes the development of a 1D-multiD integrated approach for the simulation of complex muffler configurations such as reverse chambers with inlet and outlet pipe extensions and perforated silencers with the addition of sound absorbing material.

The present work describes the numerical modeling of medium-speed marine engines, operating in a fumigated dual-fuel mode, i.e. with the second fuel injected in the ports. This engine technology allows reducing engine-out emissions while maintaining the engine efficiency and can be fairly easily retrofitted from current diesel engines. The main premixed fuel that is added can be a low-carbon one and can additionally be of a renewable nature, thereby reducing or even completely removing the global warming impact. To fully optimize the operational parameters of such a large marine engine, computational fluid dynamics can be very helpful. Accurately describing the combustion process in such an engine is key, as the prediction of the heat release and the pollutant formation is crucial. Auto-ignition of the diesel fuel needs to be captured, followed by the combustion and flame propagation of the premixed fuel.

SI engines fueled with hydrogen represent a promising powertrain solution to meet the ambitious target of carbon-free emissions at the tailpipe. Therefore, fast and reliable numerical tools can significantly support the automotive industry in the optimization of such technology. In this work, a 1D-3D methodology is presented to simulate in detail the combustion process with minimal computational effort. First, a 1D analysis of the complete engine cycle is carried out on the user-defined powertrain configuration. The purpose is to achieve reliable boundary conditions for the combustion chamber, based on realistic engine parameters. Then, a 3D simulation of the power-cycle is performed to mimic the combustion process. The flow velocity and turbulence distributions are initialized without the need of simulating the gas exchange process, according to a validated technique.

Pursuing a sustainable energy scenario for transportation requires the blending of renewable oxygenated fuels such as alcohols into commercial hydrocarbon fuels. From a chemical kinetic perspective, this requires the accurate description of both hydrocarbon reference fuels (n-heptane, iso-octane, toluene, etc.) and oxygenated fuels chemistry. A recent systematic investigation of linear C2-C5 alcohols ignition in a rapid compression machine at p = 10-30 bar and T = 650- 900 K has extended the scarcity of fundamental data at such conditions, allowing for a revision of the low temperature chemistry for alcohol fuels in the POLIMI mechanism. Heavier alcohols such as n-butanol and n-pentanol present ignition characteristic of interest for application in HCCI engines, due to the presence of the hydroxyl moiety reducing their low temperature reactivity compared to the parent linear alkanes (i.e. higher octane number).

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.