RandomParticleDiffusion

This model is the stochastic counterpart of the deterministic particle

diffusion model. It is intended for diffusion-method development and comparisons between random-walk and deterministic transport strategies.

This model is appropriate for:

• stochastic particle diffusion benchmarks
• Monte-Carlo transport verification
• comparison against deterministic diffusion solvers

PDEs solved by model:

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

where $D>0$ is a positive diffusion coefficient.

The natural scaling is set by the diffusion coefficient:

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

1. struphy.propagators.push_random_diffusion.PushRandomDiffusion

The following scalars are tracked during simulation:

• No default scalar diagnostics are defined by this model.

Create and initialize a random diffusion model:

from struphy.models import RandomParticleDiffusion

model = RandomParticleDiffusion()

model.hydrogen.var