pravash 1.2.0

Fluid dynamics simulation — SPH, Euler/Navier-Stokes, shallow water, buoyancy, drag, vortex dynamics
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

use pravash::common::{FluidConfig, FluidMaterial};
use pravash::sph::{SphSolver, create_particle_block, max_speed, total_kinetic_energy};

fn main() {
    let mut particles = create_particle_block([0.2, 0.5], [0.3, 0.3], 0.02, 0.001);
    let config = FluidConfig::water_2d();
    let viscosity = FluidMaterial::WATER.viscosity;
    let mut solver = SphSolver::new();

    println!("Pravash — SPH fluid simulation (SphSolver with spatial hash)");
    println!("Particles: {}", particles.len());
    println!("Timestep:  {} s", config.dt);
    println!();

    for frame in 0..100 {
        solver.step(&mut particles, &config, viscosity).unwrap();

        if frame % 10 == 0 {
            let ke = total_kinetic_energy(&particles);
            let ms = max_speed(&particles);
            println!("Frame {frame:>4}  KE: {ke:>10.4}  max_speed: {ms:>8.4}");
        }
    }

    println!();
    println!("Done.");
}