4.6 Article

An augmented Cartesian grid method for Stokes-Darcy fluid-structure interactions

Journal

Publisher

WILEY
DOI: 10.1002/nme.5131

Keywords

Stokes-Darcy fluid-structure interaction; augmented immersed interface method; fast Poisson solver; BJS interface condition; analytic solution of Stokes and Darcy coupling

Funding

  1. US NSF [DMS-1522768]
  2. NIH [5R01GM96195-2]
  3. CNSF [11371199, 11471166, BK20141443]
  4. Direct For Mathematical & Physical Scien
  5. Division Of Mathematical Sciences [1522768] Funding Source: National Science Foundation

Ask authors/readers for more resources

A new finite difference method based on Cartesian meshes is proposed for solving the fluid-structure interaction between a fluid flow modeled by the Stokes equations and a porous media modeled by the Darcy's law. The idea is to introduce several augmented variables along the interface between the fluid flow and the porous media so that the problem can be decoupled as several Poisson equations. The augmented variables should be chosen so that the Beavers-Joseph-Saffman and other interface conditions are satisfied. In the discretization, the augmented variables have co-dimension one compared with that of the primitive variables and are solved through the Schur complement system. A non-trivial analytic solution with a circular interface is constructed to check second-order convergency of the proposed method. Numerical examples with various interfaces and parameters are also presented. Some simulations show interesting behaviors of the fluid-structure interaction between the fluid flow and the porous media. The computational framework can be applied to other multi-phase and multi-physics problems. Copyright (c) 2015 John Wiley & Sons, Ltd.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.6
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Mathematics, Applied

A generalized modulus-based Newton method for solving a class of non-linear complementarity problems with P-matrices

Rui Li, Zhi-Lin Li, Jun-Feng Yin

Summary: Inspired by the modulus-based Newton method for linear complementarity problems, a generalized modulus-based Newton method was proposed in this study to solve a class of non-linear complementarity problems with P-matrices. A sufficient condition for convergence was obtained and numerical experiments showed that the proposed method is efficient and outperforms existing methods.

NUMERICAL ALGORITHMS (2022)

Article Computer Science, Interdisciplinary Applications

A new FV scheme and fast cell-centered multigrid solver for 3D anisotropic diffusion equations with discontinuous coefficients

Kejia Pan, Xiaoxin Wu, Hongling Hu, Yunlong Yu, Zhilin Li

Summary: This paper introduces a new extrapolation cascadic multigrid method for solving 3D anisotropic diffusion equations, avoiding the costly solution of local linear systems and demonstrating its efficiency and robustness.

JOURNAL OF COMPUTATIONAL PHYSICS (2022)

Article Mathematics, Applied

Optimal convergence of three iterative methods based on nonconforming finite element discretization for 2D/3D MHD equations

Jiali Xu, Haiyan Su, Zhilin Li

Summary: This paper focuses on analyzing nonconforming iterative finite element methods for 2D/3D stationary incompressible magneto-hydrodynamics equations. It introduces a locally stabilization term to address the instability issue in the finite element method for velocity field and pressure, along with proposing three effective iterative methods. The theoretical analysis is verified through numerical experiments in the end.

NUMERICAL ALGORITHMS (2022)

Article Mathematics, Applied

A new patch up technique for elliptic partial differential equation with irregularities

Swarn Singh, Suruchi Singh, Zhilin Li

Summary: This paper presents a new technique based on cubic spline interpolation for solving second order elliptic equations with irregularities in one dimension and two dimensions. The proposed method provides second order accurate solutions at the interface and regular points, as well as for the first and second order derivatives. Numerical experiments demonstrate the effectiveness and accuracy of the method for 2D problems.

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS (2022)

Article Mathematics, Applied

A Variant Modified Skew-Normal Splitting Iterative Method for Non-Hermitian Positive Definite Linear Systems

Rui Li, Jun-Feng Yin, Zhi-Lin Li

Summary: We propose a variant modified skew-normal splitting iterative method to solve a class of large sparse non-Hermitian positive definite linear systems. The preconditioned version of the proposed method is also constructed using the preconditioning technique. Theoretical analysis shows the unconditional convergence of the proposed method. Numerical experiments demonstrate the efficiency and superior performance of the proposed method compared to other methods.

NUMERICAL MATHEMATICS-THEORY METHODS AND APPLICATIONS (2022)

Article Mathematics, Applied

Numerical analysis of a free boundary problem with non-local obstacles

Zhilin Li, Hayk Mikayelyan

Summary: This paper deals with the challenging obstacle-like minimization problem in a cylindrical domain, where the non-local nature of the obstacle makes it difficult to apply standard optimization techniques. The paper introduces a new algorithm that can compute the global minimum and provides numerical testing and construction of exact solutions.

APPLIED MATHEMATICS LETTERS (2023)

Article Engineering, Multidisciplinary

An efficient extrapolation multigrid method based on a HOC scheme on nonuniform rectilinear grids for solving 3D anisotropic convection-diffusion problems

Shuanggui Hu, Kejia Pan, Xiaoxin Wu, Yongbin Ge, Zhilin Li

Summary: This paper presents an efficient multigrid method combined with a high-order compact finite difference scheme on nonuniform rectilinear grids for solving 3D diagonal anisotropic convection-diffusion problems with boundary/interior layers. The method involves a fourth-order compact finite difference scheme, novel multigrid prolongation operators based on quintic Lagrange interpolation and completed Richardson extrapolation, and a SSOR-preconditioned biconjugate gradient stabilized smoother. Numerical experiments demonstrate fourth-order accuracy and improved efficiency compared to the state-of-the-art algebraic multigrid method for large linear systems arising from second order elliptic PDE discretization.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2023)

Article Mathematics, Applied

Adapting Cubic Hermite Splines to the Presence of Singularities Through Correction Terms

Sergio Amat, Zhilin Li, Juan Ruiz-Alvarez, Concepcion Solano, Juan C. Trillo

Summary: This paper presents a new technique for computing accurate first order derivatives of a piecewise smooth function close to singularities using a Hermite spline. The technique corrects the system of equations to achieve O(h(4)) accuracy even near singularities. Theoretical proofs and numerical experiments confirm the effectiveness of the proposed method in eliminating oscillations and smearing of singularities.

JOURNAL OF SCIENTIFIC COMPUTING (2023)

Article Mathematics, Applied

New Sixth-Order Compact Schemes for Poisson/Helmholtz Equations

Kejia Pan, Kang Fu, Jin Li, Hongling Hu, Zhilin Li

Summary: Some new sixth-order compact finite difference schemes for Poisson/ Helmholtz equations on rectangular domains in both two-and three-dimensions are developed and analyzed. The finite difference and weight coefficients of the new methods have analytic simple expressions, different from a few schemes in the literature. The new compact schemes are exact for 2D and 3D Poisson equations if the solution is a polynomial less than or equal to 6, and the coefficient matrices are M-matrices for Helmholtz equations, guaranteeing the discrete maximum principle and convergence.

NUMERICAL MATHEMATICS-THEORY METHODS AND APPLICATIONS (2023)

Article Mathematics, Applied

HIGH ORDER COMPACT SCHEMES FOR FLUX TYPE BCS

Zhilin Li, Kejia Pan

Summary: In this paper, new fourth order compact schemes for Robin and Neumann boundary conditions are developed for boundary value problems of elliptic PDEs in two and three dimensions. These schemes utilize carefully designed undetermined coefficient methods and can be applied to various elliptic PDEs including both flux and linear boundary conditions. The developed schemes are versatile and generally have M-matrices as coefficient matrices, ensuring well-posed problems and convergence of the methods. Examples with large wave numbers and oscillatory solutions are presented to demonstrate the performance of the new schemes.

SIAM JOURNAL ON SCIENTIFIC COMPUTING (2023)

Article Computer Science, Interdisciplinary Applications

Numerical integration rules with improved accuracy close to discontinuities

Sergio Amat, Zhilin Li, Juan Ruiz-Alvarez, Concepcion Solano, Juan C. Trillo

Summary: This work introduces a novel nonlinear technique for accurate numerical integration of any order, even with data containing discontinuities and known only at grid points. The technique includes correction terms that depend on the size of jumps and derivatives of the function at discontinuities. These correction terms allow recovering the accuracy of classical integration formulas near the discontinuities and can be added during integration or as post-processing. Numerical experiments validate the theoretical conclusions of this article.

MATHEMATICS AND COMPUTERS IN SIMULATION (2023)

Article Mathematics, Applied

Stable high order FD methods for interface and internal layer problems based on non-matching grids

Zhilin Li, Kejia Pan, Juan Ruiz-Alvarez

Summary: This paper presents alternative approaches for interface and internal layer problems using non-matching grids. The developed methods achieve high order accuracy by using two mesh sizes near and away from the interface or internal layer. For problems with straight interfaces or boundary layers parallel to one axis, a fourth order compact finite difference scheme is constructed at border grid points. For problems with curved interfaces or internal layers, a level set representation is utilized to build a fine mesh within a tube around the interface. A new super-third seven-point discretization is developed for hanging nodes to guarantee the discrete maximum principle.

NUMERICAL ALGORITHMS (2023)

Article Mathematics, Applied

Accurate derivatives approximations and applications to some elliptic PDEs using HOC methods

Jin Li, Zhilin Li, Kejia Pan

Summary: In this study, new high order compact finite difference schemes are developed to accurately approximate the derivatives of the solutions to some elliptic partial differential equations (PDEs). The convergence analysis shows that the accuracy of the computed derivatives is the same as that of the solution. The developed schemes take into account the partial differential equations, including the source term and/or the boundary conditions, and have important applications in solving incompressible Stokes equations with periodic boundary conditions.

APPLIED MATHEMATICS AND COMPUTATION (2023)

Article Mathematics, Applied

Analysis of nonconforming IFE methods and a new scheme for elliptic interface problems

Haifeng Ji, Feng Wang, Jinru Chen, Zhilin Li

Summary: In this paper, a significant discovery has been made regarding nonconforming immersed finite element (IFE) methods for solving elliptic interface problems. It has been shown that IFE methods without penalties may not converge optimally in the presence of non-zero tangential derivative of the exact solution and coefficient jump on the interface. To address this issue, a new parameter-free nonconforming IFE method with additional terms on interface edges has been developed to recover the optimal convergence rates. The method has been shown to have unisolvent IFE basis functions on arbitrary triangles.

ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS (2023)

Article Mathematics, Applied

AN L8 SECOND ORDER CARTESIAN METHOD FOR 3D ANISOTROPIC INTERFACE PROBLEMS

Baiying Dong, Xiufeng Feng, Zhilin Li

Summary: A second order accurate method is proposed for general three dimensional anisotropic elliptic interface problems. The method handles finite jumps across smooth interfaces by deriving 3D interface relations and using different discretizations. The convergence analysis is shown by splitting errors and using various numerical examples.

JOURNAL OF COMPUTATIONAL MATHEMATICS (2022)

No Data Available