Transfer Matrix Computation for Wave Action Simulation in an Internal Combustion Engine

Abstract

A new technique for simulating engine pressure waves consisting of linking pressure response and mass flow rate excitation in the frequency domain has been presented. This is achieved on the so-called “dynamic flow bench”. With this new approach, precise, fast and robust results can be obtained while taking into account all the phenomena inherent to compressible unsteady flows. The method exhibited promising results when it was incorporated in a GT-Power/Simulink coupled simulation of a naturally aspirated engine. However, today’s downsized turbocharged engines come with more stringent simulation necessities, where discontinuities such as the charge air cooler (CAC) must be correctly modeled. Simulating such engines with the transfer function methodology is quite difficult because it requires mounting the entire intake line on the bench. Modeling wave action for these engines requires an understanding in the frequency domain of the flow’s characteristics through the different elements that make up the intake line. This leads us to study the acoustic transfer matrices. In order to split the intake line into separate elements, a straight duct of 185mm length is chosen as a first reference. It is mounted on the dynamic flow bench and pressure response is measured after an impulse mass flow excitation. Transfer functions of relative pressure and mass flow rate are then identified at given points upstream and downstream of this reference tube. These functions produce the desired transfer matrix poles. The resulting matrix is validated by inserting the tube in the intake lines of two four-cylinder engines which are modeled in GT-Power. Pressure and mass flow are registered at the measurement points of the tube from the simulation. The time series data upstream of the tube is treated in the frequency domain and the transfer matrix is used to calculate the corresponding downstream values. Measured values from the native simulation and those calculated using the transfer matrix propagation are then compared. Finally, the experimental technique for identifying transfer matrices of more complex elements using two versions of the previous tube is presented.