Fractal Scaling of Moffatt Vortices in Triangular Cavity Flow

기술 요약: 삼각 공동 유동 내 모팻 와류의 프랙탈 스케일링

문제 정의
본 연구는 삼각 공동 내의 느리고 점성 있는 비압축성 유동에서 발생하는 코너 와류(모팻 와류, Moffatt vortices)의 형성, 정량화 및 프랙탈 특성을 다룬다. 스토크스 유동(Stokes flow)의 모서리 근처에서 무한한 역회전 와류 시퀀스가 존재한다는 사실은 모팻(Moffatt)과 딘 및 몽타뇽(Dean and Montagnon)에 의해 이론적으로 확립되었고 타네다(Taneda)에 의해 실험적으로 시각화되었으나, 삼각 기하 구조에서의 이러한 구조에 대한 엄밀한 수치적 특성화 연구는 여전히 제한적이다. 또한, 물리적 통찰을 얻기 위해 면적-둘레 방법(area-perimeter method)을 사용하여 이러한 특정 와류 시퀀스의 프랙탈 성질을 조사한 연구는 이전에 수행되지 않았다. 저자들은 또한 이러한 자기 유사 구조가 스토크스 영역을 넘어 다양한 기하학적 구조(삼각형 대 사각형)와 스케일에서도 지속되는지 확인하고자 한다.

연구 방법론
본 연구는 상용 솔버인 ANSYS Fluent(Release 21.1)를 이용한 유한체적법 기반의 수치적 접근 방식을 채택한다. 비정상(non-dimensional) 형태의 나비에-스토크스 방정식을 사용하여 등변 삼각 공동(높이가 밑변의 두 배인 형태) 내부의 정상 상태 2차원 비압축성 유동을 해결한다.

• 수치 스킴: 압력 기반 결합 솔버(pressure-based coupled solver)를 사용한다. 대류 항은 2차 업윈드(upwind) 스킴으로, 확산 항은 중앙 차분(central difference) 스킴으로 이산화한다.
• 격자 해상도: 격자 독립성 연구를 수행하기 위해 네 가지 서로 다른 격자 해상도(2,028개에서 127,308개의 요소 범위)를 생성한다. 주요 유동 정량화에는 가장 미세한 격자(요소 크기 0.005)를 사용하며, 격자 스케일링 효과가 프랙탈 특성에 미치는 영향을 평가하기 위해 모든 격자의 데이터를 분석한다.
• 프랙탈 분석: 면적-둘레 방법을 사용하여 와류 시퀀스의 프랙탈 차원($D$)을 계산한다. 와류의 최외곽 유선은 반타원형으로 근사한다. 둘레($P$)는 라마누잔(Ramanujan)의 근사식을 사용하여 추정하고, 면적($A$)은 기하학적으로 계산한다. 프랙탈 차원은 멱법칙 관계인 $P \propto A^{D/2}$로부터 도출된다.
• 검증: 결과는 타네다의 실험적 유동 시각화 및 모팻의 이론적 예측과 질적 및 양적으로 검증된다. 또한 다양한 레이놀즈 수($Re$)에서 사각 공동 유동과의 비교 분석이 수행된다.

주요 결과

1. 와류 정량화: $Re=1$에서, 시뮬레이션은 하단 꼭짓점 근처에서 7개의 코너 와류($V_1$~$V_7$)로 구성된 캐스케이드를 성공적으로 분해한다. 연속적인 와류의 크기($r$)와 강도($\psi$)는 급격히 감소하며, 점근적 비율은$r_{n+1}/r_n \approx 0.49$및$\psi_{n+1}/\psi_n \approx 0.0012$로 나타난다. 이 비율은 모팻의 이론적 예측과 밀접하게 일치한다.
2. 프랙탈 차원: 계산된 와류의 프랙탈 차원은 1과 2 사이에 위치하며, 이는 이들의 비정수적 프랙탈 성질을 확인시켜 준다. 역관계가 관찰되는데, 가장 큰 1차 와류($V_1$)가 가장 작은 프랙탈 차원($D \approx 1.23$)을 가지는 반면, 가장 작은 분해된 와류($V_7$)는 $D \approx 1.78$에 도달한다. 이러한 경향은 와류가 코너로 접근함에 따라 공간적 구속이 증가하는 것에 기인한다.
3. 스케일링 법칙: 1차 와류의 기하학적 스케일링을 기반으로 임의의 연속 와류($n$)의 프랙탈 차원을 추정하는 경험적 관계식을 제안한다:
   $$D(n) \approx 2 \left[ \frac{\log P_1 - (n-1)\log 2}{\log A_1 - 2(n-1)\log 2} \right]$$
   이 관계식은 임의의 격자 해상도에 대해 일반화된다. $\log P$와 $\log A$ 사이에서 강한 선형 상관관계($R^2 \approx 0.99$)가 발견되었으며, 이는 이산적 척도 불변성(discrete scale invariance)을 확인해 준다.
4. 레이놀즈 수의 영향: $Re = 100$및$Re = 500$을 조사한다. 1차 와류의 중심이 오른쪽 벽 쪽으로 이동하고 2차 와류가 확장됨에도 불구하고, 고차 와류($V_3, V_4$ 등)는 그 특유의 형상과 교차하는 와도 부호를 유지한다. 와류 시퀀스의 기하학적 자기 유사성은 이러한 중간 정도의 레이놀즈 수 영역에서도 지속된다.
5. 기하학적 보편성: 사각 공동과의 비교 분석을 통해 사각 기하 구조의 코너 와류 또한 자기 유사 시퀀스와 프랙탈 특성(예: 사각 공동의 1차 코너 와류에 대해 $D \approx 1.56$)을 보임을 밝혀냈다. 이는 프랙탈 조직이 특정 기하학에 국한된 인위적 현상이 아니라, 구속된 점성 유동의 견고한 특징임을 시사한다.

의의 및 주장
본 논문은 삼각 공동 내의 모팻 와류를 엄밀하게 규명함으로써 고전적 유체 토폴로지와 프랙탈 기하학 사이의 간극을 메운다고 주장한다. 주요 기여는 다음과 같다:

• 수치적 검증: 최대 7개의 뚜렷한 에디(eddy)를 분해함으로써 삼각 기하 구조 내 모팻 와류 캐스케이드에 대한 고해상도 수치적 확인을 제공한다.
• 프랙탈 특성화: 이 코너 와류들이 와류의 크기 및 강도와 체계적으로 연결된 비정수 프랙탈 차원을 보유하고 있음을 입증한다.
• 방법론적 프레임워크: 1차 와류의 기하학적 구조를 기반으로 연속 와류의 프랙탈 차원을 추정할 수 있는 실용적인 경험적 관계식을 도입하였으며, 이는 임의의 격자 해상도에 적용 가능하다.
• 견고성: 코너 와류의 자기 유사적, 프랙탈적 성질이 스토크스 유동이나 삼각형 기하에 국한되지 않고, 중간 레이놀즈 수 영역 및 사각 공동에서도 지속됨을 입증한다.

저자들은 면적-둘레 방법이 특히 가장 작고 깊은 곳에 위치한 와류에 대해 격자 민감도를 가진다는 점을 인정하며, 이들은 신뢰할 수 있는 차원 추출을 위해 충분히 미세한 공간 해상도를 필요로 한다. 결론적으로, 본 연구는 더 복잡한 기하 구조 및 더 높은 레이놀즈 수 영역에서의 다중 스케일 와류 구조에 대한 향후 연구를 위한 기초를 제공한다.