259
Views
8
CrossRef citations to date
0
Altmetric
Research Articles

Execution of novel explicit RKARMS(4,4) technique in determining initial configurations of extra-solar protoplanets formed by disk instabilityFootnote

&
Pages 1-8 | Received 03 Jul 2014, Accepted 22 Nov 2015, Published online: 08 May 2019

Abstract

Implementation of a novel embedded Runge–Kutta fourth order four stage arithmetic root mean square technique to determine initial configurations of extra-solar protoplanets formed by gravitational instability is the main goal of this present paper. A general mathematical framework for the introduced numerical technique is described in addition to error estimation description. It is noticed that the numerical outputs through the employed novel RKARMS(4,4) method are found to be more effective and efficient in comparison with the results obtained by the classical Runge–Kutta technique.

1 Introduction

From literature review it is noticed that the ever-increasing advances in high performance computer technology have enabled several researchers towards science and engineering to employ novel numerical techniques to simulate physical phenomena. Intensive techniques are frequently required for the solution of real time practical problems and they often need the systematic application of a range of elementary techniques. In the development of new numerical methods, simplifications required to be made to progress towards an optimal solution. As a result, numerical algorithms do not usually give the exact answer to a given problem, or they can only tend towards a solution getting closer and closer with each iteration. Numerical techniques exhibit certain computational characteristics during their real time implementation. It is significant to consider these characteristics while selecting a specific technique for implementation. The characteristics which are critical to the success of implementation are accuracy, rate of convergence, numerical stability and efficiency. Numerical algorithms must review the factors such as, determination of the correctness of various steps, reduction of the number of steps, if necessary, and increase in the speed of solving the problem, respectively.

To solve many different problems under signal processing, communication, electronic and transistor circuits, Runge–Kutta (RK) method is being applied to obtain the required solution (CitationAlexander and Coyle, 1990). CitationShampine and Gordon (1975) discussed the normal order of a RK algorithm having the approximate number of leading terms of an infinite Taylor series, which calculates the trajectory of a moving point. CitationYaakub and Evans (1993) presented a new fourth order RK method based on the root mean formula for solving initial value problems (IVPs) in numerical studies. A new embedded fourth order RK method which is actually two different RK methods but of the same order p = 4 has been introduced by CitationEvans and Yaakub (1995). CitationBader (1987, Citation1998) introduced the RK–Butcher algorithm for finding the truncation error estimates, intrinsic accuracies and the early detection of stiffness in coupled differential equations arising in theoretical chemistry problems. CitationYaakub and Evans (1999) introduced a new fourth order RK technique for IVPs with error control. CitationButcher (1987, Citation1990, Citation2003) derived the best RK pair along with an error estimates and by all statistical measures it appeared as RK–Butcher algorithms. In order to overcome step-size constraint imposed by numerical stability, many new techniques have been developed in recent past. To confirm this, recently CitationPonalagusamy and Senthilkumar (2009) proposed a novel fourth order embedded RKARMS(4,4) technique based on RK arithmetic mean and root mean square with error control in detail to solve the real time application problems efficiently in image processing under CNN model. A detailed illustration related to the local truncation error (LTE), the global truncation error (GTE), error estimates and control for fourth order and four stage RK numerical algorithms is eventually addressed by CitationSenthilkumar (2009).

The formation of planetary systems has been a topic of interest to the mankind ever since the dawn of civilization. However, scientific theories for the formation of the system largely date from CitationDescartes (1644) when he proposed his vortex theory in this regard. Since that time many theories have been advanced. In most cases these theories were primarily speculative because of the lack of observational characteristics of the system. Fortunately, for the theorists of today, there are some convenient observational constrains of the system. The two end mechanisms, namely core accretion and disk instability, in principle, can form gas giant protoplanets. Though the core accretion mechanism has been, so far, adopted as the main theory of planetary formation both in our solar system and elsewhere, it fails to explain properly the recently discovered extrasolar protoplanets by direct imaging (see CitationDodson-Robinson et al., 2009). With this difficulty encountered by the core accretion models, the disk instability model, once in vague, has been reformulated with fragmentation from massive protoplanetary disks and has been advanced through the investigations of many authors (e.g., CitationBoss, 1997; CitationMayer et al., 2002, Citation2004; CitationBoley et al., 2010; CitationCha and Nayakshin, 2011). But this model is also criticized by some investigations with the argument that disk instabilities are unable to lead to the formation of self-gravitating dense clumps (CitationPickett et al., 2000; CitationCai et al., 2006; CitationBoley et al., 2007). Although some questions arise as to whether stable protoplanets could be formed or not by disk instability, the idea is believed to be a promising route to the rapid formation of giant planets in our solar system and elsewhere (CitationBoss, 2007). Unfortunately, the initial structures of the protoplanets formed via gravitational instability are still unknown and different numerical models can be found to report different configurations (CitationHelled and Schubert, 2008; CitationHelled and Bodenheimer, 2011).

It is pertinent to point out here that depending upon opacities of grains, present in protoplanets, as well as on initial conditions different investigators assumed different heat transports at different regions of protoplanets at different stages of their evolution. CitationDeCampli and Cameron (1979), in their investigation, assumed initial protoplanets to be fully convective with a thin outer radiative zone as a consequence of higher opacity and much work has since then been devoted to the evolution of planetary system including our own from such types of initial protoplanets (e.g., CitationBodenheimer et al., 1980; CitationWuchterl et al., 2000; CitationHelled et al., 2008; CitationHelled and Bodenheimer, 2011). It is well-known that depending on obeying the law L/4πR2=(surface opacity)-1, where L represents the luminosity and R is the protoplanetary radius, or on slow contraction, initial protoplanets may be fully convective (see CitationDeCampli and Cameron, 1979), which is consistent with CitationHelled et al. (2005). To investigate planetary evolution from such types of protoplanets, a series of studies were conducted by CitationPaul et al. (2008, Citation2012, Citation2013) and the obtained results were found to be in good agreement with the estimates by other investigations (see e.g., CitationHelled and Schubert, 2008; CitationHelled et al., 2008). However, recently CitationBoss (1998, Citation2002, Citation2007) in his investigations assumed the protoplanets to be in radiative equilibrium, which is consistent with earlier investigation by CitationBodenheimer (1974) who calculated a completely radiative Jovian mass structure assuming a constant grain opacity of 0.14 cm2 g−1. It is of interest to note here that in the case of radiative heat transfer, conduction is also taken part (CitationBöhm-Vitense, 1997). Based on the idea, CitationPaul et al. (2008) investigated initial structure of a Jovian mass protoplanet, which was further extended by CitationPaul and Bhattacharjee (2013) for investigating initial structures of extra solar protoplanets and the obtained results were found to be consistent with the results reported in some studies with rigorous treatment of the problem (see CitationPaul et al., 2013).

In this communication we intend to reinvestigate the model of CitationPaul and Bhattacharjee (2013) assuming heat transport following them to be conductive-radiative by a novel explicit RKARMS(4,4) method in order to test its validity and efficiency and to see how our computed results compare the estimates obtained with other investigations.

The rest of the article is structured as follows. The theoretical foundation of the problem in addition to boundary conditions is presented in Section 2. Numerical technique is adopted in Section 3. A brief description of the explicit RKARMS(4,4) technique along with local truncation error and error control is addressed in Section 4 and in its subsection. In Section 5, discussion of the obtained results as well as conclusion is presented.

2 Theoretical foundation

As in CitationPaul et al. (2008)Citationand Paul and Bhattacharjee (2013), the structure of a protoplanet, assuming the heat transport to be conductive-radiative, is given by the following set of equations:

The equation of hydrostatic equilibrium,(1) dP(r)dr=-GM(r)r2ρ(r).(1)

