입자 사슬에서 파동 방정식으로

두 질량에서 연속적인 장으로

결합 진동에서는 두 질량의 위치를 두 좌표로 기술했다. 그런데 줄은 모든 위치에서 변위를 갖는 것처럼 보인다. 파동 방정식은 이 두 설명을 연결한다. 같은 물리적 재료와 같은 초기 모양을 유지하면서 작은 요소의 수를 늘리면, 좌표의 목록이 연속적인 장의 기술로 이어진다.

핵심 질문은 이산 기계 모델이 언제 같은 연속체 실험을 근사하는가이다. 애니메이션에 점을 더 그리는 것만으로는 답할 수 없다. 질량·결합·공간 파장·초기속도·에너지를 정하고 실제 운동을 비교해야 한다. 이 단원은 서로 독립적인 초기조건을 갖는 여러 모드의 주기 사슬을 사용한다. 앞의 단일 Fourier 모드 PDE 비교에 물리적 격자 세분화 실험을 더한다.

물리적 가정과 SI 단위

평형 위치 $x_j=ja$의 횡변위를 $q_j(t)$라 하자. 전체 길이는 $L=Na$이며 인덱스는 $q_N=q_0$으로 순환한다. 주기 영역은 반복되는 패턴을 뜻한다. 자유단이나 고정단을 가진 줄과 경계 조건이 다르다.

선밀도 $\mu$의 단위는 kg/m, 장력 $\tau$의 단위는 N이다. 입자 질량과 이웃 결합 강성은

$$m=\mu a,\qquad \kappa=\frac{\tau}{a},\qquad m\ddot q_j=\kappa(q_{j+1}-2q_j+q_{j-1})$$

로 둔다. $q$는 m, $m$은 kg, $\kappa$는 N/m이다. 일정한 장력과 작은 횡기울기가 줄의 힘을 선형화할 근거다. 굽힘 강성·감쇠·장력 변화·기하 비선형성은 포함하지 않는다. 큰 기울기를 입력해도 선형 방정식의 계산값은 나오지만 실제 줄을 정확히 나타낸다는 뜻은 아니다.

이웃 입자의 힘을 유도하기

오른쪽 줄 조각의 작은 기울기는 $(q_{j+1}-q_j)/a$다. 장력의 횡방향 성분은 약 $\tau(q_{j+1}-q_j)/a$이고, 왼쪽 조각은 $-\tau(q_j-q_{j-1})/a$를 준다. 두 힘을 더하면 위의 이차 차분식이 나온다.

한 입자만 위로 올리고 이웃을 그대로 두었을 때 힘은 아래로 향해야 한다. 부호를 뒤집으면 진동하는 용수철 대신 불안정한 힘이 된다. 또한 $x_j$는 평형 좌표다. 횡변위가 증가한다고 질량이 그림의 가로축을 따라 이동하는 것은 아니다.

그림이 아닌 같은 재료를 세분화하기

$L,\mu,\tau$를 고정하고 $N$을 두 배로 늘린다. 간격과 입자 질량은 절반, 결합 강성은 두 배가 된다.

$$a\mapsto a/2,\quad m\mapsto m/2,\quad \kappa\mapsto2\kappa,\quad \frac{\kappa a^2}{m}=\frac{\tau}{\mu}=c^2.$$

예를 들어 $L=8$ m, $\mu=1$ kg/m, $\tau=1$ N에서 $N=16$이면 $a=0.5$ m, $m=0.5$ kg, $\kappa=2$ N/m이다. $N=32$이면 각각 $0.25$ m, $0.25$ kg, $4$ N/m이다. 두 경우 모두 총질량은8 kg이고 연속체 파동 속도는 $c=1$ m/s다.

$m,\kappa$를 고정한 채 $a$만 줄이면 선밀도와 장력이 바뀐다. 그것은 다른 재료의 열이며 이 실험의 연속체 수렴을 입증하지 못한다. 무엇을 고정했는지를 수치 실험의 조건으로 먼저 적어야 한다.

Taylor 전개와 근사의 조건

매끄러운 연속 변위 $u(x,t)$에 대해

$$u(x+a)-2u(x)+u(x-a) =a^2u_{xx}+\frac{a^4}{12}u_{xxxx}+\cdots$$

이다. 운동방정식에 넣고 $\mu a$로 나누면

