Fluid model

IncompressibleNavierStokesSPH

Incompressible Navier-Stokes equations discretized with smoothed particle hydrodynamics (SPH).

Overview

IncompressibleNavierStokesSPH is the particle-based SPH fluid model...

This model is appropriate for:

Incompressible Navier-Stokes:

ρ (\partial_{t} u + u \cdot \nabla u) = - \nabla p - \nabla \cdot σ \nabla \cdot u = 0 .

where $σ$ is the viscous stress tensor.

The viscous stress tensor for a Newtonian fluid is given by:

σ = - μ (\nabla u + (\nabla u)^{T} - \frac{2}{3} (\nabla \cdot u) I)

where $μ$ is the dynamic (shear) viscosity and $I$ is the identity tensor.

The characteristic speed is thermal:

u^= v^_{th} .

Time integration is performed by the following propagators (in sequence):

struphy.propagators.push_eta.PushEta
struphy.propagators.push_vxb.PushVxB (if with_B0 is True)
struphy.propagators.push_vin_viscous_potential.PushVinViscousPotential (if with_viscosity is True)
struphy.propagators.poisson_solve.PoissonSolve
struphy.propagators.push_vin_efield.PushVinEfield

Create and initialize the viscous Euler SPH model:

from struphy.models import ViscousEulerSPH

model = ViscousEulerSPH()

model.euler_fluid.var