Equação de Helmholtz

Predefinição:Mais notas A equação de Helmholtz é um tipo de equação diferencial parcial que é expressa da seguinte forma:

${\mathrm{\nabla }}^{2}A+k^{2}A=0$

onde ∇² é o Laplaciano, k é o número de onda, e A é a amplitude.

A equação, que recebeu o nome de Hermann von Helmholtz, surge em vários domínios da física e engenharia, tipicamente para descrever fenómenos físicos que são dependentes do tempo. Ela corresponde a um caso geral da Equação de Laplace.

Motivação e usos

Predefinição:Sem-fontes A equação de Helmholtz frequentemente surge no estudo de problemas físicos envolvendo equação diferencial parcial (EDP) no espaço e no tempo. A equação de Helmholtz, que representa uma forma independente do tempo da equação de onda, resulta da aplicação da técnica de separação de variáveis para reduzir a complexidade da análise.

Por exemplo, considere a equação de onda $$({\mathrm{\nabla }}^2-\frac{1}{c^2}\frac{{\partial }^2}{\partial t^2})u(𝐫,t)=0.$$

A separação de variáveis começa assumindo que a função de onda Predefinição:Math é separável: $$u(𝐫,t)=A(𝐫)T(t).$$

Substituindo esta forma na equação de onda e simplificando, obtemos a seguinte equação: $$\frac{{\mathrm{\nabla }}^{2}A}{A}=\frac{1}{c^{2}T}\frac{{\mathrm{d}}^{2}T}{\mathrm{d}{t}^2}.$$

Note que a expressão do lado esquerdo depende apenas de Predefinição:Math, enquanto que a expressão do lado direito depende apenas de Predefinição:Mvar. Como resultado, esta equação é válida no caso geral se, e somente se, ambos os lados da equação forem iguais ao mesmo valor constante. Este argumento é fundamental na técnica de resolver equações diferenciais parciais lineares por separação de variáveis. Desta observação, obtemos duas equações, uma para Predefinição:Math, e outra para Predefinição:Math $$\frac{{\mathrm{\nabla }}^{2}A}{A}=-k^2$$ $$\frac{1}{c^{2}T}\frac{{\mathrm{d}}^{2}T}{\mathrm{d}{t}^2}=-k^2,$$

onde escolhemos, sem perda de generalidade, a expressão Predefinição:Math para o valor da constante. (É igualmente válido usar qualquer constante Predefinição:Mvar como constante de separação; Predefinição:Math é escolhida apenas para conveniência nas soluções resultantes.)

Rearranjando a primeira equação, obtemos a equação de Helmholtz (homogénea): $${\mathrm{\nabla }}^{2}A+k^{2}A=({\mathrm{\nabla }}^2+k^2)A=0.$$

Da mesma forma, após fazer a substituição Predefinição:Math, onde Predefinição:Mvar é o número de onda e Predefinição:Mvar é a frequência angular (assumindo um campo monocromático), a segunda equação torna-se

$$\frac{{\mathrm{d}}^{2}T}{\mathrm{d}{t}^2}+{\omega }^{2}T=(\frac{{\mathrm{d}}^2}{\mathrm{d}{t}^2}+{\omega }^2)T=0.$$

Agora temos a equação de Helmholtz para a variável espacial Predefinição:Math e uma equação diferencial ordinária de segunda ordem no tempo. A solução no tempo será uma combinação linear de funções seno e cosseno, cuja forma exata é determinada pelas condições iniciais, enquanto a forma da solução no espaço dependerá das condições de fronteira. Alternativamente, transformadas integrais, como a Laplace ou transformada de Fourier, são frequentemente usadas para transformar uma equação diferencial parcial hiperbólica numa forma da equação de Helmholtz.

Devido à sua relação com a equação de onda, a equação de Helmholtz surge em problemas em áreas da física como o estudo da radiação eletromagnética, sismologia e acústica.

