Applicable and Comparative Research of Compressor Mass Flow Rate and Isentropic Efficiency Empirical Models to Marine Large-Scale Compressor

: A compressor is an indispensable component of marine large two-stroke diesel engines. For this type of engine, the compressor mass ﬂow rate and isentropic e ﬃ ciency empirical models are preferred for both the working cycle dynamic simulation research and the design and testing of control and diagnostics algorithms due to their compact and simple structures, and satisfactory prediction accuracy. Due to absence of comprehensive applicable and comparative research on compressor mass ﬂow rate and isentropic e ﬃ ciency empirical models for large-scale marine compressors in the literature, two marine compressors with di ﬀ erent size, ﬂow rate range, and speed range were selected as research objects in this paper


Introduction
Due to the vital role of the turbocharger in improving power density, lowering fuel consumption, and reducing emissions for diesel engines, researchers have conducted extensive theoretical and experimental studies to exploit its capability further.In the field of automotive engines, several new types of turbocharging systems have emerged, such as electrically assisted turbocharging systems (E-Turbo), variable geometry turbochargers (VGTs), and variable geometry compressors (VGCs).It was revealed by many pieces of research that, compared with traditional turbocharging technology, these advanced technologies can not only effectively improve turbocharger performance at low loading conditions but also substantially extend the operating range [1,2].Nowadays, large marine two-stroke diesel engines are generally equipped with a turbocharging system but they still lag behind automotive engines to some extent.Nevertheless, the turbocharging system substantially contributes to an improvement in power, economic and emission performance for marine diesel engines, due to the fact that they can effectively increase the density of scavenging air, which not only contributes to the complete burning of injected fuel but also allows more fuel to be injected into the combustion chamber in a single working cycle [3].
Due to the important role of turbochargers for both automotive and marine diesel engines, it is indispensable to develop corresponding mathematical models with satisfactory predictive accuracy for not only the working cycle dynamic simulation research of turbocharged diesel engines but also the design and testing of control and fault diagnosis algorithms.Therefore, researchers have conducted extensive studies on modeling the methodologies of turbochargers.The modeling methodologies in the existing literature exhibit different levels of complexity, accuracy, and computational requirements, which can be roughly classified into the look-up table method, the curve fitting method, the neural network method, and the theoretical modeling method [4].It should be noted that most research focuses on small-scale automotive turbochargers, and relevant research on large-scale marine turbochargers are scarce.The reason is that the large size and high flow range of marine turbochargers requires special and expensive experimental equipment, making it very costly to carry out relevant experimental studies; in addition, the manufacturers of marine turbochargers rarely provide turbocharger performance maps, which makes it very difficult for model development and tuning.Nevertheless, some relevant studies on large-scale marine turbochargers, especially for compressors, can be found in the literature.In [5], the compressor mass flow rate was expressed as a polynomial function of rotational speed and pressure ratio.Karlson [6] utilized the exponential function to fit the relationship between the dimensionless flow coefficient and the head coefficient for a compressor, and the isentropic efficiency was expressed as a polynomial function of the mass flow rate and rotational speed.Hansen et al. [7] used a two-dimensional fourth-order polynomial to model the mass flow rate and isentropic efficiency of a compressor.Sun et al. [8] developed a compressor mathematical model that can predict the surging characteristics and proposed a method to determine the characteristics curve of a compressor under unstable working conditions.Tu and Chen [9] attempted to model compressor mass flow rates based on the Kriging algorithm.Wang et al. [10] developed a theoretical model for a marine centrifugal compressor based on thermodynamic laws.In [11], a marine turbocharger model based on the average line method was proposed, the highlight of which was that the model could be tuned with a relatively small amount of measured data.
By analyzing the above-mentioned researches, it can be found that most modeling methodologies of marine compressors for both the mass flow rate and isentropic efficiency fall into the curve fitting method category.This is due to the fact that relative to their lookup table counterparts, stronger physical fundamentals can be embedded into the models developed with this method.In addition, they are much more computationally efficient than those developed with theoretical modeling methods.As a result, the curve fitting method can yield time-efficient mathematical models with compact and simple structures, which can not only satisfy the requirements of control-and/or diagnostic-oriented dynamic simulations, but also possibly apply to design, analysis, and performance simulation [12].It should be noted that the models developed with the curve fitting method are often referred to as empirical, semi-empirical, or mean value models in the literature.
By extensive literature survey, we have not found any paper providing a comprehensive applicable and comparative study of compressor mass flow rates and empirical isentropic efficiency models to a large-scale marine compressor.Two similar papers can be found in the literature.However, their research object was the centrifugal compressor for jet engines and refrigeration systems, respectively, and in addition, they only focused on the design operating area [12,13].To fill in this gap, two large-scale marine compressors with different sizes, flow rates, and speed ranges were selected as the research object in this paper, and a comprehensive applicable and comparative analysis was conducted to compare and analyze the predictive ability of several classical models, and several recently proposed compressor mass flow rate and empirical isentropic efficiency models.The comparative range includes both the design and off-design operating area.
The rest of the paper is organized as follows: The characteristic of compressor performance map and the definitions of several relevant dimensionless numbers is first introduced in Section 2.Then, various compressor mass flow rate and isentropic efficiency models are described briefly in Section 3.After that, the main technical parameters of two marine compressors, comparative methodology, and error evaluation criteria are presented in Section 4. Results of the comparisons and analysis results of mass flow rate and isentropic efficiency models are given in Sections 5 and 6, respectively.Section 7 summarizes the comparative and analysis results and discusses some limitations of existing empirical compressor models and feasible solutions.The last section of this paper is dedicated to conclusions and several guidelines are summarized, which can be followed when selecting or developing models for large-scale marine compressors.

Compressor Performance Map
The working characteristic of a compressor is usually represented by its performance map, as shown in Figure 1, with mass (or volume) flow rates and pressure ratios as the horizontal and vertical coordinates, respectively.In the performance map, measured data points with the same rotational speed are connected, forming the iso-speed lines, and those with the same isentropic efficiency are connected forming the iso-efficiency lines.Essentially, the compressor mass flow rate and isentropic efficiency empirical models can be considered as a mathematical description of the performance map [12].To cover different operating environments, corrected quantities including the corrected rotational speed N corr and corrected mass flow rate .m corr are usually utilized in the performance map with their definitions as follows
where T ref and p ref are the reference temperature and pressure, respectively, T in and p in are the temperature and pressure at the compressor inlet, respectively.The derivation of Equations ( 1) and ( 2) is based on the dimensional analysis of the compression process in a compressor [14].For a compressor, its whole operating area is normally divided into three zones, as shown in Figure 1, i.e., design operating zone, surging zone, and choking zone.In the design operating zone, the compressor is able to work stably and achieve high efficiency.In the surging zone, the momentum of airflow cannot stably overcome the adverse pressure gradient across the compressor, and reverse flow occurs periodically, resulting in a mild flutter to a wildly fluctuating boost, and the resultant noise known as "surging" [12,15].To prevent premature turbo failure, the compressor should avoid working in this zone continuously.The curve connecting the surging point on each iso-speed line is referred to as a surging line.When entering the choking zone, the flow reaches a state of sonic velocity, and the mass flow rate will no longer increase despite the decrease in pressure ratio, which is characterized by a steeply descending iso-speed line at the lower right side of the performance map, as shown in Figure 1.Similar to the surging line, the curve connecting the choking point on each iso-speed line is referred to as the choking line.
As shown in Figure 1, most of the compressor performance maps provided by turbocharger manufacturers only contain discrete measured data points in the design operating zone; however, for marine compressors, its actual rotational speed may be lower than the lowest rotational speed showed in the performance map, or higher than the highest one, and its pressure ratio may be lower than one under certain working conditions, such as main engine low-load running, activation of the auxiliary blower, and ship maneuvering [16].Therefore, it is required that the developed compressor model is capable of extrapolating to these off-design operating zones accurately and reasonably.In addition, since the performance map only provides several discrete iso-speed and iso-efficiency lines, the developed model is also required to interpolate within the unknown zones between these discrete lines with sufficient accuracy.Based on the above analysis, the range of applicable and comparative analysis conducted in this paper includes both the predictive accuracy in the design operating area and the extrapolation ability in the off-design operating areas.The off-design operating areas include the area with rotational speed lower than the lowest one available in the performance map, the area with rotational speed higher than the highest one, and the area to the right of the curve that connects the measured operating point of each iso-speed line with highest flow, which are termed LS, HS and LP area in this paper for convenience of expression, respectively It should be noted that by observing and comparing the performance map of small-scale automotive and large-scale marine compressors, it is known that the latter has a narrower shape, and its choking and surging lines are much closer than the former [16].Due to these differences, some compressor mass flow rates and isentropic efficiency empirical models that are well suited to small-scale automotive compressors may not be suitable for large-scale marine compressors.

Dimensionless Numbers
When developing compressor mass flow rate and isentropic efficiency empirical models, several dimensionless numbers are often utilized mainly including head coefficient Ψ, flow coefficient φ, and Mach number M, and their definitions are as follows [17] where c p , k, R and ρ are the constant pressure specific heat, specific heat ratio, gas constant, and density of air, respectively;d and U is the impeller outlet diameter and the blade tip speed, respectively.By observing Equations ( 3)-( 6), it is found that most of the physical quantities involved in the compression process of a compressor are covered by these dimensionless numbers, which greatly facilitates the model development process.Essentially, the theoretical foundation of many compressor empirical models existing in the literature is the identification of the functional relationship among these dimensionless numbers [18][19][20].

Compressor Mathematical Models
Researchers around the world have proposed various types of compressor modeling methodologies with different levels of complexity and accuracy, which can be roughly classified into four categories as following:

•
Look-up table method.This method utilizes the measured data points available in the compressor performance map to establish two-dimensional interpolation mesh, as shown in Figure 2, with which the mass flow rate and isentropic efficiency can be calculated at given rotational speed and pressure ratio conditions by using an appropriate interpolation algorithm.Since it is easy to realize, this method has been adopted in many numerical simulation studies [21,22].Nevertheless, it should be noted that because the standard interpolation algorithm is not continuously differentiable, simulation results with obvious discontinuous points may be obtained.
In addition, the extrapolation results obtained with this method are not reliable due to lacking strong physical fundamentals [23].

•
Neural network method.Due to the strong nonlinear mapping and generalization ability of the neural network, many researchers have attempted to apply this method to the prediction of the working characteristics of a compressor, and satisfactory results are attained [24,25].It was pointed out by Gholamrezaei and Ghorbanian [25] that BP (Back-Propagation) neural network with the Levenberg-Marquardt algorithm can achieve satisfactory prediction accuracy.The difficulty of this method lies in how to determine the structure of the neural network appropriately.The under-fitting phenomenon may occur if the structure is too simple and the quantity of nodes is too less; on the other hand, although the prediction capability of the neural network may improve to some extent by increasing its structural complexity and the number of nodes, its generalization ability may degrade, i.e., over-fitting phenomenon may occur; • Curve fitting method (empirical, semi-empirical, or mean value models).This method is widely used for predicting the working characteristics of compressors, mainly owing to its simpler and much more compact model structure than their theoretical counterparts, as well as their satisfactory predictive and extrapolation ability [3,6,[16][17][18][19][20][26][27][28][29].The curve fitting method can be further divided into two sub-categories: Some studies directly fit the relationship between the rotational speed, mass flow rate, pressure ratio, and isentropic efficiency with appropriate functions; while some researchers utilize suitable functions to fit the relationship between different dimensionless numbers derived from measured quantities.It should be noted that the difficulty of this method lies in how to select appropriate fitting functions and the determination of initial guesses of model parameters.

•
Theoretical method.The parameterization of compressor models developed with the above mentioned three methods strongly depends on the performance map; thus, their extrapolation ability is limited to some extent.Aiming at this problem, some researchers have attempted to develop theoretical models to predict the working characteristics of compressors without the help of performance maps.These theoretical models can be divided into 3-D, 2-D, and 1-D models in order of decreasing complexity [2,11,[30][31][32].Due to the large computational requirements, the 3-D and 2-D theoretical models are not suitable for control-and diagnostic-oriented dynamic simulations.Moreover, detailed compressor geometric dimension information is always needed.Compared with 3-D and 2-D models, the structure of 1-D theoretical models is relatively simpler and more compact, and it only requires basic geometric information.For developing 1-D theoretical models, it was first assumed that the compression process of a compressor is the isentropic compression process of compressible fluid, based on which the working characteristics of a compressor can be preliminarily estimated with the Euler equation, continuity equation, and energy conservation equation.Then, various types of losses (friction loss, incidence loss, shock loss, slip, wake loss, etc.) occurring during the actual compression process are modeled with the help of well-known empirical and semi-empirical formulas.Finally, the working characteristics of a compressor can be predicted by combing the results of the above two steps.It should be noted that although the theoretical foundation of 1-D theoretical models is the physical description of the actual compression process, many model parameters still need to be parameterized with the measured data provided in the performance map; in addition, some of the geometric information required by the 1-D theoretical model is still difficult to obtain, especially for large-scale marine compressors.
(a) (b) Because the ultimate purpose of this paper is to find a compressor model that is most suitable for the modeling of marine two-stroke diesel engine used for the marine engineering simulator, which has a relatively high requirement on both simulation speed and predictive accuracy, so only the empirical models developed with the curve fitting method are compared and analyzed in this paper owing to their simplicity, superior time-efficiency, and satisfactory predictive ability.

Compressor Mass Flow Rate Models
This section will briefly introduce several classical models and recently proposed compressor mass flow rate empirical models.In terms of the selection of the input and output parameters, these models can be divided into two categories with the model structure, as shown in Equations ( 7) and (8), respectively.Moraal and Kolmanovsky [23] pointed out in their research that the models manifested different strengths and weaknesses in terms of their sensitivity to input errors and model rigidity.Nevertheless, most researchers still select Equation (7) to develop compressor mass flow rate sub-models in the working cycle dynamic simulation model of turbocharged diesel engines, i.e., pressure ratio and rotational speed as the input, and mass flow rate as the output [33,34].
It should be noted that there are several models that are capable of working in both directions, such as the Leufvén and Llamas Ellipse Model introduced in Section 3.1.7[3,16,28,29].

The Guan Cong Model
Guan et al. [34] proposed a zonal modeling approach for marine compressors, which is referred to as the Guan Cong model in this paper.The basic idea of the Guan Cong model is to divide the whole operating area into several zones according to the iso-speed lines available in the performance map, and then utilize Equation (9) to fit the relationship between φ, Ψ, and M in each zone: The six model parameters k 1 to k 6 in Equation ( 9) can be determined through the least square method with the measured data points at the upper and lower boundaries of the iso-speed line of the current zone.Due to lacking lower boundary iso-speed line for the first zone, the model parameters corresponding to the lowest iso-speed line are used when extrapolating into this zone, similarly, for the last zone, the model parameters corresponding to the highest iso-speed line are utilized.Undoubtedly, the extrapolation strategy adopted by the Guan Cong model will inevitably influence its extrapolation accuracy, especially for the HS area, where the compressibility of air changes with compressor rotational speed significantly.

The Karlson Exponential Model
When developing the mass flow rate model for the TCA55 marine compressor, Karlson [6] proposed two types of exponential models with a different structure to fit the relationship between φ and Ψ, which are referred to as the Karlson-I and Karlson-II models, respectively.Equation ( 10) presents the Karlson-I model [6] φ where a, b, and c are all functions of the Mach number.Karlson-II model is defined as follows [6] where, α, β, and Ψ th are also expressed as polynomial functions of Mach number, similar to the Karlson-I model.

The Malkhede Model
When developing the MVEM (mean value engine model) for a marine turbocharged diesel engine, Malkhede et al. [5] directly adopted the polynomial function as shown in Equation (12) to calculate the mass flow rate of the compressor at given rotational speed and pressure ratio conditions:

The Kolmanovsky Model
Similar to the Karlson model, the Kolmanovsky model utilizes the exponential function to describe the relationship between the flow coefficient and the pressure ratio [23].
Polynomial functions, as shown in Equations ( 13) and ( 14) are used to describe the changing trend of flow coefficient and pressure ratio with compressor rotational speed at the surging line: For the area to the right side of the curve represented by Equations ( 13) and ( 14), the relationship between φ and Π can be described with the exponential function as follows: For the left side, the relationship can be described with the following linear function: where A, B, and C are functions of U.
In turn, the pressure ratio can be derived from the ∆h is obtained by Equation ( 18) based on the isentropic compression process as follows:

The Müller-II Model
Besides the Müller-I model, Müller et al. [35] also proposed another compressor mass flow rate model with simpler structure as shown in Equation (20), which is referred to as the Müller-II model in this paper: The theoretical foundation of the Müller-II model is that ∆h is it is mainly influenced by U 2 , meanwhile, .m has relatively less influence on ∆h is under low and medium rotational speed conditions, which happens to be a primary operating area for a compressor.

The Leufvén and Llamas Ellipse Model
The idea of the Leufvén and Llamas ellipse model is showed graphically in Figure 3 [3,16,28,29].By using the choking line (ChL) and the zero-slope line (ZSL), the whole operating area is divided into three zones: design operating zone, choking zone, and surging zone.ChL is the curve connecting the choking point of each iso-speed line, and ZSL is a curve that connects to the surging point.The iso-speed line between ZSL and ChL is represented by a super-ellipse with a curvature of CUR and crossing both the choking point ( .m ch , Π ch , and surging point ( .m zs , Π zs ).
m ch , Π zs , Π ch , and CUR are all functions of normalized compressor rotational speed N n (N n = N/N max ).
• Choking zone ( .m > .m ch ) When entering the choking zone, the mass flow rate will no longer change with the pressure ratio.Therefore, a straight line with a slope of infinity and crossing the choking point can be used to describe the changing trend of the iso-speed line in this zone: .m = .m ch (22) • Surging zone (0 m zs ) To reduce the model complexity, it was pointed out by Llamas and Eriksson [16] that a straight line with a slope of .m max /0.15Π max crossing the zero slope point ( .m zs , Π zs ) can be used to describe the changing trend of the iso-speed line in this zone: Actually, the ellipse model can also predict the pressure ratio at a given mass flow rate and rotational speed conditions, and it can thus simulate the surging phenomenon [3,16,28,29].

The Karlson Model
Karlson [6] put forward an isentropic efficiency model for a large-scale marine compressor based on a fifth-order polynomial function of .m and Π with a total number of 21 model parameters.It should be noted that the theoretical foundation of the Karlson isentropic efficiency model is based on multiple regression analysis; therefore, its extrapolation accuracy cannot be guaranteed to some extent.

The Guan Cong Model
In addition to the mass flow rate model for the marine compressor, Guan et al. [34] also proposed a zonal isentropic efficiency model, the theoretical foundation of which is the approximately linear relationship existing between the dimensionless torque coefficient Γ and flow coefficient φ.The definition Γ is as follows where ω (ω = 2πN/60) is the compressor angular velocity.
As similar to the mass flow rate model proposed by Guan et al. [34], the Guan Cong isentropic efficiency model also adopts the zonal modeling approach.As shown in Equation ( 25), the model parameters p i for the operating points within the current zone can be estimated through linear interpolation where N upper and p i,upper are the rotational speed, and corresponding model parameters for the upper boundary of the current zone; p i,lower and N lower are for the lower boundary.
To further improve the prediction accuracy, a correction procedure is carried out by the Guan Cong model based on the fact that the highest measured isentropic efficiency η max of each iso-speed line can be approximated with a second-order function of rotational speed where η non−corr is the predicted isentropic efficiency without correction; η non−corr,max is the maximum value of predicted uncorrected isentropic efficiency at a rotational speed of N; η corr is the corrected isentropic efficiency.

The Zeng Tao Model
Zeng et al. [36] proposed a mathematical model for predicting the compressor absorbed power m), are defined.It was revealed by Zeng et al. [34] that a strong second-order relationship exists between C power and C speed , and a strong linear relationship between the logarithm of C power and C speed .Consequently, the isentropic efficiency can be estimated at given conditions of mass flow rate, angular velocity and pressure ratio according to its definition.

The Hadef Model
Based on the compressor isentropic efficiency model proposed by Martin and Talon [20], Hadef et al. [19] made appropriate improvements and assumptions and proposed the Hadef model.The theoretical basis of the Hadef model is the application of the Euler equation in the velocity triangle at the inlet and outlet of the compressor blade where ∆h act is the actual specific enthalpy change across the compressor blade; ρ out , β out , and A out are the air density, outlet angle, and flow area at the outlet of compressor blade, respectively.Based on the assumptions made by Hadef et al. [19], ∆h act can be regarded as a linear function of mass flow rate at given rotational speed conditions, and thus Equation ( 27) can be rewritten as: After obtaining ∆h act at given conditions of rotational speed and mass flow rate and the specific isentropic enthalpy change ∆h is at given pressure ratio condition, the isentropic efficiency can be calculated with its definition: η = ∆h is /∆h act .

Llamas Model
Llamas and Eriksson [28] further improved the Hadef model and proposed the Llamas model.Compared with the Hadef model, the Llamas model takes into account the change of air density at the inlet and outlet of the compressor blade as well as the influence of friction loss on ∆h act .The relationship between ∆h act and .m and N n is expressed as follows where c is friction loss correction coefficient; ∆h act,max , .m max , and N max is the maximum value of the measured ∆h act , .m, and N, respectively, which is used to normalize the base functions a(N n ), b(N n ), and N. The derivation process of the base functions a(N n ) and b(N n ) is too long, which will not be introduced in this paper, and interested readers can refer to the original paper [28].

Test Compressors, Comparative Methodology, and Error Evaluation Criteria
To ensure unambiguous conclusions, the performance of various models was compared and analyzed for two marine compressors with different impeller outlet diameter, rotational, and flow ranges.The model numbers were A270-L59 and TCA88-25070, respectively, with their main technical parameters as shown in Table 1.It should be noted that a larger number of measured sample data exists in the performance map of the A270-L59 compressor relative to that of the TCA88-25070 compressor, which was 289 for the formal and 54 for the later, and in addition, the distribution of iso-speed lines in the former was more intensive with an interval of 600 RPM.Consequently, the parameterization and the predictive result were expected to be better when the A270-L59 compressor was selected as the research object and the relevant comparative and analysis result was thus believed to be more representative.Nevertheless, the results obtained based on the TCA88-25070 compressor can be used to validate the comparative and analysis result obtained based on the A270-L59 compressor, to some extent.As mentioned in Section 2, the whole operating area of a compressor is divided into design, LS, HS, and LP area, so the predictive and extrapolation ability of various models was investigated, compared, and analyzed in these four areas separately.For the design operating area, all the measured sample data available in the performance map were used to parameterize the models and then each iso-speed line was predicted by the models to compare with the measured ones.For the LS or HS area, the measured sample data at the lowest or highest iso-speed line in the performance map were removed when parameterizing the models, and then the removed iso-speed line was extrapolated by the models for comparison with the measured one.It should be noted that the choking point on each iso-speed line was only measured and provided in the performance map of the TCA88-25070 compressor, so the analysis on the extrapolation ability of each model in the LP area was only performed on this compressor.For a certain mass flow rate model, its LP area extrapolation ability lies in whether the choking flow can be predicted accurately; on the other hand, the LP area extrapolation ability of an isentropic efficiency can be evaluated by investigating whether efficiency will drop to zero when the pressure ratio across the compressor decreases to one.In addition, it should be noted that the measured mass flow rate rather than predicted one by mass flow rate model was used when comparing the predictive ability of each isentropic efficiency model.
In this paper, the Curve Fitting toolbox in Matlab was adopted to determine the model parameters with the least square method, and the values given in the original papers were referred to when determining the initial guesses on these model parameters.In order to quantify the predictive ability of various compressor mass flow rate and isentropic efficiency model with their model parameters determined with the same measurement data set, the following four criteria for error evaluation were adopted in this paper, including corrected coefficient of determination R 2 c , mean absolute percentage error MAPE, PEB ±5% , and PEB ±10% , that is percentage of total data points within the error bounds ± 5% and ± 10%, respectively, with their definitions as follows [12,13] where RD i is the relative deviation of ith measured data point; y meas is the measured value; y pred is the predicted value; y mean is the mean measured value; n is the number of measured data points; and m is the number of predictors (m = 4 in the Hadef model, for instance).
The reason why the corrected coefficient of determination R 2 c was chosen rather than the coefficient of determination R 2 as the criteria for error evaluation is that R 2 is influenced by the number of measured sample data and predictors [12,13].

Design Operating Area
Figures 4 and 5 presents the comparison between the measurement sample data and the predicted iso-speed lines for the two compressors, respectively.Tables 2 and 3 give the corresponding error evaluation results for each model.By observing the comparative and error evaluation results shown in Figures 4 and 5 and Tables 2  and 3, it can be found that in general, all the investigated models were capable of characterizing the changing trend between mass flow rate and pressure ratio for both compressors in the design operating area, but with different predictive accuracy.However, by comparing Tables 2 and 3, it can be seen that the predictive accuracy of each mass flow rate model dropped moderately when the research object changed from the A270-L59 to the TCA88-25070 marine compressor.This is perhaps due to the fact that a relatively smaller number of measured sample data exists in the performance map of the TCA88-25070 compressor relative to the A270-L59 compressor, which is 54 for the former, and 289 for the later, as mentioned in Section 4. Consequently, the overall predictive accuracy of the compressor models could be polluted to some extent in the existence of several predictive points with relatively large errors, like for the TCA88-25070 compressor.Nevertheless, in general, the ranking of these models in terms of predictive accuracy was found to be consistent for the two different marine compressors by analyzing the obtained error evaluation results, which demonstrates the validity of the results given in this paper, to some extent.
Among these models, the Guan Cong, Karlson-I, and Karlson-II models were the top three mass flow rate models with MAPEs around 1% for the A270-L59 compressor and around 2.5% for the TCA88-25070 marine compressor; additionally, they are also capable of characterizing the changing trend of the iso-speed line when approaching the choking point.Of the top three models, the Guan Cong model achieved the best predictive accuracy, mainly owing to its zonal modeling methodology.
Although it was revealed by Leufvén and Eriksson [3] that desirable predictive accuracy can be achieved with the Leufvén and Llamas ellipse model, the error evaluation results shown in Tables 2  and 3 indicate that the predictive accuracy of this model is undesirable for the two marine compressors investigated in this paper.This phenomenon is caused mainly due to its theoretical foundation: A super-ellipse is utilized to approximate the iso-speed line, and the maximum pressure ratio in the design operating area should not exceed the pressure ratio at ZSL, i.e., Π zs , otherwise, the straight line as shown in Equation ( 23) will be used to calculate the mass flow rate.Nevertheless, it cannot be guaranteed that the pressure ratio of all the measured sample data would be less than Π zs at a given rotational speed condition.As a result, for those sample data with a pressure ratio larger than Π zs the predicted mass flow rates will obviously deviate from the measured ones, which impairs the model predictive accuracy significantly.If those measured sample data are removed, the predictive accuracy of the Leufvén and Llamas ellipse model will improve greatly.Based on the above analysis, it can be concluded that the Leufvén and Llamas ellipse model has a relatively high requirement for the quality of measured sample data, thus weakening its practicability in engineering practices.
Despite the relatively lower predictive accuracy in the whole design operating area, it was found that for both the Müller-I and -II models, satisfactory predictive accuracy could be achieved under low rotational speed conditions, and the reason for this phenomenon is mainly owed to their theoretical foundation as introduced in Sections 3.1.5and 3.1.6.For the Malkhede model, its PEB ±10% was larger than 90%, and MAPE less than 5.5% for both compressors, so its predictive ability was considered to be satisfactory if the requirement on predictive accuracy is not extremely high.However, it can be observed in Figures 4 and 5, that after crossing a certain operating point, the mass flow rate will decrease with a decrease in pressure ratio, which is unrealistic and violates the actual changing trend.Although the Kolmanovsky model is capable of characterizing the changing trend of the iso-speed line when approaching the choking point, its predictive accuracy was relatively low among the investigated models, especially for the TCA88-25070 compressor.

Low-Pressure Ratio (LR) Area
Figure 6 presents the changing trend of choking pressure ratio Π ch and choking volume flow rate .
V ch with rotational speed for the TCA88-25070 marine compressor.As observed in Figure 6a, Π ch increased monotonously with the rotational speed and all the measured values of Π ch were greater than one.With a further decrease in rotational speed Π ch will decrease monotonously, and at a certain rotational speed, the condition of Π ch will become less than one.The research carried out by Leufvén and Eriksson [29] revealed that Π ch is approximately equal to 0.5 when the compressor rotational speed decreases to zero.Similarly, as observed in Figure 6b, .
V ch also increased monotonously with compressor rotational speed, but the increasing magnitude decreased gradually.Although the predictive accuracy of the Leufvén and Llamas ellipse model is relatively undesirable in the design operating area for both compressors, special base functions are defined in this model to characterize the changing trend of choking pressure ratio and choking mass (or volume) flow rate with the rotational speed.As seen in Figures 4 and 5, the Leufvén and Llamas ellipse model was able to describe the changing trend of the iso-speed line when approaching the choking point, its slope gradually approached infinity.In addition to the Leufvén and Llamas ellipse model, the Karlson-I, Karlson-II, and Kolmanovsky models utilize the curve characteristic of the exponential function to characterize the compressor choking phenomenon.As can be seen in Figures 4 and 5, for the Kolmanovsky model, the predicted iso-speed line entered into the choking area either too early or late; on the other hand, the Karlson-I and Karlson-II models appeared to have similar LP area extrapolation capabilities.So, in this section, only the Leufvén and Llamas ellipse models and the Karlson-I model were compared and analyzed in terms of LR area extrapolation ability.
Besides the measured data points, the choking pressure ratio predicted by the Leufvén and Llamas ellipse model is also depicted in Figure 6a.In Figure 6b, the choking volume flow rate predicted by the Leufvén and Llamas ellipse model and the Karlson-I model is depicted along with the measured sample data.As can be observed in Figure 6, both of the two models were able to accurately characterize the changing trend of choking volume flow rate at rotational speeds; in addition, the Leufvén and Llamas ellipse model was also capable of describing the changing trend of choking pressure ratio with speed, which gave it unique superiority relative to the Karlson-I model.V ch,exp by the Karlson-I model as well as the corresponding relative error RD ell and RD exp , under each speed condition.As observed in Table 4, RD ell is generally lower than RD exp , and the MAPE of the Leufvén and Llamas ellipse model is also lower than that of the Karlson-I model.In addition, since the Leufvén and Llamas ellipse model was capable of predicting the choking pressure ratio accurately and reasonably, the validity of the changing trend of the iso-speed line was guaranteed to some extent when approaching the choking point.In conclusion, although the predictive accuracy of the Leufvén and Llamas ellipse model is undesirable in the design operating area, its LR area extrapolation ability was the best among the compressor mass flow rate models investigated in this paper.

Low Speed (LS) Area
In Figures 7 and 8, the HS area extrapolation results of each model are compared with the A270-L59 and the TCA88-25070 compressor as modeling objects, respectively, and Tables 5 and 6 provide corresponding error evaluation results.It can be found from the comparative and error evaluation results that the LS area extrapolation ability of the Guan Cong, Karlson-I, Karlson-II, Müller-I, and Müller-II models is obviously better than the other models when the A270-L59 compressor was selected as modeling object.However, it dropped significantly for both the Karlson-I and Karlson-II models for the TCA88-25070 compressor, which is maybe due to the relatively smaller of measured data points and the sparse iso-speed lines available in the performance map of this particular compressor.
Among these models, the Guan Cong model manifests a satisfactory LS area extrapolation accuracy for both compressors with a MAPE of less than 5.5%, which is mainly attributed to its zonal modeling methodology as well as the fact that air is effectively incompressible under low-speed conditions, and thus, the changing trend between mass flow rate and pressure ratio will not change significantly.It is interesting to note that the Müller-I model achieved the best extrapolation ability for both compressors with a MAPE of 1.5994% for the A270-L59 and 3.5147% for the TCA88-25070, which was mainly due to its theoretical foundation.However, it should be noted that this model belongs to the second type of mass flow rate model as shown in Equation ( 8), which is not suitable for the dynamic simulation modeling of turbocharged diesel engines.It can be observed from Figures 7 and 8 that several predicted values by the Leufvén and Llamas ellipse model significantly deviate from the measured ones because of the greater pressure ratio than that at ZSL.For the Kolmanovsky and Malkhede models, their LS area extrapolation accuracy was undesirable, which will be not discussed further.

High Speed (HS) Area
It can be observed from the comparative and error evaluation results, as shown in Figures 9 and 10 and Tables 7 and 8, that all the investigated models achieved satisfactory HS area extrapolation accuracy except for the Guan Cong and Ellipse models.It is interesting to note that relative to its LS area extrapolation accuracy, the Guan Cong model dropped significantly in the HS area with a MAPE as high as 28.7346% for the A270-L59 compressor and 17.5797% for the TCA88-25090 compressor; in addition, for the TCA88-25090 compressor, the predicted mass flow rate decreased with a decrease in pressure ratio, which is unrealistic and violates the actual changing trend.The reason causing the inferior HS area extrapolation accuracy for the Guan Cong model may be due to the fact that when extrapolating to the HS area, the model parameters corresponding to the highest iso-speed lines are still adopted, which cannot capture the influence of rotational speed change on the fluid physical properties, especially the compressibility.Similar to the LS area, several predicted values by the Leufvén and Llamas ellipse model significantly deviated from the measured values.In contrast to the LS area, satisfactory extrapolation accuracy was achieved by both the Karlson-I and Karlson-II models in the HS area, with a MAPE less than 2.6% for the A270-L59 compressor and 5% for TCA88-25070 compressor, which demonstrates their stable and robust HS area extrapolation capability to some extent.

Design Operating Area
Figures 11 and 12 present the prediction results for each compressor isentropic efficiency model along with the measured data points for the A270-L59 and TCA88-25070 compressors, respectively.Tables 9 and 10 provide the corresponding error evaluation results.In order to clearly display and compare the prediction results, only part of the iso-speed lines are depicted in Figures 11 and 12     By observing and analyzing the comparative and error evaluation results, it was found that in general, all the investigated models are able to characterize the changing trend of isentropic efficiency with the mass flow rate for both the A270-L59 and TCA55-25070 marine compressors.Meanwhile, satisfactory prediction accuracy was attained for all the models and compressors with a PEB ±10% equal to 100%, PEB ±5% greater than 89%, and a MAPE less than 2.8%.Among these models, the Guan Cong model achieved the best prediction accuracy for both compressors, which is mainly attributed to its zonal modeling methodology and the correction procedure adopted, as introduced in Section 3.2.2.Although the predictive accuracy of both the Hadef model and the Llamas model were slightly inferior to the Guan Cong model, they were still satisfactory, with a MAPE of less than 1% and a PEB ±5% approaching 100% for both compressors, and this is attributed to the excellent parameterization results achieved with R 2 approaching one, and the strong physical fundamentals embedded in the two models.It should be noted that although the Llamas model can be regarded as an improved version of the Hadef model, the improvement level is not very evident for the two marine compressors investigated in this research, as can be observed in Tables 9  and 10.Despite the excellent parameterization results achieved for both the first-and second-order Zeng Tao models, with R 2 approaching one, their predictive accuracy was inferior compared to the Guan Cong, Hadef, and Llamas models with a MAPE around 1.6% for the TCA55-25070 compressor, and around 2.5% for the A270-L59 compressor.Nevertheless, it should be noted that the predictive accuracy of the Zeng Tao model was still satisfactory if the requirement on predictive accuracy was not extremely high, because the PEB ±10% of this model was still 100% and the MAPE was between 1.5% and 2.8%, for both compressors.In addition, the Zeng Tao model presents strong practicability in engineering practice, as it can be parameterized with a relatively small number of measured sample data.The predictive accuracy of the Karlson model dropped moderately when the modeling object switched from the A270-L59 to the TCA55-25070 compressor; in addition, the R 2 c of this model was only 0.3939 for the TCA55-25070 compressor, and this is due to the fact that the number of model parameters existing in the Karlson model is too much with respect to the amount of measured sample data existing for this particular compressor (21 model parameters and 54 measured sample data).

Low-Pressure Ratio (LR) Area
Similar to the compressor mass flow rate model, the isentropic efficiency model is also required to extrapolate to the LR area.In general, only operating points with relatively high isentropic efficiency are measured and depicted in the compressor performance map, for instance, the lowest isentropic efficiency value in the performance map of A270-L59 marine compressor is as high as 65%, and it is 74% for the TCA55-25070 compressor.Due to the absence of corresponding measured data points, it is difficult to study the LR area extrapolation ability for each isentropic efficiency model.
For the isentropic efficiency model, the compressor rotational speed, pressure ratio, and mass flow rate are generally selected as the input variables, and it should be noted that η will be necessarily zero when Π is equal to one according to its definition.Although the mass flow rate when Π is equal to 1 is generally not provided in the performance map, its value .m Π=1 can be estimated with the help of Leufvén and Llamas ellipse model that manifests superior LR area extrapolation ability.Based on the above analysis, the LR area extrapolation ability of each model can be investigated by comparing η = 0 and the predicted isentropic efficiency η Π=1 with Π = 1 and corresponding .m Π=1 as model input at each rotational speed condition.
Since the definition of isentropic efficiency is directly utilized to estimate the isentropic efficiency of the Zeng Tao, Hadef, and Llamas models, therefore the predicted isentropic efficiency will necessarily be zero when Π = 1; in addition, as the term "Π (k−1)/k − 1" is contained in the definition of the dimensionless torque coefficient Γ,as shown in Equation ( 24), a similar result would also be obtained by the Guan Cong model.As a result, only the LR area extrapolation ability of the Karlson model was required to be investigated.As observed in Tables 11 and 12, η Π=1 predicted by the Karlson model deviated from η = 0 significantly for both compressors, and some of the predicted values were even greater than one or negative, which violates the physical meaning of the isentropic efficiency.The inferior LR area extrapolation ability of the Karlson model is mainly due to its theoretical foundation, which is based on multiple regression and lacks strong physical fundamentals.In Figures 13 and 14, the LS area extrapolation results of each model are compared with the A270-L59, and the TCA88-25070 compressors as modeling objects, respectively, and Tables 13 and 14 provide corresponding error evaluation results.fundamentals.Although the HS area extrapolation accuracy of the Guan Cong model dropped slightly when compared to its LS area extrapolation accuracy, it still achieved the best predictive accuracy among these models in the HS area.
By observing the error evaluation results shown in Tables 9-16, it can be seen that for the Zeng Tao, Hadef and Llamas models, their predictive accuracy was relatively stable and satisfactory in both the design and off-design areas.This is the consequence of the relatively strong physical fundamentals embedded in these models.

Summary and Discussion
In this paper, a comprehensive applicable and comparative study of compressor mass flow rates and empirical isentropic efficiency models was carried out and performed on two marine compressors with different impeller outlet diameters, flow ranges, and rotational speed ranges.In general, consistent comparative analysis results were obtained for the two different compressors, and they are summarized as follows: Compressor mass flow rate model: • In general, In the design operating area, all of the investigated models were capable of characterizing the changing trend between mass flow rate and pressure ratio, among which the Guan Cong model achieved the best predictive accuracy, mainly owing to its zonal modeling methodology; • Besides the design operating area, superior LS and LP extrapolation accuracy were also achieved by the Guan Cong model.However, it dropped significantly in the HS area for both compressors due to the considerable influence of rotational speed on the fluid physical properties under high-speed conditions; • A satisfactory predictive result was achieved by both the Müller-I and -II models under low rotational speed conditions, which is mainly attributed to its theory foundation; • The Leufvén and Llamas ellipse model was capable of describing the compressor choking phenomenon accurately and reasonably, including both the choking flow and choking pressure.
In addition, the Guan Cong, Karlson-I, Karlson-II, and Kolmanovsky models could also describe the choking phenomenon, but they were inferior to the ellipse model;

•
It is interesting to note that satisfactory predictive and extrapolation accuracy was achieved by both the Karlson-I and -II models for the A270-L59 compressor.However, the HS area extrapolation ability of the two models deteriorated significantly when the TCA88-25070 compressor was selected as the modeling object, which indicates that their HS area extrapolation ability is unstable and may vary with compressors.
Compressor isentropic efficiency model: • When extrapolating to LP and HS area, the Karlson model manifested an undesirable extrapolation ability mainly due to its theoretical foundation, which is based on multiple regression and lacks strong physical fundamentals; • In the design operating area, the predictive accuracy of the Guan Cong model was the highest, mainly owing to its zonal modeling methodology and the adopted correction procedure.Additionally, this model also had satisfactory LP, LS, and HS extrapolation abilities; • Owing to the strong quadratic relationship existing between the power and speed coefficients, and the strong linear relationship between the logarithm of power and speed coefficients, the Zeng Tao model was very practical in engineering practice, as it can be parameterized readily with a relatively small number of sample data points.

•
For the whole operating area, both the Hadef model and the Llamas model provided satisfactory predictive and extrapolation ability for both compressors investigated in this paper, and this is because that the theoretical foundation of the two models are based on basic physical laws (the Euler equation for turbomachinery and the laws of thermodynamics, for instance) and the definition of isentropic efficiency is directly utilized.In addition, the predictive and extrapolation ability of the Llamas model was not improved, although obviously, it is regarded as an improved version of the Hadef model.
It can be found from the applicable and comparative research that all the models investigated in this paper belong to the empirical model category, and this is mainly due to the fact that this type of model has a simpler structure and is easier to tune, and its predictive and extrapolation ability was satisfactory to some extent.Consequently, empirical compressor models are widely adopted in the study of performance, energy analysis and control-or diagnosis-oriented simulations.Nevertheless, due to relatively weaker physical fundamentals embedded in these empirical models, relative to their theoretical counterparts, satisfactory predictive and extrapolation ability cannot be guaranteed for an arbitrary compressor.For example, the LS area extrapolation accuracy of both the Karlson-I and -II models dropped significantly when the modeling object switched from the A270-L59 to the TCA88-25070 compressor, as shown in Section 5.2.2.Undoubtedly, the research carried out in this paper will be more meaningful if these empirical models are additionally compared with theoretical models in terms of simulation speed, predictive and extrapolation accuracy, and complexity in model tuning.However, a detailed compressor geometry dimension is always required for developing theoretical models, which is usually difficult to obtain for large-scale marine compressors including the two compressors investigated in this paper.Actually, the ultimate purpose of this paper was to find a compressor model that is most suitable for the modeling of marine two-stroke diesel engines used for marine engineering simulation, which has a relatively high requirement on both simulation speed and predictive accuracy, and it can be seen from the comparative analysis results that some of the empirical models satisfied this requirement very well.
Another point that should be noted is that none of these empirical compressor models appear to achieve satisfactory predictive and extrapolation accuracy in the whole operating area simultaneously.For instance, although extremely satisfactory LR area extrapolation accuracy was achieved by the Leufvén and Llamas ellipse mass flow rate model, its predictive ability in design operating area and extrapolation ability in LS and HS area was undesirable.Actually, this limitation can be broken with a zonal modeling approach, i.e., selecting the model with the best predictive accuracy for each operating area.For example, for the A270-L59 compressor, the design and LS operating area can be modeled by the Guan Cong model, LP area modeled by the ellipse model, and HS area modeled by the Karlson-II model.However, it should be noted that an appropriate blending algorithm is required to guarantee a smooth transition at the border between different operating areas.The authors will follow this path and construct a set of empirical models, which can achieve satisfactory predictive accuracy in the whole operating area.
Generally, for achieving satisfactory predictive results with these empirical models, a large number of measured sample data is required for model tuning and parameterization.For a small-scale automotive compressor, its performance map can be readily obtained from the manufacturers, or even be found on the internet.However, as large-scale marine compressor turbocharger manufacturers rarely provide its performance maps, it makes the model parameterization procedure very difficult.Nevertheless, the shop test report of the marine diesel engine usually provides compressor operating performance at several engine loading points, including rotational speed, pressure and temperature at the inlet and outlet, and sometimes mass flow rate.Therefore, it will be very practical to develop a kind of compressor model that can be parameterized only with these measured performance data provided by the shop test report.An approximate methodology, following this idea proposed by Theotokatos [37], is often adopted by authors when developing simulation models for marine diesel engines in cases where the compressor performance map is lacking.This method is based on the fact that the engine is loaded according to the propeller law in a marine propulsion plant, and the operating points of the compressor lie on a single curve on the compressor map under steady-state conditions; in addition, this curve is close to the optimum compressor efficiency parabola owing to the fine matching between the turbocharger and engine.Consequently, the compressor pressure ratio and isentropic efficiency can be modeled as second-order polynomial functions of the turbocharger rotational speed with only three model parameters at most, which can be then parameterized readily with the compressor operating performance at several loading points.Similarly, as strong quadratic relationship exists between the power and speed coefficient and strong linear relationship between the logarithm of power and speed coefficient in the Zeng Tao isentropic efficiency model [36], the R-square of fitting results of which was more than 0.99 for both the A270-L59 and the TCA88-25070 marine compressors.Thus, this model can be parameterized readily only with the compressor performance information at the five to six operating points provided by shop test reports.In general, with fewer model parameters existing in the empirical models, the model parameterization procedure becomes much easier.However, it should be noted that the parameterization result is expected to be more reliable with more measured sample data available.

Conclusions
A comprehensive applicable and comparative study of compressor mass flow rate and isentropic efficiency empirical models was carried out on two large-scale marine compressors.Based on the comparative analysis results obtained in this research, several guidelines are summarized, which can be followed when developing simulation models for both compressor mass flow rate and isentropic efficiency:

•
Different types of empirical compressor models manifest different predictive accuracy and extrapolation ability in different operating areas; therefore, it is necessary to select the most suitable modeling methodology for the particular area concerned;

•
The zonal modeling approach proposed and implemented by Guan et al. [34] is able to effectively improve the predictive accuracy for both compressor mass flow rate and isentropic efficiency models.However, its LS and HS area extrapolation strategy requires special attention as the change in compressor rotational speed has influences on fluid physical properties especially compressibility under high-speed conditions; • To ensure robust extrapolation ability and reliable extrapolation results, it is indispensable to enhance the physical fundamentals of the developed compressor models; • When developing compressor models, it is necessary to pay special attention to their sensitivity to the errors in model parameters, and this is due to the fact that if the sensitivity is high, satisfactory predictive accuracy cannot be guaranteed even though the parameterization result is satisfactory, with a high value of R 2 c for instance; • Since the mass flow rate is normally selected as the input variable for the isentropic efficiency model, the predictive errors of the mass flow rate model will propagate to the isentropic efficiency model.So, it will be beneficial to parameterize the mass flow rate and isentropic efficiency model simultaneously with an appropriate parameterization algorithm in order to achieve satisfactory predictive accuracy in the three dimensions, including the pressure ratio, mass flow rate, and isentropic efficiency.
Both the A270-L59 and TCA88-25070 compressors are very representative of the air boosting of large two-stroke marine diesel engines, and consistent comparative analysis results were obtained for the two different compressors.Thus, the research conducted in this paper is believed to be very helpful for researchers who are interested in the dynamic simulation study of large two-stroke marine diesel engines to select the best mathematical models for their application and research.
Author Contributions: H.S. contributed to the parameterization of mathematical models, analysis of simulation results, and the writing of the manuscript.C.Z. provided language help and some technical support.J.Z. provided research conditions, team support, and financial support.B.Y. and B.J. gave some useful suggestions.All authors have read and agreed to the published version of the manuscript.Funding: This work was supported by the National Natural Science Foundation of China, grant number 51479017, and the Fundamental Research Funds for the Central Universities, grant number 3132019315.

Figure 2 .
Figure 2. Two-dimensional interpolation meshes constructed by look-up table method: (a) mass flow rate model; (b) isentropic efficiency model.

. 1 . 5 .
The Müller-I Model Müller-I model was developed by appropriately simplifying the compressor 1-D theoretical model developed by Müller et al. [35].For the Müller-I model, the specific isentropic enthalpy change ∆h is is expressed as a function of compressor blade tip speed U and mass flow rate .m

Figure 3 .
Figure 3. Sketch of the main characteristics of the Leufvén and Llamas ellipse model.• Design operating zone ( .m zs ≤ .m ≤ .m ch ) In the design operating zone, a super-ellipse with a curvature of CUR is used to describe the relationship between .m and Π

.W 3 )
with the feature of low complexity.The theoretical basis of the Zeng Tao model is based on the Euler equation for turbomachinery.Meanwhile, various types of power losses caused by such as slip, windage and friction are also taken into account.In the Zeng Tao model, the two performance parameters power coefficient C power (C power = .and speed coefficient C speed (C speed = ω/ .

Figure 6 .
Figure 6.Measured and predicted choking pressure and volume flow rate: (a) choking pressure; (b) choking volume flow rate.
Figures 11 and 12  present the prediction results for each compressor isentropic efficiency model along with the measured data points for the A270-L59 and TCA88-25070 compressors, respectively.Tables9 and 10provide the corresponding error evaluation results.In order to clearly display and compare the prediction results, only part of the iso-speed lines are depicted in Figures11 and 12.

Figure 11 .
Figure 11.Prediction results of compressor isentropic efficiency model in design operating area for the A270-L59 compressor: (a) the Karlson model; (b) the Guan Cong model; (c) the Second-order Zeng Tao model; (d) the First-order Zeng Tao model; (e) the Hadef model; (f) and the Llamas model.

Figure 12 .
Figure 12.Prediction results of each compressor isentropic efficiency model in the design operating area for the TCA88-25070 compressor: (a) the Karlson model; (b) the Guan Cong model; (c) the Second-order Zeng Tao model; (d) the First-order Zeng Tao model; (e) the Hadef model; (f) and the Llamas model.

Table 1 .
Main technical parameters of the A270-L59 and the TCA88-25070 marine compressors.

Table 2 .
Error evaluation results of each compressor flow rate model in the design operating area for the A270-L59 compressor.

Table 3 .
Error evaluation results of each compressor mass flow rate model in the design operating area for the TCA88-25070 compressor.

Table 4
presents the measured choking volume flow rate .V ch,meas , the predicted values .V ch,ell by the Leufvén and Llamas ellipse model, and .

Table 4 .
Measured and predicted choking volume flow rates and prediction errors at each rotational speed condition for the TCA88-25070 compressor.

Table 5 .
Error evaluation results of the LS area extrapolation results of each compressor mass flow rate model for the A270-L59 compressor.

Table 6 .
Error evaluation results of LS area extrapolation results of each compressor mass flow rate model for the TCA88-25070 compressor.

Table 7 .
Error evaluation results of HS area extrapolation results of each compressor mass flow rate model for the A270-L59 compressor.

Table 8 .
Error evaluation results of HS area extrapolation results of each compressor mass flow rate model for the TCA88-25070 compressor.

Table 9 .
Error evaluation results of each compressor isentropic efficiency model in the design operating area for the A270-L59 compressor.

Table 10 .
Error evaluation results of each compressor isentropic efficiency model in the design operating area for the TCA88-25070 compressor.

Table 11 .
Predicted isentropic efficiency value from the Karlson model when pressure ratio equals to 1 for the A270-L59 compressor.

Table 12 .
Predicted isentropic efficiency value from the Karlson model when pressure ratio equals to 1 for TCA55-25070 compressor.