Kinetic Equilibrium Prediction at TCV using RAPTOR and FBT

Este artigo apresenta um novo fluxo de trabalho de previsão de equilíbrio cinético para o TCV, que acopla simulações de transporte RAPTOR com cálculos de equilíbrio inverso FBT para gerar previsões rápidas e precisas de parâmetros críticos e correntes de bobinas, otimizando assim a preparação de disparos do tokamak.

O essencial

Imagine que você é o capitão de um navio futurista, o Tokamak, que tenta viajar através do espaço usando um "sol em uma garrafa" (plasma superaquecido) para gerar energia. O problema é que esse sol é instável, quente demais e quer escapar o tempo todo. Se você não controlar perfeitamente a forma e a posição desse sol, ele bate nas paredes do navio e se apaga (ou pior, causa uma explosão).

Para evitar isso, os cientistas precisam planejar a viagem com antecedência. É aqui que entra o TCV (o navio deles na Suíça) e o novo sistema que o artigo apresenta.

Vamos simplificar o que eles fizeram usando uma analogia de cozinhar um bolo perfeito:

1. O Problema: A Receita Antiga

Antes, quando os cientistas queriam fazer um "bolo de plasma" (um pulso de energia), eles usavam uma receita baseada apenas na intuição e em experiências passadas. Eles diziam: "Coloque X quantidade de farinha (corrente elétrica) e Y quantidade de ovos (aquecimento)".

• O risco: Eles não sabiam exatamente como a massa (o plasma) ia crescer ou mudar de forma enquanto o bolo assava. Se a massa crescesse mais do que o esperado, ela colapsava. Se crescesse menos, o bolo ficava sem graça. Eles precisavam adivinhar a "forma interna" do bolo para saber quanto calor aplicar.

2. A Solução: O "Previsor de Bolos" (RAPTOR + FBT)

Os autores criaram um novo sistema de duas partes que trabalham juntas, como um chef de cozinha e um arquiteto conversando em tempo real antes de começar a cozinhar.

• O Chef (RAPTOR): É um simulador super-rápido. Ele pega a receita (o plano de energia, aquecimento e gás) e diz: "Se fizermos assim, a temperatura no centro vai subir, a massa vai ficar mais densa aqui e mais fina ali". Ele prevê como o "bolo" vai se comportar quimicamente e fisicamente.
• O Arquiteto (FBT): É o especialista em estrutura. Ele olha para o que o Chef previu e diz: "Ok, se a massa vai ficar desse jeito, precisamos ajustar os suportes magnéticos (as bobinas) de tal forma para segurar essa forma específica".

A Mágica (KEP):
Antes, o Chef e o Arquiteto trabalhavam separados. O Chef fazia uma previsão, e o Arquiteto tentava adivinhar como segurar. Agora, eles conversam em loop.

1. O Chef diz: "A massa vai ficar assim".
2. O Arquiteto diz: "Então eu preciso ajustar os suportes para isso".
3. O Chef diz: "Ok, com esses novos suportes, a massa vai mudar um pouco...".
4. Eles repetem isso algumas vezes (em segundos) até chegarem a um plano perfeito onde a massa e os suportes estão em equilíbrio.

Isso é chamado de Previsão de Equilíbrio Cinético (KEP). É como se você pudesse simular o bolo inteiro no computador antes de ligar o forno, garantindo que ele não queime nem desmorone.

3. O Que Eles Descobriram?

Ao testar esse novo sistema em 211 "bolos" (experimentos reais) no TCV, eles viram coisas incríveis:

• Precisão Cirúrgica: O sistema conseguiu prever com muita precisão quantos "amperes" (corrente) cada bobina precisava para segurar o plasma. Antes, eles erravam em algumas centenas de amperes; agora, acertam muito mais.
• Formas Difíceis: Eles conseguiram planejar formas de plasma muito estranhas e complexas (como um "floco de neve" ou formas com triângulos invertidos) que antes eram muito difíceis de controlar. O sistema avisou: "Cuidado, essa forma é instável, precisamos ajustar os suportes aqui e ali".
• Segurança: Ao saber exatamente como o plasma vai se comportar, eles evitam que o "sol" bata nas paredes, o que protege o equipamento e permite que a experiência dure mais tempo.

