varFEM: variational formulation based programming for finite element methods in Matlab

Yue Yu

DOI: https://doi.org/10.48550/arXiv.2206.06918

2022-06-14

Abstract:This paper summarizes the development of varFEM, which provides a realization of the programming style in FreeFEM by using the Matlab language.

Numerical Analysis,Mathematical Software

What problem does this paper attempt to address?

Object Oriented Finite Element Method and Its Application in C++

王春江,钱若军,董石麟,赵金城

DOI: https://doi.org/10.3321/j.issn:1006-2467.2004.06.024

2004-01-01

Abstract:This paper totally concluded the basic theory of object oriented finite element method (OOFEM) and its realization procedure, and also compared the same and the different between the traditional finite element method (FEM) and the OOFEM. For spatial structures analysis software it presented the system architecture by using unified modeling language (UML) and the software programs by using the object-oriented language C++, which proves that the OOFEM is advanced and the software architecture written by C++ is applicable.
A General Finite Element Method: Extension of Variational Analysis for Nonlinear Heat Conduction with Temperature-Dependent Properties and Boundary Conditions, and Its Implementation As Local Refinement.

Xin Yao,Yihe Wang,Jianxing Leng

DOI: https://doi.org/10.1016/j.camwa.2021.08.024

IF: 3.218

2021-01-01

Computers & Mathematics with Applications

Abstract:In simulation of heat conduction with temperature-independent physical properties and boundary conditions (BCs), Galerkin residual analysis and variational analysis yield equivalent finite element method (FEM), the conventional FEM. However, if the properties and BCs are temperature-dependent, it is discovered that their derivatives further induce nonlinearity of FEM which consequently generates divergence between the two analyses. A general FEM, extension of variational analysis, is derived as general form of conventional FEM modeling nonlinear heat conduction. Numerical examples demonstrate that the general FEM produces results with considerably higher accuracy and stability and also possesses higher performances on conforming with both two analyses. Since general FEM degenerates to conventional FEM if derivatives are of small-amplitude or zero and its direct implementation to the entire domain is costive, general FEM is alternatively utilized as local refinement of governing equation only to points with significant derivatives. The strategy of local refinement is optimized to enhance efficiency.
Vector Form of Intrinsic Finite Element Method for Incompressible Fluids

Akram Samy,Shu Li,Xingfei Yuan,Chengwei Liu,Yongcan Dong

DOI: https://doi.org/10.1016/j.compfluid.2024.106319

IF: 3.077

2024-01-01

Computers & Fluids

Abstract:Vector form of intrinsic finite element (VFIFE) is a numerical method widely used in solid mechanics. However, it's hard to extend the VFIFE method to fluid mechanics since the traditional VFIFE method fails to reflect the analytical equilibrium of multiple variables in the continuum. Therefore, under the framework of analytical mechanics, this paper proposes Lagrange's equation of the second kind in fluid mechanics with the extremum condition of Lagrange power functional. And a vectorized motion equation of incompressible viscous fluids is deduced from Lagrange's equation. By using several efficient algorithms in the finite difference method (FDM) and the finite element method (FEM), the NS equation is decomposed into four governing equations of vector form for fluid mechanics. In addition, with the application of the classic Smagorinsky sub-grid scale model in large eddy simulation (LES), this paper puts forward turbulence modelling with VFIFE procedure, and a corresponding MATLAB program is developed. Two typical examples are given to demonstrate the applicability and efficiency of the proposed large eddy simulation with VFIFE method. The proposed algorithm can effectively eliminate the non-physical oscillation of the pressure, and obtain much accurate results with a small number of grids.
A general finite element method

Xin Yao,Yihe Wang,Jianxing Leng

DOI: https://doi.org/10.1016/j.camwa.2021.08.024

2021-01-01

Abstract:• A general FEM for nonlinear heat conduction with temperature-dependent properties and boundary conditions is developed. • The internal nonlinearity of FEM induced by derivatives of temperature-dependent terms is revealed in general FEM. • The induced nonlinearity terminates the equivalence of Galerkin residual and variational analyses. • The general FEM is adopted as local refinement of governing equation in pursuit for efficiency. In simulation of heat conduction with temperature-independent physical properties and boundary conditions (BCs), Galerkin residual analysis and variational analysis yield equivalent finite element method (FEM), the conventional FEM. However, if the properties and BCs are temperature-dependent, it is discovered that their derivatives further induce nonlinearity of FEM which consequently generates divergence between the two analyses. A general FEM, extension of variational analysis, is derived as general form of conventional FEM modeling nonlinear heat conduction. Numerical examples demonstrate that the general FEM produces results with considerably higher accuracy and stability and also possesses higher performances on conforming with both two analyses. Since general FEM degenerates to conventional FEM if derivatives are of small-amplitude or zero and its direct implementation to the entire domain is costive, general FEM is alternatively utilized as local refinement of governing equation only to points with significant derivatives. The strategy of local refinement is optimized to enhance efficiency.
Unified Implementation of Finite Element Methods involving Jumps and Averages in Matlab

