Leveraging unstructured grids for direct numerical simulations of wall turbulence

본 논문은 국소 콜모고로프 척도에 따라 격자 크기를 조정하여 기존 직교 격자와 동등한 정확도를 달성하면서도 고 레이놀즈 수 및 복잡한 리블릿 기하학적 구조에서 계산 비용을 크게 절감하는 벽면 난류의 직접 수치 시뮬레이션을 위한 비정렬 격자 생성 프레임워크인 {\eta}-그리드를 소개한다.

핵심

자동차 위를 흐르는 공기나 선체 옆을 지나가는 물의 흐름을 시뮬레이션해 보라고 상상해 보세요. 이를 컴퓨터로 정확하게 수행하기 위해 과학자들은 **직접 수치 시뮬레이션 (DNS)**이라는 기법을 사용합니다. DNS 를 유체 (공기나 물) 를 수백만 개의 작고 보이지 않는 정육면체 (격자) 로 분해하는 거대한 3 차원 디지털 현미경을 만드는 것이라고 생각하세요. 컴퓨터는 그런 다음 각 정육면체가 어떻게 움직이고 이웃과 상호작용하는지 계산합니다.

문제는 표면 (예: 배의 측면) 근처의 유체 흐름이 매우 혼란스럽고 세밀하다는 점입니다. 선명한 그림을 얻으려면 표면 바로 옆에 이러한 작은 정육면체가 엄청나게 많이 필요합니다. 하지만 표면에서 멀어질수록 혼란은 완화되므로 그렇게 작은 정육면체가 필요하지 않습니다.

구식 방법: "단단한 벽돌 벽"

전통적으로 과학자들은 **직교 격자 (Cartesian grid)**를 사용했습니다. 동일한 단단한 벽돌로 벽을 쌓는다고 상상해 보세요.

• 문제: 표면 근처의 미세한 세부 사항을 보려면 바닥의 벽돌을 매우 작게 만들어야 합니다. 하지만 이러한 벽돌은 단단하고 직선으로 연결되어 있기 때문에, 세부 사항이 중요하지 않은 곳에서도 벽의 꼭대기까지 동일한 작은 벽돌을 사용해야 합니다.
• 결과: 결국 대부분의 벽돌이 불필요한 수조 개의 작은 벽돌로 된 벽이 됩니다. 이로 인해 컴퓨터 시뮬레이션이 매우 느리고 비용이 많이 들게 됩니다. 마치 조수를 측정하기 위해 해변의 모래 알갱이 하나하나를 세어 보려는 것과 같습니다.

새로운 해결책: "똑똑하고 늘어나는 그물"

이 논문은 **$\eta$-격자 (eta-grid)**라는 새로운 방법을 소개합니다. 단단한 벽돌 대신 똑똑하고 늘어나는 낚시 그물을 상상해 보세요.

• 작동 원리: 저자들은 필요한 세부 사항의 양에 따라 그물의 구멍 크기가 자동으로 변하도록 시스템을 설계했습니다.
  • 표면 근처 ("내부 층"): 그물은 유체의 작고 혼란스러운 소용돌이를 잡기 위해 매우 작고 조여진 구멍을 가지고 있습니다.
  • 더 멀리 ("외부 층"): 유체가 더 차분해지면 그물이 자동으로 늘어나 구멍이 훨씬 커집니다.
• 비밀 재료: 이러한 구멍의 크기는 **콜모고로프 척도 (Kolmogorov scale, $\eta$로 표기)**라는 것에 기반합니다. 이를 유체 내 특정 높이에서 존재할 수 있는 "가장 작은 소용돌이"로 생각하세요. 새로운 격자는 단순히 이렇게 말합니다. "이 특정 높이에서 가장 작은 소용돌이를 잡을 만큼만 구멍 크기를 크게 하되, 그 이상은 크게 하지 마라."

이것이 중요한 이유

저자들은 이 "똑똑한 그물"을 스펙트럴 요소법과 유한 체적법과 같은 두 가지 다른 유형의 컴퓨터 코드에서 테스트하고, 구식인 "단단한 벽돌" 방법과 비교했습니다.

1. 정확도: 결과는 거의 동일했습니다. "똑똑한 그물"은 마찰과 속도 같은 주요 측정값에서 1% 미만의 차이로 "단단한 벽돌"만큼 물리학을 정확하게 포착했습니다.
2. 엄청난 절감: 여기서 마법이 일어납니다.
   • 매끄러운 표면 (평평한 벽과 같은) 의 경우, 새로운 격자는 고속에서 필요한 "벽돌" (격자점) 수를 약 90% 줄였습니다.
   • 거친 표면 (항력을 줄이도록 설계된 "리블릿"이라는 작은 홈이 있는 벽과 같은) 의 경우, 절감 효과는 더욱 극적이어서 격자점이 97% 적게 필요했습니다.

"홈이 파진 벽"의 비유

