Partiendo con la ecuacion del calor con conductividad constante $k = k_0$ $$ \frac{\partial T}{\partial t} = k_0\frac{\partial^2T}{\partial x^2} $$
$\implies$ Problema adimensionalizado: $$ \dfrac{\partial T'}{\partial t'} = \dfrac{\partial^2 T'}{\partial x'^2} $$
Se define $r = \dfrac{\epsilon}{h^2}$ y se tiene:
$$T_j^{n+1} = rT_{j+1}^n + (1-2r)T^n_j+ rT^n_{j-1}$$Para analizar la estabilidad del algoritmo se asume una solucion de tipo $T^n_j = \sigma^ne^{ikx_j}$.
Reemplazando se tiene y usando que $x_j = x_0 + hj$:
$$\begin{matrix} \sigma^{n+1}e^{ikx_j} = r\sigma^ne^{ikx_{j+1}} + (1- 2r)\sigma^ne^{ikx_j} + r\sigma^ne^{ikx_{j-1}} \\ \\ \sigma^{n+1} = r\sigma^ne^{ikh} + (1-2r)\sigma^n+r\sigma^ne^{-ikh} \\ \\ \sigma = re^{ikh}+(1-2r)+re^{-ikh} \\ \\ = 2r\cos(kh)+(1-2r) \\ \\ \sigma = 1-2r(1-\cos(kh)) \end{matrix}$$Para que no diverja, $\sigma \in (-1, 1) \implies \vert\sigma\vert < 1$
$-1 < 1-2r(1-\cos(kh)) < 1$ a la vez. $\rightarrow -2r(1-\cos(kh)) > -2 \implies r(1-\cos(kh)) < 1$
En el peor caso $\cos(kh) = -1\quad r = \dfrac{\epsilon}{h^2} < \dfrac{1}{2}$
Si bien el algoritmo es estable, resulta ser caro para los pasos temporales y muy restrictivo para alanzar buena precision.
Renombrando $s = \dfrac{\epsilon}{2h^2}$
$$\begin{matrix} T_j^{n+1} - T^n_j = s\left(T^{n+1}_{j+1} - 2T^{n+1}_j + T^{n+1}_{j-1} + T^n_{j+1} - 2T^n_j + T^n_{j-1}\right) \\ \\ T_j^{n+1} - s\left(T^{n+1}_{j+1} - 2T^{n+1}_j + T^{n+1}_{j-1}\right) = T_j^n + s\left(T^n_{j+1} - 2T^n_j + T^n_{j-1}\right) \\ \\ \implies - sT^{n+1}_{j+1} + (2s+1)T^{n+1}_j - sT^{n+1}_{j-1} = sT^n_{j+1} + (1-2s)T^n_j + sT^n_{j-1} \end{matrix}$$Asumiendo una solucion con la forma $T = \sigma e^{ikx}$ y recordando que $x_j = x_0 + hj$
$$\begin{matrix} -s\sigma^{n+1}e^{ikx_{j+1}} + (2s+1)\sigma^{n+1}e^{ikx_j} - s\sigma^{n+1}e^{ikx_{j-1}} = s\sigma^{n}e^{ikx_{j+1}} + (1-2s)\sigma^{n}e^{ikx_{j}} + s\sigma^{n}e^{ikx_{j-1}} \\ \\ \implies -s\sigma e^{ikh} + (2s+1)\sigma - s\sigma e^{-ikh} = se^{ikh}+(1-2s)+se^{-ikh} \\ \\ \implies -2s\sigma\cos(kh) + (2s+1)\sigma = 2s\cos(kh)+(1-2s) \\ \\ \implies \sigma = \frac{1-2s(1-\cos(kh))}{1+2s(1-\cos(kh))} \implies \vert\sigma\vert < 1 \end{matrix}$$Siempre es incodicionalmente estable.
Esto se escribe como una matriz tri-diagonal: $$\left[\begin{matrix} -s & (2s+1) & -s & 0 & 0 & 0 & 0\\ 0 & -s & (2s+1) & -s & 0 & 0 & 0\\ 0 & 0 & -s & (2s+1) & -s & 0 & 0\\ 0 & 0 & 0 & -s & (2s+1) & -s & 0\\ 0 & 0 & 0 & 0 &-s & (2s+1) & -s \end{matrix}\right]\left(\begin{matrix} \vdots \\ T_{j-1} \\ T_j \\ t_{j+1} \\ \vdots \end{matrix}\right)^{n+1}$$
Un problema de tipo $A\vec{\phi} = \vec{b}$
In [ ]: