Are you an EPFL student looking for a semester project?
Work with us on data science and visualisation projects, and deploy your project as an app on top of Graph Search.
Concept
Stiffness matrix
Summary
In the finite element method for the numerical solution of elliptic partial differential equations, the stiffness matrix is a matrix that represents the system of linear equations that must be solved in order to ascertain an approximate solution to the differential equation.
For simplicity, we will first consider the Poisson problem
on some domain Ω, subject to the boundary condition u = 0 on the boundary of Ω. To discretize this equation by the finite element method, one chooses a set of basis functions {φ_1, ..., φ_n} defined on Ω which also vanish on the boundary. One then approximates
The coefficients u_1, u_2, ..., u_n are determined so that the error in the approximation is orthogonal to each basis function φ_i:
The stiffness matrix is the n-element square matrix A defined by
By defining the vector F with components the coefficients u_i are determined by the linear system Au = F. The stiffness matrix is symmetric, i.e. A_ij = A_ji, so all its eigenvalues are real. Moreover, it is a strictly positive-definite matrix, so that the system Au = F always has a unique solution. (For other problems, these nice properties will be lost.)
Note that the stiffness matrix will be different depending on the computational grid used for the domain and what type of finite element is used. For example, the stiffness matrix when piecewise quadratic finite elements are used will have more degrees of freedom than piecewise linear elements.
Determining the stiffness matrix for other PDEs follows essentially the same procedure, but it can be complicated by the choice of boundary conditions. As a more complex example, consider the elliptic equation
where is a positive-definite matrix defined for each point x in the domain. We impose the Robin boundary condition
where ν_k is the component of the unit outward normal vector ν in the k-th direction. The system to be solved is
as can be shown using an analogue of Green's identity. The coefficients u_i are still found by solving a system of linear equations, but the matrix representing the system is markedly different from that for the ordinary Poisson problem.
Official source
About this result
This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.