KR20120013501A - Optimization design method for casing grooves of an axial compressor and the optimized casing grooves - Google Patents

Abstract

PURPOSE: An optimization design method for a casing groove of an axial compressor is provided to optimize the shape of a casing groove using an optimization theory. CONSTITUTION: A design variable of a casing groove is selected in consideration of the shape of an axial compressor. An objective function is determined(S10). A design area is set while the design variable is limited(S20). Numerical analysis is performed in the selected design area(S30). An optimization point is searched in the design area based on the numerical analysis result(S40).

Description

Optimal design method for casing grooves of an axial compressor and the optimized casing grooves}

The present invention relates to a casing groove optimum design method and optimally designed casing groove of the axial compressor, and more specifically, to the casing groove through the numerical definition of the shape of the axial compressor, lattice generation, boundary conditions and internal flow field Performance of the grooved axial compressor and the basic grooved axial compressor and the grooved axial compressor optimized by the present technology to verify the maximum stall margin in the optimized grooved axial compressor. The present invention relates to a casing groove optimum design method and optimally designed casing groove of an axial compressor, which enables the design of a casing groove to continuously improve the stall margin.

Axial compressors are turbomachines operating in a high efficiency and high flow rate area, and are mainly used in gas turbines such as jet engines and marine engines. When the axial compressor reaches the low flow rate region, severe vibration occurs due to stall and surge phenomena, which causes many problems not only in performance such as efficiency but also in the operational stability of the compressor. The design should take account of the conditions.

Recently, through the three-dimensional numerical analysis based on Computational Fluid Dynamics (CFD), it is possible to predict the flow mechanism such as stall and surge occurring in the low flow region of the axial compressor. Actively done.

As a result of the study through the numerical analysis, the casing groove is formed by digging a groove on the case surface of the axial compressor to stall stall by delaying stall phenomenon caused by the interaction of tip leakage vortex and shock wave occurring near the case of the axial compressor. It is known that margins can be improved and surge margins can be increased.

With this technique, U.S. Patent Publication No. US2003 / 0138317 discloses a configuration of casing treatment for a fluid compressor for forming one or more grooves to prevent the growth of potentially unstable vortex near the tip of the compressor wing. The above configuration merely describes the formation of grooves in the case, but does not appear at all in the configuration for optimizing the casing groove by limiting the shape, size and number of grooves.

In addition, in recent years, studies have been conducted to numerically calculate the shape of the casing groove by trial and error by various methods, but it is possible to design a casing groove for continuously improving the stall margin through an optimal design technique. The method is not yet proposed.

The present invention has been made to solve the problems described above, the object of the present invention is to optimize the shape of the casing groove in the entire design area through a systematic optimization design using optimization theory rather than a limited test through trial and error The present invention provides an optimal casing groove design method and an optimal casing groove for an axial compressor.

According to an aspect of the present invention,

Formulating a problem of selecting a design variable of the casing groove in consideration of the shape of the axial compressor and determining an objective function; A design area selection step of setting a design area by limiting a range of design variables; A numerical analysis step of performing numerical analysis in the selected design area; And an optimum point searching step of searching for an optimum point in the design area based on the numerical analysis result.

At this time, the design variable is characterized in that the width and width of the casing groove dimensioned by the length of the tip cord.

으로 정의되는 축류압축기의 스톨마진인 것을 특징으로 한다.In addition, the objective function is

It is characterized in that the stall margin of the axial compressor defined by.

: 최대효율점에서의 유량,

:스톨인근점에서의 유량,

: 최대효율점에서의 전압력,

: 스톨인근점에서의 유량.)(At this time, SM: Stall margin,

= Flow rate at the maximum efficiency point,

Flow rate near the stall

= Voltage force at the maximum efficiency point,

: Flow rate near stole.)

In the search for the optimum point, the optimal design parameter value of the casing groove may be determined from a neural network technique.

On the other hand, the optimally designed casing groove according to the present invention,