Resolução da equação de Helmholtz usando separação de variáveis

A solução para a equação de Helmholtz espacial: $${\mathrm{\nabla }}^{2}A=-k^{2}A$$ pode ser obtida para geometrias simples usando a separação de variáveis.

Membrana vibrante

Predefinição:Sem-fontes O análogo bidimensional da corda vibrante é a membrana vibrante, com as bordas fixadas para serem imóveis. A equação de Helmholtz foi resolvida para muitas formas básicas no século XIX: a membrana retangular por Siméon Denis Poisson em 1829, o triângulo equilátero por Gabriel Lamé em 1852 e a membrana circular por Alfred Clebsch em 1862. A cabeça de tambor elíptica foi estudada por Émile Mathieu, levando à equação diferencial de Mathieu.

Se as bordas de uma forma são segmentos de linha reta, então uma solução é integrável ou conhecível em forma fechada apenas se for expressável como uma combinação linear finita de ondas planas que satisfazem as condições de fronteira (zero na fronteira, ou seja, membrana fixada).

Se o domínio é um círculo de raio Predefinição:Mvar, então é apropriado introduzir coordenadas polares Predefinição:Mvar e Predefinição:Mvar. A equação de Helmholtz toma a forma $$A_{rr}+\frac{1}{r}{A}_r+\frac{1}{r^2}{A}_{\theta \theta }+k^{2}A=0.$$

Podemos impor a condição de fronteira de que Predefinição:Mvar desaparece se Predefinição:Math; assim $$A(a,\theta )=0.$$

A técnica de separação de variáveis leva a soluções de teste da forma $$A(r,\theta )=R(r)\mathrm{\Theta }(\theta ),$$ onde Predefinição:Math deve ser periódica de período Predefinição:Math. Isto leva a

$${\mathrm{\Theta }}^{″}+n^2\mathrm{\Theta }=0,$$ $$r^2{R}^{″}+r{R}^{\prime }+r^2{k}^{2}R-n^{2}R=0.$$

Segue-se da condição de periodicidade que $$\mathrm{\Theta }=\alpha \cos n\theta +\beta \sin n\theta ,$$ e que Predefinição:Mvar deve ser um inteiro. A componente radial Predefinição:Mvar tem a forma $$R(r)=\gamma J_n(\rho ),$$ onde a função de Bessel Predefinição:Math satisfaz a equação de Bessel $${\rho }^2{J_n}^{″}+\rho {J_n}^{\prime }+({\rho }^2-n^2)J_n=0,$$ e Predefinição:Math. A função radial Predefinição:Math tem infinitas raízes para cada valor de Predefinição:Mvar, denotadas por Predefinição:Math. A condição de fronteira de que Predefinição:Mvar desaparece onde Predefinição:Math será satisfeita se os números de onda correspondentes forem dados por $$k_{m,n}=\frac{1}{a}{\rho }_{m,n}.$$

A solução geral Predefinição:Mvar então assume a forma de uma série de Fourier generalizada de termos envolvendo produtos de Predefinição:Math e o seno (ou cosseno) de Predefinição:Math. Estas soluções são os modos de vibração de uma cabeça de tambor circular.

Soluções tridimensionais

Predefinição:Sem-fontes Em coordenadas esféricas, a solução é:

$$A(r,\theta ,\phi )=\sum _{\ell =0}^{\mathrm{\infty }}\sum _{m=-\ell }^{\ell }(a_{\ell m}{j}_{\ell }(kr)+b_{\ell m}{y}_{\ell }(kr))Y_{\ell }^m(\theta ,\phi ).$$

Esta solução surge da solução espacial da equação de onda e da equação de difusão. Aqui Predefinição:Math e Predefinição:Math são as funções de Bessel esféricas, e Predefinição:Math são os harmónicos esféricos (Abramowitz e Stegun, 1964). Note que estas formas são soluções gerais e requerem condições de fronteira para serem usadas em qualquer caso específico. Para domínios exteriores infinitos, uma condição de radiação pode também ser necessária (Sommerfeld, 1949).

