Application of the Double Sumudu-Generalized Laplace Transform Decomposition Method to Solve Singular Pseudo-Hyperbolic Equations

: In this study, the technique established by the double Sumudu transform in combination with a new generalized Laplace transform decomposition method, which is called the double Sumudu-generalized Laplace transform decomposition method, is applied to solve general two-dimensional singular pseudo-hyperbolic equations subject to the initial conditions. The applicability of the proposed method is analyzed through demonstrative examples. The results obtained show that the procedure is easy to carry out and precise when employed for different linear and non-linear partial differential equations.


Introduction
Nonlinear partial differential equations are capable of explaining several phenomena, for instance, diffusion, electrodynamics, fluid dynamics, elasticity, etc.The majority of problems in physics and engineering can be represented by nonlinear differential equations.Pseudo-hyperbolic equations are a mathematical model that appear mostly in the theory of longitudinal and sideways vibrations of pseudo hyperbolic elastic bars (beams).The pseudo-hyperbolic equation is a partial differential containing mixed partial derivatives with respect to time and space.However, it is difficult to find clear solutions to linear and nonlinear pseudo-hyperbolic equations in general.There are several methods that have been discussed to solve some pseudo-hyperbolic equations, such as the residual power series method [1], the H1-Galerkin mixed finite element method in [2], the finite difference method [3], and the modified double Laplace transform method [4].A sufficient condition for the nonexistence of weak solutions to the nonlinear pseudo-hyperbolic equation has been offered in the Heisenberg theory [5].The definition of the parabolicity of systems of equations with partial derivatives, which generalizes the concept of Petrovsky parabolicity, has been developed in [6].Established on this explanation, a priori evaluations of pseudo hyperbolic PDEs with constant coefficients are conducted and the solvability of the Cauchy problem is formulated in [7].The modified Laplace Adomian decomposition is a famous method that was obtained by combining the Laplace transform and Adomian decomposition, see [8].This variety of equations was discussed in general non-fixed waves in a layered and spinning liquid, see [9].The approximate solution of multidimensional diffusion problems under an ABC fractional order derivative has been discussed by Nadeem et al. in [10].The authors of [11] developed a new method which is called a fractional analytical scheme in order to obtain the approximate results of fourth-order parabolic fractional partial differential equations.The authors of [12] have examined the solution of singular linear and nonlinear one-dimensional pseudo-hyperbolic equations through the double Laplace decomposition method.Lately, a strong and useful method has been proposed to find approximate analytical solutions to partial differential equations, which is named the double Laplace decomposition method [13].The Adomian decomposition method and series extension method have effectively been applied to obtain the solution for diffusion equations problems, see [14].In [15], the authors discussed the solution of the heat equation of fractional order in two dimensions by recommending the so-called triple Laplace transform method.The approximate analytical solution of the one-dimensional coupled nonlinear sine-Gordon equation was obtained by using the double Sumudu transform iterative method in [16].The authors of [17] studied the solution of nonlinear dynamical models with non-integer order derivatives by using a Laplace-type integral transform coupled with Adomian's approach.Recently the double Sumudu-generalized Laplace decomposition method was applied to obtain the solution of a 2 + 1-pseudo parabolic equation, see [18].The main aim of this work is to present a novel theorem for the double Sumudu-generalized Laplace transform.Moreover, we offer the application of this new transformation to partial derivatives.Finally, by examining examples, we discover that (DSGLTDM) is a strong tool for the solution of the linear, nonlinear, and coupled singular (2 + 1-D) pseudo-hyperbolic equations in comparison to the Adomian decomposition method (ADM), the homotopy analysis method (HAM), and the variational iteration method (VAM).This technique opens up new ways for solving complicated linear and nonlinear partial differential equations related to different scientific and engineering fields.
The rest of this work is arranged as follows.In Section 2, we start with some basic definitions and theorems of the existence condition of the (DSGLT) and the (DSGLT) of function xy f (x, y, t), respectively.In Section 3, we study two theorems that are considered the main contribution in this work.Section 4 presents how to obtain an approximate analytical solution for the 2 + 1-dimensional linear and nonlinear singular pseudo-hyperbolic equation using the (DSGLTDM), and for each problem we offer an example.In Section 5, we study the approximate analytical solution for the singular (2 + 1-dimensional) coupled pseudo-hyperbolic equation by using the (DSGLTDM).Finally, the concluding remarks are given in Section 6.

Some Fundamental Ideas of the (SGLT)
In this section, we introduce the essential definitions of the (DST), (GLT), (SGLT) and (DSGLT).

Definition 1 ([19]
).The double Sumudu transform of the function f (x, t) in the positive quadrant of the xt-plane is given by The (GLT) of the function f (t) is denoted by G α in the following definition.
Definition 2 ([20]).If f (t) is an integrable function defined for all t ≥ 0, its (GLT) G α is the integral of f (t) times s α e − t s from t = 0 to ∞.It is a function of s, say F(s), and is denoted by where s ∈ C and α ∈ Z.
The (SGLT) of the function f (x, t) is determined by F(u 1 , s) in the following definition.
Definition 3 ([21]).The (SGLT) of the function f (x, t) of two variables x and t defined in the first quadrant of the xt-plane by the double integral in the form of S x G t indicates the (SGLT) and the symbols u 1 and s denote transforms of the variables x and t in Sumudu and generalized Laplace transforms, respectively.
The inverse (SGLT) ) is determined by the following formula: The definition of the (DSGLT) is given in the following.

Definition 4 ([18]
).Let f (x, y, t) be a continuous function of three variables x, y and t.Then, the (DSGLT) is defined by where S x S y G t denotes the (DSGLT) and u 1, u 2 and s represent transforms of the variables x, y and t in the (DSGLT), respectively.
The existence condition of the (DSGLT) is given in the next theorem.
Theorem 1 ([18]).The function f (x, y, t) is defined on (0, X), (0, Y) and (0, T) of exponential order (x, y, t), then the double Sumudu-generalized Laplace transform of f (x, y, t) exists for all Re 1 The following theorem is useful in this study.
Theorem 2 ([18]).If the (DSGLT) of the function f (x, y, t) is presented by S x S y G t ( f (x, y, t)) = F(u 1 , u 2 , s), then the (DSGLT) of the functions xy f (x, y, t), determined by In the following, we address the (DSGLT) of the function ψ x , ψ xx , ψ tt , which is defined in [18] as follows

The Idea of the (DSGLT) to the Functions yψ tt , xψ tt and xyψ tt
In this part, we study two theorems which are the main contributions of this paper; in the first theorem, we discuss how to apply the (DSGLT) to the functions yψ tt and xψ tt , and in the second theorem, we apply the (DSGLT) to the function xyψ tt .Theorem 3. The (DSGLTs) of the partial derivatives yψ tt and xψ tt are denoted by Proof.On utilizing the partial derivatives with respect to u 2 for Equation (2), one can obtain calculating the derivative in the above equation according to u 2 , we yield By placing Equation (11) into Equation ( 10), we get Equation ( 12) becomes consequently By organizing Equation ( 14), we will obtain the proof for Equation (8) as follows In the same way, we can prove Equation (8).
Theorem 4. The (DSGLT) of the partial derivatives xyψ tt is specified by Proof.By taking the partial derivative according to u 1 for Equation (2), we have we calculate the partial derivative inside the brackets as follows: Placing Equation ( 17) into Equation ( 16), we yield the partial derivative with respect to u 2 for Equation ( 18) is computed as the following Therefore, we rewrite Equation ( 19) as follows One can rearrange Equation ( 21) and apply Equation (7) to prove Equation ( 15)
Allow us to consider a general singular (2 + 1-dimensional) pseudo-hyperbolic equation in the following form to the following initial conditions: where 1 x (x ψ x ) x and 1 y y ψ y y are known as Bessel's operators and f (x, y, t), f 1 (x, y) and f 2 (x, y) are familiar functions.The solution of Equation ( 22), is given by the following steps: Step 1: Multiplying Equation ( 22) by xy, we obtain Step 2: By employing the (DSGLT) first for both sides of Equation ( 24), applying the differentiation property of the (DST) and the initial condition stated in Equation ( 23) and employing Theorems 1 and 3, we gain Step 3: By multiplying Equation (25) by s 2 u 1 u 2 , using integral, from 0 to u 1 and 0 to u 2 with respect to u 1 and u 2 , respectively, and dividing the new equation by u 1 u 2 , we have Step 4: The series solution of the singular (2 + 1-dimensional) pseudo-hyperbolic equation denoted by: where Ψ(u 1 , u 2 , s), F 1 (u 1 , u 2 ) and F 2 (u 1 , u 2 ) are the (DSGLT) and (DST) of the functions ψ(x, y, t), f 1 (x, y) and f 2 (x, y), respectively.
Step 5: On using the inverse (DSGLT) for both sides of Equation ( 26), and with the assistance of Equation ( 27), we find In light of the first approximation, we have the rest of the components, ψ n+1 , n ≥ 0, are described by Assume that the inverse of the (DSGLT) with respect to u 1 , u 2 , and s for Equations ( 28) and (29) exists.To offer an obvious outline of this method, we provide the following example.
Example 1.The singular (2 + 1-D) pseudo-hyperbolic equation is given by: subject to the initial ψ(x, y, 0) = 0, ψ t (x, y, 0) = x 2 − y 2 (31) By performing the above-mentioned method, the elements of the sequence are computed as follows: (32) As stated by the (DSGLTDM), we gain the following elements at n = 0: Likewise, at n = 1, we obtain that In the same way, let n = 2 As the sum of all terms, we have hence, the approximation solution of Equation ( 30) is confirmed by We obtain the exact solution: In the forthcoming problem, we apply the proposed method to a nonlinear pseudohyperbolic equation.
with initial conditions By performing the previous analysis, the primary approximation is described by the remaining terms are identified by where nonlinear terms A n and B n are decomposed as and the nonlinear terms ψ ∂ψ ∂x 1 and ψ ∂ψ ∂x 2 are in the following forms: and To show the applicability of this method for nonlinear problems, we study the following example.

