Abstract

Although both battery and super-capacitor are important power sources for hybrid electric vehicles, there is no accurate configuration theory to match the above two kinds of power sources which have significantly different characteristics on energy and power storage for the goal of making good use of their individual features without size wasting. In this paper, a new performance is presented that is used for analysis and optimal design method of battery and super-capacitor for hybrid energy storage system of a parallel hybrid electrical vehicle. In order to achieve optimal design with less consumption, the power-energy function is applied to establish direct mathematical relationship between demand power and the performance. During matching process, firstly, three typical operating conditions are chosen as the basis of design; secondly, the energy and power capacity evaluation methods for the parameters of battery and super-capacitor in hybrid energy storage system are proposed; thirdly, the mass, volume, and cost of the system are optimized simultaneously by using power-energy function. As a result, there are significant advantages on mass, volume, and cost for the hybrid energy storage system with the matching method. Simulation results fit well with the results of analysis, which confirms that the optimized design can meet the demand of hybrid electric vehicle well.

1. Introduction

The energy storage system (ESS) is taken as the auxiliary power source to alleviate load shifting and improve fuel economy and emission on hybrid electric vehicles (HEVs) [13]. Li-ion battery is widely used in these applications owing to its high energy density, proven safety, and competitive cost [4, 5]. Super-capacitor (SC) with long lifecycle and high energy efficiency is evaluated in automotive industry and academia for hybrid energy storage system (HESS) [6], in cooperation with battery and DC/DC converter [79]. According to above literatures, the first potential benefit of HESS is represented by the power loss reduction in the energy storage. In fact, the energy efficiency of SC is higher than that of battery, especially at significant currents. Moreover, SC allows regeneration even when the vehicle is working in critical ambient conditions (i.e., at low temperatures) [1012], in which the battery cannot operate in regenerative mode [13]. The key issue with HESS optimization is how to effectively protect the battery by using SC in order to improve system efficiency and prolong service life of battery [14]. Three main factors should be considered: the HESS topology, battery and SC size, and control strategy (CS) [15].

Regarding the HESS topology, there are three major types: passive, semiactive, and fully active. The passive HESS is the simplest topology because the battery and SC package are connected in parallel and coupled directly to the DC bus [11, 16]. In this topology, the SC works as a low-pass filter for battery. Although no additional electronic components are needed, the filtering effect inherently depends on the resistance of the battery/SC as well as the capacity because of the difference of voltage characteristics between battery and SC. Nevertheless, from the costs and the actual effects, passive HESS is rarely used in practical engineering. On the contrary, because the fully active HESS employs two DC/DC converters and additional control circuit, this topology demands compromise in terms of cost, mass, volume, efficiency, and simplicity [15, 17]. At the same time, the losses in DC/DC should not be ignored [18]. Based on the results of studies [19, 20], the semiactive HESS is favored to be the most popular topology, which only employs one DC/DC converter. It is a good tradeoff between performances, cost, and complexity. In additional, most control strategies can be implemented in this topology. Taking into account the need to charge/discharge frequently and the losses in DC/DC, the topology similar to the one used in [21] with a smaller DC/DC in comparison with the conventional HESS, is regarded as a case study in this paper. The SC provides energies directly to the motor, and battery with DC/DC will participate into work according to the state of SC when necessary.

For optimal energy management, a variety of online CSs for HESS have been empirically designed and discussed in recent literatures, such as the “all or nothing” [1], rule-based [19], filtration based [20], model predictive [22], fuzzy logic CSs [18, 23], efficient CS [24], and minimizing total fuel consumption [25]. Furthermore, under certain driving cycle and configurations, the offline CSs such as dynamic programming (DP) and Pontryagin’s minimum principle algorithms could find the global optimization solution. Santucci et al. [21] propose a DP algorithm including a simplified battery aging model. Vinot et al. [26] use Pontryagin’s minimum principle to reduce the RMS of the battery current.

