이동 제어체적과 질량 수송

움직이는 사진 속 질량은 왜 달라질까?

강물의 한 구간을 사진으로 찍는다고 생각해 봅시다. 같은 자리에 카메라를 두면 새로운 물이 계속 들어오고 나갑니다. 드론이 흐름을 따라가면 사진 안에 같은 물을 더 오래 담을 수 있습니다. 사진 속 물의 양이 변했다고 해서 물이 생성된 것은 아닙니다. 먼저 무엇을 따라 움직이는 관측인지 정해야 보존 법칙을 올바르게 읽을 수 있습니다. 이 단원은 밀도의 이동과 관측 구간의 이동을 독립적으로 조절하는 매끄러운 예제로 그 차이를 계산합니다.

선수 개념과 실제 단위

미분·정적분, 속도, 밀도, 바깥 법선의 부호와 연속방정식을 복습하세요. 흐름은 x축 방향이며 단면적 A는 일정합니다. 시간은 실제 초, 위치와 반폭은 미터입니다. 밀도에 부피를 곱하면 질량이고 밀도에 속도·면적을 곱하면 초당 통과 질량입니다. 경계 질량 유량의 그래프와 저장 질량의 그래프는 단위가 다르므로 같은 원시 축에서 빼지 않습니다.

$$[\rho]=\mathrm{kg/m^3},\quad [A]=\mathrm{m^2},\quad [M]=\mathrm{kg},\quad [A\rho u]=\mathrm{kg/s}.$$

물질 집합과 제어체적

물질 집합은 처음 선택한 동일 입자들을 따라갑니다. 제어체적은 우리가 정한 공간의 창이며 유체와 다른 속도로 움직여도 됩니다. 유체의 속도를 v, 창의 속도를 w라고 합시다. 입자가 창을 통과하는 속도는 v 자체가 아니라 v−w입니다. v=w인 특수한 경우에만 이 예제의 평행 이동 창이 같은 입자 집합을 따라갑니다. 단순히 화면이 움직인다는 이유로 물질 집합이라고 부르지 마세요.

$$\dot x_{\mathrm{particle}}=v,\qquad \dot x_{\mathrm{window}}=w,\qquad u_{\mathrm{relative}}=v-w.$$

수송 정리는 경계의 이동을 센다

움직이는 영역 안의 적분을 미분하면 두 효과가 있습니다. 이미 영역 안에 있던 위치에서 장이 시간에 따라 바뀌는 효과와 경계가 이동해 다른 위치를 포함하는 효과입니다. 수송 정리는 이 기하학적 관계를 표현합니다. G는 단위 부피당 양이고 경계 속도는 w입니다. 이 식 자체는 물리 보존 법칙이 아닙니다. 어떤 양이 어떻게 수송되고 생성되는지에 대한 물리식과 결합해야 합니다.

$$\frac{d}{dt}\int_{V(t)}G\,dV=\int_{V(t)}\partial_tG\,dV+\int_{\partial V(t)}G\,\mathbf w\cdot\mathbf n\,dS.$$

왼쪽과 오른쪽의 부호

창의 왼쪽 경계를 xL, 오른쪽을 xR이라 합시다. 왼쪽 바깥 법선은−x, 오른쪽은+x입니다. 양의 x 방향 흐름이 왼쪽으로 들어오면 저장 질량을 늘리고, 오른쪽으로 나가면 줄입니다. 경계가 함께 움직이는 만큼을 빼면 다음 식이 됩니다. 두 경계의 속도가 서로 다른 일반 구간에도 같은 원리가 적용되지만, 이번 밀도 실험에서는 두 경계가 같은 w로 움직입니다.

$$\dot M=A[\rho(x_L,t)(v-w)-\rho(x_R,t)(v-w)],\quad I_L^{\rm out}=-A\rho_L(v-w),\quad I_R^{\rm out}=A\rho_R(v-w).$$

양의 밀도가 이동하는 해 가족

기준 밀도 ρ0에 코사인 대비 α를 더합니다. |α|≤0.8이므로 모든 위치의 밀도는 양수입니다. k가 양수이면 ρ0는 한 주기의 공간 평균입니다. k=0에서는 코사인이1이므로 실제 균일 밀도는 ρ0(1+α)입니다. ρ0라는 입력 이름을 모든 경우의 관측 평균으로 읽지 마세요. 무한한 x축에 양의 밀도가 있으므로 전체 축의 총질량은 유한하지 않습니다. 우리는 유한 창의 질량을 관측합니다.

$$\rho(x,t)=\rho_0[1+\alpha\cos(k(x-vt))],\qquad \rho_0>0,\quad |\alpha|\le0.8.$$

