J. Astrophys. Astr. (2019) 40:27 © Indian Academy of Sciences https://doi.org/10.1007/s12036-019-9587-0
Numerical study of singular fractional Lane–Emden type equations arising in astrophysics
ABBAS SAADATMANDI∗ , AZAM GHASEMI-NASRABADY and ALI EFTEKHARI
Department of Applied Mathematics, University Kashan, Kashan 87317-53153, Iran. ∗Corresponding author. E-mail: [email protected]
MS received 20 December 2018; accepted 13 April 2019; published online 17 June 2019
Abstract. The well-known Lane–Emden equation plays an important role in describing some phenomena in mathematical physics and astrophysics. Recently, a new type of this equation with fractional order derivative in the Caputo sense has been introduced. In this paper, two computational schemes based on collocation method with operational matrices of orthonormal Bernstein polynomials are presented to obtain numerical approximate solutions of singular Lane–Emden equations of fractional order. Four illustrative examples are implemented in order to verify the efficiency and demonstrate solution accuracy.
Keywords. Fractional Lane–Emden type equations—orthonormal Bernstein polynomials—operational matrices—collocation method—Caputo derivative.
1. Introduction hydraulic pressure at a certain radius r. Due to the contribution of photons emitted from a cloud in radi- The Lane–Emden equation has been used to model ation, the total pressure due to the usual gas pressure many phenomena in mathematical physics, astrophysics and photon radiation is formulated as follows: and celestial mechanics such as the thermal behavior of 1 RT a spherical cloud of gas under mutual attraction of its P = ξ T 4 + , 3 V molecules, the theory of stellar structure, and the theory ξ, , of thermionic currents (Chandrasekhar 1967). This is a where the parameters T R and V stand for the radia- singular nonlinear ordinary differential equation (ODE) tion constant, the absolute temperature, the gas constant ( ) which was first introduced by U.S. astrophysicist, and the specific volume respectively. Let M r be the , ,φ Jonathan Homer Lane (Lane 1870) who was interested mass of a sphere of radius r and G g respectively in computing both the temperature and the density of denote the constant of gravitation, acceleration of grav- mass on the surface, and 37 years later was studied ity and gravitational potential of gas. By having the in more detail by Emden (1907). In astrophysics, this relation GM(r) dφ ODE is a dimensionless version of Poisson’s equation g = =− , for the gravitational potential of a simple stellar model r 2 dr (Momoniat & Harley 2006). We now briefly explain one can determine P and φ by meeting the following mathematical modeling of the thermal behavior of a three conditions: spherical cloud which leads to the classic Lane–Emden dP =−gρdr = ρdφ, (1) equation. See Chandrasekhar (1967) and Parand et al. 2 (2010) for a comprehensive coverage of the subject. ∇ φ =−4πGρ, (2) γ P = Kρ . (3) 1.1 Modeling of the thermal behavior of a spherical Here, ρ is the density of gas at a distance r from the cen- cloud of gas tre of the spherical cloud and the constants γ = 1+1/m and K are experimentally determined. The constant m For the sake of simplicity, we consider a spherical gas is called the polytropic index and is related to the ratio cloud as shown in Figure 1. Let P(r) denote the total of specific heats of gas comprising the star. Using the 27 Page 2 of 12 J. Astrophys. Astr. (2019) 40:27