On the efficiency and robustness of implicit methods in computational astrophysics

Astronomy and Astrophysics – Astronomy

Scientific paper

Rate now

  [ 0.00 ] – not rated yet Voters 0   Comments 0

Details

19

Scientific paper

Most of the numerical methods used in astrophysical fluid dynamics rely on explicit time-stepping schemes, whereas the higher robustness of implicit methods which constitute the core of modern computational fluid dynamics is rarely explored. In this paper, we survey some modern implicit solvers which are specially adapted to multi-dimensional problems and discuss their potential and range of application in comparison to common explicit methods. Special emphasis is put on the aspect of efficiency and robustness. Our reference set of equations are those corresponding to radiative magneto-hydrodynamics (MHD) modeling self-gravitating and partially and/or fully ionized flows. Explicit methods may be viewed as a very special class of so-called `defect-correction iterations' for solving an implicit discretization. Within this context one can design various implicit methods, ranging from weakly to fully implicit, which allow to follow evolutionary phases on much longer time scales than the dynamic one. We particularly present a new three-stages implicit numerical method for searching strongly time-dependent, quasi-stationary and steady-state solutions for the above-mentioned equations. Preconditioned Krylov-space and multilevel techniques are employed for enhancing the efficiency and robustness of the computation. The spatial discretization is on highly nonuniform tensor-product meshes and uses cartesian, cylinder or spherical coordinates depending on the geometrical structure of the problem. The accuracy is of second-order in space and time and can easily be increased without modifying the structure of the scheme. The algebraic solver consists of a pre-conditioned transpose-free Krylov iteration for the conservation equations and optimized multigrid algorithms for solving the Poisson equation for the gravitational potential and the transport-diffusion equation for the radiation density in flow regions with dynamically varying optical depths. The existing implementation of the proposed method employs axis-symmetry in order to reduce the problem to two dimensions.

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

On the efficiency and robustness of implicit methods in computational astrophysics 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 On the efficiency and robustness of implicit methods in computational astrophysics, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and On the efficiency and robustness of implicit methods in computational astrophysics will most certainly appreciate the feedback.

Rate now

     

Profile ID: LFWR-SCP-O-1257592

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