Una técnica numérica variacional que aproxima soluciones de problemas de contorno al reemplazar espacios de ensayo infinitodimensionales por espacios polinomiales por piezas de dimensión finita definidos en una malla, conduciendo a sistemas algebraicos dispersos que resuelven con solvers lineales o no lineales.

Principio

Discretizar la forma débil: elegir una malla y grado polinómico, ensamblar las contribuciones locales de los elementos en una matriz global de rigidez (y masa) mediante integrales variacionales y resolver el sistema algebraico resultante; la convergencia sigue de las propiedades de aproximación y estabilidad de los espacios discretos.

Demostración

Resolver −Δu=f con Dirichlet homogéneo en un dominio poligonal: triangular el dominio, elegir funciones base P1 (lineales por piezas) de soporte local, ensamblar la matriz de rigidez A_{ij}=∫Ω ∇ϕ_i·∇ϕ_j y el vector carga b_i=∫Ω fϕ_i, luego resolver Au=b para los valores nodales de u_h.

Aplicación incorrecta

Usar elementos no conformes sin corrección, omitir refinamiento de malla donde aparecen singularidades, o ignorar condiciones de estabilidad en problemas mixtos (inf-sup) produce aproximaciones inexactas o no convergentes.

Consecuencia

Conduce a aproximaciones prácticas y convergentes con estimaciones de error computables, estrategias adaptativas de refinamiento y solvers eficientes; permite manejar geometrías complejas y simulaciones a escala ingenieril cuando se implementa correctamente.

Inversión

Abandonar bases polinomiales locales y estructura de malla (por ejemplo, usar bases globales espectrales o métodos sin malla) cambia propiedades de aproximación, condicionamiento y compensaciones de implementación y puede ser mejor o peor según la regularidad y la clase de problema.

Límite

Se aplica a problemas que admiten una formulación variacional en dominios donde se pueden construir mallas; excluye enfoques puramente sin malla, algunos operadores no locales sin localización adecuada y casos donde la elección de elemento debe adaptarse (coeficientes singulares, dominios rugosos).

Tensión semántica

A menudo contrastada con métodos de diferencias finitas, volúmenes finitos y espectrales: FEM enfatiza la consistencia variacional, el soporte local y la adaptividad de la malla, mientras que las alternativas priorizan simplicidad, leyes de conservación o convergencia exponencial para soluciones suaves.

Síntesis

El método de elementos finitos convierte una EDP variacional en un sistema algebraico finito mediante aproximación polinómica local sobre una malla, equilibrando orden de aproximación, resolución de malla y estabilidad para obtener soluciones numéricas fiables y control de error.