Commit 6a8cbb14 authored by Matthieu Constant's avatar Matthieu Constant
Browse files

inject wbnd

parent 1e659ed5
......@@ -68,7 +68,7 @@ rhop = 2500
# numerical parameters
dt = .0005 # time step
#genInitialPosition(outputdir, N, r, lx, ly, rhop)
genInitialPosition(outputdir, N, r, lx, ly, rhop)
friction=0.3
p = scontact.ParticleProblem(2)
......@@ -84,10 +84,10 @@ fluid = fluid.FluidProblem(2,g,[nu0*rho0,nu1*rho1],[rho0,rho1])
fluid.load_msh("mesh.msh")
fluid.set_weak_boundary("Bottom","wall")
fluid.set_weak_boundary("Lateral","wall")
fluid.set_weak_boundary("Top","pressure",[0, 1])
fluid.set_weak_boundary("Injection","velocity",[0, outerBndV, 1])
fluid.set_wall_boundary("Bottom")
fluid.set_wall_boundary("Lateral")
fluid.set_open_boundary("Top",pressure=0, concentration=1)
fluid.set_open_boundary("Injection",velocity=[0, outerBndV], concentration=1)
fluid.set_particles(p.mass(), p.volume(), p.position(), p.velocity())
......
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