$$u_{tt}=c^2u_{xx}+\frac{c^2a^2}{12}u_{xxxx}+\cdots, \qquad c=\sqrt{\tau/\mu}$$

를 얻는다. 선두 항만 남긴 식이 파동 방정식이다. 공간 파수 $\xi$의 성분에서 무차원 조건은 $|\xi a|\ll1$이다. 간격을 미터로 썼을 때 숫자가 작다는 사실만으로 충분하지 않다. 초기 모양이 그 간격에 비슷한 크기로 급변하면 높은 공간 미분 항이 중요해진다.

이산 스펙트럼과 연속체 스펙트럼

$q_j(t)=\exp(i\xi ja-i\omega_at)$를 이산 방정식에 넣고 이웃의 위상 인자를 정리하면

$$\omega_a^2=\frac{4c^2}{a^2}\sin^2\!\left(\frac{\xi a}{2}\right), \qquad \omega=c|\xi|$$

다. 여기서 사용하는 음이 아닌 해상 가능한 모드에서는

$$\frac{\omega_a}{\omega} =\frac{\sin(\xi a/2)}{\xi a/2} =1-\frac{(\xi a)^2}{24}+O((\xi a)^4)$$

가 된다. 같은 물리적 파수에서 격자의 진동수는 연속체보다 낮다. 격자의 분해능 한계 근처에서는 차이가 크다. 두 정확한 해 중 하나가 수치적으로 불안정하다는 뜻이 아니라 두 방정식의 분산 관계가 다르다는 뜻이다. 일반적인 파동 묶음과 군속도는 뒤의 별도 주제에서 다룬다.

초기 변위와 초기속도를 독립적으로 주기

$\xi_n=2\pi n/L$인 유한 실수 Fourier 전개를 사용하자.

$$f(x)=\sum_n[A_n\cos(\xi_nx)+B_n\sin(\xi_nx)], \qquad g(x)=\sum_n[C_n\cos(\xi_nx)+D_n\sin(\xi_nx)].$$

변위 계수의 단위는 m, 속도 계수는 m/s다. 줄을 휘게 놓는 것과 초기 충격을 가하는 것은 서로 다른 모양일 수 있다. 사슬에는 $q_j(0)=f(x_j)$와 $\dot q_j(0)=g(x_j)$를 주고, 연속체에는 같은 물리적 함수 $f,g$를 준다.

이 실험에서는 서로 다른 정수 모드를 $0\le n<N/2$로 제한한다. 해상하지 못하는 모드를 조용히 다른 모드로 바꾸는 aliasing과 Nyquist 모드의 특별한 사인 좌표를 피하기 위해서다. 영모드는 $\sin0=0$이므로 코사인 계수만 가진다.

영모드를 포함하는 정확한 시간 진화

각진동수 $w\ne0$, 초기값 $A$, 초기속도 $C$인 모드 좌표는

$$z(t)=A\cos(wt)+C\frac{\sin(wt)}w, \qquad \dot z(t)=-Aw\sin(wt)+C\cos(wt)$$

다. 사슬에서는 $w=\omega_a$, 연속체에서는 $w=\omega$를 사용한다. 구현은 $\sin(wt)/w=t\operatorname{sinc}(wt)$와 $\operatorname{sinc}(0)=1$을 사용한다. $w=0$이면 $A+Ct$다. 계산이 특이해지거나 무조건 정지하는 것이 아니다.

상수 초기속도는 줄 전체의 횡방향 병진을 만든다. 변형 에너지는 없지만 운동에너지와 운동량은 있다. 그림을 가운데 놓으려고 이 운동을 빼 버리면 초기값 문제가 달라진다. 평균 위치뿐 아니라 평균 속도도 확인해야 한다.

점과 매끄러운 곡선이 뜻하는 것

사슬에는 실제 좌표 $N$개가 있다. 이 점들을 통과하는 매끄러운 삼각함수 보간선은 표시를 돕는 도구이며 추가 질량이 아니다. 연속체 장은 표본 사이에도 정의된다. 입자 위치에서 연속체를 평가하면 같은 위치의 변위를 비교할 수 있다.

시각 $t$에서 사용할 수 있는 오차는

$$\varepsilon_N(t)=\left[\frac1N\sum_{j=0}^{N-1} |q_j(t)-u(x_j,t)|^2\right]^{1/2}$$

