The spatiotemporal evolution of temperature of leaves-on-branchlet carbon based electric double layer capacitors (EDLCs) under imposed constant current was studied using a continuum thermal model. The hot spot aggregated at the tips of graphene petals (GPs), particularly at the high concave surface, at the beginning of the charging step. As the charging proceeded, the overall temperature rose continuously, and the temperature distribution was likely uniform throughout the graphene petals due to an increasingly uniform distribution of ions on GPs surfaces. To elucidate the effects of electrode geometry on the change of temperature, several simple two-dimensional structures were also simulated in the charging step. Concave and planar structures contributed to high temperature change, while a convex structure tended to alleviate the hot spot. An insight into geometric effects on the thermal behavior may lead engineers to develop a new class of nanomaterials for supercapacitors.

Electric double layer capacitors (EDLCs) are an important class of electrical energy storage devices which store energy in the form of electric double layers. The charging mechanism is highly reversible physical adsorption of ions into the porous electrodes, which empower these devices to show a remarkable power performance (15kW/kg) and greater life expectancy (> 1 million cycles). However, they store a small amount of energy (5Wh/kg) when compared with batteries. Optimization of the performance of EDLCs based on porous activated carbons is highly challenging due to complex charging process prevailing in the Nano pores of electrodes. Molecular simulations provide information at the molecular scale which in turn can be used to develop insights that can explain experimental results and design improved EDLCs. The conventional approach to simulate EDLCs places both the electrodes and electrolyte region in a single simulation box. With present day computers, however, this one-box method limits us to system sizes of the order of nanometres whereas the size of a typical EDLC is at least of the order of micrometres. To overcome this system size limitation, a Gibbs-ensemble based Monte Carlo (MC) method was recently developed, where the electrodes are simulated in a separate simulation boxes and each box is subjected to periodic boundary conditions in all the three directions. This allows us to eliminate the electrode-electrolyte interface. The simulation of the bulk electrolyte is avoided through the use of the grand canonical ensemble. The electrode atoms in the electrode are maintained at an equal constant electric potential likewise the case in a pure conductor with the use of the constant voltage ensemble. In this thesis, the Gibbs-ensemble based MC simulations are performed for an EDLC consisting of porous electrodes. The simulations are performed with aqueous electrolytes of type MX and DX2 (where M=Na+, K+; D=Ca+2; X=Cl , F ) for a wide variety of operating conditions. The water is modelled as a continuum background with a dielectric constant value of 30. The electrodes are silicon carbide-derived carbon, whose microstructure generated from reverse MC technique, is used in the simulations. The results from these simulations help us understand the charge storage mechanism, the effect of size and valence of ions on the performance of nonporous carbon based EDLCs when the hydration effects are indignant. The thesis first demonstrates the presence of finite size effects in the simulations performed with the one-box method for KCl electrolyte. The capacitance (ratio of the charged stored on the positive electrode to the voltage applied) values obtained for KCl electrolyte with the one-box method are significantly higher than the corresponding values obtained from the Gibbs-ensemble method. This shows the presence of finite size effects in the one-box method simulations and justices the use of the Gibbs-ensemble based method in our simulations. The fundamental characteristics of aqueous electrolytes in the EDLC are analyzed with the simulation results for KCl electrolyte. In agreement with experiments and modern mean held theory, the capacitance monotonically decreases with voltage (bell-shaped curve) due to overcrowding of ions near the electrode surface. The charge storage mechanism in both the electrodes is mainly a combination of countering (ions oppositely charged to that of the electrode) adsorption and ion exchange, where coins (ions identically charged to that of the electrode) are replaced with countering. However, at higher voltages, the mechanism is predominantly counter ion adsorption because of the scarcity of coins in the electrodes. The mechanism is preferentially more ion exchange for the positive electrode because of its relatively bulky countering, Cl . The shifting of mechanism towards counter ion adsorption at higher voltages and preferential ion exchange process for the positive electrode are in qualitative agreement with the recent experimental results. The constraint of equal electric potential on all the electrode atoms of the amorphous electrode in the simulations resulted in a non-uniform average charge distribution on the electrodes. It shows that the Gibbs-ensemble simulation approach can account for the polarization effects which arises due to a complex topology of the electrodes. In agreement with earlier experiments and simulation studies, the local structure analyses of the electrodes shows that the highly conned ions store charge more efficiently. On the application of voltage difference between the electrodes, the electrolyte ions move towards higher degree of con ned regions of the electrodes indicating the charging process involves local rearrangement and rescuing of electrolyte ions. The thesis also discusses the effect of temperature and bulk concentration on the performance of EDLCs. The Gibbs-ensemble based simulations are performed for the EDLC with varying temperature and bulk concentration for the KCl electrolyte independently. In agreement with the Guo -Chapman theory and experiments, the capacitance decreases with the temperature and increases with the bulk concentration. This is because the concentration of countering in the electrodes decreases with an increase in the temperature but increases with an increase in the bulk concentration. Lastly, the effect of ion size and valency on the performance of EDLCs is analyzed. The capacitance monotonically decreases with voltage (bell-shaped curve) for all the electrolytes, except for NaF, where a maximum is observed at a non-zero finite voltage (camel-shaped curve). The capacitances of NaCl and NaF are greater than that for KCl and KF, respectively. This is because the smaller Na+ ions have more accessibility to narrow con ned regions, where the charge storage efficiency is high. As expected, the capacitance for CaCl2 and CaF2 are highest among their monovalent counterparts, NaCl and KCl; NaF and KF, respectively. This is attributed to the relatively smaller double layer thickness of the bivalent Ca+2 ions. Interestingly, at higher voltages, the capacitance for the bivalent electrolytes approaches the capacitance for the monovalent electrolytes because the concentration of Ca+2 ions in the negative electrode increases sluggishly with voltage due to a strong electrostatic repulsion between Ca+2 ions.