초기 파형·반사·국소 에너지 흐름

파동은 초기값 문제다

줄의 모양을 한 장 찍은 사진만으로 미래를 정할 수 없다. 같은 모양도 순간적으로 정지해 있거나 왼쪽·오른쪽으로 이동하거나 별도의 초기 충격을 받은 상태일 수 있다. 파동 방정식에는 초기 변위와 초기속도가 모두 필요하다. 특성선이 끝점에 도달하면 경계에서 제공하는 정보도 필요해진다.

이 단원은 국소적인 매끄러운 파형을 따라간다. 변위·속도·공간 기울기·에너지 밀도·부호 있는 에너지 흐름을 계산한다. 움직이는 모양과 실제 물질의 운동, 에너지의 이동을 구분하는 것이 목표다. 앞의 주기 사슬에는 반사 끝점이 없었다. 여기서는 반직선의 고정단 또는 자유단을 선택할 수 있다.

매질과 방정식 및 물리적 범위

일정한 선밀도 $\mu>0$의 단위는 kg/m, 장력 $\tau>0$의 단위는 N이다.

$$u_{tt}=c^2u_{xx},\qquad c=\sqrt{\tau/\mu}.$$

변위 $u$는 m, 속도 $u_t$는 m/s, 기울기 $u_x$는 무차원이다. 작은 횡기울기를 가진 균일하고 감쇠 없는 선형 줄을 가정한다. 분포 외력과 시간 의존 경계 구동은 포함하지 않는다. 유한 구간에서의 여러 번 반사·Robin 조건·공간적으로 변하는 매질도 이 실험의 범위 밖이다.

물리 영역은 전체 직선 $x\in\mathbb R$ 또는 반직선 $x\ge0$이다. 유한한 그림 창은 관측 구간이지 두 번째 물리 경계가 아니다. 창을 벗어나는 파형은 그 위치에서 자동으로 반사되지 않는다.

두 개의 초기 장이 필요한 이유

초기조건은

$$u(x,0)=f(x),\qquad u_t(x,0)=g(x)$$

로 지정한다. $g=0$이면 국소 변위가 두 방향 파형으로 나뉜다. $g=-cf'$이면 같은 변위가 오른쪽으로 이동하고, $g=cf'$이면 왼쪽으로 이동한다. 처음의 변위 사진이 같더라도 서로 다른 물리적 준비다.

$c$를 곱해야 $g$가 올바른 속도 단위를 가진다. 그래프를 어느 방향으로 움직여 그렸는지만으로 이 초기조건을 대신할 수 없다. 이 실험에서는 독립적인 국소 속도 bump도 더하므로 변위와 속도가 반드시 하나의 진행파 관계를 가질 필요가 없다.

d’Alembert 공식을 유도하기

전체 직선의 일반해를 $u=R(x-ct)+L(x+ct)$라 쓰자. $t=0$에서

$$R+L=f,\qquad -cR'+cL'=g$$

다. $H'=g$인 원시함수를 택하면 $R=f/2-H/(2c)$, $L=f/2+H/(2c)$로 둘 수 있다. 원시함수의 임의 상수는 두 성분을 반대 방향으로 이동시키며 실제 합에서는 사라진다. 따라서

$$u(x,t)=\frac{f(x-ct)+f(x+ct)}2+ \frac1{2c}\int_{x-ct}^{x+ct}g(s)\,ds$$

를 얻는다. 적분을 생략하면 초기속도를0으로 바꾼 셈이다. $t=0$에서 적분 구간은 한 점으로 줄지만, 움직이는 끝점을 시간 미분하면 $g(x)$가 복원된다. 적분값이 처음에0이라는 것과 그 시간 미분이0이라는 것은 다르다.

실제 속도와 모양의 기울기는 다르다

$\ell=x-ct$, $r=x+ct$라 쓰고 직접 미분하면

$$u_t=\frac c2[f'(r)-f'(\ell)]+\frac12[g(r)+g(\ell)],$$ $$u_x=\frac12[f'(\ell)+f'(r)]+\frac1{2c}[g(r)-g(\ell)]$$

