선형 Euler 음향과 세 특성

소리 신호는 유체와 같은 속도로만 움직이지 않는다

천천히 흐르는 공기 옆에서 말을 한다고 생각해봅시다. 공기는 물질을 운반하지만, 압력의 변화는 공기에 대해서도 전파됩니다. 실험실의 관찰자가 보면 흐름을 따라가는 정보, 흐름보다 빠르게 내려가는 정보, 흐름에 거슬러 올라가는 정보가 구분됩니다. 신경망이 다음 유동 프레임을 예측할 때도 중요합니다. 매끄러운 영상을 만들더라도 압력이 잘못된 속도로 움직이면 물리적 예측은 틀릴 수 있습니다.

이 단원은 손으로 미분할 수 있는 작은 섭동의 매끄러운 예를 만듭니다. 선형 방정식의 정확한 해와, 그 방정식을 비선형 유체의 근사로 사용하는 일을 구분합니다. 학습 손실이 작거나 애니메이션이 자연스럽다는 사실만으로 이 구분이 없어지지는 않습니다.

선수 지식과 물리량

편미분, 질량 보존, 진행파와 행렬의 고유값을 복습합니다. 체적력과 열전달이 없는 일차원 비점성 이상기체를 생각합니다. 균일한 배경의 밀도, 속도, 압력을 각각 rho-zero, U, p-zero로 표시합니다. 비열비 gamma가 압력 방정식을 닫습니다. 음속은 양의 배경 물성이며 물질의 이동 속도와는 다릅니다. 아래의 프라임은 미분이 아니라 섭동입니다. 미분은 별도의 아래첨자로 씁니다.

물리량	단위	의미
밀도와 밀도 섭동	kg/m³	실제 체적당 질량
U, 음속, 속도 섭동	m/s	구분해서 읽어야 하는 속도
압력과 압력 섭동	Pa	면적당 힘
파수	rad/m	공간에 따른 위상 변화율

$$\rho=\rho_0+\rho',\quad v=U+u',\quad p=p_0+p',\qquad c^2=\frac{\gamma p_0}{\rho_0},\quad D_0=\partial_t+U\partial_x.$$

일차 항을 남겨 선형계를 유도하기

연속 방정식, 원시변수 운동량 방정식과 이상기체 압력 방정식에서 출발합니다. 각 변수에 배경과 섭동의 합을 대입합니다. 균일한 배경의 미분은 사라집니다. 연속 방정식의 곱을 펼치면 배경 밀도에 속도 섭동의 기울기를 곱한 항과, 섭동끼리 곱한 항이 함께 생깁니다. 앞의 일차 항은 남기고 뒤의 이차 항은 생략합니다. 진폭이 작다는 근사이지, 불편한 미분을 마음대로 없애도 된다는 뜻이 아닙니다. 짧은 파장은 큰 기울기를 만들 수 있으므로 차원이 있는 잔차를 읽을 때도 주의합니다.

$$\rho_t+v\rho_x+\rho v_x=0,\qquad v_t+v v_x+\frac{p_x}{\rho}=0,\qquad p_t+v p_x+\gamma p v_x=0.$$

아래의 선형 방정식들은 서로 결합되어 있습니다. 압력 기울기는 속도를 바꾸고 속도 발산은 밀도와 압력을 바꿉니다. 이 연결 중 하나를 지우면 올바른 음향 전파가 나오지 않습니다.

$$D_0\rho'+\rho_0u'_x=0,\qquad D_0u'+\frac{p'_x}{\rho_0}=0,\qquad D_0p'+\rho_0c^2u'_x=0.$$

변수 대신 특성 정보를 추적하기

속도 방정식에, 음향 임피던스로 나눈 압력 방정식을 더합니다. 빼면 반대쪽 음향 조합을 얻습니다. 압력 방정식에서 음속 제곱을 곱한 연속 방정식을 빼면 세 번째 조합이 생깁니다. 각 조합은 이제 하나의 수송 방정식을 따릅니다. 이것이 플럭스 행렬의 대각화를 실제로 읽는 방법입니다. 원래 물리량들은 결합되어 있지만, 적절한 조합으로 표현한 정보는 독립적으로 전파됩니다.

$$w_+=u'+\frac{p'}{\rho_0c},\quad w_-=u'-\frac{p'}{\rho_0c},\quad s=p'-c^2\rho'.$$

s의 단위는 압력입니다. 이 선형화에서 엔트로피와 관련된 접촉 변수이지만 J/(kg·K)로 측정하는 열역학적 엔트로피 자체는 아닙니다.

