고전 및 라그랑주 역학
Advanced Mechanics
3단계: 비선형 진동·카오스·연속체 역학·다체 시뮬레이션
이론적 기초: 우리가 지금까지 무엇을 가정하고 있었는가
2단계에서 배운 라그랑주 역학과 해밀턴 역학을 잠시 돌이켜 보자. 우리는 일반화 좌표 qᵢ를 도입하고, 라그랑지안 L = T − V를 통해 오일러-라그랑주 방정식을 유도했다. 해밀토니안 H는 위상 공간(phase space)에서 에너지 보존량이었고, 뇨터 정리(Noether's theorem)는 대칭성이 보존량을 낳는다는 아름다운 원리였다. 1단계에서는 스프링-질량 계처럼 단순 조화 진동을 배웠고, 복원력이 F = −kx처럼 변위에 비례하는 경우를 다뤘다.
여기서 질문을 하나 던져보자: F = −kx라는 복원력은 과연 항상 성립하는가? 스프링을 아주 살짝 당기면 훅의 법칙이 꽤 잘 맞는다. 그런데 스프링을 극단적으로 잡아당기면 어떻게 될까? 혹은 진자를 아주 작은 각도가 아니라 90도, 심지어 180도 가까이 당기면? 고등학교 물리 시간에 배운 단진자의 주기 공식 T = 2π√(L/g)는 θ가 작을 때만 성립하는 근사식이다. sin θ ≈ θ라는 가정, 즉 선형화(linearization) 가정 하에서 얻은 결과다. 현실 세계의 대부분은 사실 비선형(nonlinear)이다. 3단계는 바로 그 "근사"의 껍데기를 벗겨내고 현실의 복잡한 얼굴을 마주하는 단계다.
비선형 시스템이 왜 그렇게 특별한지 직관적으로 이해하기 위해 아주 간단한 예를 들어보자. 7살짜리 아이에게 설명한다면 이렇게 말할 수 있다: "나비 한 마리가 브라질에서 날갯짓을 하면 텍사스에서 토네이도가 생길 수 있다." 이것이 **카오스(Chaos)**의 통속적 표현, 이른바 **나비 효과(Butterfly Effect)**다. 에드워드 로렌츠(Edward Lorenz)가 1963년에 기상 예측 모델을 돌리던 중 발견한 것으로, 초기 조건을 0.506127 대신 0.506으로 반올림했더니 시뮬레이션 결과가 완전히 달라졌다. 중학생 수준에서 이를 이해하면 "같은 방정식인데 시작점이 아주 조금 달라도 결과가 완전히 달라지는 것"이고, 고등학생-대학생 수준에서는 **초기 조건에 대한 민감 의존성(sensitive dependence on initial conditions)**으로 정의된다. 그리고 현장 연구자 수준에서는 이를 **양의 리아푸노프 지수(positive Lyapunov exponent)**로 정량화한다.
[노트 기록] 선형 시스템 vs 비선형 시스템: 선형 시스템은 중첩 원리(superposition)가 성립한다. 즉 두 해의 합도 해다. 비선형 시스템은 중첩 원리가 깨진다. 방정식에 x², sin x, x·ẋ 등 변위나 속도의 곱 또는 비선형 함수가 포함되면 비선형이다.
연속체 역학과 탄성 이론은 또 다른 방향으로의 확장이다. 1단계와 2단계에서 우리는 줄곧 질점(point mass) 또는 **강체(rigid body)**를 다뤘다. 질점은 크기가 없고, 강체는 변형이 없다. 하지만 현실에서 금속 빔은 하중을 받으면 휘고, 고무공은 충돌하면 찌그러진다. 건물은 지진파에 흔들리고, 뼈는 충격을 받으면 금이 간다. 이처럼 **변형 가능한 물체(deformable body)**를 다루는 것이 연속체 역학이다. 마지막으로 다체 문제(N-body problem)는 1단계의 뉴턴 역학에서 출발하지만, N ≥ 3이 되는 순간 해석적 해가 존재하지 않는다는 놀라운 사실과 마주한다. 이는 앙리 푸앵카레(Henri Poincaré)가 1889년에 수학적으로 증명했다. 이 세 가지 주제—카오스, 연속체 역학, 수치 다체 시뮬레이션—은 표면적으로 달라 보이지만, 모두 "실제 세계는 단순한 방정식으로 완전히 기술되지 않는다"는 하나의 진실을 향해 수렴한다.
1. 비선형 진동과 카오스 이론
선형 진동자에서 비선형 진동자로
2단계에서 단순 조화 진동자의 라그랑지안은 L = ½mẋ² − ½kx²였다. 이로부터 유도된 운동 방정식은 ẍ + ω₀²x = 0이고, 해는 x(t) = A cos(ω₀t + φ)로 깔끔하게 떨어진다. 진폭 A가 달라져도 진동 주기는 ω₀에 의해 고정된다. 이것이 선형 진동자의 핵심 특성이다. 그런데 실제 스프링은 너무 많이 당기면 "빳빳해지는" 경향이 있다. 이를 경화(hardening) 스프링이라 하고, 반대로 더 늘어지는 것을 연화(softening) 스프링이라 한다. 이를 수식으로 표현하면 복원력은 F = −kx ± βx³으로 나타난다. 이로부터 **더핑 방정식(Duffing equation)**이 탄생한다.
더핑 방정식: ẍ + δẋ + αx + βx³ = γ cos(ωt)
여기서 δ는 감쇠 계수, α와 β는 선형 및 비선형 복원력 계수, γ cos(ωt)는 외부 주기 구동력이다. δ = β = γ = 0이면 단순 조화 진동자로 환원된다. α < 0, β > 0이면 이중 우물 포텐셜(double-well potential) 형태가 되는데, 포텐셜 V(x) = −½αx² + ¼βx⁴는 x = 0에서 불안정 평형점을 가지고, x = ±√(−α/β)에서 두 개의 안정 평형점을 가진다. 여기에 주기적 구동력을 가하면 어떻게 될까? 스스로 잠깐 생각해보라: 안정 평형점이 두 개 있는 계에 주기적으로 흔들리는 힘을 가하면 입자는 어느 쪽 우물에 머무를까?
[노트 기록] 더핑 방정식: ẍ + δẋ + αx + βx³ = γ cos(ωt). 매개변수 α < 0, β > 0일 때 이중 우물 포텐셜 형성. 비선형 항 βx³가 카오스 거동의 핵심이다.
위상 공간에서 카오스를 '보다'
2단계에서 위상 공간(phase space)의 개념을 배웠다. 상태 (x, ẋ)로 이루어진 2차원 공간에서 시스템의 궤적을 그리면 선형 조화 진동자는 닫힌 타원을 그린다. 감쇠 진동자는 나선형으로 수렴해서 원점으로 빨려들어간다. 이 수렴하는 점을 **점 어트랙터(point attractor)**라 한다. 만약 수렴하는 목적지가 점이 아니라 닫힌 곡선이라면 **한계 사이클(limit cycle)**이라 한다. 이제 더핑 진동자처럼 외부에서 주기적인 힘이 가해지는 계에서 카오스가 발생하면 어트랙터는 어떤 모양이 될까? 놀랍게도 점도 아니고, 한계 사이클도 아닌, **분수 차원(fractal dimension)**을 가지는 기묘한 구조가 나타난다. 이것이 **이상한 어트랙터(strange attractor)**다. 2차원 공간에 갇혀 있는 것처럼 보이지만 무한히 많은 층으로 이루어진 구조다. 가장 유명한 것이 로렌츠 어트랙터(Lorenz attractor)인데, 나비 모양의 3차원 구조를 가진다. 로렌츠 방정식은 다음과 같다.
ẋ = σ(y − x), ẏ = x(ρ − z) − y, ż = xy − βz
(σ, ρ, β는 매개변수로, 로렌츠가 기상 대류를 모델링할 때 쓴 계수들이다. 표준 카오스 거동은 σ = 10, ρ = 28, β = 8/3에서 나타난다.)
위상 공간에서 카오스 궤적은 이상한 어트랙터 위를 영원히 돌아다니지만, 절대로 같은 점을 두 번 통과하지 않는다. 아주 가깝게 시작한 두 궤적은 지수적으로 빠르게 벌어진다. 바로 이 지수적 발산의 속도를 **리아푸노프 지수(Lyapunov exponent)**로 정량화한다.
리아푸노프 지수: 카오스의 정량적 지표
두 초기 조건이 아주 작은 거리 δ₀만큼 떨어져 시작했다고 하자. 시간이 지남에 따라 그 거리의 변화를 살펴보면 δ(t) ≈ δ₀ · e^(λt) 로 표현된다. 여기서 λ가 리아푸노프 지수다. λ > 0이면 두 궤적이 지수적으로 발산하며 카오스의 필요 조건이다. λ < 0이면 수렴, λ = 0이면 주기적 또는 준주기적 거동을 의미한다. 실제로 리아푸노프 지수를 계산하는 방법은 수치적으로 이루어진다. 초기 조건 x₀를 잡고, 아주 살짝 다른 x₀ + ε₀에서 시뮬레이션을 동시에 돌린다. 시간이 지나면서 두 궤적 사이의 거리 ε(t)를 추적하되, 거리가 너무 커지면 재규격화(renormalization)를 해준다. 이 과정을 반복하면 평균 발산율 λ를 계산할 수 있다. λ의 단위는 [1/시간]이며, 카오스 계의 예측 가능 시간은 대략 1/λ 정도다.
[노트 기록] 리아푸노프 지수 λ: δ(t) ≈ δ₀ · e^(λt). λ > 0 → 카오스. 예측 가능 지평(prediction horizon) ≈ 1/λ. 대기 예측의 경우 λ⁻¹ ≈ 5일, 즉 일기예보가 일주일 이상 되기 힘든 물리적 근거다.
비프르케이션과 카오스로의 경로
카오스가 갑자기 터지는 것은 아니다. 하나의 매개변수(예: 구동력의 진폭 γ)를 서서히 증가시키면 계는 단계적으로 변한다. 처음에는 단순한 주기 운동을 하다가, 매개변수가 임계값을 넘으면 갑자기 주기가 두 배로 늘어난다(주기 배가 분기, period-doubling bifurcation). 더 증가시키면 또 두 배, 그리고 또 두 배... 이것이 페이겐바움(Feigenbaum) 수열이다. 미첼 페이겐바움(Mitchell Feigenbaum)은 1975년에 연속적인 분기점 사이의 비율이 δ ≈ 4.6692...라는 보편 상수로 수렴함을 발견했다. 이 상수는 계의 세부 내용에 무관하게 나타나는 **보편성(universality)**의 한 예다. 생각해보면 놀랍다: 더핑 진동자든, 로지스틱 사상(logistic map)이든, 특정 유형의 비선형 계는 모두 같은 수 4.6692로 카오스로 진입한다.
2. 연속체 역학과 탄성 이론
질점에서 연속체로
1단계에서 강체를 다룰 때 우리는 물체를 무수한 질점의 집합으로 보고, 그 집합의 회전 관성 모멘트를 I = ∫r²dm으로 계산했다. 연속체 역학은 이 아이디어를 극단까지 밀고 간다. 물체는 무한히 많은 무한소 질점들의 연속체이며, 각 점은 변형에 의해 자신의 위치가 변할 수 있다. 핵심 아이디어는 이렇다: 어떤 물체의 작은 조각에 힘이 가해지면 그 조각은 변형된다. 이 변형을 두 가지 측면에서 기술해야 한다. 하나는 물체 내부에서 힘이 어떻게 전달되는가(응력, stress)이고, 다른 하나는 물체가 실제로 얼마나 변형되었는가(변형률, strain)이다.
응력 텐서
벽돌 위에 책을 올려놓는다고 상상하자. 책의 무게는 벽돌의 윗면을 수직으로 누른다. 이 힘을 단위 면적으로 나눈 것이 **응력(stress)**이다. 수직 방향의 응력을 수직 응력(normal stress), 면에 평행한 방향의 응력을 **전단 응력(shear stress)**이라 한다. 이제 3차원으로 일반화하면, 물체 내부의 임의의 점에서 x, y, z 각각의 면을 잘라낼 수 있고, 각 면에는 3개 방향의 응력 성분이 존재한다. 총 9개의 응력 성분이 나오며, 이것이 응력 텐서 σᵢⱼ를 이룬다. σᵢⱼ는 j방향 면에 i방향으로 작용하는 응력을 의미한다. 각운동량 보존에 의해 σᵢⱼ = σⱼᵢ가 성립하므로 (2단계에서 배운 뇨터 정리와 연결해 생각해보라), 독립적인 성분은 사실 6개다.
[노트 기록] 응력 텐서 σᵢⱼ: 3×3 대칭 행렬. σᵢᵢ = 수직 응력 (i = j), σᵢⱼ (i ≠ j) = 전단 응력. 단위: Pa (Pascal = N/m²). 대칭성 σᵢⱼ = σⱼᵢ에 의해 독립 성분 6개.
변형률 텐서
응력이 힘의 관점이라면, 변형률은 변형의 관점이다. 물체 내 한 점의 변위를 u = (u₁, u₂, u₃)라 할 때, 선형 변형률 텐서는 다음과 같이 정의된다.
εᵢⱼ = ½(∂uᵢ/∂xⱼ + ∂uⱼ/∂xᵢ)
이것도 대칭 텐서다. 대각 성분 εᵢᵢ는 각 방향으로의 단순 신장률, 비대각 성분 εᵢⱼ (i ≠ j)는 전단 변형에 해당한다. 이 정의가 왜 저렇게 생겼는지 스스로 생각해보라: 만약 물체 전체가 강체 회전을 하면 변형이 없어야 한다. 그런데 단순히 ∂uᵢ/∂xⱼ만 쓰면 강체 회전도 잡아버린다. 대칭화한 형태 ½(∂uᵢ/∂xⱼ + ∂uⱼ/∂xᵢ)가 바로 순수한 변형만을 포착한다.
일반화된 훅의 법칙과 구성 방정식
1단계에서 배운 스프링의 훅의 법칙 F = −kx는 사실 이것의 1차원 특수 경우다. 3차원에서 일반적인 선형 탄성 물체는 σᵢⱼ = Cᵢⱼₖₗ εₖₗ 로 기술된다. 여기서 Cᵢⱼₖₗ은 **탄성 텐서(elasticity tensor)**로 4차 텐서다. 원리적으로 3⁴ = 81개의 성분이 있지만, 대칭성에 의해 등방성(isotropic) 재료에서는 단 두 개의 독립 매개변수로 완전히 기술된다. 이 두 매개변수가 라메 상수(Lamé constants) λ, μ이고, 이들로부터 우리에게 친숙한 영률(Young's modulus) E = μ(3λ + 2μ)/(λ + μ) (재료가 당겨지거나 압축될 때 저항하는 정도)와 푸아송 비(Poisson's ratio) ν = λ/[2(λ + μ)] (재료를 한 방향으로 당길 때 수직 방향으로 수축하는 비율)가 유도된다. 철의 경우 E ≈ 200 GPa, ν ≈ 0.3이다. 고무는 E ≈ 0.01~0.1 GPa, ν ≈ 0.5에 가깝다(거의 비압축성).
[노트 기록] 등방성 탄성체 구성 방정식: σᵢⱼ = λδᵢⱼεₖₖ + 2μεᵢⱼ (아인슈타인 합산 규약). δᵢⱼ는 크로네커 델타. E = μ(3λ+2μ)/(λ+μ), ν = λ/[2(λ+μ)]. 역으로 λ = Eν/[(1+ν)(1−2ν)], μ = E/[2(1+ν)].
**나비에 방정식(Navier's equations)**은 이 구성 방정식을 뉴턴의 운동 방정식에 대입해 얻은 변위 u에 대한 편미분 방정식이다. 정적 평형 상태(가속도 = 0)에서 (λ + μ)∇(∇·u) + μ∇²u + f = 0 이고, 여기서 f는 물체력(중력 등)이다. 이 방정식의 해가 각 점의 변위 u를 알려주고, 그로부터 변형률→응력 순으로 계산할 수 있다. 이것이 구조 해석(structural analysis)과 유한 요소법(Finite Element Method, FEM)의 이론적 토대다.
3. 다체 문제와 수치 적분 기법
두 물체에서 세 물체로: 해가 사라지다
1단계에서 뉴턴의 중력 법칙 F = −GMm/r²를 배웠고, 케플러 문제(Kepler problem)는 2체 문제를 환원 질량으로 1체 문제로 바꿔 해석적으로 풀 수 있다. 해석적 해가 존재한다는 것은 모든 초기 조건에 대해 임의의 시간 t에서의 위치를 공식으로 계산할 수 있다는 뜻이다. N = 3이 되면 어떻게 될까? 푸앵카레는 1889년에 이 문제에 매달렸다. 그는 원래 스웨덴 국왕이 제시한 상금 문제 ("태양계는 안정한가?")에 답하려 했고, 처음에는 안정하다는 증명을 제출했다. 그러나 검토 과정에서 자신의 증명에 오류가 있음을 발견했고, 이를 수정하는 과정에서 사실상 3체 문제의 카오스적 거동을 발견했다. 3체 문제에는 해석적 일반해가 없다. 특수한 초기 조건(라그랑주 해, 8자형 주기 궤도 해 등)을 제외하면, N ≥ 3 중력 다체 계는 원칙적으로 카오스적이다. 우리는 수치 적분에 의존해야 한다.
수치 적분: 미분 방정식을 컴퓨터로 푸는 법
N개 물체에 대해 운동 방정식은 다음과 같다.
drᵢ/dt = vᵢ, dvᵢ/dt = (1/mᵢ) Σⱼ≠ᵢ Gmᵢmⱼ(rⱼ − rᵢ)/|rⱼ − rᵢ|³
이를 수치적으로 풀기 위해 시간을 작은 간격 Δt로 나눈다. **오일러 방법(Euler Method)**은 가장 단순한 1차 방법으로, x(t+Δt) ≈ x(t) + v(t)·Δt, v(t+Δt) ≈ v(t) + a(t)·Δt 로 다음 상태를 직선 근사로 예측한다. 오차가 O(Δt)로 축적되며, 역학 시스템에서 에너지를 보존하지 못한다는 치명적 단점이 있어 행성 궤도 시뮬레이션에 쓰면 궤도가 나선형으로 발산하거나 수렴한다.
**4차 룽게-쿠타(RK4)**는 현재 상태에서 Δt 간격 동안 4번의 기울기를 계산하고 이를 가중 평균해 다음 상태를 예측하는 4차 정확도(오차 O(Δt⁴)) 방법이다. k₁ = f(tₙ, yₙ), k₂ = f(tₙ + Δt/2, yₙ + Δt·k₁/2), k₃ = f(tₙ + Δt/2, yₙ + Δt·k₂/2), k₄ = f(tₙ + Δt, yₙ + Δt·k₃)라 하면 yₙ₊₁ = yₙ + (Δt/6)(k₁ + 2k₂ + 2k₃ + k₄)다. 훨씬 정확하지만 RK4도 일반적으로 에너지를 완벽히 보존하지 못한다.
**심플렉틱 적분기(Symplectic Integrator)**는 2단계에서 배운 해밀턴 역학의 구조, 즉 위상 공간의 **심플렉틱 구조(symplectic structure)**를 보존하도록 설계된 방법이다. 가장 간단한 예인 속도 베를레(Velocity Verlet) 알고리즘은 다음과 같다: x(t+Δt) = x(t) + v(t)·Δt + ½a(t)·Δt², 이어서 a(t+Δt) = F[x(t+Δt)]/m를 새로 계산하고, v(t+Δt) = v(t) + ½[a(t) + a(t+Δt)]·Δt로 속도를 업데이트한다. 이 방법은 시간 반전 대칭성을 가지며 에너지를 평균적으로 보존한다. 분자 동역학(MD), 천체역학 시뮬레이션 모두 이 계열의 심플렉틱 적분기를 표준으로 사용한다.
[노트 기록] 수치 적분기 비교: 오일러는 1차 정확도, O(Δt) 오차, 에너지 비보존(부적합). RK4는 4차 정확도, O(Δt⁴) 오차, 단기 시뮬레이션 적합. 속도 Verlet은 2차 정확도이지만 시간 반전 대칭, 장기 에너지 보존 → 천체역학 표준.
실제 태양계 시뮬레이션에서 소행성대, 카이퍼 벨트까지 포함하면 N이 수백만에 달한다. 연산 복잡도는 O(N²)이므로 이를 효율적으로 처리하기 위해 반스-헛(Barnes-Hut) 알고리즘(트리 구조를 이용해 O(N log N)으로 줄임)이나 입자-메쉬(Particle-Mesh) 방법(푸리에 변환 활용)이 사용된다. 현재 우주 구조 형성 시뮬레이션(Illustris TNG 등)은 수십억 개의 입자를 다룬다.
프로젝트: 예제 문제 모음
(문제만 제시한다. 풀이와 정답은 없다. 각 문제를 스스로 사고하며 접근하라. 총 풀이 시간 약 40분으로 설계했다.)
프로젝트 A: 더핑 진동자와 카오스 분석 (약 15분)
[문제 A-1] 더핑 방정식 ẍ + 0.3ẋ − x + x³ = 0.5 cos(1.2t) 가 있다. 이 계의 포텐셜 에너지 V(x) = −½x² + ¼x⁴를 스케치하라. 평형점은 어디에 몇 개 있으며, 각각 안정인가 불안정인가? 각 평형점이 안정/불안정인 이유를 V''(x)의 부호로 설명하라. (단서: V''(x) > 0이면 포텐셜이 위로 오목, 즉 복원력이 평형점으로 향한다.)
[문제 A-2] 위의 더핑 계에서 초기 조건 (x₀, ẋ₀) = (1.0, 0.0)과 (1.001, 0.0)으로 각각 시뮬레이션을 100초 동안 돌린다고 하자. 어떤 수치 방법을 쓸지 선택하고 그 이유를 써라. 두 궤적 사이의 거리 |Δx(t)|를 반로그 그래프(y축이 log 스케일)에 그리면 어떤 형태가 나타날 것으로 예상하는가? 이 그래프의 기울기와 리아푸노프 지수 λ는 어떤 관계인가? 수식으로 도출하라.
[문제 A-3] 구동 진폭 γ를 0.1에서 0.5까지 0.05 간격으로 변화시키면서, 각 γ값에서 500번의 진동 이후(과도 응답이 사라진 후) x(t)를 기록하라. 구동력의 주기마다 x값을 점으로 찍어 **푸앵카레 단면(Poincaré section)**을 γ에 대한 함수로 그려라. γ가 증가함에 따라 점의 개수가 어떻게 변하는가? 주기 배가 분기(period-doubling bifurcation)가 일어나는 γ의 임계값을 대략적으로 찾아라. 이 전체 그림을 **분기 다이어그램(bifurcation diagram)**이라 한다.
프로젝트 B: 탄성 이론 기초 계산 (약 10분)
[문제 B-1] 단면적 A = 1 cm², 길이 L₀ = 1 m인 강철 봉(E = 200 GPa, ν = 0.3)에 인장력 F = 10 kN이 가해진다. (a) 수직 응력 σ₁₁을 계산하라. (b) 수직 변형률 ε₁₁ = σ₁₁/E를 계산하라. (c) 봉의 길이 변화 ΔL을 계산하라. (d) 푸아송 비를 이용해 봉의 단면 방향(ε₂₂ = ε₃₃ = −ν·ε₁₁)의 변형률을 계산하고, 단면 반지름(r₀ = 5.64 mm)의 변화량을 구하라.
[문제 B-2] 재료 A(E_A = 70 GPa, 알루미늄)와 재료 B(E_B = 5 GPa, 나일론)가 직렬로 연결되어 있다(같은 힘이 두 재료 모두를 통해 전달됨). 각각 길이 0.5 m, 단면적 동일. 총 신장량을 계산하고, 각 재료의 기여도 비율(%)을 구하라. 결과로부터 "직렬 연결에서 변형의 대부분은 어떤 특성의 재료에서 발생하는가?"에 대한 일반 원리를 언어로 서술하라.
프로젝트 C: 3체 중력 문제 수치 시뮬레이션 설계 (약 15분)
[문제 C-1] 동일한 질량 m = 1, G = 1인 단위 계에서 세 물체의 초기 조건이 다음과 같다. 물체 1: r₁ = (−1, 0), v₁ = (v_x, v_y). 물체 2: r₂ = (1, 0), v₂ = (−v_x, v_y). 물체 3: r₃ = (0, 0), v₃ = (0, −2v_y). 이 계가 질량중심 기준에서 전체 운동량이 0이 되려면 v_x와 v_y 사이에 어떤 관계가 있어야 하는가? 총 선운동량 P = Σ mᵢvᵢ = 0 조건으로부터 유도하라.
[문제 C-2] 위 계에 속도 베를레(Velocity Verlet) 알고리즘을 구현해 시뮬레이션하는 의사코드(pseudocode)를 작성하라. 각 시간 스텝에서 수행해야 할 연산 순서(위치 업데이트 → 힘 재계산 → 속도 업데이트)를 명시하라. 또한 시뮬레이션의 에너지 보존 여부를 검증하기 위해 어떤 물리량을 매 스텝 계산해 기록해야 하는지 서술하라.
[문제 C-3] 세 물체 중 하나의 초기 위치에 δ = 10⁻⁶의 교란을 주어 두 개의 시뮬레이션을 동시에 돌린다. t = 0, 5, 10, 20, 50 시간 단위에서 두 시뮬레이션 사이의 위치 차이 |Δr(t)|를 계산하라. 이 값들을 반로그 그래프에 찍어 리아푸노프 지수 λ를 추정하라. 만약 λ > 0이면 이 계가 카오스적이라고 결론 내릴 수 있는가? 그 판단의 한계(단순히 λ > 0만으로 카오스라 단정하기 어려운 이유)는 무엇인가?
평가 기준
이번 3단계 최종 결과물[연구: 태양계 N-body 시뮬레이터]은 다음 기준으로 평가한다. 시뮬레이션 정확도(40점): 에너지와 각운동량 보존이 장기 시뮬레이션(t = 1000 단위 이상)에서 허용 오차(0.1%) 이내로 유지되는가, 속도 Verlet과 RK4를 비교했는가. 카오스 분석(40점): 리아푸노프 지수를 실제로 수치 계산했는가, 분기 다이어그램을 구현했는가, 결과를 문헌값과 비교했는가. 연구 리포트(20점): 수식 유도의 논리적 완결성, 시각화의 명료성, 결과의 물리적 해석의 깊이.
이것이 3단계의 전체 지형이다. 카오스 이론에서는 결정론적 방정식에서 예측 불가능성이 탄생하는 역설을, 연속체 역학에서는 텐서로 기술되는 변형의 세계를, 수치 역학에서는 컴퓨터가 어떻게 자연을 흉내 내는가를 배웠다. 이 세 기둥은 현대 공학과 물리학 전반에 뿌리를 내리고 있다—기상 예측, 재료 설계, 우주 탐사, 분자 동역학, 유체 시뮬레이션까지. 프로젝트를 진행하면서 막히는 부분이 생기면 다시 이 글로 돌아와서 해당 절을 반복해 읽어보라. 이해는 한 번에 오지 않는다. 같은 개념을 세 번 다른 각도에서 보면 어느 순간 전체 그림이 뚜렷하게 보이기 시작한다.