Discrete and Continuous Models and Applied Computational ScienceDiscrete and Continuous Models and Applied Computational ScienceDiscrete and Continuous Models and Applied Computational Science2658-46702658-7149Peoples' Friendship University of Russia named after Patrice Lumumba (RUDN University)8377ArticlesСтатьиResearch ArticleMathematical and Computational Oil Spills ModellingМатематическое и компьютерное моделирование эволюции нефтяного пятнаTrepachevaA VТрепачеваАлина ВикторовнаЮжно-Российский региональный центр информатизацииatrepacheva@sfedu.ruBurtykaPh BБуртыкаФилипп БорисовичЮжно-Российский региональный центр информатизацииfburtyka@sfedu.ruComputer Center of Southern Federal UniversityЮжный федеральный университет150220142NO2 (2014)№2 (2014)28128608092016Copyright ©; 2014, Трепачева А.В., Буртыка Ф.Б.2014Трепачева А.В., Буртыка Ф.Б.http://creativecommons.org/licenses/by/4.0https://journals.rudn.ru/miph/article/view/8377The model of oil spill evolution on the water surface based on quasi-linear convection-diffusion equation is considered. To solve the latter the method of finite differences is used. Two-dimensional problem is reduced to one-dimensional by alternating direction method. For one-dimensional case traditional difference approximations are briefly discussed. And also difference schemes obtained by computer algebra based method of automatic generation of difference schemes are considered. For generated schemes the order of approximation and linear numerical dissipation are estimated. For explicit schemes stability conditions are briefly discussed. In the absence of convection a numerical comparison of traditional implicit difference scheme and one implicit scheme generated automatically with similarity solution of quasi-linear diffusion equation is carried out. This comparison shows that generated implicit scheme allows to obtain the relative error less than for traditional scheme. Using obtained implicit scheme evaluation of oil slick thickness changes in the presence of oil evaporation is done. Two different models were used to estimate oil evaporation rate. The first one is based on hypothesis that oil evaporation is regulated by air boundary layer. The second assumes that oil evaporation is regulated by oil diffusion. Calculations show that the choice of model essentially influences on oil slick thickness and oil total volume changes in time.В работе рассматривается модель переноса нефтяного пятна по поверхности воды, базирующаяся на квазилинейном уравнении конвекции-диффузии, составленном относительно толщины пятна. Для решения последнего используется метод конечных разностей. Двумерная задача сводится к одномерной с помощью метода переменных направлений. Для одномерного случая рассматриваются традиционные разностные аппроксимации, а также разностные схемы, которые можно получить с помощью метода автоматической генерации разностных схем, основанного на методах компьютерной алгебры. Для сгенерированных схем оценивается порядок аппроксимации и схемная вязкость. Для явных схем рассматривается вопрос об условиях устойчивости. В случае отсутствия конвекции проводится численное сравнение автомодельного решения квазилинейного уравнения теплопроводности с традиционной неявной разностной схемы и одной неявной схемой, сгенерированной автоматически. Сравнение показало, что сгенерированная неявная схема позволяет получить относительную погрешность решения меньшую, чем для традиционной схемы. В соответствии с полученной неявной схемой проводится расчёт изменения толщины пятна при условии наличия испарения нефти с поверхности пятна. При этом были опробованы две различные модели испарения, одна из которых предполагает, что испарение нефти регулируется в основном приграничным воздушным слоем, а вторая исходит из предположения, что испарение нефти носит в большей степени диффузионный характер. Расчёты показали, что выбор модели испарения существенно влияет на изменения толщины пятна и суммарного объёма нефти с течением времени.quasilinear advection-diffusion equationfinite difference methodGrorbner basisdifference schemes generationoil spillsoil evaporationквазилинейное уравнение конвекции-диффузииметод конечных разностейгенерация разностных схембазис Грёбнераиспарение нефти1.Fingas M. Oil Spills Science and Technology: Prevention, Response and Cleanup. - Burlington: Elseivier, 2011.2.Tkalich P. A CFD Solution of Oil Spill Problems // Environ. Modeling and Software. - 2006. - Vol. 21. - Pp. 271-282.3.Роуч П. Вычислительная гидродинамика. - М.: Мир, 1980.4.Zadeh E., Hejazi K. Eulerian Oil Spills Model using Finite-Volume Method with Moving Boundary and Wet-Dry Fronts // Model. Simul. Eng. - 2012. - Vol. 2012. - P. 33.5.Тихонов А.Н., Самарский А.А. Уравнения математической физики. - М.: Наука, 1977.6.Gerdt V.P., Blinkov Y.A., Mozzhilkin V.V. Gr¨obner Bases and Generation of Difference Schemes for Partial Differential Equations // Symmetry, Integrability and Geometry: Methods and Applications. - 2006. - Vol. 2. - P. 26.7.Gerdt V.P., Robertz D. Maple Package for Computing Gr¨obner Bases for Linear Recurrence Relations // Nuclear Instruments and Methods in Physics Research. - 2006. - Vol. 2. - P. 26.8.Gerdt V.P., Blinkov Y.A. On Computer Algebra-Aided Stability Analysis of Difference Schemes Generated by Means of Gr¨obner Bases // Computer Algebra and Differential Equations, Acta Academiae Aboensis, Ser. B. - 2007. - Vol. 67, No 2. - Pp. 168-177.9.Stiver W., Mackay D. Evaporation Rate of Spills of Hydrocarbons and Petroleum Mixtures // Computer Algebra and Differential Equations, Acta Academiae Aboensis, Ser. B. - 1984. - Vol. 18. - Pp. 834-840.