Abstract

This paper presents some novel solutions to the family of the Helmholtz equations (including the constant forced undamping Helmholtz equation (equation (1)) and the constant forced damping Helmholtz equation (equation (2))) which have been reported. In the beginning, equation (1) is solved analytically using two different techniques (direct and indirect solutions): in the first technique (direct solution), a new assumption is introduced to find the analytical solution of equation (1) in the form of the Weierstrass elliptic function with arbitrary initial conditions. In the second case (indirect solution), the solution of the undamping (standard) Duffing equation is devoted to determine the analytical solution to equation (1) in the form of Jacobian elliptic function with arbitrary initial conditions. Moreover, equation (2) is solved using a new ansatz and with the help of equation (1) solutions. Also, the evolution equations (equations (1) and (2)) are solved numerically via the Adomian decomposition method (ADM). Furthermore, a comparison between the approximate analytical solution and approximate numerical solutions using the fourth-order Runge–Kutta method (RK4) and ADM is reported. Furthermore, the maximum distance error for the obtained solutions is estimated. As a practical application, the Helmholtz-type equation will be derived from the fluid governing equations of quantum plasma particles with(out) taking the ionic kinematic viscosity into account for investigating the characteristics of (un)damping oscillations in a degenerate quantum plasma model.

1. Introduction

The ordinary and partial differential equations have played an important role in explaining many natural phenomena, in addition to their applications in many engineering and physical problems. Due to the great role played by these equations, many authors focused their efforts on finding some solutions to these equations [110]. The undamping Helmholtz equation and undamping Duffing equation in addition to their family (including friction/damping force in addition to excitation/perturbation force) are among the most famous differential equations in dynamic, electrical, and engineering systems [1116]. A lot of physical and engineering problems such as the human eardrum oscillations, the dynamics of the ships, the electrical circuits signal oscillations, heavy symmetric gyroscope, and microperforated panel absorber [1721] have been investigated using different solutions of the Helmholtz-type oscillator. The Helmholtz-type oscillator is a second-order differential equation with a quadratic nonlinear term in addition to some other terms. For realistic physical situations, we cannot ignore both the friction/dissipation force and excitation/perturbation forcing. Accordingly, the general form of the Helmholtz equation reads [2224]where denotes the displacement of the system, is the natural frequency, is a nonlinear system parameter, represents the damping factor, and F is a constant/excitation force. The first equation in system (1) is called the constant forced damping Helmholtz equation (equation (2)). If the coefficient of the damping term is neglected, equation (2) reduces to the constant forced Helmholtz equation (equation (1)). Also, if both the coefficient of the damping term and the constant force are neglected, the traditional form of Helmholtz equation is covered . The initial value problem (IVP) (1) and its family have many applications in several fields, starting from analyzing the signals that propagate in electrical circuits, plasma physics, general relativity, betatron oscillations, vibrations of shells, vibrations of the acoustically driven human eardrum, solid-state physics, etc. [2633].

It is well known that, in the absence of both friction and the excitation forces from the IVP (1), the unforced and undamping Helmholtz equations are covered. The exact analytic solutions to the unforced and undamping Helmholtz equation have been derived in detail in the literature in terms of the Weierstrass elliptic function [3438] and Jacobi elliptic functions [3841]. Generally, to solve any quadratic or cubic nonlinear second-order differential equation, firstly, we should transform it to an elliptic integral and then we solve it [24]. It is known that the unforced and undamping Helmholtz equations are completely integrable, so they have exact solutions, but if the friction force (damping term) is included, then the unforced damping Helmholtz equation becomes nonintegrable and cannot support an exact solution for arbitrary values of its coefficients . Thus, under certain condition, the unforced damping Helmholtz equation has been solved analytically in terms of the Jacobi elliptic functions by Johannessen [24]. Also, Almendral and Sanjuán [27] derived an exact solution to the unforced damping Helmholtz equation using the Lie theory under certain conditions for the coefficients .

In many realistic physical models, both the damping and the excitation/external terms are very important to be included, and thus the problem becomes more complicated to find its analytical solutions. In this paper, we will derive some analytical solutions to equation (1) in the terms of the Weierstrass and Jacobian elliptic functions. Also, an approximate analytical solution to equation (2) for arbitrary values to the coefficients and the initial conditions will be derived in detail. Moreover, the problem under consideration will be solved numerically via using the RK4 and ADM to make a comparison between the obtained solutions and the approximate numerical solutions. Furthermore, the maximum distance error between the approximate analytical solution and the approximate numerical solutions will be estimated. Also, the dynamics of nonlinear oscillations that can be generated in the RLC electronic circuits and quantum plasma will be investigated using the solution of equation (2).

