
We introduce a novel constructive approach to define time evolution of
critical points of an energy functional. Our procedure, which is different from
other more established approaches based on viscosity approximations in infinite
dimension, is prone to efficient and consistent numerical implementations, and
allows for an existence proof under very general assumptions. We consider in
particular rather nonsmooth and nonconvex energy functionals, provided the
domain of the energy is finite dimensional. Nevertheless, in the infinite
dimensional case study of a cohesive fracture model, we prove a consistency
theorem of a discretetocontinuum limit. We show that a quasistatic evolution
can be indeed recovered as a limit of evolutions of critical points of finite
dimensional discretizations of the energy, constructed according to our scheme.
To illustrate the results, we provide several numerical experiments both in one
and two dimensions. These agree with the crack initiation criterion, which
states that a fracture appears only when the stress overcomes a certain
threshold, depending on the material.

We introduce a new machinery to study the large time behavior for general
classes of HamiltonJacobi type equations, which include degenerate parabolic
equations and weakly coupled systems. We establish the convergence results by
using the nonlinear adjoint method and identifying new long time averaging
effects. These methods are robust and can easily be adapted to study the large
time behavior of related problems.

Characterization results for equality cases and for rigidity of equality
cases in Steiner's perimeter inequality are presented. (By rigidity, we mean
the situation when all equality cases are vertical translations of the
Steiner's symmetral under consideration.) We achieve this through the
introduction of a suitable measuretheoretic notion of connectedness and a fine
analysis of barycenter functions for sets of finite perimeter having segments
as orthogonal sections with respect to an hyperplane.

We consider a numerical scheme for the one dimensional time dependent
HamiltonJacobi equation in the periodic setting. We present a new and simple
proof of the rate of convergence of the approximations based on the adjoint
method recently introduced by Evans.

We provide a geometric characterization of rigidity of equality cases in
Ehrhard's symmetrization inequality for Gaussian perimeter. This condition is
formulated in terms of a new measuretheoretic notion of connectedness for
Borel sets, inspired by Federer's definition of indecomposable current.

The adjoint method, recently introduced by Evans, is used to study obstacle
problems, weakly coupled systems, cell problems for weakly coupled systems of
HamiltonJacobi equations, and weakly coupled systems of obstacle type. In
particular, new results about the speed of convergence of common approximation
procedures are derived.

The adjoint method introduced in [Eva] and [Tra] is used, to construct
analogs to the AubryMather measures for non convex Hamiltonians. More
precisely, a general construction of probability measures, that in the convex
setting agree with Mather measures, is provided. These measures may fail to be
invariant under the Hamiltonian flow and a dissipation arises, which is
described by a positive semidefinite matrix of Borel measures. However, in the
important case of uniformly quasiconvex Hamiltonians the dissipation vanishes,
and as a consequence the invariance is guaranteed.

A new necessary minimality condition for the MumfordShah functional is
derived by means of second order variations. It is expressed in terms of a sign
condition for a nonlocal quadratic form on $H^1_0(\Gamma)$, $\Gamma$ being a
submanifold of the regular part of the discontinuity set of the critical point.
Two equivalent formulations are provided: one in terms of the first eigenvalue
of a suitable compact operator, the other involving a sort of nonlocal capacity
of $\Gamma$. A sufficient condition for minimality is also deduced. Finally, an
explicit example is discussed, where a complete characterization of the domains
where the second variation is nonnegative can be given.