All issues
- 2024 Vol. 16
- 2023 Vol. 15
- 2022 Vol. 14
- 2021 Vol. 13
- 2020 Vol. 12
- 2019 Vol. 11
- 2018 Vol. 10
- 2017 Vol. 9
- 2016 Vol. 8
- 2015 Vol. 7
- 2014 Vol. 6
- 2013 Vol. 5
- 2012 Vol. 4
- 2011 Vol. 3
- 2010 Vol. 2
- 2009 Vol. 1
-
О построении линейно неявных схем, LN-эквивалентных неявным методам Рунге–Кутты
Компьютерные исследования и моделирование, 2012, т. 4, № 3, с. 483-496В работе предложен новый класс безитерационных схем (явно-неявных), который позволяет получать методы, повторяющие на линейных неавтономных задачах свойства лучших неявных жестко-точных методов Рунге–Кутты [Хайрер, Ваннер,1999] – RadauIIA и LobattoIIIC. Для этого используется понятие LN-эквивалентности методов [Ширков, 2012]. С использованием среды аналитических вычислений получены уравнения порядка и затухания таких методов и найдены коэффициенты некоторых схем до 3-го порядка включительно. Проводится численное исследование новых методов на классических тестах, применяемых для проверки схем, разрабатываемых для жестких систем.
Ключевые слова: жесткие системы обыкновенных дифференциальных уравнений, жестко- точные методы Рунге–Кутты и схемы Розенброка, устойчивость и L-эквивалентность.
Constructing of linearly implicit schemes which are LN-equivalent to implicit Runge–Kutta methods
Computer Research and Modeling, 2012, v. 4, no. 3, pp. 483-496Views (last year): 2. Citations: 2 (RSCI).New family of linearly implicit schemes are presented. This family allows to obtain methods which are equivalent to stiffly accurate implicit Runge–Kutta schemes (such as RadauIIA and LobattoIIIC) on nonautonomous linear problems. Notion of LN-equivalence of schemes is introduced. Order conditions and stability conditions of such methods are obtained with the use of media for computer symbolic calculations. Some examples of new schemes have been constructed. Numerical studying of new method have been done with the use of classical tests for stiff problems.
-
Численно-аналитическое интегрирование уравнений свободного движения тяжелой точки вблизи звукового пика показателя степенного сопротивления
Компьютерные исследования и моделирование, 2013, т. 5, № 5, с. 785-798Показано, что для различных баллистических профилей во всем скоростном пространстве сила сопротивления изменяется со скоростью V по закону R(V)=Mg·w(V/WT)n(V), где WT — близкая к звуковой пороговая скорость, w=R(WT), n(V) — значение показателя в кусочно-степенной формуле. Методом, базирующимся на преобразованиях Лежандра, найдена отражающая пик n(V) поправка к невозмущенной резольвентной функции f(b)=abb'', a(b) — подкасательная к траектории, b=tgθ — ее наклон.
Ключевые слова: баллистический профиль, сопротивление, скорость, показатель степени, звуковой пик, преобразование Лежандра, резольвентная функция, поправка теории возмущений.
Numerical-analytical integrating the equations of a point mass projectile motion at the velocities close to sonic peak of air drag exponent
Computer Research and Modeling, 2013, v. 5, no. 5, pp. 785-798It is shown that the relative air drag force for many different ballistic profiles obeys the law as follows R(V)=Mg·w(V/WT)n(V) with V being the velocity, WT — some threshold velocity close to that of sound, w equals to R(WT) and n(V) is the exponent in broken power Gȃvre formula. Using the Legendre transformation and in frames of perturbation approach received was the expression for addition δabb''(b) to resolvent function abb''(b), where a(b) is an intercept and b=tgθ, θ — inclination angle.
-
Маршевый алгоритм решения задачи переноса излучения методом коротких характеристик
Компьютерные исследования и моделирование, 2014, т. 6, № 2, с. 203-215В работе изложена процедура построения численных решений для задачи переноса излучения. В этом подходе численное решение строится последовательно от границы области вдоль направления распространения излучения. Проведено тестирование алгоритма задаче распространения излучения нагретого шара.
Space-marching algorithm for solving radiative transfer problem based on short-characteristics method
Computer Research and Modeling, 2014, v. 6, no. 2, pp. 203-215Views (last year): 10. Citations: 3 (RSCI).A procedure of approximate solving of the radiation transfer problem is presented. The approximated solution is being built successively from the domain border along the direction of radiation propagation. The algorithm was tested for model problem of hot ball radiation.
-
Алгоритм численного интегрирования потенциально-потоковых уравнений в сосредоточенных параметрах с контролем корректности приближенного решения
Компьютерные исследования и моделирование, 2014, т. 6, № 4, с. 479-493Данная работа посвящена разработке алгоритма численного интегрирования системы дифференциальных уравнений потенциально-потокового метода моделирования неравновесных процессов. Этот метод был разработан автором в опубликованных им ранее работах. В настоящей работе рассмотрение ограничивается системами с сосредоточенными параметрами. Также ранее была разработана автором методика анализа корректности приближенного решения системы потенциально-потоковых уравнений для систем в сосредоточенных параметрах. Целью настоящей статьи является объединение этой методики с современными численными методами интегрирования систем обыкновенных дифференциальных уравнений и разработка методики численного интегрирования систем уравнений потенциально-потокового метода, позволяющей гарантировать корректность приближенного решения.
Ключевые слова: потенциально-потоковый метод, уравнения потенциально-потокового метода, численное интегрирование уравнений, анализ корректности приближенного решения.
Numerical integration algorithm potentially-streaming equations in lumped parameters to control the correctness of the approximate solution
Computer Research and Modeling, 2014, v. 6, no. 4, pp. 479-493Views (last year): 4. Citations: 3 (RSCI).This work is devoted to development of an algorithm for numerical integration of differential equations potentially-streaming method simulation of non-equilibrium processes. This method was developed by the author in his earlier published works. In this paper, consideration is limited to systems with lumped parameters. Also previously developed method for analyzing the correctness of the author of the approximate solution of the system potentially-streaming equations for systems in lumped parameters. The purpose of this article is to combine this technique with modern numerical methods for integrating systems of ordinary differential equations and the development of methods of numerical integration of systems of equations potentially-streaming method that allows to guarantee the correctness of the approximate solution.
-
Неявный итерационный полинейный рекуррентный метод в применении к решению задач динамики несжимаемой вязкой жидкости
Компьютерные исследования и моделирование, 2015, т. 7, № 1, с. 35-50В работе рассматриваются результаты применения неявного итерационного полинейного рекуррентного метода решения систем разностных эллиптических уравнений, возникающих при численном моделировании динамики несжимаемой вязкой жидкости. Исследование проводится на примере решения задачи о стационарном течении в плоской каверне с подвижной крышкой, сформулированной в естественных переменных ($u, \,v, \,p$) при больших значениях чисел Re (до 20 000) и сеточных разрешений (до 2049×2049). Демонстрируется высокая эффективность метода при расчете полей поправки давления. Анализируются проблемы решения задачи при больших числах Re.
The implicit line-by-line recurrence method in application to the solution of problems of incompressible viscous fluid dynamics
Computer Research and Modeling, 2015, v. 7, no. 1, pp. 35-50Views (last year): 3. Citations: 3 (RSCI).In the paper the results of applying the implicit line-by-line recurrence method for solving of systems of elliptic difference equations, arising, in particular, at numerical simulation of dynamics of incompressible viscous fluid are considered. Research is conducted on the example of the problem about a steady-state two-dimensional lid-driven cavity flow formulated in primitive variables ($u,\, v,\, p$) for large Re (up to 20 000) and grids (up to 2049×2049). High efficiency of the method at calculation of a pressure correction fields is demonstrated. The difficulties of constructing a solution of the problem for large Rе are analyzed.
-
Классификация динамических режимов переключения намагниченности в трехслойной ферромагнитной структуре в зависимости от спин-поляризованного тока инжекции и внешнего магнитного поля. I. Продольная анизотропия
Компьютерные исследования и моделирование, 2016, т. 8, № 4, с. 605-620В приближении однородной намагниченности построена математическая модель ячейки памяти MRAM c осью анизотропии, расположенной в плоскости запоминающего ферромагнитного слоя ячейки и ориентированной параллельно ее краю (продольная анизотропия). Модель базируется на уравнении Ландау–Лифшица–Гильберта с токовым членом в форме Слончевского–Берже. Выведена система обыкновенных дифференциальных уравнений в нормальном виде, описывающая динамику намагниченности в трехслойной вентильной структуре Co/Cu/Co в зависимости от величины тока инжекции и внешнего магнитного поля, параллельного оси анизотропии магнитных слоев. Показано, что при любых токах и полях система имеет два основных состояния равновесия, расположенных на оси, совпадающей с осью анизотропии. Проведен анализ устойчивости этих состояний равновесия. Выписаны уравнения для определения дополнительных состояний равновесия. Показано, что в зависимости от величины внешнего магнитного поля и тока инжекции система может иметь всего два, четыре и шесть симметричных относительно оси анизотропии положений равновесия. Построены бифуркационные диаграммы, характеризующие основные типы динамики вектора намагниченности свободного слоя. Проведена классификация фазовых портретов на единичной сфере в зависимости от управляющих параметров (тока и поля). Изучены особенности динамики вектора намагниченности в каждой из характерных областей бифуркационной диаграммы и численно построены траектории переключения. Для построения траекторий использовался метод Рунге–Кутты. Найдены параметры, при которых существуют неустойчивые и устойчивые предельные циклы. Установлено, что неустойчивые предельные циклы существуют вокруг основного устойчивого равновесия на оси, совпадающей с осью анизотропии, а устойчивые циклы — вокруг неустойчивых дополнительных равновесий. Граница области существования устойчивых предельных циклов рассчитана численно. Обнаружены новые типы динамики под влиянием внешнего магнитного поля и спин-поляризованного тока инжекции: случайное и неполное переключение намагниченности. Аналитически определены значения пороговых токов переключения в зависимости от внешнего магнитного поля. Численно выполнены оценки времени переключения в зависимости от величин управляющих параметров.
Ключевые слова: память MRAM, одноосная анизотропия, намагниченность, свободный слой, закрепленный слой, уравнение Ландау–Лифшица–Гильберта, переключение намагниченности.
Classification of dynamical switching regimes in a three-layered ferromagnetic nanopillar governed by spin-polarized injection current and external magnetic field. I. Longitudinal anisotropy
Computer Research and Modeling, 2016, v. 8, no. 4, pp. 605-620Views (last year): 2. Citations: 6 (RSCI).The mathematical model of the magnetic memory cell MRAM with the in-plane anisotropy axis parallel to the edge of a free ferromagnetic layer (longitudinal anisotropy) has been constructed using approximation of uniform magnetization. The model is based on the Landau–Lifshits–Gilbert equation with the injection-current term in the Sloncžewski–Berger form. The set of ordinary differential equations for magnetization dynamics in a three-layered Co/Cu/Cu valve under the control of external magnetic field and spin-polarized current has been derived in the normal coordinate form. It was shown that the set of equations has two main stationary points on the anisotropy axis at any values of field and current. The stationary analysis of them has been performed. The algebraic equations for determination of additional stationary points have been derived. It has been shown that, depending on the field and current magnitude, the set of equations can have altogether two, four, or six stationary points symmetric in pairs relatively the anisotropy axis. The bifurcation diagrams for all the points have been constructed. The classification of the corresponding phase portraits has been performed. The typical trajectories were calculated numerically using Runge–Kutta method. The regions, where stable and unstable limit cycles exist, have been determined. It was found that the unstable limit cycles exist around the main stable equilibrium point on the axis that coincides with the anisotropy one, whereas the stable cycles surround the unstable additional points of equilibrium. The area of their existence was determined numerically. The new types of dynamics, such as accidental switching and non-complete switching, have been found. The threshold values of switching current and field have been obtained analytically. The estimations of switching times have been performed numerically.
Indexed in Scopus
Full-text version of the journal is also available on the web site of the scientific electronic library eLIBRARY.RU
The journal is included in the Russian Science Citation Index
The journal is included in the RSCI
International Interdisciplinary Conference "Mathematics. Computing. Education"