Semi-analytical solution of non-homogeneous Duffing oscillator equation by the Padé differential transformation algorithm

Abstract

The Duffing oscillator equation is one of important equations that model several nonlinear phenomena in science and engineering. The differential transform method (DTM) is applied to obtain the solutions of homogeneous and non-homogeneous Duffing oscillator equations under the influence of different initial conditions. In fact, the DTM can only find the solutions of the Duffing oscillator equation in a small domain, therefore the method has been improved using Padé approximation. We obtain the solutions by the Padé-DTM (PDTM) in long domain. This article proves the validity, simplicity and applicability of the PDTM method. In addition, the accuracy of the PDTM and the instability of the solution are discussed.

Keywords

Differential transformation method the Padé approximation forced Duffing oscillator equation

Introduction

Over the recent decades, many scientific phenomena and applications have been described by nonlinear ordinary differential equations (ODEs) and Partial differential equations (PDEs).^1–6 One of these significant equations is called the Duffing oscillator equation which is utilized in many sciences such as physics and engineering.^7,8 It was discovered in 1918 by electrical engineer Ger man Duffing.⁹ The general form of Duffing equation¹⁰ is given by the following form

\ddot{q} (t) + α \dot{q} (t) + β q (t) + γ q^{3} (t) = f (t),

(1)

and it subjects to following initial conditions (ICs)

q (0) = w, \dot{q} (0) = z,

(2)

where α, β, γ, w and z are real constants. The parameters α, β and γ control the size of damping, stiffness and amount of non-linearity in the restoring force, respectively. If β > 0, the system is hardening, dealing with a softening system if β < 0 and if β = 0, the system is softening stiffness system. The Duffing equation has many applications such as pendulum application. In nonlinear dynamics, the pendulum is the quintessential dynamical system. The equation of motion describing the angular displacement θ of the harmonically excited pendulum is

m l^{2} \frac{d^{2} θ}{d t^{2}} + m g l \sin θ = M \cos ω t,

(3)

where l is the length of the pendulum, g is the acceleration due to gravity and M is the amplitude of the applied moment.¹¹ Because of the negative cubic term in the stiffness moment, the pendulum exhibits a softening stiffness characteristic.¹²

The Duffing equation has been solved in the literature by several methods such as Variation iteration method,^13–15 Homotopy perturbation transform method,^16,17 an effective approach,¹⁸ a new approach method,¹⁹ Homotopy analysis method^9,20 and modified differential transform method.^2,21 The mathematician continuously have improved the computational methods to obtain fast and high accurate results. In this paper, we will turn the light on the differential transformation method (DTM).

The DTM was first introduced by Pukhov to compute linear and nonlinear initial value problems in electric circuit analyses.²² Chen and Ho invented the DTM for solving PDE and discovered closed-form series solutions for various linear and nonlinear initial value problems.²³ Hassan showed that the DTM is usable on a large selection of PDEs and is able to effortlessly get closed-form solutions.^24–26 The results by the DTM is more accurate than some other numerical and semi-analytical methods for solving ODEs.²⁷ So the advantages of the two-dimensional DTM are ability to find accurate solutions and unaffected by tiny or big numbers unlike perturbation approaches. In addition, the DTM can be used regardless of whether the governing equations and boundary/initial conditions of a given nonlinear problem have modest or big numbers and the DTM solves equations without the use of starting estimations or an auxiliary linear operator. On the other hand, one of the DTM disadvantages is solving the solution in small domain for some problems in physics.^28–30 Therefore, the DTM needs to be improved. Recently, the Padé approximation has been used to improve some iterative methods. The Padé approximation was introduced by Her-Moths and Padé. It has been utilized in physics for over 40 years and there are numerous notable examples in the field. Therefore, we aim to use the Padé approximation technique in order to improve the DTM to solve the Duffing oscillator equation.³¹ However, the DTM does not require a small parameter such as the traditional perturbation methods or a small embedded parameter such as the Homotopy perturbation method (HPM). In addition, the HPM fails to solve the damping nonlinear oscillator.³² In fact, the DTM does not required integration such as Adomian decomposition method (ADM) or Laplacian Adomian decomposition method (LADM). The ADM and LADM fail to solve non-integrable equation because their numerical scheme includes integration operators.

The novelty of this article is presenting the semi-analytical solution of the non-homogeneous Duffing oscillator equation by the Padé differential transformation method (PDTM) and finding the relation between the number of iteration and the order of the Padé approximation.

This article is organized as follows: Padé-Differential Transform Method presents steps of the algorithm of PDTM. Stability introduces some applications for forced Duffing oscillator equation and shows the solutions by using PDTM. Applications is the study of the PDTM accuracy and last section is the conclusion of the work.

Padé-Differential transform method

This section explains the steps of a logarithm of the PDTM. In general, the idea starts by finding a solution as a power series by the DTM. However the radius of convergence of this power series may not be big enough to hold the two limits. Therefore, the Padé approximation is used to modify the produced series as rational functions since the Padé approximation delivers the best estimate by converting a function’s approximate polynomial into rational functions of polynomials of a specific order.^33–35

Assume the power series $\sum_{j = 0}^{\infty} a_{j} y^{j}$ representing a function f(y)

f (y) = \sum_{j = 0}^{\infty} a_{j} y^{j} .

(4)

The Padé approximate is a rational function that has the form

[P / N] = \frac{G_{P} (y)}{O_{N} (y)},

(5)

where G_P(y) and O_N(y) are a polynomial of degree P and N at most, respectively. Thus, we have

\begin{array}{l} f (y) = & b_{0} + b_{1} y + b_{2} y^{2} + b_{3} y^{3} + \dots .., \\ G_{P} (y) = & g_{0} + g_{1} y + g_{2} y^{2} + g_{3} y^{3} + \dots + g_{P} y^{P}, \\ O_{N} (y) = & o_{0} + o_{1} x + o_{2} y^{2} + o_{3} y^{3} + \dots + o_{N} y^{N} . \end{array}

The numerator and denominator coefficients are P + 1 and N + 1, respectively, in equation (5) as can be seen. We can multiply the numerator and denominator by a constant while keeping [P/N] intact and also apply the normalization constraint O_N(0) = 1. Because of the P + 1 independent numerator and N independent denominator coefficients, there are a total of P + N + 1 unknown coefficients. This value suggests that the power series should generally fit [P/N] equation (4) in the orders 1, y, y², y³, ….., y^1+N. As a result, we have a formal power series

\sum_{j = 0}^{\infty} b_{j} y^{j} = \frac{g_{0} + g_{1} x + g_{2} y^{2} + \dots + g_{P} y^{L}}{o_{0} + o_{1} y + o_{2} y^{2} + \dots + o_{N} y^{N}} + O (y^{P + N + 1}) .

We achieve that by cross multiplying

(b_{0} + b_{1} y + b_{2} y^{2} + b_{3} y^{3} + \dots) (o_{0} + o_{1} y + o_{2} y^{2} + \dots + o_{N} y^{N}) = g_{0} + g_{1} y + g_{2} y^{2} + \dots + g_{P} y^{P} + D (x^{P + N + 1}) .

We get the following sets of equations by combining the coefficients of 1, y, y², y³, …, y^P and the coefficients of y^P+1, y^P+2, …, y^P+N

{\begin{array}{l} b_{0} = g_{0}, \\ b_{1} + b_{0} o_{1} = g_{1}, \\ b_{2} + b_{1} o_{1} + b_{0} o_{2} = g_{2}, \\ ⋮ \\ b_{P} + b_{P - 1} o_{1} + b_{P - 2} o_{2} + \dots + b_{0} o_{P} = g_{P} . \end{array}

(6)

and

{\begin{array}{l} b_{P + 1} + b_{P} o_{1} + \dots + b_{P - N + 1} o_{N} = 0, \\ b_{P + 2} + b_{P + 1} o_{1} + \dots + b_{P - N + 2} o_{N} = 0, \\ ⋮ \\ b_{P + N} + b_{P + N - 1} o_{1} + \dots + b_{P} o_{N} = 0. \end{array}

(7)

If m < 0 is true, then b_m = 0 is true for consistency, and o_i = 0 is true for accuracy, i > N.

By solving equations (6) and (7), we obtain

[P / N] = \frac{| \begin{matrix} b_{P - N + 1} & b_{P - N + 2} & \dots & b_{P + N} \\ ⋮ & ⋮ & ⋱ & ⋮ \\ b_{P} & b_{P + 1} & \dots & b_{P + N} \\ \sum_{i = N}^{P} b_{i - N} y^{i} & \sum_{i = N - 1}^{P} b_{i - N + 1} y^{i} & \dots & \sum_{i = 0}^{L} b_{i} y^{i} \end{matrix} |}{| \begin{matrix} b_{P - N + 1} & b_{P - N + 2} & \dots & b_{P + 1} \\ ⋮ & ⋮ & ⋱ & ⋮ \\ b_{P} & b_{P + 1} & \dots & b_{P + N} \\ y^{N} & y^{N - 1} & \dots & 1 \end{matrix} |} .

(8)

The Padé equations are equations (6) and (7), which generally determine the Padé numerator and denominator. Equation (8) is used to construct the [P/N] Padé approximant.³⁶

Definition (1): Differential transformation in one-dimension

If q(t) is analytic function in the domain T that is continuously differentiated with respect to time t, then

\frac{d^{h} q (t)}{d t^{h}} = ϕ (t, h), \forall t \in T,

(9)

for t = t_j, ϕ(t, h) = ϕ(t_j, h), where h belongs to the set of non-negative integers, denoted as H domain. Therefore, the differential transformation is defined as

Q_{j} (h) = ϕ (t_{j}, h) = {[\frac{d^{h} q (t)}{d t^{h}}]}_{t = t_{j}} \forall h \in H,

(10)

in the domain H, where Q_j(h) is the spectrum of q(t) at t = t_j.

Definition (2): Differential inverse transformation in one-dimension

Assume the following is a Taylor series for a function q(t)

q (t) = \sum_{h = 0}^{\infty} \frac{{(t - t_{j})}^{h}}{h!} Q (h) .

(11)

The inverse transformation of Q(h) is known as equation (11), and Q(h) is defined using definition (1)

Q (h) = N (h) {[\frac{d^{h} u (t) q (t)}{d t^{h}}]}_{t = t_{j}},

(12)

where h = 0, 1, 2, …, ∞. The function q(t) can thus be described as follows

q (t) = \frac{1}{u (t)} \sum_{h = 0}^{\infty} \frac{{(t - t_{j})}^{h}}{h!} \frac{Q (h)}{N (h)},

(13)

where N(h) ≠ 0, u(t) ≠ 0. The weighting factor is defined as N(h), and u(t) is treated as a kernel corresponding to q(t). Equations (11) and (13) are equal if N(h) = 1 and u(t) = 1. As a result, equation (11) can be considered as a special case of equation (13). The transformation with

N (h) = \frac{1}{h!}

and u(t) = 1 is used in this paper. Then, equation (12) becomes

Q (h) = \frac{1}{h!} [\frac{d^{h} q (t)}{d t^{h}}], where h = 0,1, \dots, \infty .

(14)

A differential equation in the domain of interest can be translated into an algebraic equation in the H domain using the differential transform, and q(t) can be derived as a finite-term Taylor series plus a remainder

\begin{array}{l} q (t) & = \frac{1}{u (t)} \sum_{h = 0}^{m} \frac{{(t - t_{0})}^{h}}{h!} \frac{Q (h)}{N (h)} + R_{m_{+} 1} (t) \\ = \sum_{h = 0}^{m} {(t - t_{0})}^{h} Q (h) + R_{m + 1} (t) . \end{array}

Stability

In this section, we study the stability of the Duffing equation

\ddot{q} + α \dot{q} + β q + γ q^{3} = 0,

(15)

where α, β, γ are constants. The equation is rewritten as following system of (q, q₁)

\dot{q} = q_{1},

(16)

q_{1}^{'} = - α q_{1} - β q - γ q^{3} .

(17)

The equilibrium points for the system are found by setting q′ = 0 and $q_{1}^{'} = 0$ . Thus, the equilibrium point for the system (16) is p₀ = (0, 0). The Jacobin matrix for the system at p₀ is

J (p_{0}) = (\begin{matrix} 1 & 0 \\ - α & - β \end{matrix}) .

The eigenvalues of the J(p₀) are λ₁ = 1 and λ₂ = −β, so p₀ is saddle point and unstable, but it can be stable for some initial conditions which means the neighbouring trajectories approach the equilibrium point only for certain initial conditions.^37,38

Applications

Example 1

Assume the damped Duffing oscillator equation with ICs as follows²¹

\ddot{q} + α \dot{q} + β q + γ q^{3} = 0, q (0) = w, \dot{q} (0) = z .

(18)

Applying the DTM and using Tables 1 gives the following scheme

\begin{array}{l} (h + 2) (h + 1) Q (h + 2) + α (h + 1) Q (h + 1) \\ + β Q (h) + γ S = 0, \\ Q (0) = W, Q (1) = Z, \end{array}

where

S = \sum_{x = 0}^{h} \sum_{m = 0}^{h - x} Q (x) Q (m) Q (h - x - m) .

The recursive equations of equation (18) for different values of h are as follows

\begin{array}{l} h = 0 : 2 Q (2) + α Q (1) + β Q (0) + γ (Q {(0)}^{3}) = 0, \\ h = 1 : 6 Q (3) + 2 α Q (2) + β Q (1) + γ (3 Q {(0)}^{2} Q (1)) = 0, \\ h = 2 : 12 Q (4) + 3 α Q (3) + β Q (2) + γ {(Q {(0)}^{2} Q (2) + 3 Q (0) Q (1))}^{2} = 0, \\ ⋮ \end{array}

Let consider the following value of parameters

α = 0.5, β = γ = 25, w = 0.1, z = 0.

(19)

For the given values in equation (19), q(t) is obtained as follows

\begin{array}{l} q (t) = 0.1 - 1.2625 t^{2} + 0.2104 t^{3} + 2.6828 t^{4} - 0.5392 t^{5} \\ - 2.6563 t^{6} + 0.6152 t^{7} . \end{array}

(20)

Table 1.

The one-dimensional differential transformation for basic operations.²⁹

Original function	Transformed function
q(t) ± v(t)	Q(h) ± V(h)
$\frac{d^{m} q (t)}{d t^{m}}$	$\frac{(h + m)!}{h!} Q (h + m)$
q(t)v(t)	$\sum_{l = 0}^{h} U (l) V (h - l)$
y ^m	$δ (h - m) = {\begin{array}{l} 1 & h = m \\ 0 & Otherwise \end{array}$
exp(λt)	$\frac{λ^{h}}{h!}$
sin(αt + β)	$\frac{α^{h}}{h!} \sin (\frac{h π}{2} + β)$
cos(αt + β)	$\frac{α^{h}}{h!} \cos (\frac{h π}{2} + β)$

To improve the result, we use the Laplace transform (L) on the series solution equation (20) as follows

\begin{array}{l} L [q (t)] = \frac{0.1}{s} - \frac{2.525}{s^{3}} + \frac{1.2625}{s^{4}} + \frac{64.3875}{s^{5}} \\ - \frac{64.7031}{s^{6}} - \frac{1912.5}{s^{7}} + \frac{3100.5}{s^{8}} . \end{array}

(21)

Replacing s by $\frac{1}{t}$ and then applying the Padé approximation of order three and four ( $P_{[\frac{3}{3}]}$ and $P_{[\frac{4}{4}]}$ ). Then, letting $t = \frac{1}{s}$ gives the following²⁸

P_{[\frac{3}{3}]} = \frac{0.1 s^{2} + 0.05 s + 0.05}{s^{3} + 0.5 s^{2} + 25.75 s},

(22)

P_{[\frac{4}{4}]} = \frac{0.1 s^{3} + 0.1667 s^{2} + 22.834 s + 11.4208}{s^{4} + 1.7 s^{3} + 253.6 s^{2} + 143.7 s + 5738.1} .

(23)

Applying the inverse Laplace transform (L⁻¹) into equations (22) and (23), we get the solutions by

P_{[\frac{3}{3}]}

and

P_{[\frac{4}{4}]}

, respectively, as

\begin{array}{l} q (t) = & 1.94 \times 10^{- 3} + 2.38 \times 10^{- 4} e^{- 2.5 \times 10^{- 1} t} \\ (411 \cos (5.068 t) + 20.273 \sin (5.068 t)) . \end{array}

\begin{array}{l} q (t) = & W e^{((- 06.0107 \times 10^{- 1} - 15.0816 i) t)} \\ + Z e^{((- 06.0107 \times 10^{- 1} + 15.0816 i) t)} \\ + A e^{((- 2.4894 \times 10^{- 1} - 5.0125 i) t)} \\ + B e^{((- 02.4894 \times 10^{- 1} + 5.0125 i) t)}, \end{array}

(24)

where

\begin{array}{l} W = 1.6932 \times 10^{- 5} - 1.3567 i \times 10^{- 4}, \\ Z = 1.6932 \times 10^{- 5} + 1.3567 i \times 10^{- 4}, \\ A = 4.9983 \times 10^{- 2} + 2.5633 i \times 10^{- 3}, \\ B = 4.9983 \times 10^{- 2} + 2.5633 i \times 10^{- 3} . \end{array}

The Padé approximation of $P_{[\frac{4}{4}]}$ yielded the real part of the solution in equation (24) as

\begin{array}{l} Re (y (t)) = e^{- 6.0107 \times 10^{- 1} t} (3.3846 \times 10^{- 5} \cos (15.0816 t) - 2.7134 \times 10^{- 4} \sin (15.0816 t)) \\ + 9.9966 \times 10^{- 2} e^{- 0.24894 t} \cos (5.0125 t) . \end{array}

Figure 1 presents the solution by PDTM for third and fourth order. The best result was found by applying PDTM of fourth order. Figure 2 shows the effect of the parameters α, β and γ in the solution. We realized that the frequency is changed by different α and β but does not change by γ.

Figure 1.

The comparisons of the results of $P_{[\frac{3}{3}]}$ , $P_{[\frac{4}{4}]}$ , and fourth-order Runge–Kutta.

Figure 2.

The effect parameters in Duffing equation.

Example 2

Consider the Duffing oscillator equation in the following type¹⁰

\ddot{q} + 3 q - 2 q^{3} = \cos (t) \sin (2 t),

(25)

that subjects to ICs

q (0) = 0, \dot{q} (0) = 1.

The exact solution of the given equation is

q (y) = \sin (y) .

(26)

By applying the DTM and using Tables 1, we have the following scheme

\begin{array}{l} (h + 2) (h + 1) Q (h + 2) + 3 Q (h) \\ - 2 \sum_{x = 0}^{h} \sum_{m = 0}^{h - x} Q (x) Q (m) Q (h - x - m) \\ = \frac{1^{h}}{h!} \cos (\frac{h π}{2}) \frac{2^{h}}{h!} \sin (\frac{h π}{2}) . \end{array}

Then, we apply the Padé approximation of order [3/3], [4/4] and N = 8 to obtain the solution as we see in Figure 3. The best solution is obtained by P_[4/4].

By plotting the solution of example (2) by DTM and by PDTM in Figures 4 and 5, respectively, we proved that the solution by PDTM converges to the exact solution in long domain, but the solution by DTM diverges after t = 3.

Figure 3.

The comparisons of the results of the $P_{[\frac{3}{3}]}$ , $P_{[\frac{4}{4}]}$ , and fourth-order Runge–Kutta.

Figure 4.

The comparisons of the solution by differential transform method with exact solution at t = 5, N = 8.

Figure 5.

The comparisons of the results of the $P_{[\frac{3}{3}]}$ , $P_{[\frac{4}{4}]}$ and exact solution.

Example 3

Let the non-dimensional Duffing oscillator equation^10,39

\ddot{q} + β q + γ q^{3} = | β - 1 | \cos (t),

(27)

which subjects to ICs

q (0) = 0, \dot{q} (0) = 0,

where β = 30 and γ = 0.1. By using the DTM, we obtain

\begin{array}{l} (h + 2) (h + 1) Q (h + 2) + β Q (h) + \\ γ \sum_{x = 0}^{h} \sum_{m = 0}^{h - x} Q (x) Q (m) Q (h - x - m) \\ = | β - 1 | \frac{1^{h}}{h!} \cos (\frac{h π}{2}) . \end{array}

Then, we apply P_[3/3] as well as as P_[4/4] and the results are presented in Figure 6.

Figure 6.

We note that $P_{[\frac{4}{4}]}$ is better than $P_{[\frac{3}{3}]}$ by the comparison with fourth-order Runge–Kutta.

However, we noticed that when the number of iteration is N = 8, then the best results are obtained when the order of the Padé approximation is [4/4]. The example 2 and 3 are non-damped, thus they have a periodic property. The following formula is the relationship between the frequency and the amplitude^40,41

ω^{2} = β + \frac{3}{4} γ A^{2},

where (A) is the amplitude. The frequency formulation of example 2 is

ω^{2} = 3 + \frac{3}{4} (- 2) A^{2} = 3 - \frac{3}{2} A^{2},

and in example 3 is

ω^{2} = 30 + \frac{3}{4} (0.1) A^{2} = 30 + \frac{3}{40} A^{2} .

In addition, for damped Duffing equation α ≠ 0

q^{''} (t) + α q^{'} + β q (t) = F_{0} \cos (ω t),

(28)

the practical resonance is

\sqrt{β - α^{2} / 2, β - α^{2} / 2 > 0}

and the pure resonance is

\sqrt{β}, β > 0

. Therefore, the practical and pure resonance are not equal and the general solution is the sum of two harmonic oscillations, hence it is bounded. In case of non-damped Duffing equation (α = 0), the practical and pure resonance are equal

\sqrt{β}, β > 0

. Thus, the solution is unbounded oscillatory solution as t → ∞.

Accuracy of the PDTM

In order to determine which order of the Padé approximation that satisfies the smallest error and since we do not have exact solution, the relative error takes the following formula

Relative Error = \frac{‖ Q - q ‖}{‖ Q ‖},

where q is the solution by the PDTM and Q is the solution by RK or the exact solution if it is available. Tables 2 and 3 present the errors, for example, 1 and 2, respectively. We noticed that the accurate results for homogeneous or non-homogeneous Duffing oscillator equation are obtained when the order of the Padé approximation is floor[N/2] where N is the number of iteration.⁴²

Table 2.

The accuracy of homogeneous Duffing oscillator equation (example 1).

N	$P_{[\frac{\frac{N}{2}}{\frac{N}{2}}]}$	$P_{[\frac{\frac{N}{2} - 1}{\frac{N}{2} - 1}]}$	$P_{[\frac{\frac{N}{2} - 2}{\frac{N}{2} - 2}]}$
8	2.9601 × 10⁻²	1.8208 × 10⁻¹	4.5741 × 10⁻²
9	2.9601 × 10⁻²	1.8208 × 10⁻¹	4.5741 × 10⁻²
10	2.9589 × 10⁻²	2.9601 × 10⁻²	1.8208 × 10⁻¹
11	2.9589 × 10⁻²	2.9601 × 10⁻²	1.8208 × 10⁻¹

Table 3.

The accuracy of non-homogeneous Duffing oscillator equation (example 2).

N	$P_{[\frac{\frac{N}{2}}{\frac{N}{2}}]}$	$P_{[\frac{\frac{N}{2} - 1}{\frac{N}{2} - 1}]}$	$P_{[\frac{\frac{N}{2} - 2}{\frac{N}{2} - 2}]}$
8	1.7173 × 10⁻²	1.7173 × 10⁻²	1.7173 × 10⁻²
10	1.7173 × 10⁻²	1.7173 × 10⁻²	1.7173 × 10⁻²
18	1.7173 × 10⁻²	7.2836 × 10⁻²	1.7173 × 10⁻²
20	1.7146 × 10⁻²	1.7173 × 10⁻²	7.2836 × 10⁻²

Conclusion

The article introduced the Padé differential transform approach to find a numerical solution for a non-homogeneous Duffing oscillator equation. The scheme of PDTM does need normalization or discretization and is simple and powerful method to find the numerical solution for the non-homogeneous class of ODEs. The accurate results are obtained when the order of the Padé approximation takes about half of the number of iteration. A comparison of the results with 4th-Runge–Kutta’s results reveals that the solution is extremely accurate. In future, the PDTM can be applied for variety of classes of ODEs.

Footnotes

Acknowledgements

The Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, Saudi Arabia has founded this preoject, under grant No. (KEP-MSc: 36-665-1443). The authors, therefore, acknowledge with thanks DSR technical and financial support.

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the The Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, Saudi Arabia has founded this project, under grant No.(KEP-MSc: 35-665-1443)

ORCID iD

Noufe H Aljahdaly

References

Salas

El-Tantawy

Aljahdaly

. An exact solution to the quadratic damping strong nonlinearity duffing oscillator. Math Probl Eng 2021; 2021: 1–8.

Aljahdaly

El-Tantawy

. On the multistage differential transformation method for analyzing damping duffing oscillator and its applications to plasma physics. Math 2021; 9(4): 432.

Aljahdaly

. New application through multistage differential transform method. AIP Conf Proc 2020; 2293: 420025.

Aljahdaly

Alharbey

. Fractional numerical simulation of mathematical model of hiv-1 infection with stem cell therapy. AIMS Math 2021; 6(7): 6715–6725.

Alqudah

Aljahdaly

. Global stability and numerical simulation of a mathematical model of stem cells therapy of hiv-1 infection. J Comput Sci 2020; 45: 101176.

Pirmohabbati

Sheikhani

Najafi

, et al. Numerical solution of fractional mathieu equations by using block-pulse wavelets. J Ocean Eng Sci 2019; 4(4): 299–307.

Aljahdaly

El-Tantawy

. Novel anlytical solution to the damped kawahara equation and its application for modeling the dissipative nonlinear structures in a fluid medium. J Ocean Eng Sci 2021.

Behera

Aljahdaly

Virdi

. On the modified (g’/g²)-expansion method for finding some analytical solutions of the traveling waves. J Ocean Eng Sci 2021.

Al-Jawary

Abd-Al-Razaq

. Analytic and numerical solution for duffing equations. Int J Basic Appl Sci 2016; 5(2): 115–119.

10.

Yusufoğlu

. Numerical solution of duffing equation by the laplace decomposition algorithm. Appl Math Comput 2006; 177(2): 572–580.

11.

Lenci

Rega

. Forced harmonic vibration in a duffing oscillator with negative linear stiffness and linear viscous damping. In: Kovacic

Brennan

(eds). The duffing equation: nonlinear oscillators and their behaviour. Hoboken, NJ: Wiley, 2011, pp. 219–276.

12.

Salas

Castillo H

. Exact solution to duffing equation and the pendulum equation. Appl Math Sci 2014; 8(176): 8781–8789.

13.

Khalid

Sultana

Arshad

, et al. A comparison between new iterative solutions of non-linear oscillator equation. Int J Computer Appl 2015; 128(4): 1–5.

14.

Nadeem

Yao

. Solving system of partial differential equations using variational iteration method with he’s polynomials. J Math Computer Sci 2019; 19(3): 203–211.

15.

Nadeem

. He–laplace variational iteration method for solving the nonlinear equations arising in chemical kinetics and population dynamics. J Math Chem 2021; 59(5): 1234–1245.

16.

Khan

Vázquez-Leal

Faraz

. An efficient new iterative method for oscillator differential equation. Scientia Iranica 2012; 19(6): 1473–1477.

17.

Habib

Islam

Batool

, et al. Numerical solutions of the fractal foam drainage equation. GEM-International J Geo Math 2021; 12(1): 1–10.

18.

Turkyilmazoglu

. An effective approach for approximate analytical solutions of the damped duffing equation. Physica Scripta 2012; 86(1): 015301.

19.

Bülbül

Sezer

. Numerical solution of duffing equation by using an improved taylor matrix method. J Appl Math 2013; 2013: 1–6.

20.

Sheu

Chen

, et al. Chaotic dynamics of the fractionally damped duffing equation. Chaos, Solitons & Fractals 2007; 32(4): 1459–1468.

21.

Nourazar

Mirzabeigy

. Approximate solution for nonlinear duffing oscillator with damping effect using the modified differential transform method. Scientia Iranica 2013; 20(2): 364–368.

22.

Pukhov

. Computational structure for solving differential equations by taylor transformations. Cybernetics 1978; 14(3): 383–390.

23.

Chen

. Solving partial differential equations by two-dimensional differential transform method. Appl Math computation 1999; 106(2–3): 171–179.

24.

Hassan

IAH

. Application to differential transformation method for solving systems of differential equations. Appl Math Model 2008; 32(12): 2552–2559.

25.

Hassan

IAH

. Comparison differential transformation technique with adomian decomposition method for linear and nonlinear initial value problems. Chaos, Solitons & Fractals 2008; 36(1): 53–65.

26.

Hassan

IAH

. Differential transformation technique for solving higher-order initial value problems. Appl Math Comput 2004; 154(2): 299–311.

27.

Hatami

Jing

. Differential transformation method for newtonian and non-newtonian nanofluids flow analysis: Compared to numerical solution. Alexandria Eng J 2016; 55(2): 731–739.

28.

Momani

Ertürk

. Solutions of non-linear oscillators by the modified differential transform method. Comput Math Appl 2008; 55(4): 833–842.

29.

Abazari

Soltanalizadeh

. Reduced differential transform method and its application on kawahara equations. Thai J Math 2012; 11(1): 199–216.

30.

Tamboli

Tandel

. Solution of the time-fractional generalized burger–fisher equation using the fractional reduced differential transform method. J Ocean Eng Sci 2021.

31.

Yamada

Ikeda

. A numerical test of padé approximation for some functions with singularity. Int J Comput Math 2014; 2014: 1–17.

32.

Moatimid

Zekry

. Forced nonlinear oscillator in a fractal space. Facta Universitatis, Ser Mech Eng 2022; 20(1): 001–020.

33.

Zheng

. Approximate solutions to mhd falkner-skan flow over permeable wall. Appl Math Mech 2011; 32(4): 401–408.

34.

Baker

Jr . Essentials of Padé approximants. Amsterdam, Netherlands: Elsevier, 1975.

35.

Nadeem

Yao

. Solving the fractional heat-like and wave-like equations with variable coefficients utilizing the laplace homotopy method. Int J Numer Methods Heat Fluid Flow 2020; 31: 273–292.

36.

Sayyed

Singh

Bano

. Analytical solution of mhd slip flow past a constant wedge within a porous medium using dtm-padé. Appl Math Comput 2018; 321: 472–482.

37.

Tian

Ain

Anjum

, et al. Fractal n/mems: from pull-in instability to pull-in stability. Fractals 2021; 29(02): 2150030.

38.

Tian

. From inner topological structure to functional nanofibers: Theoretical analysis and experimental verification. Membranes 2021; 11(11): 870.

39.

Buchanan

Turner

. Numerical methods and analysis. New York, NY: McGraw-Hill, 1992.

40.

. On the frequency-amplitude formulation for nonlinear oscillators with general initial conditions. Int J Appl Comput Mathematics 2021; 7(3): 1–5.

41.

Liu

. A modified frequency-amplitude formulation for fractal vibration systems. Fractals 2021; 30: 2250046.

42.

Zhang

Nadeem

Rauf

, et al. A novel approach for the analytical solution of nonlinear time-fractional differential equations. Int J Numer Methods Heat Fluid Flow 2020; 31: 1069–1084.