Investigation on performance of SOFC in hydrocarbon fuel S. Senthil Kumar, Sharad Chauhan, B. Shriprakash, S. T. Aruna CSIR National Aerospace Laboratories Bangalore-560017 Presentation at COMSOL 2016, Bangalore, Oct 20 21 2016
Solid Oxide Fuel Cell (SOFC) At the anode: H 2 + 2O 2 2H 2 O + 4e. At the cathode: O 2 + 4e 2O 2 The overall cell reaction: O 2 + 2H 2 2H 2 O Roles of Electrolyte: Oxygen ion conduction Physically seperates the fuel from oxidant Roles of Electrode: Hosts triple phase boundary to support electrochemical reactions Provides path for O 2- ions/electrons Provides channels for gas diffusion Gives mechanical support to system
Hydrocarbon Compatible SOFC Most promising clean energy source Expansion of potential fuel range Eliminates fuel preprocessing Increases overall system efficiency Bottle Neck: Carburization of SOFC anodes
Genesis of Present work State of art anode cermet is prone to carburization in hydrocarbon fuel An optimum feed/ operating conditions have to be established to minimize carbon deposition. In the typical SOFC anode, carbon deposition is an unsteady state progress Analysis of unsteady variation of porosity and catalytic activity of SOFC anode would be helpful in establishing the durability of SOFC in given fuel feed conditions Objectives To simulate the performance of SOFC with the reformed feed (CH 4 +Steam) using COMSOL To predict the performance degradation over the long duration (3000 h)
Physics Involved Secondary current Distribution: Determines current profile. Accounts for the effect of the electrode kinetics and losses due to resistance Transport of Concentrated Species: Determines species flux across electrode. Involves flow of species across the porous electrodes via diffusion and transport of oxide ion Free and Porous Media Flow: Determines flow profile. Accounts for flow in channel and porous media
Reactions Involved Methane-steam reforming reaction rate: R r = a(k f,r p CH4 p H2 O k b,r p CO p 3 H2 ) CO water-gas shift reaction rate: R s = a(k f,s p CO p H2 O k b,s p CO2 p H2 ) Methane cracking reaction rate: Boudouard reaction rate: R C = a 1 3600M C k f,c p CH4 p 2 H 2 K p,c 1 + k H p H2 2 R B = a 1 k f,bk CO p CO 3600M C 1 + K CO p CO + 1 K B p CO2 p CO 1 K CO2 K CO p CO2 p CO 2
Mathematical models for Parameter study Porosity variation rate: Where r C = (R C + R B ) Catalyst activity variation rate:. dε dt = εr CM C ρ C da dt = k ar C 2 c C a Permeability: κ = κ 0 ε ε 0 3.55 Anode Exchange current densities : RT 11 i 0,H2 = 2.1 10 F p H2 O 1.78 10 9 p H2 0.266 exp 1.2 105 RT RT 11 i 0,CO = 0.84 10 F p CO2 1.63 10 9 p CO 0.266 exp 1.2 105 RT
Mole Fraction Efficiency in % Effect of Steam 0.3 0.25 0.2 S/C=0.1 H2 98 96 94 92 800K 900K 1000K 0.15 0.1 0.05 CH4 C(s) CO2 CO 90 88 86 1100K 1200K 0 400 600 800 1000 1200 1400 Temperature (K) 84 82 1 2 3 4 5 S/C ratio Complete utilization of methane after >1023 K Efficiency decreases with increasing S/C Ƞ decreases steeply with increasing temperature
Fuel Utilization Hydrogen mole fraction in anode Carbon monoxide mole fraction in anode S/C =1, V cell =0.7 V and T=1073 K Good utilization of H 2 in fuel CO consumption was sluggish
Effect of fuel velocity Fuel utilization decreases with increasing velocity CO consumption was affected more than H 2 with increasing velocity Fuel velocity of 0.4 m/s had good depletion profile and reasonable current density
Carbon activity H 2 reduction reaction Boudouard reaction Carbon activity was <1 up to the fuel velocity of 0.9 m/s At higher velocity, the difference in carbon activity between above reactions was much pronounced
C in mol/m 3 Porosity Permeability (m 2 ) a Transient Studies Carbon concentration vs Time Catalytic activity vs Time 60000 40000 20000 0 0 1000 2000 3000 4000 Time in hours 1.5 1 0.5 0 0 1000 2000 3000 4000 Time in hours 0.42 0.4 0.38 0.36 0.34 0.32 Porosity vs Time 0 1000 2000 3000 4000 Time in hours After 120 days of Operation Permeability vs Time 1.50E-14 1.00E-14 5.00E-15 0.00E+00 0 1000 2000 3000 4000 Time in hours 15 % reduction in porosity 50 % reduced permeability
Carbon deposition vs. Time 720 h 1440 h 2160 h 2880 h Carbon builds-up quickly during the initial days Carbon deposit nears saturation in long duration
Power (W/cm 2 ) i (A/cm 2 ) V_cell Performance Degradation 1.8 Current vs Time 1.5 Polarization curve 1.6 1 1.4 1.2 0.5 1 0 1000 2000 3000 4000 Time in hrs 0 0 0.5 1 1.5 2 i (A/cm 2 ) 1.2 1 0.8 0.6 0.4 0.2 0 Power curve 0 0.5 1 1.5 2 i (A/cm 2 ) 0 h 1440 h 2880 h 29% drop in performance after 3000 h 0 h 1440 h 2880 h Concentration polarization become much pronounced with time
Conclusions S/C ratio of 1 was found to be suitable for SOFC operating temperature of 1073 K Fuel velocity of 0.4 m/s found to be suitable in the perspective of fuel utilization and carbon activity Transient studies showed 29 % drop in performances over the period of 3000 h due to carbon deposition
Acknowledgements The Director, NAL Head, SED Division CSIR-Network project (CERMESA )-Funding M/s. Carborundum Universal Ltd. - Funding THANK YOU