Structure of Analytical and Numerical Wave Solutions for the Nonlinear (1 + 1)-Coupled Drinfel’d–Sokolov–Wilson System Arising in Shallow Water Waves

: In this article, we successfully obtain novel solutions for the coupled Drinfel’d–Sokolov–Wilson DSW system utilizing various methods. These include soliton solutions characterized by hyperbolic, rational, and trigonometric functions. Speciﬁcally, the generalized exponential rational function method (GERFM) and a modiﬁed version of the new Kudryashov method (MVNK) are employed to derive diverse soliton solutions for the system. Additionally, we demonstrate numerical solutions for the coupled Drinfel’d–Sokolov–Wilson system using adaptive moving mesh and uniform mesh methods. Also, we study the stability and error analysis of the numerical schemes. To validate the accuracy and reliability of the exact solutions obtained through analytical methods, we compare them with the numerical solutions both analytically and graphically. The techniques presented in this article are deemed suitable and acceptable and can be effectively applied to solve other nonlinear evolution systems.


Introduction
Partial differential equations (PDEs) are powerful tools for modeling and analyzing various phenomena in nature and physics.These include plasma physics, population dynamics, fluid and gas flow, electromagnetic fields, wave propagation in liquids, radiation, optical fibers, heat transfer, and other processes [1][2][3][4][5][6].PDEs offer an effective and successful approach to studying and comprehending these phenomena.Understanding soliton propagation is critical for grasping important oceanic phenomena, like how nonlinear waves move in shallow or deep seas, how transverse waves move in shallow water, how magneto-hydrodynamic waves move in plasma, and how phonon packets behave in nonlinear crystals.The study of PDEs allows us to better understand the future behavior of various phenomena accurately by utilizing exact and numerical solutions [7][8][9][10].It is essential to explore the precise solutions of nonlinear phenomena to gain insight into their long-term behavior.As a result, the development of systematic methods for deriving analytical solutions to PDEs has become a captivating and widely studied area of research among scholars.There are various techniques available to obtain analytical solutions to partial differential equations.Some of these techniques include the improved Kudryashov approach, the generalized direct algebraic strategy, the first integral approach, and many others.To learn more about these techniques and the analytical solutions they generate, one can refer to references [11][12][13][14][15][16].
In this study, we investigate a coupled system of nonlinear partial differential equations representing a nonlinear (1 + 1)-coupled Drinfel'd-Sokolov-Wilson system [17] Γ t + aΨΨ x = 0, Ψ t + bΨ xxx + cΓΨ x + dΨΓ x = 0, (1) where x and t are the space variable and time variable, both of which are, respectively, independent, and the dependent variables are Γ, Ψ. Finding solutions for the shallow water equations, also known as DSW equations, has been a significant area of research.These equations are versatile and can be used to model various water flow situations that involve gravity and shear stress.They are crucial for simulating significant events in engineering and physics, such as floods, dam breaks, and flows through vegetated regions.Therefore, much effort has been put into understanding and solving these equations, which are critical in various scientific and engineering applications.Various studies have analyzed the the coupled Drinfel'd-Sokolov-Wilson system for its numerical and exact solutions.For example, direct algebraic techniques, the Adomian decomposition method, the variational approach, and others, have been used to thoroughly study the the coupled Drinfel'd-Sokolov-Wilson system [18][19][20][21].Recently, there has been much interest in investigating the fractional-order coupled Drinfel'd-Sokolov-Wilson system using a variety of techniques, including the use of conformable derivatives, the homotopy perturbation method, the G G expansion, the F-expansion method, tanh and extended tanh methods, the truncated Painlevé method, and the exp-function method.Furthermore, the discrimination system for polynomial and Jacobi elliptical functions has been explored [22][23][24][25][26].In [27], the finite difference method was used to obtain numerical solutions for the coupled Drinfel'd-Sokolov-Wilson system.Previous studies on Equation (2) have predominantly concentrated on obtaining numerical solutions without actively addressing or reducing the resulting errors.In contrast, our study successfully reduced errors and received accurate numerical solutions for system (1).
This paper is motivated by the recent advancements discussed in the literature review.The primary objective of this study is to utilize the GERFM technique and the MVNK technique to develop various traveling wave solutions for the coupled Drinfel'd-Sokolov-Wilson system.The proposed methods have distinct advantages, which are listed below.These methods present a diverse range of reliable traveling wave solutions expressed in trigonometric, hyperbolic, and rational forms.These solutions act as dependable tools in interpreting intricate phenomena and gaining insights into their underlying dynamics.In addition, this paper aims to apply the adaptive moving mesh and uniform mesh method to system (1) in order to obtain its numerical solutions.It is worth noting that the initial condition for the numerical scheme is derived from the constructed exact solutions.The primary concept behind the employed numerical method is to distribute the mesh points in the solution curvature regions.While many researchers focus on analytically finding traveling wave solutions for the coupled Drinfel'd-Sokolov-Wilson system, only a few scientists have explored the numerical solutions to this problem with a high level of precision, aiming to minimize the error.The adaptive moving mesh approach ensures that the points are evenly distributed in areas with high error, effectively reducing the error.This technique, which is not commonly available in most numerical algorithms, significantly improves the accuracy of the results.To ensure the accuracy of the solutions, comparing the exact and numerical solutions is crucial.While some studies focus solely on finding exact solutions, this study goes a step further by comparing the exact and numerical solutions to ensure their accuracy and correctness.
The structure of this paper is as follows: Section 2 introduces the analytic solution for Equation (5).Section 3 focuses on presenting the numerical solution for system (1) using both fixed mesh and adaptive moving mesh techniques.This section also includes the results and subsequent discussion.Finally, Section 4 highlights the most significant findings discovered in this study.

