Analisis Numerico

Ecuaciones Diferenciales en Derivadas
Parciales: La ecuación de Laplace ,Técnicas

de solución , Condiciones de frontera y
Volumen de control
 Introducción
 - Se denominan ecuaciones diferenciales parciales
(EDP) a aquellas ecuaciones que involucran derivadas
parciales de una función desconocida con dos o más
variables independientes.
Ejemplo:
 La ecuación
 La ecuación dede Laplace es
Laplace es una ecuación
una ecuación en
en derivadas
derivadas parciales
parciales de de
segundo
segundo orden
orden de
de tipo
tipo elíptico,
elíptico, que
que recibe
recibe ese
ese nombre
nombre en en honor
honor al
al
físico
físico yy matemático
matemático Pierre
Pierre –– Simón
Simón Laplace.
Laplace.
 La
 La ecuación
ecuación dede Laplace
Laplace aparece
aparece enen muchas
muchas otras
otras ramas
ramas de de la
la
Física
Física teórica
teórica como
como la Astronomía,
la Astronomía, la Electrostática,
la Electrostática, la
la Mecánica
Mecánica
de
de fluidos
fluidos o
o la
la Mecánica
Mecánica Cuántica.
Cuántica.
Ecuaciones Elípticas
 Para abordar la resolución numérica de las ecuaciones

elípticas utilizaremos como caso de estudio a la ecuación de
Laplace por ser utilizada en diversas áreas de ingeniería
donde se trata con problemas que involucran la determinación
de un potencial. Por ser un problema simple de plantear y
resolver utilizaremos el caso de flujo de calor en régimen
estacionario en una placa delgada. La expresión es:
Para abordar la solución numérica, trataremos a la placa como una malla de puntos discretos
(nodos) donde plantearemos la representación en diferencias finitas de la ecuación
diferencial, lo cual transforma a la EDP en una ecuación algebraica en diferencias. Utilizando
diferencias finitas centradas de segundo orden, entonces podemos escribir:
Reemplazando estas expresiones en la EDP, queda:

EJEMPLO:
La placa delgada de espesor ∆z. Se muestra un elemento, con
el cual se hace el balance de calor.
Donde
Ecuación
= flujo de calor en la dirección de la

qi
dimensión i [cal/(cm2 · s)], k = coeficiente de
difusividad térmica (cm2 /s), = densidad del
material (g/cm3 ), C = capacidad calorífica del
material [cal/(g · °C)] y T = temperatura (°C),
que se define como
Sustituyendo la ecuación, se obtiene
H = calor (cal) y V = volumen (cm3 )

 Técnica de solución :
a) La ecuación Laplaciana en Malla usada para la solución por diferencias

diferencias finitas de las EDP elípticas en dos variables
independientes, como la ecuación de Laplace.
 Las diferencias centrales basadas en el esquema de malla
 Sustituyendo estas expresiones en la ecuación

 En la malla cuadrada de la figura , ∆x = ∆y, y reagrupando
términos, la ecuación se convierte en
 En la malla cuadrada de la figura, ∆x = ∆y, y reagrupando términos, la ecuación se

convierte en
 Esta relación, que se satisface por todos los puntos interiores de la placa, se conoce como
ecuación Laplaciana en diferencias. Además, se deben especificar las condiciones de
frontera en los extremos de la placa para obtener una solución única. El caso más simple es
aquel donde la temperatura en la frontera es un valor fijo
 Sin embargo, = 75 y = 0, y, por lo tanto,

 Ecuaciones similares se pueden desarrollar para los otros puntos interiores. El
resultado es el siguiente conjunto de nueve ecuaciones simultáneas con nueve
incógnitas:
b) Método de Liebman
 Debido al gran número de ceros que aparecen en la matriz A (Ax=b), los métodos
aproximados representan un mejor procedimiento para obtener soluciones de EDP
elípticas. El método comúnmente empleado es el de Gauss-Seidel, (método de
Liebmann)
 Con esta técnica, la ecuación laplaciana en diferencias se