In the casing groove of the axial compressor, the depth and width of the casing groove dimensioned by the length of the blade tip chord are 0.177 and 0.164, respectively.

According to the present invention, the neural network model is used to optimize the shape of the casing groove, and compare the performance of the axial compressor without grooves, the axial compressor with basic grooves, and the axial compressor with grooves optimized by the present technology. By verifying the maximum stall margin in axial compressors with optimized grooves, the casing groove can be designed to continuously improve the stall margin.

1 is a flow chart showing a casing groove optimum design method of the axial compressor according to the present invention.
2 schematically shows a casing groove of an axial compressor;
Figure 3 (a), (b) is a view showing the three-dimensional shape and the meridian shape of NASA Rotor 37 used in the present invention.
Figure 4 (a), (b) is a diagram showing a hexagonal lattice system of the single channel used in the numerical analysis step of the present invention shown in Figure 1, and a hexagonal lattice system of the groove.
5 is a view showing a comparison with the existing performance test to verify the analysis results in the numerical analysis step of the present invention shown in FIG.
6 is a view showing a three-dimensional search results by the neural network model in the optimum point search step of the present invention shown in FIG.
7 shows a comparison of an optimized casing groove and a reference shape casing groove according to the invention.
8 is a diagram illustrating a sensitivity analysis of design variables of an objective function.
FIG. 9 shows a comparison of compressor performance with grooves of reference and optimum shape grooves. FIG.
(A)-(c) is a figure which shows the grooveless shape, the reference shape, the optimal shape voltage force, and tip leakage vortex at the stall vicinity point of a shape without a casing groove.
(A)-(c) is a figure which shows the Mach number with respect to the shape without a groove, a reference shape, and an optimal shape in the vicinity of the stall of a shape without a groove.
12 (a) and 12 (b) are diagrams showing Mach number distributions of a reference shape and an optimal shape at a stall neighborhood of the reference shape.
13 (a) and 13 (b) show velocity vectors of a reference shape and an optimal shape in the AB section in FIGS. 12 (a) and 12 (b).

Hereinafter, with reference to the accompanying drawings will be described in detail a preferred embodiment of the casing groove optimal design method and optimally designed casing groove of the axial compressor according to the present invention.

1 is a flow chart showing a casing groove optimum design method of the axial compressor according to the present invention, Figure 2 is a schematic view showing a casing groove of the axial compressor, Figure 3 (a), (b) is used in the present invention Figure 3 shows a three-dimensional shape and a meridional shape of NASA Rotor 37, Figure 4 (a), (b) is a hexagonal lattice system of a single channel used in the numerical analysis step of the present invention shown in Figure 1, the groove FIG. 5 is a diagram showing a hexahedral lattice system, and FIG. 5 is a diagram showing a comparison result with an existing performance test to verify an analysis result in a numerical analysis step of the present invention illustrated in FIG. 1, and FIG. Figure 3 is a view showing a three-dimensional search results by the neural network model in the optimum point search step, Figure 7 is a view showing a comparison of the casing groove of the optimized casing groove and the reference shape according to the present invention, 8 is a diagram illustrating the sensitivity analysis of the design variables of the objective function. FIG. 9 is a diagram showing the performance of a compressor having a grooveless shape and a groove having a reference shape and an optimal shape, and FIGS. c) is a diagram showing the grooveless shape, the reference shape, the optimal shape of the voltage force and the tip leakage vortex at the vicinity of the stall of the shape without the casing groove, and (a) to (c) of FIG. Fig. 12 (a) and 12 (b) show the Mach number distribution of the reference shape and the optimal shape at the point near the stall of the reference shape. 13 (a) and 13 (b) are diagrams showing velocity vectors of the reference shape and the optimum shape in the AB section in FIGS. 12 (a) and 12 (b).