Escrevendo Predefinição:Math a função Predefinição:Math tem assimptótica

$$A(r_0)=\frac{e^{ik{r}_0}}{r_0}f(\frac{{𝐫}_0}{r_0},k,u_0)+o\left(\frac{1}{r_0}\right)\text{ à medida que }{r}_0\to \mathrm{\infty }$$

onde a função Predefinição:Mvar é chamada amplitude de dispersão e Predefinição:Math é o valor de Predefinição:Mvar em cada ponto da fronteira Predefinição:Math

Soluções tridimensionais dada a função num plano bidimensional

Dado um plano bidimensional onde A é conhecido, a solução para a equação de Helmholtz é dada por: $$A(x,y,z)=-\frac{1}{2\pi }\underset{-\mathrm{\infty }}{\overset{+\mathrm{\infty }}{\iint }}{A}^{\prime }(x^{\prime },y^{\prime })\frac{e^{ikr}}{r}\frac{z}{r}(ik-\frac{1}{r}),d{x}^{\prime }d{y}^{\prime },$$

onde

$A^{\prime }(x^{\prime },y^{\prime })$ é a solução no plano bidimensional, $r=\sqrt{(x-x^{\prime }{)}^2+(y-y^{\prime }{)}^2+z^2},$ À medida que z se aproxima de zero, todas as contribuições da integral desaparecem exceto para r=0. Assim $A(x,y,0)=A^{\prime }(x,y)$ até um fator numérico, que pode ser verificado como 1 transformando a integral em coordenadas polares $(\rho ,\theta )$.

Esta solução é importante na teoria da difração, por exemplo, na derivação da difração de Fresnel.

Aproximação paraxial

Predefinição:Mais informações Na aproximação paraxial da equação de Helmholtz, a amplitude complexa Predefinição:Mvar é expressa como $$A(𝐫)=u(𝐫)e^{ikz}$$ onde Predefinição:Mvar representa a amplitude complexa que modula a onda plana sinusoidal representada pelo fator exponencial. Então, sob uma suposição adequada, Predefinição:Mvar aproximadamente resolve $$\underset{\perp }{\overset{2}{\mathrm{\nabla }}}u+2ik\frac{\partial u}{\partial z}=0,$$ onde ${\mathrm{\nabla }}_{\perp }^2\stackrel{\text{ def }}{=}\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}$ é a parte transversal do Laplaciano.

Esta equação tem aplicações importantes na ciência da ótica, onde fornece soluções que descrevem a propagação de ondas eletromagnéticas (luz) na forma de ondas parabólicas ou feixes gaussianos. A maioria dos lasers emite feixes que assumem esta forma.

A suposição sob a qual a aproximação paraxial é válida é que a derivada em Predefinição:Mvar da função de amplitude Predefinição:Mvar é uma função lentamente variável de Predefinição:Mvar:

$$\left|\frac{{\partial }^{2}u}{\partial z^2}\right|\ll \left|k\frac{\partial u}{\partial z}\right|.$$

Esta condição é equivalente a dizer que o ângulo Predefinição:Mvar entre o vetor de onda Predefinição:Math e o eixo óptico Predefinição:Mvar é pequeno: Predefinição:Math.

A forma paraxial da equação de Helmholtz é encontrada substituindo a expressão acima para a amplitude complexa na forma geral da equação de Helmholtz como segue:

$${\mathrm{\nabla }}^2(u(x,y,z)e^{ikz})+k^{2}u(x,y,z)e^{ikz}=0.$$

A expansão e o cancelamento produzem o seguinte:

$$(\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2})u(x,y,z)e^{ikz}+(\frac{{\partial }^2}{\partial z^2}u(x,y,z))e^{ikz}+2(\frac{\partial }{\partial z}u(x,y,z))ik{e}^{ikz}=0.$$