Yue Yu

DOI: https://doi.org/10.48550/arXiv.2209.10800

2022-12-12

Abstract:We provide unified implementations of the finite element methods involving jumps and averages in Maltab by combing the use of the software package varFEM, including the adaptive finite element methods for the Poisson equation and the $C^0$ interior penalty methods for the biharmonic equation. The design ideas can be extended to other Galerkin-based methods, for example, the discontinuous Galerkin methods and the virtual element methods.

Numerical Analysis
Programming of Finite Element Methods in MATLAB

Long Chen

DOI: https://doi.org/10.48550/arXiv.1804.05156

2018-04-14

Abstract:We discuss how to implement the linear finite element method for solving the Poisson equation. We begin with the data structure to represent the triangulation and boundary conditions, introduce the sparse matrix, and then discuss the assembling process. We pay special attention to an efficient programming style using sparse matrices in MATLAB.

Numerical Analysis
Programming of linear virtual element methods in three dimensions

Yue Yu

DOI: https://doi.org/10.48550/arXiv.2202.08951

2022-02-18

Abstract:We present a simple and efficient MATLAB implementation of the linear virtual element method for the three dimensional Poisson equation. The purpose of this software is primarily educational, to demonstrate how the key components of the method can be translated into code.

Numerical Analysis
Implementation of Object-oriented Program for Local P-version Adaptive FEM

夏小和,王建华,陈震东,沈为平

DOI: https://doi.org/10.3969/j.issn.0254-0053.2002.04.004

2002-01-01

Abstract:In this paper, C++ langrage with Object-oriented features was used design a P-version adaptive FEM program, which contains the modules of high-order and Serendipity shape function, high-order element, adaptive course and error estimate. In the end, besides a series of class and their functions, the input and output modules were constructed, which make the program become whole adaptive FEM software. Examples show that the errors decrease with the increasing of p refinement. The local P-version a-daptive zone is upgraded concentrate in stress-singularity zone. The abouve results show that the formulate derivation and FEM program implementation are correct.
Some researches on finite element method

Da-wei YANG,Zong-lin WANG,Xiu-hui JIN

DOI: https://doi.org/10.3321/j.issn:0367-6234.2006.02.020

2006-01-01

Abstract:After reconstructing with an object-oriented technique, traditional FEM programme's readability, maintainability, and extendibility are improved. This paper presents and resolves some problems in object-oriented FEM programs, including FEM program organization, code style, dynamic identification and creation of element object, as well as two dimensional various band width of array scheme of general stiffness. These works indicate that in order to achieve the final aim of object-oriented FEM programming, it is necessary to consider more profoundly in both general terms and details in object-oriented views.
On a vectorized basic linear algebra package for prototyping codes in MATLAB

Alexej Moskovka,Talal Rahman,Jan Valdman,Jon Eivind Vatne

2024-03-16

Abstract:When writing high-performance code for numerical computation in a scripting language like MATLAB, it is crucial to have the operations in a large for-loop vectorized. If not, the code becomes too slow to use, even for a moderately large problem. However, in the process of vectorizing, the code often loses its original structure and becomes less readable. This is particularly true in the case of a finite element implementation, even though finite element methods are inherently structured. A basic remedy to this is the separation of the vectorization part from the mathematics part of the code, which is easily achieved through building the code on top of the basic linear algebra subprograms that are already vectorized codes, an idea that has been used in a series of papers over the last fifteen years, developing codes that are fast and still structured and readable. We discuss the vectorized basic linear algebra package and introduce a formalism using multi-linear algebra to explain and define formally the functions in the package, as well as MATLAB pagetime functions. We provide examples from computations of varying complexity, including the computation of normal vectors, volumes, and finite element methods. Benchmarking shows that we also get fast computations. Using the library, we can write codes that closely follow our mathematical thinking, making writing, following, reusing, and extending the code easier.

Mathematical Software
Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements

Immanuel Anjam,Jan Valdman

DOI: https://doi.org/10.1016/j.amc.2015.03.105

2015-05-11

Abstract:We propose an effective and flexible way to assemble finite element stiffness and mass matrices in MATLAB. We apply this for problems discretized by edge finite elements. Typical edge finite elements are Raviart-Thomas elements used in discretizations of H(div) spaces and Nedelec elements in discretizations of H(curl) spaces. We explain vectorization ideas and comment on a freely available MATLAB code which is fast and scalable with respect to time.

