Abstract

Through this article, we focus on the extension of travelling wave solutions for a prevalent nonlinear pseudo-parabolic physical Oskolkov model for Kevin-Voigt fluids by using two integral techniques. First of all, we explore the bifurcation and phase portraits of the model for different parametric conditions via a dynamical system approach. We derive smooth waves of the bright bell and dark bell, periodic waves, and singular waves of dark and bright cusps, in correspondence to homoclinic, periodic, and open orbits with cusp, respectively. Each orbit of the phase portraits is envisaged through various energy states. Secondly, with the help of a prevalent unified scheme, an inventive version of exact analytic solutions comprising hyperbolic, trigonometric, and rational functions can be invented with some collective parameters. The unified scheme is an excitably auspicious method to procure novel interacting travelling wave solutions and to obtain multipeaked bright and dark solitons, shock waves, bright bell waves with single and double shocks, combo waves of the bright-dark bell and dark-bright bell with a shock, dark bell into a double shock wave, and bright-dark multirogue type wave solutions of the model. The dynamics of the procured nonlinear wave solutions are also presented through 2-D, 3-D, and density plots with specified parameters.

1. Introduction

Recently, discoveries of novel explicit and analytical travelling wave solutions of nonlinear models have attracted great attention from blooming scientists [13]. A series of travelling wave exact solutions with various nonlinear dynamics can be illustrated as unique complex phenomena in diverse scopes of nonlinear sciences, for instance, optics [13], fluid mechanics [4, 5], shallow water wave and magneto-sound transmission in plasmas [6], ion-acoustic plasma waves [7, 8], shallow water waves [9, 10], solid-state physics and chemical kinetics [11], visco-elastic Kelvin Voigt fluid [12], shallow water propagation in the harbor and coastal region [13], and nonlinear quantum field theory [14] . Due to the significance of such field of studies, theoretical probes on physical action emulation of NEEs are sought day by day. Hence, much effort has been devoted to securing exact and analytical travelling wave solutions of the NEEs. Many scholars are compelled by NEEs to raise travelling wave solutions by implementing a number of integral approaches. These approaches that are accomplished in current literature include the extended Kudryashov method [1], the first-integral method [2], the extended trial equation method [3], Bäcklund transformation [4, 5], the extended simple equation method [6], modified simple equation method [7, 8], Hirota bilinear form [9, 10], new generalized -expansion method [11], -expansion method [14], the extended -expansion method [15], Cole-Hopf transformation method [16], Adomian decomposition method [17], exp(−φ(ξ))-expansion method [1820], planner dynamical system scheme [21, 22] and many more.

Among the above procedures, the planner dynamical scheme is the most important, stable, and reliable approach as it provides various types of phase orbits. Depending on the parametric conditions, each type of orbit (periodic, heteroclinic, and homoclinic) yields periodic, kink, and bell wave solutions, respectively, under selective initial/boundary conditions. This type of approach is still least explored on numerous nonlinear models; hence, it breeds a vast challenge for proximate ascertainment of the models’ dynamical behaviors. The (1 + 1)-dimensional Oskolkov model for Kevin-Voigt fluids is one of the nonlinear models [2325] with huge applications in contemplating incompressible visco-elastic Kelvin-Voigt fluid flows. This category of the pseudo-parabolic model, which involves time derivative in the highest differential term, can also be exemplified for wider purviews such as fluid flows of fissured rock, second order shear of fluids, thermodynamics, consolidation of dust, and transmission of long wave propagations with small amplitudes. Hence, the Oskolkov model is anticipatorily examined via diverse schemes by various researchers across the globe [2325]. Several fundamental and featured procedures have also been established to procure the exact solutions of the Oskolkov model.

The modified simple equation method [12] was used to derive more exact and explicit solutions of the Oskolkov model. The tanh-coth method [23] was applied to accomplish exact solutions of the nonlinear pseudo-parabolic model. Asymptotic stable and unstable solutions to the Oskolkov model were analyzed graphically too [24]. The nonlinear shock wave, bell wave, and various dynamical motions in the presence of an exotic periodic force were discussed briefly for the generalized Oskolkov model [25]. Although the enormous investigation and profound dynamical solutions to the model have been scrutinized by many researchers, there are still some unresolved problems like bounded travelling wave solutions with their corresponding phase orbit, solutions with boundary conditions, and rogue wave solutions. All these updates spark inquisitions on how the results are exhibited and how they vary with the changing parameters. To retort such accounts, we need to derive solutions involving some of the prominent parameters with suitable initial/boundary conditions.

