Industrial Centrifugal Compressors Design Validation, using Rapid Method Evolutionary Algorithms to Predict Map Characteristics Evaluation.

Industrial Centrifugal Compressors Design Validation, using Rapid Method Evolutionary Algorithms to Predict Map Characteristics Evaluation.

J. C. Statharas*, N. W. Vlachakis

Department of Aircraft Engineering, Central Greece Technological University (TEI),

Psachna, Evias, 34400, GREECE.

Abstract

:- The purpose of this work is to model a radial compressor with a new method of predicting its map. This method calculates the air mass flow rate and the maximum pressure at each constant velocity curve of compressor operation. Two modeling approaches are proposed and applied to commercial compressors.

Keywords:- Centrifugal compressor, turbocharger, compressor map, compressor modeling, empirical model

1 INTRODUCTION

1. Air compressors

   Air compressors are mechanical devices used to compress air for many technical applications. Two of the main categories of air compressors are the axial and radial (or centrifugal) compressors. Axial and centrifugal compressors transfer energy from a rotating device to the air (or a gas), using shaft torque to impart the moment to the gas.

   Axial compressors are mainly used in large scale turbomachine applications such as in the power production industry and the aviation and marine propulsion.

2. Centrifugal compressors

   Radial (or centrifugal compressors) on the other hand have a wider use in small (also in micro) or medium scale applications such as in internal combustion supercharging or turbocharging and aviation propulsion (eg. Turboshafts).

   Really a large market exists for the internal combustion engine turbocharger, which is used in diesel and gasoline engines. The internal combustion engine turbocharger compresses the gas entering the manifold of the engine and is driven by a turbine running on the engine exhaust gases. In turbocharging a wide range of stable operation is necessary. This is an advantage of a centrifugal compressor with its generally wider operational area. Other applications of the compressor are for ventilation, process and chemical industries.

   In the majority of cases, the so-called radial flow impeller is used. The main reasons for establishing this technology are the simplicity and robustness of the construction, the low weight and size, and finally, with even a single stage, that is,

   only a compression step, a satisfactory high pressure ratio that meets the turbocharging needs is achieved.

   The performance of centrifugal compressors is due to the fast or very fast rotation of the impeller. The impeller rotation gives both kinetic energy and pressure to the air induced. The diffuser, as it guides the air peripherally to the spiral shell, slows it down by increasing its pressure. The diffuser may have fixed vanes or consist of single parallel walls. The spiral shell is a duct with continuously growing cross section. The duct surrounds the diffuser in the form of a spiral and collects the air from all the peripheral positions of the diffuser, leading it to a volute. The volute function is to duct flow away and to decelerate the flow further until the exit of the compressor. A detailed presentation of the main parts of the compressor is given in section 2.

3. Compressor map

   It is very useful for design purposes to know the behavior of the compressor during steady state operation. The steady state performance of the compressor is often specified in the form of a compressor map by curves that relate the rotational speed, the pressure rise across the compressor and the mass flow through the compressor. Examples of compressor maps are shown in section 3. Compressor maps usually are provided by the manufactures. Although compressor maps contain plenty of information for the compressor behaviour it is difficult for use in analytical studies. So modelling aproaches are necessary.

