4.7 Article

Mimetic finite difference method for the Stokes problem on polygonal meshes

Journal

JOURNAL OF COMPUTATIONAL PHYSICS
Volume 228, Issue 19, Pages 7215-7232

Publisher

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

Keywords

Incompressible Stokes equations; Mimetic discretization; Polygonal mesh

Ask authors/readers for more resources

Various approaches to extend finite element methods to non-traditional elements (general polygons, pyramids, polyhedra, etc.) have been developed over the last decade. The construction of basis functions for such elements is a challenging task and may require extensive geometrical analysis. The mimetic finite difference (MFD) method works on general polygonal meshes and has many similarities with low-order finite element methods. Both schemes try to preserve the fundamental properties of the underlying physical and mathematical models. The essential difference between the two schemes is that the MFD method uses only the surface representation of discrete unknowns to build the stiffness and mass matrices. Since no extension of basis functions inside the mesh elements is required, practical implementation of the MFD method is simple for polygonal meshes that may include degenerate and non-convex elements. In this article, we present a new MFD method for the Stokes problem on arbitrary polygonal meshes and analyze its stability. The method is developed for the general case of tensor coefficients, which allows us to apply it to a linear elasticity problem, as well. Numerical experiments show, for the velocity variable, second-order convergence in a discrete L 2 norm and first-order convergence in a discrete H-1 norm. For the pressure variable, first-order convergence is shown in the L 2 norm. (C) 2009 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 Engineering, Multidisciplinary

Curvilinear Virtual Elements for 2D solid mechanics applications

E. Artioli, L. Beirao da Veiga, F. Dassi

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2020)

Article Mathematics, Applied

The Stokes complex for Virtual Elements in three dimensions

L. Beirao da Veiga, F. Dassi, G. Vacca

MATHEMATICAL MODELS & METHODS IN APPLIED SCIENCES (2020)

Article Mathematics, Applied

Polynomial preserving virtual elements with curved edges

L. Beirao da Veiga, F. Brezzi, L. D. Marini, A. Russo

MATHEMATICAL MODELS & METHODS IN APPLIED SCIENCES (2020)

Review Mathematics, Applied

An adaptive curved virtual element method for the statistical homogenization of random fibre-reinforced composites

E. Artioli, L. Beirao da Veiga, M. Verani

FINITE ELEMENTS IN ANALYSIS AND DESIGN (2020)

Article Engineering, Multidisciplinary

The arbitrary-order virtual element method for linear elastodynamics models: convergence, stability and dispersion-dissipation analysis

Paola F. Antonietti, Gianmarco Manzini, Ilario Mazzieri, Hashem M. Mourad, Marco Verani

Summary: The study introduces the conforming virtual element method for numerical approximation of two-dimensional elastodynamics problem, proves stability and convergence of the method, and derives optimal error estimates under different refinements. Experimental results demonstrate the method's effectiveness on various computational meshes and show exponential convergence under p-refinement. Dispersion-dissipation analysis reveals that polygonal meshes exhibit similar properties to classical simplicial/quadrilateral grids.

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING (2021)

Article Engineering, Multidisciplinary

A virtual element method for the miscible displacement of incompressible fluids in porous media

L. Beirao da Veiga, A. Pichler, G. Vacca

Summary: This paper presents a virtual element (VE) discretization for a time-dependent coupled system of nonlinear partial differential equations, aiming to investigate the capabilities of virtual element methods (VEM) for complex fluid flow problems. By combining VEM with a time stepping scheme, a theoretical analysis of the method was developed under the assumption of a regular solution. The scheme was then tested on both regular and realistic test cases to validate its effectiveness.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2021)

Article Mathematics, Applied

Virtual elements for Maxwell's equations

L. Beirao da Veiga, F. Dassi, G. Manzini, L. Mascotto

Summary: We introduce a low order virtual element discretization for time dependent Maxwell's equations, which allows for the use of general polyhedral meshes. Both the semi- and fully-discrete schemes are considered. We derive optimal a priori estimates and validate them through numerical experiments. As key findings, we discuss novel inequalities associated with de Rahm sequences of nodal, edge, and face virtual element spaces.

COMPUTERS & MATHEMATICS WITH APPLICATIONS (2022)

Article Mathematics, Applied

Stability and Conservation Properties of Hermite-Based Approximations of the Vlasov-Poisson System

Daniele Funaro, Gianmarco Manzini

Summary: This study provides spectral approximation based on Hermite-Fourier expansion of the Vlasov-Poisson model with high-order artificial collision operators. The analysis considers the necessary conditions related to the artificial collision term, number of spectral modes, and time-step in order to ensure stability in appropriate norms. The study starts with a Hermite discretization of a simple linear problem in one dimension and extends partially to cover the complete nonlinear Vlasov-Poisson model.

