Search Struphy

AdhocTorus

Ad hoc tokamak MHD equilibrium with circular concentric flux surfaces.

AdhocTorusQPsi

Ad hoc tokamak MHD equilibrium with circular concentric flux surfaces.

AdiabaticPhi

Electrostatic potential for adiabatic electrons, computed from … Not used by a model yet.

Building a simulation

How models, domains, grids, and time integration fit together.

CanonicalMaxwellian2D

Canonical Maxwellian distribution function in (η₁, η₂, η₃, v∥, μ) coordinates.

CircularTokamak

Tokamak MHD equilibrium with circular concentric flux surfaces.

Citation and reproducibility

How to cite Struphy and document reproducible simulations.

CoaxialWaveguideElectric_r

Initial-condition perturbation in the norm basis.

CoaxialWaveguideElectric_theta

Initial-condition perturbation in the norm basis.

CoaxialWaveguideMagnetic

Initial-condition perturbation in the norm basis.

Code and parallelization

The architecture and parallel implementation of Struphy.

ColdPlasma

Cold plasma model: electron-fluid current coupled with Maxwell's equations via a cold-plasma Ohm's law.

ColdPlasma

Base class for a distribution as a Dirac-delta in velocity (vth = 0).

ColdPlasmaVlasov

Hybrid cold-plasma model: a cold electron fluid coupled with a hot kinetic species and Maxwell's equations.

Colella

Mapped geometry available for Struphy simulations.

ConstantVelocity

Constant-velocity background equilibrium. Represents a simple fluid equilibrium with a spatially-constant bulk velocity ( ux , uy , uz ) and configurable density/pressure…

Contact us

Ask questions, report bugs, contribute, or contact the maintainers.

Cuboid

Mapped geometry available for Struphy simulations.

CurlCurlSolve

Advances e. Not used by a model yet.

CurrentCoupling5DCurlb

Advances u and energetic_ions. Used by 1 model.

CurrentCoupling5DDensity

Advances u. Used by 1 model.

CurrentCoupling5DGradB

Advances u and energetic_ions. Used by 1 model.

CurrentCoupling6DCurrent

Advances ions and u. Used by 1 model.

CurrentCoupling6DDensity

Advances u. Used by 1 model.

CurrentSheet

Current sheet equilibrium.

DESCunit

Mapped geometry available for Struphy simulations.

DeterministicParticleDiffusion

Diffusion equation discretized with a deterministic particle method.

DriftKineticElectrostaticAdiabatic

Electrostatic drift-kinetic model for a single ion species with adiabatic electrons.

EfieldWeightsCoupling

Advances e and ions. Used by 2 models.

EQDSKequilibrium

Interface to `EQDSK file format `_.

Erf_z

Initial-condition perturbation.

Examples

Reproducible Struphy simulations and benchmark problems.

FaradayExtended

Equations: Faraday's law … Not used by a model yet.

FEEC basics

An introduction to compatible finite element spaces and operators.

GenericCartesianFluidEquilibrium

Generic Cartesian fluid equilibrium with callable fields. This class extends CartesianFluidEquilibrium to allow user-defined callable functions for velocity, pressure, an…

GenericCartesianFluidEquilibriumWithB

Generic Cartesian fluid equilibrium with magnetic field and callable fields. This class extends GenericCartesianFluidEquilibrium to include magnetic field and its gradien…

GuidingCenter

Guiding-center equation for a single species in a static background magnetic field.

GVECunit

Mapped geometry available for Struphy simulations.

GyroMaxwellian2D

A gyrotropic Maxwellian depending on two velocities (v∥, μ), n=2, where v∥ = 𝐯 · b₀ and μ = v⟂²/(2B₀) is the magnetic moment, with B₀ the background magnetic field strength.

GyroMaxwellian2Dvperp

A gyrotropic Maxwellian depending on two velocities (v∥, v⟂), n=2, where v∥ = 𝐯 · b₀ and v⟂ is the radial component of a polar coordinate system perpendicular to the magentic direction b₀.

Hall

Advances b. Used by 1 model.

HasegawaWakatani

Hasegawa-Wakatani equations in 2D for drift-wave turbulence.

HasegawaWakataniStep

Advances n and omega. Used by 1 model.

HollowCylinder

Mapped geometry available for Struphy simulations.

HollowTorus

