Introducción

La forma en la que las EDP se presentan habitualmente en la modelización de fenómenos de la

Ciencia y Tecnología es la de modelos de evolución en los que se describe la dinámica a lo
largo del tiempo de determinada cantidad o variable (también a veces denominada estado) que
puede representar objetos y fenómenos de lo más diversos que van desde la posición de un
satélite en el espacio hasta la dinámica de un átomo, pasando por los índices bursátiles o el
grado en que una enfermedad afecta a la población. En otras palabras, los modelos dinámicos o
de evolución son los más naturales en la medida que reproducen nuestra propia concepción del
mundo: un espacio tridimensional que evoluciona y cambia en el tiempo.
Cuando el estado o variable de un modelo o sistema de evolución es finito-dimensional, el
modelo más natural es un sistema de EDO, cuya dimensión coincide precisamente con el del
número de parámetros necesarios para describir dicho estado. Así, por ejemplo, para posicionar
una partícula en el espacio necesitamos de tres variables dependientes del tiempo y para
describir su dinámica un sistema de tres ecuaciones diferenciales en la que la variable
independiente es el tiempo.
Pero en muchas ocasiones, como es el caso sistemáticamente en el contexto de la Mecánica de
Medios Continuos, la variable de estado es infinito-dimensional. Esto ocurre por ejemplo
cuando se pretende describir la deformación de cuerpos elásticos o la temperatura de un cuerpo
sólido en los que la deformación o temperatura de cada uno de los puntos de ese medio continuo
constituye una variable o incógnita del sistema.
Los modelos matemáticos naturales en este caso son las EDP. En la teoría clásica de EDP ´estas
se clasifican en tres grandes grupos: elípticas, parabólicas e hiperbólicas. El modelo elíptico por
excelencia involucra el operador de Laplace

IMPORTANCIA DE LAS ECUACIONES DIFERENCIALES EN LA INGENIERIA

AMBIENTAL

- Realización de proyectos ambientales

- Interpretación de los fenómenos de la naturaleza por medio de expresiones o modelos
matemáticos
- Tratamiento de aguas residuales
- Sistema de recolección y tratamiento de residuos solidos
- Estudios de contaminación
- Diagnóstico, evaluación y monitoreo de ecosistemas

LA ECUACION DE LAPLACE
En tres dimensiones, el problema consiste en hallar funciones reales doblemente diferenciables,
una función scriptstyle u de variables reales x, y, y z, tal que

En coordenadas cartesianas

En coordenadas cilíndricas,

En coordenadas esféricas,

Muchas veces se escribe de la siguiente manera:

Donde es el operador de Laplace o "laplaciano", que también se escribe como:

Donde es la divergencia, y es el gradiente o sino, algunas veces la notación puede ser:

Donde Delta también es el operador de Laplace.

Las soluciones de la ecuación de Laplace se denominan funciones armónicas.
Si del lado derecho de la igualdad se especifica una función, f(x, y, z), es decir, si la ecuación se
escribe como:

Entonces se tiene la "ecuación de Poisson", por lo que la ecuación de Laplace es un caso

particular de esta. La ecuación de Laplace también es un caso particular de la ecuación de
Helmholtz.
La ecuación de Laplace, así como también la ecuación de Poisson, son los ejemplos más
simples de ecuaciones en derivadas parciales elípticas.
Condiciones de contorno o frontera

Problema de Dirichlet
El problema de Dirichlet para la ecuación de Laplace consiste de hallar una solución en algún
dominio D tal que U sobre su contorno o frontera es igual a una función determinada:

Como el operador de Laplace aparece en la ecuación del calor, una interpretación física de este
problema es lo siguiente: fijar la temperatura sobre el contorno del dominio de acuerdo a una
especificación determinada de la condición de contorno. La temperatura fluye hasta que alcanza
un estado estacionario en el que dicha temperatura en cada punto del dominio no cambia más.
La distribución de la temperatura en el interior será entonces la solución correspondiente al
problema de Dirichlet.

Problema de Neumann
Las condiciones de contorno de Neumann para la ecuación de Laplace no especifica la función
u en sí mismo sobre el contorno , pero sí su derivada normal. Físicamente, esto corresponde
a la construcción de un potencial para un campo vectorial cuyo efecto es conocido en el
contorno de D:

Las soluciones de la ecuación de Laplace son funciones armónicas; son todas analíticas dentro
del dominio donde la ecuación se satisface. Si cuales quiera de dos funciones son soluciones a la
ecuación de Laplace (o de cualquier ecuación diferencial homogénea), su suma (o cualquier
combinación lineal) es también una solución. Esta propiedad, llamada principio de
superposición, es muy útil, por ejemplo, las soluciones de problemas complejos pueden
construirse simplemente sumando las soluciones.
Ecuación de Laplace en dos dimensiones
La ecuación de Laplace en dos variables independientes:

Funciones analíticas
Las partes reales e imaginarias de un función analítica en los complejos satisfacen la ecuación
de Laplace. Es decir, si , y si

Entonces la condición necesaria para que sea analítica es que se satisfagan las ecuaciones
de Cauchy-Riemann:

Donde ux es la primera derivada parcial de u con respecto a x.

La ecuación de Laplace para varphi implica que la condición de integrabilidad para psi se
satisface:

Asi puede definirse con una integral de línea. La condición de integrabilidad y el teorema de
Stokes implican que el valor de la integral de línea que conecta dos puntos es independiente del
camino. El par de soluciones resultante de la ecuación de Laplace se denominan funciones
armónicas conjugadas. Esta construcción solo es válida localmente, o siempre que el camino no
esté rodeando una singularidad. Por ejemplo, si r y θ son coordenadas polares y entonces una
función analítica correspondiente es

Sin embargo, el ángulo θ es univaluada solo en una región que no incluye el origen.
La estrecha relación entre la ecuación de Laplace y las funciones analíticas establece que
cualquier solución de la ecuación de Laplace tiene derivadas en todos los órdenes, y puede
expandirse en series de potencias, al menos dentro de un círculo que no incluya una
singularidad. Esto está en contraste con las soluciones de la ecuación de onda, que por lo
general tiene menor regularidad.
Hay una íntima conexión entre las series de potencias y las series de Fourier. Si expandimos una
función f en series de potencias dentro de un círculo de radio R, esto significa que con
coeficientes definidos adecuadamente cuyas parte real e imaginaria están dadas por:
Entonces

La cual es una serie de fourier de f.

Flujo de fluido
Sean las cantidades u y v las componentes horizontal y vertical del campo de velocidad del flujo
incompresible estacionario e irrotacional en dos dimensiones. LA condición de que el flujo es
incompresible es que

Y la condición de que el flujo es irrotatacional es que

Si definimos el diferencial de ψ como

Entonces la condición de incompressibilidad es la condición de integrabilidad para este

diferencial: la función resultante se llama función de corriente porque es constante a lo largo de
las líneas de flujo. Las primeras derivadas de ψ son

Y la condición de irrotacionabilidad establece que ψ satisface la ecuación de Laplace. La

función armónica φ que es el conjugado de ψ se denomina potencial de velocidad. Las
ecuaciones de Cauchy-Riemann establecen que

Así, para cada función analítica corresponde a un flijo de fluido incompresible estacionario e
irrotational en el plano. Las parte real es el potencial de velocidad, y la parte imaginaria es la
función de corriente.

Electrostática
De acuerdo a las ecuaciones de Maxwell, un campo eléctrico (u,v) en un espacio de dos
dimensiones que es independiente del tiempo satisface

donde rho es la densidad de carga. La primera ecuación de Maxwell es la condición de

integralidad para el diferencial

así el potencial eléctrico φ puede construirse para satisfacer

La segunda ecuación de Maxwell establece que

La cual es la ecuación de Poisson.

Es importante observar que la ecuación de Laplace puede usarse en problemas de tres
dimensiones en electroestática y flujo de fluido así como en dos dimensiones

