New article: A fully variational numerical method for structural topology optimization based on a Cahn-Hilliard model