Finite Element Method

The Finite Element method is a powerful tool for finding approximate solutions to continuous problems. The methodology is similar to other numerical techniques that approximate continuous partial differential equations with discrete algebraic equations.

Simcenter STAR-CCM+ currently uses the Finite Element method in Solid Mechanics, Electromagnetism, and Viscous Flow applications.

Discretization of the Continuous Domain

Mathematical models describing physics phenomena usually consist of a set of partial differential equations, which, together with constitutive equations, form a closed set of equations. The first step for reducing the continuous set to a discrete set of algebraic equations is the discretization of the continuous space domain into a finite number of subdomains, which are called elements. Typical element shapes are triangles, quadrilaterals, tetrahedra, and hexahedra. Tetrahedra and hexahedra can be connected using transition elements like prisms and pyramids.

Weak Form of PDEs

The partial differential equations (PDE) are converted to an integral or weak form that is satisfied over each element domain. Typically, the weak form of a PDE is constructed as follows:
  • The PDE is multiplied by an arbitrary test function and integrated over the element domain.
  • The order of the derivatives is reduced using integration by parts.
  • Appropriate Dirichlet and Neumann boundary conditions are applied.

Discretization of the Unknowns

In each element, the distribution of the dependent variables is constructed from the discrete values that the variables assume at specific locations, for example, at the element nodes or edges. Shape functions interpolate the discrete values to any other point in the element, providing the local (per-element) distribution of the dependent variables. Mapping operations between the local element domain and the global physical domain ensure that the variables are continuous across adjacent elements.

All integrals are approximated with appropriate numerical techniques, leading to the final set of discrete algebraic equations.

Solution of the Linear Algebraic Equations

The system of algebraic equations is solved using a sparse direct solver.