Mach-type soliton in the Novikov-Veselov equation
Using the reality condition of the solutions, one constructs the Mach-type soliton of the Novikov-Veselov equation by the minor-summation formula of the Pfaffian. We study the evolution of the Mach-type soliton and find that the amplitude of the Mach stem wave is less than two times of the one of th...
Gespeichert in:
| Veröffentlicht in: | Symmetry, Integrability and Geometry: Methods and Applications |
|---|---|
| Datum: | 2014 |
| 1. Verfasser: | |
| Format: | Artikel |
| Sprache: | English |
| Veröffentlicht: |
Інститут математики НАН України
2014
|
| Online Zugang: | https://nasplib.isofts.kiev.ua/handle/123456789/146342 |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Zitieren: | Mach-type soliton in the Novikov-Veselov equation/ Jen-Hsu Chang // Symmetry, Integrability and Geometry: Methods and Applications. — 2014. — Т. 10. — Бібліогр.: 31 назв. — англ. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| id |
nasplib_isofts_kiev_ua-123456789-146342 |
|---|---|
| record_format |
dspace |
| spelling |
Jen-Hsu Chang 2019-02-09T08:55:25Z 2019-02-09T08:55:25Z 2014 Mach-type soliton in the Novikov-Veselov equation/ Jen-Hsu Chang // Symmetry, Integrability and Geometry: Methods and Applications. — 2014. — Т. 10. — Бібліогр.: 31 назв. — англ. 1815-0659 2010 Mathematics Subject Classification: 35C08; 35A22 DOI:10.3842/SIGMA.2014.111 https://nasplib.isofts.kiev.ua/handle/123456789/146342 Using the reality condition of the solutions, one constructs the Mach-type soliton of the Novikov-Veselov equation by the minor-summation formula of the Pfaffian. We study the evolution of the Mach-type soliton and find that the amplitude of the Mach stem wave is less than two times of the one of the incident wave. It is shown that the length of the Mach stem wave is linear with time. One discusses the relations with V-shape initial value wave for different critical values of Miles parameter. The author thanks the referees for their valuable suggestions. This work is supported in part by the National Science Council of Taiwan under Grant No. NSC 102-2115-M-606-001. en Інститут математики НАН України Symmetry, Integrability and Geometry: Methods and Applications Mach-type soliton in the Novikov-Veselov equation Article published earlier |
| institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| collection |
DSpace DC |
| title |
Mach-type soliton in the Novikov-Veselov equation |
| spellingShingle |
Mach-type soliton in the Novikov-Veselov equation Jen-Hsu Chang |
| title_short |
Mach-type soliton in the Novikov-Veselov equation |
| title_full |
Mach-type soliton in the Novikov-Veselov equation |
| title_fullStr |
Mach-type soliton in the Novikov-Veselov equation |
| title_full_unstemmed |
Mach-type soliton in the Novikov-Veselov equation |
| title_sort |
mach-type soliton in the novikov-veselov equation |
| author |
Jen-Hsu Chang |
| author_facet |
Jen-Hsu Chang |
| publishDate |
2014 |
| language |
English |
| container_title |
Symmetry, Integrability and Geometry: Methods and Applications |
| publisher |
Інститут математики НАН України |
| format |
Article |
| description |
Using the reality condition of the solutions, one constructs the Mach-type soliton of the Novikov-Veselov equation by the minor-summation formula of the Pfaffian. We study the evolution of the Mach-type soliton and find that the amplitude of the Mach stem wave is less than two times of the one of the incident wave. It is shown that the length of the Mach stem wave is linear with time. One discusses the relations with V-shape initial value wave for different critical values of Miles parameter.
|
| issn |
1815-0659 |
| url |
https://nasplib.isofts.kiev.ua/handle/123456789/146342 |
| citation_txt |
Mach-type soliton in the Novikov-Veselov equation/ Jen-Hsu Chang // Symmetry, Integrability and Geometry: Methods and Applications. — 2014. — Т. 10. — Бібліогр.: 31 назв. — англ. |
| work_keys_str_mv |
AT jenhsuchang machtypesolitoninthenovikovveselovequation |
| first_indexed |
2025-11-27T00:42:59Z |
| last_indexed |
2025-11-27T00:42:59Z |
| _version_ |
1850789185714126848 |
| fulltext |
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 10 (2014), 111, 14 pages
Mach-Type Soliton in the Novikov–Veselov Equation
Jen-Hsu CHANG
Department of Computer Science and Information Engineering,
National Defense University, Tauyuan County 33551, Taiwan
E-mail: jhchang@ndu.edu.tw
Received September 18, 2014, in final form December 10, 2014; Published online December 18, 2014
http://dx.doi.org/10.3842/SIGMA.2014.111
Abstract. Using the reality condition of the solutions, one constructs the Mach-type soliton
of the Novikov–Veselov equation by the minor-summation formula of the Pfaffian. We study
the evolution of the Mach-type soliton and find that the amplitude of the Mach stem wave
is less than two times of the one of the incident wave. It is shown that the length of the
Mach stem wave is linear with time. One discusses the relations with V -shape initial value
wave for different critical values of Miles parameter.
Key words: Pfaffian; Mach-type soliton; Mach stem wave; V -shape wave
2010 Mathematics Subject Classification: 35C08; 35A22
1 Introduction
Recently, the resonance theory of line solitons of KP-(II) equation (shallow water wave equation)
∂x(−4ut + uxxx + 6uux) + 3uyy = 0
has attracted much attractions using the totally non-negative Grassmannian [1, 3, 5, 6, 16,
22, 23], that is, those points of the real Grassmannian whose Plucker coordinates are all non-
negative. For the KP-(II) equation case, the τ -function is described by the Wronskian form
with respect to x. The Mach reflection problem describes the resonant interaction of solitary
waves appearing in the reflection of an obliquely incident waves onto a vertical wall. John
Miles discussed an oblique interaction of solitary waves and found a resonant interaction to
describe the Mach reflection phenomena [26]. In this end, he predicts an extraordinary fourfold
application of the wave at the wall. The Miles theory in terms of the KP equation and the
Mach-type solution in KP observed experimentally are investigated in [16, 17, 18, 20, 31] (and
references therein). The point is that irregular reflection can be described by the (3142)-type
soliton and the stem in the middle part should be a Mach stem wave. Inspired by their works,
one can consider the Novikov–Veselov equation similarly.
One considers the Novikov–Veselov (NV) equation [4, 11, 15, 28, 30] with real solution U :
Ut = Re
[
∂3
zU + 3∂z(QU)− 3ε∂zQ
]
, (1.1)
∂̄zQ = ∂zU, t ∈ R,
where ε is a real constant. The NV equation (1.1) is one of the natural generation of the famous
KdV equation and can have the Manakov’s triad representation [24]
Lt = [A,L] +BL,
where L is the two-dimension Schrödinger operator
L = ∂z∂̄z + U − ε
jhchang@ndu.edu.tw
http://dx.doi.org/10.3842/SIGMA.2014.111
2 J.H. Chang
and
A = ∂3
z +Q∂z + ∂̄3
z + Q̄∂̄z, B = Qz + Q̄z̄.
It is equivalent to the linear representation
Lφ = 0, ∂tφ = Aφ. (1.2)
We remark that when ε → ±∞, the Veselov–Novikov equation reduces to the KP-I (ε → −∞)
and KP-(II) (ε→∞) equation respectively [10]. To make a comparison with KP-(II) equation,
we only consider ε > 0.
Let φ1, φ2 be any two independent solutions of (1.2). Then one can construct the extended
Moutard transformation using the skew product [2, 13, 25]
W (φ1, φ2) =
∫
(φ1∂φ2 − φ2∂φ1)dz −
(
φ1∂̄φ2 − φ2∂̄φ1
)
dz̄ +
[
φ1∂
3φ2 − φ2∂
3φ1
+ φ2∂̄
3φ1 − φ1∂̄
3φ2 + 2
(
∂2φ1∂φ2 − ∂φ1∂
2φ2
)
− 2
(
∂̄2φ1∂̄φ2 − ∂̄φ1∂̄
2φ2
)
+ 3Q(φ1∂φ2 − φ2∂φ1)− 3Q̄
(
φ1∂̄φ2 − φ2∂̄φ1
)]
dt, (1.3)
such that
Û(t, z, z̄) = U(t, z, z̄) + 2∂∂̄ lnW (φ1, φ2), Q̂(t, z, z̄) = Q(t, z, z̄) + 2∂∂ lnW (φ1, φ2),
is also a solution of the NV equation (1.1).
For fixed potential U0(z, z̄, t) and Q0(z, z̄, t) of the NV equation (1.1), we can take any 2N
wave functions φ1, φ2, φ3, . . . , φ2N (or their linear combinations) of (1.2). Then the 2N -step
successive extended Moutard transformation can be expressed as the Pfaffian form [2, 27] (also
see [12, 29])
U = U0 + 2∂∂̄[ln Pf(φ1, φ2, φ3, . . . , φ2N )], Q = Q0 + 2∂∂[ln Pf(φ1, φ2, φ3, . . . , φ2N )],
where Pf(φ1, φ2, φ3, . . . , φ2N ) is the Pfaffian defined by
Pf(φ1, φ2, φ3, . . . , φ2N ) =
∑
σ
ε(σ)Wσ1σ2Wσ3σ4 · · ·Wσ2N−1σ2N ,
and Wσiσj = W (φσ(i), φσ(j)) is the extended Moutard transformation (1.3), σ being some per-
mutations.
To construct the N -solitons solutions, we take V = U = 0 in (1.1) and then (1.2) becomes
∂∂̄φ = εφ, φt = φzzz + φz̄z̄z̄, (1.4)
where ε is non-zero real constant. The general solution of (1.4) can be expressed as
φ(z, z̄, t) =
∫
Γ
e
(iλ)z+(iλ)3t+ ε
iλ
z̄+ ε3
(iλ)3
t
ν(λ)dλ, (1.5)
where ν(λ) is an arbitrary distribution and Γ is an arbitrary path of integration such that the
r.h.s. of (1.5) is well defined. One takes νm(λ) = δ(λ − pm), where pm is a complex number.
Define
φm =
φ(pm)√
3
=
1√
3
eF (pm),
Mach-Type Soliton in the Novikov–Veselov Equation 3
where
F (λ) = (iλ)z + (iλ)3t+
ε
iλ
z̄ +
ε3
(iλ)3
t.
Plugging (φm, φn) into the extended Moutard transformation (1.3), we obtain
W (φm, φn) = i
pn − pm
pn + pm
eF (pm)+F (pn). (1.6)
To study resonance, we introduce the real Grassmannian (or the 2N×M matrix) to construct
N solitons. To this end, one considers linear combination of φn. Let
~Ψ = (φ1, φ2, φ3, . . . , φM )T
and H be an 2N ×M (2N ≤M) of real constant matrix (or Grassmannian). Suppose that
H~Ψ = ~Ψ∗ = (Ψ∗1,Ψ
∗
2,Ψ
∗
3, . . . ,Ψ
∗
2N )T ,
that is,
Ψ∗n = hn1φ1 + hn2φ2 + hn3φ3 + · · ·+ hnMφM , 1 ≤ n ≤ 2N.
Then one has by the minor-summation formula [14, 19]
τN = Pf(Ψ∗1,Ψ
∗
2,Ψ
∗
3, . . . ,Ψ
∗
2N ) = Pf
(
HWMH
T
)
=
∑
I⊂[M ], ]I=2N
Pf
(
HI
I
)
det(HI), (1.7)
where the M ×M matrix WM is defined by the element (1.6) and HI
I denote the 2N ×M
submatrix of H obtained by picking up the rows and columns indexed by the same index
set I. By this formula, the resonance of real solitons of the Novikov–Veselov equation can
be investigated just like the resonance theory of KP-(II) equation [16, 21, 23]. Finally, the
N -solitons solutions are defined by [7, 9]
U(z, z̄, t) = 2∂∂̄ ln τN (z, z̄, t), V (z, z̄, t) = 2∂∂ ln τN (z, z̄, t).
To obtain the real potential U , the following reality conditions [9] for resonance have to be
considered
|pk|2 = |qk|2 = ε > 0, k = 1, 2, 3, . . . ,m,
given m pairs of complex numbers (p1, q1), (p2, q2), . . . , (pm, qm).
Letting pm =
√
εeiαm and removing i factor from (1.6) afterwards, one has
W (φm, φn) = − tan
αn − αm
2
eφmn , (1.8)
where
φmn = F (pm) + F (pn) = −2
√
ε[x(sinαm + sinαn) + y(cosαm + cosαn)]
+ 2tε
√
ε(sin 3αm + sin 3αn). (1.9)
Therefore, given a 2N ×M matrix H, the associated τH -function can be written as by (1.7) [8]
τH =
∑
I⊂[M ], ]I=2N
ΓIΛI(x, y, t),
4 J.H. Chang
where
ΛI(x, y, t) = Pf(W2N ) = (−1)N
2N∏
i=2, i>j
tan
αi − αj
2
e
2N∑
m=1
F (pm)
,
ΓI being the 2N × 2N minor for the columns with the index set I = {i1, i2, i3, . . . , i2N}. Also,
to keep τH totally positive (or totally negative), we assume that the matrix H belongs to the
totally non-negative Grassmannian [22, 23] and the angle αn satisfies the following condition:
−π
2
≤ α1 < α2 < α3 < · · · < αM−1 < αM ≤
π
2
.
For one-soliton solution, we have, −π
2 ≤ αi < αj < αk ≤ π
2
τ1 = tan
αi − αj
2
eφij + a tan
αi − αk
2
eφik
= aeφik tan
αi − αk
2
[
1 +
1
a
tan
αi−αj
2
tan αi−αk
2
eF (pj)−F (pk)
]
= aeφik tan
αi − αk
2
[
1 + eF (pj)−F (pk)+θjk
]
,
where a is a constant and the phase shift
θjk = ln
1
a
tan
αi−αj
2
tan αi−αk
2
= ln
tan
αi−αj
2
tan αi−αk
2
− ln a.
Hence the real one-soliton solution is [8]
U = 2∂z∂z̄ ln aeφik tan
αi − αk
2
[
1 + eF (pj)−F (pk)+θjk
]
= 2∂z∂z̄
[
1 + eF (pj)−F (pk)+θjk
]
=
1
2
∣∣pk − pj∣∣2 sech2
[
F (pj)− F (pk) + θjk
2
]
= 2ε sin2
(
αk − αj
2
)
sech2
[
F (pj)− F (pk) + θjk
2
]
= A[j,k] sech2 1
2
(
~K[j,k] · ~x−Ω[j,k]t+ θjk
)
. (1.10)
From (1.9) the amplitude A[j,k], the wave vector ~K[j,k] and the frequency Ω[j,k] are defined by
A[j,k] = 2ε sin2
(
αk − αj
2
)
,
~K[j,k] = 2
√
ε(− sinαj + sinαk,− cosαj + cosαk),
Ω[j,k] = 2ε
√
ε[− sin 3αj + sin 3αk], (1.11)
The direction of the wave vector ~K[j,k] =
(
Kx
[j,k],K
y
[j,k]
)
is measured in the clockwise sense from
the y-axis and it is given by
Ky
[j,k]
Kx
[j,k]
=
− cosαj + cosαk
− sinαj + sinαk
= − tan
αj + αk
2
,
that is,
αj+αk
2 gives the angle between the line soliton and the y-axis in the clockwise sense. In
addition, the soliton velocity V[j,k] is [8]
V[j,k] =
ε
4
sin 3αk − sin 3αj
sin2 αj−αk
2
(sinαk − sinαj , cosαk − cosαj) . (1.12)
Mach-Type Soliton in the Novikov–Veselov Equation 5
The paper is organized as follows. In Section 2, one investigates Mach-type or (3142)-type
soliton for the Novikov–Veselov equation. One shows the evolution of the Mach-type soliton and
obtains the relation of the amplitude of the Mach stem wave ([1,4]-soliton) with the one of the
incident wave ([1,3]-soliton). Furthermore, the length of the Mach stem wave is linear with time.
In Section 3, we discuss the relations with V -shape initial value wave for different critical value
of Miles parameter κ. It is shown that the amplitude of the Mach stem wave is less than two
times of the one of the incident wave. In Section 4, we conclude the paper with several remarks.
2 Mach type soliton
In this section, we investigate the Mach-type or (3142)-type soliton. The corresponding totally
non-negative Grassmannian is the the matrix [16]
HM =
[
1 a 0 −c
0 0 1 b
]
,
where a, b, c are positive numbers. When c = 0, one has the O-type soliton for the Novikov–
Veselov equation. For V -shape initial value wave, one can introduce parameter κ to determine
the evolution into Mach-type or O-type soliton (see next section). We remark that the Y -
shape, O-type, and P -type solitons for the Novikov–Veselov equation are investigated in [8].
Now,
HM
1√
3
[φ(p1), φ(p2), φ(p3), φ(p4)]T =
1√
3
[
φ(p1) + aφ(p2)
φ(p3) + bφ(p4)
]
=
[
Ψ∗1
Ψ∗2
]
.
A direct calculation yields by (1.7), (1.8) and (1.9)
τM = W (Ψ∗1,Ψ
∗
2) = W (φ1, φ3) + bW (φ1, φ4) + aW (φ2, φ3) + abW (φ2, φ4) + cW (φ3, φ4)
= tan
α1 − α3
2
e−2
√
ε[x(sinα1+sinα3)+y(cosα1+cosα3)]+2tε
√
ε(sin 3α1+sin 3α3)
+ b tan
α1 − α4
2
e−2
√
ε[x(sinα1+sinα4)+y(cosα1+cosα4)]+2tε
√
ε(sin 3α1+sin 3α4)
+ a tan
α2 − α3
2
e−2
√
ε[x(sinα2+sinα3)+y(cosα2+cosα3)]+2tε
√
ε(sin 3α2+sin 3α3)
+ ab tan
α2 − α4
2
e−2
√
ε[x(sinα2+sinα4)+y(cosα2+cosα4)]+2tε
√
ε(sin 3α2+sin 3α4)
+ c tan
α3 − α4
2
e−2
√
ε[x(sinα3+sinα4)+y(cosα3+cosα4)]+2tε
√
ε(sin 3α3+sin 3α4), (2.1)
where
−π
2
≤ α1 < α2 < α3 < α4 ≤
π
2
.
To investigate the asymptotic behavior for |y| → ∞, we use the notation [16], considering the
line x = −cy,
ηm(c) = −c sinαm + cosαm.
When ηm(c) = ηn(c), one gets
c =
cosαm − cosαn
sinαm − sinαn
= − tan
αm + αn
2
.
6 J.H. Chang
Since
[ηm(c)− ηi(c)]
∣∣∣
c=− tan
αi+αj
2
= cosαm − cosαi + tan
αi + αj
2
(sinαm − sinαi)
= (sinαm − sinαi)
[
tan
αi + αj
2
− tan
αi + αm
2
]
,
we have the following order relations among the other ηm(c)′s along c = − tan
αi+αj
2{
ηi = ηj < ηm if i < m < j,
ηi = ηj > ηm if m < i or m > j.
Then by a similar argument in [16], one knows that by (1.10):
(a) For y � 0, there are two unbounded line solitons, whose types from left to right are
[1, 3], [3, 4].
(b) For y � 0, there are two unbounded line solitons, whose types from left to right are
[4, 2], [2, 1].
It can be verified by the Maple software.
Now, we can discuss the relations between the parameters a, b, c and phase shifts of these line
solitons. Let us first consider the line solitons in x� 0. There are two solitons which are [3,4]-
soliton and [2,1]-soliton. The [3,4]-soliton is obtained by the balance between the exponential
terms W (φ1, φ3) and bW (φ1, φ4), and the [2,1]-soliton is obtained by the balance between the
exponential terms W (φ1, φ3) and aW (φ2, φ3). Therefore, the phase shifts of [3,4]-soliton and
[2,1]-soliton for x� 0 are given by
θ+
[3,4] = ln
tan α3−α1
2
tan α4−α1
2
− ln b, θ+
[2,1] = ln
tan α3−α1
2
tan α3−α2
2
− ln a.
For the line solitons in x � 0, there are two solitons, which are [1,3]-soliton and [4,2]-soliton.
The [1,3]-soliton is obtained by the balance between the exponential terms cW (φ3, φ4) and
bW (φ1, φ4), and the [4,2]-soliton is obtained by the balance between the exponential terms
cW (φ3, φ4) and aW (φ2, φ3). Therefore, the phase shifts of [1,3]-soliton and [4,2]-soliton for
x� 0 are given by
θ−[1,3] = ln
tan α4−α1
2
tan α4−α3
2
+ ln
b
c
, θ−[4,2] = ln
tan α3−α2
2
tan α4−α3
2
+ ln
a
c
.
So one can see that
θ−[1,3] + θ+
[3,4] = θ−[4,2] + θ+
[2,1] = total phase shift = ln
tan α3−α1
2
tan α4−α3
2
− ln c.
We define the parameter s (representing the total phase shift) by s = e
−θ−
[4,2]
−θ+
[2,1] , which leads
to
a =
tan α3−α1
2
tan α3−α2
2
se
θ−
[4,2] , b =
tan α3−α1
2
tan α4−α1
2
se
θ−
[1,3] , c =
tan α3−α1
2
tan α4−α3
2
s.
Hence we know that the three parameters a, b, c can be used to determine the locations of three
asymptotic line solitons, that is, two in x� 0 and one in x > 0. The s-parameter represents the
Mach-Type Soliton in the Novikov–Veselov Equation 7
Figure 1. The middle portion, having maximum amplitude, is the [1,4]-soliton (stem wave). The y-axis
is slightly enlarged to make the middle portion longer.
Figure 2. α1 = − 23
50π, α2 = − 1
5π, α3 = 1
5π, α4 = 23
50π, a = b = c = 1, ε = 5.
relative locations of the intersection point of the [1,3]-soliton and [3,4]-soliton with the x-axis.
Especially, when s = 1, θ−[4,2] = 0, θ−[1,3] = 0, all of the four solitons will intersect at (0, 0) when
t = 0. One remarks that the bounded line soliton [1,4] (Mach stem wave), obtained by the
balance between the exponential terms W (φ1, φ3) and cW (φ3, φ4), has the maximal amplitude
among all the solitons by (1.11) (Fig. 1) and the velocity is obtained by (1.12). Furthermore,
when t < 0, there is a bounded line [2,3]-soliton (Fig. 2, the left side of the triangle), obtained by
the balance between the exponential terms abW (φ2, φ4) and cW (φ3, φ4).
8 J.H. Chang
Figure 3. Initial wave.
Now, we consider the case α3 = −α2 ≥ 0, α4 = −α1 ≥ 0, and the amplitude
A = A[1,3] = A[4,2] ≤ 2ε (2.2)
is fixed. Then one can see that [1, 3]-soliton and [4, 2]-soliton is symmetric to the x-axis and
similarly for [3, 4]-soliton and [2, 1]-soliton. By (1.11), one knows
α3 + α1
2
≤ α3 − α1
2
=
α3 + α4
2
= arcsin
√
A
2ε
.
Therefore the angle between the [1,3]-soliton and the y-axis (counter-clockwise) is less than the
the angle between the [3,4]-soliton and the y-axis (clockwise). We see that given A and 2ε
there is a critical angle ϕC = arcsin
√
A
2ε for the angle between the [1,3]-soliton and the y-axis
(counter-clockwise). Then one can introduce the following Miles-parameter [5, 8, 26] to describe
the interaction for the Mach-type solution, noticing α3+α1
2 ≤ 0,
κ =
| tan α3+α1
2 |
tan α3+α4
2
=
| tan α3+α1
2 |
tanϕC
=
| tan α3+α1
2 |√
A
2ε−A
≤ 1. (2.3)
From (1.11), we have thus using κ
A = A[1,3] = A[4,2] =
2ε(tanϕC)2
1 + (tanϕC)2
, A[3,4] = A[2,1] =
2ε(tanϕC)2
1
κ2
+ (tanϕC)2
≤ A,
A[1,4] = 2ε sin2 α4 − α1
2
= 2ε
[
sin
(
α4 − α3
2
+
α3 − α1
2
)]2
=
2ε(tanϕC)2(κ+ 1)2
[1 + (tanϕC)2][1 + κ2(tanϕC)2]
= A
(κ+ 1)2
[1 + κ2(tanϕC)2]
< 4A. (2.4)
Remark. To make a comparison with KP-(II), we see that
2ε−A = 2ε
(
1− sin2 α3 − α1
2
)
= 2ε cos2 α3 − α1
2
.
Mach-Type Soliton in the Novikov–Veselov Equation 9
Figure 4. The stem wave moves to the left.
When ε→∞, α1 → −π
2 and α3 → π
2 such that
ε cos2 α3 − α1
2
=
1
4
.
Then
κ→
| tan α3+α1
2 |
√
2A
,
which is the Miles parameter of KP-(II) to describe the interactions of water wave solitons [16,
17, 18].
Since the [1,4]-soliton (Mach stem wave) is increasing its length with time but its end points
will lie in a line (see Figs. 4 and 5), we can obtain them as follows.
We choose s = 1, θ−[4,2] = 0, θ−[1,3] = 0 such that [1,3]-soliton and [1,4]-soliton will intersect at
(0, 0) when t = 0 (see Fig. 3). From (1.10)), the ridges of [1,3]-soliton and [1,4]-soliton are given
by F (p1)− F (p3) = 0, F (p1)− F (p4) = 0, which lead to
x(− sinα1 + sinα3) + y(− cosα1 + cosα3) + tε(sin 3α1 − sin 3α3) = 0,
x(− sinα1 + sinα4) + y(− cosα1 + cosα4) + tε(sin 3α1 − sin 3α4) = 0.
Noticing that α3 = −α2 ≥ 0, α4 = −α1 ≥ 0, one gets
x =
tε sin 3α4
sinα4
= tε(4 cos2 α4 − 1), (2.5)
y =
x(sinα1 − sinα3) + tε(− sin 3α1 + sin 3α3)
− cosα1 + cosα3
= tε
(4 cos2 α4 − 1)(sinα1 − sinα3) + (− sin 3α1 + sin 3α3)
− cosα1 + cosα3
(2.6)
= 4tε
sinα3(sinα3 + sinα4)(sinα4 − sinα3)
cosα3 − cosα4
= 4tε sinα3(sinα3 + sinα4) cot
α3 + α4
2
.
10 J.H. Chang
Figure 5. The length of the stem wave is increasing.
Using (2.3), one has
sinα3 = sin
(
α3 + α4
2
− α4 − α3
2
)
=
(1− κ) tanϕC√
[1 + (tanϕC)2][1 + κ2(tanϕC)2]
,
cosα4 = cos
(
α3 + α4
2
+
α4 − α3
2
)
=
1− κ(tanϕC)2√
[1 + (tanϕC)2][1 + κ2(tanϕC)2]
,
sinα3 + sinα4 = 2 sin
α3 + α4
2
cos
α3 − α4
2
=
2 tanϕC√
[1 + (tanϕC)2][1 + κ2(tanϕC)2]
,
4 cos2 α4 − 1 =
4[1− κ(tanϕC)2]2
[1 + (tanϕC)2][1 + κ2(tanϕC)2]
− 1
=
3 + (tanϕC)2[3κ2(tanϕC)2 − κ2 − 8κ− 1]
[1 + (tanϕC)2][1 + κ2(tanϕC)2]
. (2.7)
A simple calculation yields using (2.4)
y = 8tε tanϕC
1− κ
[1 + (tanϕC)2][1 + κ2(tanϕC)2]
= 4tA[1,4]
1− κ
(1 + κ2)(tanϕC)
,
tanχ =
y
x
=
8(1− κ) tanϕC
3 + (tanϕC)2[3κ2(tanϕC)2 − κ2 − 8κ− 1]
.
Hence one knows that the length of [1,4]-soliton is linear with time and its end points will lie
in a line having slope ± tanχ (see Figs. 4 and 5). Furthermore, from (2.5), one gets that the
[1,4]-soliton moves to the right if α4 <
π
3 , and moves to the left if α4 >
π
3 . In particular, if
α4 = π
3 or by (2.7)
3 + (tanϕC)2
[
3κ2(tanϕC)2 − κ2 − 8κ− 1
]
= 0. (2.8)
then [1,4]-soliton’s length is increasing along the y-axis. When κ = 1 (or α3 = 0), one has A = ε
2
by (2.8) and α4 = π
3 . In this special case, the soliton is fixed. It is different from the KP-(II)
case [18, 20].
Mach-Type Soliton in the Novikov–Veselov Equation 11
3 Relations with V -shape initial value waves
In this section, we investigate some relations with the V -shape initial value wave for the Novikov–
Veselov equation (1.1), ε being fixed, as compared with the KP-(II) case [16, 17, 18, 20]. The
main purpose is to study the interactions between line solitons, especially for the meaning of
the critical angle ϕC .
Recalling the one-soliton solution (1.10) and (1.11), one considers the initial data given in
the shape of V with amplitude A and the oblique angle ϕI < 0 (measured in the clockwise sense
from the y-axis):
A sech2
[√
2A cosϕI(x− |y| tanϕI)
]
. (3.1)
For simplicity, one considers A ≤ 2ε. We notice here the V -shape initial wave is in the negative x-
region. The main idea is that we can think the initial value wave as a part of Mach-type
soliton (2.1) or O-type soliton [8], that is, c = 0 in (2.1). In order to identify those soliton
solutions from the V -shape (3.1), we denote them as [i+, j+]-soliton for y � 0 and [i−, j−]-
soliton for y � 0. Solitons for y → ±∞ have by (1.10)
A = 2ε sin2 αj+ − αi+
2
= 2ε sin2 αi− − αj−
2
, (3.2)
ϕI =
αj+ + αi+
2
= −
(
αi− + αj−
2
)
.
Assume that i+ < j+ and i− > j−. Then symmetry gives
αi+ = −αi− , αj+ = −αj− . (3.3)
Using the parameter (2.3) [5, 8, 26]
κ =
| tanϕI |√
A
2ε−A
=
| tanϕI |
tanϕC
,
one can yield, noticing that ϕC =
αj+−αi+
2 =
αi−−αj−
2 = arctan
√
A
2ε−A from (3.2),
• κ ≥ 1⇒ |ϕI | ≥ ϕC ⇒ −π
2 ≤ αi+ < αj+ < αj− < αi− ≤ π
2 (O-type),
• 0 < κ < 1⇒ |ϕI | < ϕC ⇒ −π
2 ≤ αi+ < αj− < αj+ < αi− ≤ π
2 (Mach-type).
We remark here that if κ = 1 (or α3 = 0), then it is of O-type by (2.4) and (2.6). One can
see that if the angle ϕI is small, then an intermediate wave called the Mach stem ([1,4]-soliton)
appears. The Mach stem, the incident wave ([1,3]-soliton) and the reflected wave ([3,4]-soliton)
interact resonantly, and those three waves form a resonant triple. It is similar to the KP-(II)
case [16].
Let’s compute the maximal amplitude of the Mach stem ([1,4]-soliton) for fixed amplitude A
and ε. By (2.4), a simple calculation shows that
dA[1,4]
dκ
= A
2(κ+ 1)[1− κ(tanϕC)2]
[1 + κ2(tanϕC)2]2
. (3.4)
Hence one can see that when κ = 1/(tanϕC)2, that is,
(tanϕC)(| tanϕI |) = 1,
12 J.H. Chang
the Mach stem has the maximal amplitude. Consequently, if
ϕC + ϕI =
π
2
, (3.5)
then one obtains by (3.4), recalling that 0 < κ < 1 (or tanϕC > 1, i.e., A > ε),
Amax
[1,4] = A
(
1 +
1
(tanϕC)2
)
= 2ε < 2A. (3.6)
Therefore one sees that from (3.4), A and ε being fixed,
• 0 < κ < 1
(tanϕC)2
, the amplitude A[1,4] (stem wave) is increasing;
• κ = 1
(tanϕC)2
(or (3.5)), the amplitude A[1,4] has the maximal value 2ε;
• 1
(tanϕC)2
< κ < 1, the amplitude A[1,4] is decreasing.
It is noteworthy that the maximal amplitude is independent of A. Also, we know that the
maximal amplitude of Mach stem for NV equation is less than twice of the incident wave’s
one; however, for the KP equation (shallow water waves), the Mach stem’s amplitude can be
four times of the incident wave’s one [18]. This is the different point from the case of the KP
equation.
On the other hand, one can see that for κ > 1 (O-type) we have 0 ≤ αj+−αi+
2 ≤ π
4 by (3.3),
that is, A ≤ ε. Thus, if we choose A such that
ε < A ≤ 2ε, (3.7)
we get π
2 < αi− ≤ π; therefore, under the condition (3.7), the initial value wave (3.1) would
develop into a singular O-type soliton by (2.1) (c = 0) when ε is fixed. On the other hand, when
|ϕI | ≤ π
2 , A and κ are fixed, one can choose
ε =
A
2
[
1 +
(
κ
tanϕI
)2
]
≥ A
2
.
Then we can obtain regular soliton solutions.
Finally, from (2.5) one remarks that the [1,4]-soliton (stem wave) moves to the right if α4 <
π
3 ,
and moves to the left if α4 >
π
3 . The former case is different from the KP equation (shallow water
waves); i.e., the stem wave moves with the same side of incident wave for the KP equation. On
the other hand, if we replace the condition (2.2) by A = A[3,4] = A[2,1] ≤ 2ε, then by (1.11) the
[3,4]-soliton (the incident wave) has smaller amplitude than the [1,3]-soliton’s one (the reflected
wave). But this is not physically interesting.
4 Concluding remarks
One investigates the Mach-type (or (3142)-type) soliton of the Novikov–Veselov equation. The
Mach stem ([1,4]-soliton), the incident wave ([1,3]-soliton) and the reflected wave ([3,4]-soliton)
form a resonant triple. From (3.6), we see that the amplitude of Mach stem is less than two
times of the one of the incident wave, which is different from the KP equation [18]; moreover, the
length of the Mach stem is computed and show it is linear with time (2.6). On the other hand,
one uses the parameter κ (2.3) to describe the critical behavior for the O-type and Mach-type
solitons and notices that it depends on the the fixed parameter ε. We see that the amplitude A
of the incident wave is small than 2ε; furthermore, if ε < A < 2ε, then the soliton will be
singular. Now, a natural question is: what happens if A > 2ε when ε is fixed in (1.1)? Another
Mach-Type Soliton in the Novikov–Veselov Equation 13
question is the minimal completion [20]. It means the resulting chord diagram has the smallest
total length of the chords. This minimal completion can help us study the asymptotic solu-
tions and estimate the maximum amplitude generated by the interaction of those initial waves.
A numerical investigation of these issues will be published elsewhere.
Acknowledgements
The author thanks the referees for their valuable suggestions. This work is supported in part
by the National Science Council of Taiwan under Grant No. NSC 102-2115-M-606-001.
References
[1] Ablowitz M.J., Baldwin D.E., Nonlinear shallow ocean-wave soliton interactions on flat beaches, Phys.
Rev. E 86 (2012), 036305, 5 pages, arXiv:1208.2904.
[2] Athorne C., Nimmo J.J.C., On the Moutard transformation for integrable partial differential equations,
Inverse Problems 7 (1991), 809–826.
[3] Biondini G., Chakravarty S., Soliton solutions of the Kadomtsev–Petviashvili II equation, J. Math. Phys.
47 (2006), 033514, 26 pages, nlin.SI/0511068.
[4] Bogdanov L.V., Veselov–Novikov equation as a natural two-dimensional generalization of the Korteweg–de
Vries equation, Theoret. and Math. Phys. 70 (1987), 219–223.
[5] Chakravarty S., Kodama Y., Soliton solutions of the KP equation and application to shallow water waves,
Stud. Appl. Math. 123 (2009), 83–151, arXiv:0902.4423.
[6] Chakravarty S., Lewkow T., Maruno K.I., On the construction of the KP line-solitons and their interactions,
Appl. Anal. 89 (2010), 529–545, arXiv:0911.2290.
[7] Chang J.H., On the N -solitons solutions in the Novikov–Veselov equation, SIGMA 9 (2013), 006, 13 pages,
arXiv:1206.3751.
[8] Chang J.H., The interactions of solitons in the Novikov–Veselov equation, arXiv:1310.4027.
[9] Dubrovsky V.G., Topovsky A.V., Basalaev M.Y., New exact multi line soliton and periodic solutions with
constant asymptotic values at infinity of the NVN integrable nonlinear evolution equation via dibar-dressing
method, arXiv:0912.2155.
[10] Grinevich P.G., The scattering transform for the two-dimensional Schrödinger operator with a potential
that decreases at infinity at fixed nonzero energy, Russ. Math. Surv. 55 (2000), 1015–1083.
[11] Grinevich P.G., Manakov S.V., Inverse problem of scattering theory for the two-dimensional Schrödinger
operator, the ∂̄-method and nonlinear equations, Funct. Anal. Appl. 20 (1986), 94–103.
[12] Hu H.-C., Lou S.-Y., Construction of the Darboux transformaiton and solutions to the modified Nizhnik–
Novikov–Veselov equation, Chinese Phys. Lett. 21 (2004), 2073–2076.
[13] Hu H.-C., Lou S.-Y., Liu Q.-P., Darboux transformation and variable separation approach: the Nizhnik–
Novikov–Veselov equation, Chinese Phys. Lett. 20 (2003), 1413–1415, nlin.SI/0210012.
[14] Ishikawa M., Wakayama M., Applications of minor-summation formula. II. Pfaffians and Schur polynomials,
J. Combin. Theory Ser. A 88 (1999), 136–157.
[15] Kazeykina A.V., Novikov R.G., Large time asymptotics for the Grinevich–Zakharov potentials, Bull. Sci.
Math. 135 (2011), 374–382, arXiv:1011.4038.
[16] Kodama Y., KP solitons in shallow water, J. Phys. A: Math. Theor. 43 (2010), 434004, 54 pages,
arXiv:1004.4607.
[17] Kodama Y., KP solitons and Mach reflection in shallow water, arXiv:1210.0281.
[18] Kodama Y., Lectures delivered at the NSF/CBMS Regional Conference in the Mathematical Sciences “Soli-
tons in Two-Dimensional Water Waves and Applications to Tsunami” (UTPA, May 20–24, 2013), available
at http://faculty.utpa.edu/kmaruno/nsfcbms-tsunami.html.
[19] Kodama Y., Maruno K.-I., N -soliton solutions to the DKP equation and Weyl group actions, J. Phys. A:
Math. Gen. 39 (2006), 4063–4086, nlin.SI/0602031.
[20] Kodama Y., Oikawa M., Tsuji H., Soliton solutions of the KP equation with V -shape initial waves,
J. Phys. A: Math. Theor. 42 (2009), 312001, 9 pages, arXiv:0904.2620.
http://dx.doi.org/10.1103/PhysRevE.86.036305
http://dx.doi.org/10.1103/PhysRevE.86.036305
http://arxiv.org/abs/1208.2904
http://dx.doi.org/10.1088/0266-5611/7/6/005
http://dx.doi.org/10.1063/1.2181907
http://arxiv.org/abs/nlin.SI/0511068
http://dx.doi.org/10.1007/BF01039213
http://dx.doi.org/10.1111/j.1467-9590.2009.00448.x
http://arxiv.org/abs/0902.4423
http://dx.doi.org/10.1080/00036810903403343
http://arxiv.org/abs/0911.2290
http://dx.doi.org/10.3842/SIGMA.2013.006
http://arxiv.org/abs/1206.3751
http://arxiv.org/abs/1310.4027
http://arxiv.org/abs/0912.2155
http://dx.doi.org/10.1070/rm2000v055n06ABEH000333
http://dx.doi.org/10.1007/BF01077263
http://dx.doi.org/10.1088/0256-307X/21/11/001
http://dx.doi.org/10.1088/0256-307X/20/9/301
http://arxiv.org/abs/nlin.SI/0210012
http://dx.doi.org/10.1006/jcta.1999.2988
http://dx.doi.org/10.1016/j.bulsci.2011.02.003
http://dx.doi.org/10.1016/j.bulsci.2011.02.003
http://arxiv.org/abs/1011.4038
http://dx.doi.org/10.1088/1751-8113/43/43/434004
http://arxiv.org/abs/1004.4607
http://arxiv.org/abs/1210.0281
http://faculty.utpa.edu/kmaruno/nsfcbms-tsunami.html
http://dx.doi.org/10.1088/0305-4470/39/15/012
http://dx.doi.org/10.1088/0305-4470/39/15/012
http://arxiv.org/abs/nlin.SI/0602031
http://dx.doi.org/10.1088/1751-8113/42/31/312001
http://arxiv.org/abs/0904.2620
14 J.H. Chang
[21] Kodama Y., Williams L., KP solitons, total positivity, and cluster algebras, Proc. Natl. Acad. Sci. USA 108
(2011), 8984–8989, arXiv:1105.4170.
[22] Kodama Y., Williams L., The Deodhar decomposition of the Grassmannian and the regularity of KP solitons,
Adv. Math. 244 (2013), 979–1032, arXiv:1204.6446.
[23] Kodama Y., Williams L., KP solitons and total positivity for the Grassmannian, Invent. Math. 198 (2014),
637–699, arXiv:1106.0023.
[24] Manakov S.V., The method of the inverse scattering problem, and two-dimensional evolution equations,
Russian Math. Surveys 31 (1976), no. 5, 245–246.
[25] Matveev V.B., Salle M.A., Darboux transformations and solitons, Springer Series in Nonlinear Dynamics,
Springer-Verlag, Berlin, 1991.
[26] Miles J.W., Resonantly interacting solitary waves, J. Fluid Mech. 79 (1977), 171–179.
[27] Nimmo J.J.C., Darboux transformations in (2 + 1)-dimensions, in Applications of Analytic and Geometric
Methods to Nonlinear Differential Equations (Exeter, 1992), NATO Adv. Sci. Inst. Ser. C Math. Phys. Sci.,
Vol. 413, Kluwer Acad. Publ., Dordrecht, 1993, 183–192.
[28] Novikov S.P., Veselov A.P., Two-dimensional Schrödinger operator: inverse scattering transform and evolu-
tional equations, Phys. D 18 (1986), 267–273.
[29] Ohta Y., Pfaffian solutions for the Veselov–Novikov equation, J. Phys. Soc. Japan 61 (1992), 3928–3933.
[30] Veselov A.P., Novikov S.P., Finite-gap two-dimensional potential Schrödinger operators. Explicit formulas
and evolution equations, Sov. Math. Dokl. 30 (1984), 588–591.
[31] Yeh H., Li W., Kodama Y., Mach reflection and KP solitons in shallow water, Eur. Phys. J. ST 185 (2010),
97–111, arXiv:1004.0370.
http://dx.doi.org/10.1073/pnas.1102627108
http://arxiv.org/abs/1105.4170
http://dx.doi.org/10.1016/j.aim.2013.06.011
http://arxiv.org/abs/1204.6446
http://dx.doi.org/10.1007/s00222-014-0506-3
http://arxiv.org/abs/1106.0023
http://dx.doi.org/10.1007/978-3-662-00922-2
http://dx.doi.org/10.1017/S0022112077000093
http://dx.doi.org/10.1007/978-94-011-2082-1_18
http://dx.doi.org/10.1016/0167-2789(86)90187-9
http://dx.doi.org/10.1143/JPSJ.61.3928
http://dx.doi.org/10.1140/epjst/e2010-01241-0
http://arxiv.org/abs/1004.0370
1 Introduction
2 Mach type soliton
3 Relations with V-shape initial value waves
4 Concluding remarks
References
|