Through this article, we shed light on the progress of nonlinear travelling waves such as bright and dark bell waves, periodic waves, and singular bright and dark cusp waves as corresponding to each orbit of the phase portraits on the dynamical Oskolkov model by utilizing the planner dynamical strategy [20, 21]. Besides, rogue waves, multisoliton, interaction of solitons, and breather wave solutions, which are no less fascinating, have also received gradually increased attention in the field of nonlinear dynamical systems [2628]. The long wave limit of the breather localized wave was derived by Yue et al. [29]. Abundant research on these topics has been performed and different types of nonlinear localized waves have been revisited for strategic recipes of various nonlinear physical models materialized in mathematical physics and engineering applications [30, 31]. The Haar wavelet technique has been used to analyze the dynamical phenomena of the fractional-order nonlinear model [32, 33]. Nowadays, investigations on the fractional model are studied by analytical techniques [34, 35] and numerically [36, 37], exposing huge efforts. Notwithstanding, we also aim to construct the solutions to multirogue type waves, shock waves, interactions of bright-dark and dark-bright bell solitons, interactions of shock-bell waves, and singular multisoliton to the (1 + 1)-dimensional Oskolkov model by taking benefit of the lavish unified method [38, 39].

The creation of this article is designed as follows: in Section 2, we deliberate the bifurcation and phase portrait analyses of the Oskolkov model. In Section 3, we exhibit parametric expressions of explicit solutions of the Oskolkov model. The exact wave solutions and their graphical illustrations of the model via the unified method are presented in Section 4. Finally, the perspectives of the current study on this kind of dynamical wave phenomenon are discussed in the last section.

2. Bifurcations and Phase Portraits of the Oskolkov Model

The regular form of the (1 + 1) dimensional Oskolkov model (1) iswhere and in which is a constant whereas is the wave speed. Now, we reinstate (1) into the form of (ODE) as follows:where the primes stand for the derivative with respect to

Equation (2) in a concise form can be reformulated aswhere , and

This (3) can further be written into a prototype system of differential formwhich yields a Hamiltonian

Next, we attempt to detect the bifurcation of the phase orbits of the dynamical system in (4) with various conditions on the parameters From the observation, it is clear that a smooth homoclinic orbit of the system in (4) arises from the smooth solitary waves in equation (1). When the solution of of the system guarantees and , is said to be homoclinic and heteroclinic orbits for and respectively. Normally, a homoclinic orbit of the (4) correlates with a solitary wave solution of the (1) and a heteroclinic orbit of the (4) correlates with a king (antiking) wave solution of the (1). Alike, a periodic orbit of the system in (4) keeps up a correspondence with a periodic travelling wave solution of (1). Therefore, we can explore to locate all possible phase portraits of the dynamical system in (4) with the help of different constraints on the parameters

For critical points in an equilibrium situation, we have to consider and then that prototype (4) will provide two equilibrium points and , if . Furthermore, the dynamical prototype (4) yields one equilibrium point at . Now, let us assume, to be the coefficients’ matrix of the linearized prototype (4) at equilibrium points and let . Therefore, we have

With the help of the bifurcation theorem [21, 22] and the above inspection, we will get the following remarks:Group-1: for the nature at the origin is a saddle point and is a center point, whereas the corresponding bifurcations of phase portraits of the prototype (4) are visualized in Figures 1(a) and 1(b), respectively.Group-2: for the nature at origin is a center point and is a saddle point, whereas the corresponding bifurcations of phase portraits of the prototype (4) are depicted in Figures 2(a) and 2(b), respectively.Group-3: for the nature at origin is a cusp point, whereas the corresponding bifurcations of phase portraits of the prototype (4) are depicted in Figures 3(a) and 3(b), respectively.

3. Explicit Expressions of Solution for the Model (1)

