Numerical methods and programming
 RUS  ENG JOURNALS   PEOPLE   ORGANISATIONS   CONFERENCES   SEMINARS   VIDEO LIBRARY   PACKAGE AMSBIB
 General information Latest issue Archive Search papers Search references RSS Latest issue Current issues Archive issues What is RSS

 Num. Meth. Prog.: Year: Volume: Issue: Page: Find

 Num. Meth. Prog., 2016, Volume 17, Issue 4, Pages 365–379 (Mi vmp843)

On the correctness of numerical simulation of gravitational instability with the evolution of multiple gravitational collapses

V. N. Snytnikova, O. P. Stoyanovskayab

a Boreskov Institute of Catalysis SB RAS, Novosibirsk
b Novosibirsk State University

Abstract: We study the correctness of numerical models based on the SPH-method for nonstationary problems of gravitational gas dynamics with the development of gravitational instability, including with the formation of multiple collapses of a gas in a circumstellar disk. The original differential initial boundary value problems are ill-posed because of their instability with respect to variations of input data. It is shown that the following regularization is performed in the SPH-based numerical method of solving such problems: (1) if the solution to the original unstable problem exists on the entire time axis, then the condition of the method's stability with respect to small variations of input data is satisfied due to a bounded range of the variables; (2) if the solution of the original unstable problem exists only on a bounded time interval, as in the case of multiple collapses, then a stable numerical method can be developed on a class of functions bounded by a fixed constant chosen from physical considerations. On this class of bounded functions, the original problem becomes well-posed. A combination of the SPH method with a grid-based method of calculating the gravitational forces allows one to ensure the boundedness of numerical solutions. In order to clarify a meaning of the approximate numerical solutions obtained by computer simulations, it is necessary to use the integral functions weakly sensitive to the details of the numerical algorithm in use.

Keywords: well-posedness, numerical methods, gravitational gas dynamics, instability, stability with respect to the input data, stochastics, Smoothed-Particle Hydrodynamics (SPH), circumstellar disc, gravitational collapse, gas clump.

Full text: PDF file (1195 kB)
UDC: 523.4-52; 573.552

Citation: V. N. Snytnikov, O. P. Stoyanovskaya, “On the correctness of numerical simulation of gravitational instability with the evolution of multiple gravitational collapses”, Num. Meth. Prog., 17:4 (2016), 365–379

Citation in format AMSBIB
\Bibitem{SnySto16} \by V.~N.~Snytnikov, O.~P.~Stoyanovskaya \paper On the correctness of numerical simulation of gravitational instability with the evolution of multiple gravitational collapses \jour Num. Meth. Prog. \yr 2016 \vol 17 \issue 4 \pages 365--379 \mathnet{http://mi.mathnet.ru/vmp843}