4. Por Que Isso Importa?

Imagine que você está pilotando um foguete para Marte.

• Sem o sistema: Você decola e torce para os sensores não falharem, ajustando as asas às pressas enquanto o foguete treme.
• Com o sistema: Você roda uma simulação no computador que diz exatamente como o foguete vai reagir a cada turbulência, e você programa os motores para fazerem os ajustes antes mesmo de sair do chão.

Para a fusão nuclear (a energia do futuro), isso é vital. Para construir usinas de energia limpa e infinita (como o projeto ITER), precisamos saber exatamente como controlar o plasma. Este trabalho no TCV é um "laboratório de testes" que está ensinando a todos os outros como fazer essa previsão com antecedência.

Resumo da Ópera:
Eles criaram um "oráculo digital" que une a previsão do comportamento do plasma com o cálculo dos ímãs que o seguram. Isso permite que os cientistas planejem experimentos com muito mais segurança e precisão, transformando a arte de "adivinhar" a física do plasma em uma ciência de "previsão exata". É um passo gigante para tornar a energia de fusão uma realidade prática.

Resumo técnico

Título: Previsão de Equilíbrio Cinético no TCV utilizando RAPTOR e FBT

Autores: C. Contré et al. (Equipe TCV e EUROfusion)
Instituição: École Polytechnique Fédérale de Lausanne (EPFL), Suíça.

---

1. Problema e Motivação

A otimização do desempenho de tokamaks enfrenta limites operacionais rigorosos; ultrapassá-los pode levar a sistemas de proteção ou interrupções (disruptions) que encerram prematuramente os pulsos.

• Desafio Atual: As ferramentas de preparação de pulsos (shot preparation) tradicionais dependem frequentemente de estimativas empíricas ou simplificações dos perfis internos de plasma (pressão e densidade de corrente) para calcular as correntes necessárias nas bobinas de campo poloidal (PF).
• Limitação: A falta de previsão precisa dos perfis cinéticos internos (temperatura, densidade, corrente) antes do experimento resulta em estimativas imprecisas de grandezas críticas como a indutância interna ($l_i$) e o beta poloidal ($\beta_{pol}$). Isso pode causar desalinhamento da forma do plasma, instabilidades verticais e dificuldades em atingir configurações magnéticas complexas (como snowflake ou triangularidade negativa).
• Objetivo: Desenvolver um fluxo de trabalho de "Previsão Primeiro" (Predict-First) que simule a evolução cinética e de transporte do plasma antes do disparo, utilizando essas previsões para refinar o cálculo do equilíbrio magnético e das correntes das bobinas.

2. Metodologia

O trabalho propõe um acoplamento entre duas ferramentas computacionais principais, formando o fluxo de trabalho KEP (Kinetic Equilibrium Prediction):

A. Simulação de Transporte (RAPTOR)

• Função: O RAPTOR é um código de transporte 1.5D rápido que resolve equações de difusão para densidade de corrente, temperaturas (elétrons e íons) e densidade de partículas.
• Modelos Utilizados:
  • Modelo Baseado em Gradiente: Utiliza um regime crítico onde os gradientes logarítmicos de temperatura permanecem "rígidos" (stiff) em uma região central, permitindo prever perfis estáveis.
  • Transição L-H: Implementação de uma lei de escala estendida (Martin scaling) para prever o tempo de transição do modo de baixa para alta confinamento (L para H), adaptada para configurações favoráveis e desfavoráveis (incluindo Double Null e Negative Triangularity).
  • Temperatura de Íons: O modelo resolve as equações de calor para íons e elétrons simultaneamente ou utiliza uma lei de escalonamento empírica, dependendo da precisão necessária.
  • Densidade: Utiliza regras heurísticas baseadas no programa de gás e aquecimento para prever a densidade média de linha ($n_{e,l}$) antes do disparo.

B. Cálculo de Equilíbrio Inverso (FBT)

• Função: O FBT (Free Boundary Tokamak) é um solucionador de equilíbrio estático usado rotineiramente no TCV para determinar as correntes nas bobinas necessárias para manter uma forma e posição desejadas.
• Acoplamento: Em vez de usar perfis de pressão e corrente simplificados (polinomiais) como entrada padrão, o FBT recebe os perfis preditos pelo RAPTOR ($p'$ e $TT'$).
• Iteração: O processo é iterativo. O RAPTOR gera perfis cinéticos baseados nas correntes das bobinas; o FBT usa esses perfis para recalcular as correntes das bobinas necessárias para manter a forma. O sistema converge em poucas iterações (geralmente 2-3), garantindo consistência entre o transporte e o equilíbrio magnético.

3. Contribuições Principais

1. Fluxo de Trabalho KEP: Primeira implementação bem-sucedida de um ciclo de previsão de equilíbrio cinético integrado ao sistema de preparação de pulsos do TCV.
2. Extensão do Modelo de Transporte: Adaptação do modelo de transporte baseado em gradiente para plasmas de Triangularidade Negativa (NT), permitindo prever o desempenho de confinamento superior observado em modos L de NT, sem transição para H-mode.
3. Previsão de Transição L-H: Desenvolvimento de um critério de limiar de potência ($P_{LH}$) que considera a geometria do divertor (configurações favoráveis vs. desfavoráveis e Double Null), melhorando a precisão da previsão do momento da transição.
4. Validação Estatística: Benchmarking do modelo em 211 disparos do TCV, cobrindo uma vasta gama de formas (LSN, SN, Snowflake), modos (L e H) e regimes de aquecimento (Ohmico e NBI).

4. Resultados

• Convergência Rápida: O acoplamento FBT-RAPTOR converge em poucos segundos/minutos por disparo, permitindo sua utilização no tempo real de preparação de experimentos (entre pulsos).
• Precisão de Energia: As previsões de energia térmica total (elétrons e íons) estão dentro de 20% dos valores experimentais reconstruídos para a maioria dos disparos, utilizando apenas o cronograma de pulsos e estimativas de densidade.
• Correção de Correntes das Bobinas: A integração dos perfis preditos altera significativamente as correntes de alimentação (feedforward) calculadas para as bobinas PF. As correções variam de dezenas a centenas de Amperes, especialmente durante fases de alta energia (H-mode) e em configurações complexas.
• Melhoria no Controle de Forma:
  • Caso Padrão (LSN): Ajuste mais preciso da posição do ponto X, reduzindo erros de alinhamento.
  • Caso Complexo (NT + Snowflake): Em configurações de Triangularidade Negativa com divertor Snowflake (notoriamente difíceis de controlar verticalmente), o uso do KEP resultou em uma estabilidade significativamente maior, mantendo os pontos X primário e secundário próximos aos alvos até o final do disparo, evitando instabilidades verticais que ocorreram em preparações tradicionais.
• Indutância Interna ($l_i$): O método fornece estimativas mais realistas de $l_i$ e $\beta_N$, superando as limitações de reconstruções puramente magnéticas que dependem de medições de loops diamagnéticos (que podem ser ambíguas).

5. Significado e Impacto

• Operacional: O KEP permite aos operadores do tokamak ajustar o plano do pulso com base em previsões físicas mais realistas, aumentando a probabilidade de sucesso do experimento e reduzindo o risco de interrupções.
• Científico: Demonstra que a previsão de perfis cinéticos internos é essencial para o cálculo correto de equilíbrio, especialmente em regimes de alto desempenho e configurações magnéticas avançadas.
• Futuro (ITER e além): A metodologia apresentada é um passo crucial para a preparação de experimentos em dispositivos de próxima geração, como o ITER, onde a margem de erro é menor e a complexidade dos cenários exige simulações integradas precisas antes do disparo. O trabalho também aponta para futuras melhorias, como o acoplamento mais rígido com modelos de aquecimento (NBI, ECRH) e transporte de impurezas.

Em resumo, o artigo estabelece um novo padrão para a preparação de experimentos em tokamaks, substituindo estimativas empíricas por uma simulação física integrada que melhora a estabilidade, a precisão da forma e a confiabilidade dos pulsos de plasma.