Grupo3 - Ondas2

Introdução

Equações diferenciais parciais (EDP's) hiperbólicas geralmente podem ser formuladas a partir de teoremas de conservação. Um exemplo é a equação do tipo:

$\frac{\partial \overrightarrow{U}}{\partial t}+\mathrm{\nabla }\cdot \overrightarrow{F}(U)=\overrightarrow{S}(U)$,

onde $\overrightarrow{U}$ é o vetor de densidades da quantidade conservada, i.e., $\overrightarrow{U}=(U_1,...,U_n)$, $\overrightarrow{F}$ é o fluxo de densidade e $\overrightarrow{S}$ é um termo genérico representando fontes ou sumidouros.

Uma classe especial de equações hiperbólicas são as chamadas equações de adveção, na qual a derivada temporal da quantidade conservada $\overrightarrow{U}(\overrightarrow{x},t)$ é proporcional à sua derivada espacial. Nesses casos, $\overrightarrow{F}(U)$ é diagonal e dada por:

$\overrightarrow{F}(U)=v\overrightarrow{I}\cdot \overrightarrow{U}$,

onde $\overrightarrow{I}$ é a matriz identidade.

Considerando apenas uma dimensão e com $\overrightarrow{U}\equiv u$, temos a equação de adveção:

$\frac{\partial u}{\partial t}+v\frac{\partial u}{\partial x}=0$,

onde $v$ é a velocidade de propagação do pulso gerado. A equação admite uma solução analítica da forma $u=f(x-vt)$, representando uma onda se movendo na direção $x$.

A equação da onda em uma dimensão é uma EDP hiperbólica de segunda ordem dada por

$\frac{{\partial }^{2}u}{\partial t^2}=v^2\frac{{\partial }^{2}u}{\partial x^2}.$

E admite duas soluções, representadas por pulsos, $f(x+vt)$ e $f(x-vt)$.

Assumindo que $v\ne v(x)$ na equação da onda, nos restringimos a problemas lineares. Além disso, se escrevermos

$k=v\frac{\partial u}{\partial x}$, $s=\frac{\partial u}{\partial t},$

então a equação da onda pode ser escrita como um sistema de três equações diferenciais de primeira ordem:

$\{\begin{array}{l}\frac{\partial k}{\partial t}=v\frac{\partial s}{\partial x}\\ \frac{\partial s}{\partial t}=v\frac{\partial k}{\partial x}\\ \frac{\partial u}{\partial t}=s\end{array}$

Em notação vetorial, o sistema acima pode ser reescrito na forma conservativa como: $\frac{\partial U}{\partial t}+\frac{\partial F(U)}{\partial x}=0$,

onde $U=\left(\begin{array}{c}k\\ s\end{array}\right),\text{e}F(U)=\left(\begin{array}{cc}0& -v\\ -v& 0\end{array}\right)$

O Problema Físico

O Modelo de Corda Ideal

Para uma primeira abordagem da equação da onda, podemos primeiro dividir o comprimento $L$ da corda em $K$ intervalos de comprimentos iguais, dessa forma $\mathrm{\Delta }x=\frac{L}{K}$. Cada intervalo é discretizado, portanto, como $x_i$, $i=0,1,...,K$. Também podemos dividir o tempo em intervalos iguais $\mathrm{\Delta }t$ e denotá-los como $t_n$, $n=0,1...$.

Tendo feita a discretização das variáveis, podemos aproximar a equação da onda por diferenciação finita:

$\frac{u(i,n+1)-2u(i,n)+u(i,n-1)}{(\mathrm{\Delta }t{)}^2}=v^2\frac{u(i+1,n)-2u(i,n)+u(i-1,n)}{(\mathrm{\Delta }x{)}^2}$.

Sabendo que essa discretização da equação da onda pode ser verificada como sendo o método Leapfrog (ver seção do método de Leapfrog), podemos resolver a equação para $u(i,n+1)$ para sabermos o deslocamento de uma partição da corda no momento de tempo seguinte, assim obtendo

$u(i,n+1)=2(1-r^2)u(i,n)+r^2[u(i+1,n)+u(i-1,n)]-u(i,n-1)$,

onde $r=v\frac{\mathrm{\Delta }t}{\mathrm{\Delta }x}.$

Um Quadro Mais Realístico - O Modelo de Corda Rígida

Para nos aproximarmos de um modelo mais real, podemos adicionar um termo à equação original da onda que corresponde ao efeito de fricção em uma corda. De acordo com [1], a equação da onda mais geral com efeito de fricção pode ser escrita como

$\frac{{\partial }^{2}y}{\partial t^2}=v^2(\frac{{\partial }^{2}y}{\partial x^2}-ϵL^2\frac{{\partial }^{4}y}{\partial x^4}),$

onde $v$ é a velocidade transversal de propagação do pulso na corda, dada pela relação $v=\sqrt[]{\frac{T}{\rho }}$ (sendo $T$ a tensão na corda e $\rho$ a densidade linear da mesma), $ϵ$ é um parâmetro adimensional de fricção que representa a rigidez da corda e $L$ o comprimento da corda.

O parâmetro $ϵ$ é dado por

$$,

onde $\kappa$ é o raio da corda, $E$ é o Módulo de Young e $S$ a área da secção da corda.

Ao discretizarmos a equação da onda em uma corda com fricção e a resolvendo para $u_i^{n+1}$ obtemos:

$u_i^{n+1}=[2-2r^2-6ϵr^2{K}^2]u_i^n-u_i^{n-1}+r^2[1+4ϵK^2][u_{i+1}^n+u_{i-1}^n]-ϵr^2{K}^2[u_{i+2}^n+u_{i-2}^n].$

O fato de essa discretização depender do deslocamento da corda em posições $i-2$ e $i+2$ implica em precisarmos simular "pontos fantasmas" quando integramos os extremos das cordas. Para fazermos isso, podemos ou utilizar a aproximação $u_{-1}^n=-u_{+1}^n$ ou podemos considerar esses "pontos fantasmas" como pontos presos e, portanto, sempre iguais a zero.

Os Métodos Utilizados

Foi realizada uma abordagem ao problema da corda real a partir de três métodos diferentes de integração numérica. Os três são métodos para fins de resolução de equações diferenciais parciais da forma apresentada anteriormente.

O método mais básico é chamado de FTCS (Forward-Time-Centered-Space) e consiste em duas expansões de Taylor ao redor do ponto $x_j$:

$$

$$

Subtraindo as duas expressões, encontramos a expressão

$$,

A qual podemos substituir na equação da onda, juntamente com a discretização da derivada parcial temporal. Temos então que, para um sistema linear de equações hiperbólicas:

${\text{U}}_j^{n+1}={\text{U}}_j^n-\frac{\mathrm{\Delta }t}{2\mathrm{\Delta }x}[{\text{F}}_{j+1}^n-{\text{F}}_{j-1}^n]+{𝒪}(\mathrm{\Delta }{t}^2,\mathrm{\Delta }{x}^2\mathrm{\Delta }t)$

Visto que essa última notação é mais genérica, ela será utilizada para a explicação dos métodos posteriores.

O Método de Lax-Friedrichs

O método de Lax-Friedrichs consiste em substituir o termo ${\text{U}}_j^n$ com sua respectiva média espacial, i.e., $u_j^n=(u_{j+1}^n+u_{j-1}^n)/2$. Logo, temos a seguinte equação de recorrência:

${\text{U}}_j^{n+1}=\frac{1}{2}({\text{U}}_{j+1}^n+{\text{U}}_{j-1}^n)-\frac{\mathrm{\Delta }t}{2\mathrm{\Delta }x}[{\text{F}}_{j+1}^n-{\text{F}}_{j-1}^n]+{𝒪}(\mathrm{\Delta }{x}^2)$

O Método de Lax-Wendroff

O método de Lax-Wendroff é a extensão do método de Lax-Friedrichs de segunda ordem. Calculamos o vetor $U$ a partir de um passo médio de Lax-Friedrichs:

${\text{U}}_{j+\frac{1}{2}}^{n+\frac{1}{2}}=\frac{1}{2}[{\text{U}}_{j+1}^n+{\text{U}}_j^n]-\frac{\mathrm{\Delta }t}{2\mathrm{\Delta }x}[{\text{F}}_{j+1}^n-{\text{F}}_j^n]+{𝒪}(\mathrm{\Delta }{x}^2)$,

${\text{U}}_{j-\frac{1}{2}}^{n+\frac{1}{2}}=\frac{1}{2}[{\text{U}}_j^n+{\text{U}}_{j-1}^n]-\frac{\mathrm{\Delta }t}{2\mathrm{\Delta }x}[{\text{F}}_j^n-{\text{F}}_{j-1}^n]+{𝒪}(\mathrm{\Delta }{x}^2)$,

E encontramos os fluxos ${\text{F}}_{j\pm \frac{1}{2}}^{n+\frac{1}{2}}$ a partir dos valores de ${\text{U}}_{j\pm \frac{1}{2}}^{n+\frac{1}{2}}.$

Logo, com um meio passo de Leapfrog, temos a expressão final do método:

${\text{U}}_j^{n+1}={\text{U}}_j^n-\frac{\mathrm{\Delta }t}{\mathrm{\Delta }x}[{\text{F}}_{j+\frac{1}{2}}^{n+\frac{1}{2}}-{\text{F}}_{j-\frac{1}{2}}^{n+\frac{1}{2}}]+{𝒪}(\mathrm{\Delta }{x}^2)$