En este trabajo, consideramos aproximaciones numéricas del modelo de cristal de campo de fase tipo Cahn-Hilliard y construimos un esquema de elemento finito completamente discreto para él. El esquema es la combinación del método de elementos finitos para la discretización espacial y un método de cuadratización de energía invariante para el avance en el tiempo. No solo es lineal y de segundo orden en el tiempo, sino que también es incondicionalmente estable en energía. Demostramos rigurosamente la estabilidad energética incondicional y llevamos a cabo varios ejemplos numéricos adicionales para demostrar la estabilidad y precisión del esquema desarrollado numéricamente.