Perturbative modelling of Arrhenius reactive transport in non-isothermal Gordon– Schowalter Couette–Poiseuille flow

기술 요약: 비등온 Gordon–Schowalter Couette–Poiseuille 유동에서의 Arrhenius 반응성 수송에 대한 섭동 모델링

문제 정의
본 연구는 제1차 균일 Arrhenius 반응을 겪는 점탄성 유체의 평면 Couette–Poiseuille 유동에서 발생하는 결합된 열 및 질량 전달 문제를 다룬다. 구체적인 물리적 시스템은 거리 $H_0$로 떨어진 두 평행판 사이를 흐르는 Gordon–Schowalter 유체를 포함한다. 판들은 동일한 속도 $V$로 움직이며(판의 프레임 내에서 순 Couette 전단력을 제거함), 유동은 외부 압력 구배에 의해 전적으로 구동된다. 벽면은 온도 $T_0$인 등온 상태이며 완벽하게 반응하는 조건(Dirichlet 경계 조건)을 갖는다. 유체는 분산된 열원으로 작용하는 점성 소산(viscous dissipation)을 나타내며, 반응 속도는 온도에 지수적으로 민감하게 반응한다. 분석은 작은 Brinkman 수($Br \ll 1$)와 약한 점탄성 유동($We < 1$) 영역에서 수행되며, 점탄성, 열 소산, 화학 반응성이 결합된 "삼중 결합(triple coupling)"—기존 문헌에서 다뤄지지 않은 시나리오—에 초점을 맞춘다.

방법론
저자들은 합성 파라미터 $\varepsilon \equiv (1-a^2)We^2$($a$는 이송 미분 파라미터, $a \in [-1, 1]$, $We$는 Weissenberg 수)를 기반으로 한 섭동법을 사용한다.

1. 구성 방정식 모델링: Gordon–Schowalter 구성 방정식을 사용하며, 이는 정적 단순 전단(steady simple shear)에서 폐쇄형 전단 응력 해를 허용한다. 전단 점도는 $|a| < 1$일 때 속도 저하(rate-thinning) 특성을 보이지만, 상향 및 하향 이송 Maxwell 극한($a = \pm 1$)에서는 일정하다.
2. 섭동 계층 구조: 속도, 온도, 농도장은 $\varepsilon$과 $Br$의 거듭제곱에 따라 전개된다. 분석은 세 가지 뚜렷한 물리적 기여를 분리한다:
   • 기본 뉴턴 역학적 반응성 Graetz 장(field).
   • 뉴턴 유체의 열–반응 수정항(점성 가열이 포함된 뉴턴 유체).
   • 혼합 점탄성 수정항: 이는 탄성, 소산, 그리고 운동학적 민감성이 동시에 존재해야 하는 항이다.
3. 해석적 해: 속도 수정항은 고정된 응력에서 구성 법칙을 역산하여 도출된다. 온도장은 에너지 평형의 폐쇄형 적분을 통해 얻어진다. 반응성 장은 Sturm–Liouville 고유값 문제로 취급되며, 여기서 점탄성 섭동은 축 방향 감쇠율 및 전환 길이를 결정하는 기본 고유값($\lambda_1^2$)의 이동을 유발한다.
4. 수치적 검증: 해석적 결과는 슈팅 방법(shooting methods)과 유한 차분 이산화를 사용하여 에너지 평형 및 전체 반응성 수송 방정식의 직접 수치 해와 비교 검증된다.

주요 기여 및 결과

1. 혼합 그룹의 특징(The Mixed Group Signature): 열–반응 응답에 대한 주요 점탄성 수정항은 혼합 그룹 $Br \cdot We^2 (1-a^2)$에 비례한다. 이 항은 전단 점도가 일정하게 유지되는 상향 이송 Maxwell($a=1$) 및 하향 이송 Maxwell($a=-1$) 극한에서 동일하게 0이 된다. 효과는 코로테이션(corotational, $a=0$) 케이스에서 최대가 된다. 결과적으로, 이송 미분 파라미터 $a$는 단순히 응력장에만 머무는 것이 아니라, 온도 상승, 축 방향 감쇠율, 그리고 전환 길이에 측정 가능한 흔적(signature)을 남긴다.

2. 이중 경쟁 경로(Dual Competing Routes): 점탄성 섭동은 서로 반대되는 부호를 가진 두 가지 가산적인 1차 경로를 통해 농도장에 영향을 미친다:

   • 이송 경로(Advective Route): 고정된 압력 구배에서, 속도 저하 응답은 관통 유속을 증가시킨다. 이는 반응물이 소비되기 전에 반응물을 더 멀리 운반하여, 전환 길이를 증가시킨다(고유값의 음의 이동).
   • 열–반응 경로(Thermal–Reactive Route): 동일한 속도 저하 현상은 고정된 구동 응력 하에서 전단율을 높이고, 이에 따라 점성 가열을 증가시킨다. 이는 온도를 높여 Arrhenius 운동학을 가속화하고, 전환 길이를 감소시킨다(고유값의 양의 이동).

3. 교차수(The Crossover Number, $K$): 이러한 경로 간의 경쟁은 단일 무차원 교차수 $K$에 의해 지배된다.

   • 만약 $K < 1$이면, 이송 경로가 지배적이다: 즉, 점탄성이 뉴턴 설계보다 더 긴 반응기 길이를 요구하게 된다.
   • 만 만약 $K > 1$이면, 열–반응 경로가 지배적이다: 즉, 점성 가열이 반응을 충분히 가속하여 필요한 반응기 길이를 단축시킨다.
   • 반응기 크기에 미치는 순 점탄성 효과의 부호는 $K=1$에서 정확히 반전된다.

4. 정량적 정확도: 속도 및 온도 수정항에 대한 폐쇄형 표현식이 도출되었다. 1차 섭동 이론이 예측하는 고유값 이동은 직접 수치 해와 비교했을 때 1.9% 미만의 상대 오차를 보였다. 두 경로의 가산성 또한 0.13% 이내의 오차로 검증되었다.

의의 및 주장
본 논문은 열적 결합이 존재할 때 점탄성 흔적이 반응성 수송에 어떻게 나타나는지에 대한 최초의 구조적 분석을 제공한다고 주장한다.

• 구조적 통찰: 본 연구는 표준 Maxwell 모델($a=\pm 1$)이 직선형 전단 유동에서 특정 점탄성 흔적에 대해 "구조적으로 맹목적(structurally blind)"임을 입증한다. 이는 법선 응력이 해당 운동학에서 일을 수행하지 않기 때문이다. 오직 혼합 이송 미분을 가진 유체($|a| < 1$)만이 이 효과를 나타낸다.
• 공학적 유용성: 본 연구는 반응성 코팅, 폴리머 가공, 채널 반응기 응용 분야를 위한 간결한 설계 규칙을 제공한다. 도출된 교차수 $K$를 통해 엔지니어는 점탄성 유체를 다룰 때 뉴턴 기반 설계가 반응기 길이를 과소 예측할지 혹은 과대 예측할지를 예측할 수 있다.
• 설계 맵: $(Da_{II}, Br \gamma_A)$ 평면상의 파라메트릭 맵을 제시하여 점탄성 수정이 무시할 만한 수준인지, 이송 지배적인지, 혹은 운동학 지배적인지를 구분하는 영역을 나타낸다.

저자들은 자신들의 모델이 점성 가열, 반응 운동학, 그리고 탄성이 상호작용하는 시스템을 위한 완전히 재현 가능한 벤치마크 역할을 한다고 결론지으며, 효과의 부호 반전은 Arrhenius 항의 엄격한 선형화 한계를 벗어나서도 나타나는 견고한 특징임을 명시하였다.