A multiple shooting approach for the numerical treatment of stellar structure and evolution

Computer Science – Numerical Analysis

Scientific paper

Rate now

  [ 0.00 ] – not rated yet Voters 0   Comments 0

Details

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.

No associations

LandOfFree

Say what you really think

Search LandOfFree.com for scientists and scientific papers. Rate them and share your experience with other people.

Rating

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.

Rate now

     

Profile ID: LFWR-SCP-O-1255274

  Search
All data on this website is collected from public sources. Our data reflects the most accurate information available at the time of publication.