Mostrar código-fonte de Invariantes de Riemann

'''Os invariantes de Riemann''' são transformações matemáticas feitas em um sistema de equações conservativas para as tornarem mais fáceis de serem resolvidas. As invariantes de Riemann são constantes ao longo de curvas características de equações diferenciais parciais onde obtêm o nome de invariantes. Elas foram obtidas pela primeira vez por [[Bernhard Riemann]] em seu trabalho com ondas planas nas dinâmicas dos gases. 

== Teoria matemática ==
Considere o conjunto de [[Lei de conservação|equações]] de [[Lei de conservação|conservação]] : 

: <math> l_i\left(A_{ij} \frac{\partial u_j}{\partial t} +a_{ij}\frac{\partial u_j}{\partial x} \right)+l_j b_j=0 </math>

Onde <math>A_{ij} </math> e <math>a_{ij} </math> são os [[Elemento (matemática)|elementos]] das [[Matriz (matemática)|matrizes]] <math>\mathbf{A}</math> e <math>\mathbf{a}</math>e <math>l_{i} </math> e <math>b_{i} </math> são elementos de vetores. Será perguntando se é possível reescrever essa equação para 

: <math> m_j\left(\beta\frac{\partial u_j}{\partial t} +\alpha\frac{\partial u_j}{\partial x} \right)+l_j b_j=0 </math>

Para fazer isso essas curvas serão introduzidas no plano (x,t) definido pelo campo vectorial <math>(\alpha,\beta)</math>. O termo em parênteses será reescrito em termos da derivada total x,t são parametrizados como<math>x=X(\eta),t=T(\eta)</math> 

: <math> \frac{d u_j}{d \eta}=T'\frac{\partial u_j}{\partial t}+X'\frac{\partial u_j}{\partial x}</math>

comparando as duas últimas equações, encontramos 

: <math> \alpha=X'(\eta), \beta=T'(\eta) </math>

que agora pode ser escrito em forma característica 

: <math> m_j\frac{du_j }{ d \eta }+l_jb_j = 0</math>

onde devemos ter as condições 

: <math>l_iA_{ij}=m_jT'</math>
: <math>l_ia_{ij}=m_jX'</math>

onde <math>m_j</math> pode ser eliminado para dar a condição necessária 

: <math>l_i(A_{ij}X'-a_{ij}T')=0</math>

portanto, para uma [[Trivialidade (matemática)|solução não trivial]], o determinante é igual a zero 

: <math>|A_{ij}X'-a_{ij}T'|=0</math>

Para os invariantes de Riemann nós estamos preocupados com o caso quando a matriz <math>\mathbf{A}</math> é uma [[matriz identidade]] da forma 

: <math> \frac{\partial u_j}{\partial t} +a_{ij}\frac{\partial u_j}{\partial x}=0 </math>

observe que isso é homogêneo devido ao vetor <math>\mathbf{n}</math> ser zero. Na forma característica o sistema é 

: <math> l_i\frac{du_i }{dt }=0</math> com <math>\frac{dx }{dt }=\lambda </math> 

Onde <math>l</math> é o [[Autovalores e autovetores|vetor próprio]] esquerdo da matriz <math>\mathbf{A}</math> e <math>\lambda 's</math> é a velocidade característica dos [[Autovalores e autovetores|autovalores]] que a matriz <math>\mathbf{A}</math> possui, os quais satisfazem 

: <math> |A -\lambda\delta_{ij}|=0 </math>

Para simplificar essas equações características, nós podemos fazer a transformação <math> \frac{dr}{dt}=l_i\frac{du_i}{dt}</math> 

que chega a isso 

: <math> \mu l_idu_i =dr </math>

Um fator integrante <math>\mu</math> pode ser multiplicado para ajudar a integrar nesse caso. Assim o sistema agora esta na forma característica 

: <math> \frac{dr}{dt }=0</math> em <math>\frac{dx}{dt}=\lambda_i</math> 