Ecuación de Laplace en tres dimensiones

Solución fundamental
Una solución fundamental de la ecuación de Laplace satisface:

Donde la función delta de Dirac delta es una fuente unitaria concentrada en un punto. No es una
función en sí, sin embargo puede pensarse como el límite de funciones cuya integral sobre todo
el espacio es unitaria, y cuya región donde la función es distinta de cero es solo en un punto (ver
solución débil). Es común elegir una convención de signos diferente para esta ecuación, esto se
hace cuando se define la solución fundamental. Frecuentemente la elección de este signo es
conveniente para trabajar con un que es un operador positivo. Así la definición de la
solución fundamental implica que, si el laplaciano de u es integrado sobre cualquier volumen
que encierra el punto de la fuente, entonces:

La ecuación de Laplace no cambia bajo un cambio de coordenadas, y entonces podemos esperar

que la solución fundamental puede obtenerse entre soluciones que solo dependen de la distancia
r del punto de la fuente. Si elegimos el volumen de una bola de radio a alrededor del punto de la
fuente, entonce por el teorema de la divergencia de Gauss:

Entonces

Sobre una esfera de radio r que tiene como centro al punto de la fuente y por lo tanto

Un argumento similar muestra que en dos dimensiones:

Función de Green
Una función de Green es una solución fundamental que también satisface una condición
adecuada de contorno S de un volumen V. Por ejemplo, satisface

Ahora si u es cualquier solución de la ecuación de Poisson en V:

Y u toma valores de contorno g sobre S, entonces podemos aplicar la identidad de Green, (una
consecuencia del teorema de la divergencia) el cual satisface

Las notaciones un y Gn se refieren a derivadas normales a S. En vista de que las condiciones

satisfacen u y G, este resultado simplifica a

Así la función de Green describe la influencia en scriptstyle de f y g. Para el caso del interior de
una esfera de radio a, la función de Green puede obtenerse por medio de la reflexión el punto de
la fuente P a distancia ρ del centro de la esfera se refleja a lo largo de la línea radial al punto P'
que es en una distancia

Se observa que si P está dentro de la esfera, entonces P' estará fuera de la esfera. La función de
Green está dada entonces por

Donde R es la distancia al punto de la fuente P y R' es la distancia al punto reflejado P'. Una
consecuencia de esta expresión para la función de Green es la fórmula integral de Poisson. Sea
ρ, θ, y φ las componentes de coordenadas esféricas del punto de la fuente P. Aquí θ es el ángulo
con el eje vertical, la cual es contraria a la notación matemática estadounidense, pero cumple
con el estándar europeo y la práctica de la física . Entonces la solución de la ecuación de
Laplace equation dentro de la esfera está dado por
Donde

Una consecuencia simple de esta fórmula es que si u es una función armónica, enotnces el valor
de u en el dentro de la esfera es el valor medio de los valores sobre la esfera. Esta propiedad de
valor medio implica inmediatamente que funciones armónicas no constantes no pueden tomar
su valor máximo en un punto interior.
Electroestática
En el espacio libre la ecuación de laplace de cualquier potencial electroestático debe ser igual a
cero ya que rho (densidad de carga volumétrica) es cero en el espacio libre.
A partir del gradiente del potencial se obtiene el campo eléctrico

Tomando la divergencia del campo eléctrico se obtiene la ecuación de Poisson, que relaciona el
potencial eléctrico con la densidad de carga

En el caso particular del espacio libre ( ) la ecuación de Poisson se reduce a la de

Laplace.
Usando el teorema de la unicidad y mostrando que un potencial satisface la ecuación de Laplace
(la segunda derivada de V debería ser cero en el espacio libre) y el potencial tiene los valores
correctos en el contorno, el potencial entonces está unívocamente definido.
Un potencial que no satisface la ecuación de Laplace junto con la condición de contorno es un
potencial electroestático inválido.

En el valor de   en la región   es determinada por el valor de la función en el borde   
se les denomina condiciones de borde de Dirichlet que viene dado por:

En la figura, la región   está definida como: