Investigation of the Surface Brightness Model in the Milky Way via Homotopy Perturbation Method ()
1. Introduction
The Ambartsumian equation is used in the theory of surface brightness in the Milky Way [1]. It is a linear delay differential equation given by [2]
(1)
where q is a constant and
(2)
where
is also a constant. Existence and uniqueness were discussed in [3]. Although the Adomian decomposition method (ADM) was effective to deal with various types of equations [4] - [18], the HPM55 [19] [20] is preferred here to analyze (1-2). Details of the HPM were introduced by Ayati and Biazar [20]. Moreover, it will be shown that the present analysis posses more accuracy over the previous method in the literature.
2. Application of the HPM
First, Equation (1) is rewritten as
(3)
where
(4)
On substituting (4) into (3), we have
(5)
which yields
(6)
and
(7)
Hence,
(8)
From (7) and (8), the 1st-order system is given as
(9)
Therefore
(10)
The 2nd-order system is given by
(11)
By solving the system (11) for
, we have
(12)
Proceeding as above we obtain the 3rd-order system as
(13)
with the corresponding solution
(14)
The calculated higher-order solutions are obtained by MATHEMATICA and then implemented to producing the results in Section 5. From Equation (4), the HPM gives the series solution of Equation (1) as
(see Ayati and Biazar [20]) by
. This infinite series is approximated by replacing infinity with n-term, hence, the approximate solution, denoted by
, is given by [19] [20]:
(15)
The residual
is given as
(16)
and the approximate solution in Ref. [2] is
(17)
3. Analytic Solution
Equation (1) can be written as
(18)
Assume that
(19)
Accordingly, we have
(20)
and
(21)
Inserting Equations (19)-(21) into Equation (18), we get
(22)
i.e.,
(23)
and this yields
(24)
The last equation implies that
(25)
Accordingly,
(26)
Hence
(27)
The initial condition (2), gives
by
(28)
Thus
(29)
and the convergence of series in (29) is discussed in the next section in detail.
4. Convergence Analysis
Definition 1: Let
be a sequence of real functions, each function of which is defined for all x on a real interval
. For each particular x such that
consider the corresponding sequence of real numbers
. Suppose that the sequence
converges for every x such that
, and let
. Then we say that the sequence of real functions
converges pointwise on the interval
, and the function
thus defined i called the limit function of the sequence.
Definition 2: Let be a sequence of real functions, each function of which is defined for all x on a real interval. The sequence is said to converge uniformly to on if, given any, there exists (which depends only upon) such that .
Theorem 1: Let be a sequence of real functions converges uniformly to on and suppose that each function is continuous on, then the limit function is continuous on.
Theorem 2: Weierstrass M-Test
1) Let be a sequence of positive constants such that the series of constants converge.
2) Let be a series of real functions such that for each.
Then the series converges uniformly on.
Theorem 3: From (29), Let, then is a sequence of positive constants and the series converges.
Proof: Since and, then is a positive constant and. Also, the expression is always positive, where. Therefore, is a sequence of positive constants. To prove convergence of the series, we have from the ratio test that
(30)
which proves the convergence of the series.
Theorem 4: The solution given by (29) converges uniformly on the whole domain.
Proof: From the previous theorem, it we showed that is a sequence of positive constants and the series converges. This meets the first requirement of the Weierstrass M-Test in theorem 2 above. Since converges, we suppose that its sum equals, i.e.,. In order to satisfy the second requirement of theorem 2, we rewrite the solution (29) in terms of a new variable x, where, as:
(31)
To prove convergence of in (31), it is sufficient to prove the convergence of the series. At this stage, we have from that
(32)
where , and this completes the proof.
5. Validation of Numerical Results
The HPM and another direct approach have been applied in the previous sections to obtain the approximate solutions and the analytic solution, respectively, in terms of exponential functions with negative powers. The convergence of the analytic solutions was discussed in Section 4. The obtained approximate solutions by the HPM are to be analyzed here in view of the analytic solution (29) to stand on their accuracy. This can be achieved via performing comparisons with the results in the literature. The comparisons between the present results and those of Ref. [2] and [21] are presented in Table 1. The present HPM is of higher accuracy as observed from Table 1. This indicates the advantages of the current approach over those in the literature when analyzing the Ambartsumian equation.
The behavior of is displayed in Figure 1 (), Figure 2 (), Figure 3 (), and Figure 4 () for. It can be seen from these figures that the maximum values of are, , , and, respectively. This proves the efficiency of the HPM over the previous ones in [2] and [21].
6. Conclusion
The HPM was applied to solve the Ambartsumian equation in terms of exponential functions. The obtained solution was valid in the whole domain, while the corresponding solution in the literature [2] was only valid in sub-domains.
Table 1. Comparison of the present results with the corresponding results in the literature.
Moreover the obtained residual tends to zero as the q increases. In view of references [2] and [21], the present HPM is of higher accuracy when compared with the obtained exact solution.