A numerical study to analyze the interplay of Weissenberg number and viscosity ratio in a log-strain tensorial model for viscoelastic fluids

기술 요약: 점탄성 유체를 위한 로그-변형 텐서 모델의 수치 해석

문제 정의
본 연구는 점탄성 유체의 시뮬레이션과 관련된 계산적 난제, 특히 바이센베르크 수($Wi$)와 점도 비($\xi$) 사이의 상호작용에 초점을 맞춘다. 점도 비는 용매 점도($\eta$)에 대한 유효 점도(고분자 기여분$\eta_p = \kappa\tau_r$)의 비율로 정의된다. 점탄성 유동에서$Wi$의 역할은 잘 알려져 있으나,$\xi$의 독립적인 영향은 종종 간과된다. 또한, 높은$Wi$에서의 수치 시뮬레이션은 변형 텐서의 양의 정치성(positive definiteness) 상실과 높은 변형률 영역에서의 수치적 붕괴를 특징으로 하는 고 바이센베르크 수 문제(High Weissenberg Number Problem, HWNP)에 자주 직면한다. 저자들은 로그-변형 텐서 모델(탄성 응력과 탄성 변형 사이의 로그 관계를 활용함)을 사용하여 이러한 매개변수들을 체계적으로 탐구한다.

방법론
저자들은 로그-변형 모델의 지배 방정식을 풀기 위해 안정화된 혼합 유한 요소 정식화(stabilized mixed finite element formulation)를 채택한다. 방법론은 다음과 같은 핵심 구성 요소를 포함한다:

• 지배 방정식: 이 모델은 $\log B_{el}$을 포함하는 비선형 이완 항이 결합된 상향 유도 미분(upper-convected derivative)을 사용하여 탄성 변형 텐서($B_{el}$)의 진화를 기술한다. 탄성 응력은 $T_{el} = \kappa \log B_{el}$로 정의된다. 이 정식화는 변형 텐서의 대칭성과 양의 정치성을 본질적으로 보존한다.
• 시간 이산화: 1차 후방 오일러(Backward Euler) 방식의 준암시적(semi-implicit) 기법을 사용한다. 운동학적 장(field)과 비선형성을 분리하기 위해, 대류 및 구성 이완 항은 이전 시간 단계의 기지량을 사용하여 선형화된다.
• 이류 처리: 변형 텐서의 진화는 일반화된 리 미분(Lie derivative) 정식화를 통해 처리된다. 프레임 불변성(frame-invariance)을 보장하고 안정성을 높이기 위해 준 라그랑주 특성 매핑(semi-Lagrangian characteristic mapping)이 텐서 필드 업데이트에 사용된다. $B_{el}$의 행렬식(determinant)과 양의 정치성 제약 조건을 엄격히 강제하기 위해 각 증분 단계에서 투영 단계(projection step)가 적용된다.
• 공간 이산화: 공간 영역은 속도와 압력을 위한 동일 차수 연속 선형 라그 Lagrange 요소(P1-P1)와 변형 텐서를 위한 조각별 상수 불연속 갈레르킨(DG0) 요소로 이산화된다. inf-sup 조건을 회피하고 해를 안정화하기 위해 속도-압력 시스템에 변분 다중 스케일(Variational Multiscale, VMS) 프레임워크가 적용된다. 변형 텐서 필드는 시간 이산화 및 리 미분 처리를 통해 안정화되며, 테스트된 영역에서는 추가적인 VMS 안정화가 불필요하다.
• 벤치마크 케이스: 연구는 평행 채널 유동, 실린더 주변 유동, 4:1 평면 수축 유동이라는 세 가지 전형적인 압력 구동 유동을 분석한다. 시뮬레이션은 탄성과 소산의 효과를 격리하기 위해 낮은 레이놀즈 수($Re \ll 1$)에서 수행된다.

주요 기여 및 결과

1. 점도 비($\xi$)의 역할:
   본 연구는 $\xi$가 $Wi$가 증가함에 따라 비뉴턴 거동의 정도를 결정하는 데 있어 결정적인 매개변수임을 입증한다.

• 평행 채널 유동에서 $\xi$가 낮을 때 속도 프로파일은 대략적으로 포물선 형태를 유지한다. $\xi$가 증가함에 따라, 프로파일은 전단 희박(shear-thinning) 유체의 특징인 중심선 근처의 "플러그(plug-like)" 영역을 발달시킨다. 이 플러그 영역은 $Wi$가 증가함에 따라 줄어들지만, 이 현상은$\xi$가 충분히 클 때만 관찰 가능하다.
• 결과는 $Wi$단독으로는 유동 특성을 결정하지 못하며,$\xi$의 크기가 전단 희박의 정도와 그에 따른 뉴턴 유체 프로파일로부터의 편차를 결정한다는 것을 나타낸다.

2. 항력 감소 및 실린더 주변 유동:
   실린더 주변 유동 시뮬레이션은 전단 희박 거동과 일치하는 항력 감소 효과를 보여준다. 항력 계수($C_D$)는 유효 변형률이 증가함에 따라 감소한다. 연구는 높은 $\xi$에 대해 모델이 국부 유동 유형에 강하게 의존함을 강조한다. 이 모델은 단순 전단(simple shear)에서는 전단 희박성을 보이지만, 신장 유동(extensional flow)에서는 속도 의존성을 보이지 않는데, 이는 더 단순한 모델들과 구별되는 특징이다.

3. 4:1 평면 수축:
   수축 유동에서, 모델은 $Wi$가 증가함에 따라 코너 와류(corner vortex)의 크기가 비단조적으로 변화할 것을 예측한다. 특히, 조사된 매개변수 범위 전체에서 이차 립 와류(secondary lip vortices)의 핵 생성은 완전히 나타나지 않는다. 이러한 부재는 재진입 모서리(re-entrant corner) 근처의 높은 변형률이 점성 저항을 감소시켜, Oldroyd-B와 같은 일정 점도 모델에서 흔히 보이는 이차 순환 구역의 형성을 방지하는 전단 희박 모델(예: PTT, Giesekus)의 거동과 일치한다.

4. 일반화 뉴턴 유체(GNF) 모델과의 비교:
   저자들은 로그-변형 모델의 해를 단순 전단에서의 전단 희박 거동에 맞게 점도를 맞춘 GNF 모델과 비교하였다.

• 평행 채널: GNF와 로그-변형 모델은 거의 동일한 속도 프로파일을 산출하였으며(상대 오차 ~2.6%), 이는 단순 전단 지배 유동에서는 속도 의존적 점도만으로도 충분함을 확인시켜 준다.
• 실린더 주변 유동: 두 모델 사이에 상당한 차이가 발생했다. 속도 프로파일은 단순히 스케일링 인자 차이가 아니라 형태 자체가 달랐다. 저자들은 이를 실린더 상류 및 하류의 신장 유동 영역 존재 때문으로 보며, 여기서 GNF 모델의 점도 의존성(전단으로부터 유도된)은 전체 텐서 로그-변형 모델의 소산 메커니즘을 포착하는 데 실패한다.

의의 및 주장
본 논문은 로그-변형 모델이 제안된 안정화된 수치 기법과 결합되어, HWNP의 전형적인 수치적 붕괴 없이 상대적으로 높은 바이센베르크 수($Wi \sim 35$)까지 점탄성 유동을 성공적으로 시뮬레이션한다고 주장한다. 본 연구는 $Wi$에만 집중할 때 흔히 간과되는 미묘한 부분인, 점도 비$\xi$가 비뉴턴 유동 프로파일을 생성하는 독립적이고 중요한 요소임을 강조한다.

또한, GNF 모델과의 비교는 비탄성 근사의 한계를 분명히 한다. GNF 모델은 단순 전단 지배 유동에는 적합할 수 있으나, 로그-변형 모델은 유동 유형(전단 vs 신장)에 따른 필수적인 레올로지 특징을 포착하며, 이는 복잡한 기하학적 구조에서의 유동을 정확하게 예측하는 데 필요하다. 저자들은 변형 텐서의 양의 정치성을 보존하는 로그-변형 모델의 내재적 수학적 구조가 수치적 안정성과 물리적 타당성을 확보하는 데 필수적이라고 결론짓는다.

향후 연구 방향
저자들은 향후 연구가 3차원 유동 및 더 복잡한 기하학적 구조를 처리할 수 있는 계산 방법을 개발하고, 더 높은 레이놀즈 수에서 관성과 탄성의 상호작용을 조사하는 데 집중해야 한다고 제안한다. 또한, 다중 이완 모드 또는 상태 의존적 이완 시간을 포함하도록 구성 모델을 확장할 것을 제안한다.