Fisher–KPP 성장 전선과 원천 수지

성장하는 전선은 퍼지는 구름과 다르다

작은 개체군이 비어 있는 영역으로 퍼지면서 동시에 증식한다고 생각하자. 경계 수송을 모두 계산한 뒤에도 농도나 점유 영역이 증가할 수 있다. 확산만으로는 이 증가를 설명할 수 없으므로 원천항이 필요하다. Fisher–KPP 방정식은 확산에 의한 퍼짐과 수용 농도 근처에서 느려지는 성장을 결합한다. 균일한 배경 유동을 더하면 상대 전파 과정은 유지한 채 이류까지 함께 관찰할 수 있다.

여기서는 하나의 정확한 이동 프로필을 다룬다. 매끄러운 식으로 실제 장을 미분하고 유한 관측 창에서 적분할 수 있다는 장점이 있다. 임의의 초기조건을 푸는 범용 해법은 아니다. 특히 이 특수 프로필의 속도를 모든 초기조건에서 선택되는 보편 속도로 해석하면 안 된다. 해석적 예제가 무엇을 증명하고 무엇을 남겨 두는지 구분하는 것도 학습의 일부다.

가정과 단위, 열린 계의 의미

공간은 전 직선이다. 확산계수 D와 성장률 r은 엄격하게 양수인 상수이고 수용 농도 K와 단면적 A도 양수다. 배경 속도 U는 균일하게 지정한다. 농도는 음수가 아니며 정확한 프로필에서는 0과 K 사이에 있다. 성장을 가능하게 하는 자원·기질·에너지 공급은 별도 변수로 모델링하지 않는다. 따라서 원천은 모델 밖의 변수와의 교환이나 전환을 나타내며, 닫힌 계에서 물질이 무에서 생긴다는 주장이 아니다.

양	단위	의미
c, K	kg/m³	모델 물질 농도와 수용 농도
D	m²/s	확산계수
r	1/s	낮은 농도에서의 성장률
U	m/s	지정한 배경 유속
j	kg/(m² s)	부호 있는 전체 플럭스
S	kg/(m³ s)	성장 원천

$$c_t+Uc_x=Dc_{xx}+rc\left(1-\frac{c}{K}\right),\qquad j=Uc-Dc_x,\qquad c_t+j_x=S,\quad S=rc\left(1-\frac{c}{K}\right).$$

원천은 0과 K 사이에서 양수이고 양 끝에서는 0이다. 절반 농도에서 가장 크다. 하지만 국소 성장률 하나만으로 전선 전체의 이동을 결정할 수는 없다. 수송과 주변 농도 모양도 함께 작용한다. 같은 위치의 농도 변화율과 반응 원천을 비교하면서 그 차이를 확산·이류 항으로 설명해보자.

먼저 방정식의 척도를 정한다

배경 유동과 함께 움직이는 좌표를 사용하면 균일한 이류를 제거할 수 있다. 거리는 확산·성장 길이로, 시간은 성장 시간으로 측정한다. 농도를 K로 나누면 무차원이 된다. 변환된 방정식에는 물리 계수가 없어 특수 이동해를 확인하고 나중에 단위를 복원하기 쉽다.

$$u=\frac{c}{K},\qquad \xi=\sqrt{\frac rD}(x-Ut),\qquad \tau=rt,\qquad u_\tau=u_{\xi\xi}+u(1-u).$$

길이척도는 D/r 자체가 아니라 그 제곱근이다. 대응하는 속도척도는 Dr의 제곱근이다. 무차원 논문에서 본 속도 숫자를 m/s 축에 그대로 넣는 오류는 차원 점검으로 예방할 수 있다. 단위가 맞는다는 사실만으로 해가 맞는 것은 아니지만, 단위가 틀리면 물리적 해석은 이미 실패한다.

logistic 제곱 형태의 정확한 프로필

x-zero는 초기 프로필의 편리한 기준 위치다. 이 프로필은 배경 유동에 대해 오른쪽으로 전파한다. 역길이 lambda, 상대 속도 v, 실험실 속도를 다음과 같이 정의한다.