The equation of conservation of mass,(2) dM(r)dr=4πr2ρ(r).(2)

The equation of conductive-radiative heat flux,(3) 8σH3×10-24T3(r)ρ(r)+ηdT(r)dr=-CR4πRGM2(r)r3.(3)

The gas law,(4) P(r)=kμHρ(r)T(r).(4)

with the following boundary conditions(5) T(r)=0,P(r)=0atr=RsurfaceM(r)=Matr=R(surface)M(r)=0atr=0center.(5)

3 Numerical approach in structure determination

For the solution of structure equations, we have nondimensionalized them in addition to the boundary conditions following CitationPaul and Bhattacharjee (2013), which can be given by(6) dpdy=pqt(1-y)2,(6) (7) dqdy=-p(1-y)2t,(7) and(8) dtdy=CRcpq2(1-y)3(at4+bp),(8) wherea=8σH3×10-24μHGMkR3,b=Mη4πR3,andc=M2k16π2R5μH,as by means of the above transformations, ρ is reduced to the form(9) ρ=M4πR3pt,(9) where the boundary conditions are given by(10) t(y)=0,p(y)=0aty=0q(y)=1aty=0q(y)=0aty=1.(10)

Now, analytic solution of the system specified by Eqs. (Equation(6)Equation(8)) with the boundary conditions specified by Eq. (Equation10), as they stand, is impossible (CitationPaul et al., 2013), resort to be taken to numerical technique. But because of the existence of vanishing denominators in the basic equations, the integration cannot be started right from either of the boundaries, and hence the boundary conditions should be developed. In our investigation, we used developed surface boundary conditions, which are available in CitationPaul and Bhattacharjee (2013). With the developed boundary conditions, inserting the values of the required parameters involved, we have solved Eqs. (Equation(6)Equation(8)) numerically by the newly introduced explicit RKARMS(4,4) method from y=0.01 downwards to the point 0.99 to get the distributions of p,q, and t. The values of the required parameters of the present study are similar to those used in the study of CitationPaul and Bhattacharjee (2013). With the distributions of p and t, the density distribution is then easily obtained by Eq. (Equation9). The structures of the protoplanets are found to be dependent on a parameter CR. The best values of CR for the prescribed protoplanetary masses 0.3MJ,1MJ,3MJ,5MJ,7MJ and 10MJ satisfying the third condition of Eq. (Equation10) can be found to be 0.026, 0.2, 1.27, 2.43, 4.03 and 8.4, respectively, which can be found to be the same with the ones found in CitationPaul and Bhattacharjee (2013). The results of our calculation are shown in diagrammatic forms through .

Figure 1 Initial temperature–pressure profiles of some protoplanets.

Figure 2 Density distribution inside some initial protoplanets.

Figure 3 Error estimation of the proposed method through p, q and t with references to the protoplanets with masses 1MJ and 10MJ; (a) for 1MJ and (b) for 10MJ.

4 A description about explicit RKARMS(4,4) numerical technique

The sstage RK(4,4) technique for solving the IVP y=f(x,y(x)),x0 x xn subject to y(x0)=y0 can be given by(11) yn+1=yn+hi=1sbiki,(11) whereki=f(xn+cihyn+hi=1saijki),ci=j=1saij,i=123,,swith s dimensional vectors c and b and the s×s matrix A(aij). A general s-stage RK pair can be written in an array form as

The symbols C,A and bT have order s and that C,A and b̂T have order (s+1). Using the second method, the value of y at x=xn+1 can be expressed as(12) ŷn+1=yn+hi=1sb̂iki,(12) whereas the same for the first method is expressed by Eq. (Equation11).

From the embedded form, the LTE may be computed from the formula: LTE=yn+1-ŷn+1. It is of interest to note here that LTE leads to control step size.

The four stage method with the Butcher array form is written as follows:

