Fully developed flow at the inlet

1

Hello everyone,

I have a 2D flow in a rectangular channel. How can I impose a fixed flow rate at the inlet in the sense that the flow will be fully developed at the inlet?

Thank you

2

hello giannoko, you can set a variable of the sine function as the positive flow velocity on your injection boundary, like this:

# Inflow boundary condition for velocity
v_prof = Expression(("2.5e-2*sin(x[1]*pi/0.03)"), degree=2)
bc2 = DirichletBC(W.sub(0).sub(0), v_prof, boundary_markers, 1)
bc3 = DirichletBC(W.sub(0).sub(1), zero,   boundary_markers, 1)