Exact Solutions of the Coupled Drinfel'd-Sokolov-Wilson System
This section introduces the MVNK and the GERFM methods for finding soliton solutions of nonlinear evolution equations.
where Ψ = Ψ(x, t) is an unknown function, P 1 is a polynomial of Ψ(x, t), and various partial derivatives are involved in Equation (2).The wave transformation is applied as follows: into Equation (2) to reduce Equation (2) into an ordinary differential equation (ODE) expressed as after substituting the traveling wave into the system (1), we get the following expression: from Equation (3), we have by substituting Equation (4) into system (1), the system is converted into a single equation since ca 2s + ad s = α, Equation ( 5) reduces to an ODE given by −sΨ balancing Ψ ζζ with Ψ 3 in Equation ( 6) calculates the value of N = 1.

The MVNK Method
To apply MVNK [28],Equation (2) can be expressed as follows: , z j = 0 or c j = 0, (7) where the values of z 0 , z j , and c j for j = 1, 2, 3, . . ., N are to be later determined.The following function is satisfied by σ(ζ), and N is the homogeneous balance of numbers the Jacobi equation is satisfied by the following expression: the exact solutions of (5) are the constants z 0 , z 1 , and c 1 in Equation ( 9) are to be determined, such that z 1 = 0 or c 1 = 0, and σ(ζ) satisfies Equations ( 7) and (8).By substituting Equation ( 9) into Equation ( 6) and rearranging the terms, we obtain a system of algebraic equations.Solving this system gives us the following soliton solutions.Family 1: If we set Y = 0 in Equation ( 8), we have Cluster-1: , c 1 = 0, s = b, exact traveling wave solutions are provided by , If we set U 2 = 1 in Equation ( 8), then we have If we set U 2 = −1 in Equation ( 8), then we have , s = −2b, exact traveling wave solutions are provided by , If we set U 2 = 1 in Equation ( 8), then we have If we set U 2 = −1 in Equation ( 8), then we have Family 2: If we set Y = −1 2U in Equation ( 8), we obtain Cluster-1: , s = 4b, exact traveling wave solutions are provided by If we set U 2 = 1 2 in Equation ( 8), then we have If we set U 2 = − 1 2 in Equation ( 8), then we have

