Equação de Ginzburg-Landau complexa

De Física Computacional
Revisão de 20h57min de 27 de abril de 2024 por Joaovgm (discussão | contribs) (→‎Dedução)
Ir para navegação Ir para pesquisar

A equação de Ginzburg-Landau complexa (CGLE) surgiu inicialmente em 1969 como um modelo para o inicio de instabilidades em problemas de convecção de fluídos. A partir de então, ela se tornou uma das equações não lineares mais estudadas da física, descrevendo uma variedade enorme de fenômenos como:

  • Ondas não lineares;
  • Transições de fase de segunda ordem;
  • Supercondutividade;
  • Superfluidez;
  • Condensado de Bose-Einstein.

A equação de Ginzburg-Landau complexa, quando escrita de modo a minimizar o número de constantes, é dada pela equação abaixo:

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial A}{\partial t} = (1+ic_1)\nabla^2 A + A - (1-ic_3) A|A|^2. }

É possível deduzir a CGLE a partir do oscilador linear harmônico por meio de argumentos de simetria, encontrando a equação de Stuart-Landau, e, em seguida, considerando um sistema estendido no espaço.

Dedução

Espaço de fase do oscilador harmônico

A energia de um oscilador harmônico é expressa pela equação abaixo, onde Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle E} é a energia, Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle q} e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle p} a coordenada e seu respectivo momento, Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle m} é a massa e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega_0} a frequência angular

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle E = \frac{p^2}{2m} + \frac{1}{2} m \omega_0^2 q^2. }

Ao realizar as seguintes mudanças de variáveis, Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle q \rightarrow q/m^{1/2}} e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle p \rightarrow p m^{1/2}} , a equação da energia produz trajetórias circulares no espaço de fase de e

Essa é uma importante simetria do oscilador harmônico linear, resultando que a sua energia é proporcional ao quadrado da amplitude de oscilação, não dependendo da fase. Isso sugere uma motivação, qual é o menor termo não linear que pode ser adicionado de modo a preservar essa simetria. Para tanto, o estado do sistema será descrito em coordenadas polares, onde é a amplitude e a fase

Define-se, então, a variável complexa , portanto a equação acima pode ser reescrita como

Ao realizar a transformação de variável , com , a equação acima permanece inalterada. Ou seja, a equação é invariante a rotações. Então, busca-se uma função não linear tal que

também seja invariante a rotações.

Então, perante às transformações e , a função deve satisfazer

para que seja possível fatorar o termo responsável pela rotação e obter novamente a equação original.

Considerando pequenas oscilações, é possível expandir Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(A, A^*)} em potências de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle A} e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle A^*} até a menor ordem possível que satisfaça a condição e que introduza uma não linearidade à equação. Com isso, obtém-se

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(A, A^*) = \alpha_1 A + \alpha_2 |A|^2 A, \quad \alpha_1 = \alpha_{1r} + i \alpha_{1i}, \quad \alpha_2 = \alpha_{2r} + i \alpha_{2i} }

Utilizando o resultado encontrado e expressando em coordenadas polares por meio de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle A = R e^{i\phi}}

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \dot{R} = \alpha_{1r} R + \alpha_{2r} R^3, \quad \dot{\phi} = \omega_0 + \alpha_{1i} + \alpha_{2i} R^2. }

Em seguida, muda-se para o referencial que gira com a mesma frequência do oscilador harmônico por meio da definição de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \phi = \varphi + \omega_0 t} . As novas equações obtidas são

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \dot{R} = \alpha_{1r} R + \alpha_{2r} R^3, \quad \dot{\varphi} = \alpha_{1i} + \alpha_{2i} R^2. }


Método FTCS