Numerical models for design and optimization of fuel cell-based power generation systems require an accurate representation of the SOFC stack performance over a wide range of operating conditions. Response surface techniques were used with detailed SOFC stack model results to create a computationally efficient reduced order model (ROM) of the stack that retains desirable information about its internal state, which can be important in system optimization studies to preclude operation under off-design conditions that can adversely impact overall reliability. A tool was developed to randomly sample selected input parameters, obtain stack solutions for these sampled cases, perform regression to obtain response surfaces for stack performance and other metrics of interest, and export the response surfaces as a ROM. The ROM was then successfully implemented in an Aspen Plus® system model for a natural gas fuel cell (NGFC) power system. Additional work to characterize the ROM’s approximation error is ongoing.