리블릿 부분을 이해하기 위해 골프공의 질감이나 상어 피부처럼 작은 평행 홈으로 덮인 벽을 상상해 보세요.

• 구식 방법: 이를 시뮬레이션하기 위해 단단한 벽돌 방법은 홈이 격자를 전체적으로 미세하게 만들었기 때문에 어디서나 벽돌을 작게 유지해야 했습니다. 이는 직물에서 멀리 떨어진 부분까지 스웨터의 실 하나하나를 세어 보려는 것과 같습니다.
• 신식 방법: "똑똑한 그물"은 그 작은 홈들에서 몇 인치만 위가 되면 흐름이 다시 매끄러워진다는 것을 알고 있습니다. 따라서 그물은 홈 바로 위의 구멍을 즉시 늘려 더 이상 중요하지 않은 미세한 세부 사항을 무시합니다.

결론

저자들은 유체 시뮬레이션을 위한 스마트 줌 렌즈처럼 작동하는 프레임워크를 개발했습니다. 이는 필요한 곳 (벽 근처) 에는 계산 능력을 집중시키고 필요하지 않은 곳에서는 이를 완화합니다.

• 매끄러운 벽의 경우: 시뮬레이션이 커질수록 노력의 증가 속도가 훨씬 느려집니다.
• 거친 벽의 경우: 더욱 효율적으로 확장되어, 이전에는 처리할 수 없었던 복잡한 항력 감소 표면의 시뮬레이션을 컴퓨터에서 가능하게 합니다.

간단히 말해, 그들은 동일한 고품질 작업을 슈퍼컴퓨터가 한 달 걸려 수행할 수 있는 작업을 며칠 만에 완료할 수 있는 수준으로 줄이는 계산 능력의 일부만으로 수행할 수 있는 방법을 찾았습니다.

기술 요약

기술 요약: 벽면 난류의 직접 수치 시뮬레이션을 위한 비정렬 격자 활용

문제 제기
벽면 난류의 직접 수치 시뮬레이션 (DNS) 은 전체 영역에 걸쳐 콜모고로프 길이 척도 ($\eta$) 를 해석해야 한다는 요구 사항으로 인해 계산 비용이 매우 높습니다. 매끄러운 표면과 비매끄러운 표면 (예: 침수 경계 방법을 통한) 에 널리 사용되는 기존 직교 격자 솔버는 일반적으로 전체 영역에 걸쳐 고정되거나 천천히 변화하는 흐름 방향 ($\Delta x^+$) 과 횡방향 ($\Delta z^+$) 격자 해상도를 강제합니다. 이러한 접근법은 콜모고로프 척도가 커지는 외부 영역에서 과도한 해석을 초래하여 계산 병목 현상을 야기합니다. 이는 리블릿과 같은 복잡한 표면을 포함하는 시뮬레이션에서 특히 심각하여, 격자가 미세 홈을 해석하기 위해 표면 근처에서 극도로 정밀해야 하지만, 직교 프레임워크에서는 이 정밀한 해상도가 전체 경계층 두께에 걸쳐 불필요하게 유지됩니다.

방법론
저자들은 $\eta$-격자라고 명명된 새로운 비정렬 격자 생성 프레임워크를 제안합니다. 핵심 원리는 벽면 수직 ($\Delta y^+$) 과 횡방향 ($\Delta z^+$) 격자 크기를 국소 콜모고로프 척도 ($\eta^+$) 에 비례하도록 설정하여, 난류 척도가 더 큰 외부 유동 영역에서 격자를 조밀하게 만드는 것입니다.

1. 격자 구성:

   • 내부 층: 벽면 인접한 얇은 층 (약 50 점성 단위) 은 기존 DNS 와 유사한 점성 스케일 격자 크기 ($0.3 \lesssim \Delta y^+ \lesssim 4$, $\Delta z^+ \simeq 5$) 를 사용합니다. 비매끄러운 표면 (예: 리블릿) 의 경우, 격자는 최소 표면 파장 ($\ell^+$) 을 $\ell^+/30 \lesssim \Delta y^+, \Delta z^+ \lesssim 4$로 해석합니다.
   • 외부 영역: 내부 층 위에서는 격자 크기가 국소 콜모고로프 척도에 비례하여 증가합니다 ($\Delta y^+ \simeq \Delta z^+ \simeq 2\eta^+$).
   • 스케일링 법칙: 저자들은 난류 채널 유동과 제압력 구배 (ZPG) 난류 경계층 (TBL) 에 대한 로그 영역과 웨이크 영역에서의 $\eta^+$에 대한 반경험적 피팅을 정교화했습니다. 그들은 $y^+ \approx \delta^+/2$에서 전환되는 $\eta^+$에 대한 구간별 멱함수 피팅을 제안합니다.
   • 리블릿 적응: 리블릿 기하학의 경우, 이 프레임워크는 2 차 유동과 켈빈 - 헬름홀츠 롤러를 해석하기 위해 정사각형 셀로 채워진 "리블릿 하층"($\delta_r^+$) 을 도입하며, 그 위에서는 $\eta$ 기반 스케일링으로 부드럽게 전환됩니다.

