[논문 리뷰] A Modified Smagorinsky Subgrid Scale Model for the Large Eddy Simulation of Turbulent Flow

Tommy Kunhung Kim|ArXiv.org|2002. 09. 26.

Fluid Dynamics and Turbulent Flows참고 문헌 26인용 수 4

한 줄 요약

이 논문은 대규모 난류 시뮬레이션(Large Eddy Simulation, LES)에서 벽면 근처에서 발생하는 잘 알려진 문제인 과도한 난류 점성도 문제를 해결하기 위해 수정된 스마고르스키 미소스케일(Smagorinsky subgrid scale, SGS) 모델을 제안한다. 유한요소법과 필터링 분석을 활용해 고전적 스마고르스키 모델에 벽면 적응 기능을 통합함으로써, 새로운 모델은 벽면 근처에서 SGS 점성도를 감소시키면서도 계산 효율성을 유지한다. 이는 원래의 스마고르스키 모델을 능가하고, 채널 유동과 후방 경사면 유동에서 동적 SGS 모델의 성능과 맞먹는다.

ABSTRACT

In the field of Large Eddy Simulation, the Smagorinsky subgrid scale model (in some form) is the most commonly accepted and used subgrid scale model. The purpose of this paper is to address the main weakness of the Smagorinsky model, its poor performance near the wall. The goal is to establish a model that corrects the Smagorinsky model near the walls while at the same time minimizing the computational overhead. A version the Dynamic Subgrid Scale model is also incorporated into the finite element code to facilitate comparisons with the new model near the walls. One of the unique characteristics of Large Eddy Simulations as compared to other methods of dealing with turbulent flows is the idea of filtering. In this paper we define what a filter is and also address an issue related to filters; the error that results when the filtering and differential operations are interchanged. This error is studied under the context of the Finite Element Method which allows us to focus on the function being filtered rather than the filter kernal function, which has been the usual approach in studying this error.

연구 동기 및 목표

• 스마고르스키 SGS 모델의 근본적 한계인 고체 벽면 근처에서의 과도한 난류 점성도 문제를 해결한다.
• 계산 오버헤드를 증가시키지 않고 경계 근처의 수치적 소산을 줄이는 벽면 적응형 SGS 모델을 개발한다.
• 기존 LES 시뮬레이션에서 사용되는 유한요소 프레임워크와의 호환성을 유지한다.
• 표준 난류 유동에서 신규 모델을 원래의 스마고르스키 모델 및 동적 SGS 모델과 비교한다.
• 유한요소법의 맥락에서 필터-도함수 교환 오차를 연구하며, 필터 커널의 성질이 아닌 해의 형태에 중점을 둔다.

제안 방법

• 벽면 감쇠 함수를 도입하여 고전적 스마고르스키 SGS 모델을 수정함으로써, 경계 근처에서 난류 점성도를 감소시킨다.
• 유한요소법(Finite Element Method, FEM)을 활용해 국소 해의 형태 정보를 이용함으로써, 특정 조건 하에서 필터-도함수 교환 오차의 두 번째 차수 정확도를 달성한다.
• McCallen의 이전 작업을 확장하여, 수정된 SGS 모델을 Hydra 유한요소 코드에 구현한다.
• 벽면 행동과 수치 안정성의 기준으로 동적 SGS 모델(Germano 등)을 통합한다.
• 스케일 분리와 필터링이 나비에-스토크스 방정식에 미치는 영향을 연구하기 위해 토포-햇 및 가우시안 필터를 적용한다.
• 난류 유동 시뮬레이션에서 반복적 해법의 안정성 및 수렴 문제로 인해 압력 푸아송 방정식에 직접 해법을 사용한다.

실험 결과

연구 질문

• RQ1 수정된 스마고르스키 모델은 원래의 스마고르스키 모델에 비해 벽면 근처에서 난류 점성도를 어떻게 감소시키는가?
• RQ2 유한요소법은 필터-도함수 교환 오차를 어느 정도 감소시키며, 어떤 조건에서 두 번째 차수 정확도를 달성하는가?
• RQ3 수정된 스마고르스키 모델의 성능은 난류 통계 및 수치 안정성 측면에서 동적 SGS 모델과 비교해 어떻게 되는가?
• RQ4 이중 필터링을 요구하거나 음수의 난류 점성도를 도입하지 않고도 수정된 SGS 모델이 벽면 적응 행동을 달성할 수 있는가?
• RQ5 메쉬 해상도(예: 88×32×16)는 채널 및 후방 경사면 유동에서 수정된 모델의 정확도에 어떤 영향을 미치는가?

주요 결과

• 수정된 스마고르스키 모델은 벽면 근처에서 난류 점성도를 성공적으로 감소시켰으며, 벽면 방향 난류 강도 및 레이놀즈 전단 응력 프로파일을 통해 확인되었다.
• 채널 유동 시뮬레이션에서의 난류 점성도 데이터는 수정된 모델이 벽면 경계를 감지하고 이를 바탕으로 점성도를 감쇠시킴을 보여주며, 동적 SGS 모델과 유사한 행동을 보였다.
• 수정된 모델은 채널 유동과 후방 경사면 구성 모두에서 난류 통계를 예측하는 데 원래의 스마고르스키 모델보다 뛰어난 성능을 보였다.
• 유한요소법을 사용할 경우, 해의 국소 형태에 대한 지식 덕분에 특정 조건 하에서 필터-도함수 교환 오차가 두 번째 차수 정확도를 확보함을 입증하였다.
• 동적 SGS 모델은 벽면 근처에서 효과적이지만, 음수의 난류 점성도로 인해 수치적 불안정성이 발생하는 반면, 수정된 스마고르스키 모델은 이를 피할 수 있었다.
• 이중 필터링과 추가적인 동적 계산을 피하기 때문에, 수정된 모델는 낮은 계산 비용으로 동적 SGS 모델과 맞먹는 성능을 달성하였다.