케플러 시간과 타원의 면적 법칙

타원에는 시계가 필요하다

궤도 방정식은 가능한 위치를 알려 주지만 각 위치에 언제 도달하는지는 알려 주지 않습니다. 타원 위의 점을 그림의 각도에 따라 일정하게 움직이면 속도·각운동량·힘이 잘못될 수 있습니다. 이 단원에서는 케플러 방정식을 풀어 역제곱 인력의 타원에 실제 시간을 부여합니다.

위성이 원천 가까이를 빠르게 지나고 멀리에서는 천천히 움직이는 모습을 떠올리세요. 같은 시간에는 힘의 중심에 대해 같은 면적을 쓸고 지나갑니다. 같은 거리를 가거나 타원의 기하학적 중심에 대해 같은 각도만큼 움직인다는 뜻이 아닙니다. 위치와 속도를 함께 계산해 이 차이를 검사합니다.

충돌하지 않는 결합 궤도 계열을 고른다

고정 원천 모델 q̈=−μq/r³에서 μ와 m은 양수입니다. 긴반지름 a>0과 이심률 0≤e≤0.95를 고릅니다. 초기 위치는 초점에서 가장 가까운 근점입니다. e=0이면 모든 점이 같은 거리에 있으므로 초기점은 특별한 근점이 아니라 선택한 위상입니다.

$$r_p=a(1-e)>0,\qquad r_a=a(1+e),\qquad b=a\sqrt{1-e^2}.$$

이 계열은 특이한 원점에 도달하지 않습니다. 이심률 상한은 모든 제어 조합을 엄격한 타원 범위에 두기 위한 교육용 범위입니다. 더 큰 이심률을 임의로 0.95로 바꾸어 계산하는 근사가 아닙니다. 포물선 비행·쌍곡선 비행·방사 충돌 시간은 다른 매개화와 정의역 처리가 필요하며 여기서 계산하지 않습니다.

고정 방향각은 실험실에서 타원을 회전시킵니다. 순행·역행 선택은 회전 방향을 바꿉니다. 두 조작 모두 a, e, 주기와 단위 질량당 에너지를 바꾸지 않습니다.

초점과 타원의 중심을 구분한다

근점 방향을 +x로 두면 이심근점이각 E로 위치를 나타낼 수 있습니다.

$$x=a(\cos E-e),\qquad y=b\sin E.$$

힘의 중심은 초점 (0,0)이고 타원의 기하학적 중심은 (−ae,0)입니다. E는 보조 원으로 정의하는 매개변수이며 초점에서 본 실제 극각과 다릅니다. x좌표의 −ae 이동을 빼먹으면 그림이 타원처럼 보여도 힘 중심과의 관계가 달라집니다.

방향각을 바꿀 때는 위치와 속도에 같은 회전행렬을 적용해야 합니다. 위치만 회전시키고 속도를 그대로 두면 접선 조건과 각운동량이 바뀝니다. 좌표 변환과 새로운 물리적 초기 조건을 구분하세요.

면적 법칙으로 시간 방정식을 유도한다

양의 평균 운동 n과 주기 T를 다음과 같이 정의합니다.

$$n=\sqrt{\frac\mu{a^3}},\qquad T=\frac{2\pi}{n}.$$

순행은 s=+1, 역행은 s=−1로 두고 t=0에서 근점을 출발하면 평균근점이각은 M=snt입니다. 케플러 방정식은

$$E-e\sin E=M$$

입니다. 왜 사인 보정이 들어갈까요? q를 E로 미분하면 q×dq/dE=ab(1−e cos E)입니다. E=0부터 적분한 초점 기준 부호 면적은 A(E)=ab(E−e sin E)/2입니다. 면적 속도를 일정하게 요구하면 위 시간 관계가 나옵니다.

따라서 e>0일 때 E가 시간에 비례하도록 움직이면 안 됩니다. 1−e cos E가 일정하지 않기 때문입니다. 평균근점이각은 시간에 비례하지만 이심근점이각은 일반적으로 그렇지 않습니다.

단조 방정식을 안전하게 푼다

