On mKdV and associated classes of moving boundary problems: reciprocal connections

A class of Stefan-type moving boundary problems for the canonical modified Korteweg–de Vries (mKdV) equation of soliton theory is solved via application of a similarity reduction to Painlevé II which involves Airy’s equation. A reciprocal transformation is applied to derive a linked class of solvable moving boundary problems for a basic Casimir member of a compacton hierarchy. Application of a class of involutory transformations with origin in an autonomisation procedure for the Ermakov–Ray–Reid system is then used to isolate novel solvable moving boundary problems for Ermakov-modulated mkdV equations.


Introduction
The mKdV equation is, importantly, known to be connected to the canonical Korteweg-de Vries equation of soliton theory via a transformation due to Miura [1,2]. It has physical application in the analysis of nonlinear Alfvén waves in a collisionless plasma [3] and of acoustic wave propagation in an Threshold melting conditions as previously derived by Tarzia [17] and Solomon et al. [18] for analogous types of moving boundary problems for the classical heat equation can be thereby extended.
Reciprocal-type transformations as originally introduced by Bateman [19] were applied in connection with lift and drag phenomena in two-dimensional homentropic gasdynamics. Canonical reduction in subsonic régimes may thereby be obtained via Kármán-Tsien model pressure-density gas laws to the classical Cauchy-Riemann system of hydrodynamics. The gasdynamic reciprocal transformations of [19] which characteristically depend on admitted conservation laws were subsequently shown to constitute particular Bäcklund transformations in [20].
The study of moving boundary problems in continuum mechanics has been and continues to be a subject of significant research interest (see e.g., [21][22][23] and work cited therein). By constrast, the literature on moving boundary problems for solitonic equations is sparse indeed. However, an intriguing solitonic connection of physical importance occurs in the analysis in [24] of the classical Saffman-Taylor problem with surface tension [25]. Therein, a one-parameter class of solutions was isolated in a description of the motion of an interface between a viscous and nonviscous liquid. This class, remarkably, was shown to be associated with travelling wave solutions of the canonical Dym equation of modern soliton theory [26]. Importantly, the latter may be linked to the potential mKdV equation via a reciprocal transformation. Novel integrable extensions of the Dym equation were derived in [27] in the geometric context of the binormal motion of an inextensible curve of constant curvature. Connection was made, in particular, with a Camassa-Holm equation of shallow water hydrodynamics.
Reciprocal transformations were applied in [28] to obtain exact solution to moving boundary problems such as arise in the context of the percolation of liquids through a porous medium such as soil.
In [29], a Painlevé II symmetry reduction was used to derive exact solutions to a class of moving boundary problems for the Dym equation. The latter type of reduction had its origin in a study of the evolution of the interface in a Hele-Shaw problem in [30]. A Bäcklund transformation was applied in [29] to generate iteratively a novel sequence of analytically solvable moving boundary problems for the Dym equation with interface of the type x = t n . Exact solutions were obtained therein in terms of Yablonski-Vorob'ev polynomials [31,32] with n corresponding to a sequence of values of the Painlevé II parameter. The procedure was subsequently applied in [33] to generate such exact solutions involving Yablonski-Vorob'ev polynomials for an extended Dym equation.
Here, a class of Stefan-type moving boundary problems for the mKdV equation is solved by application of a similarity reduction to Painlevé II involving Airy's equation. A novel reciprocal transformation is then introduced to isolate a linked class of solvable moving boundary problems for a canonical base Casimir member of a compacton hierarchy. In addition, a class of involutory transformations with origin in the autonomisation procedure for the Ermakov-Ray-Reid system is used to construct a wide class of solitonic Ermakov-modulated mKdV equations together with associated solvable moving boundary problems of Stefan-type.

A class of moving boundary problems for the mKdV equation
The canonical mKdV equation of modern soliton theory adopts the form and constitutes the base member of the mKdV soliton hierarchy (see e.g., [6] and literature cited therein).
Here, a class of Stefan-type moving boundary problems for the mKdV equation is considered, namely together with In the preceding, L m , P m and H 0 ℝ while i, j, k are indices to be determined by admittance of a viable symmetry reduction. In the latter connection, an ansatz is introduced into the mKdV equation (1)

