Energetic Stable Discretization for Non-Isothermal Electrokinetics Model

We propose an edge averaged finite element(EAFE) discretization to solve the Heat-PNP (Poisson-Nernst-Planck) equations approximately. Our method enforces positivity of the computed charged density functions and temperature function. Also the thermodynamic consistent discrete energy estimate which resembles the thermodynamic second law of the Heat-PNP system is prescribed. Numerical examples are provided.