27.28 도심 난류의 전산유체역학(CFD) 해석

27.28 도심 난류의 전산유체역학(CFD) 해석

1. 도심 난류 CFD 해석의 필요성

도심 환경의 난류는 건물 군집의 복잡한 기하, 다양한 스케일의 유동 현상, 그리고 열 효과의 결합으로 인해 분석적 접근으로는 정확히 기술하기 어렵다. 전산유체역학(Computational Fluid Dynamics, CFD)은 이러한 복잡 유동을 수치적으로 해석할 수 있는 강력한 도구이다. 본 절에서는 도심 난류의 CFD 해석을 위한 지배 방정식, 난류 모델, 격자 구성, 경계 조건, 그리고 검증 방법을 체계적으로 서술한다.

2. 지배 방정식과 가정

도심 난류의 CFD 해석은 비압축성 Navier-Stokes 방정식을 기반으로 한다.

\nabla \cdot \mathbf{u} = 0

\rho \frac{D\mathbf{u}}{Dt} = -\nabla p + \mu \nabla^2 \mathbf{u} + \rho \mathbf{g} + \mathbf{f}_{ext}

여기서 \mathbf{u}는 속도 벡터, p는 정압, \rho는 밀도, \mu는 동점성 계수, \mathbf{g}는 중력 가속도, \mathbf{f}_{ext}는 외부 체적력이다. 열 대류를 고려하는 경우 에너지 방정식과 부력 항이 추가된다. 난류는 레이놀즈 수 Re \sim 10^7 이상의 매우 큰 값을 가지므로 난류 모델링이 필수적이다.

3. 난류 모델링 접근

도심 난류의 CFD 해석에 사용되는 주요 난류 모델은 다음과 같다. 첫째, Reynolds Averaged Navier-Stokes(RANS) 모델은 시간 평균된 유동을 계산한다. 대표적 모델로는 표준 k-\epsilon, Realizable k-\epsilon, k-\omega SST가 있다. 계산 비용이 낮고 정상 상태 해석에 적합하다. 둘째, Large Eddy Simulation(LES)은 큰 규모 와류를 직접 해석하고 작은 규모는 서브그리드 모델로 처리한다. RANS보다 정확하나 계산 비용이 크다. 셋째, Detached Eddy Simulation(DES)은 벽면 부근은 RANS, 자유 영역은 LES로 처리하는 혼합 접근이다. 각 모델은 응용 목적과 계산 자원에 따라 선택된다.

4. 격자 구성의 핵심 요소

도심 난류 해석의 격자 구성은 다음을 고려한다. 첫째, 건물 형상을 정확히 표현하기 위해 건물 표면 주변에 세밀한 격자를 배치한다. 둘째, 벽면 함수(wall function) 또는 저-Reynolds 수 모델의 사용에 따라 벽면 부근 격자 해상도를 선택한다. 셋째, 관심 영역(예: UAM 비행 경로)에서 격자 밀도를 증가시킨다. 넷째, 전체 계산 영역의 크기는 건물 높이의 최소 5-10배 이상으로 설정하여 경계 영향을 최소화한다. 다섯째, 격자 수렴성 연구를 통해 결과의 격자 의존성을 평가한다.

5. 경계 조건의 설정

도심 CFD 해석의 주요 경계 조건은 다음과 같다. 첫째, 입구 경계에서는 대기 경계층 속도 프로파일(대수 법칙 또는 멱 법칙)을 부과한다. 둘째, 출구 경계에서는 압력 출구(pressure outlet) 또는 외삽 조건을 사용한다. 셋째, 건물 및 지면 표면에서는 점착 조건(no-slip condition)을 부과한다. 넷째, 측면 및 상부 경계에서는 대칭 조건 또는 자유 경계 조건을 사용한다. 난류 입구 조건(입구 난류 강도, 길이 스케일)의 설정도 중요하며, 종종 선행 해석의 결과를 사용한다.

6. RANS 기반 해석 방법

RANS 기반 해석은 도심 풍환경 평가의 가장 일반적 방법이다. 시간 평균된 속도장과 난류 통계를 제공하며, 수 시간 내지 수 일의 계산 시간을 요구한다. k-\omega SST 모델은 벽면 근처와 자유 유동 모두에서 양호한 성능을 보여 도심 해석에 자주 사용된다. 그러나 RANS는 비정상 유동 특성(예: 와류 흘림)을 정확히 포착하지 못한다. 이러한 한계는 비정상 RANS(URANS) 또는 LES로 보완된다.

7. LES 기반 해석

LES는 도심 난류의 시공간 구조를 상세히 해석한다. 서브그리드 스케일 모델로는 Smagorinsky 모델, dynamic 모델, WALE 모델 등이 사용된다. LES는 난류 통계뿐 아니라 순간 유동 구조도 제공하며, 오염 물질 확산이나 극한 돌풍 분석에 유용하다. 그러나 계산 자원 요구가 매우 크다. 전형적 도심 LES는 수억 격자와 수주간의 계산 시간을 요구한다.

8. 계산 자원과 병렬 계산

대규모 도심 CFD 해석은 고성능 컴퓨팅(HPC) 자원을 요구한다. 수천 개 내지 수만 개의 CPU 코어를 이용한 병렬 계산이 표준이다. 메모리 사용량도 상당하여 고성능 HPC 시스템이 필수적이다. 최근에는 GPU 가속도 도입되어 계산 효율이 향상되고 있다. 클라우드 컴퓨팅을 활용한 대규모 해석도 실용화되고 있다.

9. 대표 CFD 소프트웨어

도심 CFD 해석에 사용되는 주요 소프트웨어는 다음과 같다. 상용 소프트웨어로는 ANSYS Fluent, STAR-CCM+, OpenFOAM(오픈소스로 분류되기도 함)이 있다. 학술 연구용으로는 OpenFOAM, Code Saturne, OpenFAST 등이 있다. 풍력 공학 전용 도구로는 WAsP-CFD, Meteodyn WT 등이 있다. 각 소프트웨어는 강점과 약점이 있으며 응용 목적에 따라 선택된다.

10. 검증과 타당성 확인

CFD 해석 결과의 검증은 필수적이다. 대표적 검증 방법은 다음과 같다. 첫째, 풍동 시험 데이터와의 비교: 축척 도시 모형의 풍동 시험 결과가 기준 데이터로 사용된다. 둘째, 실측 데이터와의 비교: 기상 타워, 라이다, 드론 관측 등 실측 자료와의 비교. 셋째, 표준 벤치마크 문제와의 비교: Architectural Institute of Japan(AIJ) 가이드라인에 규정된 표준 도시 형상의 해석 결과를 기준과 비교. 이러한 검증은 해석 결과의 신뢰성을 확보하는 데 필수적이다.

11. 열 효과의 고려

도심 환경에서 열 효과는 유동 패턴에 상당한 영향을 미친다. CFD 해석에서 열 효과는 다음과 같이 반영된다. 첫째, 부력 항(Boussinesq 근사 또는 완전 부력 항)이 운동량 방정식에 포함된다. 둘째, 에너지 방정식이 풀이된다. 셋째, 건물 표면의 태양 복사 가열이 모델링된다. 넷째, 지표면 열 플럭스가 경계 조건으로 부과된다. 이러한 열 효과는 도시 열섬과 같은 현상의 재현에 필수적이다.

12. UAM 공역 평가 응용

UAM 기체의 안전한 도심 운용을 위한 공역 평가에 CFD가 활용된다. 특정 도심 영역의 3차원 풍환경 지도가 CFD로 계산되며, 평균 풍속, 난류 강도, 극한 돌풍 등의 통계가 제공된다. 이 정보는 UAM 공역 설계, 경로 계획, 그리고 기상 제한 조건 설정에 활용된다. 각국의 UAM 실증 사업에서 이러한 CFD 기반 평가가 수행되고 있다.

13. 실시간 응용을 위한 대리 모델

전통적 CFD는 계산 시간이 길어 실시간 응용에 제약이 있다. 이를 극복하기 위해 대리 모델(surrogate model)이 개발되고 있다. 주요 접근은 다음과 같다. 첫째, 다차원 룩업 테이블: 다양한 기상 조건에 대한 CFD 결과를 사전 계산하여 실시간 조회. 둘째, 감쇠 차수 모델(reduced-order model): 핵심 모드를 추출한 저차 모델. 셋째, 기계 학습 대리 모델: 신경망 등으로 CFD 결과를 학습. 이러한 접근은 실시간 UAM 공역 관리의 가능성을 제공한다.

14. 미래 발전 방향

도심 난류 CFD 해석은 다음 방향으로 발전하고 있다. 첫째, 도시 전체 규모의 초고해상도 LES. 둘째, 실시간 기상 데이터와의 결합을 통한 데이터 동화(data assimilation). 셋째, 기계 학습 기반 가속 기법. 넷째, 멀티 피지크스 결합(공기역학, 열전달, 구조, 전자기 등). 다섯째, 디지털 트윈(digital twin) 개념을 통한 도시 풍환경의 통합 관리. 이러한 발전은 안전하고 효율적인 도심 항공 운용의 기반을 제공할 것이다.

15. 출처

  • Blocken, B., “50 Years of Computational Wind Engineering: Past, Present and Future,” Journal of Wind Engineering and Industrial Aerodynamics, Vol. 129, 2014.
  • Tominaga, Y., Mochida, A., Yoshie, R., Kataoka, H., Nozu, T., Yoshikawa, M., and Shirasawa, T., “AIJ Guidelines for Practical Applications of CFD to Pedestrian Wind Environment Around Buildings,” Journal of Wind Engineering and Industrial Aerodynamics, Vol. 96, 2008.
  • Hertwig, D., Efthimiou, G. C., Bartzis, J. G., and Leitl, B., “CFD Validation of Wind Flow Around Isolated Buildings Using Wind Tunnel Measurements,” Journal of Wind Engineering and Industrial Aerodynamics, Vol. 111, 2012.
  • Ferziger, J. H., Perić, M., and Street, R. L., “Computational Methods for Fluid Dynamics,” 4th ed., Springer, 2020.
  • Spalart, P. R., “Detached-Eddy Simulation,” Annual Review of Fluid Mechanics, Vol. 41, 2009.
  • Britter, R. E., and Hanna, S. R., “Flow and Dispersion in Urban Areas,” Annual Review of Fluid Mechanics, Vol. 35, 2003.

16. 버전

  • 문서 버전: v1.0
  • 작성 기준일: 2026-04-17