다. 앞의 식은 실제 물질점의 횡속도이고 뒤의 식은 모양의 기울기다. 오른쪽 진행파에서는 $u_t=-cu_x$다. 모든 입자가 일정한 물질 속도 $c$를 가진다는 뜻이 아니다. $c$는 정보와 모양의 수평 전파 속도다. 표시한 물질점은 평형 가로 좌표를 유지하면서 횡방향으로 운동한다.

매끄러움을 통제한 유한 지지 파형

다음 함수를 사용하자.

$$b(z)=\begin{cases}(1-z^2)^3,&|z|<1,\\0,&|z|\ge1,\end{cases} \qquad f(x)=A\,b\!\left(\frac{x-x_0}{w}\right).$$

$A$는 부호 있는 변위 진폭(m), $x_0$는 중심, $w>0$는 반폭(m)이다. 함수와 일·이차 도함수는 두 지지 경계에서 모두0으로 이어진다. 내부에서는

$$b'(z)=-6z(1-z^2)^2,\qquad b''(z)=-6+36z^2-30z^4$$

다. 전체 적분은 $32/35$다. 유한한 지지 덕분에 유한 전파를 명확히 검사할 수 있다. 의존 영역 밖의0은 Gaussian 꼬리를 임의로 잘라서 만든0이 아니다. 이 함수는 $C^2$이며 지지 경계에서 모든 차수의 매끄러움이 있다고 가정해서는 안 된다.

분리·진행·충격 초기조건 만들기

실험의 초기속도는

$$g(x)=d\,cf'(x)+B\,b\!\left(\frac{x-x_0}{w}\right), \qquad d\in\{-1,0,1\}$$

로 둔다. $B$는 독립적인 속도 진폭(m/s)이다. $B=0$일 때 $d=-1$은 오른쪽 진행, $d=1$은 왼쪽 진행, $d=0$은 초기 변위의 양방향 분리다. $B\ne0$이면 방향 이름은 변위에 연결된 초기 준비를 뜻하며 모든 에너지가 그 방향으로만 흐른다는 보장이 아니다.

반직선에서는 $x_0>w$로 두어 두 초기 장이 끝점 근처에서0이 되게 한다. 계산기는 추가로0.25 m의 여유를 둔다. 이렇게 해야 초기값을 나중에 고쳐 맞추지 않고도 고정단·자유단의 호환성을 확보한다.

유한 전파와 뜻밖의 평탄부

$(x,t)$의 값은 초기 구간 $[x-ct,x+ct]$의 정보에 의존한다. 이것은 의존 구간이지 안의 모든 점이 같은 비율로 영향을 준다는 뜻이 아니다. 초기 지지에서 $ct$만큼 확장한 영역 바깥에서는 변위 항과 속도 적분이 모두0이다.

전체 직선에서 $A=0$, $B\ne0$으로 두면 “눈에 보이는 변위가 곧 그 위치의 에너지”라는 생각의 반례가 나온다. 적분 구간이 속도 bump 전체를 포함한 내부 평탄부에서는

$$u=\frac{Bw}{2c}\frac{32}{35},\qquad u_t=u_x=0$$

이다. 줄에 정적인 변위가 남지만 국소 운동·변형 에너지는 없다. 에너지는 양쪽 전면에서 이동한다. $B=0.2$ m/s, $w=1$ m, $c=1$ m/s이면 평탄부 변위는 약0.09143 m다.

반사 경계에서는 초기 장 둘 다 연장하기

반직선 고정단은 $u(0,t)=0$이다. $f,g$를 각각0에 대해 홀수로 연장한 뒤 전체 직선 공식을 적용하고 $x\ge0$만 취한다. 자유단은 $u_x(0,t)=0$이며 두 함수를 각각 짝수로 연장한다.

홀수 함수의 도함수는 짝수이고 짝수 함수의 도함수는 홀수다. 따라서 양의 반직선에서 $g=cf'$가 성립해도 $f$를 먼저 연장하고 그 전체 도함수에 $c$를 곱해 $g$로 쓰면 일반적으로 틀린다. 초기속도의 대칭이 바뀌기 때문이다. 원래 지정한 두 초기 장에 필요한 대칭을 각각 적용해야 한다.

음의 영역에 놓은 반사상은 계산 도구다. 반직선 밖의 반사상 에너지를 실제 계의 추가 에너지로 세지 않는다. 물리적 적분 구간과 수학적으로 연장한 영역을 구분한다.

