“APLICACIÓN DE DIFERENCIAS FINITAS EN LA
SOLUCIÓN DE LA ECUACIÓN DE LAPLACE
CONSIDERANDO EL MÉTODO DE LIEBMANN
(GAUSS-SEIDEL), PARA LA DISTRIBUCIÓN DE
TEMPERATURA
Oro EN UNA PLACA CALENTADA”
MÉTODOS NUMÉRICOS PARA ING.
GRUPO:
• CRUZADO ARENAS, Jhoan
•DIAZ PINO, Oswaldo David
•GIL QUIÑONES, Jaime Alberto
•HERNÁNDEZ RODRÍGUEZ, Jhunior
•ZAVALETA RODRIGUEZ, Keler
RESUMEN
El presente informe se dedica a soluciones por diferencias finitas de
ecuaciones elípticas.
Antes de poner en práctica los métodos, deducimos la ecuación de
Laplace para el problema físico de la distribución de temperatura en una
placa calentada.
Después, se describe un procedimiento estándar de solución: el método
de Liebmann (Gauss-Seidel).
Ilustraremos cómo se utiliza dicho procedimiento para calcular la
distribución de la variable escalar principal, la temperatura, así como la
de una variable vectorial secundaria: el flujo de calor.
PROBLEMA
¿Cómo usar el programa computacional de matemática
MATLAB en la solución de una ecuación diferencial
parcial elíptica sobre problemas estacionarios de
conductividad térmica mediante la ecuación de Laplace –
Método de Liebmann (Gauss-Seidel)?
OBJETIVOS
Usar el programa computacional técnico
MATLAB en la solución de ecuaciones
diferenciales parciales elípticas sobre
GENERAL conductividad térmica usando la ecuación de
Laplace considerando el método de Liebmann
(Gauss-Seidel).
- Analizar la deducción de la ecuación
diferencial parcial de tipo elíptico sobre
problemas estacionarios de conductividad
térmica.
ESPECÍFICOS
- Resolver la ecuación diferencial parcial
mediante el método en diferencias finitas
usando el programa computacional MATLAB.
FUNDAMENTO TEÓRICO
LA ECUACIÓN DE LAPLACE: La ecuación de Laplace se utiliza para
modelar diversos problemas que tienen que ver con el potencial de una
variable desconocida.
Se muestra un elemento sobre
la cara de una placa rectangular
delgada de espesor ∆z. La placa
está totalmente aislada excepto
en sus extremos, donde la
temperatura puede ajustarse a
un nivel preestablecido. El
aislamiento y el espesor de la
placa permiten que la
transferencia de calor esté
limitada solamente a las
dimensiones x y y.
En estado estacionario, el flujo de calor hacia el elemento en una
unidad de tiempo ∆t debe ser igual al flujo de salida, es decir
donde q(x) y q(y) = los flujos de calor en x y y, respectivamente [cal/
(cm2 · s)]. Dividiendo entre ∆z y ∆t, y reagrupando términos, se
obtiene:
Multiplicando el primer término por ∆x/∆x, y el segundo por ∆y/∆y se obtiene
Dividiendo entre ∆x ∆y, y tomando el límite, se llega a
TÉCNICA DE SOLUCIÓN
SOLUCIÓN DEL PROBLEMA
Temperatura de una placa calentada con condiciones de frontera
fijas.
Con el método de Liebmann (Gauss-Seidel) calcular la temperatura de
la placa delgada sometida en los bordes a condiciones.
Considerar ɛs = 0.1%.
Ta= 60 Tb= 60
Tc= 50 Td=70
L1= 5 L2= 5
∆x= 0.5 y ∆y= 0.5
ALGORITMO
Paso 1: Ingresar Ta, Tb, Tc, Td, L1, L2, h1 y h2
Paso 2: Calcular n=L1/h1, m=L2/h2
Paso 3: Ingresamos máximo de iteraciones
Paso 4: Ingresamos “e” (error de truncamiento relativo requerido)
Paso 5: Incorporando condiciones de contorno
u (i,1) = Tc , i = 1: n+2
u (i, m+2) = Td , i = 1: n+2
u (1, j) = Ta , i = 1: m+2
u (n+2, j) = Tb , i = 1: m+2
Paso 6: Incorporando valor inicial interior promedio de los bordes
p= 0.25*(Ta+Tb+Tc+Td) ; u (i, j) = p
Para i = 2: n-1
j = 2: m-1
Paso 7: Generando condición de conteo de iteraciones y señal de
convergencia
k= 0 , conv = 0
si k< máximo de iteraciones y conv ==0 , k = k+1 , t=u
Paso 8: Generando aproximaciones
para i = 2:n+1
j = 2: m+1
u (i, j) = 0.25*(u(i-1, j)+u(i+1, j)+u(i, j+1)+u(i, j-1))
Paso 9: Imprimir u
Paso 10: Graficar u
SCRIPT
RESULTADOS
Solución numérica calculada en los puntos de la malla
Representación gráfica en tres dimensiones de la solución calculada
CONCLUSIONES
- Se logró resolver el problema de conductividad térmica aplicando
diferencias finitas en la ecuación de Laplace, usando el método de
Liebmann (Gauss-Seidel) en el programa computacional MATLAB.
- Se pudo analizar y deducir la ecuación diferencial de tipo elíptico
sobre problemas estacionarios de conductividad térmica.
- Se pudo obtener información y adecuarnos a la técnica de solución
para el método de diferencias finitas, en la ecuación de Laplace.
- Se logró hacer el script que sirvió para solucionar el problema de
conductividad térmica y también se hizo el algoritmo respectivo.
¡GRACIAS!