The GERFM Method
To apply GERFM [29,30], Equation (2) can be expressed as follows: the Υ(ζ) function satisfies the following differential equation: where v 1 , v 2 , v 5 , v 6 , k 1 , k 2 , k 3 , and k 4 are complex (or real) constants.According to the GERFM, with N = 1, the solutions of Equation ( 5) are the constants r 0 , r 1 , and a 1 in Equation ( 12) are to be determined, such that r 1 = 0 or a 1 = 0, and Υ(ζ) satisfies Equations ( 10) and (11).By substituting Equation ( 12) into Equation ( 6) and rearranging the terms, we obtain a system of algebraic equations.Solving this system gives us the following soliton solutions.11) can be expressed as The algebraic system is obtained by substituting Equation (13) into Equation ( 12) and then inserting the result into Equation ( 6).The solutions are as follows: , s = −2b, the exact traveling wave solutions are provided by , s = 4b, the exact traveling wave solutions are provided by , s = −8b, the exact traveling wave solutions are provided by Cluster-4: When r 0 = 0, , a 1 = 0, s = −2b, exact traveling wave solutions are provided by Family 2: The algebraic system is obtained by substituting Equation ( 15) into Equation ( 12) and then inserting the result into Equation ( 6).The solutions are as follows: , s = −4b, the exact traveling wave solutions are provided by , s = 2b, the exact traveling wave solutions are provided by , s = 2b, the exact traveling wave solutions are provided by , s = 8b, the exact traveling wave solutions are provided by

Numerical Results
This section uses two techniques (uniform mesh and adaptive moving mesh methods) to obtain the numerical solutions to the system (1).

Numerical Solutions Using a Uniform Mesh
In this subsection, we use a fixed mesh technique on a physical domain of size [0, L] to obtain numerical results for system (1).The domain is divided into N x subintervals of uniform width h = L/N x , denoted by [x j , x j+1 ], where x j = (j − 1)h for all x j ∈ [0, L] and j = 1, 2, . . ., N x + 1.To discretize system (1), we utilize the Crank-Nicolson method in the following manner: where the associated boundary conditions of system (1) are the initial condition is selected by calculating Equation ( 14) at t = 0 as follows: where a, b, c, d are constants.

Stability
In this section, we will assess the stability of the numerical solution by utilizing Fourier's stability technique.To apply this technique, we first need to linearize Equation ( 5) as follows: where b and D 0 are constants, such that D 0 = αΨ 2 .Fully discretize Equation ( 16) is given by

Let
Ψ n j = Y n e ιξ jh and then The following is the result of placing Equation ( 18) into the scheme (17): by dividing the two sides of Equation ( 19) by Ψ n j and taking Y as a common factor, we arrive at according to the analysis and as demonstrated in Equation ( 20), the numerical scheme remains stable without any conditions as long as the absolute value of Y is not greater than one.

Error Analysis
In this section, we use the Taylor series to look at the order of accuracy, which necessitates the initial discretization of Equation ( 5) Next, we determine the order by evaluating the truncation error.Assume that where e n+1 j represents the error, while Ψ n+1 j and Ψ(x j , t n+1 ) = Ψ represent the approximate and analytical solutions, respectively.Now, we substitute Equation (22) where T n j represents the truncation error and is written as consequently, the numerical scheme's truncation error is

Convergence
We need to refine two meshes with h and k values equaling zero to compute a sequence of computations using the initial data.Additionally, we must ensure that a convergent numerical scheme is used for every fixed point (x * , t * ) within a selected domain [a, b] and [0, T e ], x j → x * , t n → t * implies Ψ n j = Ψ(x * , t * ).As previously mentioned, the numerical scheme is unconditionally stable.Now, we will demonstrate that these schemes are also unconditionally convergent.To do so, let us assume that the error e is defined as e n j = Ψ n j − Ψ(x j , t n ).Now, Ψ n j exactly satisfies the scheme described by Equation (21) where T j is the truncation error; see Equation (24).If we suppose that the maximum error for the time step is given by then substituting Equation ( 25) into Equation ( 23) yields since the above inequality holds for each j = 1, 2, . . ., N x − 1, we have since the given initial data are used, we can identify E 0 = 0. Hence, the inequality is given by E n ≤ n × kT n j , hence, the scheme described by Equation ( 21) is convergent as h, k → 0.

Numerical Solutions Using an Adaptive Mesh
To obtain the numerical results of system (1), we employ the adaptive mesh method.The first step is to transform the physical domain [0, L] into the computational domain [0, 1] using the mapping x = x(η, t) : [0, 1] → [0, L], t > 0. This transformation allows us to apply the proposed numerical technique.By using the physical coordinate x and the computational coordinate η, we can represent the solution as Ψ = Ψ(x, t), Γ = Γ(x, t) where x = x(η, t).When the coordinate x is rearranged, it becomes x j (η) = x(η j , t) where η j = (j − 1)/N x , j = 1, 2, 3, . . ., N x + 1.The MMPDE that is most frequently used is as follows: assuming that τ and κ are constants, we introduce the function Q(Γ, Ψ, x), which is known as a monitor function.This function identifies regions in which the solution changes significantly, such as areas with large curvatures or high variations in the solution, and assigns more grid points to those regions.A commonly used monitor function, as proposed by Alharbi et al. [31], is given by The arc length : To apply the MMPDE5, we use the initial condition x j = (j − 1)L/N x and the boundary conditions x(0, t) = 0 and x(L, t) = L. Additionally, the computational coordinate is denoted by η j = (j − 1)h η , where h η = L/N x and j = 1, 2, . . ., N x + 1.To discretize the physical domain, we divide it into N x equal subintervals We then apply the chain rule to system (1) and obtain: with boundary conditions the initial condition is selected by calculating Equation ( 14) at t = 0 as follows: where a, b, c, d are constants.To clarify, the spatial derivative is semi-discretized while the temporal derivative remains continuous.This means that the problem described in system (1) becomes a set of ordinary differential equations, which can be solved using line methods.We use the MATLAB ODE solver (ode15i) to numerically integrate the resulting system.The discretization of Equation ( 28) must be achieved using finite differences in the following manner: where

Results and Discussion
Significant results have been achieved through the successful use of the GERFM and MVNK methods in extracting multiple exact solutions for Equation ( 5) in this study.It is important to mention that Equation ( 5) was obtained by converting the coupled Drinfel'd-Sokolov-Wilson system (1).Figures 1 and 2 display the exact solutions for Ψ 2 (x, t) and Γ 2 (x, t) for Equation (5).In Figure 1a,b, we used the MVNK method to create 3D plots of the analytical solution for Equation (5).Meanwhile, in Figure 2c,d, we also utilized the MVNK method to obtain 3D plots of the analytical solution for Equation (5).The parameters utilized are b = 1, a = 0.4, c = 0.5, d = 0.7, and N = 1000.Figure 3a,b demonstrate the behavior of waves when one parameter is changed while keeping others constant.The parameter being modified in this case is b, which affects both the direction and amplitude of the waves.It is apparent that when b is negative, the wave direction is always negative.
This article additionally investigates the numerical solutions of the coupled Drinfel'd-Sokolov-Wilson system (1) through the utilization of two numerical methodologies: the uniform mesh and the adaptive moving mesh methods.The results obtained are precise and efficient.The results in Figures 4-6 exhibit oscillations, but the adaptive moving mesh method effectively minimizes errors and produces acceptable results without oscillations.The study presents numerical results from time t = 0 → 10.As shown in Figure 4a,b, the adaptive moving mesh approach demonstrates more significant agreement with the analytical solution than the uniform approach.The insets in Figure 4a,b provide a closer look at the steep front regions of the plot.These regions are characterized by a rapid change in values and require a higher density of points for accurate representation.In the insets, these areas are indicated by yellow circles, indicating the need for more data points in those specific regions.As a result, the outcomes for both Ψ(x, t) and Γ(x, t) seem almost identical in these regions.Figure 4c depicts the obtained values of x = (η, t) using the MMPDE5 Equation ( 26) along with the monitor function Equation (27).The parameters employed for this calculation are as follows: b = 1, a = 1, c = 1, d = 1, x 0 = −10, with N x = 800, L = 20, and t = 10.These figures demonstrate that the monitor function sends additional points to areas with ∆x = 5 × 10 −5 curvatures while reducing the mesh points in other regions with ∆x = 5 × 10 −2 (where ∆x is the distance between points).The initial uniform mesh is set to ∆x = 1 × 10 −3 across all areas.Furthermore, it can be observed from Figure 4 that the numerical solutions obtained using a uniform mesh exhibit excellent agreement with the analytical solution.However, it is worth noting that achieving this level of accuracy requires a very small value of ∆x.Also, Figure 4 shows that numerical solutions achieved through a uniform mesh require a value of N x = 10,000 to display excellent agreement with the analytical solution.On the other hand, adaptive moving mesh only requires N x = 1000 to exhibit excellent agreement with the analytical solution.It is important to mention that the mesh is redistributed through the use of a monitor function that regulates its evaluation.The monitor function choice has a big impact on how the mesh moves.To prevent mesh tangles that often occur with mesh distribution, a special function called spatial smoothing is used.
Figures 5 and 6 show the results of an analysis that compared the solitary waves of an analytical solution Equation ( 5) and a numerical solution that uses an adaptive moving mesh.In Figure 5, 3D surface plots depict the analytical solutions for Ψ 13 and Γ 13 in (a) and (c), while (b) and (d) show the numerical solutions.The analysis used the following parameters: b = 1, a = 1, c = 1, d = 1, x 0 = −10, with N x = 10,000, and t = 0 → 10. Figure 6 shows 2D plots of the same solutions.The analytical solutions for Ψ 13 and Γ 13 are depicted in (a) and (c), while the numerical solutions are shown in (b) and (d).These figures offer sufficient evidence that the numerical and analytic solutions are highly similar.The clear implication is that the adaptive moving mesh method surpasses the uniform mesh approach in terms of accuracy, reliability, and convergence.This study's methods can be used to investigate other NPDEs.

