<p>In this paper, a second-order generalized Riemann problem (GRP) solver is developed for a two-layer thin film model. Extending the first-order Godunov approach, the solver is used to construct a temporal-spatial coupled second-order GRP-based finite-volume method. Numerical experiments including comparisons to MUSCL finite-volume schemes with Runge-Kutta time stepping confirm the accuracy, efficiency and robustness of the higher-order ansatz. The construction of GRP methods requires to compute temporal derivatives of intermediate states in the entropy solution of the generalized Riemann problem. These derivatives are obtained from the Rankine-Hugoniot conditions as well as a characteristic decomposition using Riemann invariants. Notably, the latter can be computed explicitly for the two-layer thin film model, which renders this system to be very suitable for the GRP approach. Moreover, it becomes possible to determine the derivatives in an explicit, computationally cheap way.</p>

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

A Generalized Riemann Problem Solver for a Hyperbolic Model of Two-Layer Thin Film Flow

  • Rahul Barthwal,
  • Christian Rohde,
  • Yue Wang

摘要

In this paper, a second-order generalized Riemann problem (GRP) solver is developed for a two-layer thin film model. Extending the first-order Godunov approach, the solver is used to construct a temporal-spatial coupled second-order GRP-based finite-volume method. Numerical experiments including comparisons to MUSCL finite-volume schemes with Runge-Kutta time stepping confirm the accuracy, efficiency and robustness of the higher-order ansatz. The construction of GRP methods requires to compute temporal derivatives of intermediate states in the entropy solution of the generalized Riemann problem. These derivatives are obtained from the Rankine-Hugoniot conditions as well as a characteristic decomposition using Riemann invariants. Notably, the latter can be computed explicitly for the two-layer thin film model, which renders this system to be very suitable for the GRP approach. Moreover, it becomes possible to determine the derivatives in an explicit, computationally cheap way.