Mathematical Software,Numerical Analysis
Efficient P1-FEM for Any Space Dimension in Matlab

Stefanie Beuter,Stefan A. Funken

DOI: https://doi.org/10.1515/cmam-2022-0239

2024-01-25

Computational Methods in Applied Mathematics

Abstract:This paper deals with the efficient implementation of the finite element method with continuous piecewise linear functions (P1-FEM) in R d ( d ∈ N ). Although at present there does not seem to be a very high practical demand for finite element methods that use higher-dimensional simplicial partitions, there are some advantages in studying the efficient implementation of the method independent of the dimension. For instance, it provides additional insights into necessary data structures and the complexity of implementations. Throughout, the focus is on an efficient realization using Matlab built-in functions and vectorization. The fast and vectorized Matlab function can be easily implemented in many other vector languages and is provided in Julia, too. The complete implementation of the adaptive FEM is given, including assembling stiffness matrix, building load vector, error estimation, and adaptive mesh-refinement. Numerical experiments underline the efficiency of our freely available code which is observed to be of a slightly more than linear complexity with respect to the number of elements when memory limits are not exceeded.

mathematics, applied
MFEM: A modular finite element methods library

Robert Anderson,Julian Andrej,Andrew Barker,Jamie Bramwell,Jean-Sylvain Camier,Jakub Cerveny,Veselin Dobrev,Yohann Dudouit,Aaron Fisher,Tzanio Kolev,Will Pazner,Mark Stowell,Vladimir Tomov,Ido Akkerman,Johann Dahm,David Medina,Stefano Zampini

DOI: https://doi.org/10.1016/j.camwa.2020.06.009

2021-01-01

Abstract: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.

mathematics, applied
Finite Element Analysis using MATLAB

Rahul Yadav

DOI: https://doi.org/10.22214/ijraset.2022.39902

2022-01-31

International Journal for Research in Applied Science and Engineering Technology

Abstract:Abstract: Finite Element Analysis or FEA is the simulation of a physical phenomenon using a numerical mathematic technique referred to as the Finite Element Method or FEM. This process is at the core of mechanical engineering, as well as a variety of other disciplines. It also is one of the key principles used in the development of simulation software. Engineers can use these FEM to reduce the number of physical prototypes and run virtual experiments to optimize their designs. Finite Element Analysis is used in various fields like structural, fluid flow, heat transfer to estimate the behavior of a component in real environment. There now exists growing body of knowledge connected with the development of mathematical models and numerical simulations of physical model. There are various software packages like Ansys, OptiStruct, COMSOL, Solidworks and many more which provide a close estimate in simulation models. Apart from these software MATLAB also has a Partial Differential Equation (PDE) toolbox which enables us to perform these simulations using some built-in functions and codes. However, a computational numerical technique is not an end of design rather it just provides a great estimate of the final component for which the results are only as good as the input provided and the final component in most of the cases require a physical testing in environment which it is meant to perform as a validation. Keywords: Brake pedal, MATLAB, Partial differential equation (PDE) toolbox, Tetrahedral mesh.

English Else
Automating the Finite Element Method

Anders Logg

DOI: https://doi.org/10.1007/s11831-007-9003-9

2011-12-02

Abstract:The finite element method can be viewed as a machine that automates the discretization of differential equations, taking as input a variational problem, a finite element and a mesh, and producing as output a system of discrete equations. However, the generality of the framework provided by the finite element method is seldom reflected in implementations (realizations), which are often specialized and can handle only a small set of variational problems and finite elements (but are typically parametrized over the choice of mesh). This paper reviews ongoing research in the direction of a complete automation of the finite element method. In particular, this work discusses algorithms for the efficient and automatic computation of a system of discrete equations from a given variational problem, finite element and mesh. It is demonstrated that by automatically generating and compiling efficient low-level code, it is possible to parametrize a finite element code over variational problem and finite element in addition to the mesh.

Numerical Analysis
nlfem: a Flexible 2d FEM Python Code for Nonlocal Convection-Diffusion and Mechanics

Manuel Klar,Christian Vollmann,Volker Schulz

DOI: https://doi.org/10.1007/s42102-023-00108-6

2023-09-26

Journal of Peridynamics and Nonlocal Modeling

