Commit b8e7e86c authored by Matthieu Constant's avatar Matthieu Constant
Browse files

fix depot 3d

parent 3e767b64
Pipeline #4546 passed with stage
in 1 minute and 31 seconds
......@@ -44,7 +44,7 @@ def genInitialPosition(filename, rhop) :
# Particles structure builder
p = scontact.ParticleProblem(3)
# Load the mesh.msh file specifying physical boundaries names
p.load_msh_boundaries("mesh.msh", ["Top", "Box"])
p.load_msh_boundaries("mesh.msh", ["Top", "Bottom", "Z", "X"])
# Reading the file containing the radii
myRadius = np.genfromtxt('radius.csv', delimiter=',')
r2 = np.amax(myRadius)*1e-6
......@@ -74,7 +74,7 @@ nu = 1e-6 # kinematic viscosity
# Numerical parameters
outf = 100 # number of iterations between output files
dt = 1e-3 # time step
dt = 1e-4 # time step
tEnd = 100 # final time
#
......@@ -95,13 +95,11 @@ ii = 0
fluid = fluid.FluidProblem(3,g,[nu*rho],[rho])
#Set the mesh geometry for the fluid computation
fluid.load_msh("mesh.msh")
fluid.set_strong_boundary("Top",0,0)
fluid.set_strong_boundary("Top",1,0)
fluid.set_strong_boundary("Top",2,0)
fluid.set_strong_boundary("Top",3,0)
fluid.set_strong_boundary("Box",0,0)
fluid.set_strong_boundary("Box",1,0)
fluid.set_strong_boundary("Box",2,0)
fluid.set_strong_boundary("X",0,0)
fluid.set_strong_boundary("Bottom",1,0)
fluid.set_strong_boundary("Z",2,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())
......
......@@ -15,3 +15,4 @@ Physical Surface("Z") = {5,6};
Physical Point("PtFix") = {5};
Characteristic Length{ PointsOf{ Volume{:}; } } = lc;
Mesh.MshFileVersion = 2;
\ No newline at end of file
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