Commit e712e4d6 authored by Michel Henry's avatar Michel Henry
Browse files

Merge branch 'jimmy_release' of...

Merge branch 'jimmy_release' of https://git.immc.ucl.ac.be/fluidparticles/migflow into jimmy_release
parents 1e4098ec 4624dc2a
Pipeline #9501 passed with stages
in 3 minutes and 49 seconds
......@@ -160,7 +160,7 @@ tic = time.time()
#
while t < tEnd :
if ii >= 10:
genMorePolymers(fluid, p, joints, polymers, r_max, 1.6, 3*r_max, 40*r_max, 2e-3, n_chain, 0.2)
genMorePolymers(fluid, p, joints, polymers, r_max, 2.0, 3*r_max, 40*r_max, 2e-3, n_chain, 0.2)
joints, polymers = removeUselessParticls(fluid, p, joints, polymers, L, n_chain)
Gp = np.zeros((p.n_particles(),2))
......
......@@ -66,7 +66,6 @@ r_max = 1e-3 # particles max radius
# Numerical parameters
outf = 50 # number of iterations between output files
dt = 5e-5 # time step
dt = 5e-3 # time step
tEnd = 30 # final time
# Geometrical parameters
......@@ -80,7 +79,7 @@ lx = 0.1 # particles area widht
H = 0.1 # domain height
W = 1.0
v = 2.3 # tangent velocity
v = 1.3 # tangent velocity
epsilon = 8e-2 # boundary layer size
noise_force = 5 # Brownian motion
noise_variance = 25*pi/180
......@@ -244,8 +243,7 @@ try:
fluid.set_strong_boundary("Inner", 0, lambda x : (((x[:,1]-ic_position[1])/ic_r))*v)
fluid.set_strong_boundary("Inner", 1, lambda x : (-((x[:,0]-ic_position[0])/ic_r))*v)
# time_integration.iterate(fluid, p, dt, min_nsub=5, external_particles_forces=Gp)
time_integration.iterate(fluid, None, dt, min_nsub=5, external_particles_forces=Gp)
time_integration.iterate(fluid, p, dt, min_nsub=5, external_particles_forces=Gp)
t += dt
# Output files writting
......
Markdown is supported
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