4.7 Article

The lowest-order weak Galerkin finite element method for the Darcy equation on quadrilateral and hybrid meshes

Journal

JOURNAL OF COMPUTATIONAL PHYSICS
Volume 359, Issue -, Pages 312-330

Publisher

ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2018.01.001

Keywords

Darcy flow; Hybrid meshes; Lowest order finite elements; Quadrilateral meshes; Weak Galerkin

Funding

  1. US National Science Foundation (NSF) [DMS-1419077]
  2. NSF [DMS-1419077, DMS-1720473]
  3. Center of Interdisciplinary Mathematics and Statistics at Colorado State University
  4. Direct For Mathematical & Physical Scien
  5. Division Of Mathematical Sciences [1419077] Funding Source: National Science Foundation

Ask authors/readers for more resources

This paper investigates the lowest-order weak Galerkin finite element method for solving the Darcy equation on quadrilateral and hybrid meshes consisting of quadrilaterals and triangles. In this approach, the pressure is approximated by constants in element interiors and on edges. The discrete weak gradients of these constant basis functions are specified in local Raviart-Thomas spaces, specifically RT0 for triangles and unmapped RT[0] for quadrilaterals. These discrete weak gradients are used to approximate the classical gradient when solving the Darcy equation. The method produces continuous normal fluxes and is locally mass-conservative, regardless of mesh quality, and has optimal order convergence in pressure, velocity, and normal flux, when the quadrilaterals are asymptotically parallelograms. Implementation is straightforward and results in symmetric positive-definite discrete linear systems. We present numerical experiments and comparisons with other existing methods. (C) 2018 Elsevier Inc. All rights reserved.

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.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Mathematics, Applied

A lowest-order weak Galerkin finite element method for Stokes flow on polygonal meshes

Jiangguo Liu, Graham Harper, Nolisa Malluwawadu, Simon Tavener

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS (2020)

Article Mathematics, Applied

An efficient numerical solver for anisotropic subdiffusion problems

Jinying Tan, Jiangguo Liu

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS (2020)

Article Computer Science, Software Engineering

Error estimation and uncertainty quantification for first time to a threshold value

Jehanzeb H. Chaudhry, Donald Estep, Zachary Stevens, Simon J. Tavener

Summary: This study investigates the posterior error estimates of a quantity that cannot be represented as a linear functional in differential equations, providing two representations and estimating unknown terms through an adjoint-based approach. The combination of these representations yields accurate results.

BIT NUMERICAL MATHEMATICS (2021)

Article Mathematics, Applied

Penalty-Free Any-Order Weak Galerkin FEMs for Elliptic Problems on Quadrilateral Meshes

Jiangguo Liu, Simon Tavener, Zhuoran Wang

JOURNAL OF SCIENTIFIC COMPUTING (2020)

Article Mathematics, Applied

A locking-free solver for linear elasticity on quadrilateral and hexahedral meshes based on enrichment of Lagrangian elements

Graham Harper, Ruishu Wang, Jiangguo Liu, Simon Tavener, Ran Zhang

COMPUTERS & MATHEMATICS WITH APPLICATIONS (2020)

Article Engineering, Multidisciplinary

Coupling Arbogast-Correa and Bernardi-Raugel elements to resolve coupled Stokes-Darcy flow problems

Graham Harper, Jiangguo Liu, Simon Tavener, Tim Wildey

Summary: This paper presents a finite element method for solving coupled Stokes-Darcy flow problems by combining classical Bernardi-Raugel finite elements and the recently developed Arbogast-Correa (AC) spaces on quadrilateral meshes. A novel weak Galerkin methodology is employed for discretization of the Darcy equation, with piecewise constant approximants used to approximate the Darcy pressure. Rigorous error analysis and numerical experiments demonstrate the stability and optimal-order accuracy of the method.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2021)

Article Mathematics, Applied

Analysis of a 2-field finite element solver for poroelasticity on quadrilateral meshes

Zhuoran Wang, Simon Tavener, Jiangguo Liu

Summary: This paper introduces a novel 2-field finite element solver for linear poroelasticity on convex quadrilateral meshes. It discretizes Darcy flow for fluid pressure and linear elasticity for solid displacement, coupling them through implicit Euler temporal discretization to solve poroelasticity problems. The rigorous error analysis and numerical tests demonstrate the accuracy and locking-free property of this new solver.

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS (2021)

Article Biology

A Monte Carlo method to estimate cell population heterogeneity from cell snapshot data

Ben Lambert, David J. Gavaghan, Simon J. Tavener

Summary: Variation among cells in biological systems is common, and mathematical models can be used to investigate these differences. Researchers have introduced a computational sampling method called Contour Monte Carlo for estimating mathematical model parameters from snapshot distributions. This method is suitable for underdetermined systems and can help researchers explore cellular variation in biological systems.

JOURNAL OF THEORETICAL BIOLOGY (2021)

Article Computer Science, Software Engineering

A posteriori error analysis for Schwarz overlapping domain decomposition methods

Jehanzeb H. Chaudhry, Donald Estep, Simon J. Tavener

