VEQ: a fast parametric Grad--Shafranov solver for fixed-boundary tokamak equilibria with flexible source profiles

O artigo introduz o VEQ, um resolvedor paramétrico rápido para equilíbrios de tokamak de fronteira fixa que utiliza perfis de fonte flexíveis e expansões harmônicas para alcançar reconstrução de equilíbrio de baixa latência e alta precisão, adequado para consultas repetidas em fluxos de trabalho de acoplamento de geometria de transporte.

O essencial

Imagine um tokamak (um reator de fusão nuclear em forma de rosquinha) como um balão gigante e invisível cheio de plasma superquente. Para evitar que esse balão estoure ou colapse, os cientistas precisam saber exatamente como a "pele" magnética que o mantém unido é moldada e como as forças dentro dele estão equilibradas. Isso é chamado de encontrar o "equilíbrio".

Normalmente, calcular esse equilíbrio é como tentar resolver um quebra-cabeça 3D massivo onde cada peça muda de forma conforme você a toca. É preciso, mas leva muito tempo — tempo demais se você precisar verificar a forma centenas de vezes por segundo para controlar o reator.

Este artigo apresenta o VEQ (Veloce EQuilibrium), uma nova ferramenta projetada para ser o botão de "avançar rápido" para esses cálculos. Veja como funciona, usando analogias simples:

1. O "Transformista" vs. O "Pintor de Pixels"

Os métodos tradicionais são como um Pintor de Pixels. Eles tentam mapear todo o campo magnético calculando o valor em milhões de minúsculos pontos de grade (pixels). É detalhado, mas pesado e lento.

O VEQ é mais como um Transformista. Em vez de pintar cada pixel, ele descreve a forma do plasma usando um "esqueleto" matemático flexível feito de alguns botões e seletores principais (chamados de parâmetros).

• Pense na forma do plasma como um pedaço de argila.
• Os resolvedores tradicionais tentam esculpir cada pequena saliência na argila.
• O VEQ usa um conjunto de ferramentas pré-definidas de "esticar" e "dobrar" (harmônicos matemáticos e polinômios). Você apenas gira alguns seletores para esticar a argila até que ela ganhe a forma correta. Isso torna o cálculo incrivelmente rápido porque há muito menos seletores para girar do que pixels para pintar.

2. O "Tradutor Universal"

Um dos maiores problemas na pesquisa de fusão é que diferentes programas de computador falam língagens diferentes. Um programa pode fornecer a "pressão", outro a "corrente" e outro o "fator de segurança" (uma medida de estabilidade).

O VEQ atua como um Tradutor Universal. O artigo mostra que o VEQ possui seis rotas de entrada diferentes (como portas USB distintas). Você pode conectar dados de qualquer uma dessas fontes diferentes, e o VEQ os traduz para sua própria linguagem interna para resolver o problema.

• A Alegação: Os autores testaram isso alimentando o mesmo dado perfeito em todas as seis portas. Eles descobriram que, não importava qual porta utilizassem, o VEQ produzia exatamente o mesmo resultado. Isso prova que o tradutor funciona perfeitamente e não introduz erros apenas por você ter usado um cabo de entrada diferente.

3. O Equilíbrio entre "Velocidade vs. Precisão"

O artigo não diz apenas que "é rápido"; ele explora a troca entre velocidade e precisão, como escolher entre um esboço e uma fotografia.

• O Esboço (Baixos Parâmetros): Você usa poucos seletores. É instantâneo (milissegundos) e bom o suficiente para uma visualização rápida, mas pode perder detalhes minúsculos.
• A Foto (Altos Parâmetros): Você usa muitos seletores. Demora um pouco mais (ainda assim, muito rápido, cerca de 19 milissegundos para formas complexas), mas captura a forma com alta precisão.
• O Resultado: Os autores testaram isso em três tipos diferentes de formas de plasma (um formato em "D" padrão, um "H-mode" de alto desempenho e uma forma complexa com um "ponto X"). Eles descobriram que, mesmo com um pequeno número de seletores, o VEQ conseguia recriar as formas complexas com um erro tão pequeno que é quase invisível (menos de 0,2% do tamanho do reator).

4. O "Teste de Estresse" (Onde as rachaduras aparecem)

Os autores foram honestos sobre os limites. Eles verificaram onde o "equilíbrio de forças" (a tensão que segura o plasma) era perfeito e onde não era.