반사 부호와 도착 시각

$B=0$인 왼쪽 입사 파형의 반직선 해는

$$u(x,t)=f(x+ct)+\rho f(ct-x),\qquad x\ge0$$

로 쓸 수 있다. 고정단에서는 $\rho=-1$, 자유단에서는 $\rho=+1$이다. 중심은 $x_0/c$에 끝점에 도착하지만 앞부분은 $(x_0-w)/c$부터 먼저 도착한다. 유한 폭의 파형이 반사되는 과정 전체를 한순간에 뒤집는 애니메이션으로 대체하면 안 된다.

$A=0.1$ m, $x_0=2$ m, $w=1$ m, $c=1$ m/s이면 반사한 중심이 $t=4$ s에 다시 $x=2$에 있다. 고정단 변위는−0.1 m, 자유단은+0.1 m다. 두 반사파 모두 에너지를 오른쪽으로 운반한다. 변위의 부호가 에너지 이동 방향을 정하지 않는다.

실제 에너지와 국소 보존을 유도하기

평형 길이당 에너지와 흐름은

$$e=\frac12\mu u_t^2+\frac12\tau u_x^2, \qquad J=-\tau u_tu_x$$

다. $e$의 단위는 J/m, $J$는 W다. PDE에 $\mu u_t$를 곱하고 $\mu c^2=\tau$를 사용하면

$$\partial_t e =\tau(u_tu_{xx}+u_xu_{xt}) =\partial_x(\tau u_tu_x), \qquad e_t+J_x=0$$

를 얻는다. 이것은 국소 보존 법칙이다. 단지 전체 합이 일정하다는 진술보다 강하게, 구간 경계를 가로질러 무엇이 이동하는지 알려 준다. 고정단에서는 $u_t(0,t)=0$, 자유단에서는 $u_x(0,t)=0$이므로 이상적인 정지 경계의 전력 흐름은 모두0이다.

방향별 에너지와 흐름의 부호

오른쪽과 왼쪽의 음이 아닌 에너지 밀도를

$$e_R=\frac\mu4(u_t-cu_x)^2,\qquad e_L=\frac\mu4(u_t+cu_x)^2$$

로 정의하면 $e=e_R+e_L$, $J=c(e_R-e_L)$다. 순수 오른쪽 파형은 $J=ce$, 왼쪽은 $J=-ce$를 만족한다. 겹침에서 두 방향 성분이 같으면 에너지는 있어도 순흐름은0일 수 있다.

원시함수를 사용한 변위 성분 $R,L$에는 서로 반대인 임의 상수를 더할 수 있다. 방향별 에너지는 도함수로 정의되므로 이 선택에 영향받지 않는다. 성분 변위를 제곱하는 것만으로 에너지를 추정해서는 안 된다.

유한 관측 구간의 수지와 독립 적분

고정 관측 구간 $[a,b]$의 에너지는

$$E_{[a,b]}(t)=\int_a^b e(x,t)\,dx$$

이며 국소 법칙을 적분하면

$$E_{[a,b]}(t)-E_{[a,b]}(0) =\int_0^t[J(a,s)-J(b,s)]\,ds$$

다. 양의 $J(a)$는 왼쪽 유입, 양의 $J(b)$는 오른쪽 유출이다. 각각 부호 있는 적분이다. 계산기는 공간 에너지와 경계 전력을 별도로 적분한다. 에너지의 끝점 차를 흐름 적분으로 정의해 잔차를 강제로0으로 만들지 않는다.

유한 지지 다항식에서는 이동하는 지지 경계가 공간·시간 적분을 다항식 조각으로 나눈다. 각 조각의 에너지와 흐름 차수는 최대12다. 8점 Gauss–Legendre 공식은15차까지 정확하므로 모든 경계에서 구간을 나누면 실수 산술에서 정확하다. 부동소수점 오차는 남으며 작은 잔차를 숨기지 않는다. 그림용 표본과 이 적분의 계산점을 구분한다.

계산 예제와 물리적 한계

