Gaussian 농도 확산과 질량 보존

농도 봉우리가 낮아져도 물질의 양은 줄지 않을 수 있다

아주 긴 관 안에 녹아 있는 물질의 매끄러운 구름을 놓았다고 생각해봅시다. 구름이 퍼지면서 농도 봉우리는 낮아집니다. 균일한 배경 흐름이 있으면 구름 전체가 하류로 이동하기도 합니다. 두 변화는 다릅니다. 이류는 평균을 옮기고 확산은 분산을 늘립니다. 여기서 사용하는 모델에서는 어느 과정도 물질을 없애지 않습니다. 가장 높은 농도만 보면 질량 손실을 잘못 판단할 수 있습니다.

총질량과 관의 실제 단면적을 명시한 Gaussian 농도장을 만듭니다. 익숙한 Gaussian 확률 공식을 보존법칙에 연결하되 농도를 확률 자체로 읽거나 구름의 중심을 입자 궤적이라고 부르지 않습니다. 밀도, 플럭스, 적분과 모멘트에는 각각의 정의와 단위가 필요합니다.

선수 지식, 가정과 물리 단위

미분, Gaussian 적분, 평균과 분산, 유한 창 보존을 복습합니다. 관은 일정 단면적을 가진 무한 직선입니다. 확산계수는 음이 아닌 상수이고 지정한 이류 속도는 공간적으로 균일합니다. 반응, 열원이나 물질원, 벽, 변계수와 농도가 유동으로 되먹임되는 효과는 없습니다. 초기 표준편차는 양수이므로 시간 영도 특이한 점원이 아니라 매끄러운 장입니다.

물리량	단위	의미
M	kg	물질의 전체 질량
c	kg/m³	질량 농도
단면적	m²	관의 실제 단면적
U	m/s	지정한 이류 속도
D	m²/s	확산계수
j	kg/(m² s)	방향 있는 전체 질량 플럭스

$$M=A\int_{-\infty}^{\infty}c(x,t)\,dx,\qquad[c]=\mathrm{kg/m^3},\qquad[j]=\mathrm{kg/(m^2s)}.$$

희박한 스칼라 물질의 교육 모델입니다. 이 식만으로 혼합물 전체의 밀도·압력·운동량 방정식까지 주어지지는 않습니다.

이류와 Fick 플럭스에서 방정식 유도하기

Fick 구성법칙은 농도 기울기의 반대쪽으로 확산 플럭스를 정합니다. 이류의 기여는 속도와 농도의 곱입니다. 이 두 방향 있는 플럭스의 합을 연속 방정식에 넣습니다. U와 D가 상수이므로 플럭스 미분을 펼치면 익숙한 이류·확산 방정식이 됩니다. 이 관례에서 양의 U는 평균을 증가하는 x 쪽으로 이동시킵니다. 이류 항의 부호가 다른 교재의 식을 그대로 옮기지 않습니다.

$$j_{\rm adv}=Uc,\qquad j_{\rm diff}=-Dc_x,\qquad c_t+j_x=0,\qquad c_t+Uc_x=Dc_{xx}.$$

플럭스는 속도가 아닙니다. 농도가 양수인 곳에서 플럭스를 농도로 나누어 국소 비율을 만들 수는 있지만, 그것을 모든 분자의 결정론적 궤적 속도로 몰래 지정하면 안 됩니다. 확산은 퍼짐을 기술하며 이 연속체 장이 직접 모사하지 않는 미시적 요동과도 연결됩니다.

매끄러운 Gaussian 가족 구성하기

초기 구름의 전체 질량, 평균과 표준편차를 정합니다. 평균은 균일하게 병진하고 분산은 선형으로 증가합니다. 성장식에 들어가는 것은 표준편차가 아니라 분산입니다. 정규화 인자를 유지해야 퍼지면서 봉우리가 낮아져도 적분값이 그대로입니다. 단면적은 길이당 질량 분포를 실제 체적 농도로 바꿉니다.

$$\mu(t)=\mu_0+Ut,\qquad v(t)=\sigma_0^2+2Dt,\qquad c(x,t)=\frac{M}{A\sqrt{2\pi v(t)}}\exp\!\left[-\frac{(x-\mu(t))^2}{2v(t)}\right].$$

Gaussian에는 날카로운 유한 끝점이 없습니다. 양의 꼬리는 실제 장의 일부이며 그림의 결함이라며 지운 뒤 다시 정규화할 대상이 아닙니다. 화면에 보이는 유한 구간은 관측 창이지 무한 관의 물리적 경계가 아닙니다.

실제 장을 미분하기

현재 평균으로부터의 변위를 delta라고 씁니다. 지수 부분을 미분하면 공간 기울기를 얻습니다. 두 번째 미분에는 변위 제곱의 양의 항과 정규화 척도의 음의 항이 함께 있습니다. 시간 미분에서는 평균 이동과 분산 증가를 모두 고려해야 합니다. 변하는 정규화 인자를 무시하면 그림은 넓어져도 질량 보존은 틀립니다.

$$\delta=x-\mu,\quad c_x=-\frac{\delta}{v}c,\quad c_{xx}=\left(\frac{\delta^2}{v^2}-\frac1v\right)c,\quad c_t=\left[\frac{U\delta}{v}+D\left(\frac{\delta^2}{v^2}-\frac1v\right)\right]c.$$

직접 대입하면 이 가족이 PDE를 만족함을 확인합니다. 같은 식을 공유해서 얻은 잔차는 해석적 일관성 진단입니다. 실제 농도의 독립 차분과 공간 구적이 부호, 단위와 정규화 구현에 대한 추가 근거가 됩니다.

보존 질량과 정규화 모멘트 해석하기

농도에 단면적을 곱해 무한 직선 전체에서 적분하면 총질량을 얻습니다. M이 양수일 때 일차 질량 모멘트와 중심 이차 모멘트를 M으로 나눕니다. 그 평균과 분산이 선언한 Gaussian 형상 매개변수와 일치합니다. 전역적인 관측이지 모든 물질이 평균 위치에 모여 있다는 뜻은 아닙니다.

$$A\int c\,dx=M,\qquad\frac{A}{M}\int xc\,dx=\mu,\qquad\frac{A}{M}\int(x-\mu)^2c\,dx=v\quad(M>0).$$

형상이 퍼지면서 평균이 움직일 수 있고, 총질량이 일정해도 특정 유한 창의 질량은 늘거나 줄 수 있습니다. 신경망 예측을 보존 목표와 비교하기 전에 관측의 공간 범위를 구분합니다.

유한 질량을 넣은 손계산

질량2kg, 단면적1m², 초기 표준편차0.5m, 초기 평균0, 이류 속도0.5m/s, 확산계수0.125m²/s를 선택합니다.1초 후 분산은0.25에서0.5m²로 늘고 평균은0.5m가 됩니다. 봉우리 농도는2/√pi kg/m³입니다. 총질량은 여전히2kg입니다. 이 값들을 그래프의 축과 표에 대입해 단위와 정규화를 확인할 수 있습니다.

$$M=2,\ A=1,\ \sigma_0=0.5,\ \mu_0=0,\ U=0.5,\ D=0.125,\ t=1:\quad\mu=0.5\mathrm{\ m},\quad v=0.5\mathrm{\ m^2},\quad c_{\max}=\frac2{\sqrt\pi}\mathrm{\ kg/m^3}.$$

현재 평균에서 확산 플럭스는 대칭성 때문에 영입니다. 그래도 이류 속도와 농도가 양수이므로 이류 플럭스는 양수입니다. 확산 플럭스가 영이라는 사실이 전체 플럭스도 영이라는 뜻은 아닙니다.

높이, 폭과 플럭스 방향 비교하기

