Técnica numérica que aproxima derivadas en ecuaciones diferenciales mediante expresiones algebraicas de diferencias finitas sobre una malla discreta, convirtiendo las ecuaciones gobernantes continuas en un sistema de ecuaciones algebraicas para resolver las variables de campo; se aplica habitualmente a problemas transitorios y estacionarios en ingeniería estructural, térmica y geotécnica sobre mallas estructuradas.

Principio

Las derivadas se sustituyen por cocientes de diferencias finitas (hacia adelante, hacia atrás o centradas) sobre los espaciamientos de la malla de modo que el operador diferencial se convierte en un 'stencil' que relaciona la incógnita de un nodo con los valores de nodos vecinos; la elección del orden de discretización y del esquema temporal determina el error de truncamiento, la estabilidad y la convergencia y, por tanto, la fidelidad de la solución numérica.

Demostración

Situación: Conducción térmica unidimensional transitoria en una capa de suelo. Reconocimiento: La ecuación de calor se discretiza espacialmente con diferencias centradas y temporalmente con un paso explícito de Euler hacia adelante. Acción: El ingeniero implementa el M.D.F. explícito con Δx y Δt fijos que satisfacen la condición de estabilidad CFL; calcula la evolución de la temperatura en nodos discretos. Consecuencia: La solución aproxima con precisión las temperaturas transitorias mientras Δt cumpla la condición de estabilidad; si Δt es demasiado grande, el esquema explícito se vuelve inestable y genera oscilaciones no físicas.

Aplicación incorrecta

Usar un esquema explícito de diferencias finitas sin verificar el criterio de estabilidad (condición de Courant–Friedrichs–Lewy) para la malla y el paso de tiempo elegidos; este error parece plausible porque la implementación explícita es sencilla, pero conduce a divergencia numérica o resultados espurios en lugar de predicciones físicas.

Consecuencia

Aplicado con discretización y tratamiento de fronteras adecuados, el M.D.F. produce soluciones eficientes y precisas en mallas regulares y es sencillo de implementar; si se aplica mal (esquemas inestables, mal tratamiento de fronteras o mallas toscas) genera resultados no convergentes o inexactos que pueden inducir a errores en el diseño o análisis.

Inversión

Para geometrías complejas, dominios no estructurados o cuando se necesitan tamaños de elemento variables, la conveniencia y precisión del M.D.F. en mallas estructuradas se diluye frente a la flexibilidad del FEM o FVM; variantes implícitas del M.D.F. mitigan las restricciones de paso temporal pero implican sistemas algebraicos más grandes y mayor coste computacional.

Límite

Claramente dentro: Resolver una ecuación de difusión 1D o un PDE en un dominio rectangular 2D simple sobre una malla cartesiana estructurada usando stencils de diferencias centradas. Caso límite: Aplicar M.D.F. a un dominio con curvatura geométrica moderada mediante transformaciones de coordenadas. Claramente fuera: Problemas que requieran mallas no estructuradas alrededor de fronteras complejas o refinamiento adaptativo h, donde FEM o FVM son preferibles.

Tensión semántica

Simplicidad y eficiencia computacional en mallas estructuradas (M.D.F.) ↔ Flexibilidad geométrica y consistencia variacional (FEM/FVM): el M.D.F. es simple y rápido para dominios regulares pero tiene dificultades con geometría compleja y condiciones de contorno mixtas, mientras que FEM/FVM gestionan estos aspectos de forma más natural a costa de mayor complejidad de implementación.

Síntesis

El M.D.F. transforma operadores diferenciales en stencils localizados en la malla — potente y eficiente para dominios regulares y EDP clásicas — pero la validez de los resultados depende críticamente del orden de discretización, las condiciones de estabilidad y el tratamiento exacto de las fronteras; cuando priman la geometría o la conservación, conviene considerar otras discretizaciones.