Besides the control strategy, the configuration is another key factor contributing to the performance of HESS. Schupbach et al. [16] use / ratio as central indices throughout the design, which should be adjusted to match the ratio of the maximum power to the maximum energy of HESS. However, because of the complexity of the active control system and the different features between the batteries and the SC, the / ratio does not accurately describe the power and energy capability of the active-combination energy storage system (ACES) or the fluctuating power demand from the HEV with some arguments [27]. Considering the interaction between the CS and configuration, Miller et al. [20] present a codesign philosophy to get the optimal solution; cost and size ratios are selected as two factors to determine the HESS with maximum energy storage capacity. A set of optimal size ratio under various fixed costs is thus searched, and the cost and size ratio with the largest value of energy is globally searched to be the optimal solution of system design; the philosophy is also used by Hung and Wu [27]. Song et al. [28] utilize the dynamic programming (DP) approach to deal with the integrated optimization problem for deriving the best configuration and energy split strategies of a HESS for an electric city bus. It is shown that the life cycle cost of the HESS initially decreases rapidly with the addition of SC, though the rate of this reduction decreases as the amount of SC increases. Hung and Wu [27] develop a simple but innovative integrated optimization approach (IOA) for deriving the best solutions of component sizing and CSs of a HESS. The simulation results show that 6% of the total energy can be saved in the IOA case compared with the original system in two driving cycles: ECE and UDDS, and two vehicle weights, respectively.

By using DP, the optimal HESS solution in terms of a specific vehicle and driving cycle can be found, which has been proven crucial for saving energy, reducing cost, reaching high overall efficiency, and reducing the peak and root mean square (RMS) values of the battery current to extend its life span [1, 29]. However, this optimization problem is computationally expensive because a lot of possibilities exist when all the factors are alterable. Therefore their use in practical applications is challenging [30]. As a result, only one factor (e.g., the HESS topology or CS) is considered in most published literature [31], and only a few focus on two or more different factors simultaneously [28].

To present a widely applicable and mathematical method to analyze the design process, the PE function and PE density matrix which is extracted from the working data of HESS used in an HEV are proposed in [32]. The minimum matched ACES under optimal power split strategy which can meet the requirements is founded. In the meantime, the efficiency of the system can also be estimated. In the meantime, considering energy loss as the evaluation indices, it comes to the conclusion that the energy loss increases firstly then declines with the increasing trend of SC capacity. However, it only discusses the capacity relation between battery and SC qualitatively, which is not sufficient enough to guide the engineering practice quantitatively.

Overall, considering the existing research either time-consuming and cycle dependent or only suboptimized, there is no efficient and practical applicable method in multiobjective optimization when factors are alterable. This paper considers the factors such as charge/discharge capacity, size, life cycle, and DC/DC efficiency which could optimize the parameter configuration effectively.

The objectives of this paper are the following:(i)to extract the power and energy demands of ESS for a case study parallel HEV along conditions which are not cycle dependent;(ii)to propose appropriate indices (i.e., energy capacity for charging/discharging) as calculation basis of performance parameters;(iii)to present a novel and practical method in multiobjective optimization (i.e., mass, volume, and cost of HESS) by using PE function. Further optimization can be done combined with the research results of [21, 32];(iv)to verify the accuracy of analysis through simulation.

This paper is based on several functional requirements of HESS for HEV. The charge/discharge capacity of HESS which is needed to meet has been proposed through modeling and simulation. According to the characteristic of vehicle operating, the power and energy density of battery and SC are recalculated based on the charge/discharge capacity requirements. Multiobjective optimization of HESS is proposed including mass, volume, and cost.

2. Power Demand of HESS in HEV

In this paper the HEV is a parallel power system, and the configuration is shown in Figure 1. The HESS comprises SC package, battery package, DC/DC converter, engine, ISG motor, electronic controlling clutch, main drive motor (M), and transmission unit. Accessories such as air conditioning and vacuum booster systems are electrically driven. The battery and SC package are combined in parallel through the DC/DC converter to supply power to the main drive motor to support the regular operations of the system. The HESS is required to have the ability to drive the vehicle, supply enough power, and restore the regenerative braking energy efficiently. The parameters of the HEVs are listed in Table 1.

For the high voltage DC bus of the ESS, the design range of operating voltage is 250–400 V. In the case that the available energy of both the batteries and the SC is 50%, the demanded power under the following three conditions should be satisfied by the ESS [1820].