좌변의 도함수는 1−e cos E≥1−e>0입니다. 따라서 현재 타원 범위에서는 유한한 M마다 해가 하나 존재합니다. |e sin E|≤e이므로 해는 반드시 [M−e,M+e]에 있습니다.

수치해법은 이 구간을 유지하면서 Newton 단계를 시도합니다. 제안한 값이 구간을 벗어나면 이분법을 사용하고 잔차의 부호로 적절한 끝점을 갱신합니다. 계산 결과에는 E−e sin E−M의 잔차와 반복 횟수를 함께 둡니다.

작은 잔차는 이 음함수 방정식을 잘 풀었다는 증거입니다. 위치의 초점 이동이나 속도 미분, 시간 척도까지 모두 맞았다는 증거는 아닙니다. 그래서 별도의 고정밀 근과 데카르트 운동방정식 적분으로 위치·속도까지 대조합니다.

e=0이면 E=M을 직접 반환하며 t=0도 정확하게 처리합니다. E와 M은 주값 각도로 줄이지 않고 누적 회전수를 유지합니다. 현재 위치만 같아도 몇 바퀴 돌았는지에 따라 누적 면적은 달라집니다.

실제 시간으로 미분해 속도를 구한다

음함수 미분으로

$$\dot E=\frac{s n}{1-e\cos E}$$

를 얻습니다. 따라서 위치와 호환되는 속도는

$$v_x=-a\sin E\,\dot E,\qquad v_y=b\cos E\,\dot E$$

입니다. 접선처럼 보이는 화살표를 따로 정한 것이 아니라 표시 위치의 실제 시간 미분입니다. 근점에서는 분모가 1−e로 작아 빠르게 움직이고, 먼점에서는 1+e로 커져 느려집니다.

역행에서는 E의 진행 부호가 반대입니다. 방향각 회전을 적용하기 전에는 같은 양의 시각에서 순행 위치를 근점 축에 대해 반사한 점이 됩니다. 반지름은 r=a(1−e cos E)>0이며 가속도와 힘은 초점에 대한 −μq/r³, −mμq/r³로 계산합니다.

에너지·각운동량·속력을 각각 검사한다

이 계열에서는 모든 시각에

$$\epsilon=-\frac\mu{2a},\qquad h=s\sqrt{\mu a(1-e^2)},\qquad L=mh$$

입니다. 순간 운동에너지는 K=m|v|²/2, 퍼텐셜은 U=−mμ/r이며 합은 mε와 같아야 합니다. 속력은 활력 방정식도 만족합니다.

$$|v|^2=\mu\left(\frac2r-\frac1a\right).$$

이심률이 있는 궤도에서는 K와 U가 변하지만 합은 일정합니다. 전체 에너지만 그리면 이 교환을 놓치기 쉽습니다. 반대로 K가 변한다는 이유만으로 마찰 손실이 있다고 판단하면 안 됩니다. 이 모델에는 항력이나 열 항이 없습니다.

면적에는 방향과 누적 횟수가 있다

초점 기준 면적 속도와 누적 면적은

$$\dot A=\frac{q\times v}{2}=\frac h2,\qquad A(t)=\frac{ht}{2}$$

입니다. 앞서 얻은 A(E)=ab(E−e sin E)/2를 별도로 계산하면 시간 풀이와 기하의 일관성을 확인할 수 있습니다. 순행 한 바퀴는 πab, 역행 한 바퀴는 −πab를 더합니다. 여러 바퀴 돌면 계속 누적됩니다.

E를 먼저 2π로 나눈 나머지로 바꾸면 전체 회전의 면적을 잃습니다. 여기서 면적은 경로를 따라 적분한 부호 있는 q×dq/2입니다. 화면에서 방문한 픽셀 영역의 합집합 넓이가 아니므로 타원의 넓이를 여러 배 넘을 수 있습니다.

전체 타원 기하와 실제 시간 경로를 따로 읽는다

닫힌 전체 타원은 한 공전의 기준 기하입니다. 실제 경로는 0초부터 조회 시각까지만 표시하며, 조회 시각 0에서는 초기점만 있습니다. 조회점은 직접 푼 상태를 뜻합니다. 전체 기준 타원이 그려졌다는 이유로 이미 모든 점을 방문했다고 해석하지 않습니다.

시간 패널은 0~20초를 다룹니다. 반지름·속력·에너지·면적은 단위가 달라 별도 축이나 패널이 필요합니다. 표본에는 조회 시각과 구간 안의 근점·먼점을 포함합니다. E를 조밀하게 나누고 각 점을 실제 케플러 시간으로 바꾸면 빠르고 짧은 근점 통과에 시간 표본이 더 많이 모입니다. 높은 이심률에서 성긴 등시간 표본만 쓰면 이 특징을 놓칠 수 있습니다.

표본화는 그림의 해상도 선택이며 운동 방정식을 대신하는 근사가 아닙니다. 조회 상태는 음함수 방정식을 직접 풀기 때문에 가장 가까운 미리 계산한 프레임을 고르는 방식과 다릅니다.

예제와 해설이 있는 연습문제

1. 원궤도 극한. μ=a=1, e=0, 순행, t=1로 두세요. 해설: n=1, E=M=1이고 q=(cos 1,sin 1), 속력과 반지름은 1, 주기는 2π입니다. 면적은 1/2입니다. e=0이므로 이 경우에만 타원 중심과 초점이 일치합니다.

2. 이심률이 있는 초기 상태. μ=2, a=2, e=1/2, m=3, t=0입니다. 해설: 근점은 1, 먼점은 3, 짧은반지름은 √3이고 초기 속력은 √3입니다. 가속도는 (−2,0), 전체 에너지는 −3/2, L은 3√3입니다. 타원 중심 (−1,0)과 힘 중심 (0,0)을 구분하세요.

3. 주기 척도. μ와 e를 고정하고 a를 두 배로 하면 어떨까요? 해설: 주기는 2의 3/2제곱 배가 됩니다. 반대로 a를 고정하고 μ를 네 배로 하면 주기는 절반입니다. T²/a³=4π²/μ이며 같은 a, μ에서 e는 주기에 들어가지 않습니다.

4. 같은 면적과 다른 각도. 이심률이 큰 궤도의 근점과 먼점 부근에서 같은 짧은 시간을 비교하세요. 해설: 둘 다 hΔt/2의 면적을 쓸지만 근점 쪽에서 더 큰 각도를 더 빠르게 통과합니다. 극각을 일정 속도로 바꾸는 애니메이션은 이 법칙을 어깁니다.

5. 역행과 회전. 방향을 반전한 뒤 방향각을 π/2 바꾸세요. 해설: 역행은 h, L, 면적의 부호를 반전하고 같은 양의 시각의 에너지와 반지름은 유지합니다. 회전은 (x,y)를 (−y,x)로 보내며 속도도 같이 변환합니다. 주기는 유지됩니다.

6. 여러 바퀴의 기록. μ=6, a=1/2, e=0.95, t=20을 사용하세요. 해설: n=√48이므로 22바퀴 이상 돕니다. 순행 E는 40π보다 커야 하며 면적도 타원 한 개의 넓이를 여러 배 넘습니다. 주값 각도만 사용하면 현재 점은 비슷해 보여도 회전 이력과 총 면적이 잘못됩니다.

연구에서 무엇을 검증할 것인가

학습한 동역학을 거의 원인 그림에서만 검사하지 말고 초점, 높은 이심률, 긴 시간의 위상 누적에서 확인해야 합니다. 좌표 정확도·에너지 보존·시간 정확도·면적 누적은 서로 다른 요구입니다. 불변량 잔차가 작아도 위상 오차가 클 수 있으므로 독립적인 위치와 속도 참조값이 필요합니다.

계산 모듈은 70자리 이분법 근과 별도의 데카르트 역제곱 운동 적분을 참조값으로 사용합니다. 부동소수점 잔차를 정확한 대수 항등식으로 부르지 않습니다. 물리 가정은 고정 원천의 Newton 인력과 충돌하지 않는 타원 계열입니다. 유도의 원자료는 MIT OpenCourseWare 케플러 문제 강의 노트를 참고하세요.