4. Compressor modeling

   Compressor modeling is essential for predicting engine performance correctly. To get correct engine performance one has to understand how the compressor works and the amount of air it provides to the engines cylinders for combustion. If compressors performance is incorrectly modeled, the amount of air is wrong and the engine output torque may be wrong due to incorrect air-to-fuel ratio and combustion efficiency. Moreover uncertainty issues arise in fuel economy and in the amount of emissions in the environment. Therefore, compressors performance has a crucial impact on engine combustion, fuel consumption,

   emmisions to the environment and output torque. It is imperative therefore to build a precise compressor model ensuring the accuracy of the whole engine model.

   Usually, the static compressor map provided by the manufacturer is not in a form that is useful in engine system modeling. In literature, two distinct ways of modeling compressors can be found:

   The first way is to store the compressor performance data into lookup tables and use mathematical algorithm to interpolate data.

   The second one is to construct fitting functions using the experimentaly derived data to recreate the compressors performance map. This method has been widely used in many works [1-11] giving acceptable prediction accuracy. Of course a vast range of modeling approaches and complexity levels exists, where some models are especially designed for a certain purpose, while others are meant to be more general. An extended review of model families is given in [12].

   This paper mainly focuses on the fitting function method (second method) for the reason that the compressor map developed by interpolation is not always appropriate for extrapolation purposes. For instance, the regions which are not covered by the performance map can not be estimated by interpolation. Moreover the majority of the existing models are valid for the compressor used for validation, not any compressor. The proposed model has a more universal use as can be applied to any compressor.

   In this study, the models of compressor flow and efficiency are fitted based on exponential function and polynomial function, respectively. Results are compared with the compressor experimental data to verify the performance of the model. In addition, the flow models derived by exponential fitting function and interpolation are compared leading to the conclusion that the fitting function methodology surpasses the interpolation by providing better extrapolation capability. We do not underestimate the excellent work done by other researchers but we feel proud for the universality of our model. Section 4 of this work presents the results of the proposed models.

1. THE GEOMETRY OF A RADIAL COMPRESSOR The radial compressor consists of the compression wheel

   (or impeller) and the housing (or shell). These two parts constitute the rotating and stationary part of the compressor and work together to compress the introduced air and to move it radially by means of its internal configuration. The impeller main characteristics are shown in the following Fig.2.1.

   Fig. 2.1. Impeller main characteristics. 1. Inlet diameter d1, 2. Tip height b2, 3.Main blade shape, 4. Intermediate small blade (Splitter), 5. Main blade,

   1. Rear wheel side, 7. Outlet diameter d2, 8. Outlet tip edge, 9. Impeller nose.

      More details of an impeller are shown in the cross section figured in Fig.2.2.

      Fig. 2.2 .Impeller cross section. 1.Central hole for the operating shaft, 2.Blade basis, 3.Blade slope for air movement, 4.Impeller rear side extension,

      5.Blade root

      There are several types of impeller designs. Each has its avantages and its corresponding use.

      Pure radial design (Fig.2.3) is not used today because it does not perform well despite the high speeds achieved. If it is used with a vaned-type diffuser, the efficiency is increased but it lacks in air mass flow rate and is not suitable for the automotive industry requiring a wide range of revolutions. It is used in diesel engines and generators operating at a limited number of revolutions.

      Fig. 2.3 . Pure radial impeller

      Impellers with only single normal blades (Fig.2.4) are rarely seen today and only in applications with lower rotation speeds. There are not suitable for applications with high performance and high compression. They lack in the air supply at high speeds despite the good performance and high pressure that achieves.

      Fig. 2.4. Impeller wheel with normal blades.

      Impellers with intermediate small blades (splitters) (Fig.2.5) in between the main blades, appear more frequently today because at high speeds the vacuum created helps to increase the air mass flow rate. That is, the small blade helps to manage and compress the air as it is accelerated according to the wheel rotation.

      Fig.2.6.Impeller with backward leaned vanes.

      There are vanes that are fully or partially backward leaned (Fig. 2.7). Fully backward leaned vanes are those that have a reversed curve across the diameter of the impeller. Partially backward leaned vanes were devised to reduce the weight of the wheel so as to be faster on accelerations. But they are not used because they are not so effective.

      Fig. 2.7 Fully or partially backward leaned vanes

      Another feature of the compressors is the trim. To understand the definition of trim we need to know the inflow diameter (inducer) and the outflow diameter (exducer) of the compressor. More specifically trim is defined as the area ratio:

      (inducer area )2

      (d )2

      Fig. 2.5 Impeller wheel with splitters

      trim=

      exducer area

      100= 1

      d2

      100

      (2.1)

      Backward leaned blades (or vanes) (Fig.2.6) are the most common in today's design. They allow the compressor to operate at a greater range of revolutions that is necessary in high performance vehicle engines since air diffusion begins before entering the diffuser.

      The compressor housing (or shell) creates the internal space in which the air moves (Fig.2.8).

      Fig. 2.8. Compressor housing.1. Spiral, 2. Compressor outlet, 3. Inlet surface, 4. Inner connection diameter

      At the end of the blade we see the diffuser (Fig.2.9). The shape of the diffuser helps to convert the air velocity to pressure.

      Fig. 2.9. Diffuser

      In some cases the diffuser has blades (see Fig. 2.10). This helps to increase the efficiency but also greatly reduces the air mass flow rate operating range. So it will only be used in low- speed engines.

      Fig. 2.10. Diffuser with blades

      Another feature of the turbo compressors is the A/R ratio. It is the ratio of the surface area of the cross section of the spiral that is called A, to the distance defined by the center of the impeller to the radius of area A and is called R (Fig.2.11).

      Fig. 2.11. Definition of A/R ratio.

      This parameter has a strong effect on the performance and operating characteristics of the turbocharger. The A/R ratio of the compressor does not play an important role because its change does not affect valuably the whole turbocharger system. From the turbine side, however, the A/R is used to adjust flow capability and matching.

      On the intake side there is an opening perimetrically to the shell which allows recirculation of the air (Fig.12). This has the effect of reducing the surge problem especially at low revolutions since there is always air available for compressor operation and flow stabilization.

      M = U 2

      RT (3.5)

      where cp is the air specific heat at constant pressure, T is the ambient air temperature, is the compression ratio, is the ratio of air specific heats and R is the air specific gas constant.

      1. Compressor models

         1. According to [13], the basic parameter is expressed as a function of and M in the following way:

            1+ 2

            i

            i

            = , = k + k

            M , i= 1,2 ,3

            3

            i1 i2

            (3.6)

            Fig. 2.12. Opening for recirculation

