Combining maximal regularity and energy estimates for time   discretizations of quasilinear parabolic equations
  We analyze fully implicit and linearly implicit backward difference formula (BDF) methods for quasilinear parabolic equations, without making any assumptions on the growth or decay of the coefficient functions. We combine maximal parabolic regularity and energy estimates to derive optimal-order error bounds for the time-discrete approximation to the solution and its gradient in the maximum norm and energy norm. 
Georgios Akrivis, Buyang Li, Christian Lubich
math.NA