A least squares augmented immersed interface method for solving Navier–Stokes and Darcy coupling equations. (15th May 2018)

• Record Type: Journal Article
• Authors: Li, Zhilin; Lai, Ming-Chih; Peng, Xiaofei; Zhang, Zhiyue
• Is Part Of: Computers & fluids. Volume 167(2018)
• Page Start: 384
• Page End: 399
• Publication Date: 2018-05-15
• Subjects: Navier-Stokes and Darcy coupling -- Coupled fluid flow with porous media -- Least squares augmented IIM -- Fast Poisson/Helmholtz solver -- Analytic solution for Navier–Stokes and Darcy coupling -- Interface relations; 65M06 -- 65M85 -- 76M20; Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2018.03.032
• Languages: English
• ISSNs: 0045-7930

Abstract

Highlights: This is the probably first finite difference method based on Cartesian meshes and fast Poisson/Helmholtz solvers for Navier–Stokes and Darcy coupling which models the interactions between a fluid flow and a porous media. The new method has been applied to the Stokes and Darcy coupling as well. The equivalence of the original and transformed systems have been established with reasonable assumptions of regularities of the solutions though the sensitivity analysis of the new system is still an open question. Both of the velocity and pressure are second order accurate in the infinity norm for the Navier–Stokes and Darcys coupling. With the new method, we have obtained some interesting simulation results for time dependent problems that cannot be solved using the method for Stokes and Darcy coupling. The proposed method may provide an efficient way for other multi-phase and multiphysics problems. Abstract: A new finite difference method based on Cartesian meshes and fast Poisson/Helmholtz solvers is proposed to solve the coupling of a fluid flow modeled by the incompressible Navier–Stokes equations and a porous media modeled by the Darcy's law. The finite difference discretization in time is based on the pressure Poisson equation formulation. At each time step, several augmented variables along the interface between the fluid flow and the porous media are introduced so that the coupled equations can be decoupled into several Poisson/Helmholtz equations with those augmented variables acting as jumps of the unknown solution and some directional derivatives. The augmented variables should be chosen so that the Beavers–Joseph–Saffman (BJS) or Beavers–Joseph (BJ) and other interface conditions are satisfied. It has been tested that a direct extension of the augmented idea in [27] does not work well when the fluid flow is modeled by the Navier–Stokes equations. One of the new ideas of this paper is to enforce the divergence free condition at the interface from the fluid side. In this way, the Schur complement matrix for the augmented variables is over-determined and the least squares solution is used for the coupling problem. The new augmented approach enables us to solve the Navier–Stokes and Darcy coupling efficiently with second order accurate velocity and pressure in the L ∞ norm for tested problems. The proposed new idea in enforcing the divergence free condition at the interface from the fluid side has also been utilized to solve the Stokes and Darcy coupling equations and shown to outperform the original method in [27] . In additional to the detailed accuracy check for the present method, some interesting numerical simulations for Navier–Stokes and Darcy coupling have been conducted in this paper as well.