инженер-исследователь
Research engineer
кандидат физико-математических наук, заведующий лабораторией геомиграционного моделирования
старший научный сотрудник
candidate of physical and mathematical sciences, chief of the geomigration laboratory of the Nuclear safety institute of the RAS, Senior Scientific Researcher of the
доктор физико-математических наук, доцент, член-корреспондент Российской академии наук, заместитель директора
doctor of physical and mathematical sciences, associate professor, corresponding member of the Russian Academy of Sciences, Deputy Director
В статье рассматривается задача моделирования тепловой конвекции в пористых средах с учетом объемного тепловыделения, возникающая при оценке безопасности захоронений радиоактивных отходов (РАО). В первом разделе представлен краткий обзор известных гидрогеологических расчетных кодов, обладающих функционалом для решения тепловых задач (FEFLOW, SUTRA, SEAWAT, TOUGH2). Отмечается отсутствие возможности учета объемного тепловыделения вследствие радиоактивного распада в данных программных комплексах. Дается описание расчетного кода GeRa собственной разработки. Во втором разделе рассмотрена математическая модель сопряженных процессов фильтрации, массо- и теплопереноса, реализованная в коде GeRa. Модель описывает процессы в насыщенной пористой среде и учитывает радиоактивный распад, сорбцию на породе, зависимость плотности и вязкости жидкости от температуры. Уравнение теплопереноса записывается в предположении теплового равновесия между раствором и породой. Учитывается перенос тепла конвекцией и кондукцией-термодисперсией, источники в виде скважин и объемного тепловыделения вследствие радиоактивного распада. В третьем разделе представлен численный метод, реализованный в коде GeRa для решения сопряженной задачи. Дискретизация по пространству осуществляется с использованием методов конечных объемов (МКО). Для дискретизации по времени применяется схема полностью неявного итерационного сопряжения, на каждой итерации которой последовательно решаются задачи фильтрации, тепло- и массопереноса. В четвертом разделе приводится тестовая задача о конвекции тепловыделяющей жидкости в замкнутой двумерной полости с изотермическими стенками, заполненной пористой средой. Результаты решения кодом GeRa сравниваются с асимптотическим решением, полученным Haajizadeh. В пятом разделе представлены результаты моделирования в коде GeRa экспериментов Buretta и Berman по изучению режимов свободной тепловой конвекции жидкости с объемным тепловыделением в пористой среде. Сравниваются зависимости числа Нуссельта от числа Рэлея, полученные в натурном эксперименте и численно. В шестом разделе рассматривается тестовая задача о непрерывной закачке высокоактивных РАО в пласт-коллектор. Демонстрируется возможность совместного расчета процессов фильтрации, тепло- и массопереноса. Численное решение, полученное в коде GeRa, сравнивается с известным аналитическим решением.
This article is devoted to the problem of thermal convection in porous media with volumetric heat generation modelling, arising in practice of radioactive waste (RW) disposal safety assessment. In the first section a brief overview of widespread hydrogeological codes (FEFLOW, SUTRA, SEAWAT, TOUGH2) featuring the ability to solve thermal problems is done. We point out the lack of heat generation caused by radioactive decay model in these programs. The GeRa numerical code developed by the authors is presented. In the second section we consider the mathematical model of coupled groundwater flow, solute and heat transport, which is implemented in GeRa. The model describes these processes in saturated porous media and takes into account radioactive decay, sorption on the rock, the dependences of density and viscosity on temperature. The heat transport equation is written assuming thermal equilibrium between the fluid and the rock. The model includes heat transport by convection and conduction-thermal dispersion. The heat source terms can be wells and volumetric heat generation due to radioactive decay. The numerical scheme implemented in GeRa to solve the aforementioned coupled problem is introduced in the third section. The space discretization is done using finite volume methods (FVM). Sequential iterative coupling implicit scheme is used for temporal discretization. On each iteration of the scheme the flow, heat transport and solute transport problems are solved sequentially. The fourth section is devoted to the test problem of heat generating fluid convection in a closed two- dimensional cavern filled by porous material with isothermal walls. The results obtained using GeRa code are compared to the asymptotical solution deduced by Haajizadeh. In the fifth section we present the results of modelling with GeRa the experiments of Buretta and Berman in which they investigated the regimes of free thermal convection of fluid with volumetric heat generation in porous media. The dependences of Nusselt number on the Raley number measured in the experiments and calculated numerically are compared. In the sixth section we consider the test problem of continuous injection of high-level RW into an aquifer. Here the ability to model coupled flow, heat and solute transport processes is shown. The numerical solution obtained using GeRa is compared to a known analytical one.
The authors declare that there are no conflicts of interest present.