Revisiting creeping viscoelastic cross-slot flow: Global linear stability and structural sensitivity analyses

기술 요약: 크리핑 점탄성 크로스 슬롯 흐름의 재고찰

문제 제기
본 논문은 2차원 크리핑 크로스 슬롯 흐름 내 올드로이드-B(Oldroyd-B) 유체의 탄성 불안정성의 물리적 기원에 관한 오랫동안 해결되지 않은 난제를 다룬다. 거의 반세기 전 실험적으로 처음 관찰되고 20년 전 수치적으로 재현된 이 흐름은, 와이센베르크 수($Wi$)가 임계값을 초과함에 따라 정상 대칭 상태에서 정상 비대칭 상태로 전환되는 대칭성 깨짐(symmetry-breaking) 현상을 보인다. 광범위한 현상학적 연구에도 불구하고, 이 불안정성을 유도하는 정확한 메커니즘에 대해서는 여전히 논쟁이 이어지고 있다. 이전의 가설들은 정체점(stagnation point) 상류의 압축 흐름, 하류의 신장 흐름, 날카로운 모서리의 기하학적 특이점, 모서리 근처의 강한 전단, 또는 탄성과 전단 박화(shear-thinning) 사이의 상호작용 등을 원인으로 지목해 왔다. 저자들은 해석적 기저 상태(analytical base state)가 결여된 복잡한 기하 구조에서의 점탄성 유체 시뮬레이션과 관련된 수치적 어려움(고 와이센베르크 수 문제)으로 인해, 불안정성 발생 시점과 모드 구조에 대한 제1원리 기반의 이론적 예측이 부족했음을 지적한다.

방법론
저자들은 영점 레이놀즈 수에서의 표준 2D 크로스 슬롯 채널 내 확산형 올드로이드-B 유체에 대하여, 네 가지 분석 및 수치 기법을 결합한 포괄적이고 다각적인 접근 방식을 채택한다:

1. 직접 수치 시뮬레이션 (DNS): FreeFem++를 이용한 유한 요소법(FEM)을 사용하여, 불압축성 스토크스 방정식(incompressible Stokes equations)과 결합된 확산형 올드로이드-B 구성 방정식을 푼다. 수치적 안정화를 위해 이산 탄성-점성 분리 응력 구배(DEVSS-G) 정식화와 스트림라인 업윈드 페트롤-갤러킨(SUPG) 기법을 활용한다.
2. 전역 선형 안정성 분석 (LSA): 정체된 기저 상태 주변에서 지배 방정식을 선형화하여 복소 고유값과 고유 모드를 구한다. 이를 통해 임계 $Wi$와 섭동의 시간적 성장률을 예측할 수 있다.
3. 구조적 민감도 분석: 직접 및 부가(adjoint) 고유 모드를 모두 계산함으로써, 구조적 민감도 장(velocity-based 및 stress-based)을 정의하고, 국부적 피드백에 가장 민감한 공간적 핵심인 '웨이브메이커(wavemaker)' 영역을 식별한다.
4. 에너지 예산 분석: 섭동 에너지 진화 방정식을 분해하여 에너지 생성 또는 소실에 기여하는 구체적인 물리적 메커니즘을 식별한다. 저자들은 기저 상태와 섭동 사이의 결합 항, 특히 중합체 동력 밀도(polymeric power density)에 집중하여 불안정성의 에너지원을 찾아낸다.

또한, 본 연구는 날카로운 모서리를 제거한 둥근 모서리 변형 모델과 자유 정체점을 제거한 중앙 실린더 삽입 변형 모델을 조사하여, 식별된 메커니즘의 견고성을 테스트한다.

주요 결과

• 분기 및 검증: DNS는 전형적인 초임계 피치포크 분기(supercritical pitchfork bifurcation)를 성공적으로 재현하며, $Wi_{cr} \approx 0.36$에서 흐름이 대칭에서 비대칭으로 전환됨을 보여준다. LSA는 이 임계값을 일관되게 예측하며, DNS 데이터 및 동적 모드 분해(DMD) 결과와 일치하는 시간적 성장률을 제시한다.
• 불안정성 모드 구조: 주요 직접 고유 모드는 카이럴(chiral) 속도-응력 섭동으로 구성된다. 속도 섭동은 정체점 중심에 다이아몬드 형태의 와류 구조를 형성한다. 결정적으로, 법선 방향의 중합체 응력 섭동($\hat{\tau}_{yy}$)은 반대칭 패리티를 가지며, 이는 신장 흐름에 의해 생성된 복굴절 가닥(birefringent strand)을 기울이고 회전시킨다. 이 회전은 흐름의 분리선(separatrices)을 편향시켜 출구 분지(outlet branches) 간의 불균등한 유량 분포를 초래한다.
• 공간적 핵심 (웨이브메이커): 구조적 민감도 분석 결과, 불안정성은 모서리나 정체점 자체에 의해 구동되는 것이 아님이 밝혀졌다. 대신, 응력 기반 웨이브메이커는 정체점 하류의 신장 흐름 내에 존재하는 좁고 높은 신장률을 가진 능선(ridges)에 국부적으로 위치한다. 이 능선들은 중합체 형태 섭동의 '헤링본(herringbone)' 패턴을 보이며, 이는 측면 방향의 인장과 압축이 동시에 발생하는 특징을 갖는다.
• 에너지원: 에너지 예산 분석은 기저 상태와 섭동 사이의 결합이 지배적인 에너지원임을 식별한다. 구체적으로, 섭동 중합체 동력 밀도($P$)는 양수이며 웨이브메이커 분석에서 확인된 것과 동일한 수직 능선을 따라 최대가 된다. 이는 저장된 탄성 에너지가 변형된 중합체 사슬으로부터 섭동 운동 에너지로 국부적으로 전달되어 불안정성을 구동함을 의미한다.
• 기하학적 변형: 둥근 모서리와 중앙 실린더 삽형 모델 모두에서 불안정성이 지속된다. 피크 중합체 동력 밀도는 여전히 신장 흐름 영역의 수직 능선을 따라 국부적으로 나타난다. 이는 날카로운 모서리나 자유 중앙 정체점이 필수적인 불안정화 요소가 아님을 입증한다.

의의 및 주장
본 논문은 크로스 슬롯 탄성 불안정성의 물리적 기원에 관한 논쟁을 해결했다고 주장한다. 저자들은 이 불안정성이 기하학적 특이점(모서리)이나 정체점 그 자체에서 기인하는 것이 아니라, 크로스 슬롯 흐름의 특징인 신장 지배 영역에서의 탄성 에너지 방출에서 비롯된다고 결론짓는다.

구체적인 메커니즘은 다음과 같다:

1. 정체점 하류의 신장 흐름에서 상당한 법선 응력과 구배를 가진 복굴절 가닥을 생성한다.
2. 이 가닥을 기울이고 회전시키는 카이럴 섭동이 발생한다.
3. 높은 신장률을 가진 좁은 능선(웨이브메이커 영역) 내에서 저장된 탄성 에너지가 방출되어 불안정성의 성장을 촉진한다.

저자들은 이러한 결과가 기존의 상충하는 해석들을 화해시킨다고 주장한다. 즉, OSCER 실험이 보여주듯 날카로운 모서리가 필수적이지 않으며, 실린더 삽입 실험이 보여주듯 정체점이 유일한 동력도 아니지만, 핵심 요소는 높은 신장률의 흐름과 그에 따른 중합체 신장이라는 것이다. 본 연구는 해석적 기저 상태를 사용할 수 없는 복잡한 점탄성 흐름에서 전역 안정성, 구조적 민감도, 에너지 예산 분석을 통합하는 것이 불안정성을 이해하는 데 매우 유용함을 입증하였다.