add small-depot

paraview_utils/show_problem.py

@@ -101,7 +101,7 @@ if isfile(fluidFile):
     fluidDisplay = showFluid(fluid, renderView ,animationScene)
     if isfile(particleProblemFile):
         Hide(fluid, renderView)
-        showFluidVelocity(fluid, renderView)
+        # showFluidVelocity(fluid, renderView)
     else:
         animationScene.GoToLast()
         fluidDisplay.RescaleTransferFunctionToDataRange(False, True)

testcases/depot-2d/depot-small.py

@@ -75,7 +75,7 @@ rho = 1000                                      # fluid density
 nu = 1e-6                                       # kinematic viscosity
 
 # Numerical parameters
-outf = 100                                        # number of iterations between output files
+outf = 2                                        # number of iterations between output files
 dt = 1e-3                                     # time step
 tEnd = 100                                      # final time
 
@@ -118,7 +118,7 @@ tic = time.time()
 # COMPUTATION LOOP
 #
 while t < tEnd :
-    time_integration.iterate(fluid, p, dt, min_nsub=5, external_particles_forces=g*p.mass(), use_predictor_corrector=True)
+    time_integration.iterate(fluid, p, dt, min_nsub=5, external_particles_forces=g*p.mass(), use_predictor_corrector=False, check_residual_norm=1e-3)
     t += dt
 
     # Output files writting

testcases/depot-2d/depot.py

@@ -63,7 +63,7 @@ def genInitialPosition(filename, r, H, ly, lx, rhop) :
     p.write_vtk(filename,0,0)
 
 # Define output directory
-outputdir = "output"
+outputdir = "output_incompressible"
 if not os.path.isdir(outputdir) :
     os.makedirs(outputdir)
 
@@ -75,7 +75,7 @@ rho = 1000                                      # fluid density
 nu = 1e-6                                       # kinematic viscosity
 
 # Numerical parameters
-outf = 100                                        # number of iterations between output files
+outf = 2                                        # number of iterations between output files
 dt = 1e-3                                     # time step
 tEnd = 100                                      # final time
 
@@ -102,7 +102,7 @@ ii        =  0
 #
 # FLUID PROBLEM
 #
-fluid = fluid.FluidProblem(2,g,[nu*rho],[rho],drag_in_stab=1,solver="petsc",solver_options="-pc_type lu",usolid=True)
+fluid = fluid.FluidProblem(2,g,[nu*rho],[rho],drag_in_stab=1,solver="petsc4py",solver_options="-pc_type lu",usolid=True)
 # Set the mesh geometry for the fluid computation
 fluid.load_msh("mesh.msh")
 fluid.set_wall_boundary("Bottom")

testcases/depot-2d/mesh.geo

-L = 0.1;
-H = 0.1;
+L = 0.2;
+H = 0.3;
 y = 0;
 lc = 0.01;