격자 파동 묶음의 실제 에너지 전달

파동 묶음이 실제 에너지를 옮기는 방식

Gaussian 포락선으로 봉우리의 움직임은 설명할 수 있지만 정규화된 장의 제곱이 자동으로 물리적 에너지가 되는 것은 아니다. 이번에는 질량과 선형 용수철로 구성된 사슬을 사용한다. 운동 에너지·탄성 에너지·전력의 정의가 명확하므로, 군속도를 그림의 인상 대신 실제 전달량과 비교할 수 있다.

사슬에는 주기 경계가 있다. 오른쪽 끝으로 나간 묶음이 왼쪽에 나타나는 것은 두 위치가 이웃하기 때문이지, 경계가 흡수했다가 다시 생성했기 때문이 아니다. 실제 입자의 운동, 복소 해석장 포락선, 입자에 배분한 역학적 에너지와 링크 전력을 구분한다. 그 분포와 중심은 서로 같을 필요가 없다.

역학 변수와 가정의 범위

길이 $L$인 사슬에 $N$개 입자를 놓고 간격을 $a=L/N$으로 둔다. 선밀도는 $\rho$, 장력 매개변수는 $\mathcal T$다. 입자 질량은 $m=\rho a$, 링크 강성은 $\kappa=\mathcal T/a$이며 단위는 각각 kg와 N/m다. 변위 $q_j$는 m, 속도 $v_j$는 m/s다. 모든 단위를 고정하면 뒤의 에너지와 흐름을 실제 차원으로 비교할 수 있다.

동일한 입자와 최근접 조화 링크, 주기 경계를 가진 무손실 선형 모델을 가정한다. 횡방향 줄의 근사로 해석하려면 상대 변위와 기울기도 충분히 작아야 한다. 큰 설정 진폭에서도 조화 방정식 자체는 계산할 수 있지만, 그 그림만으로 작은 기울기라는 물리 가정이 검증되지는 않는다.

입자는32개 또는64개이며, 영모드와 Nyquist 끝점을 제외한 양의 Fourier 모드만 준비한다. 무질서·감쇠·외부 구동·비선형 용수철은 없다. 구현은 모든 위상에 대한 변위 기울기 상한 $\sum_n k_n A_n$도 반환한다. 상한이1에 가깝거나 크면 작은 기울기 줄 근사를 무조건 적용해서는 안 되지만, 이 값이 모든 재료 오차를 인증하는 것도 아니다. 제한된 모델이 어떤 질문에 답하는지 먼저 정하고 실제 모든 재료의 성질로 확대하지 않는다.

격자의 분산 관계 유도하기

Newton 방정식은

$$m\ddot q_j=\kappa(q_{j+1}-2q_j+q_{j-1})$$

다. $q_j\propto e^{i(kja-\omega t)}$를 대입하면

$$\omega^2=\frac{4\kappa}{m}\sin^2\frac{ka}{2},\qquad \omega(k)=\frac{2c}{a}\sin\frac{ka}{2},\quad c=\sqrt{\mathcal T/\rho},\quad 0<ka<\pi$$

를 얻는다. 선택한 양의 파수 대역에서 정현 함수가 양수이므로 이 진동수 가지를 쓴다. 대역 밖으로 확장하면서 절댓값이나 가지 관례를 생략하면 양의 진동수 해석이 틀릴 수 있다. 작은 $ka$에서는 $\omega\approx ck$가 되어 앞의 비분산 연속체 관계가 나온다. 격자 한계에 가까우면 이웃 입자의 진동을 매끄러운 연속체 파동으로 보는 근사가 약해진다.

위상속도와 군속도 및 곡률

정한 대역에서 세 양은

$$v_p(k)=\frac{\omega(k)}{k},\qquad v_g(k)=c\cos\frac{ka}{2},\qquad \omega''(k)=-\frac{ca}{2}\sin\frac{ka}{2}$$

다. 양의 해상 가능한 모드에서는 위상속도가 군속도보다 크다. 곡률은 음수이므로 좁은 대역 묶음의 포락선 chirp와 퍼짐도 그 부호에 맞게 나타난다. 격자 한계에 가까워질수록 진동수는0이 아니어도 군속도는 낮아진다. 높은 진동수가 곧 빠른 묶음 전달이라는 뜻은 아니다.

