Aplicaciones de las

Ecuaciones
Diferenciales
Parciales
Ecuación diferencial

Una ecuación diferencial es una ecuación en la

que intervienen derivadas de una o más
funciones.

Las ecuaciones diferenciales se dividen en:

1. Ecuaciones diferenciales ordinarias:
aquellas que contienen derivadas respecto a
una sola variable independiente.
2. Ecuaciones en derivadas parciales: aquellas
que contienen derivadas respecto a dos o
más variables.
Ejemplo

Ecuación diferencial ordinaria

dy
= ay
dx
Solución general: y (x) = c ea x , c ∈ Rn .

Ecuación diferencial parcial

∂u ∂u
+ =0
∂x ∂y
Una de sus soluciones: u(x, y ) = x − y .
El inicio de las Ecuaciones Diferenciales Parciales

El estudio de las ecuaciones diferenciales

parciales inicio en el siglo XVIII con los trabajos
de los suizos d’Alembert y Euler y los franceses
Lagrange y Laplace.

Las ecuaciones diferenciales parciales

aparecieron en el contexto de la modelización
matemática de fenómenos de la física del
medio continuo.
Modelos

I Cuerdas vibrantes
I El campo gravitacional Newtoniano
I Conducción del calor
I Dinámica de fluidos
I Elasticidad
I Electrostática
I Electricidad y magnetismo
La ecuación de onda unidimensional
d’Alembert, 1752: modelo de una cuerda
vibrante
∂ 2u 2
2 ∂ u
= c
∂t 2 ∂x 2

También se escibe

utt = c 2 uxx

c =velocidad de la onda.

Jean le Rond d’Alembert

(1717–1783)
Derivación de la ecuación de onda
4x = longitud de un tramo 4m = su masa µ = densidad.
Componente horizontal de la tensión T ≈ T1 cos(α) = T2 cos(β)
Segunda ley de Newton: F = 4m a

∂2u
T2 sen(β) − T1 sen(α) = µ 4x
∂2t
µ 4x ∂ 2 u
tan(β) − tan(α) =
T ∂2t
µ ∂2u
 
1 ∂u ∂u
(x + 4x, t) − (x, t) =
4x ∂x ∂x T ∂2t
Cuando 4x → 0, se obtiene

∂2u ∂2u
r
µ
2
= c2 c=
∂t ∂x 2 T
El Problema de Cauchy

Proporcionar el desplazamiento y velocidad

inicial

u(0, x) = f (x), ut (0, x) = g(x)

Solución (fórmula de d’Alembert):

Z x+ct
1 1
u(t, x) = [f (x − ct) + f (x + ct)]+ g(y ) dy .
2 2c x−ct
La ecuación de onda multidimensional
Euler, 1759, y D. Bernoulli, 1762: estudio de
ondas acústicas
∂ 2u
 2 2 2

∂ u ∂ u ∂ u
2
= c2 + +
∂t ∂x 2 ∂y 2 ∂z 2

Leonhard Paul Euler Daniel Bernoulli (1700–1782)

(1707–1783)
Operador Laplaciano
n
2
X ∂2
∆u = ∇ u = 2
, n = 2, 3
i
∂x i

La ecuación de onda usualmente se escribe

utt = c 2 ∆u ó utt = c 2 ∇2 u,
Aplicaciones

La ecuación de onda también describe

I Ondas de sonido en un tubo ó en una barra
I Oscilaciones torsionales de una barra
I Ondas largas en un canal recto
I Transmisión de electricidad en un cable
aislado de baja resistencia
I Flujos de aire supersónicos
Ecuación de Laplace

Laplace, aprox. 1780: estudio de campos

potenciales gravitacionales
∂ 2ϕ ∂ 2ϕ ∂ 2ϕ
+ + =0
∂x 2 ∂y 2 ∂z 2
Otras formas de escribirla son:

ϕxx + ϕyy + ϕzz = 0