• O Interior: No meio do plasma, o VEQ é excelente. As forças estão equilibradas quase perfeitamente.
• A Borda: Perto da pele externa (o limite), os erros são ligeiramente maiores. Isso ocorre porque o VEQ usa um esqueleto liso e flexível, mas as fronteiras reais do plasma podem ser irregulares ou ter cantos afiados (como um ponto X).
• A Conclusão: O VEQ é ótimo para o "meio" do reator. Se você precisar saber exatamente o que está acontecendo na borda extrema, pode precisar conferir com uma ferramenta mais lenta e detalhada. Mas para a maioria das tarefas de controle, o VEQ é rápido e preciso o suficiente.

5. O Teste de "Transporte"

Finalmente, eles testaram se os pequenos erros na forma do VEQ causariam grandes problemas se você tentasse usá-lo para prever como o calor se move através do plasma (um teste de "transporte").

• O Resultado: Os erros foram ínfimos (menos de 1%). É como se você medisse uma sala com uma régua levemente torta; o erro no tamanho da sala não mudaria sua decisão sobre o tamanho do tapete que deve comprar.

Resumo

VEQ é uma calculadora super-rápida para reatores de fusão. Em vez de mapear cada ponto individual, ele usa um esqueleto matemático flexível para descrever a forma do plasma.

• É rápido: Pode resolver formas complexas em milissegundos.
• É flexível: Aceita dados de muitas fontes diferentes.
• É confiável: Funciona bem para o interior do reator e é preciso o suficiente para a maioria dos sistemas de controle, desde que se fique atento às bordas externas.

Os autores concluem que o VEQ é perfeito para sistemas que precisam perguntar "Como o plasma está agora?" repetidamente, como no controle de reatores em tempo real ou na execução de milhares de simulações para encontrar as melhores condições de operação.

Resumo técnico

Resumo Técnico: VEQ – Um Solucionador de Grad–Shafranov Paramétrico Rápido para Equilíbrios de Tokamak de Fronteira Fixa

Definição do Problema
Fluxos de trabalho de modelagem integrada para tokamaks (ex: IMAS, RAPTOR, SuperCode) exigem consultas repetidas de geometria de equilíbrio, coeficientes métricos e perfis de diagnóstico. Esses fluxos de trabalho diferem de reconstruções de alta fidelidade feitas uma única vez; eles demandam uma representação que seja compacta o suficiente para varreduras de parâmetros e iterações orientadas ao controle, mas que retenha informações de conformação 2D e de balanço de forças suficientemente úteis para serem aplicáveis além de simples fechamentos analíticos. As ferramentas existentes frequentemente caem em dois extremos: fechamentos de engenharia de baixo nível (ex: modelos Miller), que são rápidos, mas carecem de diagnósticos residuais; ou solucionadores de equilíbrio completos (ex: CHEASE, EFIT, DESC), que fornecem saídas nativas de alta fidelidade, mas costumam ser computacionalmente caros ou pesados em termos de infraestrutura para uso em loops internos. Existe uma lacuna para uma representação contínua de fronteira fixa que preserve a conformação 2D útil e os diagnósticos residuais, permanecendo econômica para regeneração.

Metodologia
O artigo apresenta o VEQ (Veloce EQuilibrium), um framework paramétrico compacto projetado para preencher este regime intermediário. A implementação específica testada, VEQPy, é um solucionador de Grad–Shafranov axisimétrico de fronteira fixa.

• Representação Paramétrica: O equilíbrio é definido por um mapa de superfície de fluxo paramétrico $(\rho, \theta) \to [R(\rho, \theta), Z(\rho, \theta)]$. A geometria é parametrizada usando harmônicos de superfície de fluxo do tipo MXH (ângulos poloidais distorcidos) e polinômios de Chebyshev deslocados para perfis radiais. As incógnitas ativas são os coeficientes desses harmônicos e perfis, em vez de valores nodais de um mapa de fluxo poloidal baseado em grade.
• Formulação Variacional: O solucionador impõe um resíduo de Grad–Shafranov projetado via variação. Em vez de resolver a equação de forma forte ponto a ponto, o método restringe as variações admissíveis àquelas geradas pelos coeficientes ativos da representação paramétrica. Isso resulta em um sistema de dimensão finita onde o bloco de forma é uma projeção estruturada variaçãovelmente (condição de Petrov–Galerkin) e os blocos de perfil são fechados por projeções de momentos.
• Rotas de Entrada Flexíveis: Para garantir a compatibilidade com diversos inputs provenientes de estudos analíticos, arquivos de reconstrução e modelos de controle, o VEQ suporta seis rotas de entrada (PF, PP, PI, PJ1, PJ2, PQ). Essas rotas aceitam gradientes de pressão, funções de campo toroidal, gradientes de fluxo poloidal, corrente toroidal encerrada, densidade de corrente ou perfis de fator de segurança. Uma rota de fechamento específica mapeia esses inputs heterogêneos para uma representação de fonte normalizada comum, permitindo que todas as rotas alimentem o mesmo operador de resíduo de dimensão finita.
• Implementação: O solucionador utiliza um solver não linear (método híbrido de Powell do SciPy) com um backend acelerado (Numba). Ele separa a configuração (pré-computação dependente da grade, dados do operador) da fase de solução para permitir consultas de baixa latência repetidas.

