Next Article in Journal
Verifiable (2, n) Image Secret Sharing Scheme Using Sudoku Matrix
Next Article in Special Issue
Rogue Waves Generator and Chaotic and Fractal Behavior of the Maccari System with a Resonant Parametric Forcing
Previous Article in Journal
Micromagnetic Simulation of Saturation Magnetization of Nanocrystalline Soft Magnetic Alloys under High-Frequency Excitation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Parametric Resonance for the Nonlocal Hirota–Maccari Equation

Department of Physics, Istituto Superiore Piazza Resistenza, Piazza della Resistenza 1, 00015 Monterotondo, RM, Italy
Symmetry 2022, 14(7), 1444; https://doi.org/10.3390/sym14071444
Submission received: 24 June 2022 / Revised: 10 July 2022 / Accepted: 12 July 2022 / Published: 14 July 2022
(This article belongs to the Special Issue Nonlinear Dynamics and Applied Partial Differential Equations)

Abstract

:
The nonlocal Hirota–Maccari equation is considered when a parametric excitation is acting over the frequency of a generic mode. Using the well-known asymptotic perturbation (AP) method, two coupled equations for the amplitude and phase can be obtained. We discovered the existence of an infinite-period bifurcation when the parametric force increases its value. Moreover, symmetry considerations suggest performing a global analysis of the two couples, in such a way that we find an energy-like function and corroborate and verify the existence of this infinite period bifurcation.

1. Introduction

Resonances are recurrent topics in many physics problems. Usually, this means surprisingly large feedback when a little external perturbation is acting over the system. The undamped harmonic oscillator driven by an external force is a celebrated example, which is textbook material. If we study this simple linear oscillator, then one easily proves the existence of unbounded linearly increasing terms in the solution [1].
Nonlinear systems can also exhibit resonant behavior. We can consider many Hamiltonian chaotic systems in such a way that torus resonances are the reason for the origin of chaos, to cite only a few examples [2]. A similar picture can be observed when partial differential equations are investigated. Here, some results are known in simple linear cases, but in general, the study of resonances in complex nonlinear systems and the connected partial differential equations is always a difficult question [3,4,5,6,7,8].
We can think that a resonant perturbation can only excite a small number of modes in a nonlinear partial differential equation (NPDE). Usually, the nonlinearities mix different modes and the growth of one of them is forbidden by energy transfer to the others, specifically in conservative systems where the conservation laws are important constraints for the number and amplitude of active modes. Therefore, if the energy is transferred to the non-resonant modes and in cases when their growth is forbidden by nonlinear constraints, the role of the perturbation should be limited, and unbounded growth of the relevant quantities should be inhibited. Obviously, this is the common mechanism in many systems. However, other interesting behaviors can occur, and unexpected phenomena are possible even in a supposedly simple physically relevant model.
It can be demonstrated in this paper that an infinite period bifurcation can occur in the nonlocal Hirota–Maccari equation with two spatial dimensions and one temporal variable [9], if we consider a parametric resonance, when the external excitation frequency is nearly double that of a generic mode frequency.
Let us consider the following Hirota–Maccari equation
i ψ t + ψ x y + i ψ x x x + ψ ϕ i | ψ | 2 ψ x = 0
3 ϕ x + ( | ψ | 2 ) y = 0
where ψ = ψ ( x , y , t ) is a suitable complex function and ϕ = ϕ ( x , y , t ) is a real function. The nonlocality is given by the real function ϕ = ϕ ( x , y , t ) ; its values depend on the complex function ψ = ψ ( x , y , t ) values all over the plane (x,y). The Hirota–Maccari equation substitutes the Kadomtsev–Petviashvili equation when we consider phenomena where the spatial scale in the x-direction is different from the spatial scale in the y-direction, and, in this case, the real part of the function ψ is the wave amplitude and the real function ϕ describes the nonlocal effects [10,11].
For x = y, we obtain the well-known one-dimensional Hirota equation [12] that describes the spread of the femto-second soliton pulse in the single-mode fibers.
Now, we introduce a parametric excitation that can induce a parametric resonance into Equation (1) and arrive to the following model system
i ψ t + ψ x y + i ψ x x x + ψ ϕ i | ψ | 2 ψ x = f ψ ˜ e x p ( i K 1 x + i K 2 y i Ω t )
3 ϕ x + ( | ψ | 2 ) y = 0
where f is the parametric resonance amplitude, and the circular frequency of a given mode is
ω = K 1 K 2 K 1 4
The external force is in parametric resonance ( Ω 2 ω ) with the frequency of a generic mode.
In Section 2, using the asymptotic perturbation (AP) method [13], we build an approximate solution to Equations (3) and (4) and obtain two coupled equations for the solution amplitude and phase modulations.
In Section 3, we demonstrate the existence of an infinite-period bifurcation. If we increase the parametric strength from very small values, the solution oscillates with its almost natural frequency. However, when the excitation reaches a critical value, the solution frequency decreases further and further. We arrive to an infinite period or zero-frequency bifurcation.
As a consequence, if we carefully choose the initial conditions, then the approximate solution of Equations (3) and (4) becomes immobile, frozen and insensitive to the parametric excitation.
Finally, for higher values of the amplitude response, the frequency increases and we observe a modulated motion.
In Section 4, using symmetry considerations, a global analysis of the model system we obtained in Section 2 can be performed, utilizing an energy-like function, and the existence of an infinite period bifurcation can be validated.
In Section 5, the paper’s main findings are reviewed, and possible developments and directions for future work are discussed.