The present invention optimizes the shape of the casing groove by numerically defining the shape, lattice generation, boundary condition, and internal flow field of the axial compressor, and the axial compressor without grooves and the axial compressor with basic grooves and the present technology Comparison of the performance of axial compressors with grooves optimized by means of verifying that the stall margin is maximized in the axial compressors with optimized grooves, so that casing grooves can be designed to continuously improve the stall margin. The groove optimal design method and the optimally designed casing groove. First, the casing groove optimal design method (hereinafter referred to as 'optimal design method') of the axial compressor according to the present invention is the formulation of the problem (S10), the design area selection step. (S20), the numerical analysis step (S30) and the optimum point search step (S40) It is open configuration.

In more detail, the formulating step (S10) of the problem is to select the design variables of the casing groove 30 to be formed in the case in consideration of the shape of the axial compressor, so that the maximum within the design area, that is, the purpose of optimization It relates to the step of determining an objective function.

That is, a design variable is first selected to determine the shape of the casing groove 30 for maximizing the objective function. In the present invention, as shown in FIG. The width (W / C) and depth (D / C) of the original casing groove 30 were selected as design variables.

In addition, since the present invention aims at optimizing the shape of the casing groove 30 to maximize the stall margin, the Stall Margin (SM) defined as follows is used as the objective function.

... (1)

... (One)

Where m and PR represent the flow and voltage forces, respectively, and the subscript peak and stall represent the maximum efficiency point and the stall near point.

The flow rate (stall flow rate) at which the stall occurs in the low flow rate area near the stall was determined as the point where the numerical calculation finally converged when the flow rate was decreased.The flow rate change at the inlet was 0.001kg / s per 300 steps, Convergence was determined when the flow rate difference at the outlet was 0.5% or less and the change in the adiabatic efficiency was 0.03% or less per 100 steps.

Next, the design region selection step (S20) is to set the appropriate design region by limiting the range of design variables for performing the optimal design, and to configure the optimal grid system for analysis in the numerical analysis step (s30) to be described later In relation to the present invention, a design variable, that is, a range of the width (W / C) and the depth (D / C) of the dimensionless casing groove 30 is shown in Table 1 below.

Variables Lower bound Upper bound D / C 0.077 0.231 W / C 0.120 0.192

That is, the depth of the dimensionless casing groove 30 is limited to 0.077 to 0.231, and the width of the dimensionless casing groove 30 is limited to 0.120 to 0.192. The depth of the casing groove 30 is determined through prior calculation. It was determined according to the sensitivity of the objective function, and the width of the casing groove 30 was determined so that the five casing grooves 30 did not overlap each other within the length of the tip cord.

In addition, the twelve experimental points for calculating the target function, namely the stall margin value shown in Equation (1), are used for Latin Hypercube Sampling (LHS), which is useful for sampling a specific experimental point in a design area having a multidimensional distribution. Stall margin values at the 12 experimental points are obtained through 3D Reynolds-averaged Navier-Stokes (RANS) analysis.

As such, when the design variables and the design area are determined, an optimal lattice system is constructed for analysis. In the present invention, a test for removing lattice dependence is performed between 270,000 and 740,000 lattice numbers for a single channel. Results A total of 480,000 grids were used to calculate stall margins.

Next, the numerical analysis step (S30) is to determine the target function value at each experimental point by performing a numerical analysis in the selected design area, in the present invention using the ANSYS CFX-11.0 commercial software ANSYS axial flow compressor Numerical analysis was performed assuming that the internal flow field of is a compressible three-dimensional steady state.

The shape of the axial compressor, the grid generation, the boundary condition, the flow analysis, and the result theorem were analyzed using Blade-Gen, Turbo-Grid, CFX-Pre, CFX-Solver, and CFX-Post, respectively.

The RANS equation for the flow analysis is discretized by the finite volume method, the working fluid passing through the compressor is an ideal state air, the boundary condition of the inlet is given the uniform pressure and the total temperature, and the outlet condition is a single channel. Mass flow rate for.

The axial compressor of the present invention is an axial compressor equipped with NASA Rotor 37, which is a transonic axial compressor blade 10 having five casing grooves 30 and a small aspect ratio to improve stall margin. 37 is composed of 36 blades 10, but assuming that the flow field between two adjacent blades 10 are formed periodically in the rotational direction in order to shorten the calculation time and improve convergence, each in one flow path The analysis was performed by applying periodic conditions.

