Fluid
Maxwell, MHD, and multi-fluid systems using compatible finite elements.
Structure-preserving plasma physics in Python.
Struphy is an open-source framework for solving fluid, kinetic, and hybrid plasma models with geometric finite elements and Lagrangian particles.
What "structure-preserving" means
Many plasma physics PDEs carry exact invariants: energy, mass, and the divergence-free condition on magnetic fields among them. Standard discretizations can violate these invariants over long simulations, letting numerical error accumulate. Struphy instead builds its finite-element and particle discretizations from a discrete de Rham complex, so these structural properties hold exactly at the discrete level, not only in the limit of vanishing grid size.
Models
A common numerical framework supports several descriptions of plasma dynamics while preserving their underlying geometric structure.
Kinetic
Vlasov and drift-kinetic models represented by Lagrangian particles.
Hybrid
Coupled field-particle descriptions, including MHD-kinetic models.
Quick start
Install Struphy and compile its performance-critical kernels.
Continue with the installation guide andtutorials.
$ pip install struphy
$ struphy compilePublications
A selection of papers describing the methods and algorithms behind Struphy. See thefull publication list for more.
- The linearized Vlasov–Maxwell system as a Hamiltonian system
- Variational discretizations of ideal magnetohydrodynamics in smooth regime using structure-preserving finite elements
- Geometric Particle-In-Cell discretizations of a plasma hybrid model with kinetic ions and mass-less fluid electrons
- High-Order Structure-Preserving Algorithms for Plasma Hybrid Models
- Magneto-hydrodynamic eigenvalue solver for axisymmetric equilibria based on smooth polar splines
- MHD-kinetic hybrid code based on structure-preserving finite elements with particles-in-cell