The well-known fourth order RK arithmetic mean (RKAM(4,4)) method can be written in the Butcher array form as follows (see CitationSenthilkumar and Paul, 2012): (13) yn+1=yn+h3k1+k22+k2+k32+k3+k42,(13) wherek1=f(xn,yn),k2=fxn+h2,yn+hk12,k3=fxn+h2,yn+hk22,k4=f(xn+h,yn+hk3).

The fourth order RK method with Butcher array can also be written in the modified form as (see CitationPonalagusamy and Senthilkumar, 2009)

The fourth order RK root mean square (RKRMS(4,4)) method due to CitationYaakub and Evans (1993) is given by(14) yn+1=yn+h3k12+k222+k22+k322+k32+k422,(14) wherek1=f(xn,yn),k2=fxn+12h,yn+12hk1,k3=fxn+12h,yn+116hk1+716hk2,k4=fxn+h,yn+18hk1-1756hk2+3328hk3.

It is well-known that combination of RKAM(4,4) and RKRMS(4,4) (Eqs. (Equation(15)Equation(14))) leads to give a new formation of RKARMS(4,4), and is formulated by CitationPonalagusamy and Senthilkumar (2009) ask1=f(xn,yn)=k1,k2=fxn+h2,yn+hk12=k2,k3=fxn+h2,yn+hk22,k4=f(xn+h,yn+hk3),k3=fxn+12h,yn+116hk1+716hk2=k3,k4=fxn+h,yn+18hk1-1756hk2+3328hk3=k4,(15) yn+1=yn+h3k1+k22+k2+k32+k3+k42,(15) (16) yn+1=yn+h3k12+k222+k22+k322+k32+k422.(16)

The embedded RKARMS(4,4) method is expressed as(17) Hence,bT=yn+1AM=yn+h3k1+k22+k2+k32+k3+k42,(17) (18) b̂T=yn+1RMS=yn+h3k12+k222+k22+k322+k32+k422,(18) and the estimation of the LTE, ET=bT-b̂T. In the RKARMS(4,4) method, four stages are required to obtain the solution, which share the same set of vectors k1 and k2 using bT and b̂T approximately, but k3 and k4 use bT while k3 and k4 use b̂T.

4.1 Derivation and error estimation for explicit RKARMS(4,4) method

According to CitationLotkin (1951)Citation, Ralston (1957)Citationand Lambert (1973, Citation1980), the error estimate for fourth order RK schemes is given ψ(xn,yn:h)(73/720)ML4, where L and M are positive constants. From Eqs. (Equation(17)Equation(18)), we obtain an estimate of the LTE for the RKARMS(4,4) method as LTE=yn+1-yn+1, which may be used to control step size. The LTE for well-known RKAM(4,4) method is(19) yn+1AM=yn+LTEAM,(19) and the LTE for RKRMS(4,4) method is(20) yn+1RMS=yn+LTERMS,(20) where yn+1AMandyn+1RMS are numerical approximations at xn+1 obtained by RKAM(4,4) and RKRMS(4,4), respectively, and LTEAM and LTERMS are the LTEs for RKAM(4,4) and RKRMS(4,4), respectively. The difference between yn+1AMandyn+1RMS at xn+1 gives an error estimate as(21) yn+1AM-yn+1RMS=LTEAM-LTERMS.(21)

The LTE for RKAM(4,4) method is given by(22) LTEAM=h52880-24ffy4+f4fyyyy+2f3fyfyyy-6f3fyy2+36f2fy2fyy,(22) whereas the LTE for RKRMS(4,4) method can be set to the form(23) LTERMS=h5184320-429ffy4-64f4fyyyy-48f3fyfyyy-96f3fyy2-2454f2fy2fyy.(23)

The absolute difference between LTEAMandLTERMS is given by(24) LTEAM-LTERMS=h51843201107ffy4+128f4fyyyy+176f3fyfyyy+288f3fyy2+4758f2fy2fyy.(24)

As in Eq. (Equation20), substituting f,fy,fyy, etc. in Eq. (Equation24), it can be written as(25) LTEAM-LTERMS6457184320P4Qh5,(25) where P and Q are positive constants. If we let TOL=5.00×10-5, then by setting LTEAM-LTERMSTOL, the error control and step size selection can be determined by Eq. (Equation25) as(26) 6457184320P4Qh5<TOLorh<28.54764×TOLP4Q1/5.(26)

It is pertinent to point out that in the explicit RKARMS(4,4) method with error control program, we choose error estimation as the difference between the results obtained by RKAM(4,4) and RKRMS(4,4) methods. From Eq. (Equation25), the error estimation (ERREST) is expressed as (see )(27) ERREST=yAM-yRMS×6457184320.(27)

Table 1 Comparison of LTE, GTE, and error estimation for RKARMS(4,4) method.

5 Discussion on results and conclusion

We have analyzed initial configurations of protoplanets formed via disk instability in the mass range 0.3–10 Jovian masses by a newly proposed explicit RKARMS(4,4) method under approximate zero boundary conditions. Each of the protoplanets has been assumed to be a sphere of solar composition, which is in a steady state of quasi–static equilibrium where ideal gas law holds well, and the energy equation assumes the conduction–radiation heat transport. depicts initial temperature–pressure profiles of the protoplanets with assumed masses. It can be shown from the figure that after a point little depth from the surface down to the core region, the temperature increases with the increasing mass of the protoplanets, whereas the pressure also increases with their increasing mass except for the protoplanets with masses 1MJ and 10MJ. The results of our calculation agree fairly well with the estimates obtained in the investigations of CitationHelled and Schubert (2008)Citation, Paul and Bhattacharjee (2013)Citation, and Paul et al. (2013). shows initial density distribution inside the assumed protoplanets. The figure shows that mater is not distributed uniformly in the atmosphere, and there may have variation in parameters due to gravitational stratification. This is as to be expected for initial unsegregated protoplanets. It can be observed from the figure that the central density of the protoplanet with mass 1MJ is higher than the central condensations of the other considered protoplanets and the protoplanet with mass 3MJ is found to be rarer concerning both center and surface among all the protoplanets with assumed masses.

Density distribution obtained in the study can be found to be comparable with the ones obtained in the study of CitationPaul et al. (2013). But CitationHelled and Schubert (2008) showed that the surface density of such protoplanets decreases with their decreasing mass and central density increases with their increasing mass. It is pertinent to point out here that, in reality, not a single protoplanet formed by disk instability exists in the literature with its definite structures (CitationHelled and Bodenheimer, 2011; CitationPaul and Bhattacharjee, 2013). However, the system possesses a unique solution suggesting that disk instability is a reasonable hypothesis in planetary formation. The results of our calculation may be important in the study of evolution of extrasolar giant planets. A direct comparison of the results employing the proposed RKARMS(4,4) method is made with the results obtained by the classical Runge–Kutta 4th order (RK(4,4)) method. Only the results for 1MJ are presented (see ) because to void space consumption. From the table it is found that the results by the explicit RKARMS(4,4) method are comparable with the results obtained by the RK(4,4) method, which can be seen to be true for all the protoplanetary masses considered. We have evaluated error estimation (ERREST) for the proposed method through the results of the investigation. The results for 1MJ and 10MJ protoplanets are presented in for sake of brevity. It is inferred from that the error estimations are quite reasonable. We have tested our results for varying end points with all the possible step sizes, for which the RKARMS(4,4) method is valid. The results are found to be insensitive to the choice of the end points. In order to compare the computational efficiency of the RKARMS(4,4) method with that of the RK(4,4) method, both the codes were run on the same computer with step size 0.0001. The total computational time was found to be less for the RKARMS(4,4) method (5.978359 s) in comparison with that for the RK(4,4) method (6.065116 s). Thus the RKARMS(4,4) method is found to be more optimal in solving structure equations of protoplanets in comparison with the classical RK(4,4) method with respect to the central processing unit (CPU) time and accuracy.

Table 2 Comparative distribution of thermodynamic variables inside a 1 Jupiter mass protoplanet.

In our calculation, to make the work simple, the protoplanets are assumed to be spheres of gas and dust where ideal gas law holds well. But the Clapeyron equation of state (ideal gas law) is appropriate only for the gases with no high pressure. Also, in our calculations, we have neglected radiation effect from the parent star. Furthermore, the disturbances from the parent star and the mutual attraction among the protoplanets have not been considered. But future work will be concentrated on the evolution of extrasolar planets formed via disk instability including the factors mentioned above using an appropriate equation of state, where our intention is to implement parallel numerical algorithms that employ large number of processors. The processors perform various tasks independently and simultaneously, thereby, improving the speed of execution of complex programs dramatically. Parallel computers match the speed of supercomputers at a fraction of the cost.

Gour Chandra Paul W as born in Talora, Bogra District, Bangladesh on 20th November 1974. He received his B.Sc. honours degree on Mathematics in 1994, M.Sc. degree on Applied Mathematics in 1995 and Ph.D. degree in Astrophysics in 2005 from the Department of Mathematics, University of Rajshahi, Bangladesh. He was a Post-Doctoral Fellow in the School of Mathematical Sciences, Universiti Sains Malaysia, Pulau Pinang, Malaysia during 2010–2011. Presently he is working as a professor in the Department of Mathematics, Rajshahi University, Bangladesh. He has published many good research papers in peer-reviewed international journals with high impact factor. He is also reviewer and referee for many scientific international journals. His current research interest lies in evolution of Planetary Systems, Numerical Methods, Tide and Storm Surges Modeling and Simulation.

Sukumar Senthilkumar was born in Neyveli Township, Cuddalore District, Tamilnadu, India on 18th July 1974. He received his B.Sc Degree in Mathematics from Madras University in 1994, M.Sc Degree in Mathematics from Bharathidasan University in 1996, M.Phil Degree in Mathematics from Bharathidasan University in 1999 and M.Phil Computer Science & Engineering from Bharathiar University in 2000. Also, he received PGDCA and PGDCH in Computer Science and Applications and Computer Hardware from Bharathidasan University in 1996 and 1997 respectively. He has been awarded doctoral degree in the field of Mathematics and Computer Applications from National Institute of Technology, [An Institute of National Importance-formerly Regional Engineering College], Tiruchirappalli, Tamilnadu, India. He worked as a Post Doctoral Fellow at Universiti Sains Malaysia [An Apex University], School of Mathematical Sciences, Pulau Pinang, Malaysia. Also, he worked as a Post Doctoral Researcher at the Center for Advanced Image and Information Technology, School of Electronics and Information Engineering, Division of Computer Science and Engineering, Chon Buk National University, Jeonju, Chon Buk, South Korea. He worked as a Lecturer/Assistant Professor in the Department of Computer Science at Asan Memorial College of Arts and Science, Chennai, Tamilnadu, India. Currently, he is working as an Associate Professor in the School of Computing Science and Engineering at Vellore Institute of Technology-University, Vellore, Tamilnadu, India. He has published many good scientific research papers in international conferences and in peer reviewed/highly refereed international journals with high impact factor. He has made significant and outstanding contributions to various activities related to research work. He is also an editor, associate editor, editorial board member, advisory board member, reviewer and referee for many scientific international journals. His current research interests lies in advanced cellular neural networks, advanced digital image processing, advanced numerical analysis, methods and algorithms, advanced simulation, modeling and computing. Further, he is actively engaged in other related areas including advanced astronomy, advanced astrophysics and advanced oceanography etc.

Acknowledgments

We would like to thank both the anonymous referees for useful comments that improved the manuscript. The first author also wishes to thank Professor Shishir Kumer Bhattacharjee, Department of Mathematics, Rajshahi University, for discussions and helpful suggestions during the time this manuscript was being prepared. The second author would like to extend his sincere gratitude to Government of India, eternally for providing financial support during doctoral studies via Technical Quality Improvement Programme (TEQIP), Under Ministry of Human Resource Development, to National Institute of Technology [REC], Tiruchirappalli-620015, Tamil Nadu, India. URL: http://www.nitt.edu. Also, the second author wishes to thank Government of Malaysia for offering financial grants during post doctoral fellow scheme to Universiti Sains Malaysia, School of Mathematical Sciences, Pulau Pinang-11800, Penang, Malaysia. URL: http://www.usm.my; Further, the second author express his grateful thanks to Republic of Korea for offering financial grants through second stage of Brain Korea 21 programme for post doctoral researcher scheme to Center for Advanced Image and Information Technology, School of Electronics & Information Engineering, Chon Buk National University, 664-14, 1Ga, DeokJin-Dong, Jeonju, Chon Buk, 561-756, South Korea. URL: http://www.ei.jbnu.ac.kr. and URL: http://www.ei.bk21.jbnu.ac.kr.

Notes

Peer review under responsibility of National Research Institute of Astronomy and Geophysics.

