로렌츠 힘이 작용하는 정사각 덕트 내 난류유동 해석

Prediction of Turbulent Flows in a Square Duct Subjected to Lorentz Force

본 연구는 정사각형 단면을 갖는 덕트 내 완전 발달된 난류 자기유체역학(magnetohydrodynamic, MHD) 유동에서 자기장이 평균속도와 난류 특성에 미치는 영향을 이해하기 위하여 수행되었다. 결합된 형태의 나비에-스토크스(Navier-Stokes) MHD 방정식들을 기반으로 로렌츠 힘 효과를 포함한 확장된 형태의 레이놀즈응력 수송 모형을 적용하였으며, 이 해석은 전산 유체역학을 위한 자체 개발 코드로 수행되었다. 난류 MHD 유동에 로렌츠 힘 효과를 직접적으로 부여하기 위하여 레이놀즈응력과 소멸률 방정식들에 로렌츠 힘 관련 생성과 소멸에 관한 모형들을 추가하였으며, 레이놀즈응력 방정식의 속도-압력구배 상관관계 모형의 수정을 통해서도 간접적인효과를 발휘하도록 고려하였다. 또한, 본 연구에서는 타원-혼합 레이놀즈응력 모형과 연계되도록 로렌츠 힘과 관련된 모형계수들을완전 발달된 평판 MHD 자료들을 활용하여 조율하는 과정을 거친 후에 사각덕트 유동을 해석하였다. 예측된 결과는 Chaudhary et al.(2010)의 직접수치모사(direct numerical simulation, DNS) 자료와 비교하여 타당성을 검증하였으며, 평균속도 분포의 평탄화와 난류강도 및 전단응력의 억제 경향을 확인하였다. 해석 결과는 제안한 확장형 모형이 자기장 영향을 받는 유동장의 예측 정확도 향상에기여할 수 있음을 보여주고 복잡한 경계조건을 갖는 다양한 공학적 응용 분야에서 활용 가능성을 제시한다.

This study was conducted to understand the effects of a magnetic field on mean velocity and turbulent characteristics in fully developed turbulent magnetohydrodynamic (MHD) flow through a square duct. Based on the coupled Navier-Stokes MHD equations, an extended Reynolds stress transport model including the Lorentz force effect was applied, and the analysis was performed using an in-house computational fluid dynamics code. In order to directly apply the Lorenz force effect to the turbulent MHD flow, models related to the generation and dissipation rate of the Lorenz force were additionally applied to the Reynolds stress and dissipation rate equations, and the indirect effect was also considered through the modification of the velocity-pressure gradient correlation model of the Reynolds stress equation. In addition, the model coefficients of the models related to the Lorenz force were tuned using the fully developed channel MHD data so that they could be linked to the elliptic-blending Reynolds stress model and then the square duct flow was analyzed. The prediction results were validated through comparison with the direct numerical simulation (DNS) data of Chaudhary et al. (2010), confirming the flattening of the mean velocity distribution and the suppression trends of turbulence intensity and shear stress. The results of the analysis show that the proposed extended model contributes to improving the predictive accuracy of MHD-affected flow fields and suggests its applicability to various engineering applications involving complex boundary conditions.