봉우리 높이는 표준편차에 반비례합니다. 질량과 단면적을 유지하면서 표준편차를 두 배로 늘리면 높이는 절반이 되고 분산은 네 배입니다. 대칭 구름의 왼쪽에서 확산 플럭스는 왼쪽을, 오른쪽에서는 오른쪽을 향합니다. 충분히 큰 양의 이류를 더하면 왼쪽에서도 전체 플럭스는 양수일 수 있습니다.

$$c_{\max}(t)=\frac{M}{A\sqrt{2\pi v}},\qquad j_{\rm diff}=\frac{D(x-\mu)}v c,\qquad j=\left(U+\frac{D(x-\mu)}v\right)c.$$

색이나 화살표 하나로 두 성분의 분해를 숨기지 않습니다. 각 기여는 비영인데 상쇄되어 순플럭스가 작아질 수 있습니다. 그 부호는 방향 있는 면을 통한 수송의 부호이며, 음이 아닌 농도 자체의 부호가 아닙니다.

실제 유한 창의 질량 측정하기

왼쪽과 오른쪽 위치를 고정하고 실제 구간에서 농도를 적분한 뒤 단면적을 곱합니다. Gaussian 공식은 표준정규 CDF 차이로도 쓸 수 있습니다. 수학적 항등식이지 반드시 CDF를 직접 빼서 구현하라는 뜻은 아닙니다. 먼 꼬리에서 두 CDF가 모두1로 반올림되면 구간에 비영의 물질이 있어도 부동소수점 차이는 영이 됩니다.

$$M_W=A\int_a^b c\,dx=M\left[\Phi\!\left(\frac{b-\mu}{\sqrt v}\right)-\Phi\!\left(\frac{a-\mu}{\sqrt v}\right)\right].$$

실제 폭을 유지한 적분을 로그 배율 구적으로 직접 계산합니다. 표준화한 끝점이 같은 수로 반올림되어도 원래의 물리적 끝점과 차이는 남깁니다. 두 구적 차수의 비교는 진단이며 엄밀한 오차 보증이 아닙니다. 꼬리는 빠르게 변하므로 필요한 구간 분할도 달라집니다.

창의 수지를 독립적으로 확인하기

왼쪽 바깥 법선은 음수이고 오른쪽은 양수입니다. 방향 있는 플럭스 차이에 실제 단면적을 곱하면kg/s가 됩니다. 실제 농도의 시간 미분을 적분하면 저장 변화율을 얻습니다. 두 면의 플럭스와 저장 질량의 시간 차분으로 따로 확인할 수 있습니다.

$$\frac{dM_W}{dt}+A[j(b,t)-j(a,t)]=0.$$

평균 주위의 창 밖으로 구름이 퍼지면 창 내부 질량은 감소해도 총질량은 보존됩니다. 구름이 이류로 들어오면 내부 질량이 늘 수도 있습니다. 경계 플럭스를 빼먹으면 거짓 보존 위반을 만들게 됩니다. 모든 유한 창이 언제나 전체 질량M을 담고 있다고 가정하는 것도 다른 오류입니다.

정확한 극한 대조 구분하기

D가 영이면 분산은 고정되고 Gaussian은 넓어지지 않은 채 병진합니다. U가 영이면 평균은 그대로이고 분산이 증가합니다. D가 양수이고 시간이 영이면 초기 장은 그대로지만 시간 미분이 영일 필요는 없습니다. 진화가 바로 시작되기 때문입니다. 각 대조는 방정식의 서로 다른 항을 확인합니다.

$$D=0:\ c(x,t)=c(x-Ut,0),\qquad U=0:\ \mu=\mu_0,\qquad\frac{dv}{dt}=2D.$$

어느 대조도 벽이나 반응을 도입하지 않습니다. 초기 표준편차는 계속 양수입니다. 이를 영으로 바꾸면 이 구현의 정의 범위 밖에 있는 특이 초기조건이 됩니다.

질량 영에는 질량으로 정규화한 평균이 없다

M이 정확히 영이면 농도와 모든 물리 플럭스는 영입니다. 모멘트를 정규화할 질량이 없습니다. 매개변수로 정한 형상에는 여전히 중심과 분산을 지정할 수 있지만 관측된 물질의 평균과 분산은 아닙니다. ‘측정 평균’을 영으로 보고하면 부재를 표시하는 대신 관측을 만들어내는 셈입니다.

$$M=0:\quad c=j=0,\qquad\frac{A\int xc\,dx}{M}\ \text{is undefined}.$$

자료 처리에서도 빈 관측과 비어 있지 않으면서 중심이 영인 관측은 다릅니다. 정확한 물리적 영 상태와 정규화할 수 없는 관측을 분리해서 남깁니다. 선언한 형상 매개변수에 값을 주는 것과 그 값이 자료에서 측정되었다고 주장하는 것은 다릅니다.

먼 꼬리와 로그 관측 유지하기

M이 양수이면 정규화와 지수 감쇠를 함께 결합한 뒤 농도를 일반 숫자로 바꿉니다. Gaussian 지수만 먼저 계산하면 언더플로해도 배율을 결합한 농도는 표현 가능할 수 있습니다. 반대로 농도 자체는 범위 밖이어도 큰 공간 인자가 곱해진 미분은 표현 가능할 수 있습니다. 하나의 관측을 잃었다고 다른 관측까지 모두 지워서는 안 됩니다.

$$\ell=\log\!\left(\frac c{c_*}\right)=\log\!\left(\frac{M}{A\sqrt{2\pi v}\,c_*}\right)-\frac{(x-\mu)^2}{2v},\qquad c_*=1\mathrm{\ kg/m^3}.$$

기준 농도 c-star가 로그 안을 무차원으로 만듭니다. 일반 농도는 계산 불가인데 로그 농도가 유한하면 비영의 양의 꼬리라는 뜻이며 물리적 영이 아닙니다. 구간 확률이 언더플로해도 유한한 로그 적분은 쓸 수 있습니다. 로그 질량을 표시할 때도1kg 기준임을 명시합니다. 질량 영에는 유한 로그가 없으므로 별도의 정확한 경우로 처리합니다.

반올림된 같음 아래의 작은 퍼짐 보존하기

아주 작은 양의 D나 시간에서는 분산 증가량을 초기 분산에 더해도 같은 수로 반올림됩니다. 그래도 증가량 자체는 표현 가능할 수 있습니다. 표준편차 변화를 구할 때 거의 같은 제곱근 두 개를 빼지 않습니다. 유리화 항등식은 반올림한 봉우리 비율에서 사라지는 변화도 유지합니다.

$$\Delta v=2Dt,\qquad\Delta\sigma=\frac{\Delta v}{\sqrt{v_0+\Delta v}+\sqrt{v_0}},\qquad v_0=\sigma_0^2.$$

상대 분산 증가를 eta로 두면 같은 정규화 인자에서 봉우리 비율을 얻습니다. 아래의 안정적인 식은 반올림된 제곱근이1이 되어도 음의 변화를 보존합니다.

$$\eta=\frac{\Delta v}{v_0},\qquad\frac{c_{\max}(t)-c_{\max}(0)}{c_{\max}(0)}=-\frac{\eta}{\sqrt{1+\eta}\,[\sqrt{1+\eta}+1]}.$$

안정적으로 배율을 결합해도 비영 증가량을 표현하지 못하면 그 관측을 계산 불가로 남깁니다. 정확히 확산계수가 영인 대조로 바꾸지 않습니다. 작은 병진이 덧셈에서 사라지는 경우도 반올림된 평균 이동과 안정적인 평균 변위를 구분해 표시해야 합니다.

실험 순서와 AI 연구 활용

이류를 끄고 왼쪽과 오른쪽의 확산 플럭스를 봅니다. 시간을 늘려 봉우리 농도와 분산을 비교하고, 충분히 넓은 실제 적분 또는 독립 전역 계산으로 총질량을 확인합니다. 이류를 켜서 평균 병진과 퍼짐을 구분합니다. 구름이 고정 창으로 들어오고 나갈 때 내부 질량과 경계를 비교합니다. 확산계수 영, 질량 영, 먼 꼬리 창과 작은 양의 분산 증가를 각각 대조합니다.

학습 수송 모델에서는 농도, 플럭스, 총질량, 평균과 분산의 오차를 따로 보고합니다. 봉우리만 맞추는 손실은 잘못된 폭을 보상할 수 있고, 유한 화면은 빠진 꼬리 질량을 숨길 수 있습니다. 예측 프레임을 매번 강제로 정규화하면 잘못된 PDE 플럭스를 감출 수도 있습니다. 보지 않은 이류 속도와 확산계수를 평가용으로 남기고 실제 단면적과 단위를 유지합니다. 국소 질량 변화를 해석하기 전에 경계 수지를 확인합니다. 이 Gaussian 장은 상수 계수의 통제된 기준 예이며 반응, 벽, 혼합물 동역학이나 개별 분자 경로를 맞춘다는 증거는 아닙니다.

해설이 있는 연습

연습 1: 분산과 표준편차

확산계수가0.125m²/s일 때 매초0.25m²씩 늘어나는 양은 무엇인가요?

해설: 분산입니다. 분산 증가율은 확산계수의 두 배입니다. 표준편차는 분산의 제곱근이므로 일반적으로 시간에 선형으로 증가하지 않습니다.

연습 2: 낮아진 봉우리

질량과 관 단면적을 유지하고 표준편차를 두 배로 만들면 봉우리 농도와 총질량은 어떻게 되나요?

해설: 봉우리 농도는 절반이고 총질량은 같습니다. 구름의 폭이 넓어진 것입니다. 낮은 봉우리만으로 반응이나 질량 손실을 결론 내릴 수는 없습니다.

연습 3: 대칭성과 전체 플럭스

현재 평균에서 확산 플럭스가 영이면 전체 플럭스도 영인가요?

해설: 아닙니다. 전체에는 U와 농도의 곱인 이류 플럭스가 있습니다. 이류와 농도가 비영이면 확산 기여가 대칭성으로 사라져도 전체 플럭스는 비영입니다.

연습 4: 빈 구름

총질량이 영일 때 선언한 Gaussian 중심을 측정된 질량 정규화 평균으로 보고해도 되나요?

해설: 안 됩니다. 물리적 장은 영이고 정규화하려면 영으로 나누게 됩니다. 형상 매개변수와 관측할 물질이 없어 정의할 수 없는 평균을 구분합니다.

연습 5: 꼬리 CDF 뺄셈

표준정규 CDF 두 개가 모두1로 반올림되었습니다. 차이가 영이라는 사실이 두 실제 끝점 사이에 물질이 없음을 증명하나요?

해설: 아닙니다. 구간을 직접 적분하거나 안정적인 로그 적분을 남깁니다. 원래 물리적 폭을 보존하고 꼬리를 지우거나 보이는 부분을 재정규화하지 않습니다.

연습 6: 국소 질량 변화

전역 참조의 총질량은 같은데 고정 창 내부 질량이 감소합니다. 보존 위반이라고 판단하기 전에 무엇을 확인해야 하나요?

해설: 이류와 확산을 모두 포함하고 실제 단면적을 곱한 방향 있는 경계 플럭스를 확인합니다. 유출이 저장 변화의 원인일 수 있습니다. 구적 오차와 작은 관측의 계산 불가도 물리적 부재와 분리합니다.

자료와 적용 범위

Gilbert Strang의 MIT OCW 열·이류확산 노트는 Gaussian 퍼짐과 기본해를 다룹니다. MIT의 확산 시각화 설명은 작은 확산계수의 조절을 설명합니다. 여기서는 질량 정규화, 실제 관 단면적, 플럭스 부호와 산술 대조를 직접 유도했습니다. 해당 시각화의 초기 폭 매개변수는 이 단원의 표준편차와 관례가 다르며, 오래된 Java applet을 검증 수단으로 사용하지 않습니다.