ВЫЧИСЛИТЕЛЬНЫЕ МЕТОДЫ И ПРОГРАММИРОВАНИЕ

Работа посвящена нелокальным краевым задачам для многомерного уравнения параболического типа с переменными коэффициентами. Методом энергетических неравенств получены априорные оценки в дифференциальной и разностной трактовках для решений нелокальных краевых задач. Из полученных оценок следуют единственность и устойчивость решения каждой из рассмотренных задач по правой части и начальным данным, а также сходимость решения разностной задачи к решению исходной дифференциальной задачи в L2-норме со скоростью O(|h|+τ). Для каждой из рассмотренных задач построен алгоритм численного решения, роведены численные расчеты тестовых примеров.

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

We compare the error behavior of two methods used to find a numerical solution of the linear integro-differential Fredholm equation with a weakly singular kernel in Banach space C1[a,b]. We construct an approximation solution based on the modified cubic b-spline collocation method. Another estimation of the exact solution, constructed by applying the numerical process of product and quadrature integration, is considered as well. Two proposed methods lead to solving a linear algebraic system. The stability and convergence of the cubic b-spline collocation estimate is proved. We test these methods on the concrete examples and compare the numerical results with the exact solution to show the efficiency and simplicity of the modified collocation method.

В работе представлен алгоритм решения системы уравнений Аллена-Кана и Кана-Хиллиарда, которая описывает процесс спекания. Алгоритм не требует значительных по мощности вычислительных ресурсов и позволяет выполнить моделирование процесса спекания большого количества отдельных частиц на вычислительном узле с процессором Intel Xeon E5 2697 v3 и графическим ускорителем NVIDIA K40 за приемлемое время. Проведены эксперименты по моделированию спекания сорбентоподобных структур - упаковок сферических частиц, и на них показана эффективность алгоритма.