Трехмерная совместная численная модель медленного течения вязкой жидкости
Автор: Пак Владимир Васильевич
Журнал: Вычислительная механика сплошных сред @journal-icmm
Статья в выпуске: 1 т.8, 2015 года.
Бесплатный доступ
Разработана трехмерная совместная численная модель для описания медленного течения в расчетной области, состоящей из толстого вязкого слоя, покрытого тонким многослойным вязким пластом. Предполагалось, что пласт имеет меньшую плотность, чем слой. Модель соединяет в себе уравнения Стокса, описывающие течение в слое с уравнениями Рейнольдса в пласте. С помощью асимптотического метода для исследования эволюции течения на больших временах было получено обыкновенное дифференциальное уравнение относительно смещения границ пласта и скоростей на границе раздела между пластом и слоем. Это асимптотическое уравнение используется как внутреннее граничное условие, связывающее уравнения Стокса с уравнениями Рейнольдса. Численная реализация осуществлялась модифицированным методом конечных элементов в сочетании с методом проекции градиента. Дискретизация расчетной области осуществлялась с помощью прямоугольных параллелепипедов. В качестве базовых функций использовались квадратичные элементы. Представленная модель дает возможность соединить разнотипные уравнения гидродинамики без использования какого-либо итерационного уточнения. Это позволяет значительно сократить вычислительные затраты по сравнению с большинством ранее разработанных комплексных моделей. Численные результаты подтвердили хорошую точность предлагаемой трехмерной совместной модели.
Трехмерная совместная модель, уравнения стокса, уравнения рейнольдса, уравнения смазки, метод малого параметра, метод конечных элементов, метод проекции градиента
Короткий адрес: https://sciup.org/14320754
IDR: 14320754 | УДК: 550.347.62 | DOI: 10.7242/1999-6691/2015.8.1.6
Three-dimensional coupled numerical model of creeping flow of viscous fluid
A three-dimensional numerical coupled model is developed to describe creeping flow in a computational domain that consists of thick viscous layer overlaid by a thin multi-layered viscous sheet. The sheet is assumed to have the lower density than the layer. The model couples the Stokes equations describing the flow in the layer and the Reynolds equations describing the flow in the sheet. We investigate the long-time behavior of the flow in the sheet using the asymptotic method and derive an ordinary differential equation with respect to sheet boundary displacements and velocities at the interface between the sheet and layer. Applying the obtained equation as an internal boundary condition, we couple the Stokes and Reynolds equations. Numerical implementation is fulfilled by the modified finite element method combined with the projection gradient method. The examined domain is discretized with hexahedrons. Piecewise square basic functions are used. The proposed model enables different-type hydrodynamic equations to be coupled without any iterative improvements. This reduces significantly computational costs in comparison with the available coupled models. Numerical experiments confirm a good accuracy of the developed three-dimensional coupled model.
Список литературы Трехмерная совместная численная модель медленного течения вязкой жидкости
- Schubert G., Turcotte D.L., Olsen P. Mantle convection in the Earth and planets. -Cambridge: Cambridge University Press, 2001. -956 p.
- Tan E., Choi E., Nhoutireddy P., Gurnis M., Aivazis M. GeoFramework: Coupling multiple models of mantle convection within a computational framework//Geochemistry. Geophysics. Geosystems. -2006. -Vol. 7, no. 6. -Q06001.
- Lavier L.L., Manatschal G. A mechanism to thin the continental lithosphere at magma-poor margins//Nature. -2006. -Vol. 440, no. 7082. -P. 324-328.
- Bonnardot M.-A., Hassani R., Tric E. Numerical modelling of lithosphere-asthenosphere interaction in a subduction zone//Earth Planet. Sc. Lett. -2008. -Vol. 272, no. 3-4. -P. 698-708.
- Пак В.В. Применение метола проекции градиента к численному решению совместной системы уравнений Стокса и уравнений Рейнольдса//Вычисл. мех. сплош. сред. -2014. -Т. 7, № 1. -С. 23-29.
- Craster R.V., Matar O.K. Dynamics and stability of thin liquid films//Rev. Mod. Phys. -2009. -Vol. 81, no. 3. -P. 1131-1198.
- Пак В.В. Применение метола проекции градиента к численному решению трехмерной задачи Стокса//Вычисл. мех. сплош. сред. -2010. -Т. 3, № 2. -С. 93-102.
- Fortin M. Old and new finite-elements for incompressible flows//Int. J. Numer. Meth. Fl. -1981. -Vol. 1, no. 4. -P. 347-364.
- Борисов В.Г. О параболических краевых задачах с малым параметром при производных по t//Матем. сб. -1988. -Т. 131(173), № 3(11). -C. 293-308.
- Hu J., Millet S., Botton V., Ben Hadid H., Henry D. Inertialess temporal and spatio-temporal stability analysis of the two-layer film flow with density stratification//Phys. Fluids. -2006. -Vol. 18, no. 10. -104101.
- Трубицын В.П., Баранов А.А., Евсеев А.Н., Трубицын А.П. Точные аналитические решения уравнения Стокса для тестирования уравнений мантийной конвекции с переменной вязкостью//Физика Земли. -2006. -№ 7. -С. 3-11.