The rest of this work is organized in the following manner: in Section 2, we will introduce in detail our methodology for solving the family of the Helmholtz-type equations. Also, we will introduce some new approaches for solving equation (1) as well as the exact solution of the undamped and unforced Duffing equation, which will be devoted to finding an approximate analytical solution to equation (2). In Section 3, a comparison between the obtained solutions and the approximate numerical solution using the ADM will be investigated. In Section 4, some realistic applications will be introduced. Finally, our results will be summarized in Section 5.

2. Our Methodology for Solving the Family of the Helmholtz Equations

Before proceeding in solving the IVP (1), it is necessary to refer to two fundamental equations and their solutions: the first one is called the Duffing equation and the other is called the constant forced Helmholtz equation.

2.1. Duffing Equation and Its Solution

The analytical solution to the following IVP, which is called Duffing equation [30],is given by the following formula:

By inserting this relation into the IVP (2) and after several tedious but simple math operations, we finally get the values of and as follows:where , , , and are real numbers and is called the discriminant to Duffing (2):

Solution (3) could be expressed aswhere

For negative discriminant , the solution may be written in the following form:where

For a zero discriminant , the solution of Duffing equation (2) will bewithwhere .

When , the solution reads

In case , the solution becomes the constant function .

2.2. The Analytical Solution to the Constant Forced Helmholtz Equation

The solution of the constant forced Helmholtz equationmay be expressed in either one of the following forms. Note that .

2.2.1. First Formula

Suppose that the solution of system (13) has the following form:as well aswhere represents the constant force and gives the Weierstrass elliptic function, which satisfies the following relations:

Inserting equation (14) into system (13) and after tedious straightforward calculations, we can estimate the values of , , , , , and as follows:

The value of parameter represents a root to the following quartic equation:

2.2.2. Second Formula

As another form for the solution of system (13), let us assume that its solution is given by the following relationship:where is the solution to the following Duffing equation:

By following the same procedures in the above sections, we finally obtain

The value of parameter is a solution to the following equation:

Also, the value of parameter is a solution to the following equation:

2.3. The Approximate Analytical Solution for the Constant Forced and Damped Helmholtz Equation

Let us rewrite system (1) in the following traditional form:with .

Also, let us assume that

Now, suppose that the solution of system (24) is given bywhere is a solution to the following Helmholtz equation:

Inserting (26) into the first equation in system (24), , we obtain

Expression (28) suggests the following choices:giving us the values of and as follows:

The solution to the following IVP,may be expressed in different forms as we explained in the previous section.

3. A Comparison between Our Solutions and ADM Solution

There are many numerical methods that could be used to find approximate solutions to IVP (24). Here, we make use of the ADM [4244] to solve IVP (24) with(out) forcing term . According to this method, the first iteration/approximation for the unforced case readswhere , , and .

For the second iteration/approximation, we have

The approximate Adomian approximate solution reads

For and , we can compare between the approximate/semianalytical solution (26) (for ) and the ADM approximate numerical solution (34) and the RK4 approximate numerical solution as shown in Figure 1. Also, the maximum distance error according to the formula has been estimated:

It is clear that the semianalytical solution (26) (for ) gives good results as compared to both the RK4 and ADM approximate numerical solutions.

Now, let us find an approximate solution for IVP (24) in the presence of the forcing term, using the ADM. Accordingly, the first approximation is given bywhere , , , , and .

The second approximation to IVP (24) according to the ADM readswhere .

The approximate Adomian approximate solution is given by

For and , a comparison between the semianalytical solution (26) (for ) and the ADM approximate numerical solution (38) and the RK4 approximate numerical solution has been investigated as shown in Figure 2. Furthermore, the maximum distance error has been calculated as follows:

Also, the semianalytical solution (26) for and gives excellent results as compared to the ADM approximate numerical solution (38).

4. Quantum Plasma Oscillations

In this section, we will reduce the fluid governing equations of a quantum plasma model to an evolution equation using the RPT [4549]. After a suitable transformation, we will be able to convert the obtained evolution equation to a Helmholtz-type equation in order to investigate the characteristics behavior of the damping oscillations in the model under consideration. Now, let us assume that we have a collisionless and unmagnetized electron-ion quantum plasma consisting of inertialess degenerate trapped electrons which obey the Fermi–Dirac distribution and classical fluid cold positive nondegenerate ion. Thus, the basic normalized fluid equations that govern the nonlinear dynamics of various structures could be presented as [51, 52]where and represent the normalized electron and ion number densities, respectively, gives the normalized ion speed, is the normalized kinematic viscosity of the ions, and indicates the normalized electrostatic wave potential. Therefore, we shall adopt the adiabatic trapped degeneracy for electrons, by relying on notations similar to those in [52], wherein the fundamental algebra is expressed in detail. The electron normalized number density according to Fermi–Dirac distribution readswhere expresses the normalized temperature of the degenerate electron, , , , and . Note that expression (41) is obtained under the approximation for small wave amplitude.

