Support Vector Regression Methodology for Performance Prediction of an Automotive Torque Converter with Turbine Parametric Model ()
1. Introduction
The torque converter is an important component in an automatic transmission. It provides automatic torque amplification according to the different rotational speed between the input and output speeds without any active control, inherently suppressing engine torque fluctuations [1]. An impeller, a turbine, and a stator all are a part of a torque converter. The impeller energizes the working oil and is connected to the engine, while the turbine extracts energy from the fluid and drives the transmission shaft. Flow field in turbine passage is quite complex for the reason that the rotation speed of turbine changes with the speed ratio changes. To get a larger torque ratio at low speed ratio, the turbine blade angles always change drastically. Therefore, it is very important to understand the relationship between turbine geometrical parameters and the 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 computational fluid dynamics (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 [2]. In recent years, CFD has been widely used in turbomachinery design and optimization. Rutter et al. [3] 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. [4] optimized a double-channel pump’s impeller by combined using of CFD, multi-objective genetic algorithm (MOGA) and artificial neural networks (ANN). Shojaeefard et al. [5], and Bellary et al. [6], improved the performance of a centrifugal pump impeller based on CFD simulations. Hur et al. [7] 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 [8]-[10].
Recently, support vector machine (SVM) has been extensively applied to various areas to overcome the problem of non-linear relationships and predictions [11]-[13]. There are two main categories for support vector machines: support vector classification (SVC) and support vector regression (SVR). SVM is a learning system using high-dimensional features. SVR is based on statistical learning theory and a structural risk minimization principle, which has been successfully used for non-linear system modeling. Rajasekaran et al. [14] applied support vector regression methodology for storm surge predictions. Bermolen, et al. explored the use of SVR for the purpose of link load forecast. Shamshirband, et al. [15] presented the SVR methodology for sensor data fusion to improve tracking abililty. Taghavifar et al. [16] optimized the injection strategy-chamber geometry of diesel engine using DOE method incorporating with the SVR to predict output parameters with acceptable accuracy. Although several studies [17]-[20] were reported using various methods to predict torque converter performance, application of SVR method to predict torque converter performance is scarce.
A kernel function can be utilized to form qualified function that is used SVM. Researchers have demonstrated the use of SVR in predicting hydrological modeling and pointed out the positive performance of the RBF (radial basis function) [21] [22]. Asefa et al. [23] proposed different kernels in SVR to rainfall-runoff modeling and validated that the radial basis function (RBF) outperforms other kernel functions. In this paper, the SVR for prediction of the torque converter performance including stall torque ratio and peak efficiency is applied, with two kernel functions being investigated, namely RBF and polynomial kernel function.
The basic idea behind the SVR methodology is to collect input/output data pairs and to learn the proposed network from these data. In our scheme, the input data comes from design parameters of turbine geometry for an automotive torque converter. Design of experiment (DOE) together with CFD techniques is used to obtain the output data. To improve the design efficiency, a new parametric geometric design method of turbine is proposed by using parametric equations and Creo software.
2. Parametric Model of Turbine
2.1. Torus Design
The turbine 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 turbine inlet radius
and outlet radius
, along with the shell and core, the torus profile of turbine can be obtained as shown in Figure 1.
Figure 1. Torus profile of turbine.
2.2. Blade Design
Each turbine blade profile 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 the 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 inlet angle of turbine, and
is the exit angle of turbine. In the present study, four design parameters including inlet angle
, exit angle
, offset size
, 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
and offset size
, respectively.
Figure 2. Schematic representation of 2D design curve.
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 turbine blades are constant-thickness stamped sheet metal, the blade thickness
is defined. Once the torus shape and five blade design parameters including inlet angle
, exit angle
, offset size
, conic factor
, bias angle
and blade thickness
are determined, the basic blade geometry of turbine is generated by means of parametric equations and Creo software. The schematic representation of the construction of turbine blade is shown in Figure 3.
Figure 3. Schematic representation of the construction of turbine blade: (a) 2D design curve for turbine blade, (b) 3D curve of turbine, (c) blade profile of turbine, (d) blade geometry of turbine.
After the definition of the torus and a blade, the construction of the whole turbine geometry is easy and is based on the rotation of the generic blade around the axis; for doing so the number of blades
should be provided by the user. The completed turbine 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 turbine is determined.
Figure 4. Turbine parametric geometry.
Table 1. Parameters of turbine.
No. |
Description |
Parameter |
1 |
Exit angle of turbine at shell (˚) |
|
2 |
Inlet angle of turbine at shell (˚) |
|
3 |
Conic factor of turbine blade at shell |
|
4 |
Offset size of turbine blade at shell (mm) |
|
5 |
Exit angle of turbine at core (˚) |
|
6 |
Inlet angle of turbine at core (˚) |
|
7 |
Conic factor of turbine blade at core |
|
8 |
Offset size of turbine blade at core (mm) |
|
9 |
Bias angle of turbine (˚) |
|
10 |
Blade number of turbine |
|
11 |
Blade thickness of turbine (mm) |
|
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 appropriately. 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 103, 533 grid cells in total are used. The polyhedral cells with five prism layers is used on the model. The leakage between the elements and also between an element and the core flow is 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,
Figure 5. Torque converter grid model.
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.
Figure 6. Comparison of the CFD calculation analysis with experimental results.
4. Data Collection and Support Vector Regression (SVR)
4.1. Input Parameters
The ability of the SVR to make reasonable estimations is mostly dependent on input parameters selection. Adequate consideration of the factors of turbine effecting the torque converter performance studied is therefore crucial to developing a reliable network. As mentioned above, 11 design parameters are used to obtain an turbine parametric model. Blade scroll angle
is introduced as another design parameter to investigate the effect of the turbine 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. The scroll angle can be determined by four design parameters including inlet angle
, exit angle
, conic factor
, and offset size
, 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 7.
![]()
Figure 7. GUI interface used for scroll angle calculation.
The blade scroll angles at shell and core can be represented by the mean scroll angle and change simultaneously with the change of the mean value. Similarly, the inlet angles and exit angles at shell and core can also be represented by the mean value of inlet angle and mean exit angle. Finally, the 11 design parameters of turbine are translated into 6 parameters including inlet angle
, exit angle
, scroll angle
, bias angle
, blade thickness
, and blade number
. The six design parameters are to be defined as input for the learning techniques in this research study. The range of design parameters and base values of the turbine are summarized in Table 2.
Table 2. Range of design parameters and base values of the turbine.
Design parameters |
Base value |
Range of variation |
Blade number
|
29 |
25 - 33 |
Blade thickness
/mm |
1 |
0.9 - 1.3 |
Bias angle
/(˚) |
2.75 |
−1 - 5 |
Scroll angle
/(˚) |
−1.035 |
−3 - 1 |
Inlet angle
/(˚) |
34.75 |
25 - 45 |
Exit angle
/(˚) |
150.75 |
140 - 160 |
4.2. Design of Experiment (DOE)
4.2.1. Orthogonal Design Approach
Orthogonal design approach 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. Orthogonal design 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 [24]. 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 orthogonal design. Stall torque ratio
and peak efficiency
are used as the dynamic characteristic and economic characteristic, respectively, to evaluate the performance of torque converters. The above mentioned CFD simulation method is applied and the corresponding calculation results for 25 cases are presented in Figure 8. The 25 sets of calculation data are used for SVR modeling training.
![]()
Figure 8. (a) Stall torque ratio and (b) peak efficiency for all designed cases.
4.2.2. Latin Hypercube Design (LHD)
The Latin Hypercube design is a popular choice of experimental design when computer simulation is used to study a physical process. These designs guarantee uniform samples for the marginal distribution of each single input [25]. Five cases are obtained using LHD method and the corresponding calculation results are shown in Figure 9. The 5 sets of data as testing data to analyze the prediction performance of SVR.
Figure 9. Design variables and indicators.
4.3. Support Vector Regression (SVR)
SVR in its essence is closely related with SVM theory and its parameters such as kernel and operating process are not disparate. The
-insensitive loss function is the most widely used cost function. It can be assumed as a tube equal to the approximation accuracy that surrounds the training data [26]. The objective held by SVR is to search for the function with the most
deviations from the real destination vector for all training information received and which must be as flat as possible [27]. For linear regression problems, the regression model can be expressed by Vapnik [28] and is defined by the following function:
(10)
where
is an unknown target function,
denotes the dot product and
is the weight vector. The
-insensitive loss function is defined by the following function:
(11)
As a convex optimization problem, it can be written as:
Minimize
(12)
Subject to
(13)
where
and
variables are to satisfy the function constraints. The corresponding dual optimization problem is defined as [29]:
(14)
with constraints:
(15)
where
and
are Lagrange variables while
and
are denoted by the following equation and
and
are support vectors [28].
(16)
For nonlinear regression problems, a nonlinear mapping
of the input space onto a higher dimension feature space can be used, and then linear regression can be performed in this space [30], the nonlinear model is written as
(17)
where
(18)
(19)
(20)
where
and
are support vectors.
In this paper, the non-linear RBF kernel and polynomial kernel are investigated and are defined as [31]:
RBF:
(21)
Polynomial:
(22)
where
,
, and
are kernel parameters.
4.4. Model Performance Evaluation Criteria
The root mean square error
and correlation coefficient
are used to estimate the accuracy of the prediction model:
(23)
(24)
in which
is the value of simulation and
is the value of prediction.
5. Results and Discussion
A series of simulation have been conducted for SVR modeling and prediction, and the training data sets and testing data sets used in this paper are obtained using orthogonal design approach and LHD method, respectively. In order to eliminate dimension differences, min/max normalization method is used for data standardization and normalization and then all input and output data are standardized and normalized to the range [0, 1].
5.1. SVR Parameters Optimization
The SVR parameters must be set carefully in order to construct the SVR model efficiently [32]-[34]. Inappropriately chosen SVR parameters will result in over-fitting or under-fitting, and different parameter settings may also cause significant differences in performance [35]. Thus, selecting the optimal parameters is an important step in SVR design. Cross-validation together with grid-search method is applied to obtain the SVR optimal parameters. In
-fold cross-validation, we first divide the training set into
subsets of equal size. Sequentially one subset is tested using the classifier trained on the remaining
subsets. Thus, each instance of the whole training set is predicted once so the cross-validation accuracy is the percentage of data that are correctly classified [31]. In our scheme, 5-fold cross validation is selected and the values of epsilon and tolerance are 0.01 and 0.001 seem to perform well. To select optimal combinations of
and
for radial basis function kernels, a coarse grid-search is used first. The contours of coarse grid search results for stall torque ratio and peak efficiency predictions are shown in Figure 10.
![]()
Figure 10. Coarse grid search results for (a) stall torque ratio and (b) peak efficiency predictions.
After identifying a “better” region on the grid, a finer grid search on that region can be conducted. The three dimensional plot of finer grid search results for stall torque ratio and peak efficiency predictions are presented in Figure 11.
Figure 11. Finer grid search results for (a) stall torque ratio and (b) peak efficiency predictions.
It can be seen that the grid optimization process for stall torque ratio yields to a minimum
when
while the best
for peak efficiency is (4.5948, 0.26794). The similar grid-search method is also employed for SVR with polynomial kernel function. Table 3 provides the optimal values of parameters for this data set with polynomial and RBF kernel-based SVR.
Table 3. Optimal parameters for SVR with RBF and polynomial function.
SVR parameters |
Stall torque ratio
|
Peak efficiency
/% |
RBF |
Polynomial |
RBF |
Polynomial |
Epsilon value
|
0.1 |
0.1 |
0.1 |
0.1 |
Cost value
|
3.0314 |
0.125 |
8 |
0.022097 |
Tolerance value
|
0.001 |
0.001 |
0.001 |
0.001 |
Degree of kernel
|
- |
2 |
- |
2 |
Gamma value of kernel
|
0.37893 |
0.75786 |
0.18946 |
7.4643 |
Coefficient value of kernel
|
- |
0 |
- |
0 |
5.2. SVR Prediction
The training data is used to establish the RBF and polynomial kernel-based SVR model for stall torque ratio and peak efficiency predictions. Results in Figure 12 and Figure 13 compare the RBF based SVR (SVR_rbf) with polynomial kernel based SVR (SVR_poly) models for both training data and testing data for stall torque ratio. The root mean squared error (
) and correlation coefficient (
) served to evaluate the differences between the predicted and CFD simulation values for SVR_rbf and SVR_poly. A comparison between the SVR with RBF kernel results and the SVR with polynomial kernel results reveals that SVR with RBF kernel outperforms the SVR with polynomial kernel function in terms of estimation prediction for stall torque ratio.
Figure 12. SVR prediction results for training data of stall torque ratio.
Figure 14 represents the results of the peak efficiency of CFD and predicted data in the preparation stage, while Figure 15 demonstrates the results for the testing stage. Both the SVR with RBF kernel and SVR with polynomial kernel perform well in estimating the peak efficiency. For stall peak efficiency prediction, the SVR with RBF kernel function has very small
during training and this value is slightly larger in testing. The SVM with polynomial kernel have larger
in training data and smaller
in testing data than RBF. Table 4 compares the SVR_rbf and SVR_poly models. It can be concluded that the SVR with RBF kernel is better than the SVR with polynomial kernel at estimating the performance characteristics of a torque converter.
Figure 13. SVR prediction results for testing data of stall torque ratio.
Figure 14. SVR prediction results for training data of peak efficiency.
Figure 15. SVR prediction results for testing data of peak efficiency.
Table 4. SVR performance for prediction of stall torque ratio and peak efficiency.
Method |
Stall torque ratio
|
Peak efficiency
/% |
Training |
Testing |
Training |
Testing |
|
|
|
|
|
|
|
|
SVR_rbf |
0.000194 |
0.998621 |
0.062697 |
0.844790 |
0.002403 |
0.986141 |
0.024338 |
0.900186 |
SVR_poly |
0.014299 |
0.910808 |
0.072086 |
0.764089 |
0.005023 |
0.971390 |
0.017835 |
0.933769 |
6. Conclusions
An attempt is made in this study to investigate the performance of support vector regression (SVR) technique for predicting torque converter performance including stall torque ratio and peak efficiency. The proposed method is verified by comparing the predicted values with CFD values. Two SVR models are investigated: radial basis function (SVR with RBF kernel) and polynomial function (SVR with polynomial kernel). According to RMSE and
parameters, it can be concluded that the SVR with RBF kernel is better than the SVR with polynomial kernel at estimating the performance of torque converter based on the turbine design parameters.
The prediction results demonstrate that SVR method can predict torque converter performance with higher estimation accuracy and generalization ability. Also according to results SVR may be a promising alternative to existing prediction models. It can be seen that the prediction model overcomes the main drawback of artificial neural networks without defining network structure and trapping in the local optimum.