Peclet 수와 이류·확산 경계층

단위를 바꾸어도 수송의 성격은 바뀌지 않는다

곧은 통로에 염료가 흐른다고 생각해보자. 흐름은 염료를 운반하고, 확산은 농도 차이를 완만하게 만든다. 같은 속도를 0.5m/s와 50cm/s로 표시할 수 있다. 숫자가 커졌다는 이유로 이류가 강해진 것은 아니다. 이류와 확산을 비교하려면 단위 선택에 따라 달라지지 않는 비가 필요하다. Peclet 수가 그 비교를 제공한다.

여기서 다른 실험도 생각할 수 있다. 통로를 실제로 더 길게 만들고 확산계수와 속도를 조절하여 같은 무차원 방정식을 유지하는 것이다. 이때 정규화한 농도 곡선은 겹칠 수 있지만 실제 이동 시간과 질량 수송량은 달라질 수 있다. 단위를 바꾸는 일과 물리적으로 닮은 실험을 만드는 일을 구분하며 계산해보자.

그림을 해석하기 전에 모델을 선언한다

길이가 L인 균일한 통로에서 희석된 수동 물질의 정상 일차원 수송을 다룬다. 속도 U와 확산계수 D는 상수이며 D와 L은 양수다. 양 끝 농도는 각각 독립적으로 지정하고 0 이상으로 둔다. 실제 단면적 A도 양수다. 반응·체적 원천·온도 결합·농도가 속도를 바꾸는 효과는 제외한다. 일차원 가정은 단면 방향의 구조를 제외하거나 적절하게 평균했다는 뜻이다.

$$Uc_x=Dc_{xx},\qquad 0\le x\le L,\qquad c(0)=c_L,\quad c(L)=c_R.$$

농도 단위는 kg/m³, 속도는 m/s, 확산계수는 m²/s다. 면적당 질량 플럭스는 kg/(m² s)이며 여기에 A를 곱해야 실제 질량 유량 kg/s가 된다. 정상 상태는 저장량이 시간에 따라 변하지 않는다는 뜻이다. 통로에 물질이 없거나 이동이 멈췄다는 뜻은 아니다.

무차원 방정식을 직접 유도한다

두 경계 농도가 다르면 길이 L과 농도 차이를 기준 척도로 삼는다. 연쇄법칙을 적용하면 공간 미분을 한 번 할 때마다 1/L이 생긴다. 원래 방정식의 공통 계수를 나누면 다음 결과를 얻는다.

$$\xi=\frac{x}{L},\qquad \theta=\frac{c-c_L}{c_R-c_L},\qquad \mathrm{Pe}\,\theta_\xi=\theta_{\xi\xi},\qquad \mathrm{Pe}=\frac{UL}{D}.$$

경계조건은 θ(0)=0, θ(1)=1이 된다. 여기서는 U의 방향을 보존하므로 Pe에도 부호가 있다. 크기는 두 수송 효과를 비교하고 부호는 흐름 방향을 기록한다. 속력으로 양의 Pe를 정의하고 방향을 따로 표시하는 교재의 관례와 혼동하지 않으면 된다.

두 경계 농도가 같아도 물리 해가 사라지는 것은 아니다. 농도는 일정하고 전체 플럭스는 Uc다. 다만 농도 차이로 나눈 측정 정규화는 0/0이다. 경계 농도 차이를 1로 선언한 기저 함수를 별도로 제시할 수는 있지만, 그것을 상수장의 측정 정규화라고 부르면 안 된다.

시간 척도에서 비의 의미를 복원한다

확산이 길이 L을 고르게 만드는 특성 시간은 L²/D이고 이류가 그 길이를 운반하는 시간은 L/|U|다. 모든 입자가 정확히 그 시간에 도착한다는 뜻이 아니라 두 기작의 대표 시간을 비교하는 것이다.

$$t_D=\frac{L^2}{D},\qquad t_A=\frac{L}{|U|}\quad(U\ne0),\qquad |\mathrm{Pe}|=\frac{t_D}{t_A}.$$

