Все выпуски

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

We study the long-time evolution of axisymmetric free-surface two-layered creeping flow subject to the initial topography of its boundaries and bottom velocities. Each layer has uniform density and viscosity. The upper layer is assumed to have a smaller density than the lower layer. Based on lubrication approximation (the Reynolds equations) the nonlinear system of diffusion-type equations with respect to the surface and interface between the layers is obtained to describe this flow. Taking the dimensionless density contrast between the layers as a small parameter, we apply the method of asymptotic expansions to extract leading-term approximation for the slowly varying large-time evolution of the governing equations. An asymptotic equation relating both surface and interface displacement to the bottom velocities is derived. Based on this equation, we develop the algorithm to calculate velocity fields within layers for large time. Streamlines are used to visualize the flow. Numerical results reveal stability of the streamlines in the upper layer under variation of the bottom velocity. As geophysical applications, the developed algorithm is used to evaluate the velocity field in the crust (the upper layer) beneath the large-scale lunar multi-ring basins influenced by deep movements in the underlying mantle (the lower layer). To validate the results of modeling, we compare the calculated velocity fields with basin ridge systems obtained by experimental observations. The model comparison has shown proximity of radii of basin rings and critical points of the surface velocity.

На основе упрощенных уравнений Навье-Стокса в длинноволновом приближении построена нелинейная модель двухслойного течения вязкой жидкости со свободной границей, создаваемого начальным рельефом границ слоев. Используя метод малого параметра, исследуется эволюция течения на больших временах и определяется зависимость между движением поверхности и границы раздела жидкости. Полученные результаты применяются для расчета профиля границы кора-мантия под крупномасштабной кольцевой структурой на Луне.

The nonlinear model based on the long-wave approximation of the Navier-Stokes equations is developed to investigate the evolution of free-surface two-layered creeping flow subjected by the initial topography of the surface and interface between layers. Using the method of asymptotic expansions for the governing equations, we study a long-time evolution of the flow and reveal the relation between the surface and interface displacements. The obtained results are applied to calculate the profile of the crust-mantle interface beneath the large-scale lunar basin.

Предлагается осесимметрическая модель, построенная на основе уравнений Стокса, для исследования образования многокольцевой структуры в ползущем двухслойном течении с переменной толщиной слоев. Каждый слой имеет постоянную плотность и вязкость. Верхний слой имеет меньшую плотность, чем нижний. Течение создается рельефом поверхности и границы раздела слоев. Предполагается, что эффекты поверхностного натяжения пренебрежимо малы. Мы используем асимптотический метод многих масштабов для получения уравнений, описывающих неустойчивость, возникающую в виде волны в этом течении. С помощью преобразований Фурье и Лапласа мы исследуем уравнения главного приближения для этой неустойчивости в предположении малости возмущений. Асимптотическое исследование показывает, что эта неустойчивость проявляется в виде осесимметричной волны, длина которой соизмерима с толщиной слоев, и толщины слоев играют главную роль в пространственном распределении ее экстремумов. Остальные параметры модели влияют в основном на амплитуду волны. Получено уравнение, связывающее толщины слоев с распределением экстремумов, которое применяется для исследования закономерности расположения кольцевых хребтов, наблюдаемой для большинства крупномасштабных кольцевых структур на Луне. Используя параметры некоторых лунных кольцевых структур, мы определили радиусы последовательно расположенных экстремумов неустойчивости и провели сравнение модельных результатов с радиусами концентрических хребтов некоторых многокольцевых структур на Луне.

The axisymmetric model based on the Stokes equations is proposed to investigate the multi-ring pattern formation in two-layer creeping flow with variable thickness of layers. Each layer has uniform density and viscosity. The upper layer is lighter than the lower layer. The flow is generated by both surface and interface geometry. The effect of surface tension is supposed to be negligible. We apply the method of multiple scales to obtain the governing equations describing instability in the form of wave in the flow. Using the Fourier-Laplace method, we analyze the small-amplitude leading behavior of the instability. The asymptotic study reveals that this kind of instability manifests itself as axisymmetric wave which length is comparable with layer thickness; moreover, layer thicknesses play a major role in spatial distribution of wave extrema. The other model parameters alter mostly the wave amplitude. The equation relating extrema distribution to layer thicknesses is derived. We apply the obtained results to study a ring spacing rule observed for most multi-ring basins on the Moon. Using parameters of some lunar multi-ring basins we calculate the consecutive crest radii of the unstable wave and compare the results of simulation with the measured ring radii.

Разработана нелинейная модель трехслойного течения со свободной границей на основе упрощенных уравнений вязкой жидкости в длинноволновом приближении. Проведено асимптотическое исследование модели, которое показало существование двух различных режимов эволюции течения на малых и больших временах. Получено уравнение, связывающее смещения границ слоев на больших временах, не зависящее от предыстории течения. Модельные результаты используются для изучения поведения глубинной границы под крупномасштабной кольцевой структурой на Луне в зависимости от изменения геометрических физических параметров модели.

The nonlinear model based on the long-wave approximation of the Navier–Stokes equations is developed to study the free-surface three-layered creeping flow. An asymptotic study of the governing equations reveals two different modes of evolution at a short and long time. The relation between layers’ boundaries is obtained that is independent of a pre-history of the flow. The obtained results are applied to study a behavior of the deep interface beneath the large-scale lunar basin under the variation of geometrical and physical model’s parameters.