Como este blog destina-se a divulgar assuntos em diversos níveis, irei apresentar alguns posts relacionados a teoria das Equações Diferenciais Parciais. Veremos neste post algumas propriedades da equação N-dimensional do calor em domínios limitados. Seja
Considere o problema de autovalores:
![\begin{cases} -\triangle \phi_j = \lambda_j\phi_j,\quad \quad x \in \Omega\\\phi_j = 0 \quad \quad \text{sobre}\quad \partial \Omega\\\end{cases} [;\begin{cases} -\triangle \phi_j = \lambda_j\phi_j,\quad \quad x \in \Omega\\\phi_j = 0 \quad \quad \text{sobre}\quad \partial \Omega\\\end{cases};]](http://thewe.net/tex/%5Cbegin%7Bcases%7D%20-%5Ctriangle%20%5Cphi_j%20=%20%5Clambda_j%5Cphi_j,%5Cquad%20%5Cquad%20x%20%5Cin%20%5COmega%5C%5C%5Cphi_j%20=%200%20%5Cquad%20%5Cquad%20%5Ctext%7Bsobre%7D%5Cquad%20%5Cpartial%20%5COmega%5C%5C%5Cend%7Bcases%7D)
A teoria esprectal dos operadores compactos auto-adjuntos no espaço de Hilbert, garante que o problema de autovalores admite uma sequência crescente de autovalores positivos de multiplicidade finita tendendo para o infinito, isto é,
![0 \prec \lambda_1 \leq \lambda_2 \leq \lambda_3 \ldots \leq \lambda_n \leq \ldots \to +\infty [;0 \prec \lambda_1 \leq \lambda_2 \leq \lambda_3 \ldots \leq \lambda_n \leq \ldots \to +\infty;]](http://thewe.net/tex/0%20%5Cprec%20%5Clambda_1%20%5Cleq%20%5Clambda_2%20%5Cleq%20%5Clambda_3%20%5Cldots%20%5Cleq%20%5Clambda_n%20%5Cleq%20%5Cldots%20%5Cto%20+%5Cinfty)
Deste modo, as autofunções correspondentes
constitui uma base ortonormal para
. Sendo
um domínio limitado, as soluções do problema
podem ser desenvolvidas na base de autofunções acima, ou seja:
![u(x,t) = \sum_{k=1}^{\infty}c_ke^{-\lambda_k t}\phi_k(x) \quad \quad (2) [;u(x,t) = \sum_{k=1}^{\infty}c_ke^{-\lambda_k t}\phi_k(x) \quad \quad (2);]](http://thewe.net/tex/u%28x,t%29%20=%20%5Csum_%7Bk=1%7D%5E%7B%5Cinfty%7Dc_ke%5E%7B-%5Clambda_k%20t%7D%5Cphi_k%28x%29%20%5Cquad%20%5Cquad%20%282%29)
Note que
![u_0(x) = u(x,0) = \sum_{k=1}^{\infty}c_k\phi_k(x) [;u_0(x) = u(x,0) = \sum_{k=1}^{\infty}c_k\phi_k(x);]](http://thewe.net/tex/u_0%28x%29%20=%20u%28x,0%29%20=%20%5Csum_%7Bk=1%7D%5E%7B%5Cinfty%7Dc_k%5Cphi_k%28x%29)
Afirmação: Os coeficientes
são dados por
![c_k = \int_{\Omega}u_0(x)\phi_k(x)dx [;c_k = \int_{\Omega}u_0(x)\phi_k(x)dx;]](http://thewe.net/tex/c_k%20=%20%5Cint_%7B%5COmega%7Du_0%28x%29%5Cphi_k%28x%29dx)
De fato,
![\prec u_0(x),\phi_j(x)\succ \ = \ \prec \sum_{k=1}^{\infty}c_k\phi_k(x),\phi_j(x) \succ \ = \ \sum_{k=1}^{\infty}c_k\prec \phi_k(x),\phi_j(x)\succ [;\prec u_0(x),\phi_j(x)\succ \ = \ \prec \sum_{k=1}^{\infty}c_k\phi_k(x),\phi_j(x) \succ \ = \ \sum_{k=1}^{\infty}c_k\prec \phi_k(x),\phi_j(x)\succ;]](http://thewe.net/tex/%5Cprec%20u_0%28x%29,%5Cphi_j%28x%29%5Csucc%20%5C%20=%20%5C%20%5Cprec%20%5Csum_%7Bk=1%7D%5E%7B%5Cinfty%7Dc_k%5Cphi_k%28x%29,%5Cphi_j%28x%29%20%5Csucc%20%5C%20=%20%5C%20%5Csum_%7Bk=1%7D%5E%7B%5Cinfty%7Dc_k%5Cprec%20%5Cphi_k%28x%29,%5Cphi_j%28x%29%5Csucc)
Usando a ortogonalidade das autofunções
, segue que
![\int_{\Omega}u_0(x)\phi_j(x)dx = \ \prec u_0(x),\phi_j(x)\succ \ = \ c_j\prec \phi_j(x),\phi_j(x) \succ \ = c_j [;\int_{\Omega}u_0(x)\phi_j(x)dx = \ \prec u_0(x),\phi_j(x)\succ \ = \ c_j\prec \phi_j(x),\phi_j(x) \succ \ = c_j;]](http://thewe.net/tex/%5Cint_%7B%5COmega%7Du_0%28x%29%5Cphi_j%28x%29dx%20=%20%5C%20%5Cprec%20u_0%28x%29,%5Cphi_j%28x%29%5Csucc%20%5C%20=%20%5C%20c_j%5Cprec%20%5Cphi_j%28x%29,%5Cphi_j%28x%29%20%5Csucc%20%5C%20=%20c_j)
Usando a ortogonalidade das autofunções
Propriedade 1:
![\mid\mid u(\cdot, t)\mid \mid_{L^2(\Omega)}^{2} = \sum_{k = 1}^{\infty} c_{k}^{2}e^{-2\lambda_k t} [;\mid\mid u(\cdot, t)\mid \mid_{L^2(\Omega)}^{2} = \sum_{k = 1}^{\infty} c_{k}^{2}e^{-2\lambda_k t};]](http://thewe.net/tex/%5Cmid%5Cmid%20u%28%5Ccdot,%20t%29%5Cmid%20%5Cmid_%7BL%5E2%28%5COmega%29%7D%5E%7B2%7D%20=%20%5Csum_%7Bk%20=%201%7D%5E%7B%5Cinfty%7D%20c_%7Bk%7D%5E%7B2%7De%5E%7B-2%5Clambda_k%20t%7D)
Demonstração: Usando a definição de norma no espaço
, temos:
![\int_{\Omega}u(x,t)^2 dx = \int_{\Omega}\biggl (\sum_{j=1}^{\infty}c_je^{-\lambda_j t}\phi_j(x) \biggr)\biggl(\sum_{k=1}^{\infty}c_ke^{-\lambda_k t}\phi_k(x) \biggr)dx [;\int_{\Omega}u(x,t)^2 dx = \int_{\Omega}\biggl (\sum_{j=1}^{\infty}c_je^{-\lambda_j t}\phi_j(x) \biggr)\biggl(\sum_{k=1}^{\infty}c_ke^{-\lambda_k t}\phi_k(x) \biggr)dx;]](http://thewe.net/tex/%5Cint_%7B%5COmega%7Du%28x,t%29%5E2%20dx%20=%20%5Cint_%7B%5COmega%7D%5Cbiggl%20%28%5Csum_%7Bj=1%7D%5E%7B%5Cinfty%7Dc_je%5E%7B-%5Clambda_j%20t%7D%5Cphi_j%28x%29%20%5Cbiggr%29%5Cbiggl%28%5Csum_%7Bk=1%7D%5E%7B%5Cinfty%7Dc_ke%5E%7B-%5Clambda_k%20t%7D%5Cphi_k%28x%29%20%5Cbiggr%29dx)
ou seja,
![\mid\mid u(\cdot, t)\mid \mid_{L^2(\Omega)}^{2} = \sum_{j=1}^{\infty}\sum_{k=1}^{\infty} c_jc_k e^{-\lambda_j t - \lambda_k t}\int_{\Omega}\phi_j(x)\phi_k(x)dx [;\mid\mid u(\cdot, t)\mid \mid_{L^2(\Omega)}^{2} = \sum_{j=1}^{\infty}\sum_{k=1}^{\infty} c_jc_k e^{-\lambda_j t - \lambda_k t}\int_{\Omega}\phi_j(x)\phi_k(x)dx;]](http://thewe.net/tex/%5Cmid%5Cmid%20u%28%5Ccdot,%20t%29%5Cmid%20%5Cmid_%7BL%5E2%28%5COmega%29%7D%5E%7B2%7D%20=%20%5Csum_%7Bj=1%7D%5E%7B%5Cinfty%7D%5Csum_%7Bk=1%7D%5E%7B%5Cinfty%7D%20c_jc_k%20e%5E%7B-%5Clambda_j%20t%20-%20%5Clambda_k%20t%7D%5Cint_%7B%5COmega%7D%5Cphi_j%28x%29%5Cphi_k%28x%29dx)
Usando a ortogonalidade das autofunções, segue que
![\mid\mid u(\cdot, t)\mid \mid_{L^2(\Omega)}^{2} = \sum_{k=1}^{\infty} c_{k}^{2} e^{-2\lambda_k t} [;\mid\mid u(\cdot, t)\mid \mid_{L^2(\Omega)}^{2} = \sum_{k=1}^{\infty} c_{k}^{2} e^{-2\lambda_k t};]](http://thewe.net/tex/%5Cmid%5Cmid%20u%28%5Ccdot,%20t%29%5Cmid%20%5Cmid_%7BL%5E2%28%5COmega%29%7D%5E%7B2%7D%20=%20%5Csum_%7Bk=1%7D%5E%7B%5Cinfty%7D%20c_%7Bk%7D%5E%7B2%7D%20e%5E%7B-2%5Clambda_k%20t%7D)
ou seja,
Usando a ortogonalidade das autofunções, segue que
temos o resultado desejado. Uma consequência imediata desta propriedade é que
![\mid\mid u_0(x)\mid \mid_{L^2(\Omega)}^{2} = \sum_{k=1}^{\infty} c_{k}^{2} [;\mid\mid u_0(x)\mid \mid_{L^2(\Omega)}^{2} = \sum_{k=1}^{\infty} c_{k}^{2};]](http://thewe.net/tex/%5Cmid%5Cmid%20u_0%28x%29%5Cmid%20%5Cmid_%7BL%5E2%28%5COmega%29%7D%5E%7B2%7D%20=%20%5Csum_%7Bk=1%7D%5E%7B%5Cinfty%7D%20c_%7Bk%7D%5E%7B2%7D)
Propriedade 2: (Dissipação de Energia) Para todo
,
![E(t) \leq E(0)e^{-2\lambda_1 t} [;E(t) \leq E(0)e^{-2\lambda_1 t};]](http://thewe.net/tex/E%28t%29%20%5Cleq%20E%280%29e%5E%7B-2%5Clambda_1%20t%7D)
![E(t) = \int_{\Omega}u(x,t)^2dx [;E(t) = \int_{\Omega}u(x,t)^2dx;]](http://thewe.net/tex/E%28t%29%20=%20%5Cint_%7B%5COmega%7Du%28x,t%29%5E2dx)
Propriedade 2: (Dissipação de Energia) Para todo
onde
é o primeiro autovalor e
Demonstração: Multiplicando a equação do calor por
, temos:
![(u_t - \triangle u)u = 0 \quad \quad \Rightarrow \quad \quad uu_t - u\Delta u = 0 [;(u_t - \triangle u)u = 0 \quad \quad \Rightarrow \quad \quad uu_t - u\Delta u = 0;]](http://thewe.net/tex/%28u_t%20-%20%5Ctriangle%20u%29u%20=%200%20%5Cquad%20%5Cquad%20%5CRightarrow%20%5Cquad%20%5Cquad%20uu_t%20-%20u%5CDelta%20u%20=%200)
ou seja:
![\frac{1}{2}\frac{\partial}{\partial t}(u^2) - u\Delta u = 0 \quad \quad (3) [;\frac{1}{2}\frac{\partial}{\partial t}(u^2) - u\Delta u = 0 \quad \quad (3);]](http://thewe.net/tex/%5Cfrac%7B1%7D%7B2%7D%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%20t%7D%28u%5E2%29%20-%20u%5CDelta%20u%20=%200%20%5Cquad%20%5Cquad%20%283%29)
Integrando a expressão
em relação a
sobre
, temos:
![\frac{1}{2}\int_{\Omega}\frac{\partial}{\partial t}(u^2)dx = \int_{\Omega}u\Delta udx \quad \quad (4) [;\frac{1}{2}\int_{\Omega}\frac{\partial}{\partial t}(u^2)dx = \int_{\Omega}u\Delta udx \quad \quad (4);]](http://thewe.net/tex/%5Cfrac%7B1%7D%7B2%7D%5Cint_%7B%5COmega%7D%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%20t%7D%28u%5E2%29dx%20=%20%5Cint_%7B%5COmega%7Du%5CDelta%20udx%20%5Cquad%20%5Cquad%20%284%29)
ou seja:
Integrando a expressão
Usando o teorema de Leibniz sobre a derivação sob o sinal de integração no lado esquerdo e a identidade de Green no lado direito de
, segue que
![\frac{1}{2}\frac{d}{dt}\biggl(\int_{\Omega}u(x,t)^2 dx\biggr) = \int_{\partial \Omega}u\frac{\partial u}{\partial n}dS - \int_{\Omega}\mid \nabla u \mid^2dx [;\frac{1}{2}\frac{d}{dt}\biggl(\int_{\Omega}u(x,t)^2 dx\biggr) = \int_{\partial \Omega}u\frac{\partial u}{\partial n}dS - \int_{\Omega}\mid \nabla u \mid^2dx;]](http://thewe.net/tex/%5Cfrac%7B1%7D%7B2%7D%5Cfrac%7Bd%7D%7Bdt%7D%5Cbiggl%28%5Cint_%7B%5COmega%7Du%28x,t%29%5E2%20dx%5Cbiggr%29%20=%20%5Cint_%7B%5Cpartial%20%5COmega%7Du%5Cfrac%7B%5Cpartial%20u%7D%7B%5Cpartial%20n%7DdS%20-%20%5Cint_%7B%5COmega%7D%5Cmid%20%5Cnabla%20u%20%5Cmid%5E2dx)
Demonstração: Basta fazer
na propriedade anterior.
Neste post apareceu desigualdades, teoremas e propriedades que foram apenas enunciadas, mas prometo apresentá-las integralmente em futuros posts.
Por
,
sobre
, de modo que a primeira integral no segundo membro é nula. Assim,
![\frac{1}{2}\frac{d}{dt}E(t) = - \int_{\Omega}\mid \nabla \mid^2 dx \leq -\lambda_1 E(t) [;\frac{1}{2}\frac{d}{dt}E(t) = - \int_{\Omega}\mid \nabla \mid^2 dx \leq -\lambda_1 E(t);]](http://thewe.net/tex/%5Cfrac%7B1%7D%7B2%7D%5Cfrac%7Bd%7D%7Bdt%7DE%28t%29%20=%20-%20%5Cint_%7B%5COmega%7D%5Cmid%20%5Cnabla%20%5Cmid%5E2%20dx%20%5Cleq%20-%5Clambda_1%20E%28t%29)
pela desigualdade Poincaré. Assim,
![E(t) \leq E(0)e^{-2\lambda_1t} [;E(t) \leq E(0)e^{-2\lambda_1t};]](http://thewe.net/tex/E%28t%29%20%5Cleq%20E%280%29e%5E%7B-2%5Clambda_1t%7D)
pela desigualdade Poincaré. Assim,
Corolário: (Monotonicidade):
![\mid \mid u(\cdot, t)\mid \mid_{L^2(\Omega)} \leq \mid \mid u_0(x)\mid \mid_{L^2(\Omega)} [;\mid \mid u(\cdot, t)\mid \mid_{L^2(\Omega)} \leq \mid \mid u_0(x)\mid \mid_{L^2(\Omega)};]](http://thewe.net/tex/%5Cmid%20%5Cmid%20u%28%5Ccdot,%20t%29%5Cmid%20%5Cmid_%7BL%5E2%28%5COmega%29%7D%20%5Cleq%20%5Cmid%20%5Cmid%20u_0%28x%29%5Cmid%20%5Cmid_%7BL%5E2%28%5COmega%29%7D)
Demonstração: Basta fazer
Neste post apareceu desigualdades, teoremas e propriedades que foram apenas enunciadas, mas prometo apresentá-las integralmente em futuros posts.