|Pe|가 크면 통로 전체를 확산으로 평탄하게 만들기 전에 이류가 물질을 운반한다. 그러나 경계 가까이에서는 확산이 작용할 거리가 L보다 훨씬 짧다. 전체 길이에서 확산이 약하다는 말은 얇은 경계층에서도 확산을 지워도 된다는 뜻이 아니다.

U=0이면 유한한 이류 횡단 시간을 부여하지 않는다. 확산 시간은 그대로 사용할 수 있다. 시간 의존식을 복원하고 τ=Dt/L²를 택하면 cτ+Pe cξ=cξξ가 된다. 이류 시계를 택했다가 U=0에서 식이 특이해지는 것은 척도 선택의 문제이며 확산 방정식 자체의 실패가 아니다.

두 경계조건을 만족하는 해를 구한다

먼저 농도 기울기를 미지수로 놓으면 일차 선형 미분방정식이 된다. 이를 적분한 뒤 다시 농도를 적분하고 양 끝의 농도를 대입한다. U가 비영이면 지수 해, 정확히 0이면 원래 식에서 선형 해를 얻는다.

$$\theta(\xi)=\frac{e^{\mathrm{Pe}\xi}-1}{e^{\mathrm{Pe}}-1},\qquad c=c_L+(c_R-c_L)\theta,\qquad \mathrm{Pe}=0:\ \theta=\xi.$$

양의 Pe에서는 통로 대부분이 왼쪽 경계 농도에 가깝고 오른쪽 경계조건을 하류의 얇은 층에서 맞춘다. 속도를 반대로 하고 두 경계 농도를 교환하면 곡선을 중앙에 대해 반사한 해가 된다. 경계층의 대표 두께 D/|U|가 L보다 충분히 작을 때 얇은 층이라는 해석이 유용하다. 이 두께가 날카로운 물질 경계가 되는 것은 아니다.

D를 정확히 0으로 놓으면 방정식의 차수가 바뀐다. 일차 수송식에는 임의의 두 끝 농도를 동시에 강제할 수 없는 경우가 많다. 바깥쪽에서 평평해 보이는 해만으로 원래 유한 확산 문제와 두 경계조건을 모두 대체할 수 없다.

큰 Pe에서도 작은 꼬리를 보존한다

지수 함수를 그대로 계산하면 최종 비는 유한해도 분자와 분모가 먼저 넘칠 수 있다. P=Pe>0에서는 양의 경계 가중치 두 개를 사용한다.

$$w_R=e^{-P(1-\xi)}\frac{1-e^{-P\xi}}{1-e^{-P}},\qquad w_L=\frac{1-e^{-P(1-\xi)}}{1-e^{-P}},\qquad c=c_Lw_L+c_Rw_R.$$

정확한 가중치의 합은 1이다. 따로 반올림한 가중치를 강제로 재정규화하여 원래 장을 바꾸지는 않는다. 음의 Pe에서는 Q=−Pe>0을 써서 오른쪽 가중치를 ξ에 대한 감쇠 없는 비로, 왼쪽 가중치를 exp(−Qξ)가 곱해진 비로 계산한다. 가까운 지수의 차이는 expm1로 계산하여 정확도를 보완한다.

작은 양의 가중치가 부동소수점의 최소 양수보다 작아질 수 있다. 이때 일반 숫자 표시는 불가하지만 양의 부호와 로그는 남는다. 이를 실제 0인 경계값이나 확산이 없는 구간으로 바꾸어 표시하면 물리 해석이 달라진다.

실제 플럭스와 상쇄를 점검한다

정상 방정식은 전체 플럭스가 위치에 무관함을 뜻한다. 이류와 확산의 개별 기여는 위치에 따라 크게 바뀔 수 있으며 한 항만 검사하면 보존을 잘못 판단할 수 있다.

$$j=Uc-Dc_x,\qquad j_x=0,\qquad j=U\frac{c_L-c_Re^{-\mathrm{Pe}}}{1-e^{-\mathrm{Pe}}}\quad(\mathrm{Pe}>0),\qquad U=0:\ j=-\frac{D(c_R-c_L)}L.$$

cL=0, cR>0, Pe가 큰 양수라면 전체 플럭스는 작은 음수다. 내부의 어떤 위치에서는 이류와 확산의 큰 항을 직접 빼면서 이 작은 역방향 플럭스를 잃는다. 안정적인 경계 표현으로 계산한 값과 반올림된 직접 뺄셈을 구분해서 제시해야 한다. 같은 반올림 항으로 만든 잔차가 0이라는 사실만으로 독립적인 정확도 검증을 대신할 수 없다.

단위를 가진 손계산 예제

U=0.5m/s, L=2m, D=0.1m²/s이면 Pe=10, 확산 시간은 40s, 이류 시간은 4s다. 왼쪽 농도 1kg/m³, 오른쪽 농도 0.2kg/m³에서 전체 플럭스는 약 0.50001816kg/(m² s)다. 단면적이 1.5m²이면 실제 질량 유량은 약 0.75002724kg/s다.

오른쪽 끝에서 이류 기여는 0.1kg/(m² s)에 불과하다. 음의 농도 기울기가 약 0.40001816의 양의 확산 기여를 더하여 같은 전체 플럭스를 만든다. Uc만 비교하면 통로 내부에서 물질이 사라진다고 오해할 수 있다. 통로 안에 임의의 소모 원천을 추가할 필요는 없다.

Pe가 작을 때도 비영 보정을 남긴다

Pe가 0에 가까워지면 곡선은 연속적으로 직선에 접근한다. 그러나 내부의 점에서 속도가 작지만 비영이면 직선에 대한 첫 보정도 비영이다.

$$\theta=\xi+\frac{\mathrm{Pe}}2\xi(\xi-1)+O(\mathrm{Pe}^2).$$

중앙점의 선도 보정은 −Pe/8이다. 이 값을 0.5에 더한 결과가 다시 0.5로 반올림되어도 보정 자체는 별도로 계산해 남길 수 있다. 그림이 직선처럼 보인다는 이유로 속도가 정확히 0이라고 결론 내리지 않는다. 반대로 양 끝에서는 지정한 농도를 만족해야 하므로 보정이 정확히 0이다.

닮은 실험과 수치 격자를 구분한다

양의 배율 a,b에 대해 다음처럼 다른 실험을 만들면 Pe와 같은 ξ에서의 기저 곡선이 보존된다.

$$L'=aL,\qquad D'=bD,\qquad U'=\frac baU,\qquad t_D'=\frac{a^2}b t_D,\qquad j'=\frac ba j.$$

단면적도 A′=a²A로 닮게 바꾸면 전체 질량 유량은 ab배다. 곡선이 겹쳐도 실제 시간·길이·수송량이 같지는 않다. 반면 같은 실험을 cm 단위로 다시 표시할 때는 모든 단위를 일관되게 변환할 뿐이며 실제 관측량은 그대로다.

간격 h인 격자에서는 흔히 Peh=|U|h/(2D)를 사용한다. 이것은 통로 전체가 아니라 한 셀에서 수송 항을 비교한다. 물리 Pe가 커도 격자를 충분히 촘촘하게 하면 격자 Pe는 작아진다. 경계층을 해상하는 일과 물리 수송 상태를 바꾸는 일은 다르다. 여기의 정확한 연속 해 자체가 수치 오차 추정기는 아니다.

다른 무차원 수와 AI 모델의 전이

상수 밀도 Newtonian 운동량식에서 Reynolds 수는 관성과 점성을 비교한다. 동점성계수 ν, 물질 확산계수 D, 열확산율 α를 쓰면 Re=|U|L/ν, Sc=ν/D, Pr=ν/α다. 같은 속력과 길이 척도를 쓸 때 물질 Pe=Re Sc, 열 Pe=Re Pr이다. D와 α가 같은 차원을 가진다는 이유로 서로 바꾸면 안 된다.

아래 표는 양의 기준 속력 V, 밀도 ρ, 점성계수 μ, 중력 가속도 g, 음속 cs를 사용한다. 기준 속력은 선택한 척도이며 반드시 그 순간의 국소 속도와 같지는 않다.

그룹	정의	비교와 한계
Reynolds	Re = ρ V L / μ	관성과 Newtonian 점성의 비교다. 다른 유변 물질에는 그에 맞는 수지가 필요하다.
Mach	Ma = V / cs	속력과 음향 속도의 비교다. 작은 Ma 하나가 모든 강제 조건의 비압축성을 증명하지는 않는다.
Froude	Fr = V / sqrt(g L)	관성과 중력 척도의 비교다. 관성 척도로 나눈 중력 항의 계수는 1/Fr²다.
Prandtl	Pr = ν / α	선언한 물성에서 운동량 확산과 열확산의 비다.
Schmidt	Sc = ν / D	운동량 확산과 선택한 물질 확산의 비다.
Peclet	Pe의 크기 = V L / D	전체 길이에서 이류와 scalar 확산을 비교한다. 이 표에서는 방향을 따로 둔다.
Damköhler	Da = k L² / D	다음 평판 가족에서 확산 시간과 일차 반응 시간을 비교한다.

관성 압력 척도 ρV²와 이류 시계 L/V에는 양의 기준 V가 필요하다. 관측한 영 속도를 그대로 넣어 정규화가 특이해져도 물리 방정식 자체가 무효가 되는 것은 아니다. Pe만 맞춘다고 압축성·반응·기하·모든 경계조건이 맞지는 않는다. AI 대리모델의 입력을 설계할 때도 어떤 조합이 실험 사이의 전이를 뒷받침하고 어떤 생략 기작이 전이를 깨뜨리는지 구분해야 한다.

풀이가 있는 여섯 문제

문제 1 — 단위 변환

손계산 예제를 cm와 s로 표시하면 Pe가 바뀌는가?

풀이: U=50cm/s, L=200cm, D=1000cm²/s이므로 비는 여전히 10이다. 플럭스 숫자도 질량·면적 단위에 맞게 바꾸어야 한다. 실제 수송은 바뀌지 않는다.

문제 2 — 흐름 반전

속도를 반대로 하고 두 경계 농도를 교환하면 어떻게 되는가?

풀이: 새로운 ξ의 농도는 원래 1−ξ의 농도다. 기울기와 전체 플럭스의 부호는 반대가 되고 시간 척도의 크기는 같다.

문제 3 — 일정한 농도

두 끝의 비영 농도가 같으면 플럭스가 반드시 0인가?

풀이: 아니다. 농도가 일정하므로 확산 기여는 0이지만 전체 플럭스 Uc는 남을 수 있다. 대비로 나눈 측정 정규화는 미정이며 선언한 기저 곡선과 구별한다.

문제 4 — 닮음과 실제 유량

a=2, b=1로 놓고 면적도 기하학적으로 닮게 바꾸면 무엇이 달라지는가?

풀이: 속도와 플럭스 밀도는 절반, 확산 시간은 4배, 면적은 4배, 전체 질량 유량은 2배가 된다. 경계 데이터도 대응하면 정규화 곡선은 같다.

문제 5 — 직선처럼 보이는 곡선

Pe=10⁻²⁰에서 중앙값이 0.5로 반올림되면 속도가 0이라고 보고해도 되는가?

풀이: 안 된다. 선도 보정은 −1.25×10⁻²¹이다. 반올림된 장과 안정적인 작은 보정을 따로 표시하고 비영 입력 속도를 보존한다.

문제 6 — 커진 격자

격자 간격을 키우면 물리 Peclet 수가 증가하는가?

풀이: U,L,D가 같으면 물리 Pe는 그대로다. 격자 Pe가 커지고 경계층을 충분히 해상하지 못할 수 있다. 수치 설정 변경과 물리 매개변수 변경은 다르다.

자료와 적용 범위

Langtangen·Linge의 Finite Difference Computing with PDEs는 정지 일차원 이류·확산 문제와 지수 해를 제시한다. 여기서는 실제 계수·부호 있는 속도·두 경계 농도·단면적을 복원하여 플럭스와 닮음 변환을 직접 유도했다. 이 scalar 가족으로 다차원 난류, 농도 의존 물성이나 일반 유체 운동량계를 해결한 것은 아니다.