2. The Approximate Solution

The nonlocal Hirota–Maccari Equations (3) and (4) in resonance with a parametric excitation are investigated in this section; the forcing term f is scaled as ε 2 f , where ε is our small parameter, and we introduce the slow time
τ = ε 2 t
in order to consider larger time scales and find the amplitude modulation response.
Using the AP method [13], the approximate solution ψ ( x , y , t ) of Equations (3) and (4) can be written in the following form
ψ ( x , y , t ) = ε Ψ ( τ ) e x p ( i α )
α = K x x + K y y Ω 2 t
and τ in Equation (7) is given by (6), where the boundary conditions are
ψ ( 0 , y , t ) = ψ ( L , y , t )   0 y L   K x = 2 n 1 π L ,   n 1   integer
ψ ( x , 0 , t ) = ψ ( x , L , t )   0 x L   K y = 2 n 2 π L ,   n 2   integer
The proximity of the excitation frequency Ω to the system frequency ω can be characterized through a detuning parameter σ,
ω n _ = Ω 2 + ε 2 σ ,
where
n _ = ( n 1 , n 2 ) ,   ω n _ = 4 ( n 1 n 2 π 2 L 2 ) 16 ( n 1 π ) 4 L 4 .
Note that only a single leading mode is accounted for by the approximate solution. We assume that the system is excited near the natural frequency of a specific linear mode, and that mode is not involved in an in other resonances with any other mode. Following the AP method, the modes that are not directly excited by an external source or indirectly through an internal resonance will decay with time, if friction terms are present. Note that the variable change (6) implies that
d d t ( Ψ e x p ( i Ω 2 t ) ) = ( i Ω 2 Ψ + ε 2 d Ψ d τ ) e x p ( i Ω 2 t ) .
We use the temporal rescaling (6) in order to find the asymptotic behavior of the solution, when the nonlinear amplitude can be modified by the nonlinear effects.
The elimination of the predominant linear part of Equations (3) and (4) is obtained by the assumed solution (7), and it allows us to calculate the amplitude modulation given by the parametric excitation.
After using solution (7) for the complete Equations (3) and (4), we obtain the linear equation
( ω n _ K 1 K 2 + K 1 4 ) Ψ = 0 ,   ( order   O ( ε 0 ) ) .
If we consider Equations (3) and (4) at the order, then we find that the function satisfies the nonlinear differential equation
i Ψ τ σ Ψ K 1 | Ψ | 2 Ψ = f Ψ ˜ Ψ Φ
Φ x = 0
The complex-valued function is characterized by its amplitude and phase
Ψ = ρ e x p ( i θ ) ,
and we obtain at the model equations
d ρ d τ + f ρ s i n ( 2 θ ) = 0 ,
d θ d τ + K 1 ρ 2 f c o s ( 2 θ ) + σ Φ = 0 ,
where Φ is a function of y and t, but for simplicity, we consider it as a constant. This model system (17, 18) describes the amplitude and phase evolution according to our starting assumptions. We observe that the validity of the approximate solution (7) should be expected to be kept under control on time-scale t = O ( 1 ε 2 ) and on bounded intervals of the τ-variable. On the contrary, the approximate solution (7) loses its validity if we want to find solutions on larger intervals such that τ = O ( 1 ε ) . In the next section, we will demonstrate the possibility of an infinite-period bifurcation for the nonlocal Hirota–Maccari equation.

