A new method for the determination of clay swelling thermodynamics from computer simulation is discussed and evaluated. This method allows for the determination of temperature, pressure, and water chemical potential dependence of clay swelling from simulations at a single thermodynamic state point. The temperature dependence and pressure dependence of clay swelling are shown to be directly related to the composite system entropy and volume change, respectively, that accompany swelling. Expressions for the chemical potential dependence of clay swelling are used to determine constant pressure layer spacing and adsorption isotherms, quantities that are well suited for comparison with experimental measurements. This method is evaluated through grand isoshear ensemble simulations of Na-montmorillonite, a prototypical swelling clay. Approximations associated with all expressions are discussed with explicit calculations used to demonstrate their regimes of validity.