Search Results

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.

In the framework of an increasing demand for a more sustainable mobility, where the fuel consumption reduction is a key driver for the development of innovative internal combustion engines, Turbulent Jet Ignition (TJI) represents one of the most promising solutions to improve the thermal efficiency. However, details concerning turbulent jet assisted combustion are still to be fully captured, and therefore the design and the calibration of efficient TJI systems require the support of reliable simulation tools that can provide additional information not accessible through experiments. To this aim, an experimental investigation combined with a 3D-CFD study was performed to analyze the TJI combustion characteristics in a single-cylinder spark-ignition (SI) engine. Firstly, the model was validated against experiments considering stoichiometric mixture at 3000 rpm, wide open throttle operating conditions.

The results of an experimental study on the statistical structure of turbulence in an automotive engine are reported, with specific reference to the time-frequency domains. Autocorrelation and autospectral density coefficients were evaluated in consecutive crank-angle intervals throughout the induction and compression strokes. Eulerian time scales were obtained on the analogy of both the micro and integral time scales of turbulence for stationary flows. The spatial distribution of the turbulence structure was investigated in the combustion chamber of a diesel engine with a shallow in-piston bowl and two tangential intake ducts. The study was carried out for different swirl flow conditions, produced by deactivating one intake duct and/or by changing the engine speed. The velocity data were acquired using an advanced HWA technique, under motored conditions.

Computational Fluid Dynamics is a powerful instrument for PEM fuel cell systems development, testing and optimization. Considering the complication due to the multiple physical phenomena involved in the cell's operations, a good understanding of the micro-scale fluidic behavior in boundary layers is recommended: pressure drop along the reactants gas channels and the cooling channels has a sensible effect on parasite load in fuel cell systems (i.e. the power absorbed by the pump supplying the gases), as well as an important role in thermal transport. A correct thermal and fluid dynamic boundary layer prediction on the channel walls and the other contact surface with porous layers requires usually a dense finite element volumes discretization near wall, especially if laminar flows occur: therefore, the boundary layer computational cost tends to be the major one.

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.

The variation of turbulent flow quantities with engine speed has been investigated in the combustion chamber of an automotive diesel engine with a high-squish conical-type in-piston bowl and one helicoidal intake duct, at speeds covering the wide range of 600-3000 rpm, under motored conditions. The investigation had the main purpose of studying the engine speed effect on the structure of both cycle-resolved and conventional turbulence over the induction, the compression and the early stage of the expansion stroke. The low frequency component of the fluctuating motion was also investigated.

A new test bench has been set up and equipped in order to analyze the air mean motion and turbulence quantities in the combustion system of an automotive diesel engine with one helicoidal intake duct and a conical type in-piston bowl. A sophisticated HWA technique employing single- and dual-sensor probes was applied to the in-cylinder flow investigation under motored conditions. The anemometric probe was also operated as a thermometric sensor. An analytical-numerical procedure, based on the heat balance equations for both anemometric and thermometric wires, was refined and applied to compute the gas velocity from the anemometer output signal. The gas property influence, the thermometric sensor lag and the prong temperature effects were taken into account with this procedure. The in-cylinder velocity data were reduced using both a cycle-resolved approach and the conventional ensemble-averaging procedure, in order to separate the mean flow from the fluctuating motion.

The frequency spectral structure of turbulence spatial components was investigated in the cylinder of an automotive diesel engine with a high-squish reentrant in-piston bowl of the conical type and a helical inlet port. A sophisticated HWA technique using single- and dual-sensor probes was applied for instantaneous air velocity measurements along the injector axis at practical engine speeds, up to 3000 rpm, under motored conditions. The investigation was carried out for both cycle-resolved and conventional turbulence components, as were determined by different wire orientations, throughout the induction, the compression and the early stage of the expansion stroke. The anisotropy of turbulence spectral structure and its temporal evolution during the engine cycle were examined by evaluating the autospectral density functions and the time scales of each turbulence component in consecutive correlation crank-angle intervals.

The necessity for further reductions of in-cylinder pollutant formation and the opportunity to minimize engine development and testing times highlight the need of engine thermodynamic cycle simulation tools that are able to accurately predict the effects of fuel, design and operating variables on engine performance. In order to set up reliable codes for indicated cycle simulation in SI engines, an accurate prediction of heat release is required, which, in turn, involves the evaluation of in-cylinder turbulence generation and flame-turbulence interaction. This is generally pursued by the application of a combustion fractal model coupled with semi-empirical correlations of available geometrical and thermodynamical mass-averaged quantities. However, the currently available correlations generally show an unsatisfactory capability to predict the effects of flame-turbulence interaction on burning speed under the overall flame propagation interval.

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.

This paper presents calculations of external car aerodynamics by using the Partial-Averaged Navier-Stokes (PANS) variable resolution model in conjunction with the finite volume (FV) immersed-boundary method. The work presented here is the continuation of the study reported in Basara et al. [1]. In that work, it was shown that the same accuracy of predicted aerodynamic forces can be achieved by using Reynolds-Averaged Navier-Stokes (RANS) k-ζ-f model on both types of meshes, the standard body-fitted (BF), and on the immersed boundary (IB) mesh. Due to all well-known shortcomings of the steady state approach, in this work we deal with the Partially Averaged Navier-Stokes (PANS), which belongs to the hybrid RANS-LES (scale resolving / high fidelity) methods. This approach was developed to resolve a part of the turbulence spectrum adjusting seamlessly from RANS to DNS (Direct Numerical Simulation).