$$(\partial_t+(U+c)\partial_x)w_+=0,\quad(\partial_t+(U-c)\partial_x)w_-=0,\quad D_0s=0.$$

세 모드로 실제 장 만들기

오른쪽 음향, 왼쪽 음향, 접촉 섭동의 무차원 진폭을 각각 정합니다. 진폭의 부호도 허용합니다. 공통 위상은 실험을 간단히 만드는 선택이며 일반해에서 모든 위상이 같아야 한다는 뜻은 아닙니다. 파수는 양수입니다. 배경 유속을 바꾸면 세 특성 속도가 함께 이동하고, 음속을 바꾸면 특성 속도 사이의 간격이 달라집니다. 진폭의 부호를 바꾸는 것과 전파 방향을 바꾸는 것은 다릅니다.

$$A=a_+\cos(k[x-(U+c)t]+\phi),\quad B=a_-\cos(k[x-(U-c)t]+\phi),\quad C=a_0\cos(k[x-Ut]+\phi).$$

아래처럼 물리량을 복원하면 접촉 모드의 역할이 드러납니다. 접촉 모드는 압력과 속도를 바꾸지 않으면서 밀도만 바꿉니다. 음향 시각화의 유용한 대조 실험입니다.

$$u'=c(A+B),\qquad p'=\rho_0c^2(A-B),\qquad\rho'=\rho_0(A-B+C),\qquad s=-\rho_0c^2C.$$

시공간 그림을 정확히 읽기

인터랙티브 그림처럼 가로축에 시간, 세로축에 위치를 놓습니다. 표시한 파봉은 세 종류의 직선 중 하나를 따라갑니다. 이때 기울기는 단위 시간당 위치 변화, 즉 실험실에서의 특성 속도입니다. 여러 경로를 비교할 때도 이 축 관례를 유지하세요. 배경 유속이 영이면 밀도 접촉 신호가 정지할 수 있지만 압력파는 여전히 양쪽으로 움직입니다.

$$x_+(t)=x_+(0)+(U+c)t,\quad x_-(t)=x_-(0)+(U-c)t,\quad x_0(t)=x_0(0)+Ut.$$

이 선들을 입자 궤적이라고 부르면 안 됩니다. 물질 입자는 자신의 바뀌는 위치에서 평가한 전체 속도를 따라갑니다. 작은 진폭 음향에서도 다른 미분방정식입니다. 곡선이 매끄럽다는 모양만으로 무엇을 추적한 것인지는 알 수 없습니다.

단위를 포함한 손계산

배경 밀도를1kg/m³, 음속을2m/s, 배경 유속을0.5m/s로 두고 오른쪽 음향 진폭0.01만 켭니다. 파봉의 속도 섭동은0.02m/s, 압력 섭동은0.04Pa입니다. 이 특성 신호의 실험실 이동 속도는2.5m/s입니다. 계산을 읽기 쉽게 고른 교육용 기체의 척도이며, 상온 공기의 측정값을 재현하는 설정은 아닙니다. 신호 속도2.5와 국소 전체 유속0.52를 같은 것으로 표시하지 않습니다.

$$\rho_0=1,\ c=2,\ U=0.5,\ a_+=0.01:\quad u'=0.02\mathrm{\ m/s},\quad p'=0.04\mathrm{\ Pa},\quad\rho'=0.01\mathrm{\ kg/m^3}.$$

실제 유한 관측 창의 수지 확인하기

국소 방정식은 고정된 공간 구간의 수지도 예측합니다. 구간 내부의 시간 변화율을 적분하고, 오른쪽 끝에서 왼쪽 끝을 뺀 방향 있는 플럭스와 비교합니다. 임의의 구간은 정수 개의 파장을 포함하지 않을 수 있으므로 주기 상쇄를 가정하면 안 됩니다. 선형 밀도 플럭스에는 배경에 의한 섭동 밀도 수송뿐 아니라 속도 섭동이 배경 밀도를 운반하는 부분도 들어갑니다.

$$\frac{d}{dt}\int_a^b\rho'\,dx+[U\rho'+\rho_0u']_a^b=0.$$

이 일차원 적분의 단위는 면적당 질량입니다. 실제 질량을 얻으려면 명시한 단면적을 곱해야 합니다. 바로 kg이라고 표시하면 주어지지 않은 기하를 몰래 추가하는 셈입니다. 창 내부의 값이 변한다는 사실과 보존 위반을 혼동하지 않습니다.

이차 음향 에너지와 플럭스

