Development of a Throughﬂow-Based Simulation Tool for Preliminary Compressor Design Considering Blade Geometry in Gas Turbine Engine

: Gas turbine engines are highly intricate machines, and every component of them is closely associated with one another. In the traditional engine developing process, vast experiment tests are needed. To reduce unnecessary trials, a whole gas turbine engine simulation is extremely needed. For this purpose, a compressor simulation tool is now developed. Considering the inherent drawbacks of 0D analysis and 3D CFD (Computational Fluid Dynamics) calculation, the 2D throughﬂow method is an indispensable tool. Based on the circumferential average method (CAM), 3D Navier–Stokes is transformed into a 2D method. One phenomenon arising is that the lack of description about circumferential motion leads to the need for the blade force modeling in compressor simulation. Previous models are based on the assumption that ﬂow passes through the average stream surface without entropy increasing, which is not applicable in the CAM. An improved model is proposed based on the result analysis from CAM and NUMECA method in a linear cascade. Whereafter, the model is applied in a highly loaded and low-speed fan, which has been tested for its performance characteristics. Utilizing the new model, the error of the adiabatic efﬁciency between CAM and experiment decreases from 4.0% to 1.0% and the accuracy of the mass ﬂow, and pressure ratio remains unchanged. The time involved in the CAM simulation is nearly 70 times faster than that of the 3D simulation.


Introduction
When a new gas turbine engine is under development, abundant hardware builds are necessary to measure the performance characteristics, which are much too expensive [1]. Therefore, the simulation of an overall gas turbine engine is extremely needed. Thanks to both the advance of the CFD (Computational Fluid Dynamics) and the increasingly developed computer technology, fully three-dimensional (3D) computation is increasingly involved [2,3]. However, it is very time-demanding. For example, when Medic et al. [4] simulated a 20 • sector gas turbine engine (shown as Figure 1) on the coarse grids containing 14 million cells, a simulation was performed in nearly 2 weeks relying on 700 processors [5][6][7][8]. On the other hand, the zero-dimensional method (0D) commonly used in the performance evaluation cannot present the radial distribution of aerodynamic parameters. Therefore, the throughflow method (2D) is still a significant part of its performance prediction and rapid flow field acquisition [9]. To achieve this ultimate goal of the whole gas turbine engine simulation, the throughflow method applied in compressors is first explored in this paper.
The most widely used throughflow method is the streamline curvature method. Comparing with it, the circumferential average method (CAM) first, proposed by Smith [10] is appealing in the following aspects: choke mass flow can be predicted, certain types of shock can be captured, and the same code can be used for subsonic, transonic and supersonic meridional velocities. Much research has been performed based on the CAM in turbomachinery. However, one element indispensably mentioned in every throughflow method applied in turbomachinery is that the absence of the description about the circumferential motion leads to the blade force missing (shown as Figure 2), which needs the blade force modeling appended. As for the inviscid blade force model, it was first presented by Baralon [11]. In his method, the blade force modeling was proposed based on the fully time-dependent axisymmetric Euler equations. In this theory, the flow conformed to the average stream surface for each blade row, and the blade force must be orthogonal to the average stream surface in order to obtain a loss-free flow turning. Moreover, the blade force, which they thought determined the flow path, was its tangential component. Based on the theory of Baralon, Simon [12] modified the calculation formula of the inviscid blade force, where the tangential component of the inviscid blade force was substituted by the modulus of the inviscid blade force in the calculation formula. For the design mode of turbomachinery, it was first proposed by Sturmayr and Hirsch [13], where the target tangential velocity was prescribed through the flow angle at the trailing edge (TE). In their research, how the circumferential velocity evolves inside the blade passage in the streamwise direction is needed to be given. Dawes [14] adopted another method in the solution of the inviscid blade force, where the blade force was computed from the circumferential momentum and was not solved anymore in the bladed regions. The aforesaid methods depended on a formulation, which had a constant coefficient determining the convergence effect. Based on these methods above, convergence problems often arose when iteratively solving these inviscid blade force formulations, especially with an inappropriate relaxation factor. Considering this, Jin [15] presented a method similar to the large-particle method to avoid solving the inviscid blade force explicitly. This method consisted of two steps for each time step. First, the effects caused by the inviscid blade force in the momentum equations were neglected. Moreover, then only the effects caused by the inviscid blade force were taken into account in the governing equations. However, all these foregoing methods are dependent on the assumption of the average stream surface. Moreover, the average stream surface is calculated through the mean camber surface with the deviation angle. That is, when the flow passes through the average stream surface, the inviscid blade force is to deflect the flow, and no entropy increasing occurs. Therefore, it requires that the inviscid blade force be normal to the average stream surface. However, this assumption is not applicable in the CAM, which needs to take into consideration the effect of the blade geometry on the inviscid blade force. In this paper, the modeling of the inviscid blade force considering the effect of the blade geometry is discussed and verified.
However, this assumption is not applicable in the CAM, which needs to take into consideration the effect of the blade geometry on the inviscid blade force. In this paper, the modeling of the inviscid blade force considering the effect of the blade geometry is discussed and verified. This paper focuses on two parts of the CAM derived from the Navier-Stokes (N-S) equation. The first part is about the deduction of the inviscid blade force concerning the blade geometry from the equation and its modeling process. The second part is about the application of this modeling. It is successively applied in the throughflow simulation of a subsonic linear cascade, a transonic linear cascade. All these results are compared with 3D simulation results. Finally, it is involved in the simulation of a highly loaded and lowspeed fan named TA36 and calculation results are compared with the experimental data as well as CFD results.

Throughflow Analysis Tool
In this research, the throughflow model used is the CAM to derive the governing equations. For the 3D flow in aviation turbomachinery, the CAM method applied for any parameter is defined as where q is an aerodynamic parameter, q is the circumferential value, s ϕ and p ϕ are angular values for suction side and pressure side individually. In a vane passage shown in Figure 3, the integral equation starts from the suction side to the pressure side. A blockage factor due to the circumferential blade thickness takes on the form: This paper focuses on two parts of the CAM derived from the Navier-Stokes (N-S) equation. The first part is about the deduction of the inviscid blade force concerning the blade geometry from the equation and its modeling process. The second part is about the application of this modeling. It is successively applied in the throughflow simulation of a subsonic linear cascade, a transonic linear cascade. All these results are compared with 3D simulation results. Finally, it is involved in the simulation of a highly loaded and low-speed fan named TA36 and calculation results are compared with the experimental data as well as CFD results.

