Search Results

This paper presents a numerical study on the impact of washcoat diffusion on the overall conversion performance of catalytic converters. A comprehensive transient 1D pore diffusion reaction model is embedded in state-of-the-art 1D and 3D catalytic converter models. The pore diffusion model is discussed with its model equations and the applied diffusive transport approaches are summarized. The diffusion reaction model is validated with the help of two available analytical solutions. The impact of basic washcoat characteristics such as pore diameters or thickness on overall conversion performance is investigated by selected 1D+1D calculations. This model is also used to highlight the impact of boundary layer transfer, pore diffusion and reaction on the overall converter conversion performance. The interaction of pore diffusion and flow non-uniformities is demonstrated by 3D+1D CFD simulations.

Wet clutches drag loss simulation is essentially linked to the clutch friction surface patterns in addition to the main geometry and conditions of the interface (relative speed, separation, inner and outer radius, viscosity and boundary pressures). The clutch patterns promote cooling flow and micro-hydrodynamic effects to aid clutch separation but greatly complicate the simulation of drag loss during separation. These drag losses are important in understanding the system losses as well as finding the most effective clutch cooling strategy. Typical clutch models either only consider simple patterns, such as radial grooves, or require significant simulation efforts to evaluate. Additionally, many simple models require calibration to measurement of the actual clutch they try to model before they provide a useful 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 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.

The paper focuses on an optimization methodology, which uses Design of Experiments (DoE) methods to define component parameters of parallel hybrid powertrains such as number of gears, transmission spread, gear ratios, progression factor, electric motor power, electric motor nominal speed, battery voltage and cell capacity. Target is to find the optimal configuration based on specific customer targets (e.g. fuel consumption, performance targets). In the method developed here, the hybrid drive train configuration and the combustion engine are considered as fixed components. The introduced methodology is able to reduce development time and to increase output quality of the early system definition phase. The output parameters are used as a first hint for subsequently performed detailed component development. The methodology integrates existing software tools like AVL CRUISE [5] and AVL CAMEO [1].

The advancing electrification of the powertrain is giving rise to new challenges in the field of acoustics. Film capacitors used in power electronics are a potential source of high-frequency interfering noise since they are exposed to voltage harmonics. These voltage harmonics are caused by semiconductor switching operations that are necessary to convert the DC voltage of the battery into three-phase alternating current for an electrical machine. In order to predict the acoustic characteristics of the DC-link capacitor at an early stage of development, a multiphysical chain of effects has to be addressed to consider electrical and mechanical influences. In this paper, a new method to evaluate the excitation amplitude of film capacitor windings is presented. The corresponding amplitudes are calculated via an analytical strain based on electromechanical couplings of the dielectric within film capacitors.

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

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.

This paper presents calculations of external car aerodynamics by using the finite volume (FV) immersed-boundary method. The FV numerical codes primarily employ Reynolds-Averaged Navier-Stokes (RANS) models. In recent years, and due to possibility to run very large computational meshes, these models are usually used in conjunction with the advanced near-wall models. Moreover, it has been often demonstrated that the accuracy of RANS near-wall models relies on the mesh quality near the wall so by the rule, larger number of wall body-fitted cell-layers are employed. An immersed boundary (IB) method becomes an attractive alternative to the ‘standard’ FV approaches especially when applied to low quality CAD data. In general, the IB method is less investigated and validated for the car aerodynamics, particularly in conjunction with advanced near-wall turbulence models and an adaptive mesh refinement (AMR).

The parameterization of the electrochemical pseudo-two-dimensional (P2D) model plays an important role as it determines the acceptance and application range of subsequent simulation studies. Electrochemical impedance spectroscopy (EIS) is commonly applied to characterize batteries and to obtain the exchange current density and the solid diffusion coefficient of a given electrode material. EIS measurements performed with frequencies ranging from 1 MHz down to 10 mHz typically do not cover clearly isolated solid state diffusion processes of lithium ions in positive or negative electrode materials. To extend the frequency range down to 10 μHz, the distribution function of relaxation times (DRT) is a promising analysis method. It can be applied to time-domain measurements where the battery is excited by a current pulse and relaxed for a certain period.

The paper addresses aspects of modeling cavitating flows within high pressure injection equipment while considering the effects of liquid compressibility. The presented numerical study, performed using the commercial CFD code AVL FIRE®, mimics common rail conditions, where the variation in liquid density as a function of pressure may be relevant owing to very high pressure injection scenarios. The flow through the injector has been calculated and the conditions at the outlet of the nozzle orifice have been applied as inlet condition for subsequent Euler-Lagrangian spray calculations to investigate the effects of liquid compressibility treatment on spray propagation. Flows of such nature are of interest within automotive and other internal combustion (IC) related industries to obtain good spray and emissions characteristics.

Fully flexible dual fuel (DF) internal combustion (IC) engines, that can burn diesel and gas simultaneously, have become established among heavy-duty engines as they contribute significantly to lower the environmental impact of the transport sector. In order to gain better understanding of the DF combustion process and establish an effective design methodology for DFIC engines, high fidelity computational fluid dynamics (CFD) simulation tools are needed. The DF strategy poses new challenges for numerical modelling of the combustion process since all combustion regimes have to be modelled simultaneously. Furthermore, DF engines exhibit higher cycle-to-cycle variations (CCV) compared to the pure diesel engines. This issue can be addressed by employing large eddy simulation coupled with appropriate DF detailed chemistry mechanism. However, such an approach is computationally too expensive for today’s industry-related engine calculations.

Intensive R&D is currently performed worldwide on hybrid and electric vehicles. For full electric vehicles the driving range is limited by the capacity of currently available batteries. If such a vehicle shall increase its driving range some range extending backup system should be available. Such a Range Extender is a small system of combustion engine and electric generator which produces the required electricity for charging the batteries in time. Since the acoustic response of an electric motor driving the vehicle and of a combustion engine as part of a Range Extender is very different by nature an extensive acoustic tuning of the Range Extender is necessary to meet the requirements of exterior vehicle noise and passenger comfort. This paper describes the NVH (noise, vibration & harshness) development work of a range extender within the AVL approach of an electrically driven passenger car with range extender.

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.

Commonly, the spray process in Direct Injection (DI) diesel engines is modeled with the Euler Lagrangian discrete droplet approach which has limited validity in the dense spray region, close to the injector nozzle hole exit. In the presented research, a new reactive spray modelling method has been developed and used within the 3D RANS CFD framework. The spray process was modelled with the Euler Eulerian multiphase approach, extended to the size-of-classes approach which ensures reliable interphase momentum transfer description. In this approach, both the gas and the discrete phase are considered as continuum, and divided into classes according to the ascending droplet diameter. The combustion process was modelled by taking into account chemical kinetics and by solving general gas phase reaction equations.

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.

This work elaborates the transferability of electrode diffusion coefficients gained from fitting procedures in frequency domain to an electrochemical battery model run in time domain. An electrochemical battery model of an NMC622 half-cell electrode is simulated with sinusoidal current excitations at different frequencies. The current and voltage signals are analyzed in frequency domain via Nyquist and Bode plots. The frequency domain analysis of time domain simulations is applied to assess the numerical convergence of the simulation and the sensitivity on particle diameter, electrode and electrolyte diffusion coefficients. The simulated frequency spectra are used to fit the electrode diffusion coefficient by means of different electrical equivalent circuit models and the electrochemical battery model itself. The fitted diffusion coefficients from the different electrical equivalent circuit models deviate by one order of magnitude from the a priori known reference data.

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.

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.