References

  • R.K.AlexanderJ.J.CoyleRunge–Kutta methods and differential-algebraic systemsSIAM J. Numer. Anal.271990736752
  • M.BaderA comparative study of new truncation error estimates and intrinsic accuracies of some higher order Runge–Kutta algorithmsComput. Chem.111987121124
  • M.BaderA new technique for the early detection of stiffness in coupled differential equations and application to standard Runge–Kutta algorithmsTheor. Chem. Acc.991998215219
  • P.BodenheimerCalculations of the early evolution of JupiterIcarus231974319325
  • P.BodenheimerA.S.GrossmanW.M.DeCampliG.MarcyJ.B.PollackCalculations of the evolution of the giant planetsIcarus411980293308
  • E.Böhm-VitenseIntroduction to Stellar Astrophysics1997Cambridge University Press
  • A.C.BoleyT.W.HartquistR.H.DurisenS.MichaelThe internal energy for molecular hydrogen in gravitationally unstable protoplanetary disksAstrophys. J.6562007L89L92
  • A.C.BoleyT.HayfieldL.MayerR.H.DurisenClumps in the outer disk by disk instability: why they are initially gas giants and the legacy of disruptionIcarus2072010509516
  • A.P.BossGiant planet formation by gravitational instabilityScience276199718361839
  • A.P.BossFormation of extra solar giant planets: core accretion or disk instability?Earth Moon Planets8119981926
  • A.P.BossEvolution of the solar nebula. V. Disk instabilities with varied thermodynamicsAstrophys. J.5762002462472
  • A.P.BossTesting disk instability models for giant planet formationAstrophys. J.6612007L73L76
  • J.C.ButcherThe numerical analysis of ordinary differential equations: Runge–Kutta and general linear methods1987John Wiley & SonsChichester, UK
  • J.C.ButcherOn order reduction for Runge-Kutta methods applied to differential-algebraic systems and to stiff systems of ODEsSIAM J. Numer. Anal.271990447456
  • J.C.ButcherThe Numerical Analysis of Ordinary Differential Equations2003John Wiley & SonsChichester, UK
  • K.CaiR.H.DurisenS.MichaelA.C.BoleyA.C.MejiaM.K.PickettP.D’AlessioThe effects of metallicity and grain size on gravitational instabilities in protoplanetary disksAstrophys. J.6362006L149L152
  • S.-H.ChaS.NayakshinA numerical simulation of a super-earth core delivery from ∼100 AU to ∼8 AUMNRAS415201133193334
  • W.M.DeCampliA.G.W.CameronStructure and evolution of isolated giant gaseous protoplanetsIcarus381979367391
  • Descartes, R., 1644. Principia Philosphiae. Amsterdam.
  • S.E.Dodson-RobinsonD.VerasE.B.FordC.A.BeichmanThe formation mechanism of gas giants on wide orbitsAstrophys. J.70720097988
  • D.J.EvansA.R.YaakubA new Runge–Kutta RK (4,4) methodInt. J. Comput. Math.581995169187
  • R.HelledA.KovetzM.PodolakSettling of small grains in an extended protoplanetBull. Am. Astron. Soc.372005675
  • R.HelledG.SchubertCore formation in giant gaseous protoplanetsIcarus1982008156162
  • R.HelledM.PodolakA.KovetzGrain sedimentation in a giant gaseous protoplanetIcarus1952008863870
  • R.HelledP.BodenheimerThe effects of metallicity and grain growth and settling on the early evolution of gaseous protoplanetsIcarus2112011939947
  • J.D.LambertComputational Methods in Ordinary Differential Equations1973WileyNew York
  • J.D.LambertStiffness in computation techniques for ordinary differential equationsI.GladwellD.K.Sayers1980Academic PressNew York1946
  • M.LotkinOn the accuracy of RK methodsMTAC51951128132
  • L.MayerT.QuinnJ.WadsleyJ.StadelFormation of giant planets by fragmentation of protoplanetary disksScience298200217561759
  • L.MayerT.QuinnJ.WadsleyJ.StadelThe evolution of gravitationally unstable protoplanetary disks: fragmentation and possible giant planet formationAstrophys. J.609200410451064
  • G.C.PaulJ.N.PramanikS.K.BhattacharjeeStructure of initial protoplanetsInt. J. Mod. Phys. A23200828012808
  • G.C.PaulJ.N.PramanikS.K.BhattacharjeeGravitational settling time of solid grains in gaseous protoplanetsActa Astronaut.7620129598
  • G.C.PaulS.K.BhattacharjeeDistribution of thermodynamic variables inside extra–solar protoplanets formed via disk instabilityEgypt. J. Remote. Sensing Space Sci.1620131721
  • G.C.PaulM.M.RahmanD.KumarM.C.BarmanThe radius spectrum of solid grains settling in gaseous giant protoplanetsEarth Sci. Inform.62013137144
  • B.K.PickettP.CassenR.H.DurisenR.LinkThe effects of thermal energetics on three-dimensional hydrodynamic instabilities in massive protostellar disks. II. High-resolution and adiabatic evolutionsAstrophys. J.529200010341053
  • R.PonalagusamyS.SenthilkumarA new method of embedded fourth order with four stages to study raster CNN simulationInt. J. Autom. Comput.62009285294
  • Senthilkumar, S., 2009. New embedded Runge–Kutta fourth order four stage algorithms for raster and time-multiplexing cellular neural networks simulation. Ph.D Thesis. Department of Mathematics, National Institute of Technology [REC], Tiruchirappalli, Tamilnadu, India.
  • S.SenthilkumarG.C.PaulApplication of new RKAHeM(4,4) technique to analyze the structure of initial extrasolar giant protoplanetsEarth Sci. Inform.520122331
  • L.F.ShampineM.K.GordonComputer Solution of Ordinary Differential Equations – The Initial Value Problem1975Freeman W.H. & Co.San Francisco, USA
  • R.H.RalstonRunge–Kutta methods with minimum error bondsMath. Comput.161957431437
  • Yaakub, A.R., Evans, D.J., 1993. A new fourth order Runge–Kutta method based on the root mean square formula. Computer Studies Report No. 862. Loughborough University, United Kingdom.
  • A.R.YaakubD.J.EvansA fourth order Runge–Kutta RK(4,4) method with error controlInt. J. Comput. Math.711999383411
  • G.WuchterlT.GuillotJ.J.LissauerGiant planet formationV.ManningsA.P.BossS.S.RussellProtostars and Planets IV2000Univ. of Arizona PressTucson10811109