(i) Power up in EV mode is the first condition. In the start phase, the power supplied by single motor is used to start the entire HEV without engine working. When operating in city, public buses have to stop and start repeatedly. Public HEV bus could cancel idle mode by shutting the engine off automatically when the vehicle halts. However, accessories such as the vacuum booster and air conditioner still depend on the ESS to supply the energy. Therefore, the ESS is required to maintain the functionality of those accessories during the short period like idle and then still have the ability to power up the HEV using electricity only. The sample used in this study should have the ability to power up in EV mode after halting for 3 minutes and then speeds up to 20 km/h. The demanded power for the ESS during this process is shown in Figure 2.

(ii) For acceleration, HEV requires the motor to provide high auxiliary power during acceleration, meaning that the ESS must have a certain level of discharging ability. The sample vehicle in this study is required to accelerate from 0 to 50 km/h within 20 seconds. The demanded power for the ESS during this process is shown in Figure 3.

(iii) During regenerative braking, the main drive motor should supply the force for braking, which converts the vehicle’s kinetic energy to electrical energy that would be stored in the ESS. The ESS is required to have a certain level of charging ability. The sample used in this study must be able to fulfill a regenerative brake while traveling speed is reduced from 50 km/h to standstill. The demanded power for the ESS during this process is shown in Figure 4.

3. Evaluation Method for the ESS Performance of the HEV

The performance requirements for the ESS are reflected in two aspects: (i) power capacity; (ii) energy capacity during charging and discharging. The application environment of the HEVs’ ESS is discussed and a method is proposed to calculate the ability of the batteries and SCs on the two performance requirements in this section. The parameters for the Li-ion battery and the SC are listed in Table 2 [12, 32].

HEVs have higher demanded power but lower energy demand for the ESS. Li-ion battery has high energy density but low power density. Therefore, the energy of the Li-ion battery used in HEVs is usually more than required [3336]. In order to achieve good performance characteristics over lifespan, the state of charge (SOC) of the battery should be controlled in the range of 0.4–0.6. So it is generally controlled within that range to prolong the battery life [37].

Based on previous studies [38], there is only mild variation during charging/discharging when the SOC is 0.4–0.6; therefore, it can be treated as a constant. When the SOC is 0.4, the maximum continuous discharging power of the battery is the minimum value. That value is denoted as , which is the limit of the battery’s discharging power during operating. When the SOC is 0.6, the maximum continuous charging power of the battery is the smallest value. That value is denoted as , which is the limit of the battery’s charging power. Based on the charging and discharging performance requirement for ESS, 50% of the battery energy could be set as the balance point for charging and discharging [31]. And the battery SOC is chosen as 0.4–0.6 in where is the energy capacity for charging and discharging (Wh); OCVBAT is the open-circuit voltage (V); and is the battery capacity (Ah). Based on the principle above, the power and energy performance parameters of the battery in HEV (with the cell parameters in Table 2) can be ensured in Table 3.

SCs are used in HEV because of their high power density, efficiency, low internal resistance, and long operational lifespan [39]. The efficiency of SC during charging and discharging is reduced with charging/discharging power increasing. For HEVs, when calculating the power limits during charging/discharging, certain efficiency of the SC should be guaranteed [40, 41]. The efficiency can be calculated using the internal resistance model as shown in where represents efficiency; is the open-circuit voltage; is the current; and is the internal resistance. In this study, the lowest limit for efficiency is set as 0.9 to calculate the current limit of the SC during charging/discharging. This is shown in where and refer to the current limits of the SC during discharging/charging, respectively. Using the terminal voltage and current limits for charging/discharging, the power limit of the SC during charging/discharging can be calculated on the basis of the efficiency considerations for the HEV. This is shown in where and refer to the power limit of the SC during discharging/charging, respectively.

As shown in the equations above, the power limit of the SC in HEVs based on the efficiency consideration is determined by the internal resistance and open-circuit voltage [42]. The former is an intrinsic property, whereas the latter has a linear relationship with electricity. SCs are placed at both ends of the motor in parallel. Their operating range is 250–400 V, with a corresponding range of a single working voltage at 1.67–2.67 V. The open-circuit voltage of the SC can also be approximated at 1.67–2.67 V. In this study, is set at 1.67 V and fed into (4). Then the power limit of the SC under worst operating condition during charging/discharging is obtained. This value is then used as the power limit of the SC in HEV during charging/discharging [6].

Equation (5) is the energy calculation for SC. So the energy capacity of the SC can be calculated using voltage range, as shown in where represents the energy capacity for charging/discharging (Wh). Based on the principle above, the calculated performance of the SC is summarized in Table 4.

The volume, mass, and cost of per unit can be calculated by combining the data from Tables 24. The performance of the battery and SC in terms of power and energy is summarized in Table 5.

4. Demanded Power Based Optimal Design for the Performance of HESS

4.1. HESS Configuration

The HESS is the combination of Li-ion battery and SC. The overall performance is determined by the performance of the individual units, as well as the total number used in series [43]. Optimization of the HESS is to decide the cell quantity used in series which can achieve optimal cost, mass, and volume, while ensuring that the performance of the HESS can meet the demand.

As is shown in Figure 5, the SC package is constituted by SC branches which are arranged in parallel at both ends of the motor, with an operating voltage range of 250–400 V. There are two characteristics for the SC package: (i) the wider the range of the operating voltage is, the greater the energy will be; (ii) for the same quantity units, the larger number of units that are arranged in series, the smaller combined capacity of the SC package. Owing to these two characteristics, 150 single SC units are arranged in series to constitute a SC branch. SC branch is treated as a special indivisible unit. SC branches are grouped in parallel to form the SC package of N units. With a DC/DC converter, the range of operating voltage for a battery unit is 160–260 V. The configuration of (number of battery units) is flexible and can be arranged in either series or parallel according to the actual requirement.

4.2. Performance Analysis of the HESS

The performance of the HESS refers to its ability to satisfy the demanded power [44]. The performance is determined by the individual performance parameters, including energy capacity and charging/discharging power limits of the battery and SC package. The following three conditions are used to determine whether the HEES can satisfy the demanded power shown below.

(i) (demanded power) must be less than the sum of the limits of (battery power) and (capacitor power), which is(ii) SCs must satisfy any demanded power that exceeds the power limits of batteries. However, the integral of (demanded power from SCs) cannot exceed (energy capacity of SCs), which is(iii) The energy required for the entire charging/discharging process should be less than the sum of the energy capacity of (batteries) and (SCs), which is

To analyze the effect of parameters for the SC and battery package on the performance of the HESS more intuitively, the concept of PE function is introduced. The function is defined as where represents the energy (Wh) and represents the power (kW). in the function represents the integral value of the power versus time whose demand is greater than . This is the energy value of this power , as shown in where is the power-demand function in terms of time (kW). This function is known as the PE function. Equations (7)–(9) can be rewritten using this function as follows:

With regard to a particular HESS, it can be observed that the PE function can be used to accurately determine whether the demanded power for charging or discharging process can be met. The demanded power for the HESS under the condition of acceleration is used as an example and shown in Figure 6(a). The corresponding PE function is shown in Figure 6(b).

Based on (12), when the points (), (), and () shown in Figure 6(b) are not below the PE function curve, the HESS would be able to meet the demanded power. In the PE function figure, the straight line function passing through points () and () is shown below:

It is obvious that the gradient of the line is determined by the energy capacity and discharge limit of the battery. Therefore, this is known as the discharge performance line of the batteries. The ratio between the discharge limit and the energy capacity is (discharge rate of the batteries). This is shown in

Hence, the function of the discharge performance line of the battery can be shown as follows:

Similarly, the straight line through points () and (, 0) is known as the discharge performance line of the SC. The ratio between the discharge limit and the energy capacity is (discharge rate of the SCs). This is shown in

Therefore, the function of the discharge performance line of the SC can be shown as follows:

Point () is the intersection of the performance line of the battery with the -axis, and with the -axis it is (), respectively. Point () is the intersection between the performance lines of the batteries and SCs, which is also known as the parameter point for the performance of the HESS.

From the performance lines of the batteries and SCs, it can be concluded that the HESS has to satisfy three determining conditions to meet the arbitrary demanded power:(i)The intersection between the performance line of the battery and -axis must be above the intersection between the PE function of demanded power and the -axis, as shown in Figure 7(a).(ii)The intersection between the performance line of the SC and -axis should be at the right side of the intersection between the PE function of demanded power and -axis, as shown in Figure 7(b).(iii)The intersection between the performance lines of the batteries and SCs should be above the PE function of the demanded power, as shown in Figure 7(c).

As shown in Figure 7(c), each intersection between the performance lines of the batteries and SCs (i.e., each parameter point for the performance of the HESS) corresponds to a possible design proposal in which the HESS is able to meet the demanded power when discharging. Theoretically, Figure 7(c) can be used to obtain all possible methods of setting up the HESS.

Applying the same method, the demanded power and PE function required in EV mode (Figure 8) can be used to determine the range of parameter points for the performance of the HESS when discharging.

The range of parameter points for the performance of HESS when discharging under the EV mode and acceleration mode can be determined by combining the calculation results from acceleration and EV mode. In Figure 9, the region of the right side of the thick black dotted line represents the range of the parameter point () for the performances of the HESS that meet the demanded power when discharging.

With the same method, the range of parameter points () for the performance of the HESS that meet the demanded power during charging can be calculated using Figure 10, which shows the demanded power during charging and the corresponding PE function for the regenerative braking process.

The above method is similarly used to, respectively, derive the range of parameter points () and () for the performance of the HESS during discharging/charging. Furthermore, (18) are used to indicate the relationship between (i) the power limits of the batteries during charging/discharging and (ii) their energy capacities. This information is then used to derive the range of energy capacities and (for batteries and SCs, resp.) of the HESS that can simultaneously meet the demanded power during charging/discharging. This is shown in Figure 11:

In Figure 11, the blue and green lines mark the selection range of () which is determined by the parameter point for discharging/charging performance. Based on these two lines, the selection range of () that simultaneously meets the demand during charging/discharging is determined and is indicated by the white area at the right side of the thick dotted line ABCD.

It should be noted that, in Figure 11, the lines for the charging/discharging of the battery coincided after (18) are used to perform coordinate transformation. At point A, is 0, meaning that the SC is acting as the ESS of the HEV. At point D, is 0, meaning that the battery is acting as the ESS of the HEV. Any point within the white area, such as point E, would be able to meet the demanded power for both charging/discharging. There would always be a point F on line ABCD, which has the same abscissa but a smaller ordinate than that of point E. Therefore, point F could be used as the design reference for the HESS that is able to meet the demanded power; however, it simultaneously requires SCs with much smaller capacity. This illustrates that, based on the volume, mass, and cost considerations, the optimal () should be located on line ABCD. The cost, mass, and volume of the HESS when using batteries and SCs of various capacities are then obtained by using the respective energy price, energy mass, and energy volume ratios for the battery and SC as stated in Table 5. This is shown in

The results of the mass, volume, and cost of the HESS using battery of different capacity on line ABCD are calculated and are shown in Figure 12.

It can be observed from Figure 12 that when the battery capacity range is 0.6–1 kWh, the weight, volume, and cost of the system are almost minimum values. Compared with the use of either battery or SC as the ESS of the HEVs, the mass, volume, and cost of the ESS are reduced by 30%–50%. The point F (0.8 kWh, 130 Wh) is used as the basis for optimization design of parameters of the HESS. This method is described in the next subsection.

4.3. Optimal Design of the HESS

The energy capacity of a single battery unit is 2.6 Wh (Table 3). The battery package requires 0.8 kWh of energy capacity; therefore, 308 battery units are required. Based on the principle that the voltage of the battery package does not exceed the SC package, a configuration of 77 joints of units in series and four strings in parallel is adopted to form the battery package. The open-circuit voltage of the battery package is 246.4 V.

The energy capacity of a single SC unit is 0.45 Wh (Table 4). The SC package requires 130 Wh of energy capacity; therefore, 289 SC units are required. Using the same principle stated above, the design of the SC package uses 150 joints of units in series. Two strings are laid in parallel, meaning that the SC package comprises 300 units.

The cost, mass, price, and other results of the HESS based on the configuration stated are listed in Table 6. The HESS is designed using the same method and using the following points as the bases in Figure 11: A (0 kWh, 508 Wh); B (0.29 kWh, 238 Wh); C (1.6 kWh, 50 Wh); and D (3.55 kWh, 0 Wh). After optimization, there are significant advantages (30%–50%) in terms of the mass, volume, and cost for the HESS.

5. Simulation Test to Validate the Performance

The parameters of the battery and SC package are calculated using the single unit parameters of the SCs and batteries, which are listed in Table 7.

The internal resistance model for the battery and SC package is established using parameters stated in Table 7. The relationships between the DC/DC input and output voltages and currents are used to establish the DC/DC transformer model, which is shown in where and represent the DC/DC output current and voltage, respectively; and represent the DC/DC input current and voltage respectively; and represents the DC/DC transfer efficiency. The corresponding DC/DC control strategy model is established simultaneously. This model, along with the models for the DC/DC converter, battery package, and SC package, forms the model for the HESS. This is shown in Figure 13.

The load model is established on the basis of the demanded power under the above three operating conditions (EV, acceleration, and regenerative brake mode). The load model is then used to validate the performance of the HESS through simulation. The results are shown in Figure 14.

During the simulation, a DC/DC converter is used to control the current of battery in the range of −80 to 160 A. This is to ensure that the battery power could be controlled within its power limit range. Based on the demanded power for the two operating conditions: acceleration and powering up in EV mode, the voltage of the HESS at the load end decreases from the initial 340.5 V to 259 V and 279 V, respectively, whereas the battery SOC decreases from the initial 50% to 48.1% and 44.7%, respectively (Figures 14(a) and 14(b)). For regenerative braking, the voltage at the load end increases from 340.5 V to 393 V, whereas the battery SOC increases from 50% to 51% (Figure 14(c)).

The voltage of the HESS does not exceed 250–400 V during the simulation, while the battery SOC and power are within their limits. The simulation shows that the proposed method based on the PE function could meet the design requirement.

Further, the simulation results of the energy capacity for the SC and battery package during the operation can be verified with each other in Figure 11. The demand power for HESS () under the three operating conditions (acceleration, braking, and powering up in EV mode) is labeled as C1,  C2, and C3, respectively (Figure 15).

When point F (0.8 kWh, 135 Wh) is used to design the HESS, the demanded power for SC package under the three operating conditions could be observed on the position of the intersections between C1,  C2,  and C3 and the line FG, respectively. The length from the intersections between C1/C2/C3 and FG to point G occupies 91%, 92%, and 48% of the total length of FG, respectively. This shows that, to meet the demanded power during acceleration, braking, and powering up in EV mode, the SC package has to supply 91%, 92%, and 44% of the energy capacity, respectively.

Similarly, the demanded power for the battery package under the three operating conditions can be observed from the intersections between C1, C2, and C3 and the line HI, respectively. The length from the intersections between C1, C2, and C3 and HI to point I is 19%, 10%, and 51% of the total length of HI, respectively. This shows that, to meet the demanded power during acceleration, braking, and powering up in EV mode, the battery package has to supply 19%, 10%, and 51% of the energy capacity respectively.

The results of the analysis fit well with the simulation results, thereby validating the proposed design method for the HESS.

6. Conclusion

This paper adopts the concept of the PE function to establish quantitatively mathematical relationship between the battery/SC parameters and HEV demanded power. Based on the three modes (EV, acceleration, and regenerative braking), power and energy capacity are used as design gist for HESS configuration. Taking a hybrid bus as model, this paper calculates the power limit of each power source in HEVs with efficiency consideration. The energy capacity of a single battery unit is 2.6 Wh (Table 3). If battery package requires 0.8 kWh for energy capacity, 308 battery units are required. According to the mathematical methodology, transfer the time-power function into PE function. With PE function and relationship mentioned above, this paper analyzes the demanded power, mass, volume, and cost of HESS and formats mathematical model which contributes to the optimal HESS configuration. F (0.8 kWh,130 Wh) could get the best matching on aspects of mass, volume, and cost compared to A (0 kWh, 508 Wh), B (0.29 kWh, 238), C (1.6 kWh, 50 Wh), and D (3.55 kWh, 0 Wh). Then HESS can be designed for optimization based on analysis results. Battery package constitutes 77 joints of Li-ion battery cells in series and four strings in parallel, while SC package contains 150 joints of capacitor units in series and two strings in parallel. Finally, simulation verifies the optimized HESS that has met the demanded power. F (0.8 kWh, 135 Wh) could meet the power requirements in three driving modes, which is consistent with the foregoing analysis results. Furthermore, the calculation method is the HESS configuration theory which could guide the practice engineering of HESS design for new energy vehicles.

Future work will focus on the degradation model of battery and verification of real vehicle performance.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (51107052).