Related papers: Programming of Finite Element Methods in MATLAB
The aim of this work is to present the details of the finite element approach we developed for solving the Landau-Lifschitz-Gilbert equations in order to be able to treat problems involving complex geometries. There are several…
We consider scalar semilinear elliptic PDEs, where the nonlinearity is strongly monotone, but only locally Lipschitz continuous. To linearize the arising discrete nonlinear problem, we employ a damped Zarantonello iteration, which leads to…
We develop an efficient $hp$-finite element method for piecewise-smooth differential equations with periodic boundary conditions, using orthogonal polynomials defined on circular arcs. The operators derived from this basis are banded and…
We present a simple and efficient MATLAB implementation of the local refinement for polygonal meshes. The purpose of this implementation is primarily educational, especially on the teaching of adaptive virtual element methods.
We present here a Finite Element Method devoted to the simulation of 3D periodic structures of arbitrary geometry. The numerical method based on ARPACK and PARDISO libraries, is discussed with the aim of extracting the eigenmodes of…
We study a space-time finite element approach for the nonhomogeneous wave equation using a continuous time Galerkin method. We present fully implicit examples in 1+1, 2+1, and 3+1 dimensions using linear quadrilateral, hexahedral, and…
We focus here on a class of fourth-order parabolic equations that can be written as a system of second-order equations by introducing an auxiliary variable. We design a novel second-order fully discrete mixed finite element method to…
A macroscopic model to describe the dynamics of ion transport in ion channels is the Poisson-Nernst-Planck(PNP) equations. In this paper, we develop a finite-difference method for solving PNP equations, which is second-order accurate in…
We analyze the bit complexity of efficient algorithms for fundamental optimization problems, such as linear regression, $p$-norm regression, and linear programming (LP). State-of-the-art algorithms are iterative, and in terms of the number…
We consider families of finite elements on polygonal meshes, that are defined implicitly on each mesh cell as solutions of local Poisson problems with polynomial data. Functions in the local space on each mesh cell are evaluated via…
In this article, we establish a class of new accelerated modulus-based iteration methods for solving the linear complementarity problem. When the system matrix is an $H_+$-matrix, we present appropriate criteria for the convergence…
We use the work of Milton, Seppecher, and Bouchitt\'{e} on variational principles for waves in lossy media to formulate a finite element method for solving the complex Helmholtz equation that is based entirely on minimization. In…
We propose a novel iterative algorithm for solving a large sparse linear system. The method is based on the EM algorithm. If the system has a unique solution, the algorithm guarantees convergence with a geometric rate. Otherwise,…
This paper proposes a finite element method that couples mixed and Lagrange finite elements to efficiently capture stress concentrations in elasticity problems. The method employs conforming mixed finite elements in regions with stress…
The efficient solution of large sparse saddle point systems is very important in computational fluid mechanics. The discontinuous Galerkin finite element methods have become increasingly popular for incompressible flow problems but their…
We analyze the spatially semidiscrete piecewise linear finite element method for a nonlocal parabolic equation resulting from thermistor problem. Our approach is based on the properties of the elliptic projection defined by the bilinear…
The multimesh finite element method enables the solution of partial differential equations on a computational mesh composed by multiple arbitrarily overlapping meshes. The discretization is based on a continuous--discontinuous function…
A numerical scheme is presented for approximating fractional order Poisson problems in two and three dimensions. The scheme is based on reformulating the original problem posed over $\Omega$ on the extruded domain…
This article presents a $P_0$ finite element method for boundary value problems for linear elasticity equations. The new method makes use of piecewise constant approximating functions on the boundary of each polytopal element, and is…
The presented article contains a 2D mesh generation routine optimized with the Metropolis algorithm. The procedure enables to produce meshes with a prescribed size h of elements. These finite element meshes can serve as standard discrete…