본문 바로가기

CFD 해석 및 검증/한국원자력학회 2026 춘계학술발표회 · 2026

싼 계산은 어디까지 비싼 계산을 대신할 수 있나

관 다발을 가로지르는 흐름은 관 뒤에서 소용돌이를 떨어뜨리고, 그 압력 변동이 관을 흔듭니다. 흔들림을 예측하려면 LES 가 필요하지만 비쌉니다. 그보다 훨씬 싼 DES 가 어디까지 대신할 수 있는지를 표면 압력으로 따졌습니다.

  • 조범진
  • 송민섭

왜 이걸 보나

소형모듈원자로(SMR)는 좁은 공간에서 열을 많이 옮겨야 해서 나선형 코일 증기발생기(HCSG)를 자주 씁니다. 관들이 주 유동에 직각으로 놓여 서로 엇갈리므로, 냉각재가 관 다발을 가로질러 지나갑니다.

이 교차유동은 관 뒤에서 박리와 후류를 만들고, 관 표면의 압력을 끊임없이 흔듭니다. 이 흔들림이 유동유발진동(FIV)입니다. 설계 단계에서 관에 걸리는 비정상 하중을 정량적으로 따져 둬야 오래 운전할 수 있습니다.

그런데 시간 평균을 전제로 하는 RANS 는 와류 탈락과 벽면 압력 변동을 지나치게 눌러 버립니다. LES 는 그걸 살려내지만 계산이 비쌉니다. DES 는 벽 근처만 RANS 로 풀고 박리·후류는 LES 처럼 푸는 절충안입니다. 이 절충이 압력 기준으로도 통하는지가 이 연구의 질문입니다.

관 하나를 둘러싸고

각도에 따라 전혀 다른 일이 일어납니다

가운데 관을 기준으로 각도 θ 를 재면, 자리마다 흐름의 성격이 다릅니다.

2차원 비압축성 Navier–Stokes 를 직접 풀어 얻은 속도장입니다. 논문 조건이 아닙니다. 같은 것은 피치비 P/D = 1.32 뿐입니다. 이 장은 층류(관 지름 기준 Re 100, 간극 속도 기준 약 410)이고, 논문은 Re = 7.82 × 10⁴ 입니다.

손잡이를 끌어 각도를 바꿔 보세요.

θ ≈ 0° · 정체점

흐름이 관 앞면에 정면으로 부딪쳐 멈추는 자리입니다. 압력이 가장 높고, 세 모형이 모두 실험과 잘 맞는 구간이기도 합니다.

어떤 조건

풍동 실험 하나를 기준으로 삼았습니다

공개된 풍동 벤치마크에서 대표 조건 하나를 골라, 세 가지 난류 모형을 같은 조건에서 돌리고 가운데 관의 둘레 방향 평균 압력계수 Cp̄(θ) 를 실험과 직접 맞대었습니다.

피치비 P/D
1.32
입구 속도 U
7 m/s
간극 속도 Ug
29.2 m/s
레이놀즈 수
7.82 × 10⁴
관 지름 D
38 mm
통계 구간
0.04 – 0.24 s

비교한 세 가지

  • URANS (SST k–ω)평균 유동과 평균 압력장을 싸게 얻는 기준선
  • LES (WALE)벽 근처 난류와 비정상 후류를 풀어내는 고정밀 기준
  • DES (SST-IDDES)벽 근처는 RANS, 박리·후류는 LES 처럼 — 절충안

격자

IDDES 가 매끄럽게 넘어가도록

DES 격자는 벽 근처는 URANS-Medium 수준으로, 후류 구간은 LES 급 해상도로 짰습니다. 두 모드가 부드럽게 전환되도록 하기 위해서입니다.

엇갈린 관군의 DES 격자. 입구·집중·이탈·바깥 구간으로 나뉘어 있고, 관 사이 점성 구간을 확대해 보여 준다.
논문 Fig. 6

결과 하나

정체점에서는 맞고, 후류에서 갈립니다

앞쪽 정체점(θ ≈ 0°) 근처에서는 세 모형이 모두 실험과 잘 맞습니다. 갈라지는 곳은 박리와 후류입니다. θ ≈ 180° 의 압력 회복을 LES 는 따라가는 반면, DES 는 낮게 보고 URANS 와 비슷한 모양이 됩니다.

가운데 관의 둘레 방향 압력계수 Cp(θ). 실험, LES, DES, URANS(Fine) 네 곡선이 겹쳐 그려져 있다.
논문 Fig. 7

실험과의 차이 (논문 Table Ⅳ)

LESDESURANS
RMSE0.0700.17860.1645
MAE0.05500.16480.1392
Mean Error-0.0404-0.1648-0.1205
Max |Δ|0.15730.30480.3000
Angle at Max |Δ|240°240°320°

결과 둘

평균이 같아도 떨림은 다릅니다

벽면 압력의 시간 이력을 스펙트럼으로 보면, 평균 Cp 가 비슷한 자리에서도 주 진동수가 다릅니다. θ ≈ 90° 에서 DES 의 주 피크는 LES 의 절반 아래로 내려가, 저주파 성분이 더 크다는 것을 보여 줍니다.

위치LESDES
θ = 90°117.19 Hz (St 0.153)58.59 Hz (St 0.076)
θ = 180°195.31 Hz (St 0.254)205.08 Hz (St 0.267)

그래서

속도장만 보면 놓치는 것이 있습니다

DES 검증은 보통 속도장을 맞대어 봅니다. 이 연구는 표면 압력과 FIV 관점의 비정상 지표를 기준으로 삼았습니다. 설계 하중과 직접 이어지는 값이기 때문입니다.

결론은 절반의 합격입니다. 평균 압력의 큰 그림은 따라가지만, 후류 쪽 압력 회복을 체계적으로 낮게 봅니다. 지금 조건에서 DES 를 LES 의 대체재로 쓰려면 이 편차를 알고 써야 합니다.

다음에 해야 할 것 (논문에 적힌 대로)

  • 01통계 구간을 더 길게 잡아 저주파 성분과 평균값의 수렴을 확인하기
  • 02더 느린 속도·낮은 레이놀즈 수 조건(U = 2–3 m/s)에서 같은 지표를 다시 재기
  • 03후류 구간 해상도와 RANS–LES 전환 거동에 대한 민감도 확인
  • 04대표 조건에 대해 입자영상유속계(PIV) 로 후류와 간극 유동을 검증하기

출처

조범진 · 송민섭, “Pressure-Based Assessment of SST-IDDES for Surface Pressure and FIV-Relevant Metrics in Staggered Tube-Bundle Cross-Flow”, 한국원자력학회 2026 춘계학술발표회, 2026.

CFD 해석 및 검증 분야 보기