On the Approximate Solutions of the Constant Forced (Un) Damping Helmholtz Equation for Arbitrary Initial Conditions

$is 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 fourthorder 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.


Introduction
e 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 [1][2][3][4][5][6][7][8][9][10]. e undamping Helmholtz equation (€ q(t) + n i�1 k i q i (t) � 0, where i is the odd number, i.e., i � 1, 3, 5, . . .) and undamping Duffing equation (€ q(t) + n i�1 k i q i (t) � 0, where i � 1, 2, 3, . . .) 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 [11][12][13][14][15][16]. 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 [17][18][19][20][21] have been investigated using different solutions of the Helmholtztype oscillator. e Helmholtz-type oscillator is a secondorder 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 [22][23][24] € q + 2c _ q + αq + βq 2 � F, q(0) � q 0 & q ′ (0) � _ q 0 , (1) where q denotes the displacement of the system, α is the natural frequency, β is a nonlinear system parameter, c represents the damping factor, and F is a constant/excitation force. e first equation in system (1) is called the constant forced damping Helmholtz equation (equation (2)). If the coefficient c of the damping term (2c _ q) is neglected, equation (2) reduces to the constant forced Helmholtz equation (€ q + αq + βq 2 � F) (equation (1)). Also, if both the coefficient c of the damping term and the constant force F are neglected, the traditional form of Helmholtz equation is covered (€ q + αq + βq 2 � 0). e 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. [26][27][28][29][30][31][32][33].
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. e 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 [34][35][36][37][38] and Jacobi elliptic functions [38][39][40][41]. 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 (c, α, β).
us, 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 (c, α, β).
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). e 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.

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.

Duffing Equation and Its Solution.
e 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 c 1 and c 2 as follows: where R, S, η 0 , and _ η 0 are real numbers and Δ is called the discriminant to Duffing (2): Solution (3) could be expressed as where 2 Mathematical Problems in Engineering For negative discriminant (Δ < 0), the solution may be written in the following form: where For a zero discriminant (Δ � 0), the solution of Duffing equation (2) will be with where In

e Analytical Solution to the Constant Forced Helmholtz Equation.
e solution of the constant forced Helmholtz equation may be expressed in either one of the following forms. Note that a � − F.

First Formula.
Suppose that the solution of system (13) has the following form: as well as where a represents the constant force and ℘(t) 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 B, C, g 2 , g 3 , d 1 , and d 2 as follows: e value of parameter A represents a root to the following quartic equation: Mathematical Problems in Engineering

Second Formula.
As another form for the solution of system (13), let us assume that its solution is given by the following relationship: where η � η(t) is the solution to the following Duffing equation: By following the same procedures in the above sections, we finally obtain e value of parameter A is a solution to the following equation: Also, the value of parameter B is a solution to the following equation:

e Approximate Analytical Solution for the Constant Forced and Damped Helmholtz Equation.
Let us rewrite system (1) in the following traditional form: with β ≠ 0. Also, let us assume that Now, suppose that the solution of system (24) is given by where y � y(t) is a solution to the following Helmholtz equation: y ″ + py + βy 2 � 0, Inserting (26) into the first equation in system (24), Expression (28) suggests the following choices: giving us the values of p and d as follows: e solution to the following IVP, 4 Mathematical Problems in Engineering may be expressed in different forms as we explained in the previous section.

A Comparison between Our Solutions and ADM Solution
ere are many numerical methods that could be used to find approximate solutions to IVP (24). Here, we make use of the ADM [42][43][44] to solve IVP (24) with(out) forcing term (F). According to this method, the first iteration/approximation for the unforced case (F) reads where Θ 1 � sin(tF), Θ 2 � cos(tF), and F � ����� α − ε 2 √ . For the second iteration/approximation, we have e approximate Adomian approximate solution reads For (ε, α, β, F) � (0.1, 4, 1, 0) and x(0) � 0 & x ′ (0) � 0.2, we can compare between the approximate/semianalytical solution (26) (for F � 0) 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 L D � max t i ≤t≤t f |x RK4 (t) − x method (t)| has been estimated: It is clear that the semianalytical solution (26) (for F � 0) 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 by where Θ 1 � sin(tF), . e second approximation to IVP (24) according to the ADM reads

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 [45][46][47][48][49]. 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. us, the basic normalized fluid equations that govern the nonlinear dynamics of various structures could be presented as [51,52] where n e and n represent the normalized electron and ion number densities, respectively, u gives the normalized ion speed, η is the normalized kinematic viscosity of the ions, and ϕ indicates the normalized electrostatic wave potential. erefore, 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. e electron normalized number density according to Fermi-Dirac distribution reads where T expresses the normalized temperature of the degenerate electron, s 0 � (1 + T 2 ), s 1 � (3 − T 2 )/2, s 2 � 3(1 + T 2 )/8, and s 3 � − (1 + 5T 2 )/6. Note that expression (41) is obtained under the approximation ϕ ≪ 1 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: n(x, t) � 1 + εn (1) + ε 2 n (2) + ε 3 n (3) + · · · , u(x, t) � εu (1)  where ε is a small (ε ≪ 1) and real parameter and V ph 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]. us, we could set η � ε 1/2 η, where η is O(1). 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 [46][47][48][49], we finally obtain the Korteweg-de Vries Burgers (KdVB) equation [50]. with where A p , B p , and C p represent the coefficients of the nonlinear, dispersion, and dissipative terms, respectively, and ϕ ≡ ϕ (1) .
Using the traveling wave transformation ϕ(ξ + V f t) � q(ζ), where ζ � ξ + V f t and V f 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 C is the integration constant, c � − C p /(2B p ), α � V f /B p , and β � A p /(2B p ). 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., C p � 0, then the KdV and undamped Helmholtz equations could be covered. e KdV equation, i.e., equation (44) for C p � 0, 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: (q(ζ), q ′ (ζ), q ″ (ζ)) ⟶ 0 at ζ ⟶ ± ∞, so the integration constant C 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. e soliton solution to equation (47) in the form of the Weierstrass elliptic function ℘ could be written in the following manner: where (x, t) ≡ (ζ, τ).
Here, the obtained approximate analytical solution (26)  e 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.
Mathematical Problems in Engineering 7

Conclusions
e 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. e 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 timedependent 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. us, 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
e authors declare that they have no conflicts of interest.