que é equivalente ao [[Matriz diagonal|sistema diagonal]] <ref name="Não_nomeado-xwCK-1"></ref> 

: <math>r_t^k +\lambda_kr_x^k=0, </math> <math>k=1,...,N.</math>

A solução deste sistema pode ser dado pelo método holográfico generalizado 

== Exemplo ==
Considere que a equação de Euler unidimensional escrita em termos de densidade <math>\rho</math> e velocidade <math>u</math> é  

: <math> \rho_t+\rho u_x+u\rho_x=0 </math>
: <math> u_t+uu_x+(c^2/\rho)\rho_x=0 </math>

<math>c</math>, sendo a [[velocidade do som]], é introduzido devido à suposição isentrópica. Escreva este sistema em forma de matriz 

: <math> \left( \begin{matrix} \rho\\ u \end{matrix}\right)_t +\left( \begin{matrix} u&\rho\\ \frac{c^2 }{\rho }&u \end{matrix}\right) \left( \begin{matrix} \rho\\ u \end{matrix}\right)_x=\left( \begin{matrix} 0\\ 0 \end{matrix}\right) </math>

a partir da análise acima, os valores próprios e os vetores próprios precisam ser encontrados. Os valores próprios são encontrados para satisfazer 

: <math> \lambda^2-2u\lambda+u^2-c^2=0 </math>

para obtermos 

: <math> \lambda=u\pm c </math>

e os autovetores são encontrados como 

: <math> \left( \begin{matrix} 1\\ \frac{c }{\rho } \end{matrix}\right),\left( \begin{matrix} 1\\ -\frac{c }{\rho } \end{matrix}\right) </math>

onde os invariantes de Riemman são 

: <math> r_1=J_+=u+\int \frac{c}{\rho}d\rho, </math>
: <math> r_2=J_-=u-\int \frac{c}{\rho}d\rho, </math>

( <math>J_+</math> e <math>J_-</math> são as notações amplamente usadas na [[Fluxo compressível|dinâmica dos gases]] ). Para o caso de um gás perfeito com aquecimentos específicos constantes, existe a relação <math>c^2=\text{const}\, \gamma \rho^{\gamma-1}</math>, Onde <math>\gamma</math> é a [[Coeficiente de expansão adiabática|razão de calor específica]], para dar aos invariantes de Riemann <ref name="Não_nomeado-xwCK-1"/> <ref name="Não_nomeado-xwCK-1"/> 

: <math> J_+=u+\frac{2}{\gamma-1}c, </math>
: <math> J_-=u-\frac{2}{\gamma-1}c, </math>

para dar as equações 

: <math> \frac{\partial J_+}{\partial t}+(u+c)\frac{\partial J_+}{\partial x}=0</math>
: <math> \frac{\partial J_-}{\partial t}+(u-c)\frac{\partial J_-}{\partial x}=0</math>

Em outras palavras, 

: <math>
\begin{align}
&dJ_+ = 0, \, J_+=\text{const}\quad \text{along}\,\, C_+\, :\, \frac{dx}{dt}=u+c, \\
&dJ_- = 0, \, J_-=\text{const}\quad \text{along}\,\, C_-\, :\, \frac{dx}{dt}=u-c,
\end{align}</math>

Onde <math>C_+</math> e <math>C_-</math> são curvas características. Isso pode ser resolvido pela transformação holográfica. No plano holográfico, se todas as características colapsam em uma única curva, então nós obtemos uma única onda. Se a forma da matriz do sistema de pde´s está na forma 

: <math> A\frac{\partial v}{\partial t}+B\frac{\partial v}{\partial x}=0 </math>

Então, pode ser possível multiplicar através pela matriz inversa <math>A^{-1} </math> desde que a matriz determinante <math>\mathbf{A}</math> de não seja zero. 

== Veja também ==

* Onda simples 

== Referências ==
{{Referências}}
[[Categoria:Bernhard Riemann]]
[[Categoria:Equações de derivadas parciais]]
[[Categoria:Teoria de invariantes]]