Commit 7f50552a authored by Matthieu Constant's avatar Matthieu Constant
Browse files

poiseuille avec 2 viscosités diff

parent bf6c4877
......@@ -110,10 +110,6 @@ while t < tEnd :
#Fluid solver
fluid.implicit_euler(dt)
t += dt
s = fluid.solution()
x = fluid.coordinates()
vel = (s[:,0]+s[:,3]-1/(20*nu0*rho)*x[:,1]*(1-x[:, 1]))**2
print('Error', (vel.sum())**.5)
#Output files writting
if ii %outf == 0 :
ioutput = int(ii/outf) + 1
......
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