설정한 반송 모드 중심은 비정수여도 된다. 이는 스펙트럼 가중치의 중심이지, 정수 조건을 어긴 새로운 실제 모드를 추가한다는 뜻이 아니다. 실제 허용 파수는 정수 $n$에 대해 $k_n=2\pi n/L$다. 반송 중심에서 계산한 연속 미분값은 기준량이고, 실제 진화는 유한한 이산 스펙트럼을 사용한다.

파동 묶음과 정규화 정의하기

모드 계수를 다음과 같이 준비한다.

$$w_n=e^{-(n-n_c)^2/(2s_n^2)},\qquad A_n=A_{\max}\frac{w_n}{\sum_rw_r},\qquad U_j(t)=\sum_nA_n e^{i[k_n(ja-x_0)-\omega_nt+\phi_0]}.$$

합의 범위는 $n=1,\ldots,N/2-1$이고 실제 변위는 $q_j=\operatorname{Re}U_j$다. 삼각부등식에 의해 $|U_j|\le A_{\max}$이므로 초기 매개변수는 최고 진폭의 상한이다. 실제 입자 중 하나가 반드시 그 값을 갖는다는 약속은 아니다. 중심이 격자 사이에 있거나 초기 위상이 달라지면 실수 변위의 관측 최고점은 더 작을 수 있다.

Gaussian 모드 가중치는 양의 해상 가능한 대역에서 잘려 있다. 무한 실수선 Gaussian의 정확한 Fourier 변환과 같지 않다. 넓은 대역이나 격자 한계 부근에서는 잘림과 주기 복사본을 고려해야 한다. 최고 진폭 상한을 고정한 채 대역폭을 바꾸면 에너지도 달라진다. 실험마다 같은 에너지가 되도록 몰래 재정규화하지 않는다.

변위와 속도를 함께 준비하기

각 성분의 위상을 $\theta_{nj}=k_n(ja-x_0)-\omega_nt+\phi_0$라고 하면

$$q_j=\sum_nA_n\cos\theta_{nj},\qquad v_j=\sum_n\omega_nA_n\sin\theta_{nj}$$

다. 속도의 양의 부호는 $\cos(kx-\omega t)$를 시간으로 미분해서 나온다. 초기 속도를0으로 설정하면 같은 오른쪽 진행 묶음이 아니라 양방향 성분을 섞은 다른 초기 상태가 된다. 의도한 실험을 정하려면 초기 변위뿐 아니라 속도도 지정해야 한다.

이 모델은 유한 스펙트럼 합으로 정확하게 진화하므로 시간 간격을 반복하는 적분 오차는 없다. 그렇다고 식의 부호·정규화·초기 위상이 자동으로 맞는 것은 아니다. 독립적인 운동방정식 행렬의 지수함수와 대조한다. 실제 질량은 반환한 입자 위치에만 있으며, 부드러운 연결 곡선은 모드 보간의 시각화다.

역학적 에너지를 별도로 정의하기

입자 $j$의 운동 에너지와 링크 $j\to j+1$의 탄성 에너지는

$$K_j=\tfrac12mv_j^2,\qquad V_j=\tfrac12\kappa(q_{j+1}-q_j)^2$$

다. 이웃 링크의 절반씩을 입자에 배분하여

$$e_j=K_j+\tfrac12(V_{j-1}+V_j),\qquad E=\sum_je_j=\sum_j(K_j+V_j)$$

라 둔다. 이 값들의 단위는 J다. 국소 배분 방식이 달라지면 에너지 분포의 모양도 달라질 수 있으므로, 선택한 배분과 그에 대응하는 전력을 일관되게 사용해야 한다. 각 용수철을 양쪽 입자에서 전부 더하면 전체 에너지를 중복 계산한다.

해석장의 양 $a\sum_j|U_j|^2$는 $U$가 변위이므로 m³ 단위다. 스펙트럼 진화에서 이것도 보존되지만 역학적 에너지는 아니다. 같은 변위 진폭이어도 진동수가 높으면 더 큰 실제 에너지를 가진다. 그림의 포락선이 닮았다는 이유로 두 보존량을 같은 축에 놓지 않는다.

링크 전력과 국소 수지 유도하기

대칭적으로 배분한 입자 에너지에는 오른쪽 방향을 양수로 하는 링크 전력

$$J_j=-\frac{\kappa}{2}(q_{j+1}-q_j)(v_{j+1}+v_j)$$

가 대응한다. Newton 방정식과 링크 변형률을 사용해 $e_j$를 미분하면

$$\dot e_j=J_{j-1}-J_j$$

가 된다. 들어온 전력에서 나간 전력을 뺀 값이 해당 입자에 배분된 에너지의 변화다. 주기 사슬 전체에서 더하면 모든 링크가 상쇄되어 $\dot E=0$이다. 링크 전력의 단위는 W이며 $a\sum J_j/E$는 m/s다. 단위만으로도 전력과 에너지, 전달 속도를 구별할 수 있다.

여러 모드가 간섭하면 일부 위치에서 순간 전력이 음수일 수도 있다. 양의 군속도 가지로 준비했다는 이유로 모든 국소 화살표를 양수로 강제하지 않는다. 뒤의 전체 전달 항등식은 전체 공간 합과 정의한 에너지에 대한 식이지 모든 순간 화살표의 방향을 말하는 식이 아니다.

모드 에너지에서 올바른 가중치 찾기

해상 가능한 주기 입자에서 모드 직교성을 적용하면 준비한 진행 모드의 에너지는

$$E_n=\frac{mN}{2}\omega_n^2A_n^2,\qquad E=\sum_nE_n$$

다. 같은 변위 계수를 가진 정상파와 계수가 다른 이유는 진행 상태에 특정한 초기 속도가 함께 들어 있기 때문이다. 운동 에너지와 탄성 에너지가 모두 기여한다. 영모드와 Nyquist 모드를 제외한 조건도 직교성 논증에서 사용한다.

제곱 진폭 가중치와 에너지 가중치를 비교하면 좋은 음성 대조가 된다. 제곱 진폭에는 $\omega_n^2$가 없으므로 보통 다른 평균 속도가 나온다. 보기 좋은 포락선 가중치로 실제 역학적 에너지 가중치를 대신하면 전달 추정에 체계적 편향을 만들 수 있다. 무엇을 평균하는지 먼저 정하고 그 물리량의 가중치를 사용한다.

전체 에너지 전달의 정확한 항등식

이 유한 양의 모드 준비에서는

$$v_E=\frac{a\sum_jJ_j}{E} =\frac{\sum_nE_n v_g(k_n)}{\sum_nE_n}$$

가 성립한다. 실제 링크 전력을 전체 입자에서 더하면 서로 다른 해상 가능한 모드의 교차항이 상쇄되기 때문이다. 영장이 아니면 $0<v_E<c$다. 이는 정한 모델과 초기 준비의 정확한 항등식이지, 모든 매질의 군속도 봉우리가 항상 에너지 속도와 같다는 보편적 정리가 아니다.

대역이 충분히 좁으면 각 모드의 $v_g(k_n)$가 거의 같아서 $v_E\approx v_g(k_c)$다. 넓은 묶음에서는 반송 중심의 미분값과 크게 다를 수 있다. 장이 완전히0이면 분자·분모가 모두0이므로 전달 속도는 미정이다. 이것을 에너지가 정지했다고 측정한0 m/s로 바꾸지 않는다.

하나의 실제 준비를 수치로 비교하기

$L=8$ m, $N=32$, $\rho=0.8$ kg/m, $\mathcal T=1.7$ N, $A_{\max}=0.08$ m, $n_c=4$, 대역폭 $s_n=2$, $x_0=2.96$ m, $\phi_0=0.4$를 택한다. 간격은0.25 m, 입자 질량은0.2 kg, 링크 강성은6.8 N/m다. 연속체 기준 속도는 약1.4577 m/s다.

독립 운동방정식 계산에서 전체 에너지는 약0.06790960 J, 전체 에너지 전달 속도는1.28638809 m/s다. 반송 중심의 군속도는 약1.3468 m/s이므로 유한 대역과 에너지 가중치의 효과를 확인할 수 있다. 보존되는 해석장 노름은 약0.00776644 m³이며, 단위가 다른 것만으로도 같은 에너지가 아님을 알 수 있다.

