We present a new general methodology capable of modeling chemical reactions at solid-liquid interfaces called CECILIA (combined embedded cluster at the interface with liquid approach). The main idea is to combine the embedded duster molecular orbital or density functional methods for describing interactions at the surface of a solid with the dielectric continuum approach for modeling a liquid. More details are given on how to apply this methodology to model processes at the ionic solid-water interface. Geometries and adsorption binding energies of H2O, NaCl, Na+, and Cl- al the NaCl(001)-water interface are calculated using this approach and compared with those at the NaCl(001)-vacuum interface.