Finding A_n in a Fourier-series Approach to Viscous Fluid Slowing Down

  • Thread starter Thread starter Pigkappa
  • Start date Start date
  • Tags Tags
    Fluid viscous
Pigkappa
Messages
22
Reaction score
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.
 
Physics news on Phys.org
Your Bessel function integral can be easily evaluated using a Table of Integrals. IDK if you have any other difficult integrals, but consulting references like:

Abramowitz & Stegun, "Handbook of Mathematical Functions"

or

Gradshteyn & Ryzhik, "Table of Integrals, Series, and Products"

are very handy to have. After all, real physicists don't spend their time trying to remember all the rules of integrating weird functions when someone has already done the work.

A & S can be found here: http://people.math.sfu.ca/~cbm/aands/

and G & R can be found on the web or in most college libraries.
 
Thread 'Need help understanding this figure on energy levels'
This figure is from "Introduction to Quantum Mechanics" by Griffiths (3rd edition). It is available to download. It is from page 142. I am hoping the usual people on this site will give me a hand understanding what is going on in the figure. After the equation (4.50) it says "It is customary to introduce the principal quantum number, ##n##, which simply orders the allowed energies, starting with 1 for the ground state. (see the figure)" I still don't understand the figure :( Here is...
Thread 'Understanding how to "tack on" the time wiggle factor'
The last problem I posted on QM made it into advanced homework help, that is why I am putting it here. I am sorry for any hassle imposed on the moderators by myself. Part (a) is quite easy. We get $$\sigma_1 = 2\lambda, \mathbf{v}_1 = \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix} \sigma_2 = \lambda, \mathbf{v}_2 = \begin{pmatrix} 1/\sqrt{2} \\ 1/\sqrt{2} \\ 0 \end{pmatrix} \sigma_3 = -\lambda, \mathbf{v}_3 = \begin{pmatrix} 1/\sqrt{2} \\ -1/\sqrt{2} \\ 0 \end{pmatrix} $$ There are two ways...
Back
Top