2. 구현:

   • 이 프레임워크는 오픈 소스 패키지 Gmsh를 활용하여 멀티블록 비정렬 격자를 생성합니다.
   • 격자는 SOD2D(스펙트럴 요소 방법) 와 OpenFOAM(유한 체적 방법) 두 가지 다른 솔버를 사용하여 테스트되었습니다.
   • 시뮬레이션은 매끄러운 벽과 다양한 리블릿 기하학 (삼각형 홈) 에 대한 난류 채널 유동 및 ZPG TBL 을 포함하며, 마찰 레이놀즈 수 $\delta_0^+ = 1000$까지 수행되었습니다.

주요 기여

• 프레임워크 개발: 매끄러운 표면과 복잡한 거친 표면 모두에 적용 가능한, 국소 콜모고로프 척도에 동적으로 격자 간격을 적응시키는 견고한 멀티블록 비정렬 격자 생성 전략.
• 정교화된 스케일링 법칙: TBL 과 채널 유동의 웨이크 영역에서 콜모고로프 척도에 대한 반경험적 피팅을 업데이트하여 외부 층의 격자 크기 결정 정확도를 향상시킴.
• 솔버 중립성: $\eta$-격자 프레임워크가 서로 다른 수치 이산화 체계 (SEM 및 FVM) 에서 일관된 고정밀 결과를 산출함을 입증.
• 리블릿 전용 구성: 2 차 유동 및 KH 롤러와 같은 하층 물리를 해석하면서도 외부 유동에서 빠른 격자 조밀화를 가능하게 하는 리블릿을 위한 맞춤형 격자 전략으로, 직교/IBM 접근법의 한계를 극복.

결과

• 정확도: $\eta$-격자는 기존 직교 격자 및 참조 DNS/실험 데이터와 비교하여 1% 미만의 차이로 결과를 산출합니다. 지표에는 마찰 계수 ($C_f$), 평균 속도 분포 ($U^+$), 난류 응력 및 스펙트로그램이 포함됩니다.
• 격자 효율성:
  • 매끄러운 벽: 격자 점 수 ($N_\eta$) 는 $\propto \delta_0^{+2.5}$로 스케일링되는 반면, 쌍곡선 탄젠트 신장을 가진 직교 격자는 $\propto \delta_0^{+3.0}$입니다. $\delta_0^+ = 6000$에서 $\eta$-격자는 표준 직교 격자 격자 점의 약 10% 만 필요로 합니다 ($N_\eta/N_{Tanh} \simeq 0.1$).
  • 리블릿: 효율성 향상은 더욱 두드러집니다. 항력 감소 삼각형 리블릿의 경우 $N_\eta$는 $\propto \delta_0^{+2.0}$으로 스케일링되는 반면, 직교 격자는 $\propto \delta_0^{+3.0}$입니다. $\delta_0^+ = 6000$에서 $\eta$-격자는 직교 격자 격자 점의 약 3% 만 필요로 합니다 ($N_\eta/N_{Tanh} \simeq 0.03$).
• 리블릿 물리: 이 프레임워크는 가상 원점, 2 차 유동, 켈빈 - 헬름홀츠 롤러를 포함한 리블릿 위의 복잡한 유동 메커니즘을 성공적으로 포착하며, 중간 격자 매개변수 ($C_y=2.0, C_z=2.5, y_{in}^+=50$) 를 사용하여 격자 수렴을 달성했습니다.
• TBL 검증: 리블릿 위의 ZPG TBL 에 대한 DNS 는 웨이크 영역의 고유한 실험-DNS 불일치에도 불구하고 항력 감소 백분율 및 속도 분포와 관련하여 실험 데이터 (Baron & Quadrio; Choi & Orchard) 와 좋은 일치를 보였습니다.

의의
본 논문은 $\eta$-격자 프레임워크가 "방대한 격자 절감"을 제공하여 복잡한 표면에서의 고 레이놀즈 수 벽면 난류 DNS 를 계산적으로 실현 가능하게 만든다고 주장합니다. 고정된 해상도를 유지하는 대신 격자 해상도를 국소 물리 척도 (콜모고로프 척도) 에 정렬함으로써, 이 방법은 정확성을 희생하지 않으면서 계산 비용 (격자 점 및 시간 단계 제약) 을 극적으로 감소시킵니다. 저자들은 이 접근법이 이전에 직교 기반 방법의 prohibitive 한 비용으로 인해 방해받았던 실험 조건과 일치하는 레이놀즈 수에서 복잡한 표면 위의 TBL 에 대한 DNS 수행을 경제적으로 가능하게 한다고 강조합니다. 이 연구는 구조화된 직교 메시의 제약을 넘어 고정밀 난류 시뮬레이션을 위한 비정렬 격자의 유효성을 검증합니다.