Mapped geometry available for Struphy simulations.

HomogenSlab

Homogeneous MHD equilibrium.

HomogenSlabITG

Homogenous slab equilibrium with temperature gradient in x, B-field in z.

IGAPolarCylinder

Mapped geometry available for Struphy simulations.

IGAPolarTorus

Mapped geometry available for Struphy simulations.

ImplicitDiffusion

Advances phi. Not used by a model yet.

IncompressibleNavierStokesSPH

Incompressible Navier-Stokes equations discretized with smoothed particle hydrodynamics (SPH).

ITPA_density

Initial-condition perturbation in the physical basis.

JxBCold

Advances j. Used by 2 models.

LinearExtendedMHDuniform

Linear extended MHD with zero-flow equilibrium, for uniform background conditions only.

LinearMHD

Linear ideal MHD with zero-flow equilibrium for magnetohydrodynamic wave propagation.

LinearMHDDriftkineticCC

Hybrid linear ideal MHD coupled with energetic ions (5D drift-kinetic) via the current-coupling scheme.

LinearMHDVlasovCC

Hybrid linear MHD coupled with energetic ions (6D Vlasov) via the current-coupling scheme.

LinearMHDVlasovPC

Hybrid linear MHD coupled with energetic ions (6D Vlasov) via the pressure-coupling scheme.

LinearVlasovAmpereOneSpecies

Linearized Vlasov-Ampère equations for one kinetic species around a Maxwellian background.

LinearVlasovMaxwellOneSpecies

Linearized Vlasov-Maxwell equations for one kinetic species around a Maxwellian background.

Magnetosonic

Advances n, u and p. Used by 4 models.

MagnetosonicUniform

Advances n, u and p. Used by 1 model.

ManufacturedSolutionPotential

Initial-condition perturbation in the physical basis.

ManufacturedSolutionVelocity

Initial-condition perturbation in the physical basis.

ManufacturedSolutionVelocity_2

Initial-condition perturbation in the physical basis.

Maxwell

Maxwell's equations in vacuum for electromagnetic field evolution.

Maxwellian3D

A Maxwellian depending (η₁, η₂, η₃) and on three (n=3) Cartesian velocities.

MaxwellWeakAmpere

Advances e and b. Used by 5 models.

ModesCos

Initial-condition perturbation.

ModesCosCos

Initial-condition perturbation.

ModesCosSin

Initial-condition perturbation.

ModesSin

Initial-condition perturbation.

ModesSinCos

Initial-condition perturbation.

ModesSinSin

Initial-condition perturbation.

Numerical methods

The structure-preserving methods used throughout Struphy.

OhmCold

Advances j and e. Used by 2 models.

Orthogonal

Mapped geometry available for Struphy simulations.

Poisson

Weak discretization of Poisson's equation with a diffusion matrix, stabilization and an optional time-dependent right-hand side.

PoissonAdiabaticGyrokinetic

Advances phi. Used by 1 model.

PoissonSolve

Advances phi. Used by 4 models.

PoweredEllipticCylinder

Mapped geometry available for Struphy simulations.

PressureCoupling6D

Advances u and energetic_ions. Used by 1 model.

PressureLessSPH

Particle discretization of pressureless Euler flow with external forcing.

PushDeterministicDiffusion

Advances var. Used by 1 model.

PushEta

Advances var. Used by 10 models.

PushEtaPC

Advances var. Used by 1 model.

PushGuidingCenterBxEstar

Advances ions. Used by 4 models.

PushGuidingCenterParallel

Advances ions. Used by 3 models.

PushRandomDiffusion

Advances var. Used by 1 model.

PushVinForceField

Advances var. Used by 4 models.

PushVinSPHpressure

Advances fluid. Used by 1 model.

PushVinViscousPotential

Advances fluid. Used by 2 models.

PushVxB

Advances ions. Used by 10 models.

RandomParticleDiffusion

Diffusion equation discretized with a random particle method, via a Wiener process.

RestelliAnalyticSolutionPotential

Initial-condition perturbation in the physical basis.

RestelliAnalyticSolutionVelocity

Initial-condition perturbation in the physical basis.

RestelliAnalyticSolutionVelocity_2

Initial-condition perturbation in the physical basis.

RestelliAnalyticSolutionVelocity_3

Initial-condition perturbation in the physical basis.

ScrewPinch

Straight tokamak (screw pinch) MHD equilibrium for a cylindrical geometry of radius mathematical formulation and length mathematical formulation.

ShafranovDshapedCylinder

Mapped geometry available for Struphy simulations.

ShafranovShiftCylinder

Mapped geometry available for Struphy simulations.

ShafranovSqrtCylinder

Mapped geometry available for Struphy simulations.

Shear_x

Initial-condition perturbation.

Shear_y

Initial-condition perturbation.

Shear_z

Initial-condition perturbation.

ShearAlfven

Shear-Alfvén part of LinearMHD with zero-flow equilibrium.

ShearAlfvenB1

Advances u and b. Used by 1 model.

ShearAlfvenCurrentCoupling5D

Advances u and b. Used by 1 model.

ShearAlfvenPropagator

Advances u and b. Used by 4 models.

ShearedSlab

Sheared slab MHD equilibrium in a cube with side lengths mathematical formulation. Profiles depend on mathematical formulation solely.

ShearFluid

Sheared fluid equilibrium in a cube with side lengths mathematical formulation. Profiles depend on mathematical formulation solely.

Struphy kernels

Compiled numerical routines for splines, finite elements, geometry, linear algebra, and particles.

Team

The maintainers and contributors behind Struphy.

Time integration

Splitting methods and time-stepping schemes available in Struphy.

TimeDependentSource

Advances source. Not used by a model yet.

Tokamak

Mapped geometry available for Struphy simulations.

TokamakManufacturedSolutionVelocity

Initial-condition perturbation in the physical basis.

TokamakManufacturedSolutionVelocity_1

Initial-condition perturbation in the physical basis.

TokamakManufacturedSolutionVelocity_2

Initial-condition perturbation in the physical basis.

TorusModesCos

Initial-condition perturbation.

TorusModesSin

Initial-condition perturbation.

ToyDrift

Electrostatic drift toy model for a single ion species in a given background magnetic field.

TwoFluidQuasiNeutralFull

Advances u, ue and phi. Used by 1 model.

TwoFluidQuasiNeutralToy

Linearized, quasi-neutral two-fluid model with zero electron inertia.

VariationalBarotropicFluid

Barotropic fluid equations discretized with a variational method.

VariationalCompressibleFluid

Fully compressible fluid equations discretized with a variational method.

VariationalDensityEvolve

Advances rho and u. Used by 11 models.

VariationalEntropyEvolve

Advances s and u. Used by 3 models.

VariationalMagFieldEvolve

Advances u and b. Used by 1 model.

VariationalMomentumAdvection

Advances u. Used by 9 models.

VariationalPBEvolve

Advances p, u and b. Used by 3 models.

VariationalPressurelessFluid

Pressure-less fluid equations discretized with a variational method.

VariationalQBEvolve

Advances q, u and b. Used by 3 models.

VariationalResistivity

Advances s and b. Used by 7 models.

VariationalViscosity

Advances s and u. Used by 8 models.

ViscoResistiveDeltafMHD

Delta-f visco-resistive MHD equations discretized with a variational method.

ViscoResistiveDeltafMHD_with_q

Delta-f visco-resistive MHD equations, with the q variable (square root of the pressure), discretized with a variational method.

ViscoResistiveLinearMHD

Linear visco-resistive MHD equations discretized with a variational method.

ViscoResistiveLinearMHD_with_q

Linear visco-resistive MHD equations, with the q variable (square root of the pressure), discretized with a variational method.

ViscoResistiveMHD

Full (non-linear) visco-resistive MHD equations discretized with a variational method.

ViscoResistiveMHD_with_p

Full (non-linear) visco-resistive MHD equations, with the pressure variable, discretized with a variational method.

ViscoResistiveMHD_with_q

Full (non-linear) visco-resistive MHD equations, with the q variable (square root of the pressure), discretized with a variational method.

ViscousEulerSPH

Euler equations with viscosity discretized with smoothed particle hydrodynamics (SPH).

ViscousFluid

Full (non-linear) viscous Navier-Stokes equations discretized with a variational method.

Vlasov

Vlasov equation for a single species in a static background magnetic field.

VlasovAmpereCoupling

Advances e and ions. Used by 3 models.

VlasovAmpereOneSpecies

Vlasov-Ampère system for a single kinetic species in an electric field.