Variational multiscale finite element methods for a nonlinear convection-diffusion-reaction equation

M.S. Zhelnin, A.A. Kostina, O.A. Plekhov
2019 Computational Continuum Mechanics  
Работа посвящена построению вариационных многомасштабных методов конечных элементов для численного решения двумерных краевых задач с сингулярно-возмущенным нестационарным нелинейным уравнением конвекции-диффузии-реакции. Решения данных задач могут быстро изменяться в тонких слоях, что при применении стандартной расчетной схемы Галёркина приводит к возникновению в этих областях нефизических осцилляций. В вариационных многомасштабных методах выполняется разложение исходной задачи на сеточную и
more » ... сеточную, что позволяет учесть особенности задачи на масштабах, меньших размера элемента сетки. В данной работе рассматриваются два многомасштабных метода: VMM-ASA (Variational Multiscale Method with Algebraic Sub-scale Approximation) и RFB (Residual-Free Bubbles) метод. В первом из них подсеточная задача аппроксимируется с использованием невязки сеточного уравнения и стабилизирующих параметров. Во втором подсеточная задача решается приближенно на основе аппроксимационных функций специального вида. Постановки сеточной и подсеточной задач определяются посредством линеаризации исходной задачи по подсеточной компоненте. Компьютерная реализация методов выполнена в коммерческом пакете конечно-элементного моделирования. Эффективность предложенных методов исследована путем решения модельной краевой задачи с нелинейным уравнением. Рассмотрены случаи различной величины коэффициента диффузии. В результате вычислительных экспериментов показано, что по сравнению со стандартной расчетной схемой Галёркина многомасштабные методы дают возможность достигать более устойчивого численного решения как с меньшим количеством осцилляций, так и их меньшей амплитудой. При малой величине коэффициента диффузии, когда схема Галёркина расходится, стабилизированные методы обеспечивают приемлемое численное решение на достаточно грубой сетке. Ключевые слова: уравнение конвекции-диффузии-реакции, стабилизированный метод конечных элементов, вариационный многомасштабный метод, осцилляции численного решения This paper focuses on the development of finite element methods for solving a two-dimensional boundary value problem for a singularly perturbed time-dependent convection-diffusion-reaction equation. Solution to the problem can vary rapidly in thin layers. As a result, spurious oscillations in the solution occur if the standard Galerkin method is used. In multiscale finite element methods, the initial problem is split into grid-scale and subgrid-scale problems, which allows one to capture the features of the problem at a scale smaller than an element mesh size. In the study two methods are considered: VMM-ASA (Variational Multiscale Method with Algebraic Sub-scale Approximation) and RFB (Residual-Free Bubbles). In the first method, the subgrid problem is modeled by the residual of the grid equation and intrinsic time scales. In the second method, the subgrid problem is approximated by special functions. The grid and subgrid problems are formulated through a linearization procedure on the subgrid component applied to the initial problem. The computer implementation of the methods has been carried out using a commercial finite element package. The efficiency of the developed methods has been studied by solving a test boundary value problem for the nonlinear equation. Several values of the diffusion coefficient of the equation have been analyzed. On the basis of the numerical study, it has been shown that the multiscale methods allow one to increase the stability of a numerical solution and to decrease the quantity and amplitude of oscillations compared to the standard Galerkin method. In the case of a small diffusion coefficient, the developed methods can yield a satisfactory numerical solution on a sufficiently coarse mesh.
doi:10.7242/1999-6691/2019.12.2.13 fatcat:j6grvkvndvgrva4c4pikoasyiy