초기에는 두 원형 중심이 모두2.96 m다. 그러나 집중도를 나타내는 원형 모멘트 크기는 해석장 강도에서 약0.9305, 배분 에너지에서 약0.9231로 다르다. 초기 중심이 같다고 이후 분포나 속도의 정의까지 같아지는 것은 아니다. 같은 준비를 시간에 따라 비교하며 어떤 양을 측정하는지 계속 표시한다.

원형 중심에는 집중도 진단이 필요하다

0과 $L$ 근처의 입자를 선형 평균하면 가운데 $L/2$ 부근으로 잘못 옮길 수 있다. 비음수 가중치 $w_j$에 대해 첫 원형 모멘트를 사용한다.

$$M=\frac{\sum_jw_je^{i2\pi x_j/L}}{\sum_jw_j},\qquad x_c=\frac{L}{2\pi}\arg M\pmod L.$$

각도와 함께 $|M|$도 보고한다. 균일 분포는 $M=0$이므로 선호되는 중심 각도가 없다. 서로 마주 보는 같은 두 봉우리도 첫 모멘트가0일 수 있어, 미정 중심이 곧 균일 분포라는 역명제는 성립하지 않는다. 모멘트가 아주 작으면 작은 오차에도 각도가 크게 흔들린다. 구현은 크기가 $10^{-10}$보다 작을 때 중심을 미정으로 두며, 그보다 크더라도 작은 값은 약하게 조건화된 측정이다.

해석장 강도와 실제 배분 에너지는 서로 다른 가중치로 계산한다. 원형 모멘트를 미분하면 위치에 따른 복소 위상 가중치가 들어가므로 원형 중심의 미분을 단순히 $a\sum J/E$와 같다고 할 수 없다. 미정 중심을 지나 무조건 각도를 펼치면 가짜 점프나 매우 빠른 이동을 만들어 낼 수 있다.

좁은 대역 근사와 수렴 실험

반송 파수 $k_c$ 근처에서 격자 관계를 위상·군속도·곡률로 전개하면 앞의 Gaussian 포락선 단원과 연결된다. 실제 스펙트럼이 국소화되어 있고 누적된 고차 위상 오차가 작을 때 유용하다. 넓은 대역, 격자 한계에서의 잘림, 영역을 돌아온 주기 묶음은 각각 단순 Gaussian 설명을 깨뜨릴 수 있다.

$L,\rho,\mathcal T$를 고정하고 $N$을 늘릴 때는 간격·입자 질량·링크 강성을 함께 바꾼다. 질량과 강성을 모두 고정하면 같은 물질을 세분화한 것이 아니라 다른 물질을 정의하게 된다. 수렴을 검사할 때는 물리적 초기 스펙트럼이나 준비 조건 중 무엇을 고정하는지도 명시한다.

유한 격자가 언제나 오차인 것은 아니다. 여기서는 그것이 실제 조화 격자의 정의다. 그 격자의 시간 진화를 잘못 계산한 수치 오차와, 격자가 이상적 연속 줄과 다른 물리적 차이를 구별한다. 그림을 더 부드럽게 그리는 것만으로 실제 입자 수를 늘렸다고 주장하지 않는다.

AI 예측 대상과 연구 검증

AI 모델은 초기 상태와 재료 매개변수에서 미래 변위·국소 에너지·링크 전력을 예측할 수 있다. 출력끼리도 일관되어야 한다. 변위 오차가 작아도 에너지는 속도와 공간 차분의 제곱을 포함하므로 정확하지 않을 수 있다. 전체 보존과 국소 유입·유출 수지를 각각 평가한다.

영진폭, 좁은 대역, 중심 병진, 격자 한계 근처의 넓은 대역을 별도 대조군으로 사용한다. 보지 못한 재료 비율과 대역폭 영역에서도 평가한다. 한 궤적의 인접 표본을 잘 맞춘 것만으로 다른 분산 관계에 일반화했다고 주장하지 않는다. 해석장 노름을 보존하도록 훈련한 것과 실제 에너지를 보존하는 것도 구분한다.

독립 기준값은 실제 실수 운동방정식 행렬의 지수함수를 사용하고, 별도 FFT로 해석장을 복원한다. 이를 유한 스펙트럼 모델과 대조하면 제품의 정현식을 다시 평가한 결과끼리 비교하는 것보다 강한 검사가 된다. 수학 기반을 검증한 뒤에도 독자 화면·설명·조작은 별도로 통합 검증해야 한다.