선형계의 속도·압력 부분에는 양의 음향 에너지가 있습니다. 속도 방정식에 배경 밀도와 속도 섭동을 곱하고, 압력 방정식에 압력 섭동을 배경 밀도와 음속 제곱으로 나눈 값을 곱합니다. 둘을 더하고 곱의 미분법을 적용하면 아래 수지가 나옵니다. 배경 이류와 압력 일이 서로 다른 플럭스 기여를 만듭니다. 접촉 모드는 이 에너지 정의에 기여하지 않지만, 그 밀도 변화까지 사라졌다는 뜻은 아닙니다.

$$q=\frac{\rho_0u'^2}{2}+\frac{p'^2}{2\rho_0c^2},\qquad q_t+\partial_x(Uq+p'u')=0,\qquad[q]=\mathrm{J/m^3}.$$

이 양은 내부 에너지와 배경 운동까지 포함하는 비선형 기체의 전체 에너지와 구별해야 합니다. 서로 다른 정의의 에너지 곡선을 겹쳐 놓고 수지가 일치해야 한다고 요구하면 안 됩니다. 그림과 표에 먼저 정의를 표시합니다.

같은 장을 비선형 Euler에 대입하기

복원한 전체 물리량을 처음의 비선형 방정식에 넣습니다. 선형 잔차는 반올림 크기에 가까워야 하지만 비선형 잔차는 일반적으로 영이 아닙니다. 연속 방정식과 압력 방정식에는 섭동끼리의 곱이 남습니다. 원시변수 속도 방정식에는 전체 밀도로 나누는 연산이 있으므로, 유한한 모든 진폭에서 정확히 네 배 관계를 갖는 것이 아니라 작은 진폭 극한에서 이차 크기에 접근합니다.

$$R_\rho^{NL}=u'\rho'_x+\rho'u'_x,\qquad R_p^{NL}=u'p'_x+\gamma p'u'_x,\qquad R_v^{NL}=u'u'_x+\left(\frac1{\rho_0+\rho'}-\frac1{\rho_0}\right)p'_x.$$

진폭을 모두 절반으로 줄이면 앞의 두 잔차는 사분의 일이 됩니다. 파장을 바꾸면 미분의 척도도 변합니다. 이 비교는 선택한 해 가족의 진단이지 임의의 Euler 해에 대한 전역 오차 상한이 아닙니다.

양성, 부동소수점 관측과 한계

작은 진폭은 밀도와 압력의 양성을 보존합니다. 교육용 범위에서 각 무차원 진폭의 절댓값은0.03 이하이고 gamma는1.1부터1.67입니다. 배경 밀도와 음속은 양수로 유지합니다. 모드 크기의 합을 사용하면 모든 위치를 일일이 표본화하지 않고도 양성의 전역 조건을 확인할 수 있습니다. 다만 전체값이 양수라는 조건만으로 모든 실험에서 선형 근사가 충분히 정확하다고 결론 내리지는 않습니다.

$$\frac{|\rho'|}{\rho_0}\le |a_+|+|a_-|+|a_0|,\qquad\frac{|p'|}{p_0}\le\gamma(|a_+|+|a_-|).$$

아주 작은 접촉 변수는 압력 척도의 반올림된 두 관측을 빼면 사라질 수 있습니다. 접촉 모드에서 직접 얻은 안정적인 식과 이 뺄셈 결과를 함께 남깁니다. 속도 섭동은 표현 가능해도 그 제곱은 언더플로할 수 있습니다. 계산 불가한 이차 관측을 물리적으로 섭동이 없다는 증거로 표시하면 안 됩니다. 이 매끄러운 예에는 충격파 형성, 엔트로피 생성이나 벽 경계층 계산이 포함되지 않습니다.

가용성은 각 관측량에 따로 붙습니다. 예를 들어 아주 작은 모드에 배경 밀도를 곱하면 밀도 섭동이 표시 범위 아래일 수 있습니다. 그 정확한 기여를 양의 배경에 더한 전체 밀도는 여전히 표현 가능합니다. 압력과 비영 배경 유속에서도 같은 구분이 필요합니다. 선언한 모드의 기여를 먼저 결합한 뒤 전체값·미분·이차 에너지를 각각 반올림합니다. 따라서 섭동의 표시 불가가 전체 물리량의 표시 불가로 자동 전파되지는 않습니다. 유한 창 구적은 표본 관측을 사용하는 별도 근사이므로 자체 산술 한계가 남을 수 있습니다.

실험 순서와 AI 연구 연결

먼저 음향 모드 하나만 켜고 파봉의 변위와 특성 속도를 비교합니다. 반대 음향 모드를 추가해도 각 전파 속도는 그대로이며 간섭 무늬가 생깁니다. 접촉 모드만 켜서 압력 신호 없이 밀도가 이동하는지 확인합니다. 배경 유속을 바꾸어 왼쪽 음향 신호가 실험실에서 정지하는 지점을 찾습니다. 마지막으로 같은 시공간 좌표에서 모든 진폭을 절반으로 줄여 선형·비선형 잔차를 비교합니다.

물리 기반 신경망에서는 세 방정식 잔차의 단위와 정규화를 따로 기록합니다. 차원이 다른 잔차의 크기를 그대로 더한 점수는 물리적으로 불변인 척도가 아닙니다. 학습에 사용하지 않은 배경 유속이나 음속을 남겨 전파 속도를 시험합니다. 제한된 위상만 학습한 모델은 점별 보간을 잘해도 수송을 틀릴 수 있고, 압력 없는 접촉 신호를 누락된 자료로 오해할 수도 있습니다. 특성 이동, 유한 창 수지, 대조 실험을 각각 관측합니다. 이 검사들이 일반 압축성 유동 모델 전체를 인증하는 것은 아닙니다.

해설이 있는 연습

연습 1: 실험실에서 정지한 신호

배경 유속이 음속과 같을 때 어느 특성의 실험실 속도가 영인가요?

해설: 마이너스 음향 특성의 속도 U−c가 영입니다. 플러스 특성은 음속의 두 배로, 접촉 변수는 U로 이동합니다. 음향 정보가 정지했다고 유체도 정지했다고 결론 내리면 안 됩니다.

연습 2: 압력 없는 교란

두 음향 진폭을 영으로 두고 접촉 진폭만 남겼을 때 세 섭동을 설명하세요.

해설: 속도와 압력 섭동은 영입니다. 밀도 섭동은 배경 밀도에 접촉 모드를 곱한 값이고 U로 이동합니다. 정의한 음향 에너지 q는 영이지만 밀도장은 여전히 불균일합니다.

연습 3: 임피던스와 단위

오른쪽 음향 모드 하나에서, 공통 영점을 제외하고 압력 섭동을 속도 섭동으로 나누세요.

해설: 배경 밀도와 음속의 곱입니다. 단위는 kg/(m²·s), 즉 Pa를 m/s로 나눈 것입니다. 무차원 상관계수가 아니라 음향 임피던스입니다. 공통 영점에서는0/0을 계산하지 말고 모드의 관계식을 읽습니다.

연습 4: 진폭 세분화

배경과 파수를 고정하고 세 진폭을 모두 절반으로 줄이면 비선형 질량·압력 잔차는 어떻게 되나요?

해설: 남은 모든 항이 두 섭동 인자의 곱이므로 사분의 일이 됩니다. 선형 잔차는 산술 잡음 수준에 남습니다. 파수까지 함께 바꾼 실험에서는 이 비율을 그대로 요구할 수 없습니다.

연습 5: 유한한 구간

고정된 창 내부의 밀도 섭동 적분이 증가한다면, 오른쪽에서 왼쪽을 뺀 선형 질량 플럭스의 부호는 무엇인가요?

해설: 음수입니다. 나가는 것보다 들어오는 플럭스가 큽니다. 유한 창 내부의 변화는 보존과 양립합니다. 경계 기여를 확인하지 않고 영으로 놓으면 존재하지 않는 보존 위반을 만들게 됩니다.

연습 6: 잘못 읽은 에너지 관측

속도 섭동은 비영으로 표현되지만 음향 에너지를 계산할 때 제곱이 언더플로했습니다. 음향 교란이 없다고 결론 내릴 수 있나요?

해설: 아닙니다. 비영의 장 관측을 유지하고, 비영 곱을 표현하지 못하는 에너지 관측은 계산 불가로 남깁니다. 산술 범위와 물리적 부재를 구분하며 음향 에너지와 비선형 기체 전체 에너지의 정의도 분리합니다.

자료와 적용 범위

NASA Glenn의 Euler 방정식 소개는 물리량의 결합과 비점성 근사를 설명합니다. 위의 특성 복원과 진단은 명시한 균일 이상기체 선형화에서 직접 유도했습니다. NASA에 표시된 정상 이차원 방정식을 이 비정상 세 방정식 유도의 출처로 바꾸어 인용하지 않습니다. 다음 Taylor–Green 단원은 다른 가정을 사용하는 비압축성 비선형 해 가족을 다룹니다.