Provided is a method for optimization of fuel cells operating conditions using a hybrid model, and more particularly, a method for optimization of fuel cells operating conditions using a hybrid model which generates a life prediction model determined by time and temperature based on a theoretical performance model and an empirical durability model and estimates an optimal operation temperature in a target life based on the life prediction model.
Please help me write a proper abstract based on the patent claims. CLAIM: 1. A method for optimization of fuel cells operating conditions using a hybrid model, comprising: generating a theoretical performance model by obtaining a potential difference depending on a current density of a cell based on a thermodynamics reversible voltage, an activity loss, a resistance loss, and a concentration loss of a high temperature proton-exchange membrane fuel cell (PEMFC) which is based on a polybenzimidazole (PBI) membrane with which phosphoric acid is doped (S10); generating an empirical durability model by performing a long-term durability test which is performed for a predetermined period under experimental conditions (operating conditions) specifying the high temperature PEMFC based on the PBI membrane with which the phosphoric acid is doped, so as to predict the reduction in cell voltage over an operation time at a predetermined temperature (S20); generating a life prediction model determined by time and temperature based on the performance model generated in the generating of the theoretical performance model (S10) and the durability model generated in the generating of the empirical durability model (S20) (S30); and estimating an optimal operation temperature at a target life based on the life prediction model generated in the generating of the life prediction model (S30) (S40). 2. The method of claim 1, wherein in the generating of the performance model (S10), the potential difference depending on a current density of the cell is obtained by subtracting the activity loss (ηact), the resistance loss and the concentration loss (ηconc) from the thermodynamics reversible voltage (Erev) based on the following Equation. Ecell=Erev−ηact−ηohom−ηconc 3. The method of claim 2, wherein the thermodynamics reversible voltage is calculated based on the following Equation, E ref = - Δ   g rnx ref n   F (In the above Equation, Δgrnxref represents Gibbs' free energy and F represents a Faraday constant) the thermodynamics reversible voltage at a given temperature T is calculated based on the following Equation by introducing a change in entropy depending on temperature, E T = E ref + Δ   s n   F  ( T - T ref ) (In the above Equation, Δs represents the change in entropy depending on the temperature) when a concentration of chemical species is specified, the thermodynamics reversible voltage is calculated based on the following Equation by introducing activity, and E rev = E ref - RT n   F  ln   a H 2  O a H 2  a O 2 0.5 (In the above Equation, R represents an abnormal gas constant and a represents activity) at the given temperature T and concentration c, the thermodynamics reversible voltage is calculated based on the following Equation. E rev = E ref + Δ   s n   F  ( T - T ref ) - RT n   F  ln   a H 2  O a H 2  a O 2  0.5 4. The method of claim 2, wherein the activity loss is calculated by the following Equation. j gross = j 0  (  α   n   F   η c / RT -  - ( 1 - α )  n   F   η c / RT ) j gross = j 0   α   n   F   η c / RT η c = RT α   n   F  ln   j + j leak j 0 j 0 = j 0 ref  S Pt eff  ( c O 2 c o 2 ref )  γ    Δ   b act RT  ( 1 - T T ref ) (In the above Equation, α represents a transfer coefficient, a total of current density gross jgross includes an operating current density j, an internal current, and a loss current jleak occurring due to a leakage of reactants, by being focused on a high current density area in which a reverse reaction is disregarded, a second exponential term of the right terms representing the reverse reaction is disregarded, an exchange current density j0ref is corrected depending on a concentration and a temperature value, a reference exchange current density j0ref is measured based on a reference temperature Tref and a reference oxygen concentration CO2ref, SPteff represents an effective platinum surface area, CO2 represents an oxygen concentration in the catalyst layer, and γ represents a reaction order, j0reff and SPteff are estimated based on the experimental data, and an activity barrier energy Δbact uses 7.24·104 Jmol−1.) 5. The method of claim 2, wherein the resistance loss is represented by the following Equation. η okm = j   δ m κ m κ m = AB T   - b act RT b act = - 619.6  DL + 21750 A = 168  DL 2 - 6324  DL 2 + 65750  DL + 8460 B = 1 + ( 0.1432  T - 56.89 )  RH eff (In the above Equation, δm and κm each represent the thickness and ion conductivity of the membrane, the ion conductivity of the PBI membrane relies on a doping level (DL) of the phosphoric acid, relative humidity (RH), and temperature T, temperature dependency of ion conductivity relies on an Arrhenius law, A and B are pre-exponential coefficients describing an effect of the phosphoric acid doping level and the relative humidity, respectively, bact represents the activation energy of the membrane. 6. The method of claim 1, wherein the generating of the durability model (S20) includes: generating a draft model predicting a reduction in cell voltage over the operation time at the predetermined various temperatures while ruling out an effect of the operation time (S21); estimating each coefficient of the draft model over various times generated in the generating of the draft model (S21) based on a function depending on the temperature (S22); and generating a third-order non-linear function depending on the operation time and temperature using each coefficient estimated in the estimating of the coefficient (S22) (S23). 7. The method of claim 6, wherein for each of the predetermined various temperatures, the draft model in the generating of the draft model (S21) predicts the reduction in cell voltage over the operation time using the draft model corresponding to the following Equation, and VD(t)=x1t3−x2t2+x3t+x4 (The draft models are configured of four coefficients of third-order, second-order, first-order coefficients, and a constant) the final model is to fit the four coefficients x1, x2, x3, and x4 of the draft model depending on the operation temperature. 8. The method of claim 7, wherein in the generating of the draft model (S21), a pretreatment process of converting each coefficient into a log value prior to fitting is performed to improve a prediction ability of the draft model. 9. The method of claim 1, wherein in the estimating of the operating condition (S40), when the target life is maximal, the operation temperature is kept at 150 to 160° C.