이며 단위는 m이다. 특정 시점의 작은 오차는 우연한 위상 정렬 때문일 수 있다. 여러 시점과 진동수 비를 함께 본다. 그림의 픽셀이나 표시 표본만 늘려도 기계적 입자 수 $N$과 고유진동수는 바뀌지 않는다.

실제 에너지와 주기 연결

각 연결 용수철을 한 번만 세면 사슬의 에너지는

$$E_a=\sum_j\left[\frac12m\dot q_j^2+ \frac12\kappa(q_{j+1}-q_j)^2\right]$$

다. 연속체의 에너지는

$$E=\int_0^L\left[\frac12\mu u_t^2+ \frac12\tau u_x^2\right]dx$$

다. 각각의 정확한 운동에서는 모두 보존된다. 그러나 유한 간격에서 초기값까지 같을 필요는 없다. 같은 Fourier 모드에 대해 이차 차분과 공간 미분의 응답이 다르기 때문이다. 긴 파장 세분화 극한에서 변형 에너지도 수렴한다. 두 계산이 에너지를 잘 보존한다는 사실만으로 같은 물리를 풀었다고 결론내릴 수 없다.

영이 아닌 해상 가능한 실수 모드의 사슬 에너지는

$$E_{a,n}=\frac{\mu L}{4} [C_n^2+D_n^2+\omega_{a,n}^2(A_n^2+B_n^2)]$$

이며 영모드는 $\mu LC_0^2/2$다. 직교성 때문에 이 모드 에너지들의 합은 실제 질량·연결 에너지의 합과 같다.

국소 에너지와 부호 있는 전달

입자마다 양쪽 연결 에너지의 절반씩을 배분하자.

$$E_j^*=\frac12m\dot q_j^2+ \frac14\kappa(q_j-q_{j-1})^2+ \frac14\kappa(q_{j+1}-q_j)^2.$$

$j$에서 $j+1$로 향하는 흐름을 양수로 두면

$$J_{j+1/2}=-\frac\kappa2(q_{j+1}-q_j) (\dot q_{j+1}+\dot q_j), \qquad \dot E_j^*=J_{j-1/2}-J_{j+1/2}$$

다. 실제 배분 에너지를 힘 방정식으로 미분해 이 항등식을 검사한다. 흐름의 단위는 W다. 음수는 반대 방향으로의 물리적 전달이며 보존 실패가 아니다. 전체 주기 사슬에 더하면 내부 흐름이 모두 상쇄된다.

계산 예제: 같은 초기 모양의 세분화

$L=8$, $\mu=\tau=1$이고 초기 변위를 $0.05\cos(2\pi x/L)$, 초기속도를0으로 두자. $N=16$이면 $\xi a=\pi/8$이며 진동수 비는 $\operatorname{sinc}(\pi/16)\approx0.99359$다. $N=32$에서는 약 $0.99839$다.

간격을 절반으로 하면 선두 진동수 오차는 약4분의1이 된다. 해상 가능한 여러 모드와 관측 시각을 고정해 충분히 세분화하면 변위 오차도 이차 수렴 경향을 보인다. 모든 개별 시각의 오차가 단조롭게 감소한다는 보장은 아니다. 더 긴 시간에는 위상 차가 누적될 수 있다. 비교할 때 초기 함수·단위·관측 시간을 모두 고정한다.

안내 비교와 모델 범위

먼저 $N$을 두 배로 하면 $c$가 바뀔지 예측한다. 궤적보다 질량과 강성을 먼저 비교한다. 다음으로 초기 변위에 없는 속도 모드를 추가한다. 보기 좋은 위상을 임의로 택한 결과가 아니라 지정한 초기속도가 재구성되는지 본다.

낮은 모드와 해상 한계에 가까운 모드를 비교한다. 둘 다 유효한 격자 운동일 수 있지만 연속체 근사가 좋은 것은 앞의 모드다. 마지막으로 상수 속도를 더해 평균의 선형 운동과 운동량을 확인한다. 계산기는 입자 수8·16·32·64·128, 길이4–16 m, 양의 선밀도와 장력0.25–4 SI 단위, 유한 계수−1–1, 시간0–64 s를 지원한다. 이 범위는 교육용 계산의 범위이며 파동 이론 자체의 제한이 아니다.

다른 표현으로 검증하기

