입자 커널과 연속장

점 입자에서 연속장으로

멀리서 보면 강물은 매끄럽지만 가까이 보면 분자와 입자들이 움직입니다. 이 둘을 연결하려면 어디까지 평균하는지 정해야 합니다. 사진을 흐리게 만드는 것처럼 각 입자의 질량을 작은 이웃에 펼쳐 보세요. 곡선은 매끄러워져도 실제 입자 수가 늘거나 충돌이 생긴 것은 아닙니다. 이 단원은 평균으로 만든 장의 정확한 보존 성질과 그 장을 유체 모형이라고 해석하기 위해 추가로 필요한 가정을 나눠 살펴봅니다.

선수 지식과 선밀도의 단위

정적분, 미분의 연쇄 법칙, 질량 가중 평균과 이전 이동 제어체적 단원을 복습하세요. 여기서는 실제 삼차원 기체를 재현하지 않고 일차원 직선 위의 질량을 다룹니다. 따라서 밀도는 kg/m인 선밀도입니다. 부피 밀도 kg/m³와 같은 숫자로 비교하면 안 됩니다. 입자 질량은 kg, 속도는 m/s이고 커널의 적분은 무차원1입니다. 이후의 삼차원 계수 요동 실험과도 차원을 구분합니다.

$$[m_i]=\mathrm{kg},\quad [W_h]=\mathrm{m^{-1}},\quad [\rho_h]=\mathrm{kg/m},\quad [j_h]=\mathrm{kg/s}.$$

관측 커널을 먼저 정의하기

각 입자에서 거리 h 안에만 질량을 펼치는 대칭 커널을 선택합니다. h는 지지 구간의 반폭이며 표준편차나 격자 간격과 같은 이름이 아닙니다. 아래 사차 다항식의 계수는 적분을1로 맞춥니다. 지지 경계에서 값과 일차 미분은0으로 이어지므로 C1이지만 이차 미분까지 연속인 C2 함수는 아닙니다. 그림이 매끄러워 보인다는 이유로 임의 차수의 미분 가능성을 가정하지 마세요.

$$W_h(s)=\frac{15}{16h}\left(1-\frac{s^2}{h^2}\right)^2\mathbf1_{|s|<h},\qquad \int_{\mathbb R}W_h(s)\,ds=1.$$

같은 입자와 실제 궤적

입자 i는 처음 위치 Xi에서 일정한 속도 vi로 이동합니다. 이 예제는 외력·입자 간 힘·충돌이 없는 탄도 운동입니다. 서로 지나치는 입자를 반사시키거나 충돌 속도로 바꾸지 않습니다. 가스의 충돌과 국소 평형을 자동으로 얻는 모형이 아니라 평균 연산을 정확히 조사하기 위한 유한 예제입니다. 초기 위치와 현재 위치를 함께 보면 장의 이동이 어디서 왔는지 추적할 수 있습니다.

$$x_i(t)=X_i+v_it,\qquad \dot x_i=v_i,\qquad \dot v_i=0.$$

밀도와 질량 유속은 따로 평균한다

입자마다 질량에 커널 값을 곱해 더하면 선밀도가 됩니다. 속도까지 곱해 더하면 질량 유속입니다. 관측 속도는 두 합의 비로 정의합니다. 질량이 다른 두 입자를 같은 비율로 세는 단순 속도 평균과는 다릅니다. 질량이 큰 입자가 같은 커널 위치에 있으면 그 입자의 속도가 더 크게 반영됩니다. 표시할 때도 밀도·유속·정규화한 속도를 서로 다른 단위와 축으로 읽어야 합니다.

$$\rho_h=\sum_i m_iW_h(x-x_i),\qquad j_h=\sum_i m_iv_iW_h(x-x_i),\qquad u_h=j_h/\rho_h\quad(\rho_h>0).$$

보존식은 연쇄 법칙에서 나온다

폭 h를 고정하면 커널의 시간 변화는 입자 이동에서만 옵니다. 연쇄 법칙을 적용하면 시간 미분에−vi가 붙고 공간 미분과 상쇄합니다. 따라서 이 유한 입자 장은 국소 질량 보존식을 만족합니다. 이것은 입자 수가 충분히 많거나 실제 유체가 평형이라는 증거가 아닙니다. 하나의 고립된 입자를 넓게 펼친 경우에도 같은 항등식이 성립합니다. 보존식과 물리적 적용 가능성을 구분하는 핵심 예제입니다.

$$\partial_t\rho_h=-\sum_i m_iv_iW_h^{\prime}(x-x_i),\qquad \partial_xj_h=\sum_i m_iv_iW_h^{\prime}(x-x_i),\qquad \partial_t\rho_h+\partial_xj_h=0.$$

평균 속도로 놓치는 이차 모멘트

