Commit 9a36dd3e authored by Matthieu Constant's avatar Matthieu Constant
Browse files

change testcases

parent 176a493a
Pipeline #4808 passed with stage
in 22 seconds
......@@ -75,14 +75,14 @@ s = fluid.solution()
c = np.ndarray((fluid.n_nodes()))
x = fluid.coordinates()
c[:] = 0
#c[np.logical_and(np.abs(x[:,0])< R,np.abs(x[:,1])<R)] = 1
for i in range(len(x[:,0])):
z = (x[i,0])**2+(x[i,1])**2
R1 = (0.7*R)**2
R2 = (1.3*R)**2
c[:] = 1
c[np.logical_and(np.abs(x[:,0])< R,np.abs(x[:,1])<R)] = 0
#for i in range(len(x[:,0])):
# z = (x[i,0])**2+(x[i,1])**2
# R1 = (0.7*R)**2
# R2 = (1.3*R)**2
# #s[i,1] = -0.01
c[i] = min(max(0,1/(R2-R1)*z-R1/(R2-R1)),1)
# c[i] = min(max(0,1/(R2-R1)*z-R1/(R2-R1)),1)
fluid.set_concentration_cg(c)
......
......@@ -70,7 +70,7 @@ 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)
if use_lmgc:
friction=0.3
lmgc90Interface.scontactTolmgc90(outputdir, 2, 0, friction)
......
L = .136;
H = .25;
l = 2e-3;
l = 4e-3;
lc = 0.01;
f = 4;
......
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