이중 평균법

수치해석 계산물리 물리 마지막 수정: 2026-08-15 04:39:52

1. 개요[편집]

빠른 진동은 평균 내서 지우면 된다. 문제는 지우고 나서 남는 것이 진짜 물리라는 것.

이중 평균법(two-timing 기반의 double averaging)은 빠른 시간척도와 느린 시간척도가 분리된 계에서, 빠른 진동을 주기 평균으로 소거해 느린 유효 방정식을 얻되 1차 평균에서 사라지는 효과를 2차까지 밀어 붙잡는 섭동 기법이다. 1차 평균만으로는 “아무 일도 안 일어난다”가 나오는데 실험에서는 뒤집힌 진자가 멀쩡히 서 있는 상황 — 그 간극을 메우는 것이 두 번째 평균이다.

가장 유명한 응용이 카피차 진자(Kapitza pendulum)다. 지지점을 고속으로 위아래 흔들면 물리적으로 불안정해야 할 거꾸로 선 평형점이 안정화된다. 1차 평균에서 강제력의 평균은 정확히 0이라 아무 설명도 못 하고, 2차 항에서야 유효 퍼텐셜의 골짜기가 나타난다.

2. 표준 평균화 정리 — 1차[편집]

주기 강제를 받는 느린 계를 표준형(standard form)으로 쓴다.

\dot{x} = \varepsilon f(x,t) + \varepsilon^2 g(x,t,\varepsilon), \qquad 0 < \varepsilon \ll 1

여기서 $f$ 는 $t$ 에 대해 주기 $T$ 를 갖는다. 시간 평균을

\bar{f}(y) = \frac{1}{T}\int_0^{T} f(y,s)\,ds

로 정의하고 평균화된 계

\dot{y} = \varepsilon\,\bar{f}(y), \qquad y(0)=x(0)

를 세운다. 크릴로프-보골류보프-미트로폴스키 계열의 평균화 정리가 보장하는 것은 다음과 같다.

\lVert x(t) - y(t) \rVert = O(\varepsilon) \qquad \text{on } 0 \le t \le L/\varepsilon

읽는 법이 중요하다. 오차는 $O(\varepsilon)$ 이지만 그 보장이 유효한 구간이 $1/\varepsilon$ 이다. $x$ 자체가 그 구간에서 $O(1)$ 만큼 변하므로 이건 자명한 진술이 아니다. 순진하게 ” $\varepsilon$ 이 작으니 무시”하고 짧은 시간만 보는 것과는 차원이 다른 주장이며, 근사식을 인용할 때는 이 시간 스케일 딱지를 항상 함께 붙여 읽어야 한다. 그리고 이 보장에는 비공명 조건이 붙는다.¹

기하학적 뒷면도 있다. 평균화는 근항등 변환(near-identity transformation)

x = y + \varepsilon u_1(y,t) + \varepsilon^2 u_2(y,t) + \cdots, \qquad u_1(y,t) = \int_0^{t}\bigl[f(y,s) - \bar{f}(y)\bigr]ds

로 진동 성분을 좌표변환에 밀어 넣는 것이다. $u_1$ 의 피적분함수가 평균 0이므로 $u_1$ 은 유계 주기함수이고, 남은 $y$ 방정식에서 진동이 사라진다. 이 형식은 정규형 이론의 좌표변환과 사실상 같은 기계이며, 실제로 평균화는 주기적으로 강제된 계에 대한 정규형 계산이라고 봐도 된다.

3. 두 번째 평균이 필요한 순간[편집]

위 변환을 $\varepsilon^2$ 까지 밀면

\dot{y} = \varepsilon\,\bar{f}(y) + \varepsilon^2 \bar{f}_2(y) + O(\varepsilon^3), \qquad \bar f_2(y) = \overline{\;\frac{\partial f}{\partial y}u_1 + g - \frac{\partial u_1}{\partial y}\bar f\;}

가 되고, 근사 오차는 같은 $1/\varepsilon$ 구간에서 $O(\varepsilon^2)$ 로 내려간다. 이것이 2차 평균화, 이 문서가 말하는 이중 평균이다.

정확도가 한 자릿수 좋아지는 것이 전부라면 그저 편리한 개선일 뿐이다. 진짜 이야기는 $\bar f \equiv 0$ 인 경우다. 이때 1차 평균은 “느린 변수는 아무 데도 안 간다”고 말하고, 계의 실제 느린 거동은 전부 $\varepsilon^2 \bar f_2$ 가 지배한다. 시간척도도 $1/\varepsilon$ 이 아니라 $1/\varepsilon^2$ 으로 밀려난다. 고주파 진동 강제가 딱 이 상황이다 — 강제력의 주기 평균은 0이지만, 진동하는 강제가 자기가 만든 진동 응답과 상관(correlate) 되면서 0이 아닌 평균력을 남긴다. 이 잔여력이 바로 다음 절의 유효 퍼텐셜이다.

핵심 표현을 한 줄로 적으면 이렇다. 위치 $x$ 의 입자가 퍼텐셜 $U(x)$ 와 고주파 힘 $f(x,t)$ (주기 평균 0, 각진동수 $\omega$ 가 계의 고유 진동수보다 훨씬 큼)를 함께 받을 때, 느린 운동은

U_{\mathrm{eff}}(x) = U(x) + \frac{\langle f^2 \rangle}{2m\omega^2}

라는 유효 퍼텐셜(폰데로모티브 퍼텐셜, ponderomotive potential) 안의 운동으로 기술된다. 여기서 $\langle\cdot\rangle$ 은 빠른 주기에 대한 평균이다. 강제력의 제곱이 들어오는 것이 요점이다. 평균이 0이어도 제곱의 평균은 0이 아니다.

4. 카피차 진자[편집]

이름은 표트르 카피차의 1951년 논문에서 왔다.² 길이 $l$ 의 진자에서 지지점을 연직으로 $y_p(t) = a\cos\omega t$ 만큼 흔든다. 아래쪽 연직 방향에서 잰 각을 $\theta$ 라 하면, 지지점 가속도가 만드는 관성력을 실효 중력에 합쳐

\ddot{\theta} = -\frac{1}{l}\bigl(g - a\omega^2\cos\omega t\bigr)\sin\theta

를 얻는다. 조건은 $a \ll l$ 과 $\omega \gg \sqrt{g/l}$ 다. 이제 각을 느린 성분과 빠른 성분으로 쪼갠다.

\theta(t) = \Theta(t) + \xi(t), \qquad |\xi| \ll 1,\quad \langle \xi \rangle = 0

빠른 성분. 지배적인 고주파 항만 남기면 $\ddot{\xi} \approx (a\omega^2/l)\cos(\omega t)\sin\Theta$ 이고, $\Theta$ 를 상수 취급해 두 번 적분하면

\xi(t) = -\frac{a}{l}\cos(\omega t)\,\sin\Theta

느린 성분. $\sin(\Theta+\xi) \approx \sin\Theta + \xi\cos\Theta$ 로 전개하고 한 주기 평균하면, $\langle\cos\omega t\rangle = 0$ 이라 1차 항은 통째로 죽고 $\xi$ 와의 상관항만 남는다.

\langle \cos(\omega t)\,\xi(t)\rangle = -\frac{a}{l}\sin\Theta\,\langle\cos^2\omega t\rangle = -\frac{a}{2l}\sin\Theta

따라서

\ddot{\Theta} = -\frac{g}{l}\sin\Theta - \frac{a^2\omega^2}{2l^2}\sin\Theta\cos\Theta

이고, 이는 유효 퍼텐셜(단위 $m l^2$ 당)

V(\Theta) = -\frac{g}{l}\cos\Theta + \frac{a^2\omega^2}{4l^2}\sin^2\Theta

에서의 운동이다. $\ddot\Theta = -dV/d\Theta$ 를 확인하면 바로 맞는다.

거꾸로 선 평형의 안정화. $\Theta = \pi$ 에서

V''(\pi) = -\frac{g}{l} + \frac{a^2\omega^2}{2l^2} > 0 \quad\Longleftrightarrow\quad (a\omega)^2 > 2gl

즉 지지점 속도 진폭 $a\omega$ 가 $\sqrt{2gl}$ 을 넘으면 위쪽 평형점이 국소 최소가 되어 안정해진다. 이때 그 주위의 미소 진동 각진동수는

\omega_{\mathrm{eff}}^2 = \frac{a^2\omega^2}{2l^2} - \frac{g}{l}

이다. 실험 영상에서 거꾸로 선 진자가 옆으로 밀어도 느릿하게 되돌아오는 그 진동수가 이것이다. 아래쪽 평형( $\Theta=0$ )은 물론 여전히 안정이고, 두 평형 사이에는 $\cos\Theta = -2gl/(a^2\omega^2)$ 에서 불안정 평형이 새로 생긴다 — 고주파 흔들기가 분기 이론적으로 평형점 구조 자체를 바꾼 것이다.

5. 다중척도 전개와 공명 항[편집]

같은 기계를 형식적으로 정리한 것이 다중척도법이다. 독립 시간변수를 여러 개 도입하고

