Adaptation of the hybrid fictitious domain-immersed boundary method for Reynolds-averaged turbulence modeling

기술 요약: RANS 모델링을 위한 하이브리드 가상 도메인-함몰 경계법(HFD-IBM)의 적응

문제 정의
공학적 응용 분야에서는 유체 부품의 형상 또는 위상 최적화(TO)가 빈번하게 요구되며, 이때 전산유체역학(CFD)이 비용 함수 평가의 기초로 사용된다. CFD 기반의 TO 프레임워크에서 발생하는 주요 병목 현상은 형상이 변화함에 따라 빈번한 재격자 생성(remeshing)이 필요하다는 점이다. 함몰 경계법(Immersed Boundary Methods, IBM)은 단순하고 비정합적인 격자를 사용하고 스칼라 지시 함수(indicator fields)를 통해 기하 구조를 표현함으로써 해결책을 제시하지만, 레이놀즈 평균 나비에-스토크스(RANS) 난류 모델과의 통합에는 여전히 어려움이 있다. 주요 난제는 이방성 격자 세분화 없이 경계층을 모델링하는 것, 높은 레이놀즈 수에서 압력 및 벽면 전단 응력의 가짜 진동(spurious oscillations)을 방지하는 것, 그리고 모멘텀과 난류 변수 간의 긴밀한 결합을 보장하여 불안정성을 피하는 것이다. 또한, 기존의 IBM-RANS 솔버들은 오픈 소스가 아니거나 표준적인 CFD 관행을 결여하고 있는 경우가 많다.

방법론
저자들은 정상 상태(steady-state)의 비압축성 RANS 시뮬레이션을 위해 하이브리드 가상 도메인-함몰 경계(HFDIB) 방법을 적응시킨 오픈 소스 솔버인 openHFDIBRANS를 제안한다. 이 방법론은 OpenFOAM 라이브러리를 기반으로 구축되었으며, 다음과 같은 핵심 요소들을 통합한다:

• 지배 방정식: 솔버는 네 가지 2-방정식 난류 모델($k$-$\epsilon$, $k$-$\omega$, $k$-$\omega$ SST, realizable $k$-$\epsilon$)과 결합된 정상 상태 RANS 방정식을 다룬다. 방정식은 유체 영역으로 해를 제한하기 위해 위상 지시 함수 $\lambda$에 의해 수정되며, 경계 조건을 강제하기 위한 소스 항($f_{ib}$ 및 $S_{ib}$)을 포함한다.
• 마스킹 및 셀 분류: 계산 도메인은 $\lambda$ 필드에 기초하여 고체 내부 셀($\Omega^h_{ic}$), 경계 셀($\Omega^h_{bc}$), 자유 흐름 셀($\Omega^h_{fc}$)로 나뉜다. 마스킹 필드($\alpha_u, \alpha_k, \alpha_\omega$)는 관련 영역에 대한 소스 항의 적용을 제한한다.
• 벽 모델링 및 경계 처리: 경계층을 직접 분해(resolve)하지 않고 처리하기 위해 벽 함수(wall functions)를 사용한다.
  • 마찰 속도: 모멘텀과 난류 사이의 안정성과 긴밀한 결합을 보장하기 위해, 마찰 속도($u_\tau$)는 국부 속도 프로파일 대신 난류 변수를 사용하여 계산된다.
  • 부과된 필드(Imposed Fields): 경계층이 분해되지 않는 경계 셀($y^+_{ib} > y^+_{lam}$)의 경우, 속도는 보간을 통해 재구성되지 않고 오직 난란 점성(turbulent viscosity)에 의해서만 교정된다. 분해된 영역($y^+_{ib} \leq y^+_{lam}$)에서는 벽 함수를 기반으로 속도와 난류 운동 에너지($k$)를 부과하기 위해 다항식 보간이 사용된다.
  • 난류 변수: 특정 소멸률($\omega$)은 벽 근처에서 발산하기 때문에 $\omega$에 대한 특정 소스 항을 추가하지 않는다. 대신, 이산화 행렬을 조작함으로써 영향을 받는 셀에서 $\omega$와 $k$의 생성량($P_k$)을 직접 할당하며, 이는 효과적으로 디리클레(Dirichlet) 유사 조건을 강제한다.
