Random Number GodsRNG · Enemy Territory · Battlefield · D&D
Friday BF64 + Sunday D&D

Check the calendar. ET pubs running the new Frostbite rotation all week.

Mon solveur maison vient d’atteindre un équilibre complet de Grad–Shafranov

Permalink Preview Print Science
grim3747
Engineer
Members
46 posts 2 Feb 2020 joined
6 Sep 2026 at 15:32 #1

Salut r/Physics / r/Fusion,

Je bosse depuis plusieurs mois sur un solveur numérique maison, PowerUnit, basé à l’origine sur une architecture Yang–Mills.
Cette semaine, j’ai terminé l’intégration complète de l’équation de Grad–Shafranov, utilisée pour calculer l’équilibre MHD dans les Tokamaks.

Et aujourd’hui… j’ai obtenu mon premier équilibre complet.

🧠 Contexte rapide
L’équation de Grad–Shafranov :

Δ

𝜓
=

𝜇
0
𝑅
2
𝑝

(
𝜓
)

1
2
𝑑
𝐹
2
𝑑
𝜓
décrit l’état stationnaire d’un plasma confiné magnétiquement.
C’est une EDP elliptique non‑linéaire, très sensible aux plateaux rigides, aux instabilités UV et aux stagnations.

Les solveurs pros (EFIT, CHEASE, HELENA, PSI‑Tri) utilisent des techniques avancées pour gérer ça.

J’ai tenté de reproduire une partie de cette logique dans mon propre solveur.

⚙️ Ce que PowerUnit a accompli
Après 30 000+ itérations, le solveur a atteint un état où :

aucune direction de descente (même à 1e‑6) ne réduisait l’action MHD,

les trois directions de “controlled recovery” étaient testées,

tous les pas étaient REJECTED,

la stagnation était détectée proprement,

la procédure “smooth gauge recovery” s’est déclenchée,

et le solveur a finalement marqué :

Status = COMPLETE (MHD equilibrium)

Résidu final :

𝑅
final
=
3.421599476413
×
10
13
Build stable, aucune divergence, aucune explosion UV.

📌 Points techniques
Line search multi‑échelle (5e‑5 → 1e‑6)

Détection automatique de stagnation

Recovery directions 1/3, 2/3, 3/3

Préconditionnement inspiré de
(

Δ
)

1
)

Architecture héritée de mon solveur Yang–Mills

100 super‑cycles complets

Convergence reproductible

grim3747 - NA nights
I take the dead hours.

pupok832
Engineer
Members
37 posts 21 Oct 2018 joined
6 Sep 2026 at 15:41 #2

Slop

pupok832 - Old guard
Jaymod days. Still here.

papa_artema_ebaniy
Field Medic
Members
33 posts 6 Jul 2025 joined
6 Sep 2026 at 15:43 #3

Je t'ai rédigé ici une autre instruction que tu pourras ajouter la prochaine fois avant de publier ton message. Cela permettra à ton LLM de donner une réponse encore plus précise:

"Please summarize the theory in English so that it can be copied into a Reddit post. That means: no abrupt line breaks, no LaTeX equations (use ASCII instead), just what’s allowed in the Reddit Markdown format."

papa_artema_ebaniy - D&D nights
Nat 1s are character development.

1lyshkaqwe
Engineer
Members
43 posts 30 Jun 2022 joined
6 Sep 2026 at 15:51 #4

My hovercraft is full of chatbots.

1lyshkaqwe - NA nights
I take the dead hours.

1lyshkaqwe
Engineer
Members
43 posts 30 Jun 2022 joined
6 Sep 2026 at 15:53 #5

N'encouragez pas la merde digitale.

1lyshkaqwe - NA nights
I take the dead hours.

Sign in to reply.