NASA Rotor 37, the blade 10 used for numerical analysis, has the characteristics shown in Table 2 below.

Design flow rate, ks / s 20.19 Rotational speed, rpm 17188.7 Voltage ratio 2.106 Inlet Hub-Tip Rain 0.7 Blade aspect ratio 1.19 Tip Side Entrance Mach Number 1.48 Hub Side Entrance Mach Number 1.13 Wing tip width ratio 1.29 Number of rotor blades 37

In other words, the NASA Rotor 37 has a power ratio of 2.106 at a design flow rate of 20.19 kg / s, an efficiency of 88.9%, and a tip clearance of 0.356 mm (0.47% span). The three-dimensional shape of NASA Rotor 37 shown in Table 2 is as shown in (a) of FIG. 3, and the flow parameters such as the voltage force according to the flow rate are different from the inlet reference plane as shown in the meridion shape shown in (b) of FIG. And an inlet reference plane and an outlet reference plane are located 41.9 mm upstream and 106.7 mm downstream from the hub's leading edge (LE), respectively.

In addition, when the groove 30 is mounted on the casing 20 of the axial compressor, the interface between the casing 20 and the groove 30 is provided with a GGI (General Grid Interface) interface condition. The SST (Shear Stress Transport) model was used to predict the flow separation due to the back pressure gradient.

In the vicinity of the surface of the rotating impeller blade, a hexahedral type O grid system was used, and in other areas, a hexahedral H / J / C / L type grid system was used. As shown in (b).

In order to verify the validity of the numerical results obtained by the above process, the results of comparing the results of the RANS analysis and the performance test performed in Dunham's "CFD Validation for Through-Flow Calculations" described in the 1998 AGARD report are shown. 5 is shown.

As shown in FIG. 5, the voltage force ratio and efficiency predicted at each flow point are predicted to be somewhat lower than the test results, but the voltage force and efficiency distributions in all areas show the same tendency as the test results. Can be considered valid.

Next, the optimum point search step (S40) relates to the step of searching for the optimum point in the design area based on the result obtained in the numerical analysis step (S30), in the present invention is a kind of surrogate model (surrogate model) Optimal points are calculated using the Radial Basis Neural Networks (RBNN) model.

The neural network model is an algorithm that simulates the human's ability to learn from experience and predict from existing data. The neural network model searches for the optimal point by reflecting weights through the network's predictive ability by the basic elements of neurons.

In more detail, in the numerical analysis step (S30), after evaluating the value of the objective function for the experimental points obtained by Latin Hypercube Sampling (LHS), and constructing a neural network model based on the evaluated objective function Using SQP (Sequential Qua-dratic Programming), we find the optimal point from the neural network ratio method.

In this case, the SQP is a method for optimizing the nonlinear objective function within the nonlinear constraint, and thus a detailed description thereof will be omitted.

In addition, since the optimal point may be changed according to the initial value, the SQP may be changed several times to obtain the final optimal point of the neural network.

As such, a three-dimensional mesh plot of the optimal point optimized by the neural network model is shown in FIG. 6. As can be seen in FIG. 6, the optimum point is 0.177 and 0.164 of the depth (D / C) and the width (W / C) of the grooves dimensionless with the length of the tip cord, respectively, It can be seen that it is included in the range (depth; 0.077 ~ 0.231, width; 0.120 ~ 0.192) (optimal point comparison step (S50)).

Through the series of processes as described above, the optimum shape of the casing groove 30b was determined by the optimum design method according to the present invention. In order to check the performance improvement of the casing groove 30b having the optimum shape determined as described above, Latin Any of the 12 experimental points determined by Latin Hypercube Sampling (LHS) was chosen as the casing groove 30a of the reference shape to compare several performances with a smooth casing.

