순환 운동량과 Routh 축약

보존 운동량과 고정 좌표 구분하기

순환 좌표를 계산에서 제거한다고 그 좌표를 물리적으로 정지시키는 것은 아닙니다. 각도가 순환 좌표이면 정준 운동량이 보존됩니다. 하지만 반지름이 변하면 각속도도 변할 수 있습니다. 축약한 방정식은 이 효과를 유지해야 합니다.

이 단원에서는 부분 Legendre 변환으로 Routhian을 유도합니다. 보존 운동량을 Lagrangian에 바로 대입한 결과, 모터로 각속도를 고정한 다른 실험과 비교합니다. 세 연산은 서로 다른 질문에 답합니다. 한 반지름에서 결과가 같다는 사실만으로 같은 모델이라고 판정할 수 없습니다.

계산 모델은 정확한 유리수 상태와 반지름별 곡선을 평가합니다. 시간 궤적을 적분하지 않습니다. 그림의 반지름은 m 단위의 배치 매개변수이며 물리 시간이 아닙니다.

정규적인 역학 모델 정하기

반지름 $r>0$의 단위는 m, 방사 속도 $v_r$의 단위는 m/s, 각속도 $\omega$의 단위는 s⁻¹입니다. 다음 모델을 사용합니다.

$$L=\frac m2(v_r^2+r^2\omega^2)+\frac{\beta r^2\omega}{2}-U(r), \qquad U(r)=\frac{k r^2}{2}+\frac{a r^4}{4}.$$

양의 질량 $m$은 kg, 강성 $k$는 N/m, 사차 계수 $a$는 N/m³, $\beta$는 kg/s 단위입니다. L의 모든 항은 J입니다. 속도에 선형인 항은 자이로스코픽 결합이며 물리적 에너지에서는 상쇄됩니다. 이 모델을 사용하기 위해 전자기학의 전하 부호나 자기장 해석을 먼저 가정할 필요는 없습니다.

질량·강성·사차 계수는 정수 설정을 사용합니다. 운동량과 beta 코드는 입력/2이며 각각 Js, kg/s 단위입니다. 반지름 코드는 입력/10 m이고 방사 속도와 강제 각속도 코드는 입력/2를 해당 SI 단위로 사용합니다. 반지름0은 특이한 경계를 드러내기 위해서만 선택할 수 있습니다.

보존되는 정준 운동량 확인하기

L에는 각도 자체가 없으므로 각도에 대응하는 정준 운동량은

$$p_\theta=m r^2\omega+\frac{\beta r^2}{2}=\ell$$

로 보존됩니다. 상수 $\ell$의 단위는 Js입니다. 반면 기계적 각운동량은 $m r^2\omega=\ell-\beta r^2/2$입니다. beta가 비영이고 반지름이 변하면 순환 좌표가 있어도 기계적 각운동량은 일정하지 않을 수 있습니다.

회전 대칭이면 각운동량이 보존된다는 문장을 적용할 때 어떤 운동량을 말하는지 지정해야 합니다. 선택한 Lagrangian의 속도 미분이 정준 운동량을 정합니다. 배치에 의존하는 속도 선형 항은 정준 운동량과 기계적 운동량의 관계를 바꿉니다. 그렇다고 Noether 논증이 사라지는 것은 아닙니다.

정규적인 반지름 영역에서만 역변환하기

각속도에 대한 Hessian은 $\partial^2L/\partial\omega^2=m r^2$입니다. $r>0$에서 양수이므로 운동량 방정식은 유일하게 역변환됩니다.

$$\omega(r,\ell)=\frac{\ell}{m r^2}-\frac{\beta}{2m}.$$

운동량을 고정한 진화에서는 $\dot\omega=-2\ell v_r/(m r^3)$입니다. 정준 운동량이 그대로여도 각속도는 커지거나 작아질 수 있습니다. 전체 각도를 복원하려면 초기 각도와 이 변화율의 적분이 필요합니다. 반지름 곡선을 표시했다고 각도 궤적까지 자동으로 얻은 것은 아닙니다.

축약한 방사 속도의 Hessian은 m입니다. 각도 Hessian은 kg m², 방사 Hessian은 kg으로 단위가 다릅니다. 서로 다른 좌표의 속도에 대응하므로 하나의 수치 크기로 합치거나 좌표와 무관한 스펙트럼처럼 비교하면 안 됩니다.

운동량과 각속도의 곱을 먼저 빼기

이 단원에서는 다음 부호 규약으로 축약 Lagrangian을 정의합니다.

$$\mathcal R(r,v_r;\ell)=\left[L-\ell\omega\right]_{\omega=\omega(r,\ell)}.$$

ell을 고정한 편미분은

$$\partial_r\mathcal R=\partial_rL+ (\partial_\omega L-\ell)\partial_r\omega=\partial_rL$$

입니다. 운동량 방정식 때문에 추가 연쇄 법칙 항이 사라집니다. 반대로 $L_{\mathrm{sub}}=L(r,v_r,\omega(r,\ell))$만 미분하면 $\ell\partial_r\omega$가 남아 방사 힘이 달라집니다. 운동량 역변환의 대수 계산이 틀린 것이 아니라, 변분할 때 고정하는 조건을 바꾼 것이 문제입니다.

일부 문헌은 Routhian의 전체 부호를 반대로 정의하고 제거한 변수에 Hamilton형 부호 규약을 사용합니다. 여기서는 $\mathcal R=L-\ell\omega$를 남은 방사 좌표의 Lagrangian으로 사용합니다. 이름만 비교하지 말고 정의와 운동방정식의 규약을 함께 확인하세요.

보정 퍼텐셜과 올바른 힘 읽기

Routhian은

$$\mathcal R=\frac m2v_r^2-U_{\mathrm{eff}}(r),\qquad U_{\mathrm{eff}}=U+\frac{(\ell-\beta r^2/2)^2}{2mr^2}$$

입니다. 따라서 올바른 방사 방정식은

$$m\ddot r=\frac{\ell^2}{m r^3}-\frac{\beta^2r}{4m}-U'(r)$$

입니다. L에 바로 대입하는 잘못된 절차는

$$m\ddot r\big|_{\mathrm{naive}}=-\frac{\ell^2}{m r^3}-\frac{\beta^2r}{4m}-U'(r)$$

를 예측합니다. 힘의 차이는 $2\ell^2/(mr^3)$ N입니다. ell이 비영이면 원심 항의 부호가 바뀝니다. 이 차이는 ell의 부호에는 의존하지 않지만 복원한 각속도는 일반적으로 부호에 영향을 받습니다.

보정 퍼텐셜을 전개하면 $U+\ell^2/(2mr^2)-\beta\ell/(2m)+\beta^2r^2/(8m)$입니다. 가운데 상수는 고정한 ell에서 방사 힘에는 영향을 주지 않습니다. 하지만 선택한 에너지 규약의 절대값에는 영향을 주므로 복원한 물리 에너지와 축약 에너지를 비교할 때 유지해야 합니다.

전체 복원 항등식 확인하기

각 정규 조회에서 원래 방사 Euler-Lagrange 방정식, 운동량 복원 $p_\theta=\ell$, 변화율 $d p_\theta/dt=0$을 검사합니다. 각가속도의 복원까지 포함합니다. 각속도를 그냥 일정하게 놓으면 충분하지 않습니다.

물리적 에너지는

$$E=\frac m2v_r^2+\frac m2r^2\omega^2+U =\frac m2v_r^2+U_{\mathrm{eff}}$$

입니다. 축약 가속도와 복원한 각가속도를 대입하면 이 에너지의 순간 변화율은0입니다. 기계적 각운동량의 변화율 $-\beta r v_r$는 비영일 수 있습니다. 속도 선형 결합이 운동 성분을 재분배하면서 물리적 에너지에는 기여하지 않는다는 점에서 서로 모순되지 않습니다.

$-\partial_\ell\mathcal R=\omega$는 부분 Legendre 변환의 복원 관계입니다. 이러한 비교는 유한 차분으로 추정한 기울기가 아니라 정확한 유리식 미분으로 검산합니다.

모터가 각속도를 고정하는 경우 비교하기

모터가 $\theta(t)=\Omega t+\theta_0$를 강제할 수 있습니다. 이는 외부에서 준 시간 의존 좌표 조건이며 보존 운동량의 수준집합이 아닙니다. 정해진 상수 Omega를 L에 대입하는 것은 이 다른 방사 문제에는 올바릅니다. 그 힘은

$$F_{\Omega}=(m\Omega^2+\beta\Omega)r-U'(r)$$

입니다. 모터가 공급하는 토크와 일률은

$$\tau_{\mathrm{motor}}=(2m\Omega+\beta)r v_r, \qquad P_{\mathrm{motor}}=\Omega\tau_{\mathrm{motor}}$$

입니다. 물리적 에너지는 $\dot E_\Omega=P_{\mathrm{motor}}$를 만족합니다. 자율적인 방사 보존량 $E_\Omega-\Omega p_\theta$는 일정하지만 물리적 에너지는 일반적으로 달라집니다. 이 방사 보존량을 물리적 에너지와 같다고 부르면 모터의 일을 숨기게 됩니다.

beta 항도 토크 균형에 들어갑니다. Omega가0이면 고정 토크는 비영이어도 일률은0일 수 있습니다. 한 설정에서 일률이0이라는 사실은 각속도를 강제하는 데 외부 상호작용이 전혀 필요 없다는 일반 명제가 아닙니다.

한 점에서 일치하는 오해하기 쉬운 예제

질량1 kg, 강성1 N/m, 사차 계수0, $\ell=1$ Js, beta0, 반지름1 m, 방사 속도1 m/s, 강제 각속도1 s⁻¹로 설정합니다.

이 반지름에서는 복원 각속도와 강제 각속도가 같습니다. 올바른 고정 운동량 힘과 모터 제어 힘은 모두0입니다. 잘못 대입한 힘은−2 N입니다. Routhian은−1/2 J이고 복원한 물리 에너지는3/2 J입니다.

하지만 두 올바른 실험도 같은 모델은 아닙니다. 힘의 반지름 미분은 고정 운동량 문제에서−4 N/m, 강제 각속도 문제에서0입니다. 반지름을 조금 바꾸면 구분됩니다. 또한 모터 제어 에너지 변화율은2 W, 토크는2 N m이지만 닫힌 고정 운동량 문제의 에너지 변화율은0입니다.

축약한 학습 모델을 검증할 때도 한 상태의 일치만 보지 말고 이웃 상태의 응답·미분·고정한 보존 조건까지 비교해야 한다는 예입니다.

정준 운동량0에서도 회전하는 경우

질량1 kg, 퍼텐셜0, $\ell=0$, $\beta=1$ kg/s, 반지름1 m, 방사 속도1 m/s로 둡니다. 복원 각속도는−1/2 s⁻¹입니다. 정준 운동량은0이어도 기계적 각운동량은−1/2 Js입니다.

이 특별한 운동량 영역에서는 Routh 힘과 잘못 대입한 힘이 모두−1/4 N입니다. 이 일치가 다른 운동량에서도 그 절차가 옳다는 뜻은 아닙니다. ell을 비영으로 바꾸면 빠뜨린 원심 항의 부호가 바로 드러납니다.

물리적 에너지는5/8 J, 변화율은0, 기계적 각운동량 변화율은−1 N m입니다. $\beta r^2/2$의 변화가 반대 방향으로 상쇄하므로 정준 운동량은 일정합니다. 정준 운동량0을 각도 정지와 같은 뜻으로 해석하면 안 됩니다.

이 영운동량 예제에서 모터의 강제 각속도를−1/2 s⁻¹로 정하면 실제로 일치하는 특별한 경우도 있습니다. $(2m\Omega+\beta)=0$이므로 필요한 모터 토크와 일률이0이고 두 방사 힘은 모든 정규 반지름에서 같습니다. 이런 특별한 일치가 고정 운동량과 고정 각속도를 일반적으로 서로 바꾸어 사용해도 된다는 뜻은 아닙니다.

반지름0을 명시적으로 미정으로 남기기

반지름0에서는 각도 Hessian이0입니다. 극좌표는 각도를 지정하지 못하고 운동량 방정식은 각속도를 유일하게 결정하지 못합니다. ell이 비영이면 유한한 각속도 해가 없고, ell이0이어도 하나의 각속도를 선택하지 못합니다.

따라서 반지름0에서는 조회 상태와 조회 차원 수를 제공하지 않습니다. 작은 epsilon으로 나누거나 그림의 첫 양의 반지름을 대신 조회하지 않습니다. 곡선은 선언한 정규 영역인0.1 m부터4 m까지만 그립니다.

ell0의 어떤 식이 유한한 극한을 갖더라도 전체 각도 Legendre 사상과 극좌표계가 정규적으로 바뀌었다는 증거는 아닙니다. 한 식의 연속성과 좌표·운동량 관계 전체의 가역성을 구분하세요.

차원을 세고 제거한 정보 복원하기

정규 영역의 원래 위상 공간은 $(r,\theta,p_r,p_\theta)$의4차원입니다. $p_\theta=\ell$을 고정하면3차원 수준집합입니다. 여기서 순환 각도의 대칭 궤도로 몫을 취하면2차원 축약 위상 공간, 즉 하나의 방사 배치 자유도가 남습니다.

운동량 수준집합만으로는 차원이 홀수이며 최종 축약 심플렉틱 위상 공간이 아닙니다. 각도를 하나의 고정 숫자로 제한한 것도 아닙니다. 초기 각도와 각도 진화의 복원 정보는 방사 몫 공간 밖에 있습니다.