Summary: The study focuses on adjoint-based a posteriori error analysis for domain decomposition methods, decomposing numerical errors into different contributions and constructing a two-stage solution strategy based on this decomposition to efficiently reduce errors in the quantity of interest by adjusting the relative contributions.

BIT NUMERICAL MATHEMATICS (2021)

Article Mathematics, Applied

Robust weak Galerkin finite element solvers for Stokes flow based on a lifting operator

Zhuoran Wang, Ruishu Wang, Jiangguo Liu

Summary: This paper presents novel finite element solvers for solving Stokes flow, which are pressure-robust due to the use of a lifting operator. The solvers utilize weak Galerkin finite element schemes and local Arbogast-Correa or Arbogast-Tao spaces for construction of discrete weak gradients. The lifting operator helps in removing pressure dependence of velocity errors, ensuring the robustness of the solvers. The theoretical validation and numerical illustrations demonstrate the pressure robustness of the solvers. A comparison with the non-robust classical Taylor-Hood solver is also provided.

COMPUTERS & MATHEMATICS WITH APPLICATIONS (2022)

Article Mathematics, Applied

Penalty-Free Any-Order Weak Galerkin FEMs for Linear Elasticity on Quadrilateral Meshes

Ruishu Wang, Zhuoran Wang, Jiangguo Liu

Summary: This paper presents a family of new weak Galerkin finite element methods for solving linear elasticity in the primal formulation. These methods use vector-valued polynomials of degree k >= 0 to approximate the displacement independently in element interiors and on edges of a convex quadrilateral mesh. The new methods do not require penalty or stabilizer and are free of Poisson-locking, while achieving optimal order (k + 1) convergence rates in displacement, stress, and dilation. Numerical experiments on popular test cases demonstrate the theoretical estimates and efficiency of these new solvers. The extension to cuboidal hexahedral meshes is briefly discussed.

JOURNAL OF SCIENTIFIC COMPUTING (2023)

Article Mathematics, Applied

LOCKING-FREE CG-TYPE FINITE ELEMENT SOLVERS FOR LINEAR ELASTICITY ON SIMPLICIAL MESHES

Ruishu Wang, Zhuoran Wang, Jiangguo Liu, Simon Tavener, Ran Zhang

Summary: This paper presents numerical methods for solving linear elasticity on simplicial meshes based on enrichment of Lagrangian bilinear/trilinear finite elements. The classical 1st order Bernardi-Raugel spaces are innovatively used for this purpose, with a projection to the elementwise constant space to handle dilation in the strain-div formulation. Mixed boundary conditions are considered for error estimates in energy-norm and L-2-norms of displacement and stress, demonstrating methods free of Poisson-locking.

INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING (2021)

Article Computer Science, Interdisciplinary Applications

A new type of non-polynomial based TENO scheme for hyperbolic conservation laws

Tian Liang, Lin Fu

Summary: In this work, a new shock-capturing framework is proposed based on a new candidate stencil arrangement and the combination of infinitely differentiable non-polynomial RBF-based reconstruction in smooth regions with jump-like non-polynomial interpolation for genuine discontinuities. The resulting scheme achieves high order accuracy and resolves genuine discontinuities with sub-cell resolution.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A high-order residual-based viscosity finite element method for incompressible variable density flow

Lukas Lundgren, Murtazo Nazarov

Summary: In this paper, a high-order accurate finite element method for incompressible variable density flow is introduced. The method addresses the issues of saddle point system and stability problem through Schur complement preconditioning and artificial compressibility approaches, and it is validated to have high-order accuracy for smooth problems and accurately resolve discontinuities.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Convergence analysis and optimization of a Robin Schwarz waveform relaxation method for time-periodic parabolic optimal control problems

Gabriele Ciaramella, Laurence Halpern, Luca Mechelli

Summary: This paper presents a novel convergence analysis of the optimized Schwarz waveform relaxation method for solving optimal control problems governed by periodic parabolic PDEs. The analysis is based on a Fourier-type technique applied to a semidiscrete-in-time form of the optimality condition, which enables a precise characterization of the convergence factor at the semidiscrete level. The behavior of the optimal transmission condition parameter is also analyzed in detail as the time discretization approaches zero.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Data-driven Whitney forms for structure-preserving control volume analysis

Jonas A. Actor, Xiaozhe Hu, Andy Huang, Scott A. Roberts, Nathaniel Trask

Summary: This article introduces a scientific machine learning framework that uses a partition of unity architecture to model physics through control volume analysis. The framework can extract reduced models from full field data while preserving the physics. It is applicable to manifolds in arbitrary dimension and has been demonstrated effective in specific problems.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Higher-continuity s-version of finite element method with B-spline functions

Nozomi Magome, Naoki Morita, Shigeki Kaneko, Naoto Mitsume

Summary: This paper proposes a novel strategy called B-spline based SFEM to fundamentally solve the problems of the conventional SFEM. It uses different basis functions and cubic B-spline basis functions with C-2-continuity to improve the accuracy of numerical integration and avoid matrix singularity. Numerical results show that the proposed method is superior to conventional methods in terms of accuracy and convergence.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A cell-centred Eulerian volume-of-fluid method for compressible multi-material flows

