HasegawaWakatani

This is a reduced 2D drift-wave turbulence model. The electrostatic

potential is recovered from the vorticity at each step, and the coupled density-vorticity dynamics capture resistive drift-wave behavior without resolving full kinetic physics.

This model is appropriate for:

• 2D resistive drift-wave turbulence studies
• reduced electrostatic edge-plasma benchmarks
• testing Poisson-coupled advection-diffusion solvers

PDEs solved by model:

Density equation:

Vorticity equation:

Potential equation:

where $[\varphi ,n]={\partial }_x\varphi {\partial }_{y}n-{\partial }_y\varphi {\partial }_{x}n$, $C=C(x,y)$ and $\kappa$ and $\nu$ are constants (at the moment only $N=1$ is available).

The electrostatic potential is scaled by the thermal-speed unit,

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

1. struphy.propagators.poisson_solve.PoissonSolve
2. struphy.propagators.hasegawa_wakatani_step.HasegawaWakataniStep

The following scalars are tracked during simulation:

• No default scalar diagnostics are defined by this model.

Create and initialize a Hasegawa-Wakatani model:

from struphy.models import HasegawaWakatani

model = HasegawaWakatani()
model.em_fields.phi
model.plasma.density

model.plasma.vorticity