다음 단원의 제2종 원형 구속은 위치 조건과 그에 맞는 운동량 조건으로4차원 주변 공간의2차원 면을 고르는 다른 방식입니다. 수행하는 축약을 확인하지 않고 보존량 하나마다 차원을 하나씩 빼기보다, 조건의 수와 몫으로 제거한 대칭 궤도를 따로 세어야 합니다.

반지름 곡선을 시간 운동으로 오해하지 않기

힘·힘의 기울기·운동량·에너지·일률 곡선은 다른 설정을 고정하고 반지름을 바꾼 상태족입니다. 표시한 방사 속도로 입자가 모든 반지름을 지나가는 하나의 해가 아닙니다. 각 상태의 가속도는 해당 상태에서 운동방정식으로 구한 값입니다.

이 실험의 에너지 곡선이 반지름에 따라 달라도 올바르게 복원한 동역학 방향의 에너지 미분은0일 수 있습니다. 상태족이 그 동역학을 따라가지 않으므로 모순이 아닙니다. 모터 일률 곡선도 곧 시간에 대해 적분한 모터의 일을 뜻하지 않습니다.

학습한 축약 모델을 보고할 때는 고정 운동량, 정규 영역, 복원 방정식과 에너지 규약을 기록하세요. ell0이나 한 반지름에서 작은 힘 오차는 잘못된 축약을 숨길 수 있습니다. 몫을 생략한 차원 수나 특이한 식의 유한 극한은 빠진 변수와 잘못된 좌표계를 숨길 수 있습니다.

해설이 있는 연습문제

1. 원심 항의 부호. 첫 예제에서1 m의 올바른 힘과 잘못 대입한 힘을 구하세요.

   해설: 올바른 힘은 $1-1=0$ N, 잘못된 힘은 $-1-1=-2$ N입니다. 차이는 $2\ell^2/(mr^3)=2$ N입니다. 빠진 부분 Legendre 항 때문에 차이가 생깁니다.

2. 이웃 상태 검사. 같은 예제에서 두 올바른 모델의 힘 기울기를 구하세요.

   해설: 고정 ell에서는 $-3\ell^2/(mr^4)-k=-4$ N/m입니다. 고정 Omega에서는 $m\Omega^2-k=0$입니다. 한 반지름에서 힘이 같아도 주변 상태에 대한 응답은 다릅니다.

3. 정준·기계적 운동량. ell0, beta1의 단위 상태에서 각속도·기계적 각운동량·그 변화율을 구하세요.

   해설: 복원 각속도는−1/2 s⁻¹, 기계적 각운동량은−1/2 Js, 변화율은 $-\beta r v_r=-1$ N m입니다. 정준 운동량에서는 $\beta r^2/2$의 변화가 이를 상쇄합니다.

4. 외부의 일. 첫 조회에서 방사 가속도가0이어도 모터 제어 에너지가2 W로 증가하는 이유는 무엇인가요?

   해설: 방사 속도가1 m/s인데 모터가 각속도를 고정하므로 토크2 N m와 일률2 W를 공급합니다. 회전 운동에너지와 퍼텐셜의 증가가 이를 설명합니다. 순간 방사 가속도0이 외부의 일을 없애지는 않습니다.

5. 특이한 역변환. 반지름0, ell0에서 각속도0을 지정하면 축약이 복구되나요?

   해설: 아닙니다. 각도 Hessian이0이므로 운동량 방정식이 각속도를 결정하지 않습니다. 임의로 하나를 고른다고 역변환의 유일성이나 원점의 극각이 생기지 않습니다. 조회는 미정으로 남습니다.

6. 차원 계산. 정규4차원 위상 공간에서 순환 운동량을 고정한3차원 수준집합이 이미 방사 축약 위상 공간인가요?

   해설: 아닙니다. 각도 대칭 궤도로 몫을 취해야2차원이 됩니다. 전체 운동을 복원하려면 초기 각도와 복원 방정식이 필요합니다. 운동량 고정은 각도나 각속도 고정과 다릅니다.

읽을거리와 검증 범위

고전적인 부분 Legendre 축약은 B. Langerock, F. Cantrijn, J. Vankerschaver의 Routhian reduction for quasi-invariant Lagrangians 서론과 정리1에 정리되어 있습니다. David Tong의 Lagrangian Formalism 2.5.3–2.5.4절은 보존 각운동량을 L에 바로 대입하면 힘의 부호가 바뀔 수 있음을 설명합니다. 속도 선형 모델·모터 비교·유리수 예제는 직접 구성했습니다. 일반 비가환·준순환 축약은 추가 가정이 필요하며 여기에서 구현하지 않습니다.