We study an abstract second order inclusion involving two nonlinear single-valued operators and a nonlinear multi-valued term. Our goal is to establish the existence of solutions to the problem by applying numerical scheme based on time discretization. We show that the sequence of approximate solution converges weakly to a solution of the exact problem. We apply our abstract result to a dynamic, second-order-in-time differential inclusion involving a Clarke subdifferential of a locally Lipschitz, possibly non-convex and non-smooth potential. In the two presented examples the Clarke subdifferential appears either in a source term or in a boundary term.

In this paper we propose a time discretization of a system of two parabolic equationsdescribing diffusion-driven atom rearrangement in crystalline matter. The equationsexpress the balances of microforces and microenergy; the two phase fields are the orderparameter and the chemical potential. The initial and boundary-value problem for theevolutionary system is known to be well posed. Convergence of the discrete scheme to thesolution of the continuous problem is proved by a careful development of uniformestimates, by weak compactness and a suitable treatment of nonlinearities. Moreover, forthe difference of discrete and continuous solutions we prove an error estimate of orderone with respect to the time step.