Optimization of the Impeller Geometry for an Automotive Torque Converter Using Response Surface Methodology and Desirability Function ()
1. Introduction
Torque converters are widely used in vehicle power transmission systems and transfer the torque and power from engine to transmission by the hydraulic operation of transmission oil. Despite of its ability to provide good acceleration performance and absorb excessive vibration, its mechanical efficiency is relatively low. Nowadays, a few new techniques are in use to overcome the low efficiency. The application of the locking clutch of the turbine to the impeller and the optimized design of shape of torus and blades to enhance the hydraulic performance are the two representative methods. The first technique is very effective at the wider range of vehicle speed, but not easy to obtain an optimized shape since a large number of design parameters are involved. Typically, several torque converters are designed and tested by computational fluid dynamics (CFD) and experiment until the target performance is met [1] . In an automatic transmission, the output torque and power of the engine is imparted to the transmission oil by the impeller of the torque converter. The impeller provides the pressure rise and increases the angular momentum of the oil through the rotation and the flow turning [2] . Therefore, it is very important to understand the influence of impeller geometrical parameters on performance characteristics of a torque converter. This knowledge will greatly help designers to improve the torque converter performance.
It is generally been accepted that while the accuracy of CFD analyses has not yet achieved a level that is equivalent to experimental techniques, its ability to correctly predict the direction of any changes is reliable [3] . In recent years, CFD has been widely used in turbomachinery design and optimization. Rutter et al. [4] investigated the hydraulic performance of a centrifugal pump within the electrical submersible pump (ESP) unit in single-phase flow using CFD technique. Zhao, et al. [5] optimized a double-channel pump’s impeller by combined using CFD, multi-objective genetic algorithm (MOGA) and artificial neural networks (ANN). Shojaeefard, et al. [6] , and Bellary, et al. [7] , improved the performance of a centrifugal pump impeller based on CFD simulations. Hur, et al. [8] analyzed the flow and performance of a partially-charged water retarder by CFD. Many researchers have also studied the flows in torque converters by using CFD codes employing various methods [9] [10] [11] . In order to improve the converter performance, it is required to obtain detailed understanding and relationship between the governing parameter and its effect on the performance, including efficiency, torque ratio and the input capacity factor. Kubo, et al. [12] analyzed the effects of the turbine bias angle and the impeller flow passage contraction ratio on the internal flow characteristics and the occurrence of loss. Shin, et al. [13] investigate the effect of the flow section area, flow rate, impeller blade angle and scroll angle of impeller and turbine on the performance of a torque converter. The analytical and numerical results obtained a refined relationship between geometry and performance. Shin, et al. [14] [15] investigated the effect of reactor blade geometry with varying thickness ratios, scroll angles and slot angles on the performance of a torque converter. However, the existing studies are limited to investigate the influence of only one or two geometrical parameters individually each time, and reports in the literature regarding their combinations effect on the torque converter performance characteristics can hardly be found.
The traditional one-factor-at-a-time approach has been widely used for evaluation or optimization of these parameters [16] [17] [18] . The major disadvantage of this approach is that it fails to consider any possible interaction between factors, and therefore it may lead to misinterpretation of results. However, design of experiments (DOE) and CCD statistical approaches enable the simultaneous study of several factors and assessment of their statistical significance, as well as the evaluation of interaction effects. DOE method is widely used to find the importance level of the design parameters with respect to the optimization target and obtain the best combination of design variables. Park, et al. [19] studied a methane-fueled gas engine generator with addition of hydrogen using DOE method. Wu, et al. [20] studied the influence of volute geometrical variations on hydrodynamic characteristics of a high speed circulator pump by using CFD technology and DOE method. Taghavifar, et al. [21] applied DOE evaluation to introduce the optimum injection strategy-chamber geometry of diesel engine. CCD method is widely applied to evaluate the interaction effect and explore the optimum combination of key parameters. Sagbas [22] optimized ball burnishing process using rotatable central composite design (CCD) together with desirability function approach. Hatami, et al. [23] applied CCD approach to obtain an optimal design of the vane geometry for a variable geometry turbine. Although several studies [24] [25] [26] were reported using various methods to optimize torque converters, application of DOE combined with CCD to optimize a torque converter is scarce.
In this present paper, the main objective is to improve the overall performance of automotive torque converters by means of impeller geometry optimization using DOE and CCD methods. A new parametric geometric design method of impeller is proposed by using parametric equations and Creo software. A DOE array is used to reduce the number of analyses required at each flow area factor, with two area factors being investigated, namely 17% and 20%. The DOE post-processing analysis is also used to rank the relative importance of the geometric parameters for both area factors. Then, the CCD method is applied to in-depth evaluate the interaction effect of the three most important parameters and the optimization is done to have the best overall performance of a torque converter.
2. Parametric Model of Impeller
2.1. Torus Design
The impeller parametric design starts with the definition of the torus shape. Various design parameters including active diameter, aspect ratio and two arc radii are needed to determine the shell profile of an automotive torque converter torus. As the flow area in the circular path, in the proposed torque converter model, is assumed constant, only the design parameter area factor
is used to determine the design path and core of torus. Consequently, by defining the impeller inlet radius
and outlet radius
, along with the shell and core, the torus profile of impeller can be obtained as shown in Figure 1.
2.2. Blade Design
Each impeller blade profile (Figure 3(c)) is formed by a three dimensional (3D) curve of the shell and a 3D curve of the core. The 3D curve can be calculated from the torus and a two dimensional (2D) design curve using conformal transformation principle. Figure 2 shows the schematic representation of a 2D design curve. With the origin at the starting point of the curve, the coordinate of
is (0, 0). The curve expression can be written with the following equation:
, (1)
where A, B, C, D and E are the coefficients of the variables. It should be noted that the value of D is 1.0 and the parametric equation of the 2D design curve can be obtained as
, (2)
, (3)
, (4)
, (5)
where
is the exit angle of impeller, and
is the inlet angle of impeller. In the present study, four design parameters including exit angle
, inlet angle
, offset size d, and conic factor
are provided to calculate the 2D design curve. The coordinates of
and
can be obtained as
, (6)
, (7)
, (8)
, (9)
where
is the intersection of the two tangent lines that across the curve’s starting and ending points, respectively (Figure 2). According to the definition of 2D design curve,
value and
value equal to the length of torus l and offset size d, respectively.
The 2D design curve can be easily constructed by having the torus and the four design parameters defined. After the definition of 3D curves of shell and core, another design parameter bias angle
is used to obtain the blade profile. Given the impeller blades are constant-thickness stamped sheet metal, the blade thickness t is defined. Once the torus shape and five-blade design parameters including exit angle
, inlet angle
, offset size d, conic factor
, bias angle
and blade thickness t are determined, the basic blade geometry of impeller is generated by means of parametric equations and Creo software. The schematic representation of the construction of impeller blade is shown in Figure 3.
After the definition of the torus and a blade, the construction of the whole impeller geometry is easy and is based on the rotation of the generic blade around the axis; for doing so the number of blades z should be provided by the user. The completed impeller parametric geometry is shown in Figure 4. The design parameters considered for the parametric study are shown in Table 1, provided that the torus of impeller is determined.
Figure 2. Schematic representation of 2D design curve.
(a) (b) (c) (d)
Figure 3. Schematic representation of the construction of impeller blade: (a) 2D design curve for impeller blade, (b) 3D curve of impeller, (c) blade profile of impeller, (d) blade geometry of impeller.
3. Model Verification
3.1. Computational Approach
An automotive torque converter is selected as a reference to generate the parametric model. The reference torque converter has an active diameter of 250 mm and the number of blades in the impeller, turbine, and stator are 31, 29, and 21, respectively. STAR-CCM + software is used to generate the computational mesh and perform the internal flow calculations in an appropriate way. The computational mesh is given in Figure 5 where only one blade passage is modeled for each element to illustrate the mesh distribution in the computational field when approximately 111,000 grid cells in total are used. The polyhedral cells with five prism layers are used on the model. The leakage between the elements and also between an element and the core flow are disregarded. A cyclic boundary condition is imposed on both peripheral boundaries outside a blade passage. A no-slip wall boundary condition is also imposed on all the walls bounding the domain, with a spin applied as necessary. The interfaces between elements have been handled by using the mixing plane method. A second-order upwind differencing scheme is utilized and the SST k-w model is also used for the turbulence. Steady state simulations are performed for a range of speed ratios from 0.0 to 0.9 while maintaining an impeller speed of 2000 rpm.
3.2. Comparison between Simulation and Experimental Results
Figure 6 compares the measured and calculated overall performance including efficiency
, torque ratio
, and impeller torque factor
for the parametric model of the torque converter. As indicated here, the tendencies of the experimental data correlated relatively well with the calculated results, confirming that the parametric model and computational method are valid in general. It is of note that the CFD underestimates the impeller torque factor
at low speed ratios. The discrepancy between measurements and calculations are reasonable because of the overestimated incidence losses at low speed ratios.
4. Design of Experiments (DOE)
4.1. Influence of Flow Area
When the circuit size is determined with satisfying the performance requirements, the flow cross-section can be determined with a rule of thumb that the flow area is uniform throughout the blades. In the present study, three toque converter parametric models with different flow areas are designed and simulated to investigate the effects of flow area on their overall performance including efficiency and impeller torque factor. For comparison, the flow area is altered to 17%, 20%, and 23% of area of a circle represented by the converter diameter (fa = 17%, 20%, and 23%). The design parameters are unchanged as much as possible except the area factors and simulation models with varying area factors are shown in Figure 7.
In Figure 8, the performance analysis shows the significant increases of impeller torque factor with the flow area increases, whereas efficiency calculation does not show significant increase as shown in Figure 9. It can be concluded that the impeller torque factor is more sensitive than efficiency to the torque converter’s flow area.
Figure 6. Comparison of the CFD calculation analysis with experimental results.
(a) (b) (c)
Figure 7. Simulation models with varying area factors: (a) fa = 17%, (b) fa = 20%, and (c) fa = 23%.
Figure 8. Effects of flow area on impeller torque factor of torque converters.
Figure 9. Effects of flow area on efficiency of torque converters.
4.2. Factors and Levels in DOE
DOE is a collection of mathematical and statistical techniques to reduce the number of experiments in order to find the effect of parameters affecting a response in a process, thereby aiming for a reduction in both costs and time [27] [28] [29] . A DOE method sets out configurations (or arrangements) to be conducted using an appropriate orthogonal array; the terminology used in these arrays includes “factors”—an item that is to be varied during the simulations, “level”—the number of times a factor is to be varied during the simulations and “configuration number”—the number of simulations that are required to be run to complete the analysis [30] . The selection of an appropriate DOE array is dependent on the number of factors and the levels of the factors to be analyzed. As mentioned above, 11 design parameters are used to obtain an impeller parametric model. Blade scroll angle
is introduced as another design parameter to investigate the effect of the impeller geometry on the performance of a torque converter. The scroll angle is defined as the angle between the two planes containing the intersection of the design path and the entering and leaving edges of the blade when that blade does not lie in one axial plane. Impeller blades with varying scroll angles are shown in Figure 10.
The scroll angle can be determined by four design parameters including exit angle
, inlet angle
, conic factor
, and offset size d, provided that the torus is defined. Matlab software is used to develop a simple-to-use GUI to calculate scroll angle as shown in Figure 11.
The blade scroll angles at shell and core can be represented by the blade scroll angle on the design path. Similarly, the inlet angles and exit angles at shell and core can be represented by the inlet angle and exit angle on the design path. Finally, the 11 design parameters of impeller are translated into 6 parameters including exit angle
, inlet angle
, scroll angle
, bias angle
, blade thickness t, and blade number z. The selected factors and levels in DOE are listed in Table 2. The underline sign in the table indicates the base value of impeller design parameters.
4.3. DOE Evaluation
For this paper, six main geometrical parameters mentioned above are selected as design variables (factors) and five different values (levels) are assigned for each design parameter. So 25 (L25[56]) configurations with different combinations are generated for DOE. Stall torque ratio
and peak efficiency
are used as the dynamic characteristic and economic characteristic, respectively, to evaluate the performance of torque converters. Two area factors including 17% and 20% are considered. The final design matrix in DOE is presented in Table 3 along with the responses predicted by CFD.
One aspect of the DOE method is utilized response averages, calculated for each response and area factor in relation to a specific geometry parameter variable, to provide detail relating to the influence of the geometric factors on the performance of a torque converter. For example, to calculate the response averages relating to the 27 blade number, the average of the stall torque ratio and peak efficiency would be calculated from case number 1, 2, 3, 4 and 5. The calculated response average values for each level and for all parameters are shown for each area factor in Table 4, Table 5.
Figure 10. Impeller blades with varying scroll angles.
Figure 11. GUI interface used for scroll angle calculation.
Table 2. Factors and levels in DOE.
Table 3. Design matrix and results in DOE.
Table 4. Response averages and percentage contributions for area factor 17%.
Table 5. Response averages and percentage contributions for area factor 20%.
One other important aspect of the DOE method is the ability to calculate the percentage contribution of a geometric parameter to a specified performance characteristic. The percentage contributions are calculated from an analysis of variance that effectively measures how far the performance characteristic values for a specific geometry variable vary from the mean. The amount that the high and low parameter levels vary from the mean provides a measure of that parameter’s influence on a particular performance characteristic. This is converted into a percentage value to provide a measure of the contribution relative to the other parameters. Again percentage contributions have been calculated for each area factor. A summary of the percentage contributions for stall torque ratio
and peak efficiency
at area factors 17% and 20% is shown in Figure 12.
The percentage contributions provide a great deal of information regarding the importance of the various geometric factors of impeller. As shown in Figure 12, exit angle
has the strongest effect on peak efficiency at both 17% and 20% area factors although this effect reduces for the stall torque ratio. It should be noted that while the peak efficiency contributions for the bias angle
are relatively low, the contributions to the stall torque ratio at both area factors are surprisingly large. The inlet angle
and blade number z show a lesser contribution while the blade thickness t generally has a significantly lesser effect. It can be concluded that the inlet angle
, bias angle
, and exit angle
exert the greatest influence on the overall performance of a torque converter.
The DOE array provides useful information in the form of the percentage contributions, but not provide information on interaction effects between the geometric parameters. Later, a CCD technique will be used to gauge the interactive effect among these three dominant factors.
Figure 12. Percentage contributions of influence factors (1 - 6 correspond to z, t,
,
,
and
respectively).
5. Optimization Analysis
5.1. Central Composite Design (CCD)
The CCD consists of a two-level full or fractional factorial design (corner points), an additional design (star points) and at least one point at the center of the design space (center points) (Figure 13). If the distance from the center of the design space to a factorial point is ±1 unit for each factor, the distance from the center of the design space to a star point is
with
. The
value equals to (2m)1/4, where m is the number of factors [31] [32] . In the present study, m equals to 3 so that the
value is 1.682. Therefore, each factor is studied in five different levels (−1.682, −1, 0, +1, +1.682), which is listed in Table 6.
According to CCD, 15 torque converter cases should be modeled and their properties are presented in Table 7. It should be noted that nine points are considered for central points and two responses are considered for the optimization, stall torque ratio and peak efficiency. The corresponding calculation results for two area factors are presented in Figure 14 and Figure 15. Figure 14 reveals that case 4 has the maximum stall torque ratio, but it has a low peak efficiency and case 10 has the maximum peak efficiency. As seen in Figure 15, case 10 also has the best peak efficiency, but it has the minimum stall torque ratio and case 9 has the maximum stall torque ratio. So, an optimization study is needed which is performed using CCD technique in the present study.
5.2. Evaluation and Optimization
In CCD, a polynomial model with quadratic order is applied to responses (stall torque ratio and peak efficiency). P-value Probs are estimated to be 0.0024 and less than 0.0001 while R-squared values are 0.8002 and 0.9909 for
and
, respectively, for 20% area factor. For 17% area factor, P-values are also estimated to be less than 0.0001 while R-squared values are 0.9596 and 0.9957 for
and
, respectively, which are meaningful. By using the quadratic equation for the surfaces, the stall torque ratio and peak efficiency can be estimated by the following equations for different area factors.
(a) (b)
Figure 14. (a) Stall torque ratio and (b) peak efficiency for all designed cases at 17% area factor.
(a) (b)
Figure 15. (a) Stall torque ratio and (b) peak efficiency for all designed cases at 20% area factor.
Table 6. Factors and levels in CCD.
Table 7. Design matrix and results in CCD.
1) 17% area factor:
(10)
(11)
2) 20% area factor:
(12)
(13)
By applying the above formula one can approximate the stall torque ratio and peak efficiency for all non-simulated cases.
As mentioned above, the interaction of each parameter with other parameters should be considered at the same time on stall torque ratio and peak efficiency. Therefore, the response surfaces are presented in Figure 16 and Figure 17 for each area factor, respectively, to show this interaction on
and
. These figures show the interaction between the three parameters for the optimized design. An increase in bias angle
and inlet angle
together with an decrease in exit angle
yield higher stall torque ratio amounts although the effect of bias angle is more sensible than inlet angle and exit angle due to its high slope in the surface in Figure 16(a), Figure 16(b), Figure 17(a) and Figure 17(b). As an important outcome, exit angle has the most significant effect on peak efficiency, while bias angle and inlet angle may have optimum values that maximize the peak efficiency (Figure 16(c), Figure 16(d), Figure 17(c), and Figure 17(d)). It is of note that the minimum value of peak efficiency at 20% area factor is lower than that at 17% area factor as shown in Figure 16(c) and Figure 17(c). This phenomenon can be explained by the fact that the curvature of the core increases with the increase of flow area of a torque converter, which may result in larger flow losses.
(a) (b)(c) (d)
Figure 16. Response surface of stall torque ratio and peak efficiency for optimization at 17% area factor. (a) Exit angle and bias angle, (b) inlet angle and bias angle, (c) inlet angle and exit angle, and (d) exit angle and bias angle.
(a) (b)(c) (d)
Figure 17. Response surface of stall torque ratio and peak efficiency for optimization at 20% area factor. (a) Exit angle and bias angle, (b) inlet angle and bias angle, (c) inlet angle and exit angle, and (d) exit angle and bias angle.
With CCD, optimization is based on a parameter called “desirability”. Desirability is an objective function ranging from 0.0 outside of the limits to 1.0 at the goal. The numerical optimization finds a point that maximizes the desirability function. The characteristics of the goal may be altered by adjusting the weight or importance. For several responses and factors, all goals get combined into one desirability function [23] . In this paper, two responses are defined as stall torque ratio and peak efficiency. The goal of optimization is to find a set of conditions that meet all the goals, not to get a desirability value of 1.0. Desirability reflects the preferred ranges for each response (
). The simultaneous objective function is a geometric mean of all transformed responses:
(14)
where n is the number of responses in the measure (in this case, n = 2). If any of the responses or factors fall outside their desirability range, the overall function becomes zero. For simultaneous optimization, each response must have a low and high value assigned to each goal. In this study, the goal parameter used is “maximum” (for both stall torque ratio and peak efficiency) as follows:
(15)
where Yi is the ith response value and wt is the weight of that response. Weight adds emphasis to the goal. A weight greater than 1 (maximum weight is 10), emphasizes the goal and less than 1 (minimum weight is 0.1), deemphasizes the goal. In this paper, two responses are defined, the weights will be determined according to the designer’s performance demand.
After optimization analysis by above surface formula, CCD proposes three optimized cases for each area factor (Table 8). Optimized cases 1 and 4 have the maximum stall torque ratio while optimized cases 2 and 5 have the best peak efficiency. Optimized cases 3 and 6 have the largest desirability when the weights of responses are considered equal. It is clear that torque converters with 20% area factor have better performance characteristics than torque converters with 17% area factor.
In the present study, optimized case 6 is selected as the final optimized case which is designed and calculated. Table 9 compares the theoretical and calculation results. The CFD calculation results at the optimum parameter combination confirm the effectiveness of the response surface model for optimum design parameters. Compared to the original model, the stall torque ratio and peak efficiency of the optimized one have increased by 1.62% and 1.1%, respectively. Figure 18 shows the performance comparison of the optimized design with the original one which illustrates the improvement of performance.
6. Conclusions
A new parametric design method of the impeller for an automotive torque converter has been used to conduct a parametric study covering 11 geometric parameters. The 11 geometric parameters can be represented by 6 design parameters including the impeller blade number z, blade thickness t, bias angle
, scroll angle
, exit angle
and inlet angle
.
Table 8. Optimized cases proposed by CCD.
Table 9. Comparison of the calculation and CCD proposed values for optimized case 6.
Figure 18. Performance comparison of the optimized design with the original model.
An L25[56] DOE array has been successfully constructed for two responses including stall torque ratio and peak efficiency. The analysis of the array identified the dominant geometrical influences on the performance of the torque converter. In general, the impeller exit angle
, inlet angle
and bias angle
are the three strongest influences on the overall performance, with two area factors (17% and 20%) considered.
CCD method has been employed to investigate the interactive effect among the impeller inlet angle, exit angle and bias angle. Predictive equations are presented that can identify expected performance at specific arrangements different from those analyzed. Results show that larger impeller bias angle and inlet angle together with smaller exit angle bring about an increase of the stall torque ratio while a larger exit angle is favorable to increase the peak efficiency. A final optimized design of impeller is obtained based on desirability function approach together with the predictive equations. Compared with the original model, the stall torque ratio and peak efficiency of the optimized model are increased by 1.62% and 1.1%, respectively.
Performance prediction models with more prediction accuracy for a torque converter can be used as a further research direction.