Conclusions
In this study, we used the GERFM and MVNK methods to find some new traveling wave solutions for Equation ( 5) by expressing them in trigonometric and hyperbolic functions.We also employed the adaptive moving mesh method and uniform mesh to obtain numerical solutions while reducing the associated error.We observed that our methods yield superior outcomes compared to the cosine function approach.Previous studies on the coupled Drinfel'd-Sokolov-Wilson system mainly focused on generating numerical solutions without addressing error reduction.However, our study has successfully obtained numerical solutions while reducing the associated error.We found that the adaptive mesh approach is more reliable than the uniform mesh approach.Our results show that the adaptive mesh method effectively distributes points in regions with higher error.We also created 2D and 3D figures to illustrate the effectiveness of the employed approaches.Also, the investigation reveals that scheme (17) is second-order accurate in both k 2 and h 2 , and the scheme is unconditionally stable.Overall, the exact and numerical methods used in this study are flexible and effective in developing specific traveling wave solutions for nonlinear PDEs.
It is evident from the presented figures that the precise solutions exhibit similar characteristics to the numerical answers.Furthermore, the efficacy of the employed methodologies is clearly demonstrated in Table 1 and Figure 7.The table presents the L 2 error and CPU time required to achieve a time value of t = 10 using the adaptive approach.The initial number of points, N x , was set to 200.The L 2 error was observed to gradually converge towards a value of approximately 10 −5 , which is often considered acceptable in this context.The inaccuracy exhibits a significant decrease with the increase in the number of points.However, there is a slight increase in the CPU time as the number of points increases.For example, when the sample size is equal to N x = 1600, the L 2 error reaches a value of 1.85 × 10 −5 during a duration of t = 10 min.The marginal increase in CPU time can be attributed to the fact that the utilized functions are computed within the process.Moreover, Figure 7 illustrates a significant decrease in the relative error for Ψ and Γ as the number of points increases.The rate of error reduction is swift due to the adaptive process, which effectively provides the locations with high error with an ample and appropriate number of points.As a result, the adaptive technique demonstrates an increased processing capacity and applicability in the resolution of nonlinear partial differential equations.It has been observed that numerical systems with more mesh points are more precise than those with fewer mesh points but have more significant differences.By examining the crucial aspects of the solution's behavior, it was found that the non-uniform mesh can produce more accurate results than a fixed uniform mesh.Increasing the number of mesh points in the non-uniform mesh can improve its precision, but it requires additional time to solve partial differential equations.For example, when the number of points in the x direction is increased from 200 to approximately 1000, a non-uniform mesh is considerably faster than a fixed uniform mesh.This improvement in CPU time is more beneficial when the desired accuracy is low enough.Implicit non-uniform mesh schemes have shown promise in some problems due to their efficiency and ease of implementation, outperforming fixed uniform mesh schemes and h or hp-adaptive methods.But to keep the uniform mesh from becoming tangled up and disorganized, we need to look into different ways to make meshes in non-uniform frameworks.These techniques include the adaptive moving mesh method and optimal transport equations such as the Monge-Ampère and

Figure 7 .
Figure 7.This figures (a,b) show the convergence histories using L 2 norm as a function of the number of points N x .The parameters are given by a = 1, b = 1, c = 1, d = 1, x 0 = −10, L = 20, N = 800, and t = 10.
into Equation (21), , while Ψ(x j , t n ) omits the error indicated by the truncation error kT j

Table 1 .
The convergence histories using L 2 norm at t = 10.