Abstract:Abstract In this work, we present the mathematical foundation of an assembly code for finite element approximations of nonlocal models with compactly supported, weakly singular kernels. We demonstrate the code on a nonlocal diffusion model in various configurations and on a two-dimensional bond-based peridynamics model. Further examples can be found in D’Elia et al. (Math Models Methods Appl Sci 31(08):1505–1567, 2021). The code nlfem is published under the GNU General Public License (for details, see, e.g., https://www.gnu.org/licenses/gpl-3.0.de.html ) and can be freely downloaded at https://gitlab.uni-trier.de/pde-opt/nonlocal-models/nlfem .
Variational Principles for Rigid Finite Element Method

Xiong Zhang,Ming-wan Lu

1996-01-01

Abstract:Rigid Finite Element Method (RFEM) have been proved as a powerful method for discon- tinuum mechanics. As a kind of finite element method the variational principles for RFEM are presented and then the formulations of RFEM are derived from the minimum potential energy theory. On the basis of the variational principles, RFEM can be easily extended to other nonlin- ear problems such as elastic plastic analysis, limit analysis, etc.
Vectorized implementation of primal hybrid FEM in MATLAB

Harish Nagula Mallesham,Kamana Porwal,Jan Valdman,Sanjib Kumar Acharya

2024-01-28

Abstract:We present efficient MATLAB implementations of the lowest-order primal hybrid finite element method (FEM) for linear second-order elliptic and parabolic problems with mixed boundary conditions in two spatial dimensions. We employ the Crank-Nicolson finite difference scheme for the complete discrete setup of the parabolic problem. All the codes presented are fully vectorized using matrix-wise array operations. Numerical experiments are conducted to show the performance of the software.

Numerical Analysis
AutoFreeFem: Automatic code generation with FreeFEM++ and LaTex output for shape and topology optimization of non-linear multi-physics problems

Grégoire Allaire,Michael H. Gfrerer

2024-07-16

Abstract:For an educational purpose we develop the Python package AutoFreeFem which generates all ingredients for shape optimization with non-linear multi-physics in FreeFEM++ and also outputs the expressions for use in LaTex. As an input, the objective function and the weak form of the problem have to be specified only once. This ensures consistency between the simulation code and its documentation. In particular, AutoFreeFem provides the linearization of the state equation, the adjoint problem, the shape derivative, as well as a basic implementation of the level-set based mesh evolution method for shape optimization. For the computation of shape derivatives we utilize the mathematical Lagrangian approach for differentiating PDE-constrained shape functions. Differentiation is done symbolically using Sympy. In numerical experiments we verify the accuracy of the computed derivatives. Finally, we showcase the capabilities of AutoFreeFem by considering shape optimization of a non-linear diffusion problem, linear and non-linear elasticity problems, a thermo-elasticity problem and a fluid-structure interaction problem.

Optimization and Control
Vem++, a c++ library to handle and play with the Virtual Element Method

F. Dassi

2023-10-03

Abstract:The Virtual Element Method (VEM) is an extension of the Finite Element Method (FEM) used for handling polytopal meshes. This paper provides a brief introduction to the VEM for a two-dimensional Laplacian problem. Additionally, it highlights the differences between a VEM implementation and a FEM code, emphasising the main challenges associated with the VEM. Furthermore, this paper presents a possible approach to address these challenges: vem++, a c++ library specifically developed for working with VEM discretisation. The library is designed to handle various partial differential equations in two or three dimensions, arising from both academic and engineering problems. Its flexible design allows for the seamless integration of new features, such as novel polytopes' quadrature rules, solvers, and virtual element spaces.

Numerical Analysis

varFEM: variational formulation based programming for finite element methods in Matlab

Object Oriented Finite Element Method and Its Application in C++

A General Finite Element Method: Extension of Variational Analysis for Nonlinear Heat Conduction with Temperature-Dependent Properties and Boundary Conditions, and Its Implementation As Local Refinement.

Vector Form of Intrinsic Finite Element Method for Incompressible Fluids

A general finite element method

Unified Implementation of Finite Element Methods involving Jumps and Averages in Matlab

Programming of Finite Element Methods in MATLAB

Programming of linear virtual element methods in three dimensions

Implementation of Object-oriented Program for Local P-version Adaptive FEM

Some researches on finite element method

On a vectorized basic linear algebra package for prototyping codes in MATLAB

Fast MATLAB assembly of FEM matrices in 2D and 3D: Edge elements

Efficient P1-FEM for Any Space Dimension in Matlab

MFEM: A modular finite element methods library

Finite Element Analysis using MATLAB

Automating the Finite Element Method

nlfem: a Flexible 2d FEM Python Code for Nonlocal Convection-Diffusion and Mechanics

Variational Principles for Rigid Finite Element Method

Vectorized implementation of primal hybrid FEM in MATLAB

AutoFreeFem: Automatic code generation with FreeFEM++ and LaTex output for shape and topology optimization of non-linear multi-physics problems

Vem++, a c++ library to handle and play with the Virtual Element Method