A new weak Galerkin (WG) method is introduced and analyzed for the second order elliptic equation formulated as a system of two first order linear equations. This method, called WG-MFEM, is designed ...by using discontinuous piecewise polynomials on finite element partitions with arbitrary shape of polygons/polyhedra. The WG-MFEM is capable of providing very accurate numerical approximations for both the primary and flux variables. Allowing the use of discontinuous approximating functions on arbitrary shape of polygons/polyhedra makes the method highly flexible in practical computation. Optimal order error estimates in both discrete H^1 and L^2 norms are established for the corresponding weak Galerkin mixed finite element solutions.
This paper introduces a finite element method by using a weakly defined gradient operator over generalized functions. The use of weak gradients and their approximations results in a new concept ...called discrete weak gradients which is expected to play an important role in numerical methods for partial differential equations. This article intends to provide a general framework for managing differential operators on generalized functions. As a demonstrative example, the discrete weak gradient operator is employed as a building block in the design of numerical schemes for a second order elliptic problem, in which the classical gradient operator is replaced by the discrete weak gradient. The resulting numerical scheme is called a weak Galerkin (WG) finite element method. It can be seen that the weak Galerkin method allows the use of totally discontinuous functions in the finite element procedure. For the second order elliptic problem, an optimal order error estimate in both a discrete H1 and L2 norms are established for the corresponding weak Galerkin finite element solutions. A superconvergence is also observed for the weak Galerkin approximation.
MFEM: A modular finite element methods library Anderson, Robert; Andrej, Julian; Barker, Andrew ...
Computers & mathematics with applications (1987),
01/2021, Volume:
81, Issue:
1
Journal Article
Peer reviewed
Open access
MFEM is an open-source, lightweight, flexible and scalable C++ library for modular finite element methods that features arbitrary high-order finite element meshes and spaces, support for a wide ...variety of discretization approaches and emphasis on usability, portability, and high-performance computing efficiency. MFEM’s goal is to provide application scientists with access to cutting-edge algorithms for high-order finite element meshing, discretizations and linear solvers, while enabling researchers to quickly and easily develop and test new algorithms in very general, fully unstructured, high-order, parallel and GPU-accelerated settings. In this paper we describe the underlying algorithms and finite element abstractions provided by MFEM, discuss the software implementation, and illustrate various applications of the library.
Abstract
In the context of unfitted finite element discretizations, the realization of high-order methods is challenging due to the fact that the geometry approximation has to be sufficiently ...accurate. We consider a new unfitted finite element method that achieves a high-order approximation of the geometry for domains that are implicitly described by smooth-level set functions. The method is based on a parametric mapping, which transforms a piecewise planar interface reconstruction to a high-order approximation. Both components, the piecewise planar interface reconstruction and the parametric mapping, are easy to implement. In this article, we present an a priori error analysis of the method applied to an interface problem. The analysis reveals optimal order error bounds for the geometry approximation and for the finite element approximation, for arbitrary high-order discretization. The theoretical results are confirmed in numerical experiments.
A preasymptotic error analysis of the finite element method (FEM) and some continuous interior penalty finite element method (CIP-FEM) for the Helmholtz equation in two and three dimensions is ...proposed. H1- and L2-error estimates with explicit dependence on the wave number k are derived. In particular, it is shown that if k2p+1 h2p is sufficiently small, then the pollution errors of both methods in H1-norm are bounded by O(k2p+1 h2p), which coincides with the phase error of the FEM obtained by existent dispersion analyses on Cartesian grids, where h is the mesh size, and p is the order of the approximation space and is fixed. The CIP-FEM extends the classical one by adding more penalty terms on jumps of higher (up to pth order) normal derivatives in order to reduce efficiently the pollution errors of higher order methods. Numerical tests are provided to verify the theoretical findings and to illustrate the great capability of the CIP-FEM in reducing the pollution effect.
In the present contribution, a phenomenological constitutive model of medium manganese steels, in which both twinning-induced (TWIP) and transformation-induced (TRIP) plasticity enhancing mechanisms ...are activated, is implemented in the finite element framework. The implementation is utilized for the analysis of the full-field strain partitioning in dual-phase microstructure maps obtained from electron backscattering diffraction. The results of the finite element analysis suggest that the strain localization in the studied steel has an alternating character. Specifically, in the low strain region, most of the externally imposed deformation is accommodated by the initially softer austenite. The higher strain hardening rate of austenite due to deformation twinning (TWIP effect) and the mechanically-induced transformation to martensite (TRIP effect) results in a shift of the strain localization to ferrite. This alternating strain localization is a key feature that distinguishes the medium manganese TWIP+TRIP steel. It is shown that this alternating strain localization contributes to the superior mechanical behavior of medium manganese TWIP+TRIP steel reported in the literature.
Display omitted
Extension operators for trimmed spline spaces Burman, Erik; Hansbo, Peter; Larson, Mats G. ...
Computer methods in applied mechanics and engineering,
01/2023, Volume:
403
Journal Article
Peer reviewed
Open access
We develop a discrete extension operator for trimmed spline spaces consisting of piecewise polynomial functions of degree p with k continuous derivatives. The construction is based on polynomial ...extension from neighboring elements together with projection back into the spline space. We prove stability and approximation results for the extension operator. Finally, we illustrate how we can use the extension operator to construct a stable cut isogeometric method for an elliptic model problem.
Nektar++ is an open-source framework that provides a flexible, high-performance and scalable platform for the development of solvers for partial differential equations using the high-order ...spectral/hp element method. In particular, Nektar++ aims to overcome the complex implementation challenges that are often associated with high-order methods, thereby allowing them to be more readily used in a wide range of application areas. In this paper, we present the algorithmic, implementation and application developments associated with our Nektar++ version 5.0 release. We describe some of the key software and performance developments, including our strategies on parallel I/O, on in situ processing, the use of collective operations for exploiting current and emerging hardware, and interfaces to enable multi-solver coupling. Furthermore, we provide details on a newly developed Python interface that enables a more rapid introduction for new users unfamiliar with spectral/hp element methods, C++ and/or Nektar++. This release also incorporates a number of numerical method developments – in particular: the method of moving frames (MMF), which provides an additional approach for the simulation of equations on embedded curvilinear manifolds and domains; a means of handling spatially variable polynomial order; and a novel technique for quasi-3D simulations (which combine a 2D spectral element and 1D Fourier spectral method) to permit spatially-varying perturbations to the geometry in the homogeneous direction. Finally, we demonstrate the new application-level features provided in this release, namely: a facility for generating high-order curvilinear meshes called NekMesh; a novel new AcousticSolver for aeroacoustic problems; our development of a ‘thick’ strip model for the modelling of fluid–structure interaction (FSI) problems in the context of vortex-induced vibrations (VIV). We conclude by commenting on some lessons learned and by discussing some directions for future code development and expansion.
Program Title: Nektar++
Program Files doi:http://dx.doi.org/10.17632/9drxd9d8nx.1
Code Ocean Capsule:https://doi.org/10.24433/CO.9865757.v1
Licensing provisions: MIT
Programming language: C++
External routines/libraries: Boost, METIS, FFTW, MPI, Scotch, PETSc, TinyXML, HDF5, OpenCASCADE, CWIPI
Nature of problem: The Nektar++ framework is designed to enable the discretisation and solution of time-independent or time-dependent partial differential equations.
Solution method: spectral/hp element method
PDF
