Simulation of thin film flows with a moving mesh mixed finite element method

We present an efficient mixed finite element method to solve the fourth-order thin film flow equations using moving mesh refinement. The moving mesh strategy is based on harmonic mappings developed by Li et al. [J. Comput. Phys., 170 (2001), pp. 562-588, and 177 (2002), pp. 365-393]. To achieve a high quality mesh, we adopt an adaptive monitor function and smooth it based on a diffusive mechanism. A variety of numerical tests are performed to demonstrate the accuracy and efficiency of the method. The moving mesh refinement accurately resolves the overshoot and downshoot structures and reduces the computational cost in comparison to numerical simulations using a fixed mesh.Comment: 18 pages, 10 figure

Numerical investigations of two-phase flow with dynamic capillary pressure in porous media via a moving mesh method

{©} 2017 Elsevier Inc.Motivated by observations of saturation overshoot, this paper investigates numerical modeling of two-phase flow in porous media incorporating dynamic capillary pressure. The effects of the dynamic capillary coefficient, the infiltrating flux rate and the initial and boundary values are systematically studied using a traveling wave ansatz and efficient numerical methods. The traveling wave solutions may exhibit monotonic, non-monotonic or plateau-shaped behavior. Special attention is paid to the non-monotonic profiles. The traveling wave results are confirmed by numerically solving the partial differential equation using an accurate adaptive moving mesh solver. Comparisons between the computed solutions using the Brooks–Corey model and the laboratory measurements of saturation overshoot verify the effectiveness of our approach

Numerical investigations of two-phase flow with dynamic capillary pressure in porous media via a moving mesh method

{©} 2017 Elsevier Inc.Motivated by observations of saturation overshoot, this paper investigates numerical modeling of two-phase flow in porous media incorporating dynamic capillary pressure. The effects of the dynamic capillary coefficient, the infiltrating flux rate and the initial and boundary values are systematically studied using a traveling wave ansatz and efficient numerical methods. The traveling wave solutions may exhibit monotonic, non-monotonic or plateau-shaped behavior. Special attention is paid to the non-monotonic profiles. The traveling wave results are confirmed by numerically solving the partial differential equation using an accurate adaptive moving mesh solver. Comparisons between the computed solutions using the Brooks–Corey model and the laboratory measurements of saturation overshoot verify the effectiveness of our approach