Commit 84f39785 authored by Matthieu Constant's avatar Matthieu Constant
Browse files

depot test

parent 02481d25
Pipeline #5668 passed with stage
in 2 minutes and 9 seconds
......@@ -63,7 +63,7 @@ def genInitialPosition(filename, r, H, ly, lx, rhop) :
p.write_vtk(filename,0,0)
# Define output directory
outputdir = "output"
outputdir = "output1"
if not os.path.isdir(outputdir) :
os.makedirs(outputdir)
......@@ -109,11 +109,10 @@ fluid.set_wall_boundary("Bottom")
fluid.set_wall_boundary("Lateral")
fluid.set_wall_boundary("Top",pressure=0)
# Set location of the particles in the mesh and compute the porosity in each computation cell
fluid.set_particles(p.mass(), p.volume(), p.position(), p.velocity(),p.contact_forces())
fluid.set_particles(p.mass(), p.volume(), p.position(), p.velocity(),p.contact_forces(),init=True)
fluid.export_vtk(outputdir,0,0)
tic = time.time()
fluid.set_particles(p.mass(), p.volume(), p.position(), p.velocity(),p.contact_forces())
#
# COMPUTATION LOOP
......
......@@ -4,8 +4,8 @@ y = 0;
lc = 0.01;
Point(1) = {-L, H, 0,lc};
Point(2) = {-L, -H, 0,lc};
Point(3) = {L, -H, 0,lc};
Point(2) = {-L, 0.2, 0,lc};
Point(3) = {L, 0.2, 0,lc};
Point(4) = {L, H, 0,lc};
Line(1) = {1, 2};
Line(2) = {2, 3};
......
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