maj depot wip

testcases/depot-2d/depot.py

@@ -68,7 +68,7 @@ if not os.path.isdir(outputdir) :
     os.makedirs(outputdir)
 
 # Physical parameters
-g = -9.81                                       # gravity
+g = np.array([0,-9.81])                         # gravity
 r = 1e-3                                        # particles radius
 rhop = 1500                                     # particles density
 rho = 1000                                      # fluid density
@@ -80,8 +80,8 @@ dt = 2.5e-3                                     # time step
 tEnd = 100                                      # final time
 
 # Geometrical parameters
-ly = 7e-2                                       # particles area height
-lx = 2e-1                                       # particles area widht
+ly = 1e-1                                       # particles area height
+lx = 4e-1                                       # particles area widht
 H = 0.6                                         # domain height
 
 #
@@ -114,8 +114,7 @@ fluid.export_vtk(outputdir,0,0)
 
 tic = time.time()
 
-G = np.zeros_like(p.velocity())
-G[:,1] = p.mass()[:,0]*g
+G = p.mass()*g
 #
 # COMPUTATION LOOP
 #

testcases/depot-2d/mesh.geo

@@ -4,8 +4,8 @@ y = 0;
 lc = 0.01;
 
 Point(1) = {-L, H, 0,lc};
-Point(2) = {-L, .2, 0,lc};
-Point(3) = {L, .2, 0,lc};
+Point(2) = {-L, -H, 0,lc};
+Point(3) = {L, -H, 0,lc};
 Point(4) = {L, H, 0,lc};
 Line(1) = {1, 2};
 Line(2) = {2, 3};