76 research outputs found
A Virtual Element Method for a Nonlocal FitzHugh-Nagumo Model of Cardiac Electrophysiology
We present a Virtual Element Method (VEM) for a nonlocal reaction-diffusion
system of the cardiac electric field. To this system, we analyze an
-conforming discretization by means of VEM which can make use of
general polygonal meshes. Under standard assumptions on the computational
domain, we establish the convergence of the discrete solution by considering a
series of a priori estimates and by using a general compactness
criterion. Moreover, we obtain optimal order space-time error estimates in the
norm. Finally, we report some numerical tests supporting the theoretical
results
Mimetic Finite Difference methods for Hamiltonian wave equations in 2D
In this paper we consider the numerical solution of the Hamiltonian wave
equation in two spatial dimension. We use the Mimetic Finite Difference (MFD)
method to approximate the continuous problem combined with a symplectic
integration in time to integrate the semi-discrete Hamiltonian system. The main
characteristic of MFD methods, when applied to stationary problems, is to mimic
important properties of the continuous system. This approach, associated with a
symplectic method for the time integration yields a full numerical procedure
suitable to integrate Hamiltonian problems. A complete theoretical analysis of
the method and some numerical simulations are developed in the paper.Comment: 26 pages, 8 figure
Continuous finite element subgrid basis functions for Discontinuous Galerkin schemes on unstructured polygonal Voronoi meshes
We propose a new high order accurate nodal discontinuous Galerkin (DG) method
for the solution of nonlinear hyperbolic systems of partial differential
equations (PDE) on unstructured polygonal Voronoi meshes. Rather than using
classical polynomials of degree N inside each element, in our new approach the
discrete solution is represented by piecewise continuous polynomials of degree
N within each Voronoi element, using a continuous finite element basis defined
on a subgrid inside each polygon. We call the resulting subgrid basis an
agglomerated finite element (AFE) basis for the DG method on general polygons,
since it is obtained by the agglomeration of the finite element basis functions
associated with the subgrid triangles. The basis functions on each sub-triangle
are defined, as usual, on a universal reference element, hence allowing to
compute universal mass, flux and stiffness matrices for the subgrid triangles
once and for all in a pre-processing stage for the reference element only.
Consequently, the construction of an efficient quadrature-free algorithm is
possible, despite the unstructured nature of the computational grid. High order
of accuracy in time is achieved thanks to the ADER approach, making use of an
element-local space-time Galerkin finite element predictor.
The novel schemes are carefully validated against a set of typical benchmark
problems for the compressible Euler and Navier-Stokes equations. The numerical
results have been checked with reference solutions available in literature and
also systematically compared, in terms of computational efficiency and
accuracy, with those obtained by the corresponding modal DG version of the
scheme
The conforming virtual element method for polyharmonic and elastodynamics problems: a review
In this paper, we review recent results on the conforming virtual element
approximation of polyharmonic and elastodynamics problems. The structure and
the content of this review is motivated by three paradigmatic examples of
applications: classical and anisotropic Cahn-Hilliard equation and phase field
models for brittle fracture, that are briefly discussed in the first part of
the paper. We present and discuss the mathematical details of the conforming
virtual element approximation of linear polyharmonic problems, the classical
Cahn-Hilliard equation and linear elastodynamics problems.Comment: 30 pages, 7 figures. arXiv admin note: text overlap with
arXiv:1912.0712
Bulk-surface virtual element method for systems of PDEs in two-space dimensions
none3siIn this paper we consider a coupled bulk-surface PDE in two space dimensions. The model consists of a PDE in the bulk that is coupled to another PDE on the surface through general nonlinear boundary conditions. For such a system we propose a novel method, based on coupling a virtual element method (Beirão Da Veiga et al. in Math Models Methods Appl Sci 23(01):199–214, 2013. https://doi.org/10.1051/m2an/2013138) in the bulk domain to a surface finite element method (Dziuk and Elliott in Acta Numer 22:289–396, 2013. https://doi.org/10.1017/s0962492913000056) on the surface. The proposed method, which we coin the bulk-surface virtual element method includes, as a special case, the bulk-surface finite element method (BSFEM) on triangular meshes (Madzvamuse and Chung in Finite Elem Anal Des 108:9–21, 2016. https://doi.org/10.1016/j.finel.2015.09.002). The method exhibits second-order convergence in space, provided the exact solution is H2 + 1 / 4 in the bulk and H2 on the surface, where the additional 14 is required only in the simultaneous presence of surface curvature and non-triangular elements. Two novel techniques introduced in our analysis are (i) an L2-preserving inverse trace operator for the analysis of boundary conditions and (ii) the Sobolev extension as a replacement of the lifting operator (Elliott and Ranner in IMA J Num Anal 33(2):377–402, 2013. https://doi.org/10.1093/imanum/drs022) for sufficiently smooth exact solutions. The generality of the polygonal mesh can be exploited to optimize the computational time of matrix assembly. The method takes an optimised matrix-vector form that also simplifies the known special case of BSFEM on triangular meshes (Madzvamuse and Chung 2016). Three numerical examples illustrate our findings.Funding: Regione Puglia REFIN: Codice 901D2CAA UNISAL026openFrittelli M.; Madzvamuse A.; Sgura I.Frittelli, M.; Madzvamuse, A.; Sgura, I
- …