Struct SphSimulation

Source

pub struct SphSimulation {Show 17 fields
    pub particles: Vec<SphParticle>,
    pub h: f32,
    pub rest_density: f32,
    pub stiffness: f32,
    pub gamma: f32,
    pub viscosity: f32,
    pub surface_tension: f32,
    pub gravity: Vec3,
    pub max_dt: f32,
    pub min_dt: f32,
    pub dt: f32,
    pub cfl_factor: f32,
    pub speed_of_sound: f32,
    pub xsph_epsilon: f32,
    pub boundary_repulsion: f32,
    pub domain_min: Vec3,
    pub domain_max: Vec3,
    /* private fields */
}

Expand description

Weakly Compressible SPH (WCSPH) fluid simulation.

Reference: Becker & Teschner (2007). Uses cubic spline kernel, Tait equation of state, XSPH velocity correction, surface tension (Akinci 2013).

Fields§

§particles: Vec<SphParticle>§h: f32

Smoothing length h.

§rest_density: f32

Reference rest density (kg/m³).

§stiffness: f32

Stiffness for Tait EOS (pressure term).

§gamma: f32

Adiabatic index gamma.

§viscosity: f32

Dynamic viscosity coefficient.

§surface_tension: f32

Surface tension coefficient (Akinci).

§gravity: Vec3

Gravity.

§max_dt: f32

Maximum timestep (CFL).

§min_dt: f32

Minimum timestep.

§dt: f32

Current adaptive dt.

§cfl_factor: f32

CFL number.

§speed_of_sound: f32

Speed of sound (for Tait EOS).

§xsph_epsilon: f32

XSPH epsilon (velocity smoothing).

§boundary_repulsion: f32

Boundary repulsion force strength.

§domain_min: Vec3

Domain bounding box.

§domain_max: Vec3

Implementations§

Source §

impl SphSimulation

Source

pub fn new(h: f32, rest_density: f32) -> Self

Source

pub fn add_cube( &mut self, center: Vec3, half_extent: Vec3, spacing: f32, mass: f32, )

Add a cubic block of particles.

Source

pub fn step(&mut self)

Step the simulation by one adaptive timestep.

Source

pub fn step_with_dt(&mut self, dt: f32)

Step with a fixed timestep.

Source

pub fn advance(&mut self, total_time: f32)

Step for total_time seconds using adaptive CFL.

Source

pub fn particle_count(&self) -> usize

Source

pub fn total_kinetic_energy(&self) -> f32

Total kinetic energy.

Source

pub fn average_density(&self) -> f32

Average density.

Source

pub fn sample_pressure(&self, pos: Vec3) -> f32

Compute pressure at an arbitrary world point via SPH interpolation.

Auto Trait Implementations§

§

impl UnwindSafe for SphSimulation

Blanket Implementations§

Source §

impl<T> Any for T
where T: 'static + ?Sized,

Source §

fn type_id(&self) -> TypeId

Gets the TypeId of self. Read more

Source §

impl<T> Borrow<T> for T
where T: ?Sized,

Source §

fn borrow(&self) -> &T

Immutably borrows from an owned value. Read more

Source §

impl<T> BorrowMut<T> for T
where T: ?Sized,

Source §

fn borrow_mut(&mut self) -> &mut T

Mutably borrows from an owned value. Read more

Source §

impl<T> Downcast for T
where T: Any,

Source §

fn into_any(self: Box<T>) -> Box<dyn Any>

Convert Box<dyn Trait> (where Trait: Downcast) to Box<dyn Any>. Box<dyn Any> can then be further downcast into Box<ConcreteType> where ConcreteType implements Trait.

Source §

fn into_any_rc(self: Rc<T>) -> Rc<dyn Any>

Convert Rc<Trait> (where Trait: Downcast) to Rc<Any>. Rc<Any> can then be further downcast into Rc<ConcreteType> where ConcreteType implements Trait.

Source §

fn as_any(&self) -> &(dyn Any + 'static)

Convert &Trait (where Trait: Downcast) to &Any. This is needed since Rust cannot generate &Any’s vtable from &Trait’s.

Source §

fn as_any_mut(&mut self) -> &mut (dyn Any + 'static)

Convert &mut Trait (where Trait: Downcast) to &Any. This is needed since Rust cannot generate &mut Any’s vtable from &mut Trait’s.

Source §

impl<T> DowncastSync for T
where T: Any + Send + Sync,

Source §

fn into_any_arc(self: Arc<T>) -> Arc<dyn Any + Sync + Send>

Convert Arc<Trait> (where Trait: Downcast) to Arc<Any>. Arc<Any> can then be further downcast into Arc<ConcreteType> where ConcreteType implements Trait.

Source §