Throughflow Analysis Tool
In this research, the throughflow model used is the CAM to derive the governing equations. For the 3D flow in aviation turbomachinery, the CAM method applied for any parameter is defined as q(x, r, t) = 1 where q is an aerodynamic parameter, q is the circumferential value, ϕ s and ϕ p are angular values for suction side and pressure side individually. In a vane passage shown in Figure 3, the integral equation starts from the suction side to the pressure side.
However, this assumption is not applicable in the CAM, which needs to take into consideration the effect of the blade geometry on the inviscid blade force. In this paper, the modeling of the inviscid blade force considering the effect of the blade geometry is discussed and verified. This paper focuses on two parts of the CAM derived from the Navier-Stokes (N-S) equation. The first part is about the deduction of the inviscid blade force concerning the blade geometry from the equation and its modeling process. The second part is about the application of this modeling. It is successively applied in the throughflow simulation of a subsonic linear cascade, a transonic linear cascade. All these results are compared with 3D simulation results. Finally, it is involved in the simulation of a highly loaded and lowspeed fan named TA36 and calculation results are compared with the experimental data as well as CFD results.

Throughflow Analysis Tool
In this research, the throughflow model used is the CAM to derive the governing equations. For the 3D flow in aviation turbomachinery, the CAM method applied for any parameter is defined as where q is an aerodynamic parameter, q is the circumferential value, s ϕ and p ϕ are angular values for suction side and pressure side individually. In a vane passage shown in Figure 3, the integral equation starts from the suction side to the pressure side. A blockage factor due to the circumferential blade thickness takes on the form: A blockage factor due to the circumferential blade thickness takes on the form: where N is the number of blades. After the definition of the CAM, aerodynamic parameters can be decomposed into circumferential average and spatial fluctuation:

of 22
Based on the CAM and the main operation rule, the N-S equation in the relative frame can be derived into: where U is the vector of conservative variables, F and G are the convective flux vectors, F v and G v are the viscous flux vectors, S is the source term, F B is the inviscid blade force and F F is the viscous blade force. ρ is density, p is pressure; w x , w r ,w u are scalar values for relative velocity in x, r and ϕ direction, respectively; e and i are the relative total energy and total enthalpy per unit individually and they take on the form: where c v is the specific heat at a constant volume. In Equation (4), the first row represents the continuity equation, rows 2-4 depict the momentum equations, and the last row describes the energy balance. A time-marching finite volume method is used to solve the equation. The Edwards' LDFSS [16] scheme is applied to discretize the convective fluxes. The calculation of the viscous fluxes requires the evaluation of the primitive variables and their first derivatives at the edge middle of the considered cell. The Green-Gauss theorem is also used to estimate the gradient of the variables. An explicit Runge-Kutta method is chosen for the time discretization as it is straightforward to implement. Meanwhile, local time step and implicit residue averaging are used to accelerate the convergence.

Development of the Inviscid Blade Force
In previous research, the inviscid blade force was usually modeled as a term related to time and then solved through the iteration method. The formula proposed by Sturmayr and Hirsch [17] is: Appl. Sci. 2021, 11, 422 5 of 22 where τ is the pseudo time; C is a constant value used in iteration to regulate the convergence of iteration; w u target is the target tangential velocity; w u is the tangential velocity determined in the throughflow method.
Based on this method, researchers make several improvements upon the inviscid blade force. For example, Simon [18] put forward a model where f B is the magnitude of the inviscid blade force; n Bx ,n Br and n Bu are components of the unit normal vector of the average stream surface. The modeling of all these inviscid blade force models is based on the assumption that the effect produced by the inviscid blade force is to change the flow direction without entropy production. Therefore, the direction of the inviscid blade force must be orthogonal to the average stream surface. However, these models do not take the effect of the blade geometry into consideration, which is not suitable in the CAM.
According to Formula (4), the inviscid blade force takes the form Based on the coordinate definition of the mean camber surface: The circumferential coordinates of the suction side and the pressure side can be displayed as: and r ∂ϕ M ∂x , r ∂ϕ M ∂r , −1 T represents the normal vector of mean camber surface for a blade. Therefore, the inviscid blade force can be transformed as: That is, It can be seen that the axial and the radial components of the inviscid blade force include two parts. First, the first item is related to the circumferential inviscid blade force, and it is easy to know: where γ M and λ M are the axial flow angle and the radial flow angle, respectively, for the mean camber surface. Therefore, the axial and radial components of the inviscid blade force can be transformed into: Through this formula, it is revealed that the axial and radial components of the inviscid blade force are not only related to the circumferential component of the inviscid blade force and the flow angle of the blade camber surface but also are concerned with the blockage factor and the static pressure on the blade surface. This aspect has not been concerned in the former models. When the thickness of the blade changes rapidly along the axial or the radial direction, the inviscid blade force is no longer normal to the mean camber surface. Given this, the second component needs more consideration. However, the first item can be calculated through the circumferential component of the inviscid blade force and flow angle, whereas the second item needs to be modeled. As this study is about the directive problem of turbomachinery, the partial derivative of the blockage factor along the axial and the radial direction can be calculated. However, the pressure on the suction and pressure side cannot be acquired in CAM, which needs to be modeled.

3D CFD Method Verification
In this study, the analyses of the CAM results are based on the 3D simulation results. To perform the 3D simulation in this study, the commercial numerical software NUMECA was used to simulate the linear cascade and the TA36. The mesh generation for the calculation was completed by the IGG module of the NUMECA. The Fine module of the NUMECA is widely used to calculate the steady flow in turbomachinery. To validate the applicability of the 3D CFD method, the simulation of a controlled diffusion airfoil (CDA) linear cascade was carried out.
The experiment results of the linear cascade have been published in the works of Steinert [19]. In their papers, the coordinates of the blade geometry have been listed. Based on this, the calculation model was constructed. The H-O-H topology was adopted for the linear cascade to get a better mesh quality. After the grid independence check, the total number of grid points is nearly 7.1 × 10 5 , shown in Figure 4. Moreover, it is observed that the y+ is less than five after the calculation results are obtained. In this study, the Spalart-Allmaras (S-A) turbulence model was applied. At the inlet, the flow direction is axial, the total temperature is 288.15 K, and the total pressure is 101,325 Pa. At the outlet, the boundary condition is adjusted to make sure that the Mach number is the same as the experiment condition. No-slip condition is applied at the wall surface. total temperature is 288.15 K, and the total pressure is 101，325 Pa. At the outlet, the boundary condition is adjusted to make sure that the Mach number is the same as the experiment condition. No-slip condition is applied at the wall surface.
(a) (b) When the convergent results were acquired, the pressure distribution on the blade surface was compared with the experiment results, shown in Figure 5. From this figure, the distribution is nearly the same except for some points. Therefore, it can be reached that the NUMECA software is suitable to be applied in the further simulation of 3D models.  the distribution is nearly the same except for some points. Therefore, it can be reached that the NUMECA software is suitable to be applied in the further simulation of 3D models.  When the convergent results were acquired, the pressure distribution on the blade surface was compared with the experiment results, shown in Figure 5. From this figure, the distribution is nearly the same except for some points. Therefore, it can be reached that the NUMECA software is suitable to be applied in the further simulation of 3D models.