Principais Contribuições

1. Representação Paramétrica de Dimensão Finita: Definição de uma representação de superfície de fluxo MXH–Chebyshev compacta e a derivação do resíduo de Grad–Shafranov projetado via variação, incluindo variações convectivas associadas ao movimento da superfície de fluxo.
2. Fechamentos de Fonte Específicos por Rota: Introdução de um framework que mapeia inputs 1D heterogêneos (pressão, corrente, fator de segurança, etc.) para um sistema de resíduo comum, distinguindo a consistência da rota de entrada do erro de aproximação geométrica.
3. Abordagem de Diagnóstico Dual: Relata tanto o resíduo projetado (minimização pelo solver) quanto os resíduos de forma forte amostrados pontualmente. Isso distingue a convergência do sistema reduzido do rastreamento de balanço de forças pontual.
4. Benchmarking: Testes abrangentes em casos Solov'ev, CHEASE e derivados de EFIT G-EQDSK, quantificando a consistência da rota de entrada, erro de reconstrução de superfície de fluxo, trocas entre precisão–custo e latência de solução repetida.

Resultos

• Consistência de Rota: Testes controlados usando inputs suaves e mutuamente compatíveis gerados a partir de um equilíbrio de referência comum demonstraram que todas as seis rotas de entrada convergem para o mesmo equilíbrio de fronteira fixa paramétrico, com diferenças atribuíveis apenas aos caminhos de recuperação numérica (ex: passos de integração/diferenciação).
• Precisão e Latência: O solver foi testado em três casos representativos: um caso Solov'ev em forma de D, um caso CHEASE em modo H e um caso de ponto-X (X-point) de EFIT.
  • Forma de D: Alcançou um erro de forma normalizado pelo raio menor de $1.4 \times 10^{-3}$ com 9 parâmetros ativos e um tempo de mediana de apenas solução de 1.6 ms.
  • Modo H: Alcançou um erro de $1.1 \times 10^{-3}$ com 65 parâmetros e um tempo de mediana de 19 ms.
  • Ponto-X: Alcançou um erro de $1.9 \times 10^{-3}$ com 94 parâmetros e um tempo de mediana de 15 ms.
• Diagnósticos de Resíduo: Diagnósticos de forma forte pontuais revelaram que o enriquecimento da representação ativa melhora primariamente o balanço de forças interno. Nos casos de Modo H e Ponto-X, os valores globais de RMS e de resíduo máximo permaneceram dominados por contribuições próximas à fronteira, indicando que o ansatz de fronteira fixa não pode compensar totalmente desalinhamentos de ajuste de fronteira ou estruturas singulares (como pontos-X) através de variações de forma interna puras.
• Acoplamento de Transporte: Em um teste isolado de acoplamento geometria-transporte 1D, a resposta do perfil de temperatura aos erros de geometria do VEQ permaneceu abaixo de aproximadamente 1%.

Significância e Alegações
O artigo afirma que o VEQ é adequado para consultas repetidas de geometria de equilíbrio em fluxos de trabalho de modelagem integrada e controle, desde que diagnósticos pontuais sejam mantidos para filtrar casos que exijam refinamento de fronteira, correção local ou solucionadores de maior fidelidade.

Os autores posicionam o VE_Q como complementar aos solucionadores de reconstrução de alta fidelidade, não como um substituto universal. Seu valor principal reside em fornecer uma representação contínua e compacta com resíduos projetados explícitos e diagnósticos amostrados que são econômicos de regenerar. O artigo afirma explicitamente que os resultados de tempo refletem a latência de "apenas solução" após a configuração e não devem ser interpretados como um substituto para os tempos totais do pipeline que incluem pré-processamento, I/O e discretização. O método é destinado a preencher a lacuna entre fechamentos analíticos e pipelines nativos de solucionadores completos, ofereção um equilíbrio entre riqueza geométrica e eficiência computacional.