Scholarlink

<p>We consider a compressible Stokes problem in the quasi-stationary case coupled with a time dependent advection-diffusion equation with special emphasis on high viscosity contrast geophysical mantle convection applications. In space, we use a P2-P1 Taylor–Hood element which is generated by a blending approach to account for the non-planar domain boundary without compromising the stencil data structure of uniformly refined elements. In time, we apply an operator splitting approach for the temperature equation combining the BDF2 method for diffusion and a particle method for advection, resulting in an overall second order scheme. Within each time step, a stationary Stokes problem with a high viscosity contrast has to be solved for which we propose a matrix-free, robust and scalable iterative solver based on Uzawa type block preconditioners, polynomial Chebyshev smoothers and a BFBT type Schur complement approximation. Our implementation is using a hybrid hierarchical grid approach allowing for massively parallel, high resolution Earth convection simulations.</p>

错误：搜索内容不能为空，请输入英文关键词
错误：关键词超出字数限制，请精简
高级检索

A robust matrix-free approach for large-scale non-isothermal high-contrast viscosity Stokes flow on blended domains with applications to geophysics

  • 期刊论文

  • 出版模式：

    开放获取

  • 发表日期：

    2026年1月7日

查看全文
  • Andreas Burkhart,
  • Nils Kohl,
  • Barbara Wohlmuth,
  • Jan Zawallich

摘要

We consider a compressible Stokes problem in the quasi-stationary case coupled with a time dependent advection-diffusion equation with special emphasis on high viscosity contrast geophysical mantle convection applications. In space, we use a P2-P1 Taylor–Hood element which is generated by a blending approach to account for the non-planar domain boundary without compromising the stencil data structure of uniformly refined elements. In time, we apply an operator splitting approach for the temperature equation combining the BDF2 method for diffusion and a particle method for advection, resulting in an overall second order scheme. Within each time step, a stationary Stokes problem with a high viscosity contrast has to be solved for which we propose a matrix-free, robust and scalable iterative solver based on Uzawa type block preconditioners, polynomial Chebyshev smoothers and a BFBT type Schur complement approximation. Our implementation is using a hybrid hierarchical grid approach allowing for massively parallel, high resolution Earth convection simulations.