Commit e6396e46 authored by Matthieu Constant's avatar Matthieu Constant
Browse files

option two fluids when only fluids without grains

parent 0776c0a9
......@@ -89,6 +89,10 @@ def iterate(fluid, particles, dt, min_nsub=1, contact_tol=1e-8, external_particl
raise ValueError("external_particles_forces must have shape (number of particles,dimension!")
if particles is None and fluid is not None:
fluid.implicit_euler(dt)
# For two fluids flows, advance the concentration field with the fluid velocity.
# The number of sub-time steps for the advection is automatically computed.
if fluid.n_fluids() == 2 :
fluid.advance_concentration(dt)
if particles is not None and fluid is None:
_advance_particles(particles,external_particles_forces,dt,min_nsub,contact_tol)
if fluid is not None and particles is not None:
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment