SUBJECT: M.S. Thesis Presentation
BY: Jesus Arias
TIME: Tuesday, April 12, 2019, 9:00 a.m.
PLACE: Love Building, 295
TITLE: Numerical Modeling of a Solid Oxide Fuel Cell for use in Real-Time Simulation and Cyber-Physical Systems
COMMITTEE: Dr. Comas L. Haynes, Co-Chair (GTRI)
Dr. S. Mostafa Ghiaasiaan, Co-Chair (ME)
Dr. Yogendra K. Joshi (ME)
Dr. Aklilu T. G. Giorges (GTRI)


Cyber-physical systems provide a mechanism with which to investigate the physical phenomena and behavior of traditionally cost-prohibitive or otherwise fragile equipment. For the National Energy Technology Laboratory (NETL), this approach resulted in the Hybrid Performance (Hyper) facility which features a gas turbine-SOFC hybrid cycle utilizing real turbomachinery and a simulated SOFC stack. This allows for the investigation of combined cycle performance and control strategies, in an exhaustive manner, both without fear of destroying delicate state-of-the-art fuel cells, and with the full accuracy of real-world turbomachinery.
Issues arose between the transient response of the SOFC model being limited due to the calculation time of the SOFC model taking on average 40 milliseconds to calculate for a given timestep. In order to be able to match the speed of transients from the turbomachinery and likewise better discern transient behavior, it was determined that the SOFC model must be optimized to operate at a sample time of 5 milliseconds. To do this, both the electrochemical algorithm and the thermal algorithm used to simulate the physical behavior of the SOFC are investigated to determine where improvements can be made.
To this end the rootfinding numerical recipes of the electrochemical algorithm are investigated as the complex electrochemistry requires a highly iterative nested dual convergence loop to resolve the voltage-current relationship, and likewise the temporal discretization of the thermal algorithm is modified for the sake of higher accuracy and stability. Ultimately the new electrochemical algorithm featuring higher order rootfinding schemes proves to be efficient enough to reach the sub 5 millisecond target, and when coupled with the new temporal discretization an additional small decrease in calculation time show that a fully implicit, higher order temporal discretization can also successfully be used if desired.