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 this paper, a recently improved Computational Fluid Dynamics (CFD) methodology for virtual prototyping of the heat treatment of cast aluminum parts, above most of cylinder heads of internal combustion engines (ICE), is presented. The comparison between measurement data and numerical results has been carried out to simulate the real time immersion quenching cooling process of realistic cylinder head structure using the commercial CFD code AVL FIRE®. The Eulerian multi-fluid modeling approach is used to handle the boiling flow and the heat transfer between the heated structure and the sub-cooled liquid. While for the fluid region governing equations are solved for each phase separately, only the energy equation is solved in the solid region. Heat transfer coefficients depend on the boiling regimes which are separated by the Leidenfrost temperature.

The present work introduces a fully integrated real-time (RT) capable engine and vehicle model. The gas path and drive line are described in the time domain of seconds whereas the reciprocating characteristics of an IC engine are reflected by a crank angle resolved cylinder model. The RT engine model is derived from a high fidelity 1D cycle simulation and gas exchange model to support an efficient and consistent transfer of model data like geometries, heat transfer or combustion. The workflow of model calibration and application is outlined and base ECU functionalities for boost pressure, EGR, smoke and idle speed control are applied for transient engine operation. Steady state results of the RT engine model are compared to experimental data and 1D high fidelity simulations for 19 different engine load points. In addition an NEDC (New European Drive Cycle) is simulated and results are evaluated with data from chassis dynamometer measurements.

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

Minimizing the lubricant volume in a transmission system reduces the churning losses and overall unit costs. However, lubricant volume reduction is also detrimental to the thermal stability of the system. Transmission overheating can result in significant issues in the region of loaded contacts, risking severe surface/sub-surface damage in bearings and gears, as well as reduction in the lubricant quality through advanced oxidation and shear degradation. The increasing trend of electrified transmission input speeds raises the importance of understanding the thermal limits of the system at the envelope of the performance to ensure quality and reliability can be maintained, as well as being a key factor in the development, effecting internal housing features for the promotion of lubrication. A nodal bearing thermal model will be shown which utilizes thermal resistances and smooth particle based CFD for determining bearing lubricant feed rates during operation.

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