Комплексная численная модель медленного течения многофазной жидкости

Разработана двумерная комплексная численная модель медленного течения многофазной жидкости в расчетной области, состоящей из относительно толстого слоя двухфазной жидкости, покрытого тонким многослойным вязким пластом. На границе сопряжения разнородных подобластей происходит массобмен между легким компонентом двухфазного слоя и нижним слоем многослойного пласта. Общая система уравнений соединяет в себе уравнения вязкой компакции, описывающие течение в двухфазном слое, с уравнениями Рейнольдса - в пласте. В модели можно задавать структуру пласта с любой степенью подробности, а также поверхностные процессы, в частности эрозию, денудацию и осадконакопление. Использовалось дополнительное асимптотическое граничное условие, которое связывает разнородные уравнения гидродинамики без каких-либо процедур итерационного уточнения. Это условие дает возможность значительно сократить вычислительные затраты по сравнению с затратами, необходимыми большинству ранее разработанных комплексных моделей. Проведено численное исследование эволюции поля скоростей и границ раздела слоев. Результаты расчета показали, что на малых и больших временах наблюдаются различия в режимах эволюции поля скоростей и границ слоев. Эволюцию можно разделить, по крайней мере, на две стадии с характерными масштабами времени: быструю на малых временах и медленно изменяющуюся (квазистационарную) на больших временах. Многостадийность эволюции течения определяется не внешними факторами, а геометрическими и физическими параметрами моделируемой среды. В качестве геофизического приложения показан численный расчет процесса образования области разуплотненной мантии под земной корой в активной зоне перехода океан-континент.