div grad ϕ = 0
∇2 ϕ = 0
∆ϕ = 0 Pierre-Simon, marqués de

Laplace (1749–1827)
Funciones armónicas

Las soluciones ϕ(x, y , z) de esta ecuación se denominan funciones

armónicas y son importantes en muchos campos de la ciencia,
principalmente:

I Electromagnetismo
I Astronomía
I Dinámica de fluidos

ya que describen el comportamiento de los potenciales eléctricos,

gravitacionales y de los fluidos.

La teoría general de las soluciones de la ecuación de Laplace se

conoce como Teoría del Potencial.
Problema con valores en al frontera
Problema típico: encontrar una solución que satisface valores
arbitrarios en la frontera ∂Ω de un dominio acotado Ω.

Ejemplo: encontrar ϕ(r , θ) armónica sobre el

círculo de radio uno, tal que ϕ(1, θ) = φ(θ).
Solución (Poisson):
2π
1 − r2
Z
1
ϕ(r , θ) = u(θ0 )dθ0
2π 0 1+ r2 − 2r cos(θ − θ0 )

Siméon Denis Poisson

(1781-1840)
Poisson también es conocido por la ecuación que lleva su nombre

∆ϕ = f

la cual tiene aplicaciones en electrostática, ingeniería mecánica y

física teórica
Ecuación de calor

Fourier, 1810–1822, “Théorie analytique de la chaleur”)

∂ 2u ∂ 2u ∂ 2u
 
∂u
=k + +
∂t ∂x 2 ∂y 2 ∂z 2
Otras formas de expresarla son:
I utt = k [uxx + uyy + uzz ]
I utt = k ∇2 u
Jean Baptiste Joseph Fourier,
I utt = k ∆u
1768–1830
Derivación de la ecuación de calor
Ley de Fourier: q = −κ ∇u (flujo local de calor)

La enegía térmica total Q contenida en una región V es

ZZZ
Q= cs ρ u(x, t) dV
V

El flujo de calor que atraviesa la superficie S de la región V es

ZZ ZZ
− q · n dS = κ ∇u · n dS
S S

Conservación de la energía (suponiendo no hay fuentes)

ZZZ ZZ  ZZZ 
dQ du
= cs ρ dV = κ ∇u · n dS = ∇ · (κ ∇u) dv
dt V dt S V

∴ utt = k ∆u con k = κ/cs ρ

Problema de Cauchy
Especificar una distribución inicial u(x, 0) = f (x).
Ejemplo: El problema

ut = k uxx
u(x, 0) = cos(2x)
ux (0, t) = u(1, t) = 0

tiene solución
∞
sen(2) X 2 2
u(x, t) = + An cos(nπx) e−kn π t
2
1
Z 1
An = 2 cos(2x) cos(nπx)
0

Augustin Louis Cauchy,

1789–1857
Aplicaciones de la ecuación de calor

La ecuación de calor es de fundamental importancia en diversos

campos y temas cientificos:
I Procesos de difusión: en física, química y procesos industriales.
Difusión de partículas en pilas atómicas, difusión de vorticidad
en flujos viscosos. Movimiento lento en hidrodinámica
I Estadística: esta conectada con el movimiento Browniano
I Probabilidad: describe caminatas aleatorias. Evolución de
distribuciones en procesos aleatorios. Ecuación de Schrödinger
para una partícula libre.
I Matemática financiera: modelación de opciones
(Black–Scholes)
I Transmisión telegráfica en cables de baja conductancia o
capacitancia
Clasificación de las EDP

Entonces, los tres principales tipos de EDP de la física matemática

fueron introducidas entre la mitad del siglo XVIII y principios del siglo
XIX:

EDP lineal de 2o. orden en dos variables:

A uxx + B uxy + C uyy + . . . = 0

Secciones cónicas: A x 2 + B x y + C y 2 + . . . = 0

1. B 2 − 4AC < 0: Elíptica (Ecuación de Laplace)