For investigating the nonlinear structures and oscillations in the present model, the RPM will be employed for this purpose. Accordingly, the stretching and expansions for the independent and dependent variables are, respectively, introduced as follows:where is a small and real parameter and represents the normalized phase velocity of the unmodulated structures. It is assumed that the impact of ionic kinematic viscosity is small according to many lab experiments [53, 54]. Thus, we could set , where is . By substituting stretching (42) and expansions (43) into system (40) and by collecting the terms of different powers of , we could get a system of reduced equations. Solving the system of reduced equations for the first-two orders of by following the same procedures in [4649], we finally obtain the Korteweg–de Vries Burgers (KdVB) equation [50].withwhere , , and represent the coefficients of the nonlinear, dispersion, and dissipative terms, respectively, and .

Using the traveling wave transformation , where and gives the frame velocity, into the KdVB equation (44) and integrating once over , we get the constant forced and damped Helmholtz equation as follows:where is the integration constant, , , and . Now, we can apply the above solution of the constant forced and damped Helmholtz equation that is given in equation (26) to equation (46) for investigating the characteristics of the damped oscillations in a quantum plasma.

Note that if the ionic kinematic viscosity is neglected, i.e., , then the KdV and undamped Helmholtz equations could be covered. The KdV equation, i.e., equation (44) for , is one of the most popular soliton and cnoidal equations and has been extensively investigated. For the soliton solution, the following conditions must be fulfilled: at , so the integration constant in equation (46) must vanish. Accordingly, equation (46) could be reduced to the undamped Helmholtz equation:

It is well known that equation (47) supports many solutions such as periodic solution (see the solutions to the constant forced Helmholtz equation above) and solitons. The soliton solution to equation (47) in the form of the Weierstrass elliptic function could be written in the following manner:where .

Here, the obtained approximate analytical solution (26) to the constant forced and damped Helmholtz equation (46) will be analyzed numerically according to the quantum plasma parameters , i.e., , and , i.e., . The behavior of the quantum plasma oscillations according to the approximate analytical solution (26) and the approximate numerical solution according to the RK4 method is presented in Figure 3 for different values of quantum plasma parameters. It is clear from Figure 3 that our approximate analytical solution (26) is more accurate than the RK4 numerical solution. On the contrary, the RK4 numerical solution gives poor results and with increasing time this solution becomes unstable.

5. Conclusions

The Helmholtz-type equations including the constant forced undamping Helmholtz equation (equation (1)) and the constant forced damping Helmholtz equation (equation (2)) have been solved analytically and numerically. Two techniques were used to get the analytical solutions to equation (1). In the first technique, we used a new assumption to find an analytical solution to equation (1) in the form of Weierstrass elliptic function. In the second case, the solution of the standard Duffing equation has been utilized to find an analytical solution to equation (1) in the form of Jacobian elliptic function. However, the main goal of this paper is to solve equation (2), using the obtained solutions of equation (1). Moreover, both equation (1) and equation (2) have been solved numerically via the ADM. The analytical and approximate analytical solutions of equations (1) and (2) have been compared to the RK4 and ADM approximate numerical solutions. Furthermore, the maximum distance error between the RK4 approximate numerical solution and the approximate analytical solutions in addition to the approximate numerical solution using the ADM has been estimated. It was found that the obtained solutions are generally consistent with both RK4 and ADM solutions. Moreover, the obtained solutions have been applied for analyzing the oscillations that may arise in the quantum plasma. During the analysis, it was found sometimes that the approximate analytical solution is better than the RK4 numerical solution as shown in the quantum plasma model. Finally, these solutions may help us understand the oscillations that may arise in the different physical and engineering systems.

In future work, the similar approaches could be used for analyzing and solving higher-order nonlinear oscillator equations. Also, a damping Helmholtz–Duffing equation with time-dependent forced term is considered one of the most important and vital problems due to its great role in explaining many natural phenomena in different branches of science. Thus, in the next work, some new approaches will be devoted to find some solutions for these problems.

Data Availability

No data were used to support the findings of this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by Taif University Researchers Supporting Project (no. TURSP-2020/275), Taif University, Taif, Saudi Arabia. The authors thank Associate Professor S.A. El-Tantawy, the Chief of Research Center for Physics (RCP), Faculty of Science and Arts, Al-Baha University, Saudi Arabia, and Faculty of Science, Port Said University, Egypt, for fruitful discussions, careful reading, his contributions to applications to plasma physics, and improving their manuscript very well.