Все выпуски

В работе рассматривается новая методика моделирования реакционно-диффузионных систем на основе систем обыкновенных дифференциальных уравнений. В отличие от специализированных численных методов, таких как метод прямых, новая методика позиционируется как чистая альтернатива на модельном уровне для уравнений в частных производных. По своему описанию новый метод во многом подобен методу конечных объемов, но в отличие от него для описания диффузии применяет статистические упрощения и принципы геометрической вероятности. Главными задачами данного подхода являются упрощение качественного анализа реакционно-диффузионных систем, а также повышение эффективности численной реализации модели. Первая задача успешно решается, так как для качественного анализа динамики модели на основе систем обыкновенных дифференциальных уравнений становится возможным использовать аппарат классической теории динамических систем. Вторая задача решается лишь отчасти, так как выигрыш при сохранении приемлемой точности для численной реализации будет существенным лишь для определенных, достаточно простых, начальных распределений молекул, а также для определенных коэффициентов диффузии. При этом для формирования критериев применимости на практике мы отдельно оцениваем погрешность моделирования с использованием данной новой методики.

We consider a new technique for modelling the reaction-diffusion systems based on systems of ordinary differential equations. In contrary to the specialized numerical methods such as straight line method, this new technique is positioned at model level as a full alternative for partial differential equations. The description of this new method is quite similar to the description of finite volume method, except that it uses statistical simplifications and principles of geometric probability to describe diffusion processes. The main goal of this approach is to simplify the qualitative analysis of reaction-diffusion systems and to increase the efficiency of the numerical implementation. The first task is successfully resolved because of the fact that for the qualitative analysis of model dynamics based on ordinary differential equations it is possible to use the apparatus of the classical theory of dynamical systems. The second task is solved only partially, because the gain in efficiency while maintaining acceptable accuracy for numerical implementation will be considerable only for certain simple initial distribution of molecules, as well as for certain diffusion coefficients. To determine the criteria for practical application of this technique we also estimate the model error in general.