Homepage Inria website

Overall Objectives
Research Program
Application Domains

Section: New Software and Platforms


Finite Element Equilibrium Solver in MATLAB

Keywords: Finite element modelling - Optimal control - Plasma physics

Functional Description: FEEQS.M (Finite Element Equilibrium Solver in Matlab) is a MATLAB implementation of the numerical methods in [Heumann2015] to solve equilibrium problems for toroidal plasmas. Direct and inverse problems for both the static and transient formulations of plasma equilibrium can be solved. FEEQS.M exploits MATLAB‘s evolved sparse matrix methods and uses heavily the vectorization programming paradigm, which results in running times comparable to C/C++ implementations. FEEQS.M complements the production code CEDRES++ in being considered as fast prototyping test bed for computational methods for equilibrium problems. This includes aspects of numerics such as improved robustness of the Newton iterations or optimization algorithms for inverse problems. The latest developments aim at incorporating the resistive diffusion equation.

[Heumann2015]: Heumann, H., Blum, J., Boulbe, C., Faugeras, B., Selig, G., Ané, J.-M., Brémond, S., Grandgirard, V., Hertout, P., Nardon, E.: Quasi-static free-boundary equilibrium of toroidal plasma with CEDRES++: Computational methods and applications. In: Journal of Plasma Physics 81 (2015)