$$\lambda=\sqrt{\frac r{6D}},\qquad v=5\sqrt{\frac{Dr}{6}},\qquad X(t)=x_0+(U+v)t,\qquad z=\lambda(x-X(t)),\qquad c=\frac{K}{(1+e^z)^2}.$$

멀리 왼쪽에서는 K, 멀리 오른쪽에서는 0에 접근하며 날카로운 유한 경계는 없다. K를 바꾸면 농도 크기는 달라지지만 lambda와 v는 그대로다. 다른 계수를 고정하고 D를 늘리면 전선이 넓어지고 상대 속도가 커진다. r을 늘리면 전선은 좁아지지만 상대 속도는 역시 커진다. 두 변화가 모두 속도를 키운다고 폭까지 같은 방식으로 바뀌는 것은 아니다.

실제 공간 미분을 확인한다

서로 보완하는 logistic 인자 w와 s를 도입하자. 둘을 따로 유지하면 나중에 포화 영역의 수치 안정성도 좋아진다. 검증하려는 PDE에서 미분값을 역으로 정의하지 말고 실제 농도식을 직접 미분해야 한다.

$$w=\frac1{1+e^z},\qquad s=1-w,\qquad c=Kw^2,\qquad c_x=-2\lambda cs,\qquad c_{xx}=2\lambda^2cs(3s-1).$$

전이 영역의 1차 미분은 음수이지만 2차 미분은 부호가 바뀐다. 따라서 확산은 어떤 위치의 농도를 늘리고 다른 위치의 농도를 줄인다. D가 양수라고 D c_xx가 모든 위치에서 양수인 것은 아니다. 농도가 낮다는 이유만으로 확산이 반드시 그 위치에 유입된다고 말할 수도 없다. 주변 프로필의 곡률을 보아야 한다.

시간 미분과 방정식 대입

프로필은 모양을 바꾸지 않고 이동한다. 시간 미분은 공간 미분에 실험실 속도의 음수를 곱한 값이다. 1에서 w의 제곱을 뺀 값이 s 곱하기 1+w라는 관계를 사용하면 원천을 직접 대입하기 쉬워진다.

$$c_t=-(U+v)c_x,\qquad S=rcs(2-s),\qquad -vc_x=Dc_{xx}+S.$$

c와 s가 양수인 위치에서 마지막 식을 cs로 나누고 제안한 lambda와 v를 대입하자. 모든 위치에서 다음 항등식을 얻는다. 이는 정확한 해의 가족을 확인하는 직접 계산이다. 소프트웨어가 이 가족을 제대로 구현했는지 확인하려면 제품과 독립적으로 미분한 참조값도 필요하다.

$$2v\lambda=2D\lambda^2(3s-1)+r(2-s)=\frac{5r}{3}.$$

이류는 상대 이동 항등식에서만 사라진다. 실험실 좌표의 이동과 경계 전체 플럭스에는 여전히 남아 있다. 실험실에서 정지한 프로필에서도 확산과 성장, 이류가 활발하게 균형을 이룰 수 있다. 정지한 그림을 보고 모든 과정이 멈췄다고 해석하지 말자.

기준 위치의 농도는 4분의 1이다

z가 0이면 분모가 4이므로 X의 농도는 절반이 아니라 K/4다. 0과 1 사이의 고정 상대 농도 수준마다 서로 다른 위치가 있으며 모두 같은 실험실 속도로 움직인다.

$$c(X(t),t)=\frac K4,\qquad x_\alpha(t)=X(t)+\frac1\lambda\log\left(\alpha^{-1/2}-1\right),\qquad \frac{dx_\alpha}{dt}=U+v.$$

절반 농도에서는 로그 이동량이 음수이므로 X보다 왼쪽이다. 시각화에 4분의 1·절반·다른 농도 수준을 함께 표시하면 기준점과 절반 위치의 차이를 쉽게 볼 수 있다. K를 바꿀 때 절대 농도 임계값을 추적하는 실험은 상대 농도 수준을 추적하는 실험과 다르다. 어느 정의를 사용하는지 먼저 정해야 비교가 일관된다.

국소 수지를 손으로 계산한다

