Non-equilibrium thermodynamics of the driven resonant level model is studied using numerical simulations based on the driven Liouville von-Neumann formalism. The approach is first validated against recently obtained analytical results for quasi-static level shifts and the corresponding first order corrections. The numerical approach is then used to study far-from-equilibrium thermodynamic properties of the system under finite level shift rates. The proposed methodology allows the study of unexplored non-equilibrium thermodynamic regimes in open quantum systems.