Next Article in Journal
Cognitively Economical Heuristic for Multiple Sequence Alignment under Uncertainties
Next Article in Special Issue
Ekeland Variational Principle and Some of Its Equivalents on a Weighted Graph, Completeness and the OSC Property
Previous Article in Journal
Up and Down h-Pre-Invex Fuzzy-Number Valued Mappings and Some Certain Fuzzy Integral Inequalities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unfolding a Hidden Lagrangian Structure of a Class of Evolution Equations

School of Mathematical Sciences, Tel Aviv University, Tel Aviv 69978, Israel
Submission received: 3 November 2022 / Revised: 14 December 2022 / Accepted: 16 December 2022 / Published: 20 December 2022
(This article belongs to the Special Issue Principles of Variational Methods in Mathematical Physics)

Abstract

:
It is shown that a simple modification of the standard Lagrangian underlying the dynamics of Newtonian lattices enables one to infer the hidden Lagrangian structure of certain classes of first order in time evolution equations which lack the conventional Lagrangian structure. Implication to other setups is outlined and exemplified.

1. The Problem

Unlike the Roman Deity Janus enshrined with two equal faces, in the calculus of variations most of the attention is given to the direct venue which emerges naturally when a scientific problem is ab initio introduced via the variational principle. However, in many scientific endeavors, in particular when the underlying problem is partially phenomenological, models are derived on the level of differential equations, thus raising a natural quest for an underlying Lagrangian structure which apart from “scientific beauty” would carry many advantages, with conservation laws being one of them. The ‘upstream’ path from Equations to the underlying variational structure is known as an inverse problem of calculus of variations. Unfortunately, our ability to tackle the inverse problem and thus deduce the underlying Lagrangian structure from a given dynamical system is far more modest in spite of the long-standing attempts, c.f., [1,2,3,4,5], which date back to Helmholtz, as those efforts, rather than provide us with a definite algorithm to construct the Lagrangian, provide us the conditions for its existence or statements of equivalence rather than with a specific ‘how’ (notably, there is a far richer literature on Hamiltonian structures than on Lagrangians, though a more balanced approach may be found in texts on dispersive waves. c.f., [6,7,8]).
In the present communication, we unfold the algorithm to determine the underlying Lagrangian of a particular dynamical system as it relates to the K ( n , n ) equations
K ( n , m ) u t = ( u m ) x + ( u n ) x x x , 1 < n , m ,
and a variety of their extensions. Equation (1) were introduced by us some time ago, refs. [9,10] and for 1 < n beget the compaction, solitary wave with a finite span. In the particular case m = n , its compactons take a simple trigonometric form. For instance, when n = 2
u = 4 λ 3 cos 2 x + λ t 4 H ( 2 π | x + λ t | ) ,
with H being the Heaviside function. However, besides compactons, the K ( m , n ) equations had many other fascinating features; at the time the model was introduced, Lagrangian structure did not seem to be one of them. The apparent lack of a Lagrangian structure was the main critique directed at the K ( m , n ) models by Cooper et al. [11], who introduced an alternative, compacton-supporting model derived from an ab initio assumed Lagrangian. In what follows, for the m = n and m = n + 1 cases we shall construct the seemingly ’missing’ Lagrangian and thus, at least in part, refute their critique, but first let us note that to account for other physically relevant cases [12], the K ( m , n ) model was extended to the C ( m ; a , b ) setup [10]
C ( m ; a , b ) u t = ( u m ) x + u a ( u b ) x x x , a n d n a + b ,
which conserves
I 1 = u d x a n d I ω = u ω d x ,
where ω b + 1 a . For a = 0 , Equation (3) reduces to the K ( m , n ) setup and with model in [11] being a particular b = a + 1 case.
As shown in [10], for ω = 2 , i.e., b = a + 1 , the C ( m ; a , a + 1 ) equations follow from a conventional Lagrangian which in terms of u = ψ x reads
L = d t 1 2 ψ t ψ x + ψ x m + 1 m + 1 1 2 ψ x n 1 ψ x x 2 d x ,
adding to the conservation roster the conservation of energy. However, since in the K ( m , n ) equations a = 0 , in this setup only the semi-linear, Korteveg–deVries-like sub-cases are endowed with the Lagrangian structure (5).
Yet, insofar as the m = n and m = n + 1 cases are concerned, i.e., for the K ( n , n ) and K ( n + 1 , n ) cases, the claim in [11] was too sweeping as it refers only to the conventional Lagrangian structure (5). We have recently demonstrated [13], that both setups admit a Lagrangian structure which, however, is ‘hidden from sight’ and, as we shall see shortly, quite different from the conventional one (5). Consequently, the two conservation laws in (4) were appended with an additional, non-local, conservation law.
As we shall see shortly, in both K ( n , n ) and K ( n + 1 , n ) setups, the derivation of the non-standard Lagrangian was based on non-obvious steps/tricks. The present communication aims to re-derive the underlying Lagrangian of the K ( n , n ) setup via a more algorithmic approach, which exploits a certain formal affinity between the K ( n , n ) equations and a dynamical system describing the motion of a Newtonian mass-spring chain which is not only easily deduced from a Lagrangian setup, but is often ab initio formulated via its Lagrangian. Although tying two very different setups may seem at first to be yet another ad hoc trick, as we shall see the idea may be easily extended to other dynamical systems, and thus has potentially a much wider scope. To avoid misunderstanding, we note that we assume the considered problems to be stated on the whole line with solutions vanishing at infinity and thus satisfying the natural boundary conditions implied by the variational derivation.
Let us first briefly summarize the derivation in [13]; see also [10]. To this end, the K ( n , n ) equations are rewritten as a Hamiltonian system
u t = x L 2 δ δ u I ω w h e r e L 2 = 1 + 2 x 2 .
This step is deductive. The maneuver which leads to the breakthrough is to introduce a new variable v via
u = L [ v ] ,
with L understood in a pseudo-differential sense. In terms of L , we have
v t = x δ δ v 1 1 + n L [ v ] 1 + n d x = x L L [ v ] n .
The next, final step consists of introducing a ‘potential’ variable ψ [7,8,9], where v = ψ x that casts Equation (8) into a form which, as per ref. [4], is akin to the Lagrangian
L = d t 1 2 ψ t ψ x 1 n + 1 L [ ψ x ] n + 1 d x
of what, prior to their map, were K ( n , n , ) equations.
Two remarks are in order:
(1) As noted, in the preceding discussion we have tacitly assumed that our system is given on a line with all data vanishing at infinity, thus the natural boundary conditions assumed tacitly in the variational derivation hold trivially. If non-trivial boundary conditions are at play, the issue may become more evolved. However, recalling the very recent work of Olver [14] who further developed the old observation [1,2,3,4,5] that one may modify the variational problem without altering the corresponding Euler–Lagrange equations by adding a null Lagrangian to the integrand, we note that since this modification changes the associated natural boundary conditions, it enables one to enlarge the range of boundary value problems akin to variational techniques. Thus, the core difficulty remains to find the underlying Lagrangian, whereas its extension to handle non-natural boundary conditions may be delegated to the unfolding of a proper null Lagrangian.
(2) Note that the last step of introducing a potential in (5) is well known and apart from its more recent use in dispersive systems [6,7,8], it could be found in the classical textbooks in electrodynamics., c.f., [15], or continuum [16,17].
With a Lagrangian and its invariant properties at hand, in addition to the conservation of mass and energy, we now also have conservation of the momentum v 2 d x , which in terms of the original variables, turns into the conservation of
I 2 = u L 2 [ u ] d x .
It is easily seen that though the presented approach was applied to a very particular problem it may be quite naturally extended. Thus, for instance, one may consider
u t = x L * 2 ( F ( u ) ) ,
where L * 2 is any ‘reasonable’ linear operator and F ( u ) any nice function [4]. For instance, in Ref. [18] addressing synchronisation of oscillators F ( u ) = cos u was assumed. Using the previous notation, in terms of v we have
v t = x L * F ( L * v ) .
In terms of ψ , v = ψ x , the underlying Lagrangian reads,
L = d t { 1 2 ψ t ψ x G L * [ ψ x ] } d x , a n d F = G .
As another extension, consider the second order in time system (c.f., Equation (27))
u t t = x 2 L * 2 ( F ( u ) ) .
In terms of ψ , v = ψ x , its Lagrangian reads,
L = d t { 1 2 ψ t 2 G L * [ ψ x ] } d x ,
and F = G .
In the next section, we take a detour needed to tie our problem with the obvious Lagrangian structure of a Newtonian mass-spring lattice which will render a more natural path toward the deduction of a Lagrangian structure of the K ( n , n ) equations. Yet another take on the lattice–Lagrangian relations is briefly discussed in the Appendix A.

