Capítulo 12
Matlab
SOLUCIONES NUMÉRICAS PARA ECUACIONES DIFERENCIALES PARCIALES
Subcapítulo 1
Problema 1
Aproximar la solución a la ecuación diferencial parcial elíptica
1 ´
Utilice h k ycomparelosresultadosconlasoluci o nreal u ( x, y ) ( x y ) 2
2
Problema 2
Utilice el algoritmo 12.1 para aproximar la solución a la ecuación diferencial parcial elíptica
1 ´
Utilice h k ycomparelosresultadosconlasoluci o nrealu ( x, y ) ln( x 2 y 2 )
3
Problema 8
Una placa de plata rectangular de 6 cm por 5 cm genera calor uniformemente en cada punto a
la velocidad q = 1,5 cal /cm3 · S.
Sea x la distancia a lo largo del borde de la placa de 6 cm de longitud e y sea la distancia a lo
largo del borde de la placa de 5 cm de longitud. Suponga que la temperatura u a lo largo de los
bordes es mantenido a las siguientes temperaturas:
donde el origen se encuentra en una esquina de la placa con coordenadas (0, 0) y los bordes se
encuentran a lo largo del positivo eje x e y. La temperatura de estado estable u = u (x, y)
satisface la ecuación de Poisson:
donde K, la conductividad térmica, es 1.04 cal/cm·grados·s. Calcule la temperatura u (x, y)
usando Algoritmo 12.1 con h = 0.4 y k = 1/3.
Subcapítulo 2
Problema 1
Aproxime la solución de la siguiente ecuación diferencial parcial usando la Diferencia hacia
atrás método.
Use m = 4, T = 0.1 y N = 2, y compare sus resultados con la solución real
2
t
4
u ( x, t ) e sin x
2
Problema 2
Aproxime la solución de la siguiente ecuación diferencial parcial usando la Diferencia hacia
atrás método.
Use m = 3, T = 0.1 y N = 2, y compare sus resultados con la solución real
2
t
4
u ( x, t ) 2e sin 2 x
Problema 16
Utilice los resultados del ejercicio 15 para aproximar la solución a
con h = 0,1 y k = 0,01. Compare su respuesta en t = 0.25 con la solución real
\[u(x,t)=2{{e}^{-{{\pi }^{2}}t}}\sin \pi x+x\left( 1-x \right)\]
Problema 18
La temperatura u (x, t) de una varilla larga y delgada de sección transversal constante y
conducción homogénea. el material se rige por la ecuación de calor unidimensional. Si se
genera calor en el material, Por ejemplo, por resistencia a la corriente o reacción nuclear, la
ecuación de calor se convierte en
donde l es la longitud, ρ es la densidad, C es el calor específico y K es la difusividad térmica de
la varilla. La función r = r (x, t, u) representa el calor generado por unidad de volumen. Suponer
que
l 1.5cm, K 1.04cal / [Link].s, 10.6 g / cm3 ,
c 0.056cal / g .deg, y r ( x, t , u ) 5.0cal / cm 3 .s
Si los extremos de la varilla se mantienen a 0 ° C, entonces
\[u(0,t)=u(l,t)=0,\,\,\,\,t>0\]
Suponga que la distribución de temperatura inicial está dada por
\[u(x,0)=\sin \frac{\pi x}{l},\,\,\,\,0\le x\le l\]
Utilice los resultados del ejercicio 15 para aproximar la distribución de temperatura con h =
0,15 y k = 0.0225.
Subcapítulo 3
Problema 1
Aproxima la solución a la ecuación de onda
usando el algoritmo de diferencia finita 12.4 con m 4, N 4 y T 1.0.
´
Compare sus resultados en t 1.0 a la soluci o n realu ( x , t ) cos t sin x
Problema 2
Aproxima la solución a la ecuación de onda
usando el algoritmo de diferencia finita 12.4 con m 4, N 4 y T 0.5.
´
Compare sus resultados en t 0.5 a la soluci o n real u ( x, t ) sin t sin 4 x
Problema 3
Aproxima la solución a la ecuación de onda
\[\begin{align}
& usando\text{ }el\text{ }algoritmo\text{ }de\text{ }diferencia\text{ }finita\text{ }con\text{ h }
=\frac{\pi }{10}\,\,y\,\,k=0.05,\,\, \\
& con\text{ h}=\frac{\pi }{20}\,\,y\,\,k=0.1\,\,y\,\,\,\,\text{h}=\frac{\pi }{20}\,\,y\,\,k=0.0.5 \\
\end{align}\]
\[Compare\text{ }sus\text{ }resultados\text{ }en\text{ }t\text{ }=\text{ 0}\text{.5 }a\text{ }
la\text{ }soluci\acute{o}n\text{ }real\,\,u(x,t)=\,\cos t\,\sin x\]