Optimization of Field-scale Surface Dissolution with Thermoelastic Constraints