압축 흐름의 운동량·에너지 수지

같은 입자를 따라가면 에너지도 일정할까?

같은 기체 입자들을 따라가는 상자를 생각해 봅시다. 입자가 경계를 통과하지 않으므로 질량은 일정합니다. 그렇다면 그 안의 총에너지도 반드시 일정할까요? 팽창하는 기체는 바깥에 일을 할 수 있습니다. 질량 유출0을 에너지 교환0으로 바꾸면 이 압력 일을 놓칩니다. 이 단원에서는 각 입자의 속도가 일정한 매끄러운 압축·팽창 흐름을 골라 질량, 운동량, 내부·운동 에너지의 예산을 따로 계산합니다.

선수 지식과 모형의 범위

이동 제어체적의 상대 유속, Newton 법칙, 운동 에너지, 압력·일을 복습하세요. 일정 단면적 A의1차원 이상기체를 가정합니다. 점성, 열전도, 체적력, 화학 반응과 충격파는 포함하지 않습니다. γ는1보다 큰 일정한 비열비이며 내부 에너지 밀도와 압력을 연결합니다. 상태식의 더 일반적인 유도와 적용 범위는 열역학 과정에서 확장합니다. 지금은 지정한 매끄러운 해 가족이 세 보존식을 어떻게 만족하는지 확인합니다.

$$\mathcal I=\rho e=\frac{p}{\gamma-1},\qquad \mathcal K=\frac12\rho u^2,\qquad \mathcal E=\mathcal I+\mathcal K,\qquad \gamma>1.$$

세 가지 저장량과 단위

보존 벡터의 첫 성분은 질량 밀도, 둘째는 운동량 밀도, 셋째는 총에너지 밀도입니다. 둘째의 부호는 속도 방향을 반영하며 셋째의 운동 에너지는 속도의 제곱을 포함합니다. 밀도들을 공간 부피로 적분하면 실제 저장량이 됩니다. 에너지를 압력이나 단위 질량당 에너지와 혼동하지 마세요. 이들은 연결되어 있지만 수치와 단위가 같은 양은 아닙니다.

$$U=(\rho,\rho u,\mathcal E),\qquad [U]=(\mathrm{kg/m^3},\mathrm{kg/(m^2s)},\mathrm{J/m^3}).$$

압력은 힘과 일을 운반한다

질량 유속은ρu입니다. 운동량 유속에는 물질이 운반하는ρu²와 압력 p가 함께 들어갑니다. 총에너지 유속은 저장 밀도에 속도만 곱한 값이 아닙니다. 압력에 속도를 곱한 일률 p u가 추가됩니다. 따라서 에너지 저장은 내부·운동 에너지의 합이지만 에너지 유속에는 엔탈피에 해당하는 항이 나타납니다. 저장량의 정의와 유속의 정의를 섞지 마세요.

$$F=(\rho u,\rho u^2+p,(\mathcal E+p)u),\qquad \partial_tU+\partial_xF=0.$$

이동 경계에서는 F−wU

경계가 속도 w로 움직이면 그 경계를 통과하는 상대 유속은 F−wU입니다. 질량에는ρ(u−w)가 남고 운동량에는ρu(u−w)+p가 남습니다. 에너지는 E(u−w)+p u입니다. 마지막 항을 p(u−w)로 바꾸면 안 됩니다. 물질 경계 w=u에서도 압력 일 p u는 남습니다. 왼쪽·오른쪽 바깥 법선을 적용한 뒤 양쪽 유속 차이가 저장 변화율을 결정합니다.

$$F-wU=\bigl(\rho(u-w),\ \rho u(u-w)+p,\ \mathcal E(u-w)+pu\bigr),\qquad \frac{d}{dt}\int_{x_L}^{x_R}AU\,dx=A[(F-wU)_L-(F-wU)_R].$$

매끄러운 affine 압축·팽창

초기 위치에 선형인 속도를 줍니다. 변형 척도 s=1+at가 양수인 동안 밀도는1/s, 압력은1/s^γ로 변합니다. a가 양수이면 팽창, 음수이면 압축입니다. v0는 전체 평행 이동을 추가합니다. 제품은 a와 시간의 범위를 제한해 s가0에 접근하지 않게 합니다. s=0의 붕괴나 진공을 작은 양수로 대체하여 이 해를 계속하지 않습니다.

$$s(t)=1+at>0,\qquad \rho=\frac{\rho_0}{s},\qquad u=\frac{v_0+ax}{s},\qquad p=p_0s^{-\gamma}.$$

입자 속도는 왜 일정한가?

초기 위치 X의 입자는 x=sX+v0t로 움직입니다. 이를 한 번 미분하면 v0+aX이며 시간에 무관합니다. 현재 위치를 속도장에 대입해도 같은 값입니다. 압력의 공간 기울기가0이므로 이 가족의 물질 가속도도0입니다. 공간의 한 지점을 고정하고 측정한 속도는 시간에 따라 바뀔 수 있습니다. 고정 위치 미분과 같은 입자를 따라가는 미분을 구분하세요.

$$x(X,t)=sX+v_0t,\qquad \frac{Dx}{Dt}=v_0+aX,\qquad \partial_tu+u\partial_xu=0,\qquad \partial_xp=0.$$

연속방정식과 운동량을 대조하기

속도 기울기 H=a/s는 공간적으로 일정합니다. 밀도는 시간에 따라−Hρ만큼 변하며 질량 유속의 공간 변화가 이를 보상합니다. 운동량 밀도의 시간 미분에는−2H가 붙습니다. 이 항은ρu²+p의 공간 미분과 상쇄합니다. 단순히 ρ만 시간에 따라 바꾸고 속도나 압력을 그대로 두면 같은 해 가족이 아닙니다. 실제 장의 변화율을 함께 대조해야 합니다.

$$H=\frac as,\qquad \partial_t\rho=-H\rho,\qquad \partial_x(\rho u)=H\rho,\qquad \partial_t(\rho u)=-2H\rho u,\qquad \partial_x(\rho u^2+p)=2H\rho u.$$

내부 에너지와 운동 에너지의 국소 변화

고정 위치에서 내부 에너지 밀도는−γH I, 운동 에너지 밀도는−3H K만큼 변합니다. 두 값을 더하면 총에너지의 시간 미분입니다. 에너지 유속의 공간 미분이 이를 보상합니다. 같은 입자의 속도가 일정해도 고정 부피의 운동 에너지 밀도는 달라질 수 있습니다. 그 부피에 들어 있는 질량과 입자 구성이 바뀌기 때문입니다.

$$\partial_t\mathcal I=-\gamma H\mathcal I,\quad \partial_t\mathcal K=-3H\mathcal K,\quad \partial_t\mathcal E=-H(\gamma\mathcal I+3\mathcal K),\quad \partial_x[(\mathcal E+p)u]=H(\gamma\mathcal I+3\mathcal K).$$

고정·평행 이동·물질 구간

초기 중심 c와 반폭 h를 정한 뒤 세 관측을 비교합니다. 고정 구간은 같은 위치를, 평행 이동 구간은 일정한 관측 속도를, 물질 구간은 같은 입자를 따라갑니다. affine 흐름의 물질 구간은 폭도 변하므로 단순 평행 이동으로 대신할 수 없습니다. 모드에서 사용하지 않는 속도 입력은 실제 경계를 바꾸지 않습니다. 현재 양쪽 좌표와 각 경계 속도를 함께 확인하세요.

$$x_{L,R}^{\rm fixed}=c\mp h,\qquad x_{L,R}^{\rm translating}=c+wt\mp h,\qquad x_{L,R}^{\rm material}=s(c\mp h)+v_0t.$$

물질 구간의 질량·운동량·에너지

물질 구간의 질량은 초기 길이·단면적·밀도의 곱이며 일정합니다. 양쪽 압력이 같으므로 전체 운동량 변화도0입니다. 각 입자의 속도가 일정하므로 이 구간의 운동 에너지도 일정합니다. 하지만 팽창하면 내부 에너지가 감소합니다. 양쪽 압력 힘의 합이0이어도 서로 다른 속도로 움직이는 경계에서 하는 일의 합은0이 아닐 수 있습니다.

$$M_{\rm mat}=2Ah\rho_0,\quad \dot P_{\rm mat}=0,\quad \dot K_{\rm mat}=0,\quad I_{\rm mat}=\frac{2Ahp_0}{\gamma-1}s^{1-\gamma},\quad \dot E_{\rm mat}=-2Aahp.$$

손계산: 열교환 없이2/3J 감소

A=1m², ρ0=1kg/m³, p0=1Pa, a=0.5/s, v0=0, γ=2를 고릅니다. 초기 물질 구간은[−1,1]m입니다. t=0에서 질량은2kg, 내부 에너지는2J, 운동 에너지는1/12J입니다. t=1초에는 s=1.5, 구간은[−1.5,1.5]m, 압력은4/9Pa입니다. 내부 에너지는4/3J로 줄고 운동 에너지는 그대로여서 총에너지가2/3J 감소합니다. 그만큼 경계에서 바깥으로 압력 일을 했습니다.

$$E(0)=2+\frac1{12}=\frac{25}{12}\ \mathrm J,\qquad E(1)=\frac43+\frac1{12}=\frac{17}{12}\ \mathrm J,\qquad \Delta E=-\frac23\ \mathrm J.$$

누락한 압력 일을 음성 대조로 사용하기

같은 실제 장에서 에너지 유속의 p u만 빼 봅시다. 물질 경계에서는 대류 항도0이므로 잘못된 모형은 에너지 교환이 없다고 예측합니다. 실제 저장 변화율은 여전히−2Aahp입니다. 이 차이를 숨기려고 원천이나 장을 바꾸지 않습니다. 국소 검사에서도∂tE+∂x(Eu)=−pH가 남습니다. 손계산 예제의 최초 순간에는 올바른 순유입−1W, 잘못된 순유입0W입니다.

$$F_{E,\rm wrong}=\mathcal Eu,\qquad \partial_t\mathcal E+\partial_xF_{E,\rm wrong}=-pH,\qquad \dot E_{\rm mat}(0)=-1\ \mathrm W.$$

정확한 영·작은 변화·유효 범위

a=0이면 균일한 평행 이동 또는 정적 압력 상태가 됩니다. p0=0이면 양의 밀도를 가진 무압력 극한이며 진공과 다릅니다. 영폭 구간은 정상적인 면 관측과 구분합니다. 작은 a에서 두 끝점의 에너지가 같은 수로 반올림되어도 내부 에너지 변화는 남을 수 있습니다. log1p와 expm1으로 작은 비율 변화를 따로 계산합니다. 사용할 수 없는 산술 값과 정확한 영을 구분하고, 표시한 잔차를 임의 기하·충격파의 정확성 보증으로 확장하지 않습니다.

$$\Delta I_{\rm mat}=I_{\rm mat}(t_0)\operatorname{expm1}\!\left((1-\gamma)\operatorname{log1p}\!\frac{a(t_1-t_0)}{1+at_0}\right).$$

log1p와 expm1만 사용해도 무차원 인자가 이미0으로 반올림되었다면 작은 변화를 잃습니다. 계산은 그 인자를 분리하고 log1p(z)/z와 expm1(z)/z의 연속 극한1을 사용한 뒤, 나머지 인자를 에너지 척도와 결합하여 마지막에 표시 범위로 반올림합니다. 예를 들어 초기 내부 에너지2J, γ=1.5에서 팽창률을 binary64의 최소 양수/s로 두면1초 뒤 에너지 변화는 여전히 비영인 음수로 반올림됩니다. 두 끝점 에너지는 모두2J로 표시됩니다. 이는 선언한 모형 안의 산술 차이이며 실험 정밀도를 주장하는 결과가 아닙니다.

AI 모델의 세 예산을 검증하는 순서

먼저 물질 팽창 예제에서 질량·운동량·운동 에너지는 일정하고 내부 에너지는 줄어드는지 확인하세요. 압축으로 부호를 반전하고, 고정 구간과 평행 이동 구간에서 상대 유출을 비교합니다. 다음에는 압력 일을 누락한 대조, 정확한 a=0, 작은 비영 a를 검사합니다. 학습 모형에는 보존식뿐 아니라 상태 관계, 초기·경계·관측 데이터가 필요합니다. 각 예산의 단위·척도·허용 오차를 따로 보고하고, 같은 해석식의 항등식과 독립 구적·차분·입자 ODE 결과를 구분하세요.

풀이가 있는 연습 문제

연습 1

물질 경계 w=u이면 세 상대 유속은 모두0인가요?

해설: 질량 유속은0이지만 운동량 유속에는 p, 에너지 유속에는 p u가 남습니다.

연습 2

손계산 예제의 최초 왼쪽·오른쪽 운동량 바깥 유출과 순힘은?

해설: 왼쪽−1N, 오른쪽+1N이므로 순운동량 변화는0입니다. 두 힘의 일률은 경계 속도가 달라 상쇄하지 않습니다.

연습 3

t=1초의 총에너지 변화와 그 원인을 설명하세요.

해설: 25/12J에서17/12J로 변해−2/3J입니다. 질량 유출이나 열전도가 아니라 팽창 경계의 압력 일입니다.

연습 4

a의 부호를 음수로 바꾸면 물질 내부 에너지 변화의 부호는?

해설: s>0 범위에서 압축이 일어나 내부 에너지가 증가합니다. 각 입자 속도는 계속 일정하므로 물질 운동 에너지는 변하지 않습니다.

연습 5

a=0, v0=0, p0>0이면 모든 에너지가0인가요?

해설: 아닙니다. 정적 압력과 내부 에너지는 남습니다. 에너지 유속과 시간 변화는0입니다.

연습 6

세 보존 잔차가0인 학습 모형이면 해가 유일하게 식별되나요?

해설: 보장하지 않습니다. 상태 관계와 초기·경계·관측 정보가 필요하고, 지정 가족의 검사를 일반 해의 인증으로 바꾸면 안 됩니다.

자료와 다음 연결

이동 경계와 물질 미분의 일반 배경은 MIT의 수송 정리 강의, 응력·일과 적분형 보존식은 Stanford의 보존식 강의를 참고하세요. 이 단원의 해 가족과 손계산은 교재의 예제입니다. 다음에는 연속체 가정, 응력·변형, 일반 Euler·Navier–Stokes, 구성 법칙과 보존형 이산화로 확장합니다.