Natural element method

Last updated
20 points and their Voronoi cells Euclidean Voronoi diagram.svg
20 points and their Voronoi cells

The natural element method (NEM) [1] [2] [3] is a meshless method to solve partial differential equation, where the elements do not have a predefined shape as in the finite element method, but depend on the geometry. [4] [5] [6]

A Voronoi diagram partitioning the space is used to create each of these elements.

Natural neighbor interpolation functions are then used to model the unknown function within each element.

Applications

When the simulation is dynamic, this method prevents the elements to be ill-formed, having the possibility to easily redefine them at each time step depending on the geometry.

Related Research Articles

A discrete element method (DEM), also called a distinct element method, is any of a family of numerical methods for computing the motion and effect of a large number of small particles. Though DEM is very closely related to molecular dynamics, the method is generally distinguished by its inclusion of rotational degrees-of-freedom as well as stateful contact and often complicated geometries. With advances in computing power and numerical algorithms for nearest neighbor sorting, it has become possible to numerically simulate millions of particles on a single processor. Today DEM is becoming widely accepted as an effective method of addressing engineering problems in granular and discontinuous materials, especially in granular flows, powder mechanics, and rock mechanics. DEM has been extended into the Extended Discrete Element Method taking heat transfer, chemical reaction and coupling to CFD and FEM into account.

Topology optimization (TO) is a mathematical method that optimizes material layout within a given design space, for a given set of loads, boundary conditions and constraints with the goal of maximizing the performance of the system. Topology optimization is different from shape optimization and sizing optimization in the sense that the design can attain any shape within the design space, instead of dealing with predefined configurations.

Smoothed-particle hydrodynamics Method of hydrodynamics simulation

Smoothed-particle hydrodynamics (SPH) is a computational method used for simulating the mechanics of continuum media, such as solid mechanics and fluid flows. It was developed by Gingold and Monaghan and Lucy in 1977, initially for astrophysical problems. It has been used in many fields of research, including astrophysics, ballistics, volcanology, and oceanography. It is a meshfree Lagrangian method, and the resolution of the method can easily be adjusted with respect to variables such as density.

Mesh generation Subdivision of space into cells

Mesh generation is the practice of creating a mesh, a subdivision of a continuous geometric space into discrete geometric and topological cells. Often these cells form a simplicial complex. Usually the cells partition the geometric input domain. Mesh cells are used as discrete local approximations of the larger domain. Meshes are created by computer algorithms, often with human guidance through a GUI, depending on the complexity of the domain and the type of mesh desired. A typical goal is to create a mesh that accurately captures the input domain geometry, with high-quality (well-shaped) cells, and without so many cells as to make subsequent calculations intractable. The mesh should also be fine in areas that are important for the subsequent calculations.


Marching squares is a computer graphics algorithm that generates contours for a two-dimensional scalar field. A similar method can be used to contour 2D triangle meshes.

Fluid–structure interaction

Fluid–structure interaction (FSI) is the interaction of some movable or deformable structure with an internal or surrounding fluid flow. Fluid–structure interactions can be stable or oscillatory. In oscillatory interactions, the strain induced in the solid structure causes it to move such that the source of strain is reduced, and the structure returns to its former state only for the process to repeat.

Meshfree methods Methods in numerical analysis not requiring knowledge of neighboring points

In the field of numerical analysis, meshfree methods are those that do not require connection between nodes of the simulation domain, i.e. a mesh, but are rather based on interaction of each node with all its neighbors. As a consequence, original extensive properties such as mass or kinetic energy are no longer assigned to mesh elements but rather to the single nodes. Meshfree methods enable the simulation of some otherwise difficult types of problems, at the cost of extra computing time and programming effort. The absence of a mesh allows Lagrangian simulations, in which the nodes can move according to the velocity field.

A frontal solver, conceived by Bruce Irons, is an approach to solving sparse linear systems which is used extensively in finite element analysis. It is a variant of Gauss elimination that automatically avoids a large number of operations involving zero terms.

Johann Hadji Argyris FRS was a Greek pioneer of computer applications in science and engineering, among the creators of the finite element method (FEM), and lately Professor at the University of Stuttgart and Director of the Institute for Statics and Dynamics of Aerospace Structures.

In mathematics, in particular numerical analysis, the FETI method is an iterative substructuring method for solving systems of linear equations from the finite element method for the solution of elliptic partial differential equations, in particular in computational mechanics In each iteration, FETI requires the solution of a Neumann problem in each substructure and the solution of a coarse problem. The simplest version of FETI with no preconditioner in the substructure is scalable with the number of substructures but the condition number grows polynomially with the number of elements per substructure. FETI with a preconditioner consisting of the solution of a Dirichlet problem in each substructure is scalable with the number of substructures and its condition number grows only polylogarithmically with the number of elements per substructure. The coarse space in FETI consists of the nullspace on each substructure.

Bruce Irons (engineer)

Bruce Moncur Irons was an engineer and mathematician, known for his fundamental contribution to the finite element method, including the patch test, the frontal solver and, along with Ian C. Taig, the isoparametric element concept.

The material point method (MPM) is a numerical technique used to simulate the behavior of solids, liquids, gases, and any other continuum material. Especially, it is a robust spatial discretization method for simulating multi-phase (solid-fluid-gas) interactions. In the MPM, a continuum body is described by a number of small Lagrangian elements referred to as 'material points'. These material points are surrounded by a background mesh/grid that is used to calculate terms such as the deformation gradient. Unlike other mesh-based methods like the finite element method, finite volume method or finite difference method, the MPM is not a mesh based method and is instead categorized as a meshless/meshfree or continuum-based particle method, examples of which are smoothed particle hydrodynamics and peridynamics. Despite the presence of a background mesh, the MPM does not encounter the drawbacks of mesh-based methods which makes it a promising and powerful tool in computational mechanics.