expresa como:
y se resuelve de manera iterativa para j = 1 hasta n e i = 1 hasta m
c) Método de variables secundarias
 Como la distribución de temperatura está descrita por la ecuación de Laplace, ésta se considera la
variable principal en el problema de la placa caliente. En este caso, así como en otros problemas
donde se tengan EDP, las variables secundarias también pueden ser importantes. De hecho, en
ciertos contextos de ingeniería, la variable secundaria puede realmente ser más importante. En la
placa caliente, una variable secundaria es el flujo de calor a través de la superficie de la placa. Esta
cantidad se calcula a partir de la ley de Fourier. Las aproximaciones por diferencias finitas centradas
para las primeras derivadas se sustituyen en la ecuación para obtener los siguientes valores del
flujo de calor en las dimensiones x y y:
 El flujo de calor resultante se calcula a partir de estas dos cantidades
mediante
donde la dirección de está dada por

 CONDICIONES EN LA FRONTERA
Debido a que está libre de complicaciones, la placa rectangular con

condiciones de frontera fijas representa un ideal para mostrar cómo
se resuelven numéricamente las EDP elípticas. Ahora veremos otro
problema que ampliará nuestras habilidades para abordar problemas
más realistas. Éste considera fronteras en donde se especifica la
derivada, y fronteras que tienen forma irregular.
 Condiciones con derivada en la frontera  Un nodo frontera (0, j) en

el extremo izquierdo de
La condición de frontera fija o de Dirichlet analizada hasta
una placa caliente. Para
ahora es uno de los diferentes tipos usados en las ecuaciones
aproximar la derivada
diferenciales parciales. Una alternativa común es el caso
normal al extremo (es
donde se da la derivada, que se conoce comúnmente como
decir, la derivada x), se
una condición de frontera de Neumann. En el problema de la
localiza un punto
placa caliente, esto corresponde a especificar el flujo de calor,
imaginario (–1, j) a una
más que la temperatura en la frontera. Un ejemplo es la
distancia ∆x más allá del
situación donde el extremo está aislado. En tal caso, la
extremo
derivada es cero
Observe que para esta ecuación se necesita un punto imaginario (–1, j) que esté fuera de la
placa. Aunque este punto exterior ficticio podría parecer que representa un problema,
realmente sirve para incorporar la derivada de la condición de frontera en el problema, lo cual
se logra representando la primera derivada en la dimensión x en (0, j) por la diferencia dividida
finita
Donde se puede despejar
Ahora se tiene una relación para T–1,j que incluye la derivada. Esta relación se
sustituye en la ecuación
 Fronteras irregulares
La figura es un sistema útil para ilustrar cómo se pueden tratar las fronteras no
rectangulares. Como se muestra, la frontera inferior izquierda de la placa es
circular. Observe que tenemos parámetros adicionales (a1 , a2 , b1 , b2 ) en cada
una de las longitudes que rodean al nodo. Por supuesto que, para la placa
mostrada en la figura , a2 = b2 = 1. Conservaremos estos parámetros en la
siguiente deducción, de tal modo que la ecuación resultante sea aplicable a
cualquier frontera irregular (y no sólo a la esquina inferior izquierda de una placa
caliente). Las primeras derivadas en la dimensión x se aproximan como sigue
 Malla de una placa caliente con

una frontera en forma
irregular. Observe cómo se
utilizan los coeficientes
ponderados al considerar el
espaciamiento no uniforme en
Las segundas derivadas se obtienen a partir de estas primeras
la cercanía de la frontera no
derivadas. Para la dimensión x, la segunda derivada es
rectangular.
Sustituyendo estas ecuaciones en la ecuación obtenemos
Sustituyendo las ecuaciones obtenemos;
Es posible desarrollar una ecuación similar en la dimensión y

d) Método del volumen de control
El método por diferencias finitas o series de Taylor divide al sistema en nodos

y escribe la ecuación diferencial parcial correspondiente para cada nodo. Las
aproximaciones por diferencias finitas, entonces, se sustituyen por las
derivadas para llevar las ecuaciones a una forma algebraica. Un
procedimiento de esto es bastante simple y directo con mallas rectangulares
y coeficientes constantes.
Sin embargo, este procedimiento se vuelve un más complicado cuando se
tiene derivadas como condiciones de fronteras con forma irregular. Veamos
un ejemplo de un sistema donde se presentan dificultades adicionales. La
siguiente placa está hecha de dos materiales diferentes y los espacios en la
malla también son diferentes. Además, la mitad de su extremo superior está
sujeta a transferencia de calor convectivo; mientras que la otra mitad está
aislada. Obtener las ecuaciones para algunos nodos, ej. el nodo (4, 2)
requeriría algunas deducciones adicionales, que van más allá de los métodos
vistos.
El método del volumen de control nos da
otra manera de aproximar numéricamente
las EDP, útil en casos como el de la figura
. El método se parece a la aproximación
por puntos, donde los puntos se
determinan a través del dominio. No
obstante, en lugar de aproximar la EDP en
un punto, la aproximación se aplica al
volumen que rodea el punto. En una malla
ortogonal, el volumen está formado por las
rectas perpendiculares que pasan por el
punto medio de cada línea que une nodos
adyacentes. Un balance de calor, en
estado estacionario, se obtiene después
para cada volumen de manera similar al
anterior(sumando y restando los flujos de
calor entrantes y salientes al volumen).
EJEMPLO:
- Se tiene una placa metálica que tiene una temperatura en extremos de los lados , , yy
 Entonces vamos a dividir en cuatro partes iguales

el eje x y el eje y nos queda :
100
°
𝑇 11 𝑇 12 𝑇 13 𝑇14 𝑇 15
100
°𝑇 2 1 𝑇 22 𝑇 23 𝑇 24 𝑇 𝟎°

25
𝟓𝟎 °
100
°𝑇 3 1 𝑇 32 𝑇 33 𝑇 34 𝑇 35
𝟎°

100 °𝑇 4 1 𝑇 42 𝑇 43 𝑇 44 𝑇 45

𝟎°

𝑇51 𝑇 52 𝑇 53 𝑇 54 𝑇 55
𝟎°

𝟓𝟎 ° 𝟓𝟎 ° 𝟓𝟎 °

Entonces:
+175=0
+75=0
+75=0
+100=0
=0
=0
+150=0
+50=0
+50=0

Analisis Numerico

Cargado por

Copyright:

Formatos disponibles

Analisis Numerico

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Analisis Numerico

Cargado por

Copyright:

Formatos disponibles

Ecuaciones Diferenciales en Derivadas

Parciales: La ecuación de Laplace ,Técnicas

 Para abordar la resolución numérica de las ecuaciones

Reemplazando estas expresiones en la EDP, queda:

= flujo de calor en la dirección de la

H = calor (cal) y V = volumen (cm3 )

a) La ecuación Laplaciana en Malla usada para la solución por diferencias

 Sustituyendo estas expresiones en la ecuación

 En la malla cuadrada de la figura, ∆x = ∆y, y reagrupando términos, la ecuación se

 Sin embargo, = 75 y = 0, y, por lo tanto,

y se resuelve de manera iterativa para j = 1 hasta n e i = 1 hasta m

c) Método de variables secundarias

donde la dirección de está dada por

Debido a que está libre de complicaciones, la placa rectangular con

 Condiciones con derivada en la frontera  Un nodo frontera (0, j) en

Donde se puede despejar

 Malla de una placa caliente con

Sustituyendo las ecuaciones obtenemos;

Es posible desarrollar una ecuación similar en la dimensión y

El método por diferencias finitas o series de Taylor divide al sistema en nodos

 Entonces vamos a dividir en cuatro partes iguales

𝟓𝟎 ° 𝟓𝟎 ° 𝟓𝟎 °

También podría gustarte