Commit 3220b5fc authored by Matthieu Constant's avatar Matthieu Constant
Browse files

wbnd

parent 6a8cbb14
......@@ -71,17 +71,14 @@ fluid = fluid.FluidProblem(2,g,nu*rho,rho)
fluid.load_msh("mesh.msh")
fluid.set_strong_boundary("Box",0,0)
fluid.set_strong_boundary("Box",1,0)
fluid.set_strong_boundary("Left",0,0.001)
fluid.set_strong_boundary("Right",2,0)
fluid.set_weak_boundary("Left","Inflow")
fluid.set_weak_boundary("Right","Outflow")
fluid.set_open_boundary("Left",velocity=[0.001, 0])
fluid.set_open_boundary("Right",pressure=0)
fluid.set_wall_boundary("Box")
fluid.set_particles(p.mass(), p.volume(), p.position(), p.velocity())
fluid.export_vtk(outputdir,0,0)
tic = time.time()
while t < tEnd :
#if ii%100==0 and ii != 0:
# fluid.adapt_mesh(0.1,10,8e-4)
fluid.implicit_euler(dt, stab_param=1e-5)
forces = fluid.compute_node_force(dt)
vn = p.velocity() + forces * dt / p.mass()
......
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