Two-dimensional seepage analysis using a polygonal cell-based smoothed finite element method

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Yang, Yang, Yan, Mingjiao, Zhang, Zongliang, Yin, Yinpeng, Liu, Qiang, Li, You liang
Format: Preprint
Published: 2026
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866917278654660608
author Yang, Yang
Yan, Mingjiao
Zhang, Zongliang
Yin, Yinpeng
Liu, Qiang
Li, You liang
author_facet Yang, Yang
Yan, Mingjiao
Zhang, Zongliang
Yin, Yinpeng
Liu, Qiang
Li, You liang
contents This study develops a polygonal cell-based smoothed finite element method (CSFEM) for two-dimensional seepage analyses in porous media, covering steady-state, transient, and free-surface problems. Wachspress interpolation on convex polygonal elements is combined with cell-based gradient smoothing, so that element matrices are assembled using boundary integrals only, avoiding in-element derivatives and improving robustness on distorted and locally refined meshes. To improve efficiency, a solution-driven adaptive refinement strategy is employed to concentrate resolution near steep hydraulic gradients and evolving wet-dry interfaces. Free-surface seepage is handled by a fixed-mesh iterative scheme that updates the wetted region and boundary conditions to track the phreatic surface. Benchmark tests validate the formulation against analytical solutions and high-fidelity FEM references. In steady seepage examples, the proposed polygonal CSFEM reproduces linear hydraulic-head fields to near machine precision and yields smaller head errors than conventional FEM at the same characteristic mesh size. In transient problems, accurate head evolution and stable time responses are obtained, while adaptive refinement efficiently resolves localized high-gradient zones. For free-surface cases, the method captures the phreatic-surface profile and seepage-face development reliably without remeshing. The quadtree refinement and adaptivity provide substantial efficiency gains in degrees of freedom and runtime for a prescribed accuracy level.
format Preprint
id arxiv_https___arxiv_org_abs_2601_19357
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Two-dimensional seepage analysis using a polygonal cell-based smoothed finite element method
Yang, Yang
Yan, Mingjiao
Zhang, Zongliang
Yin, Yinpeng
Liu, Qiang
Li, You liang
Numerical Analysis
65N30, 74S05
F.2.2; I.2.7
This study develops a polygonal cell-based smoothed finite element method (CSFEM) for two-dimensional seepage analyses in porous media, covering steady-state, transient, and free-surface problems. Wachspress interpolation on convex polygonal elements is combined with cell-based gradient smoothing, so that element matrices are assembled using boundary integrals only, avoiding in-element derivatives and improving robustness on distorted and locally refined meshes. To improve efficiency, a solution-driven adaptive refinement strategy is employed to concentrate resolution near steep hydraulic gradients and evolving wet-dry interfaces. Free-surface seepage is handled by a fixed-mesh iterative scheme that updates the wetted region and boundary conditions to track the phreatic surface. Benchmark tests validate the formulation against analytical solutions and high-fidelity FEM references. In steady seepage examples, the proposed polygonal CSFEM reproduces linear hydraulic-head fields to near machine precision and yields smaller head errors than conventional FEM at the same characteristic mesh size. In transient problems, accurate head evolution and stable time responses are obtained, while adaptive refinement efficiently resolves localized high-gradient zones. For free-surface cases, the method captures the phreatic-surface profile and seepage-face development reliably without remeshing. The quadtree refinement and adaptivity provide substantial efficiency gains in degrees of freedom and runtime for a prescribed accuracy level.
title Two-dimensional seepage analysis using a polygonal cell-based smoothed finite element method
topic Numerical Analysis
65N30, 74S05
F.2.2; I.2.7
url https://arxiv.org/abs/2601.19357