Devido à desigualdade paraxial mencionada acima, o termo Predefinição:Math é negligenciado em comparação com o termo Predefinição:Math. Isto produz a equação de Helmholtz paraxial. Substituindo Predefinição:Math obtém-se então a equação paraxial para a amplitude complexa original Predefinição:Mvar:

$$\underset{\perp }{\overset{2}{\mathrm{\nabla }}}A+2ik\frac{\partial A}{\partial z}=0.$$

A integral de difração de Fresnel é uma solução exata para a equação de Helmholtz paraxial.

Equação de Helmholtz não homogénea

Predefinição:Multiple image A equação de Helmholtz não homogénea é a equação $${\mathrm{\nabla }}^{2}A(𝐱)+k^{2}A(𝐱)=-f(𝐱)\text{ em }{ℝ}^n,$$ onde Predefinição:Math é uma função com suporte compacto, e Predefinição:Math Esta equação é muito semelhante à equação de Poisson com triagem, e seria idêntica se o sinal de mais (na frente do termo Predefinição:Mvar) fosse trocado por um sinal de menos.

Para resolver esta equação de forma única, é necessário especificar uma condição de fronteira no infinito, que é tipicamente a condição de radiação de Sommerfeld

$$\underset{r\to \mathrm{\infty }}{\lim }{r}^{\frac{n-1}{2}}(\frac{\partial }{\partial r}-ik)A(𝐱)=0$$

em $n$ dimensões espaciais, para todos os ângulos (ou seja, qualquer valor de $\theta ,\varphi$). Aqui $r=\sqrt{\sum _{i=1}^n{x}_i^2}$ onde $x_i$ são as coordenadas do vetor $𝐱$.

Com esta condição, a solução para a equação de Helmholtz não homogénea é

$$A(𝐱)=\underset{{ℝ}^n}{\int }!G(𝐱,{𝐱}^{\prime })f({𝐱}^{\prime }),\mathrm{d}{𝐱}^{\prime }$$

(repare que esta integral é na verdade sobre uma região finita, já que Predefinição:Mvar tem suporte compacto). Aqui, Predefinição:Mvar é a função de Green desta equação, ou seja, a solução para a equação de Helmholtz não homogénea com Predefinição:Math igual à função delta de Dirac, então Predefinição:Mvar satisfaz

$${\mathrm{\nabla }}^{2}G(𝐱,{𝐱}^{\prime })+k^{2}G(𝐱,{𝐱}^{\prime })=-\delta (𝐱,{𝐱}^{\prime })\in {ℝ}^n.$$

A expressão para a função de Green depende da dimensão Predefinição:Mvar do espaço. Temos $$G(x,x^{\prime })=\frac{i{e}^{ik|x-x^{\prime }|}}{2k}$$ para Predefinição:Math,

$$G(𝐱,{𝐱}^{\prime })=\frac{i}{4}{H}_0^{(1)}(k|𝐱-{𝐱}^{\prime }|)$$ para Predefinição:Math, onde Predefinição:Math é uma função de Hankel, e $$G(𝐱,{𝐱}^{\prime })=\frac{e^{ik|𝐱-{𝐱}^{\prime }|}}{4\pi |𝐱-{𝐱}^{\prime }|}$$ para Predefinição:Math. Note que escolhemos a condição de fronteira de que a função de Green é uma onda que se afasta para Predefinição:Math.

Finalmente, para um n geral,

$$G(𝐱,{𝐱}^{\prime })=c_d{k}^p\frac{H_p^{(1)}(k|𝐱-{𝐱}^{\prime }|)}{|𝐱-{𝐱}^{\prime }{|}^p}$$

onde $p=\frac{n-2}{2}$ e $c_d=\frac{1}{2i(2\pi {)}^p}$.