$A=0.1$, $B=0$, $w=1$, $x_0=2$, $\mu=\tau=1$인 오른쪽 파형을 $[0,4]$에서 관측하자. 처음의 지지 $[1,3]$은 모두 안에 있다. $t=4$에는 지지가 $[5,7]$이므로 관측 에너지는0이다. 오른쪽 경계 적분은 초기 에너지와 같고 왼쪽 전달은0이다. 전체 직선의 총에너지는 계속 보존된다.

순수 진행파의 총에너지는 $E=\tau A^2\int_{-1}^1[b'(z)]^2dz/w$다. 같은 진폭에서 폭을 좁히면 에너지와 기울기가 커진다. 초기 최대 기울기는 $96|A|/(25\sqrt5\,w)$다. 이는 초기 지표일 뿐이다. 추가 속도 장과 나중의 겹침 때문에 기울기가 달라질 수 있다. 지원하는 진폭·폭·시간의 수치 범위가 모든 설정의 물리적 작은 기울기를 보장하지는 않는다.

해설이 있는 연습 문제

문제 1: 모양만으로 방향을 정할 수 있을까

같은 $f$에 대해 한 줄은 $g=0$, 다른 줄은 $g=-cf'$다. 무엇이 다른가요?

해설: 앞의 줄은 반쪽 변위 파형 둘로 나뉘고 뒤의 줄은 그 분리 없이 오른쪽으로 진행한다. 초기속도 장이 사진에 없는 정보를 제공한다.

문제 2: 속도만 준 평탄부

$A=0$, $B=0.2$, $w=1$, $c=1$이고 적분 구간이 전체 bump를 포함할 때 내부 변위를 구하세요.

해설: $Bw(32/35)/(2c)\approx0.09143$ m다. 평탄부에서는 $u_t=u_x=0$이므로 변위가 있어도 국소 에너지는0이다.

문제 3: 대칭 연장의 함정

고정단에서 $f$의 홀수 연장을 미분하는 것만으로 왼쪽 입사파의 속도를 연장할 수 없는 이유는 무엇인가요?

해설: 그 도함수는 짝수지만 고정단의 초기속도는 홀수로 연장해야 한다. 양의 반직선에 주어진 원래 $g=cf'$를 별도로 홀수 연장한다.

문제 4: 반사 부호와 에너지

고정단이 음의 변위 파형을 돌려보내면 흐름도 음수인가요?

해설: 반사 뒤에는 오른쪽으로 이동하므로 $J=ce\ge0$이다. 진폭 부호를 바꾸면 속도와 기울기 부호가 함께 바뀌어 곱과 에너지 방향은 그대로다.

문제 5: 관측 구간에서 빠져나간 에너지

$[a,b]$의 에너지가0.03 J 감소했고 왼쪽 적분 유입이0.01 J라면 오른쪽 적분 유출은 얼마인가요?

해설: $-0.03=0.01-I_b$이므로 $I_b=0.04$ J다. 전체 보존은 국소 감소를 허용하며 양쪽 경계 항을 모두 고려해야 한다.

문제 6: 지지 경계에서 적분을 나누는 이유

다항식 구적법을 적용하기 전에 이동하는 지지 경계에서 구간을 나누는 이유는 무엇인가요?

해설: 적분 함수는 각 조각에서는 다항식이지만 경계를 가로지르는 하나의 다항식은 아니다. 정확 차수 보장은 조각에 적용된다. 경계를 무시하면 그럴듯해 보여도 틀린 보존 잔차를 얻을 수 있다.

AI 검증 질문과 일차 참고문헌

학습한 파동 대리 모델은 독립적인 두 초기 장을 복원하고, 유한 전파를 지키며, 물리 반사 경계와 잘라 보인 그림을 구분하고, 국소 에너지를 부호 있는 경계 흐름과 맞춰야 한다. 단순 진행파 외에 평탄부와 독립 초기속도가 있는 반사도 검사한다. 평균 변위 오차가 작다는 사실만으로 이 성질들이 검증되지는 않는다.

Herman의 d’Alembert 강의노트1–3.2절과 Toronto의 반직선 파동 강의노트는 초기값·경계의 틀을 제공한다. 유한 지지 다항식·SI 방향별 에너지·독립 관측 구간 수지·계산기 검사는 명시한 모델에서 직접 유도했다. 이 예제만으로 회절·일반 분산 전파·모든 경계값 문제의 집필을 완료한 것은 아니다.