T_0 = t,\quad T_1 = \varepsilon t,\quad T_2 = \varepsilon^2 t, \qquad \frac{d}{dt} = \partial_{T_0} + \varepsilon\,\partial_{T_1} + \varepsilon^2\partial_{T_2} + \cdots

해를 $x = x_0(T_0,T_1,\dots) + \varepsilon x_1 + \cdots$ 로 전개한다. 이때 각 차수의 방정식에 **공명 강제 항(secular term)**이 생기는데, 이를 방치하면 $t\sin t$ 처럼 시간에 비례해 커지는 항이 나와 전개가 $t \sim 1/\varepsilon$ 에서 무너진다. 느린 변수 $T_1, T_2$ 의 의존성을 자유도로 써서 이 항들의 계수를 0으로 만드는 것이 **가해조건(solvability condition)**이고, 그 결과가 곧 평균화된 방정식이다.

교과서 예제인 무감쇠 더핑 진동자

\ddot{x} + x + \varepsilon x^3 = 0

에서 순진한 전개는 $x_1$ 에 $t\sin t$ 를 낳는다. 다중척도로 처리하면 진폭 $a$ 가 상수로 유지되고 위상만 느리게 밀려, 진동수가

\omega \approx 1 + \frac{3}{8}\varepsilon a^2

로 진폭에 의존한다는 결과가 나온다. 비선형 진동의 골자 — 진폭이 진동수를 바꾼다 — 가 공명 항 제거 한 줄에서 떨어진다. 강제·감쇠가 붙으면 같은 절차가 공명 곡선의 휘어짐과 도약(jump) 현상을 준다.

6. 같은 기계, 다른 이름들[편집]

이중 평균은 물리 분야마다 다른 이름으로 재발견되어 왔다. 아래는 전부 “빠른 위상을 평균하고 2차 잔여 효과를 취한다”의 변주다.³

분야	빠른 위상	평균이 남기는 것
역학	지지점 진동 $\omega t$	카피차 유효 퍼텐셜
이온 트랩	RF 구동 $\Omega t$	유사퍼텐셜, secular 진동
플라즈마	자이로 위상	유도중심 근사, 단열 불변량 $\mu$
천체역학	궤도 근점이각	장주기 섭동, 세속 진화

파울 트랩(Paul trap)이 특히 깔끔한 사례다. 정전기만으로는 3차원 가둠이 불가능하지만(라플라스 방정식이 안장점만 허용한다), 사중극 전극에 각진동수 $\Omega$ 의 RF를 걸면 위 공식이 그대로 적용되어

\Psi(\mathbf{r}) = \frac{q^2\lvert \mathbf{E}_0(\mathbf{r})\rvert^2}{4m\Omega^2}

라는 유사퍼텐셜(pseudopotential)이 생기고, 그 안에서 이온이 느린 secular 진동을 한다. 정확한 운동방정식은 마티외 방정식이고 안정 영역이 $q \lesssim 0.9$ 로 유한한데, 이중 평균은 그 안정 영역 깊숙한 곳( $q \ll 1$ )에서의 유효 기술이다. 안정 경계 근처에서는 평균화 전제가 무너지고 마티외 해석이 필요하다는 것이 이 근사의 유효 범위를 정직하게 말해 준다.

유도중심 근사는 규모가 다르다. 자화 플라즈마에서 입자는 자이로 반경 $\rho$ 로 자이로 진동수 $\Omega_c = qB/m$ 에 맞춰 돌면서 자기력선을 따라 흐른다. 자이로 위상에 대해 평균하면 자기 모멘트 $\mu = mv_\perp^2/(2B)$ 가 단열 불변량으로 살아남고, 1차에서 $\mathbf{E}\times\mathbf{B}$ 표류가, 2차에서 grad- $B$ 표류와 곡률 표류가 나온다. 이 평균화가 플라즈마 시뮬레이션에서 갖는 의미는 절대적이다 — 자이로운동학(gyrokinetics)은 6차원 위상공간을 5차원으로 줄이고 가장 빠른 시간척도를 제거해, 토카막 난류 시뮬레이션을 계산 가능한 문제로 만든다. 없었으면 사이클로트론 주기를 분해하느라 아무것도 못 돌린다.

7. 수치해석에서[편집]

수치 하는 사람이 이 이론에 관심을 가져야 하는 이유는 단순하다. 빠른 진동을 명시적으로 적분하면 망한다.