As shown in FIG. 7, the reference groove 30a has a depth and a width of 0.110 and 0.182, which are dimensionless by the length of the tip cord. The width of the groove is a little narrower than the reference shape, but looking at the height of the groove, since the optimum shape is much higher than the reference shape, it can be seen that each groove area is widened through optimization.

Data comparing the design variables and the objective function of a smooth casing, a casing with a reference groove 30a of a reference shape and a casing with an optimized shape groove 30b. Is shown in Table 3 below.

shape
Design variables Objective function (SM,%) Growth rate,% D / C W / C RBNN RANS SM Remarks Smooth - - - 11.63 - - Reference 0.110 0.182 - 17.23 5.60 Ref.- Smooth Optimum
0.177
0.164
19.32
18.93
7.30 Opt. -Smooth 1.70 Opt. -Ref.

As shown in Table 3, the stall margin (SM) of the compressor having the grooveless shape and the reference shape groove was estimated to be 11.63% and 17.23%, respectively, through RANS analysis. 19.32% using the neural network (RBNN) model and 18.93% through the RANS analysis.

From this, the stall margin predicted by the neural network model can be seen to be relatively accurate when compared with the calculated value through the RANS analysis.

In addition, the stall margin of the axial compressor with grooves of the reference shape is improved by 5.6% compared to the axial compressor without grooves, and the stall margin of the axial compressor with the groove of the optimum shape is 1.7%, compared to the axial compressor with grooves of the reference shape. It can be seen that 7.3% is improved compared to axial compressor without groove.

는 최적점(Optimal point)에서의 목적함수 값을 나타낸다.Next, Figure 8 shows the sensitivity analysis of the design variables of the objective function, wherein the change of the design variables (W / C, D / C) was limited to 10% of the optimal value,

Represents the objective function value at the optimal point.

As a result of the sensitivity analysis shown in FIG. 8, it can be seen that the objective function is more sensitive to the depth (D / C) than the width (W / C) of the groove.

Next, FIG. 9 compares the compressor performance with the grooveless shape and the reference shape and the groove of the optimum shape. First, the flow rate at the near stall point is 0.921 for the shape without the groove, the reference shape. 0.903 for and 0.889 for optimal shape.

In addition, the ratio of the total pressure of the grooveless shape, the reference shape, and the optimum shape in the vicinity of the stall was calculated as 2.073, 2.095, 2.076, and the total temperature ratio was calculated as 1.282, 1.286, 1.287.

Therefore, as compared with the shape without the casing groove, the casing groove is installed to reduce the flow point near the stall, and as a result, it is confirmed that the stall margin is improved, and the stall margin can be further improved compared to the reference shape through the optimization of the groove. can confirm.

In addition, when the groove is installed, the total temperature ratio is constant over the entire operating range while the total pressure rises, indicating that the high momentum flow in the groove reduces the delamination that causes pressure loss near the blade suction surface. It becomes possible.

Next, (a) to (c) of FIG. 10 show the grooveless shape, the reference shape, the optimum shape of the voltage force, and the tip leakage vortex at the vicinity of the stall of the shape without the casing groove, and one of the effects of the casing groove is It is shown that the trajectory of the tip leakage vortex is changed by changing the angle between the blade and the tip leakage vortex.

That is, the tip leakage flow due to the pressure gradient at the blade tip interacts with the inflow flow in the flow path to form a tip leakage vortex. In the case of a compressor without grooves, the tip leakage vortex is mainly influenced by the tip leakage flow. Proceed along the pressure side of the blade as shown in 10 (a).

Compared to the grooveless shape, the tip leak vortex of the grooved compressor is highly influenced by the flow in the groove, as shown in FIGS. 10 (b) and 10 (c), where tip leak vortex is formed along the groove. Will be. As a result, the diffusion of the tip leakage vortex occurs due to the pressure rise generated after the vortex passes the shock, which creates a region of low energy.