2. COMPRESSOR MODELING

There are various methods for modeling the compressor behavior which are based on the parameterization of the dimensionless air mass flow rate using several parameters. From the different parameterizations developed, a simple method has been developed in [13] and shows that a simple model gives good accuracy among the most important areas of compressor maps. The advantages of the simple model are:

• A good description of the air mass flow rate in the low-speed area, where the machines operate most of the time.

• Simple adjustment

The models are based on the parameterization of the flow factor :

where Ki coefficients to be determined.

1. According to [14], and are expressed as a quarter of an ellipsis (Fig. 3.1), which results in the flow factor being expressed only as a function of the basic parameter .

   Fig. 3.1. Variation of versus according to [14]

2. According to [15] the flow factor can also be expressed only as a function of :

   = + (1 e 2 + )

   = + (1 e 2 + )

   3

   3

   = m c

   d 2 U 1

   4

   4

   2 2

   (3.1)

   where c is the air mass flow rate, is the air density and

   U2 is the tip speed needed:

   U = d2

   2 60 (3.2)

   where N is the number of revolutions per minute (RPM). Assuming that the flow factor can be parameterized as ,

   then it is possible to invert and express the air mass flow rate as a function of :

   m = ( d 2 U )

   c

   c

   (3.7)

1. Finding the coefficients

   Irrespectively of the method used, the coefficients, which apply only to a particular compressor, must be determined. So each time we want to model a different compressor, it is necessary to repeat the procedure of determining the coefficients again and again.

   In all types, the K1, K2 and K3 parameters (if used) are determined using the least squares method in the Matlab Curve Fitting Toolbox.

   4 2 2

   (3.3)

   Ki = ki1+ ki2 M ,i= 1,2,3 (3.8)

   K1=k11+k12M

   )

   )

   We still have to model the flow factor. There are three models for this. In all the models proposed, the basic parameter and Mach-number M are used

   1

   = cP Ta 0.5U 2

   K2=k21+k22M K3=k31+k32M

   We selected the Mitsubishi TD04-09B compressor for

   modeling its map. The main characteristics of this compressor

   map is shown in Fig.3.2. Below we calculate the coefficients

   (

   (

   1

   2

   (3.4)

   are shown in Table 3.1 and the coressponding compressor

   from the Mitsubishi TD04-09B compressor map (Fig.3.2). First, from the map we find the experimental points as (mc)

   at different rotational speeds N. Then equations (3.1), (3.2), (3.4), (3.5) are used to calculate (). Table 3.2 shows an example of the results for rotational speed N=140000 RPM.

   Table 3.1 Basic characteristics of Mitsubishi TD04 -09B compressor

   Turbo model:	Compressor Trim:	d1 in/(mm)	d2 in/(mm)	Compressor Diameter Average:
   Mitsubishi TD04-09B	50	1.365 (34.671)	1.930 (49.022)	1.6475

   rotational speed N=140000 RPM
   Uc=359,17 m/sec		<>M=1,0345
   	mc
   1,986	0,037	0,044396	1,007288
   1,982	0,047	0,056395	1,004029
   1,973	0,053	0,063595	0,996681
   1,963	0,060	0,071994	0,988487
   1,949	0,067	0,080393	0,976966
   1,94	0,072	0,086393	0,969528
   1,931	0,076	0,091193	0,962066
   1,898	0,083	0,099592	0,934489
   1,87	0,087	0,104391	0,910821
   1,837	0,091	0,109191	0,882599
   1,794	0,096	0,115191	0,845277
   1,752	0,101	0,12119	0,8082
   1,695	0,107	0,128389	0,756852
   1,629	0,114	0,136789	0,695832
   1,61	0,116	0,139189	0,677939

   rotational speed N=140000 RPM
   Uc=359,17 m/sec		M=1,0345
   	mc
   1,986	0,037	0,044396	1,007288
   1,982	0,047	0,056395	1,004029
   1,973	0,053	0,063595	0,996681
   1,963	0,060	0,071994	0,988487
   1,949	0,067	0,080393	0,976966
   1,94	0,072	0,086393	0,969528
   1,931	0,076	0,091193	0,962066
   1,898	0,083	0,099592	0,934489
   1,87	0,087	0,104391	0,910821
   1,837	0,091	0,109191	0,882599
   1,794	0,096	0,115191	0,845277
   1,752	0,101	0,12119	0,8082
   1,695	0,107	0,128389	0,756852
   1,629	0,114	0,136789	0,695832
   1,61	0,116	0,139189	0,677939

   Fig. 3.2. Mitsubishi TD04-09B compressor map Table 3.2 (example N=140kRPM)

   From the calculated () using the Matlab Curve Fitting Toolbox (Table 3.3) we find the coefficients K1, K2 and K3 for different rotational speeds N (Table 3.4). For example, the coefficients for the method in [13] are shown in Fig. 3.3. We see from Figure 3.3 and in Table 3.3 the function f(x)=(a+b*x)/(c-x) corresponds to Eq.(3.6) and a, b, c are the coefficients K1, K2 and K3.

   Fig. 3.3. K1, K2 and K3 coefficients for the method in [13].

   Table 3.3 matlab curve fitting equation output

   General model: f(x) = (a+b*x)/(c-x) Coefficients (with 95% confidence bounds):

   a =	0.1917 (0.1715, 0.2119)
   b =	-1.182 (-1.253, -1.11)
   c =	0.1801 (0.1651, 0.1951)

   Table 3.4

   Jensen [13]
   N	M	K1	K2	K3
   90000	0,66504	0,2541	-1,483	0,2176
   120000	0,88672	0,2194	-1,295	0,1978
   140000	1,0345	0,1916	-1,181	0,18
   160000	1,18229	0,3534	-2.458	0,1493
   180000	1,33008	0,1867	-0,9549	0,1991

   The next step is to find the coefficients kii using the equation (3.8) that depend on the number M. The least squares method of Matlab is used. The final coefficients for all methods are shown in Table 3.5.

   Table 3.5 The coefficients of equations

   Method	k11	k12	k21	k22	k31	k32
   [13]	0,3215	-0,10135	-2,0111	0,794088	0,2361	-0,02782
   (Ellipsis) [14]	59,78	-35,261	0,587	0,302388	–	–
   [15]	786,45 3	1319,45	0,970275	2,528758	1,757268	4,14

   In [15], the ki coefficients must be carefully calculated especially at high speeds (N>140000 RPM) because in exponential terms when rounding can easily result in incorrect results. This method is proposed in large diameter turbocompressors. d2> 0.5 m.

   To calculate air mass flow rate and compression ratio according to each method, the following equations are used:

   m = d 2 U

   c 4 2

   [
   [

   2

   cp Ta

   2

   (3.9)

   ]

   =

   (3.10)

   0,5U2 + 1 1

   Figure 3.4 shows the results of the three methods applied to the TD04-09B compressor.

   Fig. 3.4. Intercomparison of the modeling methods

   The efficiency nc can also be expressed as a function of and in the following way:

   n = a 2+ + , = i1+ ai2 M

   Fig. 3.5. TD04-09B Compressor efficiency

1. AD-HOC MODELING OF THE RADIAL COMPRESSOR MAP.

It would be very useful if we could, by using some equations, to determine the characteristic constant velocity curves of the compressor map. It would also be even more useful if this resulted from simple geometric characteristics of the compressor and applied to a wide range of compressor sizes. In the preceded methods the modeling of compressor maps and the calculation of respective coefficients are made for each compressor. In our own methodology, the calculated coefficients apply to all compressors.

1. Statharas-Vlachakis calculational procedure

   The calculation of a minimum and a maximum value of the air mass flow rate and the calculation of the maximum compression ratio of the compressor will be made. The calculation is based on empirical equations and geometric characteristics of the compressor.

   The minimum value of air mass flow rate is:

   c 1

   (3.11)

   2 3 i

   ai3 M

   mmin

   = 0,08 cr1

   d 2

   1 z0,4

   4

   (4.1)

   Figure 3.5 shows the results of the efficiency nc for the TD04-09B compressor.

   The radial inlet velocity is:

   =

   =

   c

   c

   cu1

   r1

   tan a1 (4.2)

   The peripherl component of the absolute input velocity is:

   cu1= u1 (4.3)

   The peripheral velocity is:

   1

   1

   u = d1 N

   60

   The load factor is :

   = 1

   2 + 1

   The blade blockage factor of diameters ratio is :

   a= (0,53 1 0,208)( 2)

   tan

   tan a2

   (4.5)

   d d

   d 2 d 0

   0,05

   (4.8)

   The Weisner slip factor:

   cos

   where d0 =0,032 m.

   A good estimation of the radial exit velocity is suggested:

   = 1 1 c = 0,36 u

   z (4.6)

   where z is the number of blades 1 inlet flow angle

   2 outlet flow angle 1 inlet blade angle 2 outlet blade angle

   r2 2 (4.9)

   The peripheral outlet velocity is:

   2

   2

   u = d 2 N

   60 (4.10)

   The maximum pressure ratio of the compressor is :

   1

   max

   = c(1+ 0,42( 1) 2)

   (4.11)

   Where c is the correction coefficient for the outlet diameter d2:

   c= 4,17 d2+ 0,866 (4.12)

   and the tip Mach number:

   M = u2

   R T (4.13)

2. Determining the compressor map

   Fig. 4.1. Inlet-outlet compressor velocities

   Depending on the mass flow rate, the curves of the compressor map will be derived according to the following empirical equation:

   ( m)= 1,03 erfc[ 0,2 ( m i )]

   erfc max mmax

   2e

   (4.14)

   The above equation uses the complementary error function [16].

   Similarly, in relation to the mass flow rate, the curves of the compressors map can also be obtained according to the following empirical equation

   mi 2e

   m

   m

   poly ( m)= 1,02 max [1 0,3 ( )]

   max (4.15)

   The ratio of the masses i/max takes intermediate values between i=min and i=max. The coefficient of compressor morphology in both equations is given by the equation:

   e=27,78d2+1,111 (4.16)

3. Application-Examples

   	Compressor	Inlet diameter d1 (mm)	Outlet diameter d2 (mm)	trim
   1	Garrett gt06_50	22,63	32	50
   2	Garrett gt12_50	29	41	50
   3	Garrett gt15_56	32,9	43,9	56
   4	Garrett gt17_49	34,3	49	49
   5	Mitsubishi td04-09b	34,7	49	50
   6	Garrett gt25_550	48,4	60	65
   7	Garrrett gtx3077r	58	76	58
   8	Garrett gt4202r	75,2	102	54
   9	Holset hx82	98	141	48

   	Compressor	Inlet diameter d1 (mm)	Outlet diameter d2 (mm)	trim
   1	Garrett gt06_50	22,63	32	50
   2	Garrett gt12_50	29	41	50
   3	Garrett gt15_56	32,9	43,9	56
   4	Garrett gt17_49	34,3	49	49
   5	Mitsubishi td04-09b	34,7	49	50
   6	Garrett gt25_550	48,4	60	65
   7	Garrrett gtx3077r	58	76	58
   8	Garrett gt4202r	75,2	102	54
   9	Holset hx82	98	141	48

   To confirm this methodology we will apply it to various commercial compressors that cover the entire range of compressors that can be used in piston combustion engines. The compressors used have the following characteristics as shown in Table 4.1.

   Fig. 4.2. Inlet-outlet compressor angles

   The maximum mass flow rate is given by the formula :

   d 2

   Table 4.1 Compressors used

   mmax

   = ac

   2 z0,4

   r2 4

   (4.7)

   The figures 4.3-4.10 show the compressor maps predicted by using the new method. In the figures the two procedures are shown as erfc and poly. As we can see, the approach with the new prediction method is very good, since the constant velocity curves created are close to the actual operating points of the compressors given by the manufacturers. Particularly in the garrett gt06-50, gt17-49 and gt12-50 compressors the approximation of the air mass flow rate constant velocity curves but also the pressure ratio of the compressor is very good.

   3 Garett gt06-50

   2,8

   2,6

   2,4

   280000 rpm

   280000 rpm

   2,2

   260000 rpm

   260000 rpm

   2

   220000 rpm

   220000 rpm

   240000 rpm

   240000 rpm

   1,8

   200000 rpm

   200000 rpm

   1,6

   1,4

   180000 rpm

   180000 rpm

   1,2

   1

   Fig.4.5. Garrett gt17 49 compressor map

   0 0,01 0,02 0,03 0,04 0,05 0,06 0,07

   m (kg/s)

   Fig.4.3. Garrett gt06-50 compressor map

   Fig.4.4. Garrett gt12-50 compressor map

   Fig.4.6.Mitsubishi td04-09b compressor map

   Fig.4.7. Garrett gt25_550 compressor map

   Fig.4.8. Garrett gtx3077r compressor map

   Fig.4.9. Garrett gtx4202r compressor map

   Fig.4.10. Holset hx82 compressor map

   If we compare the new methodology with the old ones (Fig. 4.11) we see that we have a very good approach with respect to the older modeling approaches. Let us note that the older models do not calculate the air mass flow rate not the pressure compression ratio.

   Fig. 4.11. Intercomparison of compressor map prediction methodologies

4. Particularities of compressors

   When choosing the right compressor we have to take account of the trim . That is, the larger the ratio of the diameters square of the impeller, the greater the air mass flow rate that compressor can support.

   trim=

   trim=

   d

   d

   1

   1

   2

   100

   2

   100

   (d2 )

   Compressors with large inlet to small outlet diameter will support large amounts of air at low pressure as opposed to a small inlet to large outlet diameter compressor that supports less mass flow rate at relatively higher pressures.

   validated by comparing with the compressor performance map provided by the manufacturer. Results show that the performance predicted by the models is closely comparable with the real performance map. Moreover the error function based flow rate model shows credible extrapolation capability which the interpolation methodology does not have. The compressor model can be applied to mean value model of turbocharged diesel engine for non-linear control and state estimation.

   REFERENCES

   Fig. 4.12. Compressor maps having different trim

   In the two compressor maps (Figure 4.12) it is evident that although having the same outlet diameter d2 the larger trim '60' compressor has higher mass flow rates at lower pressures. For example, the first compressor for P2/P1 = 2.2 has an air mass flow rate of up to 44lb/min while in the second compressor for the same pressure ratio the air mass flow rate is up to 50lb / min.

1. CONCLUSIONS

The work presented in this paper address the issue of modeling the compressor performance map, proposing an error function based model for flow rate. The models are

1. L. Eriksson, L. Nielsen, J. Brugård, J. Bergström, F. Pettersson, and Per Andersson, Modeling of a turbocharged SI engine, Annual Reviews in Control, 26(1):129137, 2002.

2. .T. Gravdahl, O. Egeland, and S.O. Vatland, Active surge control of centrifugal compressors using drive torque, vol.2 of Proc. IEEE Conf. on Decision and Control, 2001.

3. L. Guzzella and A. Amstutz, Control of diesel engines, Control Systems Magazine, IEEE,18(5):5371, Oct 1998.

4. L. Guzzella and A. Amstutz, Control of diesel engines, Control Systems Magazine, IEEE,18(5):5371, Oct 1998.

5. F.K. Moore and E.M. Greitzer, A theory of post-stall transients in axial compression systems: Part i – development of equations, Trans. of ASME, J. of Engineering for Gas Turbines and Power, (108):6876, 1986. .

6. P. Moraal and I. Kolmanovsky. Turbocharger modeling for automotive control applications. In SAE World Congr., Techn. Paper 1999-01- 0908, March 1999.

7. M. Müller, E. Hendricks, and Spencer C. Sorenson, Mean value modelling of turbocharged spark ignition engines, In SAE World Congr., Techn. Paper 980784, Detroit, USA, 1998.

8. A. Plianos and R. Stobart, Modeling and control of diesel engines equipped with a two-stage turbo-system, In SAE World Congr., Techn. Paper 2008-01-1018, 2008.

9. A. Thomasson and L. Eriksson, Modeling and control of co-surge in bi-turbo engines, In Proc. of IFAC World Congr., pages 13010 13015, 2011.

10. J. Wahlström and L. Eriksson, Modelling diesel engines with a variable-geometry turbocharger and exhaust gas recirculation by optimization of model parameters for capturing non-linear system dynamics, Proc. of the IMechE, Part D: J. of Automobile Engineering, 225(7):960986, 2011.

11. . E. Ffowcs Williams and X. Y. Huang, Active stabilization of compressor surge, Journal of Fluid Mechanics, 204:245262, 6 1989.

12. O. Leufvén, Modeling for control of centrifugal compressors, PhD thesis 1516, Linköping University, 2013.

13. J. P. Jensen, A.F. Kristensen, S.C. Sorenson, N. Houbak, and E. Hendricks, Mean value modeling of a small turbocharged diesel engine, In SAE World Congr., Techn. Paper 910070, Feb 1991.

14. P. Andersson, Air Charge Estimation in Turbocharged Spark Ignition Engines, PhD thesis 989, Linköping University, 2005.

15. H. Tu, H. Chen, Modeling of a Compressor's Performance Map by Fitting Function Methodology, Advanced Materials Research Vols. 779-780 (2013) pp 1194-1198.

16. http://mathworld.wolfram.com/Erfc.html