Verlet

The Verlet integrator for SPH (Smoothed Particle Hydrodynamics) implements a velocity Verlet algorithm combined with SPH interactions. This integrator is a wrapper around the UAMMD VerletNVE integrator combined with an SPH potential.

For more details on the underlying method, please refer to the UAMMD VerletNVE documentation and UAMMD SPH documentation.

  • type: SPH, Verlet

  • parameters:

    • timeStep: real: Time step \([time]\)

    • temperature: real: Initial temperature for velocity initialization \([energy]\)

    • support: real: Support radius for SPH kernel \([distance]\)

    • viscosity: real: Kinematic viscosity of the fluid \([distance^2/time]\)

    • gasStiffness: real: Gas stiffness constant \([energy]\)

    • restDensity: real: Rest density of the fluid \([mass/distance^3]\)

Example:

"sphVerlet":{
  "type":["SPH","Verlet"],
  "parameters":{
    "timeStep": 0.01,
    "temperature": 1.0,
    "support": 1.0,
    "viscosity": 0.1,
    "gasStiffness": 1.0,
    "restDensity": 1.0
  }
}

Note

This integrator automatically initializes particle velocities according to the specified temperature if they are not already set.

Warning

This integrator assumes that particle masses are defined in the particle data. The SPH formulation also requires particle densities to be computed and stored.