In fracture mechanics, the energy release rate, , is the rate at which energy is transformed as a material undergoes fracture. Mathematically, the energy release rate is expressed as the decrease in total potential energy per increase in fracture surface area, and is thus expressed in terms of energy per unit area. Various energy balances can be constructed relating the energy released during fracture to the energy of the resulting new surface, as well as other dissipative processes such as plasticity and heat generation. The energy release rate is central to the field of fracture mechanics when solving problems and estimating material properties related to fracture and fatigue.

Gmsh

Gmsh is a finite-element mesh generator developed by Christophe Geuzaine and Jean-François Remacle. Released under the GNU General Public License, Gmsh is free software.

Nektar++ is a spectral/hp element framework designed to support the construction of efficient high-performance scalable solvers for a wide range of partial differential equations (PDE). The code is released as open-source under the MIT license. Although primarily driven by application-based research, it has been designed as a platform to support the development of novel numerical techniques in the area of high-order finite element methods.

The finite point method (FPM) is a meshfree method for solving partial differential equations (PDEs) on scattered distributions of points. The FPM was proposed in the mid-nineties in, and with the purpose to facilitate the solution of problems involving complex geometries, free surfaces, moving boundaries and adaptive refinement. Since then, the FPM has evolved considerably, showing satisfactory accuracy and capabilities to deal with different fluid and solid mechanics problems.

Sara Zahedi is an Iranian-Swedish mathematician who works in computational fluid dynamics and holds an associate professorship in numerical analysis at the Royal Institute of Technology (KTH) in Sweden. She is one of ten winners and the only female winner of the European Mathematical Society Prize for 2016 "for her outstanding research regarding the development and analysis of numerical algorithms for partial differential equations with a focus on applications to problems with dynamically changing geometry". The topic of Zahedi's EMS Prize lecture was her recent research on the CutFEM method of solving fluid dynamics problems with changing boundary geometry, such as arise when simulating the dynamics of systems of two immiscible liquids. This method combines level set methods to represent the domain boundaries as cuts through an underlying uniform grid, together with numerical simulation techniques that can adapt to the complex geometries of grid cells cut by these boundaries.

In numerical analysis, multi-time-step integration, also referred to as multiple-step or asynchronous time integration, is a numerical time-integration method that uses different time-steps or time-integrators for different parts of the problem. There are different approaches to multi-time-step integration. They are based on domain decomposition and can be classified into strong (monolithic) or weak (staggered) schemes. Using different time-steps or time-integrators in the context of a weak algorithm is rather straightforward, because the numerical solvers operate independently. However, this is not the case in a strong algorithm. In the past few years a number of research articles have addressed the development of strong multi-time-step algorithms. In either case, strong or weak, the numerical accuracy and stability needs to be carefully studied. Other approaches to multi-time-step integration in the context of operator splitting methods have also been developed; i.e., multi-rate GARK method and multi-step methods for molecular dynamics simulations.

Barna Szabó

Barna A. Szabó is a Hungarian-American engineer and educator, noted for his contributions on the finite element method, particularly the conception and implementation of the p- and hp-versions of the Finite Element Method. He is a founding member and fellow of the United States Association for Computational Mechanics, an external member of the Hungarian Academy of Sciences and fellow of the St. Louis Academy of Sciences.

Massimo Guiggiani

Massimo Guiggiani is an Italian mechanical engineer and academic with interests in vehicle dynamics, gear drives and the boundary element method (BEM). He is professor of applied mechanics at the Università di Pisa. He created an algorithm, Guiggiani's method, for the evaluation of strongly singular and hypersingular integrals. He is the author of The Science of Vehicle Dynamics.

References

  1. Sukumar, N.; Moran, B.; Belytschko, T. (21 June 1998). "The natural element method in solid mechanics". International Journal for Numerical Methods in Engineering. 43 (5): 839–887. Bibcode:1998IJNME..43..839S. doi:10.1002/(SICI)1097-0207(19981115)43:5<839::AID-NME423>3.0.CO;2-R.
  2. J. Yvonnet; D. Ryckelynck; P. Lorong; F. Chinesta (2004). "A new extension of the natural element method for non‐convex and discontinuous problems: the constrained natural element method (C‐NEM)" (PDF). International Journal for Numerical Methods in Engineering. 60 (8): 1451–1474. Bibcode:2004IJNME..60.1451Y. doi:10.1002/nme.1016. S2CID   122887431.
  3. "Large deformation analysis of elastic bodies by nonlinear Petrov–Galerkin natural element method". Advances in Mechanical Engineering. April 2019.
  4. Lu, Ping; Shu, Yang; Lu, Dahai; Jiang, Kaiyong; Liu, Bin; Huang, Changbiao (2017-01-01). "Research on Natural Element Method and the application to simulate metal forming processes". Procedia Engineering. International Conference on the Technology of Plasticity, ICTP 2017, 17-22 September 2017, Cambridge, United Kingdom. ScienceDirect. 207: 1087–1092. doi: 10.1016/j.proeng.2017.10.1135 .
  5. "What is the difference between nem (natural element method) and cnem (constrained natural element method)?". ResearchGate. Retrieved 2019-07-15.
  6. Botelho, D. P.; Marechal, Y.; Ramdane, B. (November 2016). "Vector interpolation on natural element method: Mesh sensitivity analysis". 2016 IEEE Conference on Electromagnetic Field Computation (CEFC). Institute of Electrical and Electronics Engineers. p. 1. doi:10.1109/CEFC.2016.7816353. ISBN   978-1-5090-1032-5. S2CID   27851390.