2
2. B − 4AC = 0: Parabólica (Ecuación de calor)
3. B 2 − 4AC > 0: Hiperbólica (Ecuación de onda)
Hay muchas ecuaciones de 2o. orden que no se pueden clasificar.
Por ejemplo:

Ecuación de Euler–Tricomi uxx = x uyy

útil en el estudio de flujos transónicos (aeronáutica), es:

I Hiperbólica si x > 0
I Elíptica si x < 0
I Parabólica (de tipo degenerado) sobre la linea x = 0

F/A–18 volando a velocidad transónica

Algunos métodos de solución

I Separación de variables
I Principio de superposición
I Cambio de variables
I Método de las caracterísiticas. Tipicamente
se aplica a ecs. de 1er. orden
I Series de Fourier
I Métodos de Tranformadas (Laplace,
Fourier,...)
I Funciones de Green
I Métodos de aproximación
La realidad es mucho más compleja

En aplicaciones reales las ecuaciones son más complejas y

usualmente no son lineales.

Ejemplo: en la ecuación del calor la constante de difusión térmica k

realmente depende de la solución

ut = k (u) ∆u

Esta es una EDP no–lineal

I Las EDP no lineales no tienen solución analítica, excepto en

pocos casos muy simples.

I Más aún, las EDP lineales no tienen solución analítica en

regiones no simples
Otras EDP Clásicas, 1750–1900
Ecuaciones de Euler, 1755
Flujos de fluidos inviscidos
∂ρ
+ ∇ · (ρu) = 0
∂t
∂ρu
+ ∇ · (u ⊗ (ρu)) + ∇p = 0
∂t
∂E
+ ∇ · (u(E + p)) = 0,
∂t

I ρ es la densidad del fluido

I u = (u, v , w) es el vector de velocidades
I E = ρ e + 12 (u 2 + v 2 + w w 2) es la energía por u. de volumen
I p es la presión
La segunda ecuación es la divergencia de un producto diádico y
puede expresarse en la forma (suponendo ρ constante)
 
∂
ρ + u · ∇ u + ∇p = 0
∂t
 las ondas descritas
Al igual que las olas por la ecuación de
Las ecuaciones de oceánicas Euler se “rompen” y
Euler son hiperbólicas
se forman las ondas
no lineales y sus
de choque
soluciones generales
son ondas
Ecuaciones de Navier–Stokes
Navier, 1822–1827 −→ Poisson, 1831 −→ Stokes, 1845: movimiento
de fluidos viscosos como agua y aire
 
∂u
ρ + u · ∇u = −∇p + ∇ · T + f
∂t
∂ρ
+ ∇ · (ρu) = 0
∂t
Uno de los conjuntos de EDP más útiles: describe la física de un
gran numéro de fenómenos de interés científico, académico y
económico.

Claude Louis Navier, 1785–1836 George Gabriel Stokes, 1819–1903