국소 연속방정식을 직접 확인한다

코사인의 공간 미분은 음의 사인이고 시간 미분에는 추가로−v가 붙습니다. 따라서 밀도의 시간 변화와 질량 유속의 공간 변화가 상쇄합니다. 이 관계를 손으로 계산하면 어떤 입력을 바꾸어도 이 해 가족이 연속방정식을 만족하는 이유를 볼 수 있습니다. 제품이 같은 해석식에서 만든 잔차0은 독립 검증이 아닙니다. 실제 밀도 표본의 차분과 구간 적분을 별도로 비교합니다.

$$\partial_x\rho=-\rho_0\alpha k\sin(k(x-vt)),\qquad \partial_t\rho=-v\partial_x\rho,\qquad \partial_t\rho+\partial_x(\rho v)=0.$$

창 안의 질량과 sinc

창 중심은 c+wt이고 반폭 h는 일정합니다. 코사인을 적분하면 중심의 위상과 반폭이 분리됩니다. 작은 k에서 사인 차이를 k로 바로 나누는 대신 sinc를 사용하면 정확한 k=0 극한을 자연스럽게 표현합니다. 창의 폭이 여러 파장을 덮으면 대비가 평균되어 작아질 수 있습니다. 이것은 밀도가 사라진 것이 아니라 관측 평균이 공간 구조를 덜 구분한다는 뜻입니다.

$$M(t)=2Ah\rho_0[1+\alpha\cos\phi(t)\operatorname{sinc}(kh)],\qquad \phi(t)=k[c+(w-v)t],\qquad \operatorname{sinc}(z)=\begin{cases}\sin z/z&z\ne0,\\1&z=0.\end{cases}$$

큰 배경을 먼저 소거한 유입 차이

양쪽 경계의 밀도는 큰 공통 배경을 포함합니다. 이를 반올림한 뒤 차감하면 작은 대비가 사라질 수 있습니다. 코사인 차이의 항등식으로 배경을 먼저 소거한 유입 차이를 계산합니다. 직접 차감한 값과 안정적인 가족 계산을 함께 비교하되 후자를 실제 센서가 잃은 정보를 복원한 값이라고 부르지 않습니다. 둘은 같은 입력 해 가족을 평가하는 다른 산술 방법입니다.

$$\dot M=2A\rho_0\alpha(v-w)\sin\phi(t)\sin(kh).$$

관측 기준계를 함께 바꾸기

유체 속도에 b를 더하고 창 속도에도 b를 더하면 상대 속도는 같습니다. 위치도 같은 방식으로 변환해야 같은 장면을 비교할 수 있습니다. 질량과 유입 차이는 유지되지만 실험실 기준의 질량 유속 자체는 바뀝니다. v만 바꾸고 창은 그대로 두는 조작은 기준계 변경이 아니라 다른 상대 운동입니다. AI 모델의 입력 변환을 시험할 때도 이 차이가 중요합니다.

$$x'=x+bt,\qquad v'=v+b,\qquad w'=w+b,\qquad v'-w'=v-w.$$

손계산: 저장2kg와 유입1kg/s

A=1m², ρ0=1kg/m³, α=1/2, k=π/2rad/m, v=1m/s를 고릅니다. 창은 c=1m, h=1m, w=0으로 두고 t=0을 봅니다. 구간은[0,2]m이며 왼쪽 밀도는1.5, 오른쪽은0.5kg/m³입니다. 코사인 대비의 구간 평균은0이므로 질량은2kg입니다. 그런데 순간 유입 차이는1kg/s입니다. 저장량과 변화율이 서로 다른 질문임을 확인하세요.

$$\rho_L=1.5,\quad \rho_R=0.5\ \mathrm{kg/m^3},\qquad M(0)=2\ \mathrm{kg},\qquad \dot M(0)=1\ \mathrm{kg/s}.$$

유한 시간 변화와 순간 기울기

처음과 끝의 질량 차이는 유입 차이를 시간에 대해 적분한 값과 비교해야 합니다. 마지막 순간의 기울기에 시간 길이를 곱하면 일반적으로 다릅니다. 손계산 예제에서0→0.5초의 질량 증가는 약0.450158kg입니다. 최초 기울기1kg/s에0.5초를 곱한0.5kg과도 다릅니다. 그래프의 한 점과 누적 면적을 각각 읽으세요.

$$\Delta M=\int_{t_0}^{t_1}\dot M\,dt,\qquad M(t)=2+\frac{2}{\pi}\sin\!\left(\frac{\pi t}{2}\right),\qquad M(0.5)-M(0)=\frac{\sqrt2}{\pi}\ \mathrm{kg}.$$