2. The Newtonian Chain

We start with the Hamiltonian of an unrelated mass-particle chain [19,20],
H = N + N m 2 y n ˙ 2 + P y n + 1 y n h .
Let y n = y ( t , n h ) , = h / 2 , m = ρ h , ρ 1 , and D x x , then expansion of the potential P yields
P y n + 1 y n h = P ( y x ) h 2 24 ( y x x ) 2 P ( y x ) + O ( h 4 ) .
To this order of expansion and the standard discrete-continuum association, the ab initio available Lagrangian density becomes
L = 1 2 ( y t ) 2 P ( y x ) + h 2 24 ( y x x ) 2 P ( y x ) d t d x ,
with its equations of motion
y t t = x P ( y x ) + h 2 12 P ( y x ) x ( P ( y x ) y x x ) .
Thus, if P ( s ) = s 4 / 4 and u = y x , then after one differentiation the resulting PDE reads
u t t = u 3 + h 2 8 u ( u 2 ) x x x x .
However, we shall also need an alternative description. Let
M ( ) 2 x e x e x , L D M 1 ,
and v a new variable defined via
y = M ( ) [ v ] o r v = L D [ y ] ,
then in terms of v, the Hamiltonian (13) reads
H = 1 2 ( M v t ) 2 + P ( v x ) d t d x ,
which yields
M 2 v t t = P ( v x ) x .
Acting on (20) with L D 2 and x we obtain
u t t = L D 2 P ( u ) x x w h e r e u = v x .
Actually, without invoking v, we may use y ab initio and rewrite the Hamiltonian (13) in terms of y
H = 1 2 ( y t ) 2 + P ( L D y x ) d t d x ,
with its Lagrangian density
L = 1 2 ( y t ) 2 P ( L D y x ) d t d x ,
which yields
y t t = L D P ( L D y x ) x .
Of course, since u = v x = L D y x , acting with L D x , turns Equation (24) into (21).
Note that if, say, P ( u ) = u 3 and L D = 1 + h 2 12 x 2 , obtained after carrying the expansion of L D to fourth order, Equation (21) turns into a second order in time variant of the K ( n , n ) equation
u t t = u n + h 2 12 ( u n ) x x x x .
With its Lagrangian structure given via (14). Both Equations (18) and (27) and their underlying Lagrangians provide a quasi-continuum rendition of the same lattice, with expansions centering around different nodal locations.

