Axisymmetric reduction of the adjoint operator in ideal MHD equilibrium

Resumo Técnico: Redução Axissimétrica do Operador Adjunto em Equilíbrio de MHD Ideal

Enunciado do Problema
Cálculos de equilíbrio de magnetohidrodinâmica (MHD) ideal são fundamentais para o projeto de dispositivos de confinamento magnético e para a reconstrução de estados de plasma a partir de dados de diagnóstico. Enquanto os equilíbrios axissimétricos (ex: tokamaks) são governados pela equação de Grad–Shafranov, os equilíbrios tridimensionais (3D) (ex: stellarators) requerem resolvedores numéricos. Ambos os regimes utilizam cada vez mais abordagens baseadas em otimização, onde a análise de sensibilidade em relação a inúmeros parâmetros é necessária. O método adjunto oferece um arcabouço eficiente para tal análise de sensibilidade. No entanto, uma relação matemática precisa entre a formulação adjunta do problema de equilíbrio 3D completo e o adjunto da equação de Grad–Shafranov axissimétrica não foi explicitamente demonstrada. Especificamente, permanece incerto como o operador adjunto 3D se reduz sob a suposição de axissimetria e se a estrutura resultante se alinha com as formulações adjuntas axissimétricas estabelecidas, particularmente no que diz respeito aos fatores de ponderação geométrica inerentes à energia magnética axissimétrica.

Metodologia
Os autores empregam uma abordagem de teoria de operadores para derivar o adjunto das equações de equilíbrio de MHD ideal sem assumir simetria inicialmente.

1. Formulações Adjuntas Gerais: O equilíbrio de MHD é tratado como um problema de otimização restrito onde as variáveis de estado são a densidade de corrente ($J$), o campo magnético ($B$) e a pressão ($p$). Uma Lagrangiana é construída utilizando as restrições de equilíbrio e um funcional objetivo.
2. Derivação Adjunta 3D: Ao tomar as variações da Lagrangiana em relação às variáveis de estado e aplicar o teorema da divergência sob condições de contorno fixas ($B \cdot n = 0$), os autores derivam as equações adjuntas 3D. As variáveis adjuntas ($\alpha, \beta, \gamma$) são introduzidas como multiplicadores de Lagrange duais às restrições de equilíbrio estático, distintas de deslocamentos lagrangianos físicos.
3. Redução Axissimétrica: Os autores impõem axissimetria ao sistema adjunto 3D derivado. Eles utilizam representações padrão para campos vetoriais de divergência nula em geometria toroidal, expressando o campo magnético e o campo adjunto $\alpha$ em termos de funções de fluxo poloidal ($\psi$) e ângulos toroidais ($\phi$).
4. Projeção e Simplificação: O sistema reduzido é projetado ao longo e através da direção toroidal para isolar componentes escalares. Os autores analisam especificamente a quantidade escalar $\eta = -\alpha \cdot \nabla \psi$, interpretando-a como uma perturbação do tipo fluxo.
5. Comparação de Operadores: O operador escalar resultante é comparado com o operador de Grad–Shafranov perturbado. Os autores também examinam a estrutura do produto interno, especificamente o papel do peso geométrico $1/R^2$ decorrente da variação da energia magnética poloidal.

Contribuições Principais e Resultados

• Redução Explícita: O artigo demonstra que a redução axissimétrica do sistema adjunto homogêneo 3D completo produz um operador escalar que coincide exatamente com o operador de Grad–Shafranov perturbado (Zakharov & Pletzer 1999). Este operador governa variações infinitesimais da função de fluxo de equilíbrio.
• Correspondência de Variáveis: O estudo estabelece uma correspondência direta entre as variáveis adjuntas 3D e a variável escalar reduzida. A variável adjunta escalar reduzida é identificada como $\eta = -\alpha \cdot \nabla \psi$, onde $\alpha$ é o campo adjunto incompressível associado à equação de balanço de forças.
• Estrutura de Produto Interno Ponderado: A análise esclarece que o operador escalar reduzido é formalmente autoadjunto apenas em relação a um produto interno ponderado $(u, v)_w = \int u v R^{-2} dV$. Este fator de peso $1/R^2$ mostra-se originado naturalmente da expressão axissimétrica da energia magnética poloidal, em vez de ser uma imposição arbitrária.
• Consistência Estrutural: O trabalho resolve discrepâncias com formulações anteriores (ex: Sanpei et al. 2021, 2023) que omitiram o peso geométrico no produto interno. Os autores mostram que, embora formulações não ponderadas sejam praticamente viáveis para otimização numérica, a formulação ponderada é necessária para manter a compatibilidade estrutural e a autoadjunção entre o adjunto 3D e sua redução axissimétrica.

Significância
O artigo fornece uma declaração de consistência rigorosa em nível de operador, vinculando a formulação adjunta de MHD ideal 3D completa ao problema adjunto de Grad–Shafranov axissimétrico. Ao esclarecer como a estrutura adjunta axissimétrica está incorporada no arcabouço 3D, o trabalho oferece uma base teórica unificada para análise de sensibilidade. Isto é particularmente relevante para:

• Configurações Fracamente Não-Axissimétricas: A formulação adjunta 3D serve como um "arcabouço pai" para analisar equilíbrios que são pequenas perturbações de um estado axissimétrico, permitindo a organização sistemática das respostas axissimétricas e correções 3D.
• Reconstrução e Otimização Baseadas em Adjuntos: A correspondência explícita entre variáveis 3D e axissimétricas auxilia na interpretação de resultados de sensibilidade e pode auxiliar na análise de bifurcação ou na caracterização de respostas singulares em sistemas de equilíbrio linearizados.

Os autores concluem que esta redução esclarece a relação matemática entre os problemas adjuntos de MHD ideal 3D e axissimétricos, garantindo que as propriedades estruturais (como a autoadjunção) sejam preservadas entre as dimensões quando os pesos geométicos apropriados são aplicados.