Computer Science – Numerical Analysis
Scientific paper
Apr 1994
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=1994an....315..205r&link_type=abstract
Astronomische Nachrichten (ISSN 0004-6337), vol. 315, no. 3, p. 205-234
Computer Science
Numerical Analysis
6
Abundance, Mathematical Models, Nonlinear Evolution Equations, Numerical Analysis, Stellar Atmospheres, Stellar Evolution, Stellar Models, Stellar Structure, Algorithms, Boundary Value Problems, Euler-Lagrange Equation, Partial Differential Equations, Space-Time Functions
Scientific paper
We present a new numerical method for solving the system of partial differential equations describing the structure and evolution of a spherically symmetric star. As usual, we employ the transversal method of lines in order to split the equations into a coupled spatial and temporal part. The novel features of the algorithm are the following: Instead of using the Lagrangian picture we formulate the system of partial differential equations in the Eulerian picture. We reformulate the equations of stellar structure as a multipoint boundary-value problem By means of this reformulation the rather clumsy iterative matching procedure of stellar atmosphere and interior is avoided. The multipoint boundary-value problem is solved by the multiple shooting method. This approach not only ensures a high accuracy of the stellar models calculated at each time step but also allows the free boundaries inside the star due to different energy transport mechanisms to be located exactly. The time derivatives involved in the stellar-structure equations are discretized implicitly to second order accuracy. Moreover, at each time step, the chemical abundances are determined by using a sophisticated update procedure. In this way, a high accuracy is achieved with respect to the integration in time. The algorithm has turned out to be exceedingly reliable and numerically accurate. This is shown by the evolution of a 1 Solar Mass star up to the hydrogen-shell burning phase. In this example, the virial theorem, the law of mass conservation, the law of energy conservation is fulfilled to a hitherto unattainable degree of accuracy. Since the multiple shooting method, which is at the heart of our approach, is a perfect example of a parallel algorithm, the computational speed of the algorithm might be substantially improved provided easy-to-program, high-performance parallel computers with sufficiently many processors become available in the near future.
Bulirsch R.
Pfleiderer Jorg
Reiter Johann
No associations
LandOfFree
A multiple shooting approach for the numerical treatment of stellar structure and evolution does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with A multiple shooting approach for the numerical treatment of stellar structure and evolution, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and A multiple shooting approach for the numerical treatment of stellar structure and evolution will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-1255274