Airy solution
Painlevé II admits an important class of exact solutions w = − � (z)∕ (z) corresponding to the parameter = 1∕2 where is governed by the classical Airy equation wherein � ∶= d∕dz . Iterative application of a Bäcklund transformation due to Lukashevich [34] to this class generates a chain of exact solutions of P II corresponding to a sequence of Painlevé parameters = k − 1∕2, k ℕ + . Application of this result in the context of stationary boundary value problems for the Nernst-Planck electrolytic system has been detailed in [35].
Here, the class of similarity solutions with is applied to obtain exact solution to a class of moving boundary problems (2) for the solitonic mKdV equation. In the above, Ai and Bi denote Airy functions of the first and second kind respectively.

Boundary conditions
With (13), namely in extenso it is seen that whence Thus, on use of the relation 2 2 = 1 . Accordingly, with

A Casimir connection
Under the reciprocal transformation applied to the mKdV equation (1), it is seen that whence The latter constitutes the base third order member of the compacton hierarchy of [36], namely the Casimir equation

A class of reciprocal moving boundary problems
Here, the reciprocal of the mKdV moving boundary problems determined by (1)-(3) becomes The reciprocal relations yield whence so that on use of the relation 3 = 3 . Accordingly, and the reciprocal moving boundary is determined by

Modulation
Systems which incorporate spatial or temporal modulation arise naturally in both nonlinear physics and continuum mechanics. In physics, they occur notably in nonlinear optics and the theory of Bose-Einstein condensates [33,[37][38][39][40]. In continuum mechanics they have been shown to have importance in the analysis of initial and boundary value problems in elastodynamics, visco-elastodynamics and elastostatics of both inhomogeneous and initially stressed media.
In recent work [41], modulated coupled systems of sine-Gordon, Demoulin and Manakov-type have been systematically reduced to their solitonic unmodulated counterparts via a class of involutory transformations. Here, such transformations of the type are applied to the class of moving boundary problems for the mKdV equation (1) determined by the conditions (2)-(4) to construct an associated class of solvable modulated boundary value problems. It is seen that if the relations in ℝ * are augmented by * = −1 then the involutory property ℝ * * = I results. This kind of transformation was originally applied in [42] in the autonomisation of the Ermakov-Ray-Reid system [43,44]. The latter arises notably in nonlinear optics [45][46][47][48][49][50][51], and in diverse other areas of physical application [52][53][54][55][56][57][58].

Ermakov modulation
If the modulation term (x) in the class of involutory transformations ℝ * is determined by the classical Ermakov equation in = * −1 , namely [59] then the nonlinear superposition principle is admitted in which Ω 1 , Ω 2 constitute a pair of linearly independent solutions of the auxiliary equation and with constants c together with c i , i = 1, 2 such that where W = Ω 1 Ω 2x − Ω 1x Ω 2 is the constant Wronskian of Ω 1 , Ω 2 . This result and generalisations may be obtained via Lie group methods [60]. Physical applications of the classical Ermakov equation (37) arise, in the elastodynamics of boundary-loaded neo-Hookean hyperelastic tubes of Mooney-Rivlin material as detailed in [61].
The class of Ermakov-modulated moving boundary problems with conditions (34)-(35) is connected parametrically to its unmodulated counterpart by the relations In this manner, such nonlinear Ermakov-modulated moving boundary systems may be reduced to their exactly solvable unmodulated counterparts. Modulations determined by specialisations of (39) include, inter alia, those which involve Lamé, Mathieu or Airy terms.

Conclusion
Application of Airy-type solutions of the classical Painlevé II equation has previously been made to time-independent boundary value problems associated with the Nernst-Planck electrolytic system [35]. It was established therein that a chain of exact solutions may be derived by iterated action of a Bäcklund transformation admitted by Painlevé II on the seed Airy class. Here, Airy-type solution of a Painlevé II reduction of a canonical mKdV equation is applied to solve classes of Stefan-type moving boundary problems. The application of iteration of the Painlevé II Bäcklund transformation in this context as conducted in the exact solution of moving boundary problems for the solitonic Dym equation in [29] remains to be investigated.
A moving boundary problem in plasma physics has been previously investigated in [62]. Diverse physical applications of the mKdV equation have been documented in [63]. The type of application of Airy-type representations to moving boundary problems as presented here for the mKdV equation may, in principle, be extended to the solitonic Gardner which represents a composition of the mKdV and KdV equations. This is under current investigation. The Gardner equation has important physical applications, notably in plasma physics, optical lattice theory and nonlinear wave propagation in hydrodynamics. It has recently been shown to arise in nonlinear elastodynamics [64].
Funding Open Access funding enabled and organized by CAUL and its Member Institutions.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The (43) S * (t * ) = x * | x= (t * +a) 1∕3 .
images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.