The Zhou’s Method for Solving the White-Dwarfs Equation

Abstract

In this work we apply the differential transformation method (Zhous method) or DTM for solving white-dwarfs equation which Chandrasekhar [1] introduced in his study of the gravitational potential of these degenerate (white-dwarf) stars. DTM may be considered as alternative and efficient for finding the approximate solutions of the initial values problems. We prove superiority of this method by applying them on the some Lane-Emden type equation, in this case. The power series solution of the reduced equation transforms into an approximate implicit solution of the original equation.

Share and Cite:

P. Alzate and W. Urueña, "The Zhou’s Method for Solving the White-Dwarfs Equation," Applied Mathematics, Vol. 4 No. 10C, 2013, pp. 28-32. doi: 10.4236/am.2013.410A3005.

1. Introduction

Other classical nonlinear equation which has been the object of much study is called the white-dwarf equation. This equation has the form:

(1)

with and subject to initial conditions:

(2)

The Equation (1) is one of Lane-Emden types, where

. In fact, it reduces to Lane-Emden equation with index when.

Let us consider a spherical cloud of gas (see Figure 1) and denote its hydrostatic pressure at a distance from the centre by P. Let be the mass of the spheres of radius the gravitational potential of the gas, and g the acceleration of gravity.

Then, we have the following equations:

(3)

where G is the gravitational constant.

Thus, three conditions are assumed for the determination of and P:

(4)

where ρ is the density of the gas.

(5)

and

(6)

where γ and K are arbitrary constants.

Figure 1. Spherical cloud of gas.

Solving (4) and (6) with when we have:

(7)

or

(8)

where and. If this value of ρ is replaced into Equation (5), we obtain:

(9)

where,.

This is the Lane-Emden equation. This is a second order ordinary differential equation for the density profile ϕ as a function of the radius. We also have to determine initial conditions at. Obviously, at the center the density has to be the central density, and hence at. Also, the derivative of the density has to vanish at the center, and therefore at.

Our goal is to find a function ϕ that satisfies the Lane-Emden Equation (9). In general, this has to be done numerically using the differential transformation method.

The Zhou’s Method is a semi-numerical-analytic method for solving ordinary and partial differential equations. The concept of the DTM was first introduced by Zhou in 1986 [2]. Its main application therein is to solve both linear and non-linear initial value problems in electric circuit analysis. This technique constructs an analytical solution in the form of a polynomial. The Taylor series method is computationally expensive for large orders. The differential transformation method is an alternative procedure for obtaining analytic Taylor series solution of the differential equations. The series often coincides with the Taylor expansion of the true solution at point, in the initial value case, although the series can be rapidly converged in a very small region.

Many numerical methods were developed for this type of nonlinear ordinary differential equations, specifically on Lane-Emden type equations such as the Adomian Decomposition Method (ADM) [3,4], the Homotopy Perturbation Method (HPM) [5,6], the Homotopy Analysis Method (HAM) [7] and Bernstein Operational Matrix of Integration [8], in [9] Hojjati and Parand propose a class of second derivative multistep methods (SDMM) for solving some well-known classes of Lane-Emden type equations. Finally in [10] Batiha propose the variational iteration method (VIM) for a class of Lane-Emden type equation, powerful method for the solution of linear and nonlinear equations. In this paper, we show superiority of the DTM by applying them on the some type LaneEmden type equations (white-dwarfs).

2. Differential Transformation Method

Differential transformation method of function is defined as follows:

(10)

In (10), is the original function and is the transformed function and the inverse differential transformation is defined as:

(11)

In real applications, function is expressed by a finite series and Equation (11) can be written as:

(12)

Equation (12) implies that

is negligibly small. Here, n is decided by the convergence of natural frequency in this study.

The following theorems that can be deduced from Equations (10) and (11).

Theorem 1.1. If

then

Theorem 1.2. If

then

is a constant.

Theorem 1.3. If

then

Theorem 1.4. If

, then

Theorem 1.5. If

thenwhere

Theorem 1.6. (Cárdenas A). If with then:

Proof. We proof this theorem by induction method. We assume that is true for m, then we will prove that is true also for. Using DTM to a product we have:

then,

where and Therefore,

whenever the expression reduces to:

Now, if then all terms of the sum are zero except $ obtaining:

Rearranging the expression we have:

The proofs of Theorems 1-5 are available in [1].

3. Numerical Solution for C = 0.3

To illustrate the ability of the Zhou’s method for the Lane-Emden type equation (white-dwarf), the next problem is provided for. The results reveal that this method is very effective.

Problem. To solve subject to and. First, multiplying both sides by x we have:

(13)

Here, is easy to verify that the function

has a series expansion:

(14)

where Therefore, (13) takes the form:

(15)

Using in (15) the Theorems 3, 4 and 6 we have:

(16)

or

(17)

where,

and successively. Also,

(18)

(19)

(20)

for all.

Now, from the initial conditions and we have:

(21)

Substituting Equation (21) into Equation (17) and by recursive method, the results are listed as follows:

For we have:

For we have:

and then,.

For we have:

and then,

For we have:

and then,

For we have:

and then,

Therefore using (12), the closed form of the solution can be easily written as:

A series solution obtained by Chandrasekhar [2] using series expansion was:

Table 1 shows the comparison of obtained by the DTM (method proposed in this work) and those obtained by Hojjati. The resulting graph of the whitedwarfs equation in comparison to the present method and those obtained by [5] is shown in Figure 2.

4. Conclusion

In this work, we presented the definition and handling of one-dimensional differential transformation method. Using the differential transformation, differential equations can be transformed to algebraic equations and the resulting algebraic equations are called iterative equations. This method has applied to solve some class of LaneEmden type equations as a model for the study of the gravitational potential of these degenerate stars (whitedwarfs), which are nonlinear ordinary differential equations on the semi-infinite domain (see [11]). The Figure

Table 1. Comparison between Zhou’s method and SDMM.

Figure 2. The numerical results obtained using DTM with C = 0.2.

2 and table clearly show the high efficiency of DTM to solve nonlinear equations in comparison with other analytical methods equations (see [12,13]).

5. Acknowledgements

The authors gratefully acknowledge the support of the Universidad Tecnológica de Pereira and the group GEDNOL and would like to thank the referee for his valuable suggestions that improved the presentation of the paper.

NOTES

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1] S. Chandrasekhar, “Introduction to Study of Stellar Structure,” Dover, New York, 1967.
[2] J. K. Zhou, “Differential Transformation and Its Applications for Electrical Circuits,” Huazhong University of Science & Technology Press, Wuhan, 1986.
[3] S. J. Liao, “A New Analytic Algorithm of Lane-Emden Type Equations,” Advances Applied Mathematics, Vol. 1, No. 142, 2003, pp. 1-16.
[4] A. M. Wazwaz, “A New Algorithm for Solving Differential Equations of Lane-Emden Type,” Applied Mathematics and Computation, Vol. 118, No. 2-3, 2001, pp. 287310.
http://dx.doi.org/10.1016/S0096-3003(99)00223-4
[5] R. A. Gorder, “An Elegant Perturbation Solution for the Lane-Emden Equation of the Second Kind,” New Astronomy, Vol. 2, No. 16, 2011, pp. 65-67. http://dx.doi.org/10.1016/j.newast.2010.08.005
[6] J. I. Ramos, “Series Approach to the Lane-Emden Equation and Comparison with the Homotopy Perturbation Method,” Chaos Solitons Fractals, Vol. 38, No. 2, 2008, pp. 400-408.
http://dx.doi.org/10.1016/j.chaos.2006.11.018
[7] S. Iqbal and A. Javed, “Application of Optimal, Advances,” Applied Mathematics, Vol. 1, No. 42, 2005, pp. 2948.
[8] C. M. Khalique and P. Ntsime, “Exact Solutions of the Lane-Emden Type Equation,” New Astronomy, Vol. 7, No. 13, 2008, pp. 476-480. http://dx.doi.org/10.1016/j.newast.2008.01.002
[9] G. Hojjati and K. Parand, “An Efficient Computational Algorithm for Solving the Nonlinear Lane-Emden Type Equations,” International Journal of Mathematics and Computation, Vol. 4, No. 7, 2011, pp. 182-187.
[10] B. Batiha, “Numerical Solution of a Class of Singular Second-Order IVPs by Variational Iteration Method,” International Journal of Mathematical Analysis, Vol. 3, No. 40, 2009, pp. 1953-1968.
[11] H. Davis, “Introduction to Nonlinear Differential and Integral Equations,” Dover, New York, 1962.
[12] N. Kumar and R. Pandey, “Solution of the Lane-Emden Equation Using the Bernstein Operational Matrix of Integration,” ISRN Astronomy and Astrophysics, Vol. 2011, 2011, Article ID: 351747. http://dx.doi.org/10.5402/2011/351747
[13] S. Motsa and S. Shate, “New Analytic Solution to the Lane-Emden Equation of Index 2,” Mathematics Problems Engineering, Vol. 2012, 2012, Article ID: 614796.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.