Пакет программ BPRQ предназначен для численного исследования нелинейных краевых задач на конечном отрезке [a,b] для системы из n обыкновенных дифференциальных уравнений. Векторная запись краевой задачи имеет вид:

Здесь  f (x, y, q) и g (y(a), y(b), q) - достаточно гладкие вектор-функции по совокупности аргументов, q – скалярный параметр, принадлежащий отрезку [qmin, qmax]. Предполагается, что в некоторой окрестности значения q = q0 краевая задача имеет решение y = y(t,q), дифференцируемое как по t, так и по параметру q.

  Типичной характерной особенностью нелинейных краевых задач является возможность существования областей изменения параметра, в которых одному и тому же значению параметра может соответствовать несколько решений. Возникновение множественности решений, больших градиентов и другие нелинейные эффекты учитываются в численном методе продолжения решения по параметру, который был разработан в Институте математики им. С. Л. Соболева СО РАН.

  Отметим, что в общем случае численное исследование решения нелинейной краевой задачи в зависимости от параметра, как правило, носит характер вычислительного эксперимента, для организации которого в каждом конкретном случае требуется «изобретать» индивидуальный подход к ряду неформализуемых проблем вычислительного характера, связанных с достаточно общей формулировкой рассматриваемой краевой задачи.

  Реализованный в пакете метод продолжения решения нашёл применение при численном исследовании достаточно большого числа краевых задач из различных приложений, таких как моделирование каталитических процессов, нелинейные колебания и т.д. В ряде случаев изучение множественности решений позволяет определить параметры гистерезисных явлений.

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

  Достоинства пакета BPRQ демонстрируются на моделях «Grain-Luss», «LK-Equat», «Platform», «Reactor-B», «Rele_D», «Rele_N», «Van_der_Pol», «Van_der_Pol1», «Zerno_P», моделирующих стационарные процессы каталитических реакций, в плёночной электромеханике, нелинейные колебания.

  1. Годунов С. К.
    Обыкновенные дифференциальные уравнения с постоянными коэффициентами. Издательство НГУ, 1994.-Т.1: Краевые задачи.- 264 с.
  2. Бахвалов Н. С.
    Численные методы, Москва, Наука, 1975, 632 с.
  3. Холодниок М., Клич А., Кубичек М., Марек М.
    Методы анализа нелинейных динамических моделей, Москва , Мир, 1991, 368 с.
  4. Завьялов Ю. С., Квасов Б. И., Мирошниченко В. Л.
    Методы сплайн-функций, Москва, Наука, 1980, 352 с.
  5. Фадеев С. И.
    Программа численного решения нелинейных краевых задач для систем обыкновенных дифференциальных уравнений с параметром // Вычислительные методы линейной алгебры. Новосибирск, Наука, Сибирское отделение, 1990, с. 104-200.
  6. Дятлов В. Л., Коняшкин В. В., Потапов Б. С., Фадеев С. И.
    Плёночная электромеханика.
    Новосибирск, Наука, Сибирское отделение,  1991. - 248 с.
  7. Когай В. В., Фадеев С. И.
    Применение продолжения по параметру на основе метода множественной стрельбы для численного исследования нелинейных краевых задач // Сибирский журнал индустриальной математики, т. 4, N 1(7), (2001), c. 83-101.
  8. С. И. Фадеев, В. В. Когай
    Краевые задачи для систем обыкновенных; дифференциальных уравнений. Учебное пособие. Новосибирск. НГУ. 2012, 278 с.
  9. Иоффе И. И., Письмен Л. М.
    Инженерная химия гетерогенного катализа. Ленинград, Химия, 1972.- 462 с.
  10. Зельдович Ф. Б., Баренблатт Г. И., Либрович В. Б., Махвиладзе Г. М.
    Математическая теория горения и взрыва, - Москва, Наука,1980. -  478 с.