Saltar al contenido
Introducción

Hablemos un poco de lo que haremos:

Este curso es una introducción práctica a los métodos numéricos aplicados a ecuaciones diferenciales parciales (EDPs), utilizando Python y Jupyter Notebooks. Se enfoca en la ecuación de onda, que es fundamental en muchos fenómenos de las ciencias de la Tierra. Los métodos abordados incluyen diferencias finitas, método pseudoespectral, elementos finitos y elementos espectrales.

Está diseñado para estudiantes o investigadores que buscan comprender la implementación práctica de simulaciones y los fundamentos de la computación paralela. Aunque el curso incluye matemáticas, no es teórico; se centra en cómo discretizar EDPs, implementarlas en código y garantizar resultados precisos.

En esta primera parte abordaremos conceptos como discretización, mallas computacionales, computación paralela y una introducción a Python y Jupyter Notebooks. Posteriormente, se aplicará el método de diferencias finitas a la ecuación de onda.

Resumen

Ecuaciones diferenciales parciales (EDP)s

Una ecuación diferencial parcial (EDP) es una ecuación matemática que describe cómo cambia una cantidad en función de múltiples variables independientes. A diferencia de las ecuaciones diferenciales ordinarias (EDOs), que dependen de una sola variable independiente, las EDPs involucran derivadas parciales con respecto a dos o más variables independientes.

Matemáticamente, una EDP tiene la forma general:
\( F\left(x_1, x_2, \dots, x_n, u, \frac{\partial u}{\partial x_1}, \frac{\partial u}{\partial x_2}, \dots, \frac{\partial^2 u}{\partial x_1^2}, \dots \right) = 0, \)
donde \( u = u(x_1, x_2, \dots, x_n) \) es la incógnita que queremos determinar, y \( F \) es una función que relaciona \( u \) y sus derivadas parciales de diferentes órdenes.

Algunos tipos comunes de EDPs son:

  1. Ecuación de difusión/calor
    Modela la distribución de una cantidad como el calor en un medio:
    \( \frac{\partial u}{\partial t} = \alpha \nabla^2 u, \)
    donde \( \nabla^2 u \) es el operador laplaciano que involucra derivadas espaciales.

  2. Ecuación de onda
    Describe la propagación de ondas como las sonoras o sísmicas:
    \( \frac{\partial^2 u}{\partial t^2} = c^2 \nabla^2 u, \)
    donde \( c \) es la velocidad de propagación.

  3. Ecuación de Laplace/Poisson
    Modela fenómenos estacionarios como potenciales eléctricos:
    \( \nabla^2 u = 0 \quad \text{(Laplace)}, \quad \nabla^2 u = f(x, y) \quad \text{(Poisson)}. \)

Las EDPs son esenciales en la modelación de fenómenos físicos, como el clima, la propagación de ondas o el flujo de fluidos. En la mayoría de los casos, estas ecuaciones no tienen soluciones exactas analíticas, por lo que se requieren métodos numéricos para aproximarlas. Esto las hace fundamentales en áreas como geofísica, meteorología, y simulación computacional.

Aplicación teórica

Algunas aplicaciones:

Las ecuaciones diferenciales parciales (EDP) son fundamentales en la modelización de una gran variedad de fenómenos en ciencias e ingeniería. Aquí tienes varios ejemplos de su uso en diferentes campos:

Física y Termodinámica:

Ecuación del calor: La ecuación de difusión o de calor describe cómo la temperatura u(x,t) cambia en el tiempo t y el espacio x. Es una EDP de la forma:

\[ \frac{\partial u}{\partial t} = \alpha \frac{\partial^2 u}{\partial x^2} \]

donde \(\alpha\) es la difusividad térmica. Se utiliza para modelar la conducción del calor en una barra, o la distribución de calor en un objeto.

Ecuación de Laplace: La ecuación de Laplace describe el comportamiento de campos escalarizados, como el potencial eléctrico o el campo gravitacional, en un medio sin fuentes. Es de la forma:

\[ \nabla^2 \phi = 0 \]

donde \(\nabla^2\) es el operador laplaciano y \(\phi\) es el potencial.

Mecánica de Fluidos:

Ecuaciones de Navier-Stokes: Estas ecuaciones describen el movimiento de fluidos viscosos. Son fundamentales en la dinámica de fluidos y tienen la forma general:

\[ \frac{\partial u}{\partial t} + u \cdot \nabla u = -\nabla p + \nu \nabla^2 u \]

donde \(u\) es el campo de velocidades del fluido, \(p\) es la presión, y \(\nu\) es la viscosidad cinemática. Se usan en la simulación de flujos de aire, agua y otros fluidos.

Ecuación de la onda: En la acústica y la electromagnética, la propagación de ondas se describe mediante la ecuación de ondas. Un ejemplo es la ecuación para ondas sonoras o electromagnéticas:

\[ \frac{\partial^2 u}{\partial t^2} = c^2 \nabla^2 u \]

donde \(u\) es la perturbación en el medio y \(c\) es la velocidad de propagación.

Electromagnetismo:

Ecuaciones de Maxwell: Las ecuaciones que gobiernan los campos eléctricos y magnéticos en el espacio y el tiempo son un conjunto de EDPs. Un ejemplo es la ecuación de Gauss para el campo eléctrico:

\[ \nabla \cdot E = \frac{\rho}{\epsilon_0} \]

donde \(E\) es el campo eléctrico y \(\rho\) es la densidad de carga eléctrica. Estas ecuaciones describen cómo se propagan los campos electromagnéticos.

Biología y Ecología:

Modelo de difusión y crecimiento celular: Las EDPs también se utilizan para modelar la difusión de sustancias químicas, como en el modelo de difusión de nutrientes o en la migración de células:

\[ \frac{\partial u}{\partial t} = D \nabla^2 u \]

donde \(u\) representa la concentración de la sustancia y \(D\) es el coeficiente de difusión.

Ecuaciones de reacción-difusión: Modelan fenómenos como la propagación de enfermedades o la distribución de especies. Tienen la forma general:

\[ \frac{\partial u}{\partial t} = D \nabla^2 u + f(u) \]

donde \(f(u)\) es una función que describe la reacción (como el crecimiento de una población).

Finanzas:

Ecuación de Black-Scholes: En la teoría financiera, la ecuación de Black-Scholes describe la evolución del precio de opciones financieras. Su forma es:

\[ \frac{\partial C}{\partial t} + rS \frac{\partial C}{\partial S} + \frac{1}{2} \sigma^2 S^2 \frac{\partial^2 C}{\partial S^2} = rC \]

donde \(C\) es el precio de la opción, \(S\) es el precio del activo subyacente, \(\sigma\) es la volatilidad, y \(r\) es la tasa de interés.

Ingeniería:

Ecuaciones estructurales: En el diseño de estructuras, las EDPs describen la deformación de materiales bajo fuerzas externas. La ecuación general de la elasticidad es:

\[ \nabla \cdot (\sigma) + f = 0 \]

donde \(\sigma\) es el tensor de tensiones y \(f\) es el vector de fuerzas externas.

Astronomía:

Ecuación de Poisson: En gravitación, la ecuación de Poisson es una EDP que describe cómo el potencial gravitacional \(\Phi\) se relaciona con la distribución de masa \(\rho\):

\[ \nabla^2 \Phi = 4 \pi G \rho \]

donde \(G\) es la constante de gravitación y \(\rho\) es la densidad de masa.

Geofísica:

Ecuación de Ondas Sísmicas: Las ondas sísmicas que se propagan a través de la Tierra son modeladas mediante ecuaciones diferenciales parciales. Una de las ecuaciones básicas es:

\[ \frac{\partial^2 u}{\partial t^2} = c^2 \nabla^2 u \]

donde \(u\) es el desplazamiento de las partículas en el medio y \(c\) es la velocidad de propagación de la onda.

Estos son solo algunos ejemplos, pero las EDPs tienen una aplicación vasta en muchas otras áreas, como la teoría de cuerdas, la dinámica de poblaciones, la óptica, la química, la ingeniería de materiales, y la meteorología.

Aplicación práctica

Hagamos un ejemplo sencillo para ilustar un poco.

Un ejemplo sencillo y teórico de la aplicación de las ecuaciones diferenciales parciales (EDP) es el estudio de la distribución de temperatura en una barra metálica uniforme. Este problema se describe mediante la ecuación del calor, que modela cómo varía la temperatura \( u(x,t) \) a lo largo de la barra y en el tiempo.

La ecuación del calor en una dimensión es:

\[ \frac{\partial u}{\partial t} = \alpha \frac{\partial^2 u}{\partial x^2}, \]

donde:

  1. \( u(x,t) \) representa la temperatura en el punto \( x \) de la barra en el tiempo \( t \).
  2. \( \alpha > 0 \) es la difusividad térmica del material, una constante que depende de sus propiedades físicas.
  3. \( \frac{\partial u}{\partial t} \) describe el cambio temporal de la temperatura.
  4. \( \frac{\partial^2 u}{\partial x^2} \) es el término que representa el flujo de calor debido a diferencias de temperatura a lo largo de la barra.

Contexto del problema

Imagina una barra metálica larga y delgada con extremos en \( x = 0 \) y \( x = L \). La temperatura inicial en cada punto de la barra está definida por una función \( u(x,0) = f(x) \). Además, los extremos pueden mantenerse a una temperatura fija o ser aislados térmicamente.

Condiciones

  1. Condición inicial: Define la distribución de temperatura al inicio, como \( u(x,0) = f(x) \).
  2. Condiciones de frontera:
    • \( u(0,t) \) y \( u(L,t) \) especifican las temperaturas en los extremos de la barra (pueden ser constantes o funciones del tiempo).
    • Alternativamente, puede asumirse que los extremos son aislados, lo que implica que \( \frac{\partial u}{\partial x} = 0 \) en \( x = 0 \) y \( x = L \).

Solución

La solución analítica de este problema, en casos simples, se obtiene mediante técnicas como la separación de variables, series de Fourier o transformadas. En esencia, la solución describe cómo la temperatura en cada punto de la barra evoluciona hacia un estado de equilibrio, donde la temperatura es uniforme.