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 equações de conservação :

$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$

Onde $A_{ij}$ e $a_{ij}$ são os elementos das matrizes $\mathbf {A}$ e $\mathbf {a}$e $l_{i}$ e $b_{i}$ são elementos de vetores. Será perguntando se é possível reescrever essa equação para

$m_{j}\left(\beta {\frac {\partial u_{j}}{\partial t}}+\alpha {\frac {\partial u_{j}}{\partial x}}\right)+l_{j}b_{j}=0$

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

${\frac {du_{j}}{d\eta }}=T'{\frac {\partial u_{j}}{\partial t}}+X'{\frac {\partial u_{j}}{\partial x}}$

comparando as duas últimas equações, encontramos

$\alpha =X'(\eta ),\beta =T'(\eta )$

que agora pode ser escrito em forma característica

$m_{j}{\frac {du_{j}}{d\eta }}+l_{j}b_{j}=0$

onde devemos ter as condições

$l_{i}A_{ij}=m_{j}T'$

$l_{i}a_{ij}=m_{j}X'$

onde $m_{j}$ pode ser eliminado para dar a condição necessária

$l_{i}(A_{ij}X'-a_{ij}T')=0$

portanto, para uma solução não trivial, o determinante é igual a zero

$|A_{ij}X'-a_{ij}T'|=0$

Para os invariantes de Riemann nós estamos preocupados com o caso quando a matriz $\mathbf {A}$ é uma matriz identidade da forma

${\frac {\partial u_{j}}{\partial t}}+a_{ij}{\frac {\partial u_{j}}{\partial x}}=0$

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

$l_{i}{\frac {du_{i}}{dt}}=0$ com ${\frac {dx}{dt}}=\lambda$

Onde $l$ é o vetor próprio esquerdo da matriz $\mathbf {A}$ e $\lambda 's$ é a velocidade característica dos autovalores que a matriz $\mathbf {A}$ possui, os quais satisfazem

$|A-\lambda \delta _{ij}|=0$

Para simplificar essas equações características, nós podemos fazer a transformação ${\frac {dr}{dt}}=l_{i}{\frac {du_{i}}{dt}}$

que chega a isso

$\mu l_{i}du_{i}=dr$

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

${\frac {dr}{dt}}=0$ em ${\frac {dx}{dt}}=\lambda _{i}$

que é equivalente ao sistema diagonal

$r_{t}^{k}+\lambda _{k}r_{x}^{k}=0,$ $k=1,...,N.$

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 $\rho$ e velocidade $u$ é

$\rho _{t}+\rho u_{x}+u\rho _{x}=0$

$u_{t}+uu_{x}+(c^{2}/\rho )\rho _{x}=0$

$c$, sendo a velocidade do som, é introduzido devido à suposição isentrópica. Escreva este sistema em forma de matriz

$\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)$

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

$\lambda ^{2}-2u\lambda +u^{2}-c^{2}=0$

para obtermos

$\lambda =u\pm c$

e os autovetores são encontrados como

$\left({\begin{matrix}1\\{\frac {c}{\rho }}\end{matrix}}\right),\left({\begin{matrix}1\\-{\frac {c}{\rho }}\end{matrix}}\right)$

onde os invariantes de Riemman são

$r_{1}=J_{+}=u+\int {\frac {c}{\rho }}d\rho ,$

$r_{2}=J_{-}=u-\int {\frac {c}{\rho }}d\rho ,$

( $J_{+}$ e $J_{-}$ são as notações amplamente usadas na dinâmica dos gases ). Para o caso de um gás perfeito com aquecimentos específicos constantes, existe a relação $c^{2}={\text{const}}\,\gamma \rho ^{\gamma -1}$, Onde $\gamma$ é a razão de calor específica, para dar aos invariantes de Riemann

$J_{+}=u+{\frac {2}{\gamma -1}}c,$

$J_{-}=u-{\frac {2}{\gamma -1}}c,$

para dar as equações

${\frac {\partial J_{+}}{\partial t}}+(u+c){\frac {\partial J_{+}}{\partial x}}=0$

${\frac {\partial J_{-}}{\partial t}}+(u-c){\frac {\partial J_{-}}{\partial x}}=0$

Em outras palavras,

${\begin{aligned}&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{aligned}}$

Onde $C_{+}$ e $C_{-}$ 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

$A{\frac {\partial v}{\partial t}}+B{\frac {\partial v}{\partial x}}=0$

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

Veja também

• Onda simples