We consider the use of partial differential equation-based models for energy contracts. In particular, we describe a semi-Lagrangian finite-element method solving the equations that arise in the context of one- and two-factor models. We demonstrate the application of this method to the pricing of various types of swing options, and make use of the results to explore the properties of swing contracts.