Abstract
In this paper, our aim is to further expand the use of the two-variable -expansion approach to a new coupled KdV and Z-K system, which has various significant applications in different fields of applied sciences. The KdV equation, along with shallow-water waves and long internal waves in oceans, basically explains how long, one-dimensional waves propagate in a variety of physical conditions. The study of coastal waves on the basis of the ocean is done using the Zakharov–Kuznetsov (Z-K) equation and this model is utilized to illustrate ion-acoustic wave propagation. By using this method, different forms of analytical solutions of the new coupled KdV (NCKdV) system and the new coupled Z-K (NCZ-K) system, such as solitons, multi-peak solitons, solitary waves, trigonometric, hyperbolic and rational functions and other wave solutions are constructed. The significant features of multi-peak solitons induced by the higher-order effects, including velocity variations, localization or periodicity attenuation and state transitions, are revealed. When the localization disappears then the multi-peak soliton becomes a periodic wave. The constructed solutions are also presented graphically having their applications in engineering, etc. The stability of the solution is examined by utilizing modulation instability. The results obtained show that the proposed technique is universal and efficient. In addition, this technique can also be applied to lots of other new coupled systems arising in other areas of applied sciences.
1. Introduction
Non-linear phenomena play a vital part in applied sciences. Soliton theory heavily depends on the travelling wave solutions of non-linear models in applied sciences [Citation1–4]. The generalized KdV and Zakharov–Kuznetsov (Z-K) equations are essential models for many physical phenomena in partial differential equations (PDEs), including shallow-water waves, non-linear optics, applications in astrophysical and space environments, the interaction between a water wave and a floating ice cover, gravity-capillary waves, shock waves, ion-acoustic waves, hydro-dynamic, stratified waves and many others. The study of explicit analytical solutions to new coupled systems is crucial for illuminating the properties of non-linear phenomena and soliton theory. The travelling wave solutions of these non-linear coupled systems and high-dimensional non-linear equations may enable researchers to thoroughly investigate the aforementioned natural phenomena. Therefore, it is very important to build up non-linear evolution equations and analytical solutions of new coupled systems.
In mathematical physics and the applied sciences, the study of analytical solutions in the explicit form for new coupled systems is crucial for describing the characteristics of non-linear issues and soliton theory. To gain explicit results in the form of solitons, solitary wave and elliptic function solutions, many precise and influential techniques have been developed such as Inverse scattering transform technique [Citation5], (-expansion method [Citation6], the Classical Darboux transformations [Citation7,Citation8], Bäcklund transformation [Citation9–11], expansion scheme, Hirota's bilinear [Citation12,Citation13], Painlevé analysis [Citation14,Citation15], extended tanh scheme [Citation16], auxiliary equation scheme [Citation17], direct algebraic methods [Citation18,Citation19], the mapping and extended mapping methods [Citation20], elliptic function scheme [Citation21], rational expansion method and several others [Citation22,Citation23]. A fantastic progress has been made in the analysis of the solutions, structures, interactions and other properties of solitary waves and solitons, also distinct profound results have been effectively obtained [Citation24,Citation25]. Explicit and numerical solutions were illustrated of KdV and coupled KdV equations by decomposition approaches [Citation26–28]. ‘ The concept of soliton was initiated to show that the solitary waves, which keep their characters unaffected, throughout the proliferation process and after interface, are like particles [Citation29–31]. Because of these characteristics, solitons are believed as standard data bits [Citation32]. Owing to the potential applications of optical solitons in telecommunications and ultrafast signal processing systems, it is the main target of theoretical and experimental analysis [Citation32,Citation33]. As the exact results of the PDEs, the breather plays a central task in semi-conductor quantum wells [Citation34], arrays of micro-mechanical oscillators [Citation35] and junctions of Josephson [Citation36]. Breather specifies that the solutions’ behaviour is periodic in space or time having the property of spatial or temporal localization. Two types of known breathers are the Akhmediev breather of space-periodic and the Kuznetsov-Ma soliton of time-periodic [Citation37–42]. Some other analytical results are the results of rogue waves [?, Citation42–46]. These rogue wave solutions are confined in together space and time, arise from nowhere and vanish starved of a trace [Citation43,Citation44]. These results may be acquired through the Taylor expansion of the breather results [Citation47–50].
In this paper, the two-variable -expansion approach is used to achieve the analytical solutions of the new coupled KdV and Z-K systems. As a result, distinct types of solitons and other wave solutions are obtained. The stability of the solution is examined by utilizing modulation instability (MI). The significant features of multi-peak solitons induced via the higher-order effects, including velocity variations, localization or periodicity attenuation and state transitions, are revealed. When the localization disappears, then the multi-peak soliton becomes a periodic wave. The main advantage of this method is that different forms of analytical solutions of such as solitons, solitary waves, rational solitons, trigonometric function, hyperbolic function and other wave solutions are constructed using this technique. The proposed technique is also universal and efficient. This technique can also be applied to lots of other new coupled systems arising in other areas of applied sciences.
The remaining article's main structure is as follows: Section 2 describes the two-variable technique in detail. In the third Section, the given technique is applied on the NCKdV and NCZ-K systems for producing precise wave results. The analysis of modulation instability is discussed in the fourth section to explain the stability of models. The results discussion and its physical justification are provided in the fifth section. Finally, the sixth section summarizes the entire work.
2. Proposed scheme
In order to arrive at wave results of the aforementioned equations, the -expansion technique is described in detail in this section. The following is a full introduction to this strategy from [Citation1,Citation51]:
The aforementioned approach begins with the second-order linear ODE, which is known as follows: (1) (1) by taking (2) (2) for the accuracy of the calculations. The variables Φ and Ψ have the derivatives as (3) (3) There are three types of cases in the general results of Equation (Equation1(1) (1) )
Case-I: , the solution is (4) (4) with (5) (5) Case-II: , the solution is as (6) (6) with (7) (7) Case-III: , the rational solution is as (8) (8) with (9) (9) In all cases, and are constants.
The following procedures must be taken in order to use this method to obtain precise solutions to the non-linear evolution equation. (10) (10) Step 1: Wave transformation is used to convert Equation (Equation10(10) (10) ) to ODE (ordinary differential equation); it includes as (11) (11) where is the nth order derivative with respect to ξ.
Step 2: The result could be in the form of Ψ and Φ (12) (12) here G adheres to (Equation1(1) (1) ). constant coefficients and α are to be determined. N can be determined by using the balancing principle in (Equation11(11) (11) ).
Step 3: By considering case 1 as an example, utilizing (Equation12(12) (12) ) into (Equation11(11) (11) ) while taking into account (Equation3(3) (3) ) and (Equation5(5) (5) ), a polynomial equation in Ψ and Φ is attained and produces a system of algebraic equations.
Step 4: A software package program is used to solve the system. Wave solutions in (Equation11(11) (11) ) are constructed as three different types of functions by utilizing resultant values and .
Step 5: The solution process is completed by producing results in (Equation10(10) (10) ) using (wave transformation) conversely.
3. Applications
3.1. The new coupled KdV system
Qin [Citation40] used a finite-dimensional integrable system to build a new hierarchy of non-linear evolution equations. A novel coupled KdV system [Citation52] is an interesting equation in this hierarchy is as (13) (13) where and γ are arbitrary constants. The author in [Citation17] more extended the NCKdV system into a new coupled KP system and deliberated both systems using Hirota's bilinear approach. Consider the travelling wave solution as (14) (14) where δ and ρ are wave number and frequency of results, respectively. Using (Equation14(14) (14) ) in (Equation13(13) (13) ) and the ODEs obtained as (15) (15) Using the balancing principle on (Equation15(15) (15) ) and considering the solutions as (16) (16) Utilizing Equations (Equation16(16) (16) ) with (Equation3(3) (3) ) into (Equation15(15) (15) ), we obtained system of equations in and by equating the coefficients of to zero. After solving the system, the following families of results are produced:
Family 1: (17) (17) Family 2: (18) (18) Family 3: Family 4: Case-I: , from the first family, the following soliton solution set is obtained in the hyperbolic function form as (19) (19) from the second family, the following soliton solution set is also obtained in the hyperbolic function form as (20) (20)
Case-II: , from the first and second families, the following soliton solution sets are obtained in the trigonometric function form as (21) (21) (22) (22) Case-III: , from the first and second families, thefollowing soliton solution sets are obtained in the rational function form as (23) (23) (24) (24) The more generalized resemble solitons of Equation (Equation13(13) (13) ) other families can also be constructed.
3.2. The new coupled Z-K system
The coupled KdV system (Equation13(13) (13) ) can extend to the NCZ-K system [Citation52] in the following form (25) (25) where and γ are real numbers. The authors in [Citation52] studied this system using the modified extended direct algebraic scheme. By considering the wave transformation as (26) (26) where and ν are the wave numbers and frequency, respectively. Now Equation (Equation25(25) (25) ) is transformed into ODE by using Equation (Equation25(25) (25) ) as (27) (27) On Equation (Equation27(27) (27) ) using the balancing principle, and presuming solutions as (28) (28) Utilizing Equations (Equation28(28) (28) ) with (Equation3(3) (3) ) into (Equation27(27) (27) ), and the system of algebraic equations obtained in and byequating the coefficients of to zero, the following families of solutions by solving these equations are obtained as follows:
Family 1: Family 2: Family 3: Family 4: Case-1: , from the first family, the following solution set is obtained in the hyperbolic function form as
(29) (29) From the second family, the following solution set is obtained in the hyperbolic function form as
(30) (30) Case-2: , from the first and second families, the following solution sets are obtained in the trigonometric function form as
(31) (31) (32) (32) Case-III: , from the first and second families, the following solution sets are obtained in the rational function form as
(33) (33) (34) (34) The more generalized resemble solitons of Equation (Equation25(25) (25) ) other families can also be constructed.
4. Modulation stability analysis
Various non-linear PDEs exhibit scattering and non-linear consequences as a result of uncertainty in the steady-state modulation. Utilizing linear stability analysis, the modulation instability of NCKdV and NCZ-K systems is investigated [Citation36,Citation38].
4.1. The new coupled KdV system
The following is the form of steady-state solution (SSS) of the NCKdV system (35) (35) Here, normalized optical power P. Evolution of the perturbation is verified by using the linear stability (LS) analysis. We can linearize by putting Equations (Equation35(35) (35) ) into (Equation13(13) (13) ), we have (36) (36) It is supposed that the solution of Equation (Equation36(36) (36) ) has as (37) (37) where k and ω are the wave number and frequency of perturbation, respectively. Utilizing Equations (Equation37(37) (37) ) into (Equation36(36) (36) ), the dispersion relation (DR) is acquired as (38) (38) The dispersion relations in (Equation38(38) (38) ) indicate that steady-state stability varies upon the wave number, self-phase modulation and stimulating Raman scattering. For wave numbers k, the velocity dispersion ω is real and the steady state is stable along the small perturbation if . when , it becomes unstable, i.e. ω is imaginary as the perturbation builds exponentially. It can be easily observed MI(modulation instability) when . These conditions indicate the increased rate of MI achieved spectrum can be stated as (39) (39)
4.2. The new coupled Z-K system
The following form is the SSS of the NCZ-K system (40) (40) Here, the normalized optical power is P. Evolution of the perturbation is verified by using the LS analysis. We are using Equations (Equation25(25) (25) ) and (Equation40(40) (40) ) and linearizing (41) (41) It is supposed that the solution of Equation (Equation41(41) (41) ) has as (42) (42) where ν and are the frequency of perturbation and wave numbers, respectively. By using Equations (Equation42(42) (42) ) and (Equation41(41) (41) ), the following dispersion relation is obtained (43) (43) The dispersion relations in (Equation43(43) (43) )) indicate that the steady-state stability varies upon the wave number, self-phase modulation and stimulating Raman scattering. For all wave numbers , the velocity dispersion ν is real and the steady state is stable along the small perturbation if . It turns out to be unstable if , i.e. ν is imaginary as the perturbation builds exponentially. It can easily be observed modulational instability (MI) if . These conditions indicate that the increased rate of MI gain spectrum can be conveyed as (44) (44)
5. Discussion and physical interpretation of results
The acquired results in this article are dissimilar of the gained results of several academics because Equation (Equation1(1) (1) ) is dissimilar from the existing methods. By assigning precise values of parameters, different families of solutions have been obtained for ODE (Equation1(1) (1) ). The new coupled KdV system was resolved by Hirota's bilinear approach [Citation17]. In [Citation52], modified extended direct algebraic approach has been utilized to get precise wave results of (Equation13(13) (13) ) and (Equation25(25) (25) ) in travelling wave solutions in trignometric and elliptic functions forms. In [Citation53], the researchers found Soliton solutions for a KdV equation and a generalized Hirota–Satsuma CKdV equation. By a natural decomposition method, a new solution of coupled KdV equation has been obtained in [Citation54]. Sine-Gordon and Modified Kudryashov methods were used to evaluate dual-mode Hirota-Satsuma coupled KdV equations in [Citation55]. The main advantage of this method is that different forms of analytical solutions of such as solitons, solitary waves, rational solitons, trigonometric function, hyperbolic function and other wave solutions are constructed using this technique. If or in hyperbolic function solutions, then we can obtain an analytical one-soliton solution; otherwise, we can get two soliton solutions. So, in this work, several innovative results have been accomplished, which never have been stated earlier.
In Figure , by specifying suitable values of parameters, the wave solution (Equation19(19) (19) ) in disparate constitution is represented as follows: Figure (A,B,E) describes the bright soliton and their 2D contour plot figures (Figure (B,D ,F)) correspondingly. In Figure , by specifying proper values of parameters, the wave solution (Equation20(20) (20) ) is shown in disparate forms as follows: Figure (A,C,E) are solions in multi-peak forms of dissimilar amplitude and their 2D contour plots in figures (Figure (B,D ,F)) correspondingly. In Figure , by specifying suitable values of parameters, the solution (Equation21(21) (21) ) is shown in different forms as follows: Figure (A,C,E) are solitons in periodic form of dissimilar amplitude and their 2D contourplots in Figure (B,D ,F), correspondingly.
In Figure , by specifying suitable values of parameters, the wave solution (Equation30(30) (30) ) is illustrated in disparate shapes as follows: Figure (A,C,E) are solitons in multi-peak forms and their 2D contourplot in Figure (B,D ,F) correspondingly. In Figure , by specifying suitable parameter values, the solution (Equation32(32) (32) ) is illustrated in various shapes as follows: Periodic solitons with varying amplitudes are shown in Figure (A,C, E), and their 2D contourplots are shown in Figure (B,D, F), respectively. Figure (A) depicts the DR between the ω and k of (Equation38(38) (38) ) and the DR between ν and wave numbers () of (Equation43(43) (43) ) is demonstrated in Figure (B), respectively.
6. Conclusion
We effectively applied the proposed technique to the NCKdV and the NCZ-K systems in this work. The KdV equation, along with shallow-water waves and long internal waves in oceans, basically explains how long, one-dimensional waves propagate in a variety of physical conditions. The study of coastal waves on the basis of the ocean is done using the Zakharov–Kuznetsov (Z-K) equation and this model is utilized to illustrate ion-acoustic wave propagation. As a result, different forms of analytical solutions of the new coupled KdV (NCKdV) system and new coupled Z-K (NCZ-K) system, such as solitons, multi-peak solitons, solitary waves, trigonometric, hyperbolic and rational function solutions and other wave solutions, are constructed in the explicit form using the proposed scheme. By giving the parameters, the appropriate values, novel structures of constructed solutions of these models are represented. Graphical representations of the physical structures of a few obtained results are efficient for conveying the intricate physical existences of both models. The stability of the solution is examined by utilizing modulation instability. The results obtained show that the proposed technique is universal and efficient, and can be applied on lots of other new coupled systems.
Disclosure statement
No potential conflict of interest was reported by the author(s).
References
- Shehzad K, Seadawy AR, Wang J, et al. Multi peak solitons and btreather types wave solutions of unstable NLSEs with stability and applications in optics. Opt Quantum Electron. 2023;55(1):7.
- Qian X, Lu D, Arshad M, et al. Novel traveling wave solutions and stability analysis of perturbed Kaup-Newell Schrödinger dynamical model and its applications. Chinese Phys B. 2021;30(2):020201.
- Zhang Z, Li B, Chen J, Construction of higher-order smooth positons and breather positons via Hirota's bilinear method. Nonlinear Dyn. 2021;105(3):2611–2618.
- Raza N, Jannat N, Gómez-Aguilar JF, et al. New computational optical solitons for generalized complex Ginzburg–Landau equation by collective variables. Modern Phys Lett B. 2022:2250152. doi:10.1142/S0217984922501524.
- Ablowitz MJ, Musslimani ZH. Inverse scattering transform for the integrable nonlocal nonlinear Schrödinger equation. Nonlinearity. 2016;29(3):915.
- Wang M, Li X, Zhang J. The ( G′/G)−expansion method and travelling wave solutions of nonlinear evolution equations in mathematical physics. Phys Lett A. 2008;372(4):417–423.
- Lambert F, Loris I, Springael J. Classical Darboux transformations and the KP hierarchy. Inverse Probl. 2001;17(4):1067.
- Levi D, Ragnisco O, Sym A. Dressing methods, classical Darboux transformation. Il Nuovo Cimento B (1971–1996). 1984;83(1):34–42.
- Rogers C, Schief WK. The classical Bäcklund transformation and integrable discretisation of characteristic equations. Phys Lett A. 1997;232(3–4):217–223.
- Biswas A, Konar S. Introduction to non-Kerr law optical solitons. New York: Chapman and Hall/CRC; 2006.
- Atamurotov F, Abdujabbarov A, Ahmedov B. Shadow of rotating non-Kerr black hole. Phys Rev D. 2013;88(6):064004.
- Zhang Z, Yang X, Li B, et al. Multi-lump formations from lump chains and plane solitons in the KP1 equation. Nonlinear Dyn. 2023;111:1625–1642.
- Zhang Z, Li B, Chen J, et al. Degenerate lump interactions within the Kadomtsev–Petviashvili equation. Commun Nonlinear Sci Numerical Simul. 2022;112:106555.
- Kumar S, Mohan B, Kumar A. Generalized fifth-order nonlinear evolution equation for the Sawada-Kotera, Lax, and Caudrey-Dodd-Gibbon equations in plasma physics: Painlevé analysis and multi-soliton solutions. Phys Scripta. 2022;97(3):035201.
- Islam M, Akter M, Gómez-Aguilar JF, et al. A novel study of the nonlinear Kadomtsev-Petviashvili-modified equal width equation describing the behavior of solitons. Opt Quantum Electron. 2022;54(11):1–15.
- Wu J. New explicit traveling wave solutions for three nonlinear evolution equations. Appl Math Comput. 2010;217(4):1764–1770.
- Wazwaz AM. Completely integrable coupled KdV and coupled KP systems. Commun Nonlinear Sci Numerical Simul. 2010;15(10):2828–2835.
- Abdullah, Seadawy AR, Jun W. Mathematical methods and solitary wave solutions of three-dimensional Zakharov-Kuznetsov-Burgers equation in dusty plasma and its applications. Res Phys. 2017;7:4269–4277.
- Arshad M, Seadawy AR, Lu D. Travelling wave solutions of Drinfel'd-Sokolov-Wilson, Whitham-Broer-Kaup and (2 + 1)-dimensional Broer-Kaup-Kupershmit equations and their applications. Chin J Phys. 2017;55(3):780–797.
- Arshad M, Seadawy AR, Lu D. Novel solitons and elliptic function solutions of (1 + 1)-dimensional higher order nonlinear Schrödinger equation with derivative non-Kerr nonlinear terms and its applications. Modern Phys Lett B. 2019;33(22):1950253.
- Helal MA, Seadawy AR. Exact soliton solutions of a D-dimensional nonlinear Schrödinger equation with damping and diffusive terms. Zeitschrift für Angewandte Mathematik und Physik. 2011;62(5):839–847.
- Seadawy AR, El-Rashidy K. Traveling wave solutions for some coupled nonlinear evolution equations. Math Comput Model. 2013;57(5-6):1371–1379.
- Seadawy AR. Approximation solutions of derivative nonlinear Schrödinger equation with computational applications by variational method. Eur Phys J Plus. 2015;130(9):1–10.
- Kaya D, Inan IE. Exact and numerical traveling wave solutions for nonlinear coupled equations using symbolic computation. Appl Math Comput. 2004;151(3):775–787.
- Kaya D. On the solution of a Korteweg-de Vries like equation by the decomposition method. Int J Comput Math. 1999;72(4):531–539.
- Kaya D, Aassila M. An application for a generalized KdV equation by the decomposition method. Phys Lett A. 2002;299(2-3):201–206.
- Nonlaopon K, Mann N, Kumar S, et al. A variety of closed-form solutions, Painlevé analysis, and solitary wave profiles for modified KdV–Zakharov–Kuznetsov equation in (3 + 1)-dimensions. Res Phys. 2022;36:105394.
- Kumar S, Dhiman SK. Lie symmetry analysis, optimal system, exact solutions and dynamics of solitons of a (3 + 1)-dimensional generalised BKP–Boussinesq equation. Pramana. 2022;96(1):1–20.
- Ablowitz MJ, Ablowitz MA, Clarkson PA. Solitons, nonlinear evolution equations and inverse scattering. Cambridge: Cambridge University Press; 1991.
- Gaber AA, Alsharari F, Kumar S. Some closed-form solutions, conservation laws, and various solitary waves to the (2 + 1)-D potential BK equation via lie symmetry approach. Int J Modern Phys B. 2022;36(20):2250117.
- Zhang CR, Tian B, Wu XY, et al. Rogue waves and solitons of the coherently-coupled nonlinear Schrödinger equations with the positive coherent coupling. Phys Scripta. 2018;93(9):095202.
- Serkin VN, Hasegawa A, Belyaeva TL. Solitary waves in nonautonomous nonlinear and dispersive systems: nonautonomous solitons. J Mod Opt. 2010;57(14-15):1456–1472.
- Wu XF, Hua GS, Ma ZY. Evolution of optical solitary waves in a generalized nonlinear Schrödinger equation with variable coefficients. Nonlinear Dyn. 2012;70(3):2259–2267.
- Haas F, Manfredi G, Shukla PK, et al. Breather mode in the many-electron dynamics of semiconductor quantum wells. Phys Rev B. 2009;80(7):073301.
- Sato M, Hubbard BE, Sievers AJ, et al. Observation of locked intrinsic localized vibrational modes in a micromechanical oscillator array. Phys Rev Lett. 2003;90(4):044102.
- Trias E, Mazo JJ, Orlando TP. Discrete breathers in nonlinear lattices: experimental detection in a josephson array. Phys Rev Lett. 2000;84(4):741.
- Zhang Z, Yang X, Li B, et al. Generation mechanism of high-order rogue waves via the improved long-wave limit method: NLS case. Phys Lett A. 2022;450:128395.
- Ablowitz MJ, Kaup DJ, Newell AC, et al. Method for solving the sine-Gordon equation. Phys Rev Lett. 1973;30(25):1262.
- Kibler B, Fatome J, Finot C, et al. Observation of Kuznetsov-Ma soliton dynamics in optical fibre. Sci Rep. 2012;2(1):1–5.
- Qin Z. A finite-dimensional integrable system related to a new coupled KdV hierarchy. Phys Lett A. 2006;355(6):452–459.
- Seadawy AR, Lu D. Bright and dark solitary wave soliton solutions for the generalized higher order nonlinear Schrödinger equation and its stability. Res Phys. 2017;7:43–48.
- Sun WR, Tian B, Zhen HL, et al. Breathers and rogue waves of the fifth-order nonlinear Schrödinger equation in the Heisenberg ferromagnetic spin chain. Nonlinear Dyn. 2015;81(1):725–732.
- Solli DR, Ropers C, Koonath P, et al. Optical rogue waves. Nature. 2007;450(7172):1054–1057.
- Akhmediev N, Ankiewicz A, Soto-Crespo JM. Rogue waves and rational solutions of the nonlinear Schrödinger equation. Phys Rev E. 2009;80(2):026601.
- Bludov YV, Konotop VV, Akhmediev N. Matter rogue waves. Phys Rev A. 2009;80(3):033610.
- Moslem WM, Shukla PK, Eliasson B. Surface plasma rogue waves. EPL (Europhysics Letters). 2011;96(2):25002.
- Akhmediev N, Ankiewicz A, Taki M. Waves that appear from nowhere and disappear without a trace. Phys Lett A. 2009;373(6):675–678.
- Xie XY, Tian B, Wang YF, et al. Rogue wave solutions for a generalized nonautonomous nonlinear equation in a nonlinear inhomogeneous fiber. Ann Phys (N Y). 2015;362:884–892.
- Tao Y, He J. Multisolitons, breathers, and rogue waves for the Hirota equation generated by the Darboux transformation. Phys Rev E. 2012;85(2):026601.
- Shan S, Li C, He J. On rogue wave in the Kundu-DNLS equation. Commun Nonlinear Sci Numerical Simul. 2013;18(12):3337–3349.
- Kara S, Ünsal Ö. Analytical solutions to new forms of two nonlinear partial differential equations via two variable expansion method. Partial Differ Equations Appl Math. 2022;5:100210.
- Seadawy AR, Arshad M, Lu D. Stability analysis of new exact traveling-wave solutions of new coupled KdV and new coupled Zakharov-Kuznetsov systems. Eur Phys J Plus. 2017;132(4):1–19.
- Fan E. Soliton solutions for a generalized Hirota–Satsuma coupled KdV equation and a coupled MKdV equation. Phys Lett A. 2001;282(1-2):18–22.
- Elbadri M, Ahmed SA, Abdalla YT. A new solution of time-fractional coupled KdV equation by using natural decomposition method. Abstr Appl Anal. 2020;2020:9.
- Alquran M, Jaradat I, Baleanu D. Shapes and dynamics of dual-mode Hirota-Satsuma coupled KdV equations: exact traveling wave solutions and analysis. Chinese J Phys. 2019;58:49–56.