Searching for an energy saving numerical method

Hello everybody,

i have the following problem and i would really appreciate your help.

I am searching for an energy saving numerical method for pde’s (first order).
The Newmark-Beta-Method would be really great but unfortunately i don’t have a mass matrix after the discretisation.

Thank you!

Greetings