Commit 24fb090b authored by Michel Henry's avatar Michel Henry
Browse files

compute first FluidShear solution

parent 730247af
Pipeline #8767 passed with stages
in 5 minutes and 18 seconds
......@@ -59,13 +59,14 @@ dt = 1e-3 # time step
# dt = 10 # time step
shutil.copy("mesh.msh", outputdir +"/mesh.msh")
outf = 5 # number of iterations between output files
outf = 20 # number of iterations between output files
# Object fluid creation + Boundary condition of the fluid (field 0 is horizontal velocity; field 1 is vertical velocity; field 2 is pressure)
fluid = mbfluid.FluidProblem(2,g,mu,rho,petsc_solver_type="-pc_type lu")
fluid.load_msh("mesh.msh")
fluid.set_wall_boundary("Bottom",velocity=[0,0], pressure = 0)
fluid.set_wall_boundary("Top",velocity=[1.0,0])
time_integration.iterate(fluid,None,10)
# if strong boundary on periodic line, it should be forced on both sides
# fluid.set_strong_boundary("Right",2,0)
......
......@@ -15,6 +15,4 @@ for t in all_tests :
suite.addTest(unittest.defaultTestLoader.loadTestsFromName(t))
r = unittest.TextTestRunner().run(suite)
print(r)
sys.exit(0 if r.wasSuccessful() else 1)
sys.exit(0 if r.wasSuccessful() else 1)
\ 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