Example 2. Consider the following nonlinear singular pseudo-hyperbolic equation
with the following conditions By utilizing the stated method above, we have where A n and B n are defined in Equations ( 39) and (40).The subsequent terms are denoted by By following the same method, we have Hence, as stated in Equation ( 27), we get and the exact solution of Equation ( 41) is demonstrated by

Singular (2 + 1-D) Coupled Pseudo-Hyperbolic Equation and the (DSGLTDM)
The fundamental evaluation of the (DSGLTDM) is clearly shown to indicate its ability and high accuracy by examining the general coupled singular (2 + 1-D) pseudo-hyperbolic equation of the form with the following initial conditions where f (x, y, t), g(x, y, t), f 1 (x, y), f 2 (x, y), g 1 (x, y) and g 1 (x, y) are given functions.By us- ing the (DSGLTDM), this method contains the following steps.
Step 1: Multiply both sides of Equation (45) by xy, resulting in the following equation, and Step 2: Using the (DSGLT) for both sides of Equations ( 47) and (48) and the (DST) for Equation (46), we obtain and Step 3: By operating the integral for both sides for Equations ( 50) and ( 49), from 0 to u 1 and 0 to u 2 with respect to u 1 and u 2 , respectively, we have and The series solution of Equation ( 45) is thus entirely stated by Step 4: By applying the inverse (DSGLT) for both sides of Equations ( 51) and ( 52) and applying Equation (53), one can obtain and From the above equations, we equate terms on both sides to produce recurrence relations as with initial condition ϕ(x, y, 0) = x 2 − y 2 , ϕ t (x 1 , x 2 , 0) = − x 2 − y 2 ω(x, y, 0) = x 2 − y 2 , ω t (x, y, 0) = − x 2 − y 2 , (58) By using our method above, we acquire the following components: and the remaining terms ϕ n+1 and ω n+1 , n ≥ 0, are presented by and hence the exact solution becomes ϕ(x, y, t) = x 2 − y 2 e −t , ω(x, y, t) = x 2 − y 2 e −t .

Conclusions
The (DSGLTDM) is a hybrid method that consists of the (DSGLT) and the (ADM).This method has some advantages over the previously studied methods.In this research work, we have utilized the (DSGLTDM) in order to obtain the approximate and series solutions of linear, nonlinear, and coupled singular (2 + 1-D) pseudo-hyperbolic equations.
By investigating examples, we found that the (DSGLTDM) is a powerful tool for the solution of linear, nonlinear, and coupled singular (2 + 1-D) pseudo-hyperbolic equations and compared it with the Adomian decomposition method (ADM), the homotopy analysis method (HAM), and the variational iteration method (VAM).Nevertheless, until now, there has been an open problem to investigate the rate of convergence to the exact solution for these types of problems.It can be studied by the (DSGLTDM) using an analytical solution to the other singular 2 + 1-dimensional partial differential equations, which occur in science as well as engineering, and may deliver a better knowledge of the real-world problems that are described by singular 2 + 1-dimensional partial differential equations.
Funding: The author would like to extend their sincere appreciation to Researchers Supporting Project number (RSPD 2023R948), King Saud University, Riyadh, Saudi Arabia.