Bridging between the Two Systems

The Newtonian system is second order in time. To turn it into a structurally analogous, but of first order in time system, we modify the kinetic part in (23)
y t 2 y t y x
with the resulting Lagrangian becoming
L = 1 2 y t y x P ( L D y x ) d t d x .
Identifying P with G and ψ with y, begets a Lagrangian form identical to (12), with a related dynamical system which is first order in time. Setting u = L D y x yields Equation (11). In particular, P ( s ) = 1 n + 1 s n yields the K ( n , n ) setup with the resulting sought-after Lagrangian structure.
Finally, consider again the direct expansion of the lattice and its underlying standard Lagrangian (14). Assuming again that P ( s ) = 1 n + 1 s n + 1 , m = n and repeating the y t 2 y t y x association, we recover module normalizable coefficients, the Lagrangian (5) and the underlying Equation (3); C ( m ; a , a + 1 ) , with 2 a = n 1 .
Thus, as with Equations (21) and (25), with each of the expansions being based on a different nodal location on the lattice, we obtain a different equation of motion and a different corresponding Lagrangian: one which is conventional and readily available, and the other hidden and subject of the present exploration.

3. Closing Comments

The present note is a modest take on the inverse problem of the calculus of variations to seek from a particular partial differential system its variational antecedent. At the present state of affairs, the statements regarding the underlying Lagrangian structure of a given dynamical system are affirmative rather than constructive. Every success is ad hoc and relies on certain intuitive steps which have to be inferred independently for each problem.
Exploiting a formal structural affinity between two different physico-mathematical entities wherein one, being second order in time, is endowed with a natural Lagrangian structure, the Lagrangian of the other is constructed via a natural modification. What saves the presented approach from being yet another ad hoc spiel is the fact that the approach carries to a wider set of problems wherein the first order in time dynamical system lacks an obvious underlying Lagrangian, but a better chance may await us looking at structurally similar second-order equations in time endowed with an obvious Lagrangian structure.

Funding

The research was funded in part by the ISF contract 1983/22.

Data Availability Statement

Not applicable.

Acknowledgments

The author thanks P. Olver for his constructive comments.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

Another take on the discussed problem is provided by the following first-order system [8]
d d t u n = N + N f ( u n + 1 ) f ( u n 1 ) 2 .
Let D x = x and L + 2 L D , then
u t = x L + 2 f ( u ) .
If u = L + v , then
v t = x L + f ( L + v ) .
As before, let v = ψ x , then the desired Lagrangian is
L = 1 2 ψ t ψ x Q ( L + ψ x ) d t d x ,
where f ( u ) = Q ( u ) .
We may approximate L + in two ways. In the standard approach we utilize the first two terms L + 2 L a 2 1 + D x 2 6 of the exact expression
L + 2 = 1 + D x 2 3 ! + D x 4 5 ! +
which renders L a a pseudo differential operator. But if instead we adopt
L + 2 L b 2 ( 1 + D x 2 12 ) 2 ,
with L b = 1 + D x 2 12 being a nice differential operator, and within small error covering the first three terms in (A5).

