We introduce a mathematical model, which describes the charge inversion phenomena in systems with a charged wall or boundary. This model may prove helpful in understanding semiconductor devices, ion channels, and electrochemical systems like batteries that depend on complex distributions of charge for their function. The mathematical model is derived using the energy variational approach that takes into account ion difusion, electrostatics, finite size effects, and specific boundary behavior. In ion dynamic theory, a wellknown system of equations is the Poisson-Nernst-Planck (PNP) equation that includes entropic and electrostatic energy. The PNP type of equation can also be derived by the energy variational approach. However, the PNP equations have not produced the charge inversion/layering in charged wall situations presumably because the conventional PNP does not include the finite size of ions and other physical features needed to create the charge inversion. In this paper, we investigate the key features needed to produce the charge inversion phenomena using a mathematical model, the energy variational approach. One of the key features is a finite size (finite volume) effect, which is an unavoidable property of ions important for their dynamics on small scales. The other is an interfacial constraint to capture the spatial variation of electroneutrality in systems with charged walls. The interfacial constraint is established by the difusive interface approach that approximately describes the boundary effect produced by the charged wall. The energy variational approach gives us a mathematically self-consistent way to introduce the interfacial constraint. We mainly discuss those two key features in this paper. Employing the energy variational approach, we derive a non-local partial diferential equation with a total energy consisting of the entropic energy, electrostatic energy, repulsion energy representing the excluded volume effect, and the contribution of an interfacial constraint related to overall electroneutrality between bulk/bath and wall. The resulting mathematical model produces the charge inversion phenomena near charged walls. We compare the computational results of the mathematical model to those of Monte-Carlo computations.