JOURNAL OF SCIENTIFIC COMPUTING (2021)

Article Mathematics, Applied

Equilibrium analysis of an immersed rigid leaflet by the virtual element method

L. Beirao da Veiga, C. Canuto, R. H. Nochetto, G. Vacca

Summary: This study investigates the equilibrium of a hinged rigid leaflet with an attached rotational spring in a stationary incompressible fluid within a rigid channel using theoretical and numerical methods. Sufficient conditions for the existence and uniqueness of equilibrium positions are identified based on properties of the domain functional. The proposed numerical technique utilizes the mesh flexibility of the Virtual Element Method and proves quasi-optimal error estimates through a variety of numerical experiments.

MATHEMATICAL MODELS & METHODS IN APPLIED SCIENCES (2021)

Article Engineering, Multidisciplinary

Extended virtual element method for two-dimensional linear elastic fracture

E. Benvenuti, A. Chiozzi, G. Manzini, N. Sukumar

Summary: In this paper, the authors propose an eXtended Virtual Element Method (X-VEM) for two-dimensional linear elastic fracture. The X-VEM allows for mesh-independent modeling of crack discontinuities and elastic crack-tip singularities on general polygonal meshes. The method involves an extended projector and additional basis functions constructed from standard virtual basis functions and enrichment fields. Numerical experiments demonstrate the accuracy and optimal convergence of the X-VEM formulation for mixed-mode linear elastic fracture problems.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2022)

Article Mathematics, Applied

Conforming virtual element approximations of the two-dimensional Stokes problem

Gianmarco Manzini, Annamaria Mazzia

Summary: In this paper, two conforming virtual element formulations for the numerical approximation of the Stokes problem on polygonal meshes are presented. Both formulations are inf-sup stable and have optimal convergence rates in the L2 and energy norm. The effectiveness of these numerical approximations is assessed through investigation on a representative benchmark problem.

APPLIED NUMERICAL MATHEMATICS (2022)

Article Mathematics, Applied

Stability and Interpolation Properties for Stokes-Like Virtual Element Spaces

Jian Meng, Lourenco Beirao da Veiga, Lorenzo Mascotto

Summary: In this paper, we establish stability bounds for Stokes-like virtual element spaces in both two and three dimensions. These bounds are crucial for deriving optimal interpolation estimates. In addition, we conduct numerical tests to investigate the behavior of the stability constants from a practical perspective.

JOURNAL OF SCIENTIFIC COMPUTING (2023)

Article Engineering, Multidisciplinary

A fully-discrete virtual element method for the nonstationary Boussinesq equations in stream-function form

L. Beirao da Veiga, D. Mora, A. Silgado

Summary: In this paper, a fully-coupled virtual element method is proposed for solving the nonstationary Boussinesq system in 2D. The method utilizes the stream-function and temperature fields and employs C1- and C0-conforming virtual element approaches for spatial discretization. The temporal variable is discretized using a backward Euler scheme. The well-posedness and unconditional stability of the fully-discrete problem are proved, and error estimates in H2- and H1-norms are derived for the stream-function and temperature fields. Benchmark tests are conducted to validate the theoretical error bounds and demonstrate the behavior of the fully-discrete scheme.

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING (2023)

Article Computer Science, Interdisciplinary Applications

The tensor-train mimetic finite difference method for three-dimensional Maxwell?s wave propagation equations

G. Manzini, R. Vuchkov, B. Alexandrov

Summary: Coupling the mimetic finite difference method with the tensor-train format allows for efficient low-rank numerical approximations of the solutions of the time-dependent Maxwell wave propagation equations in three dimensions. By discretizing the curl operators on the primal/dual tensor product grid complex and coupling it with a staggered-in-time second-order accurate time-marching scheme, we obtain a solver that is accurate to the second order in time and space. The use of the tensor-train format significantly improves solver performance in terms of CPU time and memory storage.

MATHEMATICS AND COMPUTERS IN SIMULATION (2023)

Article Mathematics, Applied

A VIRTUAL ELEMENT GENERALIZATION ON POLYGONAL MESHES OF THE SCOTT-VOGELIUS FINITE ELEMENT METHOD FOR THE 2-D STOKES PROBLEM

Gianmarco Manzini, Annamaria Mazzia

Summary: This study presents a new method for the numerical approximation of the Stokes problem on polygonal meshes within the framework of the virtual element method. The numerical approximation shows optimal convergence rates in various cases, except for the lowest order on triangular meshes, and in square meshes known to be unstable.

JOURNAL OF COMPUTATIONAL DYNAMICS (2022)

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)