DeterministicParticleDiffusion

This is a particle discretization of diffusion where particles follow a

deterministic drift derived from the current estimate of the density. It is intended mainly for diffusion-method development and verification rather than plasma dynamics.

This model is appropriate for:

• deterministic particle discretizations of diffusion
• transport benchmarks with positive scalar densities
• numerical comparison against stochastic diffusion methods

PDEs solved by model:

Find $u:\mathbb{R}\times \Omega \to {\mathbb{R}}^{+}$ such that

where $\mathbb{D}:\Omega \to {\mathbb{R}}^{3\times 3}$ is a positive diffusion matrix. At the moment only matrices of the form $D\ast Id$ are implemented, where $D>0$ is a positive diffusion coefficient.

The diffusion coefficient defines the normalization,

No separate plasma velocity scale is used in this model.

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

1. struphy.propagators.push_deterministic_diffusion.PushDeterministicDiffusion

The following scalars are tracked during simulation:

• No default scalar diagnostics are defined by this model.

Create and initialize a deterministic diffusion model:

from struphy.models import DeterministicParticleDiffusion

model = DeterministicParticleDiffusion()

model.hydrogen.var