3. A Reverse Infinite-Period Bifurcation

The fixed points of the model system (17, 18) are connected to periodic solutions of the starting system (3, 4), and we must consider the conditions d ρ / d t = d θ / d t = 0 . Note that the solution exists in the following cases
P 1 :   θ E = 0   ρ E = ( f σ + Φ ) K 1   f σ Φ ,
P 2 :   θ E = π 2   ρ E = ( f σ + Φ ) K 1   Φ σ f ,
We can easily determine their possible stability, and through the Jacobian matrix, we obtain the equation
λ 2 4 f K 1 ρ E 2 = 0
λ 2 = ± 4 K 1 f ρ E 2
where, in the last equation, the minus sign is for the solution P1 (elliptic point) and the plus sign is for the solution P2 (saddle point).
From Equations (21a) and (21b), we can find the small oscillation frequency around the elliptic point
Ω ˜ = 2 ρ E ( f K 1 )
Using Equations (17) and (18), we arrive at the frequency-response equation
σ = f + Φ K 1 ρ E 2 ,
where now the equilibrium point amplitude is given by (19, 20) (Figure 1).
Subsequently, we show the response-parametric excitation curve (Figure 2). For the elliptic point, the response increases with the parametric excitation.
We show that an infinite-period bifurcation occurs near the elliptic point, but it is more convenient if we consider increasing parametric excitation values.
We observe the following situation ( σ > 0 , the excitation frequency Ω is slightly smaller than the double of ω, near the equilibrium point P1):
(i) When f = 0, there is no parametric excitation and the wave amplitude is constant;
(ii) When the excitation is weak, there are no equilibrium points because of (3.1) and the wave amplitude begins to oscillate with an amplitude increasing with the excitation amplitude f and the frequency ω ˜
ω ˜ = ω K 1 ρ 2 + Φ σ ;
(iii) If the external excitation increases and reaches the critical value
f C = σ Φ ,
the oscillation frequency decreases and the wave seems to collapse, but actually it begins a very slow oscillation, meaning metaphorically that it dies and is born again, with a large period, and an infinite-period bifurcation occurs when f f C (Figure 3).
We can estimate the oscillation period T assuming in our calculations
T = 2 π Ω ˜ = 0 2 π ( d θ ) K 1 ρ E 2 σ + Φ + f c o s ( 2 θ ) = π ( 2 ) ( f c ) ( f f c )
(iv) If the initial condition for the wave amplitude is carefully tuned, that is, it corresponds to P1, then the wave amplitude seems immobile, frozen and insensitive to the external excitation with no oscillation;
(v) Suddenly, beyond the critical value (25) and for higher response amplitude values, a slow oscillation appears with the frequency
Ω ˜ = ρ E ( K 1 )

4. Global Analysis of the Model System

