Commit b17b5fe6 authored by Jonathan Lambrechts's avatar Jonathan Lambrechts
Browse files

re-organise file-structure

parent af987753
...@@ -45,7 +45,7 @@ tEnd = 100 ...@@ -45,7 +45,7 @@ tEnd = 100
#numerical parameters #numerical parameters
lcmin = 0.001 # approx r*100 but should match the mesh size lcmin = 0.001 # approx r*100 but should match the mesh size
dt = 1e-4/5 dt = 1e-4*10
alpha = 2.5e-6 alpha = 2.5e-6
epsilon = alpha*lcmin**2 /nu epsilon = alpha*lcmin**2 /nu
print('epsilon',epsilon) print('epsilon',epsilon)
...@@ -75,14 +75,12 @@ def outerBndV(x) : ...@@ -75,14 +75,12 @@ def outerBndV(x) :
print(.05/((ii-11000)%bndLoop + 1)) print(.05/((ii-11000)%bndLoop + 1))
return .05/((ii-11000)%bndLoop + 1) return .05/((ii-11000)%bndLoop + 1)
strong_boundaries = [("Top",2,0.),("TopOut",1,outerBndV),("TopOut",0,0.),("Top",1,outerBndV),("BottomOut",1,outerBndV),("BottomOut",0,0.),("Bottom",1,outerBndV),("Lateral",0,0.)] strong_boundaries = [("Top",2,0.),("TopOut",1,outerBndV),("TopOut",0,0.),("Top",1,outerBndV),("BottomOut",1,outerBndV),("BottomOut",0,0.),("Bottom",1,outerBndV),("Lateral",0,0.)]
fluid = fluid.fluid_problem("mesh.msh",g,nu*rho,rho,epsilon,strong_boundaries) fluid = fluid.fluid_problem("mesh.msh",g,nu*rho,rho,epsilon,strong_boundaries)
fluid.set_particles(p.mass(), p.volume(), p.position(), p.velocity()) fluid.set_particles(p.mass(), p.volume(), p.position(), p.velocity())
fluid.export(outputdir,0,0) fluid.export(outputdir,0,0)
tic = time.clock() tic = time.clock()
forces = g*p.mass() forces = g*p.mass()
fluid.set_particles(p.mass(), p.volume(), p.position(), p.velocity()) fluid.set_particles(p.mass(), p.volume(), p.position(), p.velocity())
......
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