Methane hydrate, a potential future energy resource, is known to occur naturally in vast quantities beneath the ocean floor and in permafrost regions. It is important to evaluate how much methane is recoverable from these hydrate reserves This article introduces the theoretical background of HydrateResSim, the National Energy Technology Laboratory (NETL) methane production simulator for hydrate-containing reservoirs, originally developed for NETL by Lawrence Berkeley National Laboratory (LBNL) It describes the mathematical model that governs the dissociation of methane hydrate by depressurization or thermal stimulation of the system. including the transport of multiple temperature-dependent components in multiple phases through a porous medium The model equations are obtained by incorporating the multiphase Darcy's law for gas and liquid into both the mass component balances and the energy conservation equations Two submodels in HydrateResSim for hydrate dissociation arc also considered a kinetic model and a pure thermodynamic model. Contrary to mole traditional reservoir simulations, the set of model unknowns or primary variables in HydrateResSim changes throughout the simulation as a result of the formation or dissociation of ice and hydrate phases during the simulation. The primary variable switch method (PVSM) is used to effectively track these phase changes. The equations are solved by utilizing the implicit time finite-difference method on the grid system, which can properly describe phase appearance or disappearance as well as the boundary conditions The Newton-Raphson method is used to solve the linear equations after discretization and setup of the Jacobian matrix. We report here the application of HydrateResSim to a three-component, four-phase flow system in order to predict the methane produced from a laboratory-scale reservoir. The first results of HydrateResSim code in a peer-reviewed publication are presented in this article The numerical solution was verified against the state-of-the art simulator TOUGH+Hydrate. The model was then used to compare two dissociation theories, kinetic and pure equilibrium Generally, the kinetic model revealed a lower dissociation rate than the equilibrium model. The hydrate dissociation patterns differed significantly when the thermal boundary condition was shifted from adiabatic to constant-temperature. The surface area factor was found to have an important effect on the late of hydrate dissociation for the kinetic model The deviation between the kinetic and equilibrium models was found to Increase with decreasing surface area factor.