Density-independent SPH solver. More...

#include <DensityIndependentSolver.h>

Inheritance diagram for DensityIndependentSolver:

Public Member Functions
	DensityIndependentSolver (IScheduler &scheduler, const RunSettings &settings)

	~DensityIndependentSolver ()

virtual void	integrate (Storage &storage, Statistics &stats) override
	Computes derivatives of all time-dependent quantities. More...

virtual void	create (Storage &storage, IMaterial &material) const override
	Initializes all quantities needed by the solver in the storage. More...

Public Member Functions inherited from ISolver
virtual void	collide (Storage &UNUSED(storage), Statistics &UNUSED(stats), const Float UNUSED(dt))
	Detects the collisions and computes new positions of particles. More...

Public Member Functions inherited from Polymorphic
virtual	~Polymorphic ()

Detailed Description

Density-independent SPH solver.

Uses solver of Saitoh & Makino [14]. Instead of density and specific energy, independent variables are energy density (q) and internal energy of i-th particle (U). Otherwise, the solver is similar to SummationSolver; the energy density is computed using direct summation by self-consistent solution with smoothing length.

Attention: Works only for ideal gas EoS!

Definition at line 24 of file DensityIndependentSolver.h.

Constructor & Destructor Documentation

◆ DensityIndependentSolver()

DensityIndependentSolver::DensityIndependentSolver	(	IScheduler &	scheduler,
		const RunSettings &	settings
	)

explicit

Definition at line 85 of file DensityIndependentSolver.cpp.

◆ ~DensityIndependentSolver()

DensityIndependentSolver::~DensityIndependentSolver ( )

default

Member Function Documentation

◆ create()

void DensityIndependentSolver::create	(	Storage &	storage,
		IMaterial &	material
	)		const

overridevirtual

Initializes all quantities needed by the solver in the storage.

When called, storage already contains particle positions and their masses. All remaining quantities must be created by the solver. The function is called once for every body in the run. The given storage is guaranteed to be homogeneous; it contains only a single material.

Parameters

storage	Particle storage that shall be modified as needed by the solver.
material	Material containing parameters of the body being created. The solver can also set up necessary timestepping parameters of the material (ranges and minimal values of quantities).

Implements ISolver.

Definition at line 203 of file DensityIndependentSolver.cpp.

◆ integrate()

void DensityIndependentSolver::integrate	(	Storage &	storage,
		Statistics &	stats
	)

overridevirtual

Computes derivatives of all time-dependent quantities.

The solver can also modify the quantities arbitrarily. It is however not recommended to perform the integration in the solver (using the time step stored in Statistics) as this is a job for timestepping. The solver can modify quantities using boundary conditions, inter-quantity relationships (such as the summation equation for density in SPH), clamping of values etc. It is also possible to add or remove particles in the storage and modify materials. Threads running concurrently with the solver must assume the solver can modify the storage at any time, so accessing the storage from different threads is only allowed before or after integrate; the is no locking for performance reasons.

All highest order derivatives are guaranteed to be set to zero when integrate is called (this is a responsibility of ITimeStepping implementation).