A global analysis of the model system (17, 18) can be performed. Considering the model system (17, 18), we obtain an energy-like function E ( ρ , θ )
E ( ρ , θ ) = K 1 2 ρ 4 + ( σ Φ ) ρ 2 + f ρ 2 c o s ( 2 θ ) d E d τ = 0
We can easily conclude that every function in the form λ E ( ρ , θ ) , with λ real numbers, can be used as an energy function. We can take λ = 1 so that the elliptic point P1 corresponds to a minimum. In Figure 4, we show the level curves for the energy function (28) near the equilibrium point (19), and we can easily observe that the infinite-period bifurcation is connected to the “bottleneck” we can see in the same figure. The bottleneck is, when in the closed curves in Figure 4, the ρ-value (y-axis) is constant and ϑ (x-axis) is varying.
The fixed point (19) downshifts the solution motion near it, and a lower frequency is born.

5. Conclusions

The nonlocal Hirota–Maccari equation with an external periodic force in parametric resonance with a generic mode has been investigated. An approximate solution with the AP perturbation method, previously used for other nonlinear partial differential equations, was suggested and analyzed qualitatively in the vicinity of critical point fC. Two coupled nonlinear equations describe the temporal evolution for the solution amplitude and phase. Frequency-response curves are shown and the existence of an infinite-period bifurcation (when the parametric excitation increases in magnitude) can be demonstrated.
Looking at Equations (19), (22) and (26), it is possible to discern certain analogies with the Hopf–Landau theory of turbulence [1,14].
This perturbation method could be applied to important resonances for other nonlinear partial differential equations in plasma physics and/or nonlinear optics.

Funding

This research received no external funding.

Data Availability Statement

Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.

Conflicts of Interest

The author declares that there are no conflicts of interest.

References

  1. Fossen, T.; Nijmeijer, H. (Eds.) Parametric Resonance in Dynamical Systems; Springer: New York, NY, USA, 2011. [Google Scholar]
  2. Kim, S.; Mackay, R.S.; Guckenheimer, J. Resonance regions for families of torus maps. Nonlinearity 1989, 2, 391–404. [Google Scholar] [CrossRef]
  3. Sahoo, S.K.; Das, H.C.; Panda, C.L. An overview of transverse vibration of axially travelling string. In Recent Trends in Applied Mathematics; Mishra, S.R., Dhamala, T.N., Makinde, O.D., Eds.; Springer: New York, NY, USA, 2021; pp. 427–446. [Google Scholar]
  4. Jackson, R.K.; Weinstein, M.I. Geometric Analysis of Bifurcation and Symmetry Breaking in a Gross–Pitaevskii Equation. J. Stat. Phys. 2004, 116, 881–905. [Google Scholar] [CrossRef] [Green Version]
  5. Susanto, H.; Cuevas, J.; Kruger, P. Josephson tunnelling of dark solitons in a double-well potential. J. Phys. B At. Mol. Opt. Phys. 2011, 44, 059003. [Google Scholar] [CrossRef] [Green Version]
  6. Susanto, H.; Cuevas, J. Josephson Tunneling of Excited States in a Double-Well Potential, in Sponaneous Symmetry Breaking. In Self-Trapping, and Josephson Oscillations; Springer: Berlin/Heidelberg, Germany, 2012; pp. 583–599. [Google Scholar]
  7. Marangell, R.; Jones, R.T.; Susanto, H. Localized standing waves in inhomogenous Schrödinger equations. Nonlinearity 2010, 23, 2059. [Google Scholar] [CrossRef]
  8. Marangell, R.; Jones, R.T.; Susanto, H. Instability of standing waves for nonlinear Schrödinger-type equations. J. Differ. Equ. 2012, 253, 1191–1205. [Google Scholar] [CrossRef] [Green Version]
  9. Maccari, A. A generalized Hirota equation in 2+1 dimensions. J. Math. Phys. 1998, 39, 6547–6551. [Google Scholar] [CrossRef]
  10. Demiray, S.T.; Pandir, Y.; Bulut, H. All exact travelling wave solutions of Hirota equation and Hirota–Maccari system. Optik 2016, 127, 1848–1859. [Google Scholar] [CrossRef]
  11. Xia, P.; Zhang, Y.; Zhang, H.; Zhuang, Y. Some novel dynamical behavior of localized solitary waves for the Hirota-Maccari system. Nonlinear Dyn. 2022, 108, 533–541. [Google Scholar] [CrossRef]
  12. Hirota, R. Exact envelope-soliton solutions of a nonlinear wave equation. J. Math. Phys. 1973, 14, 805–809. [Google Scholar] [CrossRef]
  13. Maccari, A. Coherent solutions for the fundamental resonance of the Boussinesq equation. Chaos Solitons Fractals 2013, 54, 57–64. [Google Scholar] [CrossRef]
  14. Gluzman, S.; Karpeev, D. Perturbative expansions and critical phenomenain random structured media. In Modern Problems in Applied Analysis; Drygaś, P., Rogosin, S., Eds.; Birkhauser: Basel, Switzerland, 2018; pp. 117–134. [Google Scholar]