운동량 유속에는 속도의 제곱이 들어갑니다. 따라서 평균 속도를 제곱하는 연산과 속도를 제곱한 뒤 평균하는 연산은 다릅니다. 평균으로 설명하지 못한 부분은 현재 관측 속도에서 벗어난 속도들의 제곱을 가중해 더한 값입니다. 이 중심화된 계산은 큰 두 수를 차감하는 방법보다 작은 분산을 읽기 좋고 음수가 되지 않습니다. 여기서 얻은 분산 수송을 충돌 기체의 열역학적 압력이나 온도로 바로 부르지는 않습니다.

$$Q_h=\sum_i m_iv_i^2W_h=\rho_hu_h^2+S_h,\qquad S_h=\sum_i m_i(v_i-u_h)^2W_h\ge0.$$

탄도 운동량 보존과 폐쇄 문제

입자 속도가 일정하므로 운동량 밀도 j의 시간 변화는 Q의 공간 변화와 상쇄합니다. 밀도와 평균 속도만 예측하고 S를0으로 두면 원래 입자들이 운반하던 운동량을 놓칠 수 있습니다. 이것이 더 많은 모멘트나 구성 관계가 필요한 이유를 보여주는 작은 예입니다. 다만 이 식을 곧바로 Navier–Stokes의 점성 응력이라고 읽지 마세요. 점성·충돌·열전도는 다른 물리 모형에서 유도하거나 측정해야 합니다.

$$\partial_tj_h+\partial_xQ_h=0,\qquad \partial_tj_h+\partial_x(\rho_hu_h^2)=-\partial_xS_h.$$

손계산: 평균 속도0의 두 흐름

질량1kg인 입자 두 개를 원점에 두고 속도를 각각+1과−1m/s로 정합니다. t=0, x=0, h=1m에서 커널은15/16m⁻¹입니다. 밀도는15/8kg/m, 질량 유속은0이고 평균 속도도0입니다. 그러나 운동량 유속과 분산 수송은 모두15/8N입니다. 서로 반대로 흐르는 두 집단을 정지한 한 집단으로 바꾸면 질량과 평균 속도는 같아도 이차 모멘트가 달라집니다.

$$\rho_h=\frac{15}{8}\ \mathrm{kg/m},\quad j_h=u_h=0,\quad Q_h=S_h=\frac{15}{8}\ \mathrm N.$$

유한 화면은 전체 공간이 아니다

커널을 전체 직선에서 적분하면 각 입자의 질량이 그대로 보존됩니다. 유한 관측 창에서는 커널의 일부만 포함될 수 있습니다. 화면에 보이는 질량을 전체 질량으로 다시 정규화하면 실제 유출이나 시야 밖의 질량을 숨깁니다. 두 입자 손계산에서[−0.5,0.5]m 안의 질량은203/128kg입니다. 이는 전체2kg보다 작지만 질량이 사라진 것이 아닙니다. 창을 나눠 적분한 값도 원래 창의 적분과 합이 맞아야 합니다.

$$\int_{\mathbb R}\rho_h\,dx=\sum_i m_i,\qquad M_{[a,b]}=\sum_i m_i\int_a^bW_h(x-x_i)\,dx.$$

빈 지지 구간과 영길이 창

어떤 커널도 닿지 않는 위치에서는 밀도와 유속이 모두0입니다. 그 비인 속도는0/0이므로 미정입니다. 이를 정지 유체의 속도0으로 바꾸면 빈 영역과 물질이 정지한 영역을 혼동하게 됩니다. 길이0의 창도 질량은0이지만 평균 선밀도는 미정입니다. 제한된 교육용 계산은 h를0.02~2m로, 양의 창 길이는10⁻⁸m 이상으로 한정합니다. 이보다 작은 입력은 임의로 확대하거나 다른 실험으로 바꾸지 않고 거부합니다.

밀도가 양수여도 모든 속도 모멘트를 일반 표시 범위에 담을 수 있는 것은 아닙니다. 속도가10⁻²⁰⁰m/s 정도이면 양수인 제곱도 중간 부동소수점 곱셈에서0으로 반올림됩니다. 여기서는 입력 binary64 값의 커널 모멘트와 미분을 정확한 유리수 식으로 누적한 뒤 마지막에 한 번 반올림합니다. 분산은 입자 쌍 사이 속도 차이의 제곱을 비음수 가중합으로 더하는 동치식을 사용해 큰 모멘트의 차감도 피합니다. 산술 표는 정확한0, 선형 표시 범위 밖의 비영 값과 그 부호·로그, 빈 지지 구간의 미정 평균을 구분합니다. 다만 입자 좌표와 지지 구간 판정은 일반 부동소수점 기하로 계산하므로 그 단계에서 잃은 위치 정보를 복구하지는 않습니다.

평균 폭이 변할 때 생기는 항

입자를 그대로 두고 평균 폭만 넓혀도 한 위치의 밀도는 낮아질 수 있습니다. 이것을 실제 질량 생성이나 소멸로 해석하면 안 됩니다. h가 시간에 따라 바뀌면 시간 미분에 폭의 변화율에 해당하는 항이 추가됩니다. 같은 입자·같은 시각에서 폭을 바꾸는 실험으로 이 관측 효과를 분리하세요. 공간적으로 변하는 폭이나 경계에서 다시 정규화한 커널에는 더 많은 항이 필요하며 이번 계산에는 포함하지 않습니다.

$$\partial_t\rho_h+\partial_xj_h=\dot h\sum_i m_i\partial_hW_h(x-x_i),\qquad \partial_tj_h+\partial_xQ_h=\dot h\sum_i m_iv_i\partial_hW_h(x-x_i).$$

단일 입자의 관측 변화 계산

원점의 정지한1kg 입자와 h=1m를 보세요. 원점 밀도는15/16kg/m입니다. 같은 입자를 h=2m로 관측하면15/32kg/m로 절반이 되지만 전체 적분은 여전히1kg입니다. h=1m에서 폭이 초당0.25m로 증가하는 순간의 밀도 변화는−15/64kg/(m·s)입니다. 물리적 유속은0이므로 이 변화는 오직 관측 폭에서 옵니다. 두 화면의 커널 폭이 같았는지 확인하지 않고 잔차를 비교하면 잘못된 결론을 낼 수 있습니다.

$$\rho_h(0)=\frac{15m}{16h},\qquad \partial_t\rho_h(0)=-\frac{15m\dot h}{16h^2}.$$

기준계를 바꾸어도 남는 분산

모든 입자 속도에 같은 b를 더하고 관측 좌표도 x+bt로 옮겨 보세요. 입자와 관측 위치의 상대 거리가 유지되므로 밀도는 같고 평균 속도에는 b가 더해집니다. 중심 속도 편차는 같으므로 S도 유지됩니다. 유속 자체는 기준계에 따라 변하므로 같은 수가 나와야 한다고 강제하면 안 됩니다. 변환할 양과 불변인 양을 나눠 확인하는 검사는 AI가 좌표 표현을 물리 현상과 혼동하는지 점검하는 데도 유용합니다.

$$x^{\prime}=x+bt,\quad v_i^{\prime}=v_i+b,\quad j_h^{\prime}=j_h+b\rho_h,\quad u_h^{\prime}=u_h+b,\quad S_h^{\prime}=S_h.$$

실험 순서와 AI의 식별 한계

먼저 반대 속도 두 입자에서 평균 속도0과 분산 수송의 비영 값을 비교하세요. 같은 속도로 바꾸면 분산은0이 됩니다. 다음에는 h를 줄여 빈 구간과 좁은 피크를 보고, 넓혀 이웃 입자의 정보가 합쳐지는 과정을 봅니다. 마지막으로 폭 변화율을 켜서 관측 항과 물리적 유출을 구분합니다. 같은 몇 개의 평균값을 내는 입자 배치는 여러 개일 수 있으므로 매끄러운 출력과 보존 잔차만으로 미시 상태나 국소 평형을 유일하게 찾을 수 없습니다.

풀이가 있는 연습 문제

연습 1

원점의 같은 질량 두 입자가±1m/s로 움직입니다. 평균 속도0이면 운동량 유속도0인가요?

해설: 아닙니다. h=1m인 손계산에서 Q=S=15/8N입니다. 제곱 속도는 부호 상쇄를 하지 않습니다.

연습 2

질량1kg와3kg이 같은 위치에서−1,+1m/s로 움직이면 관측 속도는?

해설: 커널이 같으므로(−1+3)/4=0.5m/s입니다. 단순 산술 평균0과 다릅니다.

연습 3

모든 커널 밖의 밀도와 유속이0일 때 속도를0으로 표시해도 될까요?

해설: 안 됩니다. 정규화는0/0이므로 미정입니다. 물질이 정지한 관측과 다릅니다.

연습 4

정지한1kg 입자를 h=1m에서2m로 관측하면 원점 밀도와 전체 질량은?

해설: 원점 밀도는15/16에서15/32kg/m로 줄지만 전체 질량은1kg입니다.

연습 5

두 입자 예제에서[−0.5,0.5]m의 질량이2kg보다 작으면 질량 보존 위반인가요?

해설: 아닙니다. 창 안에는203/128kg가 있고 나머지는 창 밖의 커널에 있습니다.

연습 6

커널이 만든 연속방정식 잔차0은 충돌·국소 평형을 증명하나요?

해설: 아닙니다. 탄도 입자 한 개도 같은 항등식을 만족합니다. 적용 척도와 분포에 대한 별도 증거가 필요합니다.

자료와 다음 연결

연속체 평균과 척도 분리의 배경은 Natalia Berloff의 Averaging in Fluid Dynamics, 평균 자유 행로와 충돌 시간의 구분은 David Tong의 Kinetic Theory §1.2를 참고하세요. 이 교재의 커널·유한 입자·선택한 세 부등식은 설명을 위한 별도 예제입니다. 계산 모형을 위 자료가 인증한다는 의미는 아닙니다.