작은 변화와 정확한 영

두 큰 질량이 같은 수로 반올림되어도 가족의 작은 비영 변화는 남을 수 있습니다. 코사인 끝점 차이를 사인 곱으로 바꾸어 따로 유지합니다. 반대로 시간 길이0, 대비0, 파수0, 물질 창의 상대 속도0은 구조적으로 정확한 영입니다. 작은 비영 양의 산술 표현 손실을0으로 바꾸지 않습니다. 수치 적분의 차수 차이와 잔차는 선택한 계산의 관측이며 엄밀한 오차 상한이 아닙니다.

$$\Delta M=-4Ah\rho_0\alpha\operatorname{sinc}(kh)\sin\!\left(\frac{\phi(t_0)+\phi(t_1)}2\right)\sin\!\left(\frac{k(w-v)(t_1-t_0)}2\right).$$

극도로 짧은 시간에서는 시간을 먼저 반으로 나누는 순간 비영 중점이0으로 반올림될 수 있습니다. 구현은 사인 항등식의 작은 인자를 질량 척도와 결합할 때까지 보존합니다. 예를 들어 원점을 중심으로 고정한 창은 최초 기울기가0이어도 양의 짧은 시간 동안 이차의 비영 질량 변화가 있습니다. 반면 같은 시간의 음수부터 양수까지 대칭으로 비교하면 변화가 정확히0입니다. 표현할 수 없는 비영 결과를 보존이나 정지의 근거로 사용하지 마세요.

길이0과 좌표 해상도

h=0이면 질량과 부피는 모두0이며 질량/부피로 관측 평균을 만들 수 없습니다. 양의 h가 중심 좌표의 해상도보다 작으면 컴퓨터의 두 끝점이 같아 보일 수도 있습니다. 이 경우 양의 폭을 임의로 넓히거나 정확한 영폭으로 바꾸지 않고 상태를 기록합니다. 닫힌 가족 적분의 폭 정보와 반올림한 끝점 좌표가 서로 다른 해상도를 가질 수 있습니다.

$$\overline\rho_{\rm observed}=\frac{M}{2Ah}\quad(h>0),\qquad h=0:\quad M=0,\quad \overline\rho_{\rm observed}\text{ undefined}.$$

AI 검증과 실험 순서

먼저 손계산 예제로 단위·부호를 맞추고 고정 창과 물질 창을 비교하세요. 다음에는 정지 밀도 위로 창만 움직여 잘못된 실험실 유속이 무엇을 놓치는지 봅니다. 마지막으로 작은 파수·영폭·작은 양의 폭을 비교합니다. 연속방정식 잔차0만으로 밀도의 초기 형태나 관측하지 않은 위치를 식별할 수 없습니다. 초기·경계 데이터, 관측량과 허용 오차를 함께 지정하고 학습 결과의 점별·구간 검사를 나누어 보고하세요.

풀이가 있는 연습 문제

연습 1

물질 창에서 v=w=1m/s이면 양쪽 상대 질량 유출과 질량 변화는 얼마인가요?

해설: 상대 속도0이므로 양쪽 상대 유출과 질량 변화는 모두0입니다. 실험실 유속 ρv는 비영일 수 있습니다.

연습 2

균일 밀도1kg/m³, A=1m², v=1, w=−0.3m/s입니다. 양쪽 바깥 유출과 순유입은?

해설: 왼쪽−1.3kg/s, 오른쪽+1.3kg/s, 순유입0입니다. 양쪽에서 물질이 통과해도 저장량이 일정할 수 있습니다.

연습 3

ρ0=1, α=1/2에서 k=0이면 실제 밀도는1인가요?

해설: 아닙니다. 코사인이1이므로 실제 밀도는1.5kg/m³입니다. 기준 밀도와 특수 극한의 실제 평균을 구분합니다.

연습 4

h=0일 때 질량0을 부피0으로 나눠 평균을0이라고 해도 될까요?

해설: 안 됩니다. 관측 평균은 미정입니다. 양의 폭에서의 해석적 극한과 별도로 표시합니다.

연습 5

손계산 예제의0→0.5초 질량 변화를 구하고1kg/s×0.5초와 비교하세요.

해설: 정확한 변화는√2/π≈0.450158kg이고 상수 기울기 근사는0.5kg입니다. 실제 기울기는 시간에 따라 달라집니다.

연습 6

서로 다른 두 균일 밀도가 모두 연속방정식을 만족합니다. 잔차0으로 둘을 구분할 수 있나요?

해설: 구분할 수 없습니다. 초기 밀도나 질량 관측 등 추가 정보가 필요합니다. 보존과 식별은 다른 질문입니다.

자료와 다음 연결

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