Pigkappa
- 22
- 1
Homework Statement
A fluid of viscosity \nu is rotating with uniform angular velocity \Omega inside a cylinder of radius a that is also rotating. At time t = 0, the cylinder is brought to a rest. The circular motion of the fluid is gradually slowed down due to the viscosity; show that u_\theta (r, t) satisfies the following relation with a Fourier-series approach:
u_\theta (r,t) = - 2 \Omega a \sum_{n \geq 1}{\frac{J_1(\lambda_n r / a)}{\lambda_n J_0(\lambda_n)} \exp{(- \lambda_n^2 \frac{\nu t}{a^2}})}
where \lambda_n is the n-th root of the Bessel function J_1.
I'll write u rather than u_\theta from now on for simplicity.
Homework Equations
Boundary condition u(a, t) = 0
Initial condition u(r, 0) = \Omega r
The fluid equations eventually give \frac{\partial u}{\partial t} = \nu (\frac{\partial^2 u}{\partial r^2} + \frac{1}{r} \frac{\partial u}{\partial r} - \frac{u}{r^2})
The Attempt at a Solution
I won't go through all the details as they take a lot to write. After a standard separation of variables u(r, t) = A(r )*B(t) an equation for A is found that can be brought in the form that has J_1 as its solution. I arrived at:
u(r, t) = \sum_{n \geq 1}{A_n J_1(\lambda_n r / a) \exp{(- \lambda_n^2 \frac{\nu t}{a^2})}}
So I need to find each A_n by using the initial condition u(r, 0) = \Omega r. I don't seem to be able to do it. I find myself doing integrals like:
\int_0^1{x^2 J_1(\lambda_n x) dx}
which I can't do. Can you please show the right steps to take to find the A_n? It might be something very easy but I haven't dabbled with Bessel's function as much as I should have in the past.
4. Where the problem comes from
This is a worked out example in Acheson's book Elementary Fluid Dynamics. He just says that the result comes from a standard Fourier-type analysis and I was trying to re-derive it.