Modeling of the Inviscid Blade Force Based on a Linear Cascade
To investigate the inviscid blade force in the CAM, a simple linear cascade is involved. Based on this model, the other relevant factors which may influence the axial inviscid blade force can be excluded. Since no change occurs in the spanwise direction, the radial blade force Br F equals 0. A CDA linear cascade is applied in the study, and its geometry parameters are shown in Table 1:

Modeling of the Inviscid Blade Force Based on a Linear Cascade
To investigate the inviscid blade force in the CAM, a simple linear cascade is involved. Based on this model, the other relevant factors which may influence the axial inviscid blade force can be excluded. Since no change occurs in the spanwise direction, the radial blade force F Br equals 0. A CDA linear cascade is applied in the study, and its geometry parameters are shown in Table 1: In this paper, when studying this model, Euler equations are adopted in the CAM, which means that viscous blade force and viscous stress are neglected. In the simulation of the 3D method, Euler equations are also involved, and they are solved using NUMECA. The boundary conditions for these two methods are the same. At the inlet, the flow direction is axial, the total temperature is 288.15 K, and the total pressure is 101,325 Pa. At the outlet, the mass flow is set to ensure the Mach number at the inlet. No-penetration condition is applied at the wall surface.

Grid Independence Verification
To prove the reliability of the calculation results, it is necessary to verify the grid independence before the study. In the CAM, three types of mesh with different numbers of grid points were generated to calculate the model. The topology of the mesh and the boundary conditions were set the same. Moreover, the grid numbers are 41 × 33, 65 × 49, 97 × 73 (axial direction radial direction), respectively. The linear cascade was simulated on the subsonic condition. Based on results from these three grids, the distribution of static pressure and the Mach number are shown in Figure 6. The curves of the case with 65 × 49 grid points coincide with that of the case with 97 × 73 grid points. Hence, in conclusion, when the number of grid points exceeds 65 × 49, the results are considered independent of the grid number. In consideration of the calculation accuracy and the requirement of saving time, the numerical method with 65 × 49 grid points was introduced into this study. 97 × 73 (axial direction radial direction), respectively. The linear cascade was simulated on the subsonic condition. Based on results from these three grids, the distribution of static pressure and the Mach number are shown in Figure 6. The curves of the case with 65 × 49 grid points coincide with that of the case with 97 × 73 grid points. Hence, in conclusion, when the number of grid points exceeds 65 × 49, the results are considered independent of the grid number. In consideration of the calculation accuracy and the requirement of saving time, the numerical method with 65 × 49 grid points was introduced into this study. The grid number of 3D grids used in NUMECA is, respectively, 160k, 310k and 550k. The boundary conditions for the three grid topologies are the same. Moreover, the boundary conditions are the same as the one in the CAM method. Based on the results from these three grids, the distribution of static pressure on the blade surface and the total pressure at the outlet are shown in Figure 7. The curves of the case with 310k grid points coincide with that of the case with 550k grid points. Hence, in conclusion, when the number of grid points exceeds 310k, the results are considered independent of the grid number. In consideration of the calculation accuracy and the requirement of saving time, the numerical method with 310k grid points was involved in this study. The grid number of 3D grids used in NUMECA is, respectively, 160k, 310k and 550k. The boundary conditions for the three grid topologies are the same. Moreover, the boundary conditions are the same as the one in the CAM method. Based on the results from these three grids, the distribution of static pressure on the blade surface and the total pressure at the outlet are shown in Figure 7. The curves of the case with 310k grid points coincide with that of the case with 550k grid points. Hence, in conclusion, when the number of grid points exceeds 310k, the results are considered independent of the grid number. In consideration of the calculation accuracy and the requirement of saving time, the numerical method with 310k grid points was involved in this study.  Therefore, grids for these two methods are displayed in Figure 8. The grid settings in CAM are 65 and 49 in the axial and radial direction, respectively. In the 3D method, the grid type is H-O-H, and the total grid number is 310k. The CAM needs a CPU time of 100 ss on a laptop, while the 3D calculation requires 2 h on 4 cores of an Intel Xeon E5-2650 V2@2.606HZ processor.   Therefore, grids for these two methods are displayed in Figure 8. The grid settings in CAM are 65 and 49 in the axial and radial direction, respectively. In the 3D method, the grid type is H-O-H, and the total grid number is 310k. The CAM needs a CPU time of 100 ss on a laptop, while the 3D calculation requires 2 h on 4 cores of an Intel Xeon E5-2650 V2@2.606HZ processor.  Therefore, grids for these two methods are displayed in Figure 8. The grid settings in CAM are 65 and 49 in the axial and radial direction, respectively. In the 3D method, the grid type is H-O-H, and the total grid number is 310k. The CAM needs a CPU time of 100 ss on a laptop, while the 3D calculation requires 2 h on 4 cores of an Intel Xeon E5-2650 V2@2.606HZ processor.

Model Analysis on the Subsonic Condition
First, a comparison of the subsonic flow in linear cascade is performed, where mass

Model Analysis on the Subsonic Condition
First, a comparison of the subsonic flow in linear cascade is performed, where mass flow at the outlet is set to ensure the inlet Mach number of M = 0.77. Figures 9 and 10 show the isolines of the pressure and the axial velocity for both methods, and the values marked on them are the absolute values. In these figures, the leading edge (LE) and the trailing edge correspond to the beginning and end of the colorful isolines. One thing worth noting is that the circumferential average results from the 3D results are circumferentially mass average except the pressure, and for the pressure, it is the circumferential area average. Through these figures, it reveals that although parameters at the inlet and outlet are consistent, contours in the internal passage are different. These differences cannot be observed in the total pressure recovery coefficient and the spanwise distribution of parameters. In Figure 9b, the pressure obtained through the 3D method increases slowly, which shows the diffusive process in the vane passage. However, in the CAM, pressure first increases rapidly and then decreases slowly, which is contradicting to real flow. The velocity contours are shown in Figure 10 also reveal the difference between these two methods. The flow passes the LE should accelerate first and then decelerate with diffusion. However, in the result of CAM, the axial velocity decreases as soon as it flows into the blade passage. Figure 11a,b show the distribution of average pressure and axial velocity at the middle span, which also reveals the aforesaid phenomenon.    Figure 12a,b shows the distribution of the tangential and the axial components of the inviscid blade force at the middle span obtained in these two methods. Through comparison, it shows that the tangential component of the inviscid blade force has the same trend, which is solved through iteration. Since the tangential component is relatively accurate, the circulation at the inlet and outlet can be guaranteed, which makes the performance parameters from CAM similar to 3D results. However, the axial component has a great difference. Therefore, it is essential to depict the influence of the blade geometry on the inviscid blade force.  Figure 12a,b shows the distribution of the tangential and the axial components of the inviscid blade force at the middle span obtained in these two methods. Through comparison, it shows that the tangential component of the inviscid blade force has the same trend, which is solved through iteration. Since the tangential component is relatively accurate, the circulation at the inlet and outlet can be guaranteed, which makes the performance parameters from CAM similar to 3D results. However, the axial component has a great difference. Therefore, it is essential to depict the influence of the blade geometry on the inviscid blade force. In Formula (15), it is shown that pressures at the suction side and the pressure side are needed to depict the inviscid blade force modeling. However, pressure on the suction and pressure side cannot be acquired in the throughflow method; the term 2 s p p p + needs to be modeled. Figure 13 shows the pressure distribution in the middle span obtained from NUMECA. PP and PS refer to pressures on the pressure side and the suction side severally, and AP means the algebraic mean value of pressure at both sides, which is 2 s p p p + . SP means the spatially averaged value from 3D simulation results. From these figures, it can be observed that the trend of AP and SP is nearly the same, especially at the middle and posterior portion of the blade. Moreover, the average spatial pressure from the 3D numerical method has the same meaning as the average pressure proposed in the CAM method. Therefore, it is appropriate to use the average pressure in the CAM method in place of the term  In Formula (15), it is shown that pressures at the suction side and the pressure side are needed to depict the inviscid blade force modeling. However, pressure on the suction and pressure side cannot be acquired in the throughflow method; the term p s +p p 2 needs to be modeled. Figure 13 shows the pressure distribution in the middle span obtained from NUMECA. PP and PS refer to pressures on the pressure side and the suction side severally, and AP means the algebraic mean value of pressure at both sides, which is p s +p p 2 . SP means the spatially averaged value from 3D simulation results. From these figures, it can be observed that the trend of AP and SP is nearly the same, especially at the middle and posterior portion of the blade. Moreover, the average spatial pressure from the 3D numerical method has the same meaning as the average pressure proposed in the CAM method. Therefore, it is appropriate to use the average pressure in the CAM method in place of the term In Formula (15), it is shown that pressures at the suction side and the pressure side are needed to depict the inviscid blade force modeling. However, pressure on the suction and pressure side cannot be acquired in the throughflow method; the term 2 s p p p + needs to be modeled. Figure 13 shows the pressure distribution in the middle span obtained from NUMECA. PP and PS refer to pressures on the pressure side and the suction side severally, and AP means the algebraic mean value of pressure at both sides, which is 2 s p p p + . SP means the spatially averaged value from 3D simulation results. From these figures, it can be observed that the trend of AP and SP is nearly the same, especially at the middle and posterior portion of the blade. Moreover, the average spatial pressure from the 3D numerical method has the same meaning as the average pressure proposed in the CAM method. Therefore, it is appropriate to use the average pressure in the CAM method in place of the term Therefore, Formula (15) turns to be: Therefore, Formula (15) turns to be: The improved inviscid blade force model is applied in the CAM, where the boundary conditions remain the same as before. The axial and tangential components of the inviscid blade force obtained through these three methods are shown in Figure 14a,b ("REVISED" means the improved CAM). Through these figures, it can be seen that the axial component from the CAM improves greatly, which is nearly the same as the results from the 3D simulation. Moreover, the tangential component improves at the same time due to the influence of the axial blade force.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 13 of 25 means the improved CAM). Through these figures, it can be seen that the axial component from the CAM improves greatly, which is nearly the same as the results from the 3D simulation. Moreover, the tangential component improves at the same time due to the influence of the axial blade force.  Figure 15 displays the flow field inside the cascade obtained through the improved inviscid blade force modeling. Through these figures, it can be seen that the contours are similar to 3D simulation results. Figure 16a,b show the distribution of average pressure and average axial velocity at the middle span, which can reveal the effect resulting from the improved modeling. Based on the original inviscid blade force modeling, the maximum error of pressure along the streamwise direction between CAM and 3D method at the middle span is 30.9%. After using the improved model, the error decreases to 9.0%. As for the axial velocity, it drops from 21.0% to 7.5% after using this model. Moreover, these differences all appear at the leading edge of the blade. These errors can be explained as follows: this improved model is not suitable for all conditions. At the LE, pressure on the suction side and the pressure side is relevant to the incidence angle and changes drastically, which needs further modification due to the incidence angle.   Figure 15 displays the flow field inside the cascade obtained through the improved inviscid blade force modeling. Through these figures, it can be seen that the contours are similar to 3D simulation results. Figure 16a,b show the distribution of average pressure and average axial velocity at the middle span, which can reveal the effect resulting from the improved modeling. Based on the original inviscid blade force modeling, the maximum error of pressure along the streamwise direction between CAM and 3D method at the middle span is 30.9%. After using the improved model, the error decreases to 9.0%. As for the axial velocity, it drops from 21.0% to 7.5% after using this model. Moreover, these differences all appear at the leading edge of the blade. These errors can be explained as follows: this improved model is not suitable for all conditions. At the LE, pressure on the suction side and the pressure side is relevant to the incidence angle and changes drastically, which needs further modification due to the incidence angle.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 13 of 25 means the improved CAM). Through these figures, it can be seen that the axial component from the CAM improves greatly, which is nearly the same as the results from the 3D simulation. Moreover, the tangential component improves at the same time due to the influence of the axial blade force.  Figure 15 displays the flow field inside the cascade obtained through the improved inviscid blade force modeling. Through these figures, it can be seen that the contours are similar to 3D simulation results. Figure 16a,b show the distribution of average pressure and average axial velocity at the middle span, which can reveal the effect resulting from the improved modeling. Based on the original inviscid blade force modeling, the maximum error of pressure along the streamwise direction between CAM and 3D method at the middle span is 30.9%. After using the improved model, the error decreases to 9.0%. As for the axial velocity, it drops from 21.0% to 7.5% after using this model. Moreover, these differences all appear at the leading edge of the blade. These errors can be explained as follows: this improved model is not suitable for all conditions. At the LE, pressure on the suction side and the pressure side is relevant to the incidence angle and changes drastically, which needs further modification due to the incidence angle.

Model Analysis on the Transonic Condition
In previous research, when the CAM is applied in the transonic condition, a special shock wave was captured [13,17,20]. However, in real 3D transonic flow, the shock wave is normal to the blade geometry other than axisymmetric. Therefore, in the CAM, the shock wave should be eliminated, where a decelerating process exists. Considering this point, the inviscid blade force involved in the transonic cascade also needs to be improved.
First, a simulation of the transonic cascade is performed through the 3D method. By changing pressure at the outlet, the Mach number at the inlet is restricted to 0.86. The isoline of Mach numbers at the middle span is shown in Figure 17, which reveals that the strong shock wave exists on the suction side. Similar to Figure 13, the distribution of pressure at the middle span is also displayed in Figure 18 for this condition. In this figure, it can be seen that the spatially averaged pressure from the 3D result is nearly the same as the algebraic mean value of the pressure side and the suction side. They have slight differences at the LE and the front half part, and its value at the shock wave of the suction side is a little different. However, they are coincident behind the shock wave. To examine the influence of these two aforesaid inviscid blade force models on the result, they are both applied in the CAM simulation. The distribution of the pressure and

Model Analysis on the Transonic Condition
In previous research, when the CAM is applied in the transonic condition, a special shock wave was captured [13,17,20]. However, in real 3D transonic flow, the shock wave is normal to the blade geometry other than axisymmetric. Therefore, in the CAM, the shock wave should be eliminated, where a decelerating process exists. Considering this point, the inviscid blade force involved in the transonic cascade also needs to be improved.
First, a simulation of the transonic cascade is performed through the 3D method. By changing pressure at the outlet, the Mach number at the inlet is restricted to 0.86. The isoline of Mach numbers at the middle span is shown in Figure 17, which reveals that the strong shock wave exists on the suction side.

Model Analysis on the Transonic Condition
In previous research, when the CAM is applied in the transonic condition, a special shock wave was captured [13,17,20]. However, in real 3D transonic flow, the shock wave is normal to the blade geometry other than axisymmetric. Therefore, in the CAM, the shock wave should be eliminated, where a decelerating process exists. Considering this point, the inviscid blade force involved in the transonic cascade also needs to be improved.
First, a simulation of the transonic cascade is performed through the 3D method. By changing pressure at the outlet, the Mach number at the inlet is restricted to 0.86. The isoline of Mach numbers at the middle span is shown in Figure 17, which reveals that the strong shock wave exists on the suction side. Similar to Figure 13, the distribution of pressure at the middle span is also displayed in Figure 18 for this condition. In this figure, it can be seen that the spatially averaged pressure from the 3D result is nearly the same as the algebraic mean value of the pressure side and the suction side. They have slight differences at the LE and the front half part, and its value at the shock wave of the suction side is a little different. However, they are coincident behind the shock wave. To examine the influence of these two aforesaid inviscid blade force models on the result, they are both applied in the CAM simulation. The distribution of the pressure and Similar to Figure 13, the distribution of pressure at the middle span is also displayed in Figure 18 for this condition. In this figure, it can be seen that the spatially averaged pressure from the 3D result is nearly the same as the algebraic mean value of the pressure side and the suction side. They have slight differences at the LE and the front half part, and its value at the shock wave of the suction side is a little different. However, they are coincident behind the shock wave.

Model Analysis on the Transonic Condition
In previous research, when the CAM is applied in the transonic condition, a special shock wave was captured [13,17,20]. However, in real 3D transonic flow, the shock wave is normal to the blade geometry other than axisymmetric. Therefore, in the CAM, the shock wave should be eliminated, where a decelerating process exists. Considering this point, the inviscid blade force involved in the transonic cascade also needs to be improved.
First, a simulation of the transonic cascade is performed through the 3D method. By changing pressure at the outlet, the Mach number at the inlet is restricted to 0.86. The isoline of Mach numbers at the middle span is shown in Figure 17, which reveals that the strong shock wave exists on the suction side. Similar to Figure 13, the distribution of pressure at the middle span is also displayed in Figure 18 for this condition. In this figure, it can be seen that the spatially averaged pressure from the 3D result is nearly the same as the algebraic mean value of the pressure side and the suction side. They have slight differences at the LE and the front half part, and its value at the shock wave of the suction side is a little different. However, they are coincident behind the shock wave. To examine the influence of these two aforesaid inviscid blade force models on the result, they are both applied in the CAM simulation. The distribution of the pressure and To examine the influence of these two aforesaid inviscid blade force models on the result, they are both applied in the CAM simulation. The distribution of the pressure and the axial velocity at the midspan are compared in Figure 19. It reveals a huge difference between these two models. When the original inviscid blade force is used, the maximum error between CAM and 3D method for pressure along the flow direction is 128.1%, and it is 47.0% for axial velocity. However, when the improved model is used, these two values are reduced to 39.3% and 19.9%, respectively. The most obvious diversity exists at the position where shock wave occurs, and at other positions, these two sets of results obtained from the revised model and NUMECA are nearly the same. Finally, one element needs to be discussed. The maximum error of the static pressure is still 39.3%, which is still large. It is observed that the maximum error exists at the location of the shock. Considering the difficulties in the simulation of the shock with the throughflow method [11,13,[20][21][22] and the improvement of the parameters in the blade passage, the method is regarded as an applicable method. To further improve the accuracy, there are two aspects that may need to consider. The first one is to add a factor in adjusting variable incidence angles, which makes this model more exact for different incidence angles. Moreover, the other one is that the circumferential fluctuation source terms [23,24] should also be considered in the governing equation since they also contribute to the flow structure. the axial velocity at the midspan are compared in Figure 19. It reveals a huge difference between these two models. When the original inviscid blade force is used, the maximum error between CAM and 3D method for pressure along the flow direction is 128.1%, and it is 47.0% for axial velocity. However, when the improved model is used, these two values are reduced to 39.3% and 19.9%, respectively. The most obvious diversity exists at the position where shock wave occurs, and at other positions, these two sets of results obtained from the revised model and NUMECA are nearly the same. Finally, one element needs to be discussed. The maximum error of the static pressure is still 39.3%, which is still large. It is observed that the maximum error exists at the location of the shock. Considering the difficulties in the simulation of the shock with the throughflow method [11,13,[20][21][22] and the improvement of the parameters in the blade passage, the method is regarded as an applicable method. To further improve the accuracy, there are two aspects that may need to consider. The first one is to add a factor in adjusting variable incidence angles, which makes this model more exact for different incidence angles. Moreover, the other one is that the circumferential fluctuation source terms [23,24] should also be considered in the governing equation since they also contribute to the flow structure.

Experimental Facility
After the improved model of the inviscid blade, force is verified in the linear cascade; it is then applied in a fan named TA36 to validate its applicability. The TA36 is a highly loaded and low aspect ratio fan, which is comprised of a rotor and a stator. Moreover, the load coefficient is up to 0.258. The designed rotational speed is 2900 rpm. The flow path of the TA36 is displayed in Figure 20. The geometrical parameters of the TA36 are shown in Table 2.

Experimental Facility
After the improved model of the inviscid blade, force is verified in the linear cascade; it is then applied in a fan named TA36 to validate its applicability. The TA36 is a highly loaded and low aspect ratio fan, which is comprised of a rotor and a stator. Moreover, the load coefficient is up to 0.258. The designed rotational speed is 2900 rpm. The flow path of the TA36 is displayed in Figure 20. The geometrical parameters of the TA36 are shown in Table 2.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 15 of 25 the axial velocity at the midspan are compared in Figure 19. It reveals a huge difference between these two models. When the original inviscid blade force is used, the maximum error between CAM and 3D method for pressure along the flow direction is 128.1%, and it is 47.0% for axial velocity. However, when the improved model is used, these two values are reduced to 39.3% and 19.9%, respectively. The most obvious diversity exists at the position where shock wave occurs, and at other positions, these two sets of results obtained from the revised model and NUMECA are nearly the same. Finally, one element needs to be discussed. The maximum error of the static pressure is still 39.3%, which is still large. It is observed that the maximum error exists at the location of the shock. Considering the difficulties in the simulation of the shock with the throughflow method [11,13,[20][21][22] and the improvement of the parameters in the blade passage, the method is regarded as an applicable method. To further improve the accuracy, there are two aspects that may need to consider. The first one is to add a factor in adjusting variable incidence angles, which makes this model more exact for different incidence angles. Moreover, the other one is that the circumferential fluctuation source terms [23,24] should also be considered in the governing equation since they also contribute to the flow structure.

Experimental Facility
After the improved model of the inviscid blade, force is verified in the linear cascade; it is then applied in a fan named TA36 to validate its applicability. The TA36 is a highly loaded and low aspect ratio fan, which is comprised of a rotor and a stator. Moreover, the load coefficient is up to 0.258. The designed rotational speed is 2900 rpm. The flow path of the TA36 is displayed in Figure 20. The geometrical parameters of the TA36 are shown in Table 2.    Figure 22 are as follows: 1-inlet flow tube, 2-measurement point of wall static pressure, 3-measurement point of inlet total pressure, 4-rotor, 5-stator, 6measurement point of outlet total and static pressure, 7-motor, 8-brackets, 9-outlet adjusting components, 10-throttling valve. The rotor is motivated through an actuating motor, and the rotational speed of the rotor is controlled by a frequency transmitter. The working point of the fan is commanded by the opening degree of the throttle valve. The throttle valve consists of a fixed conical surface and a removable annular sleeve, which can modify the opening degree of the throttle valve by adjusting the relative position of these two components. The position of the annular sleeve is controlled by a stepping motor, which can make the fan work from the maximum mass flow point to the near stall point.  Figure 22 are as follows: 1-inlet flow tube, 2-measurement point of wall static pressure, 3-measurement point of inlet total pressure, 4-rotor, 5-stator, 6-measurement point of outlet total and static pressure, 7-motor, 8-brackets, 9-outlet adjusting components, 10-throttling valve. The rotor is motivated through an actuating motor, and the rotational speed of the rotor is controlled by a frequency transmitter. The working point of the fan is commanded by the opening degree of the throttle valve. The throttle valve consists of a fixed conical surface and a removable annular sleeve, which can modify the opening degree of the throttle valve by adjusting the relative position of these two components. The position of the annular sleeve is controlled by a stepping motor, which can make the fan work from the maximum mass flow point to the near stall point.
In this test facility, two planes named 3 and 6 are chosen to measure the total pressure at the inlet and the outlet of the fan. The total pressure is measured through the five-hole total pressure comb. The test instruments and their distribution are shown in Figure 23. The static wall pressure is measured at position 2 through four static pressure measuring holes, which are uniformly distributed in the circumferential direction. The static pressure sensor involved is Honeywell-0142PC01D. The steady-state pressure signal is delivered to the data acquisition board belonging to the NI PXI system through the terminal box. Through this method, the pressure parameters of the fan system are obtained. The torque of the actuating motor in the test is measured through the torque meter. With this approach, the rotational speed and torque can be obtained from the data acquisition system.
Before the measurement, the five-hole total pressure comb is calibrated. In order to guarantee the accuracy, repeated tests are conducted. The uncertainties of the mass flow rate and efficiency in this experiment are within 0.5%.   In this test facility, two planes named 3 and 6 are chosen to measure the total pressure at the inlet and the outlet of the fan. The total pressure is measured through the five-hole total pressure comb. The test instruments and their distribution are shown in Figure 23. The static wall pressure is measured at position 2 through four static pressure measuring holes, which are uniformly distributed in the circumferential direction. The static pressure sensor involved is Honeywell-0142PC01D. The steady-state pressure signal is delivered to the data acquisition board belonging to the NI PXI system through the terminal box. Through this method, the pressure parameters of the fan system are obtained. The torque of the actuating motor in the test is measured through the torque meter. With this approach, the rotational speed and torque can be obtained from the data acquisition system. Next, the process of this experiment is described in detail. First, the area of the throttle valve is adjusted to the maximum, and the fan is accelerated to the set rotational speed. Then the opening area of the throttle valve is decreased slowly to make sure that the fan works at a quasi-steady state. With the fan working at the stall state, the throttle valve is opened to quit this stall state. In this experiment, the measurement parameters are the mass flow, total pressure ratio and efficiency. The volume flow rate of the fan is calculated as:  Therefore, the efficiency is calculated as: Before the measurement, the five-hole total pressure comb is calibrated. In order to guarantee the accuracy, repeated tests are conducted. The uncertainties of the mass flow rate and efficiency in this experiment are within 0.5%.
Next, the process of this experiment is described in detail. First, the area of the throttle valve is adjusted to the maximum, and the fan is accelerated to the set rotational speed. Then the opening area of the throttle valve is decreased slowly to make sure that the fan works at a quasi-steady state. With the fan working at the stall state, the throttle valve is opened to quit this stall state. In this experiment, the measurement parameters are the mass flow, total pressure ratio and efficiency. The volume flow rate of the fan is calculated as: where α is the discharge coefficient of the fan, and it is 0.992 in the test, A 0 is the area of the inlet, ρ 0 is the density of the atmosphere, p sF is the pressure increase of the fan. The efficiency formula is: where P u is the useful power of the fan, and P a is the shaft power. They are defined as: where U and I are the voltage and the current of the motor, cos ϕ is the power coefficient, whose range is 0.75-0.85.
where p t1 and p t2 are the total pressure at the inlet and outlet severally. Therefore, the efficiency is calculated as:

Grid Independence Verification
To verify the usability of the improved inviscid blade force model, the CAM based on N-S equations and the 3D method are both applied in the fan simulation. The computational grid topology used in both methods is shown in Figure 24.
The boundary condition for both simulations is set as follows: The total temperature and pressure are 288.15 K and 101,325 Pa severally. The inlet flow direction is axial. The pressure is restricted at the hub of the outlet, and the pressures at other positions of the outlet are obtained through a radial-equilibrium equation. At the solid wall, the no-slip boundary condition is used. The turbulence model is the S-A model.

Grid Independence Verification
To verify the usability of the improved inviscid blade force model, the CAM based on N-S equations and the 3D method are both applied in the fan simulation. The computational grid topology used in both methods is shown in Figure 24. The boundary condition for both simulations is set as follows: The total temperature and pressure are 288.15 K and 101,325 Pa severally. The inlet flow direction is axial. The pressure is restricted at the hub of the outlet, and the pressures at other positions of the outlet are obtained through a radial-equilibrium equation. At the solid wall, the no-slip boundary condition is used. The turbulence model is the S-A model.
To check the grid independence in the simulation of the TA36, three grids are made in both the CAM and NUMECA. For the CAM, the grid number is, respectively 36 × 33 + 34 × 33, 53 × 49 + 49 × 49, 59 × 55 + 55 × 55. The performance characteristics of this fan based on these three grids are obtained, shown in Figure 25. They are the pressure ratio versus the mass flow and the adiabatic efficiency versus the mass flow. The curves of the case with 53 × 49 + 49 × 49 grid points coincide with that of the case with 59 × 55 + 55 × 55 grid points. Hence, in conclusion, when the number of grid points exceeds 53 × 49 + 49 × 49, the results are considered independent of the grid number. In consideration of the calculation accuracy and the requirement of saving time, the numerical method with 53 × 49 + 49 × 49 grid points were involved in this study. The grid number used in the NUMECA for the TA36 are, respectively, 450k, 610k, 830k. The performance characteristics of this fan based on these three grids are obtained, shown in Figure 26. They are the pressure ratio versus the mass flow and the adiabatic To check the grid independence in the simulation of the TA36, three grids are made in both the CAM and NUMECA. For the CAM, the grid number is, respectively 36 × 33 + 34 × 33, 53 × 49 + 49 × 49, 59 × 55 + 55 × 55. The performance characteristics of this fan based on these three grids are obtained, shown in Figure 25. They are the pressure ratio versus the mass flow and the adiabatic efficiency versus the mass flow. The curves of the case with 53 × 49 + 49 × 49 grid points coincide with that of the case with 59 × 55 + 55 × 55 grid points. Hence, in conclusion, when the number of grid points exceeds 53 × 49 + 49 × 49, the results are considered independent of the grid number. In consideration of the calculation accuracy and the requirement of saving time, the numerical method with 53 × 49 + 49 × 49 grid points were involved in this study.
To verify the usability of the improved inviscid blade force model, the CAM based on N-S equations and the 3D method are both applied in the fan simulation. The computational grid topology used in both methods is shown in Figure 24.  The grid number used in the NUMECA for the TA36 are, respectively, 450k, 610k, 830k. The performance characteristics of this fan based on these three grids are obtained, shown in Figure 26. They are the pressure ratio versus the mass flow and the adiabatic The grid number used in the NUMECA for the TA36 are, respectively, 450k, 610k, 830k. The performance characteristics of this fan based on these three grids are obtained, shown in Figure 26. They are the pressure ratio versus the mass flow and the adiabatic efficiency versus the mass flow. The curves of the case with 610k grid points coincide with that of the case with 830k grid points. Hence, in conclusion, when the number of grid points exceeds 610k, the results are considered independent of the grid number. In consideration of the calculation accuracy and the requirement of saving time, the numerical method with 610k grid points was introduced into this study.
Therefore, the mesh for the throughflow method is made of 53 × 49 (Rotor) + 49 × 49 (Stator) nodes. The computational region for the 3D method meshes with nearly 610,000 grids, and a tip clearance of 1 mm exists at the rotor domain. efficiency versus the mass flow. The curves of the case with 610k grid points coincide with that of the case with 830k grid points. Hence, in conclusion, when the number of grid points exceeds 610k, the results are considered independent of the grid number. In consideration of the calculation accuracy and the requirement of saving time, the numerical method with 610k grid points was introduced into this study. Therefore, the mesh for the throughflow method is made of 53 × 49 (Rotor) + 49 × 49 (Stator) nodes. The computational region for the 3D method meshes with nearly 610,000 grids, and a tip clearance of 1 mm exists at the rotor domain. Figure 27a,b shows the overall characteristic lines acquired from the experiment, NU-MECA and these two CAM methods. Figure 27a is shown for the characteristic of mass flow versus total pressure ratio. Figure 27b is for mass flow versus adiabatic efficiency. From these two figures, it is observed that the performance of efficiency is improved, and the characteristic of mass flow and pressure ratio remains unchanged. The performance parameters on the design point acquired from different calculation methods are shown in Table 3. When the original inviscid blade force model is applied, mass flow and pressure ratio from CAM, experiment and 3D method is nearly the same. However, the error of adiabatic efficiency between them is more than 4.0%. This is related to the simulation accuracy of the flow field. The original model is based on the assumption that the inviscid blade force is perpendicular to the average stream surface, which leads to the wrong prediction of the field. While in fact, the inviscid blade force is not. After the improved inviscid blade force model is applied, the error of adiabatic efficiency between throughflow and experiment or 3D simulation is about 1%.  Figure 27a,b shows the overall characteristic lines acquired from the experiment, NUMECA and these two CAM methods. Figure 27a is shown for the characteristic of mass flow versus total pressure ratio. Figure 27b is for mass flow versus adiabatic efficiency. From these two figures, it is observed that the performance of efficiency is improved, and the characteristic of mass flow and pressure ratio remains unchanged.

Numerical Analysis and Discussion
points exceeds 610k, the results are considered independent of the grid number. In consideration of the calculation accuracy and the requirement of saving time, the numerical method with 610k grid points was introduced into this study. Therefore, the mesh for the throughflow method is made of 53 × 49 (Rotor) + 49 × 49 (Stator) nodes. The computational region for the 3D method meshes with nearly 610,000 grids, and a tip clearance of 1 mm exists at the rotor domain. Figure 27a,b shows the overall characteristic lines acquired from the experiment, NU-MECA and these two CAM methods. Figure 27a is shown for the characteristic of mass flow versus total pressure ratio. Figure 27b is for mass flow versus adiabatic efficiency. From these two figures, it is observed that the performance of efficiency is improved, and the characteristic of mass flow and pressure ratio remains unchanged. The performance parameters on the design point acquired from different calculation methods are shown in Table 3. When the original inviscid blade force model is applied, mass flow and pressure ratio from CAM, experiment and 3D method is nearly the same. However, the error of adiabatic efficiency between them is more than 4.0%. This is related to the simulation accuracy of the flow field. The original model is based on the assumption that the inviscid blade force is perpendicular to the average stream surface, which leads to the wrong prediction of the field. While in fact, the inviscid blade force is not. After the improved inviscid blade force model is applied, the error of adiabatic efficiency between throughflow and experiment or 3D simulation is about 1%. The performance parameters on the design point acquired from different calculation methods are shown in Table 3. When the original inviscid blade force model is applied, mass flow and pressure ratio from CAM, experiment and 3D method is nearly the same. However, the error of adiabatic efficiency between them is more than 4.0%. This is related to the simulation accuracy of the flow field. The original model is based on the assumption that the inviscid blade force is perpendicular to the average stream surface, which leads to the wrong prediction of the field. While in fact, the inviscid blade force is not. After the improved inviscid blade force model is applied, the error of adiabatic efficiency between throughflow and experiment or 3D simulation is about 1%.  Figure 28a,c shows the distribution of rotor pressure ratio and adiabatic efficiency, stator total pressure recovery coefficient, respectively. It reveals that the trend of adiabatic efficiency obtained from the improved CAM is much closer to the 3D simulation results, and the trend of rotor pressure ratio and stator total pressure recovery coefficient remains nearly consistent.  Figure 28a,c shows the distribution of rotor pressure ratio and adiabatic efficiency, stator total pressure recovery coefficient, respectively. It reveals that the trend of adiabatic efficiency obtained from the improved CAM is much closer to the 3D simulation results, and the trend of rotor pressure ratio and stator total pressure recovery coefficient remains nearly consistent. Figures 29-31 show the components of the inviscid blade force in the axial, radial, and tangential direction obtained both in the revised CAM and 3D method. It is observed that after using the improved blade force model, the calculation precision is promoted obviously. The axial and radial components are coincident with the results from the 3D simulation. On the other hand, the accuracy of the tangential component is also improved.       Figures 32 and 33 show the isolines of the axial velocity and the pressure obtained through two methods. From these pictures, it can be seen that these two isolines are nearly the same. The only difference exists at the hub of the rotor, where the flow accelerates. The phenomenon also can be seen at the isolines of pressure. These differences may arise from the reason that the inviscid blade force model still needs to be further modified to adapt to the variable incidences.   Figures 32 and 33 show the isolines of the axial velocity and the pressure obtained through two methods. From these pictures, it can be seen that these two isolines are nearly the same. The only difference exists at the hub of the rotor, where the flow accelerates. The phenomenon also can be seen at the isolines of pressure. These differences may arise from the reason that the inviscid blade force model still needs to be further modified to adapt to the variable incidences.  Figures 32 and 33 show the isolines of the axial velocity and the pressure obtained through two methods. From these pictures, it can be seen that these two isolines are nearly the same. The only difference exists at the hub of the rotor, where the flow accelerates. The phenomenon also can be seen at the isolines of pressure. These differences may arise from the reason that the inviscid blade force model still needs to be further modified to adapt to the variable incidences.  To make the comparison between the results from 2D and 3D much clearer, the distributions of these parameters along a grid line in the streamwise direction are displayed. The axial component of the inviscid blade force and the axial velocity along the gridline are analyzed as two examples. Moreover, the gridline we choose is shown in Figure 34, which is red. The aerodynamic parameters along this line are shown in Figure 35. It is observed that the axial component of the inviscid blade force is almost the same from these two methods. Moreover, for the axial velocity, the trend of these two results is the same, and the biggest error is nearly 4%. To make the comparison between the results from 2D and 3D much clearer, the distributions of these parameters along a grid line in the streamwise direction are displayed. The axial component of the inviscid blade force and the axial velocity along the gridline are analyzed as two examples. Moreover, the gridline we choose is shown in Figure 34, which is red. The aerodynamic parameters along this line are shown in Figure 35. It is observed that the axial component of the inviscid blade force is almost the same from these two methods. Moreover, for the axial velocity, the trend of these two results is the same, and the biggest error is nearly 4%. tributions of these parameters along a grid line in the streamwise direction are displayed. The axial component of the inviscid blade force and the axial velocity along the gridline are analyzed as two examples. Moreover, the gridline we choose is shown in Figure 34, which is red. The aerodynamic parameters along this line are shown in Figure 35. It is observed that the axial component of the inviscid blade force is almost the same from these two methods. Moreover, for the axial velocity, the trend of these two results is the same, and the biggest error is nearly 4%.

Conclusions
In the design process of the traditional gas turbine engines, vast experiment tests are needed. To reduce unnecessary trials, a whole engine simulation tool is extremely needed in the preliminary design period. Considering the drawbacks of 0D analysis and 3D CFD calculation, the 2D throughflow analysis method is still an indispensable tool. For this The axial component of the inviscid blade force and the axial velocity along the gridline are analyzed as two examples. Moreover, the gridline we choose is shown in Figure 34, which is red. The aerodynamic parameters along this line are shown in Figure 35. It is observed that the axial component of the inviscid blade force is almost the same from these two methods. Moreover, for the axial velocity, the trend of these two results is the same, and the biggest error is nearly 4%.

Conclusions
In the design process of the traditional gas turbine engines, vast experiment tests are needed. To reduce unnecessary trials, a whole engine simulation tool is extremely needed in the preliminary design period. Considering the drawbacks of 0D analysis and 3D CFD calculation, the 2D throughflow analysis method is still an indispensable tool. For this

Conclusions
In the design process of the traditional gas turbine engines, vast experiment tests are needed. To reduce unnecessary trials, a whole engine simulation tool is extremely needed in the preliminary design period. Considering the drawbacks of 0D analysis and 3D CFD calculation, the 2D throughflow analysis method is still an indispensable tool. For this purpose, a compressor simulation tool is first developed. Based on the average circumferential method, 3D N-S and Euler equations are transformed into 2D models. When they are involved in the simulation of the compressor, an improved inviscid blade force model considering the effect of the blade geometry is derived.
(1) The inviscid blade force is attributed to the tangential difference of pressure in the blade passage and the gradient of the blade blockage factor. The former is to deflect the flow direction, and the latter is to reduce the flow path area. However, in previous research, the assumption is that gas flows along the average stream surface of the blade, and the effect of the blade geometry is not taken into consideration, which is not applicable in the CAM.
(2) The improved inviscid blade force model is obtained from the analysis in a linear cascade on the subsonic condition. Through analyzing the results from the 3D calculation, it is observed that the spatially averaged pressure is nearly the same as the algebraic mean value of pressure on the pressure and suction side. Since the spatially average pressure has the same meaning as the pressure in CAM, it is thought to be reasonable to replace p s +p p 2 with the pressure value obtained in CAM.
(3) When compared with the 3D results, the maximum error of pressure along the streamwise direction decreases from 30.9% to 9.0%, based on the improved model on the subsonic condition. For the axial velocity, it drops from 21.0% to 7.5%. When the model is applied on the transonic condition, the maximum pressure error between CAM and 3D results in changes from 128.1% to 39.3%, and for axial velocity, the error is reduced