Next, (a) to (c) of FIG. 11 show the Mach number for the grooveless shape, the reference shape, and the optimum shape at the point near the stall of the grooveless shape. As shown, in the compressor without grooves, separation occurs at the suction side of the blade, whereas in grooved compressors, the separation at the blade is suppressed, as shown in FIGS. 11B and 11C. You can see it separated.

The phenomenon that the peeling is suppressed as described above is because the flow of the high momentum is transmitted from the pressure surface to the suction surface through the groove to suppress the peeling from the blade.

Next, (a) and (b) of FIG. 12 show Mach number distributions of the reference shape and the optimal shape at the point near the stall of the reference shape, immediately downstream of the shock line by the action of tip leakage vortex and flow path impact. It can be seen that a low speed area occurs in the area.

This low velocity region is the Vortex stagnation zone, which changes the main flow and raises the back pressure to cause delamination. Within the vortex stagnation zone the flow has no axial velocity component and only a circumferential velocity component.

Next, (a) and (b) of FIG. 13 show velocity vectors of the reference shape and the optimum shape in the AB cross-sections in FIGS. 12 (a) and (b), and refer to the circulating flow in the vortex stagnant region. Both shape and optimum shape can be found near the blade tip.

Compared with the reference shape shown in FIG. 13A, the optimum shape shown in FIG. 13B has a greater depth through optimization, which can absorb more circulating flow.

In addition, it can be seen that the circulating flow in the reference shape of the groove proceeds back to the flow path, whereas the circulating flow in the optimum shape is mitigated by the opposite flow in the groove, which reduces the flow path blocking phenomenon in the flow path and the compressor It will help to improve the stability and stall margin of.

Therefore, according to the casing groove optimum design method and optimally designed casing groove of the axial compressor according to the present invention, the shape of the casing groove can be optimized in the entire design area through systematic optimization design using optimization theory, not limited testing through trial and error. By making it possible, not only can the performance of the axial compressor be improved, but the casing groove can be designed to continuously improve the stall margin.

Although the above embodiments have been described with respect to the most preferred examples of the present invention, it is not limited to the above embodiments, and it will be apparent to those skilled in the art that various modifications are possible without departing from the technical spirit of the present invention.

The present invention relates to a casing groove optimum design method and optimally designed casing groove of the axial compressor, and more specifically, to the casing groove through the numerical definition of the shape of the axial compressor, lattice generation, boundary conditions and internal flow field Performance of the grooved axial compressor and the basic grooved axial compressor and the grooved axial compressor optimized by the present technology to verify the maximum stall margin in the optimized grooved axial compressor. The present invention relates to a casing groove optimum design method and optimally designed casing groove of an axial compressor, which enables the design of a casing groove to continuously improve the stall margin.

10: blade 20: casing
30: Groove S10: formalization stage of the problem
S20: Design Area Selection Step S30: Numerical Analysis Step
S40: Optimal point search step S50: Optimal point comparison step

Claims (5)

Formulating a problem of selecting a design variable of the casing groove in consideration of the shape of the axial compressor and determining an objective function;
A design area selection step of setting a design area by limiting a range of design variables;
A numerical analysis step of performing numerical analysis in the selected design area; And
And an optimum point search step for searching for an optimum point in the design area based on the numerical analysis result.
The method of claim 1,
The design variable is the casing groove optimum design method of the axial compressor, characterized in that the width and width of the casing groove dimensioned by the length of the tip cord.
The method of claim 1,
The objective function is

Casing groove optimum design method of the axial compressor, characterized in that the stall margin of the axial compressor.
(At this time, SM: Stall margin,

= Flow rate at the maximum efficiency point,

Flow rate near the stall

= Voltage force at the maximum efficiency point,

: Flow rate near stole.)
The method of claim 1,
In the optimum point search step, the casing groove optimum design method of the axial compressor, characterized in that for determining the optimal design parameter value of the casing groove from the neural network technique.
In the casing groove of the axial compressor,
A casing groove of an axial compressor, characterized in that the depth and width of the casing groove dimensionless with the length of the blade tip chord are 0.177 and 0.164.

Images