Physics-informed Neural Operator Learning for Nonlinear Grad-Shafranov Equation

Abstract

AI for fusion requires bridging a critical ``sim-to-real'' gap: simulation-trained models must generalize reliably under distribution shifts in safety-critical workflows. Focusing on the nonlinear Grad-Shafranov equation (GSE), we develop and analyze a physics-anchored operator-learning framework for fixed-boundary equilibrium prediction. The framework combines data anchors with PDE residual constraints and uses a physics-motivated Transformer-KAN Neural Operator (TKNO) to capture global elliptic coupling and nonlinear source response. Under multi-parameter distribution shifts, our analysis shows that data-only surrogates can develop severe OOD tails, while physics-only training may converge to incorrect solution branches; by combining sparse data anchors with PDE constraints, physics-anchored training reduces worst-tail errors on shape-driven and joint shifts. Mechanism diagnostics indicate that data anchors suppress branch-scale failures, whereas PDE residuals reduce OOD-induced high-frequency error amplification. Evaluated on EXL-50U discharge inputs against the device's operational equilibrium solver, the model achieves close agreement (RMSE $<$ 1.3%) with millisecond-level inference. These results provide a practical route toward physically reliable AI surrogates for fusion workflows. Our code will be available at https://github.com/dsqzhou/physics-anchored-gse.