  We present a natural framework for constructing energy-stable time
discretization schemes. By leveraging the Onsager principle, we demonstrate its
efficacy in formulating partial differential equation models for diverse
gradient flow systems. Furthermore, this principle provides a robust basis for
developing numerical schemes that uphold crucial physical properties. Within
this framework, several widely used schemes emerge naturally, showing its
versatility and applicability.