In this segment, we derive various explicit parametric illustrations of travelling wave solutions for the model (1). For effortlessness, the energy level of the Hamiltonian is specified as and .

3.1. Examining Group-1 Based on Figure 1 Observation (Part 2)
(I)For the state the allied homoclinic orbit at is described via where the model (1) suggests itself as a valley type smooth solitary wave solution specified by Figure 1(a). Using the relation in equation (5), yieldsCombining the first equation of the prototype (4) and (7), we attain the parametric expression of the valley-type smooth solitary wave solution as follows:where , The profile of the valley solitary wave solution (8) is shown in Figure 4(a) with (II)For the state it is equivalent to the homoclinic orbit to the stable point at identified by where the model (1) suggests itself as a valley-type smooth solitary wave solution portrayed in Figure 1(b). Using the relation in equation (5), yieldsCombining the first equation of the prototype (5) and (9), we attain the parametric expression of the valley-type smooth solitary wave solution as follows:where , The profile of the valley solitary wave solution (10) is shown in Figure 4(b) with (III)For the state or , the model (1) represents a set of smooth periodic wave solutions described by as can be viewed in Figures 1(a) and 1(b), respectively. For the state , the phase portrait of the prototype (4) is shown in Figure 1(a). In this phase, the formulation of the closed domain can be written as: , (11) where and are the meeting points of the curve identified by the orbits on the -axis that preserve the condition .Combining the first equation of the prototype equations (4) and (11), we attain the parametric formulation of the periodic solution as follows:

The shape of the periodic outcomes of (12) is displayed in Figure 4(c) with

For the state similar investigation can be used in Figure 1(b). Suppose, and are the meeting points of the curve identified by the orbits on the -axis that preserve the condition .

We acquire the parametric formulation of the periodic solution as follows:

The shape of the periodic outcomes of (13) is presented in Figure 4(d) with

3.2. Examining Group-2 Based on Figure 2 Observation (Part 2)
(I)For the state it is equivalent to the homoclinic orbit at the stable point identified by where the prototype (4) suggests itself as a smooth solitary wave solution of valley type as depicted in Figure 2(a). With the help of the relation yieldsThen, by means of the first equation of the prototype (4) and (14), we attain the parametric expression of the valley-type smooth solitary wave solution that has a similar nature as depicted in Figure 4(a) as follows:where , (II)For the state it is related to the homoclinic orbit at the stable point identified by where the model (1) represents a valley-type smooth solitary wave solution as displayed in Figure 2(b). With the help of the relation , yieldsThen, by means of the first equation of the prototype (4) and (16), we attain the parametric expression of the valley-type smooth solitary wave solution that has similar nature as depicted in Figure 4(b) as follows:where , (III)For the state or , the prototype (4) has a group of smooth periodic wave solutions identified by , which can be observed in Figures 2(a) and 2(b), respectively. For this aspect, the formulation of the periodic solution is identical to the solutions in (12) and (13).
3.3. Examining Group-3 Based on Figure 3 Observation (Part 2)
(I)For the state there is an unrestricted orbit with the alike Hamiltonian as the origin (Figure 3(a)). The open cusp orbit can be specified byThen, by means of the first equation of the prototype (4) and (18), we acquire the periodic cusp wave solution as follows:where , The shape of the singular bright cusp wave of (19) is depicted in Figure 5(a).(II)For the state there is also an unrestricted orbit with the alike Hamiltonian as the origin (Figure 3(b)). Due to this condition, we also catch out the similar periodic cusp wave solution in the form as follows:where ,

The profile of the singular dark cusp wave solution in (20) is depicted in Figure 5(b).

4. Solutions of the Oskolkov Model via the Unified Method and Graphical Illustration

Let us apply a steadfast treatment of the unified method [32, 33] to the (1 + 1) dimensional Oskolkov model (1) in this portion. Graphical presentations are also outlined here.

4.1. Solutions of the Oskolkov Model

Using the travelling wave variable, we reinstate (1) into the ODE (2). Next, we enumerate the balance number of expressions in (2) between the linear term, and the nonlinear term, by letting Through this effort, a trial solution can be formed as follows:

Here are unfamiliar constants to be formed subsequently and we consider the derivative of satisfying an ODE namely the Riccati differential equation

The solutions of the given Riccati equation are given as follows:Case-01: hyperbolic function (when ) as follows:Case-02: trigonometric function (when ) as follows:where and , are real arbitrary constants.Case-03: rational function solutions (when )

Differentiating (21) as many times as necessary while satisfying (22), these expressions are substituted back into (1). Next, we attain a polynomial of where and then by letting the coefficients of equal to zero yields as follows:

By solving the above system of equations, many sets of solutions can be attained as follows:

Set-01:

Set-02:

Set-03:

Set-04:

Set-05:

Set-06:

Now, by plugging these expressions into (21) for unfamiliar constants, different types of solutions for (1) can be acquired.

For the Set-01, if , then we obtain the hyperbolic function solutionswhere while , are arbitrary constants.

For the Set-01, if , then we obtain trigonometric function solutionswhere while and are real arbitrary constants.

If , we can search out rational function type solutions, but in this case, if we consider then the executed values of and will be undefined. So, the solution is rejected. Similarly, the other sets (Set-02 to Set-06) each provide us with eight solutions in the same way.

4.2. Numerical Illustrations of the Obtained Results and Discussion

In this section, we will deliberate the nature and behavior of all the solutions physically with 2-D and 3-D graphs, as well as density plots of the (1 + 1) dimensional Oskolkov model by instigating the most recent method, namely, the unified method. The solution conveys the multipeaked bright solitonic nature for the parametric values . Oppositely, the solution conveys the multipeaked dark solitonic nature for the parametric values as illustrated in Figure 6. The real part of presents the kink type shock wave whereas its imaginary part presents the combo bright-dark bell wave for the parameters as illustrated in Figure 7. On the other hand, the real part of conveys the combo dark-bright bell wave with a shock wave, whereas its imaginary part expresses the bright bell wave with a shock for the parameter , as illustrated in Figure 8. The real part of the solution presents the dark bell as a double shock wave but its imaginary part presents the double bright-dark bell wave with a shock for the parameters , as demonstrated in Figure 9. The real part of conveys the bright bell wave with double shocks but its imaginary part conveys the bright-dark bell wave with a shock for the parametric values , as demonstrated in Figure 10. Both the real and imaginary parts of represents the multi-rogue type waves with singularities for , as demonstrated in Figure 11. Yet, every solution obtained for the model by using this technique can be reduced to the multirogue (bright-dark) type wave with the changing signs of the parameters and . Such a dynamic is displayed in Figure 12 of the solution for the parametric values . Besides, we display its density and 2D plots of corresponding real and imaginary parts altogether.

5. Conclusions

Through this work, the key effort is to procure, test, and explore the innovative travelling wave solutions and the physical properties of the nonlinear Oskolkov model for Kevin-Voigt fluids by utilizing two reliable mathematical techniques. Based on the bifurcation theory of the dynamical scheme, we explored the bifurcation and derived a phase portrait of the model using diverse parametric conditions. We have derived smooth bright bell and dark bell waves, periodic solutions as well as singular dark cusp and bright cusp waves. These travelling waves are achieved by corresponding to homoclinic, periodic, and open orbits with the cusp of each orbit of the phase portraits, respectively. The unified method as a substantial way to derive the new travelling wave solutions has been utilized through which we have successfully achieved the solutions for multipeaked bright and dark solitons, shock waves, bright bell waves with a shock, bright bell waves with a double shock, combo bright-dark and dark-bright bell waves with a shock, dark bell into a double shock wave, and bright-dark multirogue type waves. The bifurcation analysis of this model is presented first in this article, even though most of the obtained solutions derived by direct integration from the Hamiltonian of the model according to each energy is the first step here. In particular, the bright bell waves with a shock, bright bell waves with double shock, combo bright-dark, combo dark-bright bell waves with a shock, dark bell into a double shock wave, and bright-dark multirogue type waves were not reported in the previous literature. This technique proposes solutions with free parameters that could be significant to deliberate more intricate nonlinear corporal phenomena. The obtained solutions in this paper reaffirmed that the method is very effective and easily applicable to formulate more exact travelling wave solutions than other methods of the nonlinear evolution models occurring in many mathematical physics and engineering applications.

Data Availability

The data used to support the findings of this study are included in the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.