• 솔루션 전략: 솔버는 확장된 SIMPLE 알고리 アル고리즘을 활용한다. 결정적으로, 모멘텀 예측기 이전에 부과된 필드가 계산된다. 모멘텀 방정식이 풀린 후에는 경계 조건을 직접 강제하기 위해 교정 단계($u = u + \alpha_u(u_{ib} - u)$)가 적용된다. 이 단계는 언더 릴랙세이션(under-relaxation)이 필요한 정상 상태 시뮬레이션에서 수렴을 위해 필수적이다. 비압축성 조건은 IBM 강제력과의 충돌을 피하기 위해 유체 영역에서만 강제된다.

주요 기여

1. 오픈 소스 구현: 본 연구는 널리 사용되는 OpenFOAM 프레임워크 내에서 표준 CFD 관행을 포함하는 IBM-RANS 솔버의 완전한 오픈 소스 구현을 제공한다.
2. 안정성 중심 설계: 저자들은 난류 변수를 통해 마찰 속도를 평가하고 모멘텀과 난류 수송 방정식 간의 긴밀한 연결을 강제함으로써 일반적인 IBM-RANS 불안정성 문제를 해결한다.
3. 최적화를 위한 견고성: 이 방법은 기하학적 최적화 응용 분야를 위해 특별히 설계되었으며, 재격자 생성 없이 다양한 레이놀즈 수($10^1$ ~ $10^6$)와 복잡한 기하 구조에 대해 견고함을 입증하였다.
4. 포괄적 벤치마킹: 솔버는 여섯 가지 시나리오(2D/3D 파이프 유동, 역류 단계(backward-facing step), 실린더 주위 유동, NACA-0009 에어포일, Ahmed body)에 대해 표준 체형 적합(body-fitted) 시뮬레이션(simpleFoam) 및 실험 데이터와 검증되었다.

결과
openHFDIBRANS의 성능은 표준 체형 CFD(simpleFoam) 및 실험 데이터와 비교 평가되었다:

• 정확도: 솔버는 넓은 범위의 레이놀즈 수와 서로 다른 난류 모델에 대해 simpleFoam과 일치하는 결과를 생성했다. 속도, 난란 점성 및 레이놀즈 응력에 대한 $L_2$ 노름 오차는 일반적으로 낮았으나, 난란 점성의 급격한 변화가 있는 영역(예: $k$-$\omega$ SST 및 realizable $k$-$\epsilon$ 모델)에서는 약간 더 높았다.
• 민감도: 이 방법은 격자 면에 대한 고체-유체 계면의 위치와 관련하여 만족스러운 견고성을 보였다. 격자 면과 정렬될 때 정확도가 향상되었으나, 솔버는 인터페이스가 어긋난 경우에도 기능적으로 작동하였다.
• 항력 및 양력: 실린더와 NACA-0009 에어포일의 경우, 솔버는 유동 패턴을 포착하였으며 참조 데이터와 잘 일치하는 양력 및 항력 계수를 예측하였다.
• 한계점: 특정 사례(예: 역류 단계 및 NACA-0009)에서 솔버는 고체 객체 뒤쪽의 웨이크(wake) 영역 크기를 약간 과대평가하는 경향이 있었다. 또한, 다른 RAS 모델과 마찬가지로 더 발전된 2-층 모델이 필요한 실린더의 항력 급변(drag crisis) 현상을 해결하지 못했다.

의의 및 주장
본 논문은 openHFDIBRANS가 기하학적 또는 위상 최적화가 필요한 공학 응용 분야를 위한 실행 가능한 도구라고 주장한다. 저자들은 이 방법이 웨이크 크기를 약간 과대평가할 수 있지만, 그 동작이 예측 가능하고 다양한 기하 구조에 걸쳐 일관적이라는 점을 강조한다. 이러한 일관성은 최적화 프레임워크에서 가장 중요한 요소로 식별된다. 본 연구는 표준 체형 방식의 대안으로서, 빈번한 재격사 생성의 계산 비용을 제거하면서도 공학적 설계를 위한 충분한 정확도를 유지하는 준비된 도구로서 이 솔버를 위치시킨다. 저자들이 언급한 향래의 응용 분야로는 바이오센서 및 다이 캐스팅 금형의 위상 최적화가 있다.