Shell balances in cylindrical coordinates

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
2 replies · 5K views
MexChemE
Messages
237
Reaction score
54
Hello, PF! I have some doubts about setting up shell balances in a cylindrical geometry. Consider a fluid flowing down a vertical pipe. In order to perform the momentum balance, we take a cylindrical (annular) shell of length L and width Δr. The analysis of such system can be found in chapter 2 of BSL's Transport Phenomena (section 2.3).

Now, in the book, the rate of momentum entering the shell in the axial direction through z = 0 is given by
[tex](2\pi r\Delta r) \left. \phi_{zz} \right|_{z=0}[/tex]
Where apparently [itex]2\pi r\Delta r[/itex] is the cross-sectional area of the annular shell. This confuses me; if I calculate the area by integrating over the annular cross-section, I get this
[tex]A = \int_r^{r+ \Delta r} \int_0^{2 \pi} r \ d \theta dr = \pi (2r \Delta r + (\Delta r)^2)[/tex]
I don't understand why was the [itex](\Delta r)^2[/itex] term left out from the area expression in the book.

Also, after setting up a shell balance in cartesian coordinates, we have to divide the entire balance equation by the entire volume of the shell (see section 2.2 of BSL). Apparently, this is not the case in cylindrical coordinates. Instead of dividing by 2πrΔrL, the authors divided by 2πΔrL; why did they leave out r?

Thanks in advance for any input!
 
Engineering news on Phys.org
MexChemE said:
Hello, PF! I have some doubts about setting up shell balances in a cylindrical geometry. Consider a fluid flowing down a vertical pipe. In order to perform the momentum balance, we take a cylindrical (annular) shell of length L and width Δr. The analysis of such system can be found in chapter 2 of BSL's Transport Phenomena (section 2.3).

Now, in the book, the rate of momentum entering the shell in the axial direction through z = 0 is given by
[tex](2\pi r\Delta r) \left. \phi_{zz} \right|_{z=0}[/tex]
Where apparently [itex]2\pi r\Delta r[/itex] is the cross-sectional area of the annular shell. This confuses me; if I calculate the area by integrating over the annular cross-section, I get this
[tex]A = \int_r^{r+ \Delta r} \int_0^{2 \pi} r \ d \theta dr = \pi (2r \Delta r + (\Delta r)^2)[/tex]
I don't understand why was the [itex](\Delta r)^2[/itex] term left out from the area expression in the book.
They are looking at a thin annular region where the ##(\Delta r)^2## is going to be negligible. They are going to be dividing by 2πΔrL and letting Δr approach zero, so the ##(\Delta r)^2## contribution will drop out anyway.
Also, after setting up a shell balance in cartesian coordinates, we have to divide the entire balance equation by the entire volume of the shell (see section 2.2 of BSL). Apparently, this is not the case in cylindrical coordinates. Instead of dividing by 2πrΔrL, the authors divided by 2πΔrL; why did they leave out r?

Thanks in advance for any input!
Maybe they liked the form of the resulting equation better. It doesn't really matter. If they divided by r, the resulting equation would be equivalent.

Chet
 
Chestermiller said:
They are going to be dividing by 2πΔrL and letting Δr approach zero, so the ##(\Delta r)^2## contribution will drop out anyway.
Got it. I did think that at first, but the limit had not been taken at that point, so the confusion arised.
Chestermiller said:
Maybe they liked the form of the resulting equation better. It doesn't really matter. If they divided by r, the resulting equation would be equivalent.
I checked for myself and the equations ended up being the same. Personally, I prefer dividing by the whole volume because it makes more sense to me.

Thanks for your input, Chet!