Thermal energies deposited by OB stellar clusters in starburst galaxies lead to the formation of galactic superwinds. Multi-wavelength observations of starburst-driven superwinds pointed at complex thermal and ionization structures which cannot adequately be explained by simple adiabatic assumptions. In this study, we perform hydrodynamic simulations of a fluid model coupled to radiative cooling functions, and generate time-dependent non-equilibrium photoionization models to predict physical conditions and ionization structures of superwinds using the maihem atomic and cooling package built on the program flash. Time-dependent ionization states and physical conditions produced by our simulations are used to calculate the emission lines of superwinds for various parameters, which allow us to explore implications of non-equilibrium ionization for starburst regions with potential radiative cooling.