Figure 1. Frequency (y = σ)-response (x = ρ) curve where K1 = 0.9, f = 0.02, φ = 0.03. The upper branch is the elliptic point, the lower branch corresponds to the saddle point.
Figure 1. Frequency (y = σ)-response (x = ρ) curve where K1 = 0.9, f = 0.02, φ = 0.03. The upper branch is the elliptic point, the lower branch corresponds to the saddle point.
Symmetry 14 01444 g001
Figure 2. Response (y = ρ)-parametric excitation (x = f) curve where K1 = 0.9, σ = 0.02, φ = 0.03. The upper branch is the elliptic point, the lower branch corresponds to the saddle point.
Figure 2. Response (y = ρ)-parametric excitation (x = f) curve where K1 = 0.9, σ = 0.02, φ = 0.03. The upper branch is the elliptic point, the lower branch corresponds to the saddle point.
Symmetry 14 01444 g002
Figure 3. Infinite-period bifurcation. We show the period T = y-axis and the parametric excitation f = x-axis, when fC is 0.025 or fC = 0.023.
Figure 3. Infinite-period bifurcation. We show the period T = y-axis and the parametric excitation f = x-axis, when fC is 0.025 or fC = 0.023.
Symmetry 14 01444 g003
Figure 4. Level curves in the plane x = ϑ and y = ρ. (K1 = 0.6, f = 0.02, σ = φ). Obviously, only positive values for ρ are acceptable. From the upper side, we show level curves for energy E = 0.0002, 0.0001 and 0 and we can observe a modulated motion. Subsequently, we obtain closed curves for E = −0.00001, −0.00007, −0.0001.
Figure 4. Level curves in the plane x = ϑ and y = ρ. (K1 = 0.6, f = 0.02, σ = φ). Obviously, only positive values for ρ are acceptable. From the upper side, we show level curves for energy E = 0.0002, 0.0001 and 0 and we can observe a modulated motion. Subsequently, we obtain closed curves for E = −0.00001, −0.00007, −0.0001.
Symmetry 14 01444 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Maccari, A. A Parametric Resonance for the Nonlocal Hirota–Maccari Equation. Symmetry 2022, 14, 1444. https://doi.org/10.3390/sym14071444

AMA Style

Maccari A. A Parametric Resonance for the Nonlocal Hirota–Maccari Equation. Symmetry. 2022; 14(7):1444. https://doi.org/10.3390/sym14071444

Chicago/Turabian Style

Maccari, Attilio. 2022. "A Parametric Resonance for the Nonlocal Hirota–Maccari Equation" Symmetry 14, no. 7: 1444. https://doi.org/10.3390/sym14071444

APA Style

Maccari, A. (2022). A Parametric Resonance for the Nonlocal Hirota–Maccari Equation. Symmetry, 14(7), 1444. https://doi.org/10.3390/sym14071444

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