<p>We analyze the existence of solutions for the finite-element discretized stationary and time-dependent sedimentation-consolidation partial differential equations (PDEs) using recent advances in the local-projection stabilization (LPS) method. The sedimentation model studied here couples gravitationally forced Stokes flow and a convection-diffusion equation for the solids concentration, with boundary conditions selected to mimic physical applications. This system is highly nonlinear and sensitive due to the non-constant velocity and concentration coupling terms, as well as the nonlinear flux term in the convection equation. We provide an overview of pre-existing methods already developed for solving this class of problem, establish the coercivity of the underlying operator, as well as conditions for the existence of discrete solutions. For convection-dominated regimes using non-<i>inf-sup</i> stable finite elements with LPS, we demonstrate that stabilization effectively eliminates interior and boundary layers in the velocity and pressure solutions. Numerical examples in two and three dimensions for the non-stationary case illustrate the theoretical results.</p>

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

A local-projection stabilized sedimentation model with convection-dominant flow

  • Simon Taylor,
  • Malte Braack

摘要

We analyze the existence of solutions for the finite-element discretized stationary and time-dependent sedimentation-consolidation partial differential equations (PDEs) using recent advances in the local-projection stabilization (LPS) method. The sedimentation model studied here couples gravitationally forced Stokes flow and a convection-diffusion equation for the solids concentration, with boundary conditions selected to mimic physical applications. This system is highly nonlinear and sensitive due to the non-constant velocity and concentration coupling terms, as well as the nonlinear flux term in the convection equation. We provide an overview of pre-existing methods already developed for solving this class of problem, establish the coercivity of the underlying operator, as well as conditions for the existence of discrete solutions. For convection-dominated regimes using non-inf-sup stable finite elements with LPS, we demonstrate that stabilization effectively eliminates interior and boundary layers in the velocity and pressure solutions. Numerical examples in two and three dimensions for the non-stationary case illustrate the theoretical results.