Fluidos incompresibles
Inercia
z }| {
 ∂u 
ρ | ·{z
+ u ∇u} = −∇p + µ∇2 u + |{z}
f
∂t | {z } | {z }
|{z} Aceleración Gradiente Viscosidad Otras
Aceleración convectiva de presión fuerzas
no estacionaria
∇·u=0
Coordenadas cartesianas: u = (u, v , w)
Ecuaciones de momentum
 2
∂2u ∂2u
  
∂u ∂u ∂u ∂u ∂p ∂ u
ρ +u +v +w =− +µ + + + ρgx
∂t ∂x ∂y ∂z ∂x ∂x 2 ∂y 2 ∂z 2
 2
∂2v ∂2v
  
∂v ∂v ∂v ∂v ∂p ∂ v
ρ +u +v +w =− +µ + + + ρgy
∂t ∂x ∂y ∂z ∂y ∂x 2 ∂y 2 ∂z 2
 2
∂2w ∂2w
  
∂w ∂w ∂w ∂w ∂p ∂ w
ρ +u +v +w =− +µ + + + ρgz
∂t ∂x ∂y ∂z ∂z ∂x 2 ∂y 2 ∂z 2
Ecuación de continuidad
∂u ∂v ∂w
+ + =0
∂x ∂y ∂z
Coordenadas cilíndricas
Ecuaciones de momentum
uθ2
 
∂ur ∂ur uθ ∂ur ∂ur ∂p
ρ + ur + + uz − =−
∂t ∂r r ∂θ ∂z r ∂r
1 ∂ 2 ur ∂ 2 ur
   
∂ 1 ∂rur ur 2 ∂uθ
+µ + 2 + − − + ρgr
∂r r ∂r r ∂θ2 ∂z 2 r2 r 2 ∂θ
 
∂uθ ∂uθ uθ ∂uθ ∂uθ ur uθ 1 ∂p
ρ + ur + + uz + =−
∂t ∂r r ∂θ ∂z r r ∂θ
2 2
   
∂ 1 ∂ruθ 1 ∂ uθ ∂ uθ 2 ∂ur uθ
+µ + 2 + + 2 − 2 + ρgθ
∂r r ∂r r ∂θ2 ∂z 2 r ∂θ r
 
∂uz ∂uz uθ ∂uz ∂uz ∂p
ρ + ur + + uz =−
∂t ∂r r ∂θ ∂z ∂z
1 ∂ 2 uz ∂ 2 uz
   
1 ∂ ∂uz
+µ r + 2 + + ρgz
r ∂r ∂r r ∂θ2 ∂z 2
Ecuacion de continuidad
1 ∂ 1 ∂uθ ∂uz
(rur ) + + = 0.
r ∂r r ∂θ ∂z
Aplicaciones

I Modelan el movimiento de líquidos y gases

I Modelan corrientes oceánicas
I Modelan el movimiento de estrellas dentro de las galaxias
I Modelan fenómenos complejos en biología
I Útiles studio el flujo de sangre
I Útiles en el diseño de autos y aeroplanos
I Útiles en el estudio fenómenos meteorológicos
I Útiles en el estudio de transporte de contaminantes
I Y muchas otras aplicaciones
Elasticidad lineal
Navier, 1821 y Cauchy, 1822:
Es el estudio matemático de como se
deforman los objetos sólidos y de sus
esfuerzos internos debido a cargas externas
dadas (Mecánica del medio continuo)

Ecuaciones de Navier–Cauchy
(λ + µ)∇(∇ · u) + µ∇2 u + F = 0
En notación extendida (ingeniería)
 2
∂ 2 ux ∂ 2 ux
  
∂ ∂ux ∂uy ∂uz ∂ ux
(λ + µ) + + +µ + + + Fx = 0
∂x ∂x ∂y ∂z ∂x 2 ∂y 2 ∂z 2
 2
∂ 2 uy ∂ 2 uy
  
∂ ∂ux ∂uy ∂uz ∂ uy
(λ + µ) + + +µ + + + Fy = 0
∂y ∂x ∂y ∂z ∂x 2 ∂y 2 ∂z 2
 2
∂ 2 uz ∂ 2 uz
  
∂ ∂ux ∂uy ∂uz ∂ uz
(λ + µ) + + +µ + + + Fz = 0
∂z ∂x ∂y ∂z ∂x 2 ∂y 2 ∂z 2
donde µ el el módulo de rígidez, λ es un parámetro de Lamé, u es
desplazamiento, F son las fuerzas de cuerpo.
Las ecuaciones anteriores son un caso particular de

σij,j + Fi = ρ∂tt ui Ecuaciones de movimiento (2a. ley de Newton)

σij = Cijkl εkl Ecuaciones constitutivas (Ley de Hooke)
εij = (uj,i + ui,j )/2 Ec. deformación–desplazamiento

donde

I σij = σji Tensor de esfuerzos

I ρ Densidad
I Cijkl Tensor de elasticidad
I λ, µ Constantes materiales
I εij = εji Tensor de deformación
Ecuaciones de Maxwell: Teoría electromagnética, 1864

ρ
Ley de Gauss ∇·E=
ε0
(campo eléctrico)
Ley de Gauss ∇·B=0
(campo magnético)
∂B
Ley de Faraday ∇×E=−
∂t
∂E
Ley de Ampère ∇ × B = µ0 J + µ0 ε0
∂t James Clerk Maxwell,

1831–1879

E: campo eléctrico B: campo magnético J: densidad de corriente

ρ: densidad de carga 0 : pemitividad µ0 : permeabilidad

Formas integrales
I
QS
E · dA =
ε0
IS
B · dA = 0 Teorema de la divergencia
IS
∂ΦB,S ZZZ ZZ
E · dl = − (∇ · F) dV = F·n dS
∂S ∂t V
I
∂ΦE,S
B · dl = µ0 IS + µ0 ε0
∂S ∂t

QS : Carga neta encerrada por S

R
IS = S J · dA: Corriente eléctrica que
atraviesa S
R
ΦE,S = S E · dA: Flujo eléctrico sobre S Teorema de Stokes
Z I
R ∇ × F · dΣ = F · dr
ΦB,S = S B · dA: Flujo magnético sobre S Σ ∂Σ
 Controversia sobre el papel de Maxwell
Contibrución: corrección en el término de corriente de
desplazamiento (On Physical Lines of Force 1861) Le permitió
derivar la ecuación de onda electromagnética en su artículo
Dynamical Theory of the Electromagnetic Field

1 ∂2
 
∇2 − E = 0
c 2 ∂t 2
1 ∂2
 
2
∇ − 2 2 B = 0
c ∂t

La luz es una onda electromágnetica

(confirmado experimentalmente por
Heinrich Hertz en 1887).
Se dice que las ecuaciones de Maxwell se llamaron ecuaciones de
Hertz-Heaviside pero que Einstein, por alguna razón, les llamo
ecuaciones de Maxwell–Hertz.
Johann Carl Friedrich Gauss André–Marie Ampère (1775–1836)

Michael Faraday (1791–1867

(1777–1855)
Ecuación de Korteweg–De Vries, 1896

Modelo de ondas solitarias en

aguas someras.

φt + φxxx + 6 φx = 0

I EDP no–lineal
I Soluciones analíticas y
precisas
I Apareció primero en 1877 Gustav de Vries,
Diederik Johanes Korteweg,
en un artículo de (1866–1934)
(1848–1941)
Boussinesq en Francia
Solución
Método: inverse scattering transform.

Solución con un soliton

1 1
φ(x, t) = c h√ i
2 cosh2 c
2 (x − ct − a)

c: rapidez de fase a: constante.

Solución con dos solitones:

∂2
φ(x, t) = −2 ln(1 + B1 eθ1 + B2 eθ2 + A B1 B2 eθ1 +θ2 )
∂x 2
 2
a1 − a2
θ1 = a1 x − a13 t θ2 = a2 x − a23 t A=
a1 + a2
Aplicaciones

La teoría matemática de la ecuación de KdV es

muy rica e interesante, y desde hace algunas
décadas es un tema de investigación muy
activo.
I Física y Mecánica no lineal
I Ondas en aguas someras
I Ondas internas largas en un oceano con densidad estratificada
I Ondas ion–acústicas en plasmas
I Ondas acústicas sobre un látice de cristal
I Neurociencia, conducción de señales en las neuronas
Otras ecuaciones

I La ecuación de Monge–Ampère (Monge 1755)

I La ecuación de superficie mínima (Lagrange 1760)
I La ecuaciones de Laplace y Poisson, como se aplica en
problemas de electricidad y magnetismo, comenzando por
Poisson en 1813, el libro de Green en 1828 y Gauss en 1839.
I El problema de Plateau, como modelo de burbujas de jabón, en
los 1840’s
I La ecuación de Helmholtz y el problema de valores propios para
el operador de Laplace en conexión con la acústica (1860)