Torque on the Accretion disc

Jamipat

http://www.maths.qmul.ac.uk/~rpn/ASTM735/lecture3.pdf

The diagram on page 26 is the accretion disc.

The torque acting on the inner edge of the ring (the one that has a thickness of dR in the diagram) is

RF_in = -2[itex]\pi[/itex]R³[itex]\nu[/itex]Ʃ[itex]\frac{dΩ}{dR}[/itex]

The torque acting on the outer edge of the ring is

RF_out = [2[itex]\pi[/itex]R³[itex]\nu[/itex]Ʃ[itex]\frac{dΩ}{dR}[/itex]]_R+dR

I would think that the net torque acting would be

T = [2[itex]\pi[/itex]R³[itex]\nu[/itex]Ʃ[itex]\frac{dΩ}{dR}[/itex]]_R+dR - [2[itex]\pi[/itex]R³[itex]\nu[/itex]Ʃ[itex]\frac{dΩ}{dR}[/itex]]_R

= [2[itex]\pi[/itex]R³[itex]\nu[/itex]Ʃ[itex]\frac{dΩ}{dR}[/itex]]dR

but according to equation (66) on page 27 of http://www.maths.qmul.ac.uk/~rpn/ASTM735/lecture3.pdf, the net torque is

[itex]\frac{∂}{∂R}[/itex][2[itex]\pi[/itex]R³[itex]\nu[/itex]Ʃ[itex]\frac{dΩ}{dR}[/itex]]dR

Does anyone know why [itex]\frac{∂}{∂R}[/itex] is included in the expression?

zhermes

Assume you take the limit that [itex] dR \rightarrow 0[/itex] then that is the correct expression, just based on the definition of the derivative.
Conceptually its because the torque is the difference in force between different radii.

Chronos

This sounds suspiciously like a variant on Lorentian aether, which has been thoroughly disproven.

Jamipat

Why does dR [itex]\rightarrow[/itex] 0 makes it a correct expression?

zhermes

Because that's the definition of a derivative
[tex] f'(x) \equiv \lim_{dx \to 0} \frac{ f(x+dx) - f(x) }{ dx } = \lim_{dx \to 0} \frac{ f|_{x+dx} - f|_x }{ dx } [/tex]

Jamipat

If that's the definition of the derivative, shouldn't [itex]\frac{∂}{∂R}[/itex] on the expression be [itex]\frac{d}{dR}[/itex]

zhermes

No, because you're only looking at the explicit R dependence---in your equation f(x) would be just the explicitly R dependent part. So you're just doing the partial derivative.

Jamipat

T_m = [itex]\frac{dj}{dt}[/itex]

where T_m is the torque per unit mass, and j = R²Ω(R) which is the specific angular momentum (angular momentum per unit mass)

It may be shown that

[itex]\frac{dj}{dt}[/itex] = [itex]\frac{∂j}{∂t}[/itex] + (v.[itex]\nabla[/itex])j (convective derivative)

= v_R[itex]\frac{∂j}{∂R}[/itex] [equation (68) of page 27 http://www.maths.qmul.ac.uk/~rpn/ASTM735/lecture3.pdf]

Looking at the convective derivative, I'm guessing that [itex]\frac{∂j}{∂t}[/itex] = 0 due to the derivative of j = R²Ω(R) with respect to t. If that's the case, then surely

[itex]\frac{dj}{dt}[/itex] should also be 0 which makes v_R[itex]\frac{∂j}{∂R}[/itex] equal to 0.

I don't understand why [itex]\frac{dj}{dt}[/itex] isn't 0.

zhermes

The change in specific angular momentum would be zero if there were no torque. If there is torque, its not zero.

The equation you gave, which (at least) looks right---says that [itex] \frac{dj}{dt} = (v \dot \Del ) j = v_R \frac{ \partial j}{\partial R}[/itex] (i think you're right about the partial w.r.t. time being zero---at least if we assume a steady state). But the divergence of the specific angular momentum isn't zero---so neither is the R.H.S.