가역 반응과 종별 수송 수지

구름이 움직이는 동안 물질의 표지가 바뀐다

긴 관 속의 추적 물질이 두 관측 상태 A와 B 사이를 오간다고 생각하자. 파란색 A가 빨간색 B로 바뀌면 구름의 색은 달라지지만 전체 물질은 그대로일 수 있다. A만 감지하는 장치는 신호가 줄었다고 보고한다. 이 기록만으로는 반응으로 A가 전환된 것인지, 구름이 관측 영역 밖으로 이동한 것인지 알 수 없다. 각 상태의 수지와 둘을 합친 전체 수지를 함께 세워야 한다.

여기서는 임의의 반응을 푸는 범용 수치 해법 대신 정확하게 유도할 수 있는 특수한 해의 가족을 만든다. 전체 구름은 앞 단원의 Gaussian 이류·확산 해를 따르고, 공간적으로 같은 조성 비율이 가역 일차 반응에 따라 변한다. 두 과정을 분리하면 실제 농도와 미분, 경계 플럭스를 각각 계산할 수 있다. 목표는 공식을 외우는 것이 아니라 측정된 변화 중 얼마가 이동이며 얼마가 전환인지 설명하는 것이다.

가정과 물리 단위

관은 단면적 A가 일정한 무한 직선이다. 두 표지는 같은 질량 기준을 사용하며 일정한 속도 U와 확산계수 D를 공유한다. 정반응과 역반응의 속도상수는 음수가 아니고 시간·위치·농도와 무관하다. 처음의 A 비율은 모든 위치에서 같다. 유동은 외부에서 지정하며 압력, 온도, 운동량, 농도가 유동에 주는 되먹임은 계산하지 않는다. 같은 질량의 표지 전환이라는 설정을 분자량이 다른 임의의 화학 반응으로 확대하면 안 된다.

양	단위	의미
c_A, c_B	kg/m³	각 표지의 질량 농도
j_A, j_B	kg/(m² s)	이류와 확산을 합친 부호 있는 플럭스
k_f, k_b	1/s	정반응·역반응 속도상수
S_A, S_B	kg/(m³ s)	순전환 원천
M	kg	전 직선의 유한한 전체 질량

$$j_i=Uc_i-D\partial_xc_i,\qquad \partial_tc_i+\partial_xj_i=S_i,\qquad i\in\{A,B\}.$$

양의 플럭스는 x가 증가하는 방향을 뜻한다. 양의 원천은 해당 종을 그 위치에서 늘린다. A의 원천이 양수라고 전체 물질까지 증가하는 것은 아니다. 농도가 줄어드는 현상을 해석할 때도 이동 방향과 생성·소멸의 부호를 별도로 읽어야 한다.

결합된 보존식을 유도한다

정반응 한 번은 A를 없애고 같은 질량의 B를 만든다. 역반응은 반대다. 따라서 A의 순원천은 역반응 유입에서 정반응 유출을 뺀 값이다. 두 종의 식을 더하면 원천이 정확히 상쇄된다. 이 상쇄는 가정한 화학양론의 성질이며, 어떤 수치 알고리즘이 우연히 얻어낸 결과가 아니다.

$$S_A=-k_fc_A+k_bc_B,\qquad S_B=-S_A,\qquad c=c_A+c_B,\qquad c_t+Uc_x=Dc_{xx}.$$

각 종이 빠르게 바뀌더라도 합은 원천 없는 수송 방정식을 따른다. 두 확산계수가 다르면 전체 질량 수지 자체는 여전히 세울 수 있지만, 전체 플럭스를 하나의 공통 D와 전체 농도 기울기로 표현하는 것은 일반적으로 불가능하다. 다음의 분리형 구성이 왜 가정에 의존하는지 보여주는 예다.

조성과 Gaussian 운반장을 분리한다

처음부터 양의 폭을 가진 매끄러운 구름을 사용하자. 전체 질량은 M, 초기 중심은 mu-zero, 초기 표준편차는 sigma-zero다. 같은 운반장에 시간에만 의존하는 종별 비율을 곱한다. 초기 조성이 공간적으로 일정하고 두 표지의 수송 연산자가 같기 때문에 가능한 구성이다.

$$v(t)=\sigma_0^2+2Dt,\quad \mu(t)=\mu_0+Ut,\quad c=\frac{M}{A\sqrt{2\pi v}}e^{-(x-\mu)^2/(2v)},\quad c_A=f_A(t)c,\quad c_B=f_B(t)c.$$

단면적은 길이당 질량을 부피당 농도로 바꾼다. 확률의 정규화 상수가 아니라 실제 관의 단면적이다. A 방정식에 곱 형태를 대입하면 운반장 방정식에 f_A를 곱한 항과 조성 변화에 따른 f_A' c가 남는다. 운반장 항이 상쇄되므로 반응은 하나의 상미분방정식으로 줄어든다.

$$f_A'=-k_ff_A+k_b(1-f_A),\qquad f_B=1-f_A,\qquad f_A+f_B=1.$$

마지막 두 식은 해석적 모델의 항등식이다. 컴퓨터에서 B를 계산할 때는 반올림된 A를 1에서 빼는 대신 두 양의 비율을 따로 계산해야 작은 생성량을 보존할 수 있다. 수학적 항등식과 부동소수점 계산 절차는 구분하자.

평형으로의 완화와 영 반응률 예외

두 반응률의 합이 양수이면 그 역수가 조성 완화의 시간척도다. 평형에서는 정반응과 역반응의 양이 같다. 조성은 지수적으로 평형에 접근하며, 지금의 가정에서는 구름의 위치·폭·이동 속도가 이 조성 곡선을 바꾸지 않는다.

$$\lambda=k_f+k_b>0,\qquad f_A^*=\frac{k_b}{\lambda},\qquad f_A(t)=f_A^*+(f_{A0}-f_A^*)e^{-\lambda t},\qquad \tau=\lambda^{-1}.$$

두 반응률이 모두 0이면 모든 초기 조성이 그대로 정상 상태다. 유일한 선호 평형이나 위 식으로 정의되는 유한한 완화 시간은 없다. 이때 평형을 0으로 채워 넣으면 원래 방정식에 없던 선호를 만드는 셈이다. 전체 질량 M이 0일 때도 선언된 조성 변수의 시간 변화는 수학적으로 계산할 수 있지만, 관측 질량을 전체 질량으로 나눈 측정 조성은 0/0이므로 정의되지 않는다.

총전환량과 순원천은 다르다

평형인 혼합물에서도 A와 B의 전환은 계속 일어날 수 있다. 양방향 전환량이 같으면 순원천은 0이지만 사건 자체가 사라지는 것은 아니다. 농도 변화가 아니라 사건 수나 에너지 사용을 측정하는 연구에서 특히 중요한 구분이다.

$$R_f=k_fc_A,\qquad R_b=k_bc_B,\qquad S_A=R_b-R_f,\qquad S_A+S_B=0.$$

같은 양의 반응률과 절반씩의 조성에서는 원천이 0이어도 농도가 양수인 곳의 양방향 전환량은 양수다. 반대로 물질이 전혀 없는 곳에서도 원천은 0이다. 따라서 순원천이 0이라는 숫자 하나만 보고 반응이 없다고 결론 내릴 수 없다. 시각화에서는 전체 농도, 두 방향 전환량, 부호 있는 순원천을 함께 비교하는 편이 이해에 도움이 된다.

실제 미분에서 원천항을 찾는다

곱으로 정의한 실제 농도를 미분하자. 시간 미분에는 운반장 변화와 조성 변화가 모두 들어간다. 공간 미분에는 조성 기울기가 없는데, 비율이 위치에 의존하지 않기 때문이다. 원천은 바로 추가적인 시간 변화를 담당한다.

$$\partial_tc_A=f_Ac_t+f_A'c,\qquad \partial_xj_A=f_A\partial_xj,\qquad \partial_tc_A+\partial_xj_A=f_A'c=S_A.$$

이 장을 그대로 두고 방정식에서 반응만 누락하면 잔차는 원천과 같아진다. 반응이 필요한 장에서 원천 없는 식이 실패하는지 확인하는 대조 실험이다. 장 자체를 새로운 무반응 해로 바꿔 놓고 잔차가 작다고 말하면 원래 문제의 누락을 찾을 수 없다. AI가 생성한 방정식이나 물리 제약의 필요성을 검토할 때도 같은 입력과 같은 실제 장을 유지해야 비교의 의미가 선명해진다.

조성 변화를 손으로 계산한다

정반응률 1/s, 역반응률 0.5/s, 초기 A 비율 0.8, 전체 질량 2 kg을 택하자. 1초 뒤에는 A의 비율이 평형인 3분의 1 쪽으로 상당히 이동한다. 화면에는 반올림해 표시하더라도 다음 계산에는 정확한 식을 사용한다.

$$f_A(1)=\frac13+\left(0.8-\frac13\right)e^{-1.5}\approx0.43746,\qquad M_A(1)\approx0.87492\ \mathrm{kg},\qquad M_B(1)\approx1.12508\ \mathrm{kg}.$$

A 질량은 초기 1.6 kg에서 줄었지만 전체 질량은 2 kg이다. 이때 A 비율의 시간 변화율은 약 −0.15619/s이고, 전 직선의 A 원천은 약 −0.31238 kg/s이다. B 원천은 크기가 같고 부호가 반대다. U를 바꾸면 두 구름이 함께 이동하고 D를 바꾸면 함께 퍼지지만, 분리형 가족의 조성 곡선은 바뀌지 않는다. 이 결과를 이용하면 수송 계산과 반응 계산의 역할을 따로 점검할 수 있다.

유한 관측 창에는 경계와 반응이 모두 필요하다

실제 장치는 무한한 관 전체가 아니라 유한한 영역을 본다. 국소식을 고정된 창에서 적분하고 단면적을 곱하자. 오른쪽 경계 플럭스에서 왼쪽 경계 플럭스를 뺀 것이 순유출이다. 저장량 감소는 양의 순유출, 음의 반응 원천 또는 둘의 결합으로 발생한다.

$$M_i^{[a,b]}=A\int_a^bc_i\,dx,\qquad \frac{dM_i^{[a,b]}}{dt}+A[j_i(b)-j_i(a)]=A\int_a^bS_i\,dx.$$

두 종의 창 수지를 더하면 반응 원천 적분은 사라지지만 경계 수송은 남는다. 따라서 창 안의 전체 질량은 일정할 필요가 없다. 반면 전 직선에서는 Gaussian의 무한대 경계 플럭스가 사라지므로 전체 질량이 일정하다. 그래프의 질량이 전 직선 질량인지 관측 창 질량인지 먼저 확인하자. 창을 구름 중심에 놓았는지 가장자리에 놓았는지에 따라서도 저장 변화의 해석이 달라진다.

작은 변화와 작은 꼬리를 안정적으로 계산한다

처음에 순수한 A만 있을 때, 표시된 A 비율이 1로 반올림되더라도 새로 생긴 B 비율은 양수일 수 있다. 두 비율을 양수 항의 합으로 따로 계산하고 작은 지수 차이에는 expm1 형태를 사용한다. 초기 A 비율 변화율을 d-zero라고 쓰자.

$$d_0=k_b(1-f_{A0})-k_ff_{A0},\qquad f_A'(t)=d_0e^{-\lambda t},\qquad \Delta f_A=\frac{d_0}{\lambda}(1-e^{-\lambda t}).$$

마지막 식은 lambda가 양수일 때만 쓴다. 정반응률이 10의 −20승/s이고 역반응이 없으면, 1초 뒤 B 비율은 대략 10의 −20승이다. 반올림된 A를 1에서 빼면 0이 나오지만 B가 수학적으로 없다는 뜻은 아니다. 반대로 아주 빠른 완화 뒤의 비영 변화율은 표현 가능한 최소 양수보다 작아질 수 있다. 이때 선형 관측은 제공 불가로 표시하고 부호와 로그 크기를 보존해야 한다. 구름에서 멀리 떨어진 Gaussian 꼬리 역시 그림에 보이지 않는다고 잘라내거나 다시 정규화하지 않는다.

아주 작은 반응률에서도 지수 계산 전의 정보가 중요하다. 처음의 변화율이 영과 최소 양수의 정확한 중간 크기라면, 양의 시간만큼 완화한 뒤에는 그 크기가 조금 더 작아진다. 감쇠 인자를 먼저1로 반올림하면 관측이 반올림 경계의 잘못된 쪽에 놓일 수 있다. 작은 무차원 지수에서는 다항식 보정과 입력의 곱·비를 유지한 뒤 마지막에 반올림한다. 비율 변화율을 표시할 수 없어도 질량을 곱한 전체 원천량은 표현 가능할 수 있으므로 각각 계산한다. 이 처리가 Gaussian 구적이나 모든 초월함수 계산까지 정확 산술 솔버로 바꾸는 것은 아니다.

정확한 정상 상태는 별도의 조건이다. 반응률이2/s와1/s이면 수학적 평형은3분의1이다. 하지만 이진 부동소수점에 입력한3분의1은 그 근처의 유리수이므로, 화면의 두 비율이 같아 보여도 작은 초기 변화율이 정확히 영일 필요는 없다. 같은 반응률과 초기 비율0.5는 정확한 정상 대조가 된다. 반올림된 같음을 평형 측정이라고 부르기 전에 지정한 입력, 변화율의 부호와 안정적 변화를 함께 확인하자.

시각화 실험과 연구 점검

먼저 수송을 고정하고 정반응률을 바꾸며 두 색의 농도, 합, 조성 곡선을 비교하자. 다음에는 반응률을 고정하고 속도를 바꾸자. 구름은 이동하지만 전 직선 종별 비율은 같다. 마지막으로 구름 가장자리 근처에 좁은 창을 놓고 저장 변화, 순유출, 원천 적분을 각각 관찰하자. 국소 농도 기록만으로 반응률을 식별하기 어려운 이유가 이 실험에서 드러난다. 관측 창을 바꾼 추가 데이터가 식별성에 어떤 도움을 주는지도 질문할 수 있다.

물리 제약 AI에서는 작은 학습 잔차만으로 정확성을 주장할 수 없다. 독립적으로 구성한 장의 농도와 미분, 실제 창의 원천 적분, 원천 누락 대조를 함께 확인해야 한다. 함께 제공하는 참조 데이터는 제품 계산과 별개인 고정밀 미분·적분을 사용한다. 같은 공식을 모든 항의 정의와 검증에 반복 사용하면 부호 오류도 함께 복제될 수 있다. 공간에 따라 달라지는 초기 조성, 서로 다른 확산계수, 온도 의존 반응률로 확장할 때는 새로운 유도와 검증이 필요하다.

풀이가 있는 연습

연습 1 — 평형 비율

정반응률이 2/s이고 역반응률이 1/s이면 정상 A 비율은 얼마인가?

해설: 3분의 1이다. 정반응 전환량은 2 곱하기 3분의 1 곱하기 전체 농도이고, 역반응 전환량은 1 곱하기 3분의 2 곱하기 전체 농도다. 둘은 같으며 농도가 양수이면 각각 0이 아니다.

연습 2 — 사라지는 A 신호

고정 창의 A 농도가 감소하면 전체 물질이 반드시 사라졌는가?

해설: 아니다. A가 B로 전환될 수 있고 두 종 모두 경계를 통과할 수 있다. 두 종과 경계 플럭스를 측정하거나 추정한 뒤 합의 창 수지를 비교해야 한다. 전 직선 질량 보존이 창 질량의 일정함을 뜻하지 않는다.

연습 3 — 두 반응률이 모두 0

두 반응률이 0이고 초기 A 비율이 0.7이면 조성과 평형을 어떻게 보고해야 하는가?

해설: 조성은 0.7로 유지된다. 모든 비율이 정상이므로 하나의 평형값이나 유한한 완화 시간을 선택하지 말고 평형이 유일하지 않다고 보고한다.

연습 4 — 누락된 원천

실제 반응장을 유지하고 A 식에서 반응만 누락하면 잔차는 무엇인가?

해설: f_A' c, 즉 실제 A 원천이다. 부호는 순생성인지 순소멸인지 알려준다. 같은 장을 사용하는 대조라야 수송장을 바꾸지 않고 반응항의 필요성을 확인할 수 있다.

연습 5 — 반올림된 비율

A 비율이 1로 표시되는데 안정 계산에서 작은 양의 B가 나오면 B를 강제로 0으로 만들어야 하는가?

해설: 아니다. A 쪽에서 작은 차이가 반올림으로 사라졌을 뿐이다. 독립적으로 계산한 B를 보존하고 반올림된 합과 해석적 정규화를 구분하자. 강제 절단이나 재정규화는 관찰하려던 정보를 버린다.

연습 6 — 분리형 가정 깨기

초기 A 비율이 위치에 따라 다르면 같은 곱 구성이 그대로 유효한가?

해설: 일반적으로 아니다. 조성의 공간 미분이 추가 기울기 항을 만든다. 하나의 시간 의존 비율로는 전체 조성장을 나타낼 수 없으므로 공간을 포함한 반응수송 문제를 다시 풀어야 한다.

자료와 적용 범위

Eymard·Hilhorst·Olech의 가역 반응확산 유한체적 논문은 결합된 가역 반응확산 방정식과 유계 영역 조건의 일차 자료다. 여기의 무한 직선·공통 수송 Gaussian 가족은 별도로 직접 유도한 교육용 특수화다. 논문의 수치 수렴 결과가 이 구현을 인증하거나 논문의 경계조건이 지금의 모델에 적용되는 것은 아니다.