D를 1/6 m²/s, r을 1/s, K를 4 kg/m³, U를 0으로 두자. 기준 위치에서는 lambda가 1/m이고 농도가 1 kg/m³이다. 직접 대입하면 국소 수지의 모든 항을 구할 수 있다.

$$\lambda=1,\quad v=\frac56,\quad c=1,\quad c_x=-1,\quad c_{xx}=\frac12,\quad c_t=\frac56,\quad Dc_{xx}=\frac1{12},\quad S=\frac34.$$

농도/시간 단위로 확산은 1/12, 성장은 3/4를 기여하며 합은 실제 시간 미분인 5/6이다. 이 값들은 한 위치의 값이지 적분된 질량이 아니다. 한 점의 원천에 임의의 창 폭을 곱하면 일반적으로 정확한 원천 적분이 되지 않는다. 창 안에서 원천의 크기가 바뀌기 때문이다. 이 예제는 국소 방정식과 적분 수지가 같은 법칙의 서로 다른 관측이라는 점도 보여준다.

유한 창의 저장량에는 성장 적분이 들어간다

고정된 물리 관측 창에서 보존형 식을 적분하자. 오른쪽에서 왼쪽을 뺀 경계 플럭스가 순유출이다. 원천도 저장량과 같은 창에서 적분하고 같은 단면적을 곱해야 한다.

$$M_{[a,b]}=A\int_a^bc\,dx,\qquad \frac{dM_{[a,b]}}{dt}+A[j(b)-j(a)]=A\int_a^bS\,dx.$$

저장 질량 증가는 안으로 들어오는 수송, 양의 성장 원천 또는 둘 다에서 올 수 있다. 계산 모델은 농도, 실제 시간 미분, 성장 원천을 각각 적분하고 경계에서 실제 플럭스를 평가한다. 수치 구적의 수지 일치는 진단이며 임의의 함수에 대한 인증된 오차 한계는 아니다. 별도의 고정밀 참조는 제품 함수를 호출하지 않고 실제 장을 미분·적분한다. 창의 폭과 위치를 바꿔도 같은 수지가 유지되는지 비교하자.

전 직선 전체 질량은 왜 발산하는가

왼쪽 plateau는 끝없이 이어지고 양의 K에 접근한다. 따라서 전 직선 질량은 무한하다. 확산을 사용한다는 공통점이 있어도 유한 질량 Gaussian 구름과는 다른 설정이다. 화면을 유한한 구간으로 잘라 그렸다고 해석적 장의 전 직선 질량이 유한해지는 것은 아니다.

$$\lim_{x\to-\infty}c(x,t)=K>0,\qquad A\int_{-\infty}^{\infty}c(x,t)\,dx=+\infty.$$

전체 질량은 발산이라고 표시해야 한다. 그림 구간의 유한 적분에 전체 질량이라는 이름을 붙이면 안 된다. 유한 창 수지는 여전히 잘 정의되며 유용하다. 유한한 자원 저장고나 처음에 국소화된 개체군을 연구하려면 다른 초기조건·경계조건 또는 추가 방정식이 필요하다. 지금의 정확한 전선은 그런 문제의 해를 자동으로 제공하지 않는다.

포화 여량과 작은 꼬리를 안정적으로 관측한다

왼쪽에서 농도가 K로 반올림되어도 실제 여량과 성장 원천은 양수일 수 있다. 반올림된 c/K를 1에서 빼면 원천이 거짓으로 0이 된다. 오른쪽에서는 농도가 표현 범위보다 작아져도 계수를 곱한 원천은 다시 표현 가능할 수 있다. 농도의 로그를 보존하고 보완 logistic 인자를 안정적인 softplus 형태로 계산하자.

$$\operatorname{softplus}(z)=\log(1+e^z),\quad \log w=-\operatorname{softplus}(z),\quad \log s=-\operatorname{softplus}(-z),\quad K-c=Ks(1+w),\quad S=rcs(1+w).$$

softplus의 정의를 코드에 그대로 옮겨 큰 양의 z에서 지수를 먼저 계산하면 넘침이 생길 수 있다. 양수 쪽에서는 z와 작은 로그 보정의 합으로 바꿔 계산한다. 안정적인 여량과 반올림 농도의 뺄셈은 별도 관측으로 유지하자. 비영 선형 값이 표현 범위 아래라면 정확한 0으로 바꾸지 말고 제공 불가와 로그 값을 함께 보여준다. 표시 자릿수는 계산이 실제로 구분하는 정보를 전달해야 한다.

대조 실험이 증명하는 범위

같은 실제 장에서 원천을 포함한 식과 누락한 식을 비교하자. 원천 누락 잔차는 양의 성장 원천과 같아야 한다. 이는 해당 항이 필요한지 확인한다. 정확한 프로필 매개변수에서 r을 0으로 바꾸는 것은 별개다. 길이와 속도 구성이 퇴화하므로 같은 장의 대조가 아니며, D가 0인 경우도 지금 가족의 정의 범위 밖이다.

U를 상대 속도의 음수까지 바꾸면 실험실 전선은 정지한다. 더 음수로 만들면 왼쪽으로 움직이지만 상대 전파는 여전히 오른쪽이다. 여러 상대 농도 위치를 추적하며 공통 속도를 비교해보자. 이 실험은 좌표계에 대한 직관을 주지만 특수 속도가 최소 속도이거나 임의 초기조건의 선택 속도라는 주장은 하지 않는다. 속도 선택을 연구하려면 초기 꼬리와 해당 초기값 문제를 별도로 명시해야 한다.

풀이가 있는 연습

연습 1 — 절반 농도 위치

X에서 농도가 K/2인가? 절반 농도 위치는 어디인가?

해설: X의 농도는 K/4다. 절반 위치는 X에 제곱근 2에서 1을 뺀 값의 로그를 lambda로 나눈 이동량을 더한다. 이 이동량은 음수이므로 절반 위치는 왼쪽에 있다.

연습 2 — 손계산 검산

예제에서 성장만 고려하면 국소 시간 미분을 얼마로 예측하는가?

해설: 3/4로 예측해 확산 기여 1/12를 놓친다. 실제 값은 5/6이다. 반응 법칙이 정확해도 공간 재분배를 제외하면 농도 변화 전체를 설명하지 못한다.

연습 3 — 실험실 전선 정지

어떤 U에서 실험실 전선이 정지하며 그때 반응도 멈추는가?

해설: U를 −v로 둔다. 프로필이 정지하므로 시간 미분은 0이지만 반응과 확산은 계속되고 이류가 이를 상쇄한다. 그림의 정지와 과정의 정지는 다르다.

연습 4 — 전체 적분 해석

유한한 그림 구간의 적분으로 전 직선 질량을 얻을 수 있는가?

해설: 아니다. 왼쪽의 양의 plateau 때문에 전 직선 질량은 발산한다. 유한 적분은 지정한 관측 창 안의 질량이라고 이름 붙여야 한다.

연습 5 — 거짓 영 원천

농도는 K로 반올림되는데 안정 계산은 양의 성장을 보고한다. 어느 쪽을 버려야 하는가?

해설: 둘을 서로 다른 관측으로 유지한다. 반올림 농도는 작은 여량을 잃었지만 안정 보완 인자는 이를 구분한다. 원천을 강제로 0으로 만들지 말고 차이를 표시한다.

연습 6 — 특수 속도의 의미

이 정확한 이동해를 검증하면 모든 국소 초기상태가 v로 전파함을 증명한 것인가?

해설: 아니다. 이 특정 초기 프로필이 해당 상대 속도로 평행이동한다는 사실을 확인했다. 다른 초기상태의 속도 선택은 추가 가정이 필요한 별도의 문제다.

자료와 적용 범위

Ablowitz·Zeppetella의 특수 전파 속도 Fisher 방정식 논문(1979)은 특수 무차원 속도에서의 명시적 해를 제시한다. 자료 검토에서는 출판사 초록을 확인했으며 구독 본문을 읽었다고 주장하지 않는다. 이 단원의 차원 복원과 프로필 대입은 위에서 직접 유도하고 별도로 검증한 내용이다. 접근하지 못한 본문 전체의 증명을 검토했다는 뜻이 아니다.