Timothy R. Law, Philip T. Barton

Summary: This paper presents a practical cell-centred volume-of-fluid method for simulating compressible solid-fluid problems within a pure Eulerian setting. The method incorporates a mixed-cell update to maintain sharp interfaces, and can be easily extended to include other coupled physics. Various challenging test problems are used to validate the method, and its robustness and application in a multi-physics context are demonstrated.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Two-step multi-resolution reconstruction-based compact gas-kinetic scheme on tetrahedral mesh

Xing Ji, Fengxiang Zhao, Wei Shyy, Kun Xu

Summary: This paper presents the development of a third-order compact gas-kinetic scheme for compressible Euler and Navier-Stokes solutions, constructed particularly for an unstructured tetrahedral mesh. The scheme demonstrates robustness in high-speed flow computation and exhibits excellent adaptability to meshes with complex geometrical configurations.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Multiscale sampling for the inverse modeling of partial differential equations

Alsadig Ali, Abdullah Al-Mamun, Felipe Pereira, Arunasalam Rahunanthan

Summary: This paper presents a novel Bayesian statistical framework for the characterization of natural subsurface formations, and introduces the concept of multiscale sampling to localize the search in the stochastic space. The results show that the proposed framework performs well in solving inverse problems related to porous media flows.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Constrained optimized dynamic mode decomposition with control for physically stable systems with exogeneous inputs

Jacob Rains, Yi Wang, Alec House, Andrew L. Kaminsky, Nathan A. Tison, Vamshi M. Korivi

Summary: This paper presents a novel method called constrained optimized DMD with Control (cOptDMDc), which extends the optimized DMD method to systems with exogenous inputs and can enforce the stability of the resulting reduced order model (ROM). The proposed method optimally places eigenvalues within the stable region, thus mitigating spurious eigenvalue issues. Comparative studies show that cOptDMDc achieves high accuracy and robustness.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A hybridizable discontinuous Galerkin formulation for the Euler-Maxwell plasma model

Andrea La Spina, Jacob Fish

Summary: This work introduces a hybridizable discontinuous Galerkin formulation for simulating ideal plasmas. The proposed method couples the fluid and electromagnetic subproblems monolithically based on source and employs a fully implicit time integration scheme. The approach also utilizes a projection-based divergence correction method to enforce the Gauss laws in challenging scenarios. Numerical examples demonstrate the high-order accuracy, efficiency, and robustness of the proposed formulation.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Numerical solution of the cavity scattering problem for flexural waves on thin plates: Linear finite element methods

Junhong Yue, Peijun Li

Summary: This paper proposes two numerical methods (IP-FEM and BP-FEM) to study the flexural wave scattering problem of an arbitrary-shaped cavity on an infinite thin plate. These methods successfully decompose the fourth-order plate wave equation into the Helmholtz and modified Helmholtz equations with coupled conditions on the cavity boundary, providing an effective solution to this challenging problem.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Fast and scalable computation of shape-morphing nonlinear solutions with application to evolutional neural networks

William Anderson, Mohammad Farazmand

Summary: We develop fast and scalable methods, called RONS, for computing reduced-order nonlinear solutions. These methods have been proven to be highly effective in tackling challenging problems, but become computationally prohibitive as the number of parameters grows. To address this issue, three separate methods are proposed and their efficacy is demonstrated through examples. The application of RONS to neural networks is also discussed.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A second order directional split exponential integrator for systems of advection-diffusion-reaction equations

Marco Caliari, Fabio Cassini

Summary: In this paper, a second order exponential scheme for stiff evolutionary advection-diffusion-reaction equations is proposed. The scheme is based on a directional splitting approach and uses computation of small sized exponential-like functions and tensor-matrix products for efficient implementation. Numerical examples demonstrate the advantage of the proposed approach over state-of-the-art techniques.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

A conservative semi-Lagrangian method for inhomogeneous Boltzmann equation

Sebastiano Boscarino, Seung Yeon Cho, Giovanni Russo

Summary: This work proposes a high order conservative semi-Lagrangian method for the inhomogeneous Boltzmann equation of rarefied gas dynamics. The method combines a semi-Lagrangian scheme for the convection term, a fast spectral method for computation of the collision operator, and a high order conservative reconstruction and a weighted optimization technique to preserve conservative quantities. Numerical tests demonstrate the accuracy and efficiency of the proposed method.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)

Article Computer Science, Interdisciplinary Applications

Uniqueness and numerical scheme for spherical shell-structured sources from the far field patterns with at most two frequencies

Jialei Li, Xiaodong Liu, Qingxiang Shi

Summary: This study shows that the number, centers, scattering strengths, inner and outer diameters of spherical shell-structured sources can be uniquely determined from the far field patterns. A numerical scheme is proposed for reconstructing the spherical shell-structured sources, which includes a migration series method for locating the centers and an iterative method for computing the inner and outer diameters without computing derivatives.

JOURNAL OF COMPUTATIONAL PHYSICS (2024)