References

  1. Anderson, I.A. Aspects of the inverse problem to the calculus of variations. Arch. Math. 1988, 24, 181–202. [Google Scholar]
  2. Tonti, E. Inverse problem: Its general solution. In Lecture Notes in Pure and Applied Mathematics No. 100: Differential Geometry, Calculus of Variations and Their Applications; Marcel Decker: New York, NY, USA, 1985; pp. 497–510. [Google Scholar]
  3. Giaquinta, M.; Hildebrandt, S. Calculus of Variations I. The Lagrangian Formalism; Springer: New York, NY, USA, 1996. [Google Scholar]
  4. Olver, P.J. Dirac’s theory of constraints in field theory and the canonical form of Hamiltonian differential operators. J. Math. Phys. 1986, 27, 2495. [Google Scholar] [CrossRef]
  5. Olver, P.J. On the Hamiltonian structure of evolution equations. Math. Proc. Camb. Philos. Soc. 1980, 88, 71. [Google Scholar] [CrossRef] [Green Version]
  6. Whitham, G.B. Linear and Nonlinear Waves; John Wiley and Sons: Cambridge, MA, USA, 1974. [Google Scholar]
  7. Infeld, E.; Rowlands, G. Nonlinear Waves, Solitons and Chaos; Cambridge Univesity Press: Cambridge, UK, 1990. [Google Scholar]
  8. Drazin, P.G.; Johnson, R.S. Solitons: An Introduction; Cambridge Univesity Press: Cambridge, UK, 1990. [Google Scholar]
  9. Rosenau, P.; Hyman, J.M. Solitons with finite wavelength. Phys. Rev. Letts. 1993, 70, 564. [Google Scholar] [CrossRef] [PubMed]
  10. Rosenau, P.; Zilburg, A. Compactons. J. Phys. A Math. Theor. 2018, 51, 343001. [Google Scholar] [CrossRef]
  11. Cooper, F.; Shepard, H.; Sodano, P. Solitary waves in a class of generalized Korteweg–de Vries equations. Phys. Rev. E 1993, 48, 4027. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Rubinstein, J. Evolution equations for stratified dilute suspensions. Phys. Fluids A Fluid Dyn. 1990, 2, 3. [Google Scholar] [CrossRef]
  13. Zilburg, A.; Rosenau, P. On Hamiltonian formulations of the C1(m, a, b) equations. Phys. Letts. 2017, A381, 1557. [Google Scholar] [CrossRef]
  14. Olver, P.J. Boundary Conditions and Null Lagrangians in the Calculus of Variations and Elasticity. J. Elast. 2022. [Google Scholar] [CrossRef]
  15. Stratton, J.A. Electromagnetic Theory; McGraw-Hill: New York, NY, USA, 1941. [Google Scholar]
  16. Sieniutycz, S. Conservation Laws in Variational Thermo-Hydrodynamics; Kluwer: Dordrecht, The Netherlands, 1994. [Google Scholar]
  17. Sieniutycz, S.; Berry, R.S. Canonical formalism, fundamental equation, and generalized thermomechanics for irreversible fluids with heat transfer. Phys. Rev. E 1993, 47, 1765. [Google Scholar] [CrossRef] [PubMed]
  18. Rosenau, P.; Pikovsky, A. Solitary phase waves in a chain of autonomous oscillators. Chaos 2020, 30, 053119. [Google Scholar] [CrossRef] [PubMed]
  19. Rosenau, P. Hamiltonian dynamics of dense chains and lattices: Or how to correct the continuum. Phys. Letts. 2003, A311, 39. [Google Scholar] [CrossRef]
  20. Rosenau, P.; Schochet, S. Compact and almost compact breathers: A bridge between an anharmonic lattice and its continuum limit. Chaos 2005, 15, 015111. [Google Scholar] [CrossRef] [PubMed]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rosenau, P. Unfolding a Hidden Lagrangian Structure of a Class of Evolution Equations. Axioms 2023, 12, 2. https://doi.org/10.3390/axioms12010002

AMA Style

Rosenau P. Unfolding a Hidden Lagrangian Structure of a Class of Evolution Equations. Axioms. 2023; 12(1):2. https://doi.org/10.3390/axioms12010002

Chicago/Turabian Style

Rosenau, Philip. 2023. "Unfolding a Hidden Lagrangian Structure of a Class of Evolution Equations" Axioms 12, no. 1: 2. https://doi.org/10.3390/axioms12010002

APA Style

Rosenau, P. (2023). Unfolding a Hidden Lagrangian Structure of a Class of Evolution Equations. Axioms, 12(1), 2. https://doi.org/10.3390/axioms12010002

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop