Hybrid model

LinearMHDVlasovPC

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

Overview

LinearMHDVlasovPC is the pressure-coupling counterpart to the

current-coupling hybrid model. It is targeted at linear problems where energetic-particle pressure anisotropy is the relevant feedback channel on the bulk MHD dynamics.

Use cases

This model is appropriate for:

linear pressure-coupling hybrid studies
energetic-particle pressure feedback on MHD modes
verification of PushEtaPC and pressure-coupling operators

Governing equations

PDEs solved by model:

MHD continuity:

\frac{\partial ρ ~}{\partial t} + \nabla \cdot (ρ_{0} U ~) = 0

MHD momentum:

ρ_{0} \frac{\partial U ~}{\partial t} + \nabla p ~ + \frac{A _{h}}{A _{b}} \nabla \cdot P ~_{h, ⊥} = (\nabla \times B ~) \times B_{0} + J_{0} \times B ~

J_{0} = \nabla \times B_{0}

MHD pressure:

\frac{\partial p ~}{\partial t} + \nabla \cdot (p_{0} U ~) + \frac{2}{3} p_{0} \nabla \cdot U ~ = 0

MHD induction:

\frac{\partial B ~}{\partial t} - \nabla \times (U ~ \times B_{0}) = 0

Energetic-particle Vlasov equation:

\frac{\partial f _{h}}{\partial t} + (v + U ~_{⊥}) \cdot \nabla f_{h} + [\frac{1}{ϵ} v \times (B_{0} + B ~) - \nabla U ~_{⊥} \cdot v] \cdot \frac{\partial f _{h}}{\partial v} = 0

Perpendicular pressure tensor:

P ~_{h, ⊥} = \int v_{⊥} v_{⊥}^{⊤} f_{h} d v

Normalization

Fluid and hot-particle velocities are normalized with the bulk Alfvén

speed. The kinetic pressure tensor is scaled consistently with $A_{h} m_{H} n^v^_{A}^{2}$ , and the hot cyclotron parameter is $ε$ .

Discretization

Time integration is performed by the following propagators (in sequence):

struphy.propagators.push_eta_pc.PushEtaPC
struphy.propagators.push_vxb.PushVxB
struphy.propagators.pressure_coupling_6d.PressureCoupling6D
struphy.propagators.shear_alfven_propagator.ShearAlfvenPropagator
struphy.propagators.magnetosonic.Magnetosonic

Diagnostics

The following scalars are tracked during simulation:

MHD kinetic energy: en_U
Thermal pressure energy: en_p
Magnetic energy: en_B
Energetic-particle kinetic energy: en_f
Total energy: en_tot
Lost particles: n_lost_particles

Example

Create and initialize the linear MHD-Vlasov pressure-coupling model:

from struphy.models import LinearMHDVlasovPC

model = LinearMHDVlasovPC()
model.em_fields.b_field
model.mhd.velocity

model.energetic_ions.var