For transport modelled by the time-dependent convection-diffusion equation, positivity of numerical concentrations and mass conservation are two important properties numerical solvers should respect. This paper investigates such a solver based on the implicit Euler time-marching and finite volume discretization on quadrilateral meshes. The solver uses mapped \( Q_1 \) bilinear polynomials for approximation of the concentration. A new upwinding technique is adopted to handle convection dominance. Flux correction is devised to ensure nonnegative numerical concentrations. Matlab code modules based on efficient implementation of this solver are incorporated into our package DarcyLite. Numerical experiments are presented to illustrate the performance of the new solver.

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

DarcyLite Modules for a Property-Preserving Transport Solver

  • Jiangguo Liu,
  • Boyang Yu,
  • Yingli Li

摘要

For transport modelled by the time-dependent convection-diffusion equation, positivity of numerical concentrations and mass conservation are two important properties numerical solvers should respect. This paper investigates such a solver based on the implicit Euler time-marching and finite volume discretization on quadrilateral meshes. The solver uses mapped \( Q_1 \) bilinear polynomials for approximation of the concentration. A new upwinding technique is adopted to handle convection dominance. Flux correction is devised to ensure nonnegative numerical concentrations. Matlab code modules based on efficient implementation of this solver are incorporated into our package DarcyLite. Numerical experiments are presented to illustrate the performance of the new solver.