정확한 스펙트럼 구현을 같은 모드 공식을 다시 호출해서만 검사해서는 안 된다. 이웃 힘에서 주기 이차 차분 행렬을 조립하고, 위치·속도의 일차 상태 행렬을 만든 뒤 독립적인 행렬지수로 진화시킨다. 별도로 정밀도를 높인 ODE 해도 대조한다. 마지막 에너지 합뿐 아니라 실제 좌표·힘·국소 에너지·링크 흐름을 비교한다.

연속체와 비교할 때는 같은 물리 Fourier 계수를 유지하고 분산 관계만 바꾼다. 보존은 각각의 진화를 검증하고 세분화는 두 모델의 관계를 검증한다. 특히 유한 간격에서는 이산·연속체의 초기 변형 에너지가 다르다. 같은 총합을 표시하려고 한쪽 에너지로 다른 쪽을 대체하면 이 실험에서 측정하려는 모델링 차이를 숨기게 된다.

해설이 있는 연습 문제

문제 1: 같은 줄을 세분화하기

$L=8$, $\mu=2$, $\tau=8$, $N=16$일 때 $m,\kappa,c$를 구하고 $N$을 두 배로 하세요.

해설: $a=0.5$, $m=1$ kg, $\kappa=16$ N/m, $c=2$ m/s다. 두 배로 세분화하면 $a=0.25$, $m=0.5$, $\kappa=32$이고 속도는2로 유지된다. 이 예제의 장력8은 교육용 계산기의 범위를 넘지만 유도식은 그대로 성립한다.

문제 2: 줄 전체의 횡방향 병진

모든 입자의 초기 위치가0.1 m, 속도가0.05 m/s라면12초 뒤 어떻게 되나요?

해설: 모든 변위는0.7 m다. 이웃 차이와 힘은0이다. 에너지는 $\mu L(0.05)^2/2$의 운동에너지뿐이다. 영주파수라고 초기속도까지0인 것은 아니다.

문제 3: 네 간격 길이의 파장

$\lambda=4a$일 때 $\omega_a/\omega$를 계산하세요.

해설: $\xi a=2\pi a/\lambda=\pi/2$이므로 비는 $\sin(\pi/4)/(\pi/4)\approx0.90032$다. 격자가 모드를 해상한다는 조건은 긴 파장 근사의 조건보다 약하다.

문제 4: 연결을 한 번 셀까 두 번 셀까

입자별 배분 에너지에는 연결 두 개가 있는데 전체 에너지에서는 왜 각 연결을 한 번만 세나요?

해설: 각 입자가 양쪽 연결의 절반씩을 받으므로 전체 합에서 한 연결이 두 번 등장하되 각각 가중치가 절반이다. 입자마다 연결 두 개의 전체 에너지를 넣으면 변형 에너지를 두 배로 센다.

문제 5: 늦은 시점의 작은 오차

특정한 늦은 시점에서 $\varepsilon_N(t)$가 매우 작으면 연속체 근사가 좋다는 증거인가요?

해설: 충분하지 않다. 다른 진동수도 잠시 비슷한 위상을 가질 수 있다. 여러 시점·진동수 비·같은 물리 초기조건을 유지한 세분화 열을 비교해야 한다.

문제 6: 변위 없이 주는 초기 충격

한 모드에서 $A=B=D=0$, $C\ne0$일 때 변위와 초기 에너지는 무엇인가요?

해설: 공간 코사인 성분의 변위 진폭은 $C\sin(\omega_at)/\omega_a$, 속도 진폭은 $C\cos(\omega_at)$다. 초기 에너지는 $\mu LC^2/4$다. 영모드에서는 변위가 $Ct$로 이어지고 에너지의 정규화 계수는4분의1 대신2분의1이다.

AI 모델과의 연결 및 읽을 자료

학습 시뮬레이터는 한 격자의 궤적을 잘 맞추면서도 잘못된 재료 스케일링을 학습할 수 있다. 같은 물리적 초기 장을 유지해 $N$을 바꾸고 보존·국소 전달·연속체 오차를 따로 평가한다. 스펙트럼 신경 연산자에서도 물리적 파수와 배열 인덱스를 구분해야 한다. 이 실험은 검증 질문을 제공하며 학습 모델의 일반화 정리를 대신하지 않는다.

MIT8.03SC 강의9는 사슬에서 파동 방정식으로 가는 과정을 다룬다. 여기의 다중 모드 초기값·SI 에너지·국소 흐름·세분화 계산은 명시한 기계 모델에서 직접 유도했다. 다음 단원은 반복되는 주기 영역 대신 국소 초기 파형과 실제 경계를 사용한다.