Mitochondrial cristae are dynamic invaginations of the inner membrane and play a key role in its metabolic capacity to produce ATP. Structural alterations caused by either genetic abnormalities or detrimental environmental factors impede mitochondrial metabolic fluxes and lead to a decrease in their ability to meet metabolic energy requirements. While some of the key proteins associated with mitochondrial cristae are known, very little is known about how the inner membrane dynamics are involved in energy metabolism. In this study, we present a computational strategy to understand how cristae are formed using a phase-based separation approach of both the inner membrane space and matrix space, which are explicitly modeled using the Cahn–Hilliard equation. We show that cristae are formed as a consequence of minimizing an energy function associated with phase interactions which are subject to geometric boundary constraints. We then extended the model to explore how the presence of calcium phosphate granules, entities that form in calcium overload conditions, exert a devastating inner membrane remodeling response that reduces the capacity for mitochondria to produce ATP. This modeling approach can be extended to include arbitrary geometrical constraints, the spatial heterogeneity of enzymes, and electrostatic effects to mechanize the impact of ultrastructural changes on energy metabolism.