1 Introduction
In hamiltonian mechanics a physical system is described by a state vector
and a momentum vector . The evolution in time of the system is governed by a hamiltonian function , via the equations:(1) 
where , denote derivatives with respect to time. The evolution is reversible.
We propose a dissipative modification of these equations, which extends the work from [2], [3], [4].
Definition 1.1
Given a hamiltonian and a likelihood function
a physical system described by the state vector and the momentum vector evolves such that the vector defined by:
(2) 
called the difference (or gap) from hamiltonian evolution driven by , maximizes the likelihood:
(3) 
with the notations , .
The likelihood function cannot be taken arbitrarily, otherwise the proposed modification would be too vague. We propose the following constraints on the likelihood function.
First we define the information content function associated to the likelihood as
(4) 
with the convention that .
Definition 1.2
The likelihood and the associated information content satisfy:

The information content function (4) is convex in each of the 2nd and 3rd variables and it has the needed degree of smoothness required (for example it is lower semicontinuous with respect to the relevant topologies on and )

for any , the following maxima exist
and they are either or .
The equation (3) can be rephrased as: given the hamiltonian and the information content function
, the physical system evolves such that at any moment it minimizes the information content of the gap from a hamiltonian evolution. Indeed, an evolution of the system is a curve
with the property that it minimizes the information content gap functional:(5) 
among all admissible evolution curves .
Pure Hamiltonian evolution.
Let’s pick the information content (4) to be:
This corresponds to a likelihood function:
The maximization of the likelihood (3) implies that the gap vector , therefore the evolution equations (2) reduce to the pure Hamiltonian evolution equations (1).
This example is trivial, we need a method to construct more interesting ones. One such method is based on the following observation, adapted from [6], section 2. We use the notations explained in section 3, in particular we use the duality and we suppose that we have on a topology compatible with it, so that the information content function is lower semicontinuous (lsc).
Proposition 1.3
Bipotentials were introduced in [11] as a convex analysis notion which is well adapted for applications to nonassociated constitutive laws. Bipotentials were used in soil mechanics, plasticity, damage or friction. For the theory of bipotentials see the review paper [5].
By concentrating our attention to the function , instead of the information content , we can build a host of examples. Indeed, for any lsc and convex function
the associated function
(8) 
satisfies the conditions (a), (b) from Proposition 1.3 , where is the polar, or Fenchel conjugate of . Indeed, the mentioned condition (b) is just a reformulation of the Fenchel inequality [11].
2 Examples
Viscosity, Rayleigh dissipation.
In particular, let’s pick
where is a convex, lsc function. A straightforward computation of gives:
therefore the information content has the expression:
By the corollary 1.4 we obtain the equations:
(11) 
This shows that is a Rayleigh dissipation potential.
Elastoplasticity.
Take a hamiltonian system with space state and momentum space and supplement the state and the momentum with a new pair of state and momentum spaces:
Suppose further that the pair of spaces are in duality, so that we can define
which leads us to a duality product of with itself:
In this setting, we take and
where is a convex lsc function.
There are two differences with respect to the viscosity example: there is a cartesian decomposition of the state space and momentum space, and the dissipation potential depends on (a component of) the momentum variable, while previously the dependence was on the state variable.
By a computation analoguous with the one from the previous example we obtain:
Via the corollary 1.4 we obtain the equations:
(12) 
In particular, let’s take , , Hilbert spaces, and a hamiltonian of the form:
where is the kinetic energy and is the elastic energy. We denote the elastic force by
The 2nd and the 4th equations of the system (12) give:
The system (12) reduces to the familiar equations:
(13) 
In the next example we encounter a small extension of (9).
Damage.
We consider here a standard model of damage. In [2] we introduced a formulation of hamiltonian evolution with dissipation which is the ancestor of the one proposed in this article. Then we applied that formulation to a more involved damage model, where we used an AmbrosioTortorelli functional [1] in the hamiltonian of the model.
As in the elastoplasticity example, we supplement the spaces , with a new pair . Here is a damage variable and is a conjugated variable to the damage one.
Thefore we have a duality product over :
We choose a hamiltonian
where is the kinetic energy, is a term akin to the kinetic energy, but for the variable and is the elastic energy.
The difference from the previous examples is that the damage variable and the sign of are constrained (i.e. damage cannot decrease):
We have to include these constraints in the expression of the information content:
The expression of the information content does not have the form (9), but instead it has the form
with, in this case
We obtain the equations:
(14) 
and the inequations:
(15) 
In the next example we shall need a smoothness assumption. The alternative would be to renounce at the condition 1.2 (a) and to keep only the condition (b).
Unilateral contact.
We take , to be dual, finite dimensional Banach spaces and to be a closed, nonempty set of admissible states . With the notations about Fréchet normals from [7] 1.4, we shall need that is smooth enough, in the sense that the tangent cone at any is the polar of the normal cone at x, denoted by .
There is no special form of the hamiltonian required. The information content has the form:
The smoothness assumption on is that
i.e. that the tangent cone is the polar of the normal cone. In general it is true that the normal cone is the polar of the tangent cone, but not the other way around. An alternative would be to not make any smoothness assumptions on and instead to take an information content of the form
Again, in general the tangent cone is not convex, therefore this form of information content contradicts the convexity condition 1.2 (a).
We obtain the equations:
(16) 
and the unilateral contact conditions:
(17) 
3 Notations and useful definitions
General notations.
The space of states and the space of momenta are real topological vector spaces in duality:
We suppose the usual: the duality is bilinear, continuous and for any linear and continuous functions , there exist , such that and .
The space is in duality with itself by:
The space is also symplectic, with the symplectic form defined by: for any ,
For any differentiable function the gradient of at a point is the element with the property that
and the symplectic gradient of is is defined in a similar way by the equality
If we use the partial derivatives notation
then
With the introduction of the linear conjugation
we have .
Convex analysis notations.
These are the classical ones from Moreau [8]. We add the field of reals . The addition operation is extended with for any . The multiplication with positive numbers is extended with: if then .
For any function , it’s domain is .
The set of lower semicontinuous (lsc), convex functions defined on , with nonempty domain is . The indicator function of a convex and closed set is
For any natural number , a non empty set is monotone monotone if for any collection we have the inequality:
The set is maximally monotone if it is monotone and maximal with respect to the inclusion of sets.
The set is cyclically monotone if it is monotone for any natural number . It is cyclically maximal monotone if it cyclically monotone and maximal with respect to the inclusion of sets.
The subdifferential of a function at a point is the set:
The polar of a function is
The polar is always convex and lsc.
Polars and subgradients are related by the Fenchel inequality. For any function which is convex, lsc, we define
for any , . The Fenchel inequality has two parts:
References
 [1] L. Ambrosio, V. Tortorelli, On the Approximation of Free Discontinuity Problems, Bollettino UMI 7, 6B (1992), 105123
 [2] M. Buliga, Hamiltonian inclusions with convex dissipation with a view towards applications, Mathematics and its Applications 1, 2 (2009), 228251, arXiv:0810.1419.
 [3] M. Buliga, G. de Saxcé, A symplectic BrezisEkelandNayroles principle, Mathematics and Mechanics of Solids 22, 6, (2017)
 [4] M. Buliga, A stochastic version and a Liouville theorem for hamiltonian inclusions with convex dissipation, arXiv:1807.10480
 [5] M. Buliga, G. de Saxcé, C. Vallée, Bipotentials for non monotone multivalued operators: fundamental results and applications, Acta Applicandae Mathematicae, 110, 2(2010), 955972
 [6] M. Buliga, G. de Saxcé, C. Vallée, Blurred maximal cyclically monotone sets and bipotentials , Analysis and Applications 8 (2010), no. 4, 114
 [7] A. Ya. Kruger, On Fréchet subdifferentials, Journal of Mathematical Sciences 16 (2003), no. 3, 33253358
 [8] J.J. Moreau, Fonctionnelles convexes (Rome: Istituto Poligrafico e zecca dello stato, 2003), p. 71, 97.
 [9] R.T. Rockafellar, On the maximal monotonicity of subdifferential mappings, Pacific J. Math. 33, 1 (1970), 209216
 [10] A. Oueslati, A.D. Nguyen, G. de Saxcé (2017). A Symplectic Minimum Variational Principle for Dissipative Dynamical Systems. In: Nielsen F., Barbaresco F. (eds) Geometric Science of Information. GSI 2017. Lecture Notes in Computer Science, vol 10589. Springer, Cham
 [11] G. de Saxcé, Une généralisation de l’inégalité de Fenchel et ses applications aux lois constitutives, C. R. Acad. Sci., Paris, Sér. II 314 (1992) 125129
Comments
There are no comments yet.