해설 문제

문제 1 — 같은 물질을 더 잘게 나누기

$L,\rho,\mathcal T$를 고정하고 입자 수를 두 배로 늘리면 질량과 강성은 어떻게 되는가?

해설: 간격이 절반이므로 $m=\rho a$도 절반, $\kappa=\mathcal T/a$는 두 배다. 기준 속도 $c=\sqrt{\mathcal T/\rho}$는 같다. 질량과 강성을 그대로 두면 선밀도와 장력에 해당하는 물리적 조건이 바뀌므로 같은 물질의 세분화 비교가 아니다.

문제 2 — 초기 속도의 부호

$q=A\cos(kx-\omega t)$의 $t=0$ 속도를 구하라.

해설: 시간 미분하면 $v=\omega A\sin(kx-\omega t)$이므로 초기 속도는 $\omega A\sin(kx)$다.0으로 설정하면 진행파 준비가 달라진다. 초기 변위 사진이 같더라도 속도가 다르면 이후 상태가 달라진다는 초기값 문제의 기본 원리를 확인할 수 있다.

문제 3 — 올바른 평균 속도

두 모드의 에너지가1 J와3 J이고 군속도가1.4 m/s와0.6 m/s이면 $v_E$는 얼마인가?

해설: 에너지 가중 평균은 $(1\cdot1.4+3\cdot0.6)/4=0.8$ m/s다. 단순 평균1 m/s는 어떤 모드에 실제 에너지가 더 많은지 무시한다. 진폭 가중치도 에너지 가중치와 일반적으로 다르므로 전달량의 정의에 맞게 선택해야 한다.

문제 4 — 최고 진폭과 준비 에너지

진폭 상한을 고정한 채 모드 대역폭을 넓히면 전체 에너지도 고정되는가?

해설: 아니다. 정규화는 $\sum A_n$을 고정하지만 에너지는 $\sum\omega_n^2A_n^2$에 비례한다. 진폭 배분과 진동수 가중치가 모두 바뀔 수 있다. 전달을 비교하기 전에 실제 준비 에너지를 측정하고, 하나의 정규화가 모든 물리량을 보존한다고 가정하지 않는다.

문제 5 — 주기 경계 근처의 중심

길이8 m인 고리의0과6 m에만 같은 가중치가 있으면 원형 중심은 어디인가?

해설: 두 위상은0과 $3\pi/2$이고 평균 편각은 $-\pi/4$다. 중심은8 m를 법으로7 m, 집중도는 $1/\sqrt2$다. 선형 평균3 m는 두 위치를 잇는 짧은 호에서 멀어지므로 경계 근처 묶음의 중심으로 적절하지 않다.

문제 6 — 중심 미정과 에너지 전달

진행 모드의 원형 에너지 중심이 미정이어도 전달 전력이 존재할 수 있는가?

해설: 그렇다. 순수 정현 모드의 배분 에너지는 상수와 두 배 모드 성분을 가지므로 첫 원형 모멘트가0일 수 있다. 에너지가 반드시 공간적으로 균일하다는 뜻은 아니다. 링크 전력의 전체 합은 여전히 양수일 수 있다. 중심 추적만으로 실제 흐름 정의를 대신할 수 없는 사례다.

자료와 재현 가능한 근거

MIT 8.03SC 강의13은 분산 파동의 도입을, Harvard Matthew Schwartz 강의11은 파동 묶음의 전개를 다룬다. 유한 사슬의 정규화·에너지 배분·국소 전력·전체 전달 항등식은 명시한 역학 방정식에서 직접 유도했다.

src/lib/learning-dispersion.ts의 latticeDispersivePacket이 실제 입자와 서로 다른 관측량을 계산한다. scripts/learning/generate-dispersion-references.py는 제품을 불러오지 않고 SciPy 행렬 지수함수와 FFT 복원을 사용한다. 검사는 실제 출력·힘·국소 수지·보존·병진·미정 관측량을 대조한다. 이 근거가 모든 실제 재료 근사나 아직 검사하지 않은 홈페이지 그림을 인증하는 것은 아니다.