Struct salva2d::solver::DFSPHSolver

pub struct DFSPHSolver<KernelDensity: Kernel = CubicSplineKernel, KernelGradient: Kernel = CubicSplineKernel> {
    pub min_pressure_iter: usize,
    pub max_pressure_iter: usize,
    pub max_density_error: Real,
    pub min_divergence_iter: usize,
    pub max_divergence_iter: usize,
    pub max_divergence_error: Real,
    /* private fields */
}

A DFSPH (Divergence Free Smoothed Particle Hydrodynamics) pressure solver.

Fields

min_pressure_iter: usize

Minimum number of iterations that must be executed for pressure resolution.

max_pressure_iter: usize

Maximum number of iterations that must be executed for pressure resolution.

max_density_error: Real

Maximum acceptable density error (in percents).

The pressure solver will continue iterating until the density error drops bellow this threshold, or until the maximum number of pressure iterations is reached.

min_divergence_iter: usize

Minimum number of iterations that must be executed for divergence resolution.

max_divergence_iter: usize

Maximum number of iterations that must be executed for divergence resolution.

max_divergence_error: Real

Maximum acceptable divergence error (in percents).

The pressure solver will continue iterating until the divergence error drops bellow this threshold, or until the maximum number of pressure iterations is reached.

Implementations

impl<KernelDensity, KernelGradient> DFSPHSolver<KernelDensity, KernelGradient>

where KernelDensity: Kernel, KernelGradient: Kernel,

pub fn new() -> Self

Initialize a new DFSPH pressure solver.

Trait Implementations

impl<KernelDensity, KernelGradient> PressureSolver for DFSPHSolver<KernelDensity, KernelGradient>

where KernelDensity: Kernel, KernelGradient: Kernel,

fn init_with_fluids(&mut self, fluids: &[Fluid])

Initialize this solver with the given fluids.

fn init_with_boundaries(&mut self, _boundaries: &[Boundary])

Initialize this solver with the given boundaries.

fn predict_advection( &mut self, timestep: &TimestepManager, kernel_radius: Real, contact_manager: &ContactManager, gravity: &Vector<Real>, fluids: &mut [Fluid], boundaries: &[Boundary] )

Predicts advection with the given gravity.

fn evaluate_kernels( &mut self, kernel_radius: Real, contact_manager: &mut ContactManager, fluids: &[Fluid], boundaries: &[Boundary] )

Evaluate the SPH kernels for all the contacts in contact_manager.

fn compute_densities( &mut self, contact_manager: &ContactManager, fluids: &[Fluid], boundaries: &mut [Boundary] )

Compute the densities of all the boundary and fluid particles.

fn step( &mut self, counters: &mut Counters, timestep: &mut TimestepManager, gravity: &Vector<Real>, contact_manager: &mut ContactManager, kernel_radius: Real, fluids: &mut [Fluid], boundaries: &[Boundary] )

Solves pressure and non-pressure force for the given fluids and boundaries.