Math ClubMath Club
v1 · padrão canônico

Lição 23 — Equação do Calor em 1D — Separação de Variáveis

Equação do calor 1D: método de separação de variáveis, problema de Sturm-Liouville, série de Fourier para condições iniciais gerais, decaimento exponencial dos modos.

Used in: Cálculo 3 — Unidade 3 · USP MAT0216 · ITA MA-013

ut=κuxx;u(x,t)=n=1bneκ(nπ/L)2tsin ⁣(nπxL);bn=2L0Lf(x)sin ⁣(nπxL)dxu_t = \kappa u_{xx};\quad u(x,t)=\sum_{n=1}^{\infty}b_n e^{-\kappa(n\pi/L)^2 t}\sin\!\left(\frac{n\pi x}{L}\right);\quad b_n = \frac{2}{L}\int_0^L f(x)\sin\!\left(\frac{n\pi x}{L}\right)dx

A solução da equação do calor com extremidades fixas é uma série de modos de Fourier, cada um decaindo exponencialmente. Os coeficientes bnb_n são obtidos pelo produto interno com a condição inicial f(x)f(x).

Choose your door

Rigorous notation, full derivation, hypotheses

Método de separação de variáveis

Problema modelo

ut=κuxx,0<x<L,  t>0u_t = \kappa u_{xx},\quad 0<x<L,\;t>0 u(0,t)=u(L,t)=0,u(x,0)=f(x)u(0,t) = u(L,t) = 0,\quad u(x,0) = f(x)

Passo 1: Separação

Tente u(x,t)=X(x)T(t)u(x,t) = X(x)T(t):

XT=κXTTκT=XX=λX T' = \kappa X'' T \Rightarrow \frac{T'}{\kappa T} = \frac{X''}{X} = -\lambda

(constante, pois lado esquerdo depende só de tt e lado direito só de xx).

Passo 2: Problema de Sturm-Liouville

X+λX=0,X(0)=X(L)=0X'' + \lambda X = 0,\quad X(0) = X(L) = 0

Soluções não-triviais existem apenas para autovalores λn=(nπ/L)2\lambda_n = (n\pi/L)^2 com autofunções Xn=sin(nπx/L)X_n = \sin(n\pi x/L), n=1,2,3,n = 1,2,3,\ldots

Passo 3: Equação temporal

Tn=κλnTnTn(t)=eκλnt=eκ(nπ/L)2tT'_n = -\kappa\lambda_n T_n \Rightarrow T_n(t) = e^{-\kappa\lambda_n t} = e^{-\kappa(n\pi/L)^2 t}

Passo 4: Superposição e coeficientes de Fourier

u(x,t)=n=1bnsin ⁣(nπxL)eκ(nπ/L)2tu(x,t) = \sum_{n=1}^\infty b_n \sin\!\left(\frac{n\pi x}{L}\right)e^{-\kappa(n\pi/L)^2 t}

Com u(x,0)=f(x)u(x,0) = f(x):

bn=2L0Lf(x)sin ⁣(nπxL)dxb_n = \frac{2}{L}\int_0^L f(x)\sin\!\left(\frac{n\pi x}{L}\right)dx

Comportamento assintótico

Para tt\to\infty: u0u\to 0 (esfriamento). O modo dominante é n=1n=1: ub1eκπ2t/L2sin(πx/L)u\approx b_1 e^{-\kappa\pi^2 t/L^2}\sin(\pi x/L).

Tempo de decaimento: τ=L2/(κπ2)\tau = L^2/(\kappa\pi^2) — proporcional a L2L^2 e inversamente a κ\kappa.

Exemplos resolvidos


Exercícios

Bloco A — Separação de variáveis

Bloco B — Aplicações

Bloco C — Nível ITA/USP

To continue

  • Strauss, W.A. — Partial Differential Equations, §4.1–4.3 — Wiley, 2ª ed.
  • Kreyszig, E. — Advanced Engineering Mathematics, §21.3 — Wiley, 10ª ed.
  • Churchill, R.V. & Brown, J.W. — Fourier Series and Boundary Value Problems, §2.3–2.5 — McGraw-Hill, 8ª ed.
  • Iório, V. — EDP — Um Curso de Graduação, §4 — IMPA, 2ª ed.
  • Evans, L.C. — Partial Differential Equations, §2.3 — AMS, 2ª ed.

Updated on 2026-05-28 · Author(s): Clube da Matemática

Found an error? Open an issue on GitHub or submit a PR — open source forever.