빠른 진동수 $\omega$ 를 분해하려면 $\Delta t \ll 1/\omega$ 가 필요한데, 정작 보고 싶은 물리는 $t \sim 1/\varepsilon^2$ 에서 일어난다. 필요한 스텝 수가 $\omega/\varepsilon^2$ 로 폭발한다. 강성 방정식 문서의 강성이 감쇠 스케일 분리라면, 이쪽은 진동 스케일 분리로 인한 강성이며 음함수 적분기로도 해결되지 않는다. 음함수법은 빠른 성분을 감쇠시켜 없애 버리는데, 이중 평균이 말하듯 그 빠른 성분의 상관이 정확히 우리가 원하는 유효력이기 때문이다.

실무적 처방은 세 갈래다.

해석적으로 평균해 버린다. 위 유도로 유효 방정식을 손으로 뽑고 그것만 적분한다. 가장 싸고 가장 정확하지만 문제마다 사람이 해야 한다. 자이로운동학 코드가 이 길이다.
수치적으로 평균한다. 미시 솔버를 짧은 구간만 돌려 힘의 평균을 추정하고 거시 스텝을 크게 밟는 방식 — 이종 다중척도법(HMM)과 스트로보스코픽 평균화 계열이다. 유효 방정식을 손으로 못 뽑을 때의 대안.
빠른 성분을 아예 구속으로 대체한다. 구속 동역학의 SHAKE·RATTLE이 분자동역학에서 결합 신축 진동을 없애 $\Delta t$ 를 늘리는 것이 대표적이다. 다만 이건 평균이 아니라 제거라서, 카피차형 잔여 효과가 중요한 문제에는 쓸 수 없다.

주의할 함정도 있다. 평균 방정식만 적분하면 빠른 진동의 진폭 정보가 사라진다. 물리적 관측량이 $\theta$ 인데 $\Theta$ 만 갖고 있으면 $O(\varepsilon)$ 만큼 틀린 답을 보게 되므로, 근항등 변환 $x = y + \varepsilon u_1$ 을 역으로 적용해 진동 성분을 복원해야 한다. 초기조건도 마찬가지다 — $y(0) = x(0)$ 로 그냥 두면 $O(\varepsilon)$ 오차가 초기부터 들어간다. 2차 정확도를 원하면 초기조건도 2차까지 보정해야 한다는 것이 이중 평균을 실제로 쓸 때 가장 자주 잊히는 항목이다.⁴

8. 이름에 대한 주의[편집]

“이중 평균(double averaging)“이라는 표현은 다른 분야에서 다른 뜻으로도 쓰인다. 검색할 때 헷갈리기 쉬우니 정리해 둔다.

이중 평균 나비에-스토크스(double-averaged Navier–Stokes). 거친 하상 위 개수로 유동에서 시간 평균에 이어 하상 요철 규모의 공간 평균을 추가로 취하는 하천 수리학의 기법이다. 시간척도 분리가 아니라 시간+공간 평균이라 이 문서의 내용과 다르다.
천체역학의 이중 평균. 계층적 삼체계에서 내부 궤도와 외부 궤도의 두 궤도 주기 모두에 대해 평균해 장기 세속 진화(코자이-리도프 진동 등)를 얻는 것. 이쪽은 이 문서와 같은 정신이지만, 평균 대상이 두 개의 빠른 위상이라는 점에서 스케일 하나짜리 두 번 밀기와는 구조가 다르다.

9. 관련 문서[편집]

반대로 이중 평균이 안 되는 지점도 분명하다. 빠른 진동수가 느린 계의 고유 진동수와 유리수 비를 이루면 공명이 걸려 평균화가 무너진다. KAM 정리에서 공명 원환면부터 부서지는 것과 같은 사정이며, 평균화 정리의 증명이 비공명 조건을 요구하는 이유이기도 하다. ↩
카피차(Pyotr Kapitza)는 초유동 연구로 1978년 노벨 물리학상을 받은 그 카피차 맞다. 진동 지지점 진자 자체는 그보다 훨씬 오래된 문제(스티븐슨, 1908)이지만, 유효 퍼텐셜로 깔끔하게 설명한 1951년 논문 덕에 이름이 붙었다. 좋은 설명은 발견만큼 값이 나간다는 사례. ↩
이 원리를 반대로 쓰면 진동으로 물성을 바꾸는 공학이 된다. 진동하는 체에서 분말이 유체처럼 흐르는 것, 진동 급이기(vibratory feeder)가 부품을 한 방향으로만 밀어내는 것, 진동으로 구조물의 유효 강성을 바꾸는 것 전부 같은 계열의 아이디어다. 흔들어서 되는 일이 생각보다 많다. ↩
“평균 방정식으로 바꿨는데 왜 오차가 안 줄죠”라는 질문의 8할은 초기조건 보정을 안 해서다. 이론이 약속한 $O(\varepsilon^2)$ 은 좌표변환 전체에 대한 약속이지 $y$ 방정식만의 약속이 아니다. ↩