inject

testcases/injectionWGrains/Q0.05Lyon/inject2fGrains.py

@@ -30,8 +30,7 @@ import time
 import shutil
 import random
 
-outputdir = "outputVid1"
-outputdir1 = "outputVid"
+outputdir = "output"
 if not os.path.isdir(outputdir) :
     os.makedirs(outputdir)
 
@@ -70,11 +69,12 @@ rhop = 1059
 dt = .001                                  # time step
 
 shutil.copy("mesh.msh", outputdir +"/mesh.msh")
-#genInitialPosition(outputdir, N, r, lx, ly, rhop)
+genInitialPosition(outputdir, N, r, lx, ly, rhop)
 friction=0.3
-#lmgc90Interface.scontactTolmgc90(outputdir, 2, 0, friction)
-#p = lmgc90Interface.ParticleProblem(2)
-#p.write_vtk(outputdir,0,0)
+lmgc90Interface.scontactTolmgc90(outputdir, 2, 0, friction)
+p = lmgc90Interface.ParticleProblem(2)
+p.write_vtk(outputdir,0,0)
+
 p = scontact.ParticleProblem(2)
 p.read_vtk(outputdir,0)
 # number of iterations between output files
@@ -82,8 +82,6 @@ outf = 25
 ii = 0
 t = 0
 def outerBndV(x) :
-    return 0.3
-    #print(0.265258*min((6*t**5-15*t**4+10*t**3),1))
     return 0.265258*min((6*t**5-15*t**4+10*t**3),1)
 
 fluid = fluid.FluidProblem(2,g,[nu0*rho0,nu1*rho1],[rho0,rho1],coeff_stab=0.001)