유동계수

편집 역사 토론
유체역학 전산유체역학 마지막 수정: 2026-09-06 04:44:19

1. 개요[편집]

유동계수
Discharge / Flow Coefficient
기본 정의실제 유량 ÷ 이상 유량
분해Cd = Cc · Cv (수축 × 속도)
예리연 오리피스약 0.60 ~ 0.62
벤투리 관약 0.98 ~ 0.995
디젤 분사 노즐약 0.6 ~ 0.9 (입구 모서리 형상에 의존)
주의밸브의 Cv·Kv무차원이 아니다
측정정상류 유동벤치 / 교정 유량계

“이 밸브의 유동계수가 얼마죠?” 라는 질문에 답하기 전에, 그쪽이 말하는 유동계수가 무엇인지부터 물어야 한다. 최소 네 가지가 같은 이름으로 돌아다닌다.

유동계수(discharge coefficient, CdC_d)는 어떤 유로를 실제로 통과하는 유량과, 손실이 없다고 가정했을 때의 이상 유량의 비다.

Cd  =  m˙actualm˙idealC_d \;=\; \frac{\dot m_{\text{actual}}}{\dot m_{\text{ideal}}}

정의만 보면 더할 나위 없이 간단하지만, 이 계수가 악명 높은 이유는 분모가 규약이라는 데 있다. “이상 유량”은 자연이 정해 주는 값이 아니라 사람이 고른 기준 면적과 이상화 가정으로 계산한 값이며, 기준을 바꾸면 같은 장치의 CdC_d 가 두 배로 변한다. 그래서 실무 사고의 대부분은 물리를 몰라서가 아니라 서로 다른 규약을 섞어서 일어난다.

용어도 어지럽다. 한국어로 유동계수·유량계수·토출계수·유출계수가 뒤섞여 쓰이고, 영어로도 discharge coefficient 와 flow coefficient 가 다른 뜻으로 갈라진다. 이 문서는 각 분야의 정의를 명시적으로 구분하는 것을 주 목적으로 삼는다.

2. 이상 유량은 무엇인가[편집]

비압축성 유동에서 기준은 베르누이 방정식이다. 상류 정체압과 목 정압의 차 Δp\Delta p 만으로 속도를 정하면

m˙ideal=ρAref2Δpρm˙=CdAref2ρΔp\dot m_{\text{ideal}} = \rho\,A_{\text{ref}}\sqrt{\frac{2\,\Delta p}{\rho}} \qquad\Longrightarrow\qquad \dot m = C_d\,A_{\text{ref}}\sqrt{2\rho\,\Delta p}

이다. 접근속도가 무시할 수 없는 관로 계측에서는 여기에 1/1β41/\sqrt{1-\beta^4} 가 붙고, 이 인자를 CdC_d 에 흡수한 것을 따로 유량계수 KK 라 부른다. 두 계수를 같은 기호로 쓰는 문서가 실제로 섞여 있고 β=0.7\beta = 0.7 이면 둘이 15 % 가까이 차이 나므로, 계산서를 받으면 첫 줄에서 정의식을 확인하는 것이 습관이 되어야 한다.

기체에서 압력비가 커지면 베르누이가 아니라 등엔트로피 노즐 식이 기준이 된다.

m˙=CdArefp0RT0(pp0)1/γ{2γγ1[1(pp0)(γ1)/γ]}1/2\dot m = C_d\,A_{\text{ref}}\,\frac{p_0}{\sqrt{R T_0}} \left(\frac{p}{p_0}\right)^{1/\gamma} \left\{\frac{2\gamma}{\gamma-1}\left[1-\left(\frac{p}{p_0}\right)^{(\gamma-1)/\gamma}\right]\right\}^{1/2}

압력비가 임계값 아래로 내려가면 이상 유량 쪽이 초크되어 상류 조건만의 함수가 되고(압축성 유동·라발 노즐), 그 뒤로도 CdC_d 는 대체로 완만하게 변한다. 밸브·포트처럼 압력비가 크게 흔들리는 자리에서 비압축 식으로 뽑은 CdC_d 를 압축성 조건에 그대로 쓰면 안 되는 이유가 여기 있다. 규약이 다르기 때문이지 물리가 달라서가 아니다.

3. 왜 1보다 작은가 — 수축계수와 속도계수[편집]

CdC_d 가 1 이 아닌 데는 서로 다른 두 원인이 있고, 이를 분리한 것이 고전적인 분해다.

Cd=CcCv,Cc=AvenaAref,Cv=uactualuidealC_d = C_c\,C_v, \qquad C_c = \frac{A_{\text{vena}}}{A_{\text{ref}}},\quad C_v = \frac{u_{\text{actual}}}{u_{\text{ideal}}}
  • 수축계수 CcC_c — 유선이 관성 때문에 계속 수렴해, 최소 단면(축류, vena contracta)이 기하학적 구멍보다 작아지는 효과. 점성과 무관한 순수 기하·관성 효과다.
  • 속도계수 CvC_v — 그 최소 단면에서의 실제 속도가 무손실 값보다 작은 효과. 이쪽이 점성 손실이다.

예리연 오리피스에서 두 값의 크기 차이가 극적이다. CvC_v 는 0.97~0.99 로 1 에 매우 가깝고, CcC_c 가 0.61 근처다. 즉 오리피스의 Cd0.6C_d \approx 0.6 은 마찰이 아니라 거의 전부 축류가 만든 숫자다. 자유유선 이론이 2차원 예리연 슬릿에 대해 Cc=π/(π+2)0.611C_c = \pi/(\pi+2) \approx 0.611 을 주는 것이 그 근거이며, 자세한 것은 오리피스 쪽에 있다. 반대로 벤투리 관은 수축을 완만한 원추로 처리해 축류가 생기지 않으므로 Cc1C_c \approx 1, 결국 CdC_d 가 0.99 근처로 올라간다. 두 장치의 계수 차이 전부가 CcC_c 에서 온다.

여기서 첫 번째 기호 충돌이 발생한다. CvC_v 는 속도계수이기도 하고 밸브의 유량계수이기도 하다. 둘은 무차원이냐 아니냐부터 다르다.

4. 관로 계측 — ISO 5167 계열[편집]

차압식 유량계의 CdC_d 는 세계적으로 규격화되어 있다. ISO 5167(오리피스·노즐·벤투리)이 형상·설치·탭 위치를 못 박고, 그 형상에 대해 CdC_d 를 경험식으로 준다. 오리피스는 예전 ISO 5167:1991 의 스톨츠(Stolz) 식이 2003년 개정에서 리더-해리스/갤러거 식으로 교체됐고, 지금 계산서에 찍히는 값은 후자다.

읽어야 할 것은 계수가 아니라 구조다.

  • CdC_d레이놀즈수의 함수다. ReD\mathrm{Re}_D 가 크면 상수에 수렴하지만 저유량에서는 유량의 함수가 되어 유량식이 음함수가 되고, 계산이 반복법으로 돈다.
  • CdC_d탭 위치에 의존한다. 같은 판이라도 코너 탭·플랜지 탭·D–D/2 탭이 다른 값을 쓴다. Δp\Delta p 를 어디서 재느냐가 정의의 일부이기 때문이다.
  • 규격이 보장하는 범위 밖(작은 관, 낮은 Re\mathrm{Re}, 큰 β\beta)에서는 규격값이 아니라 교정 데이터가 필요한 장치가 된다.

세부는 오리피스벤투리 관 문서에 있으므로 여기서는 반복하지 않는다.

5. 밸브의 C_v · K_v — 무차원이 아니다[편집]

프로세스 배관에서 “flow coefficient”라 하면 거의 언제나 밸브의 CvC_v 를 뜻하고, 이것은 위에서 본 무차원 CdC_d완전히 다른 종류의 양이다. 정의가 아예 실험 절차로 되어 있다.

  • CvC_v — 60 °F 물이 밸브를 지날 때 차압 1 psi 를 만드는 유량을 US gpm 으로 적은 수.
  • KvK_v — 5~40 °C 물이 지날 때 차압 1 bar 를 만드는 유량을 m³/h 로 적은 수.

실용식은 비중 GG 에 대해 Q=CvΔp/GQ = C_v\sqrt{\Delta p / G} 꼴이며, 두 정의는 단위 환산으로 이어진다. Cv=1C_v = 1 을 넣어 보면

Kv=0.2271 m3/h0.06895 bar0.865CvK_v = \frac{0.2271\ \mathrm{m^3/h}}{\sqrt{0.06895\ \mathrm{bar}}} \approx 0.865\,C_v

가 나온다. 여기서 반드시 짚어야 할 점 셋.

  1. 차원이 있다. CvC_v 는 순수한 수가 아니라 단위계에 묶인 양이다. 미터법 문서와 인치법 문서를 섞으면 15 % 가 조용히 사라진다.
  2. 면적 정보가 이미 들어 있다. SI 로 정리한 Av=Q/Δp/ρA_v = Q/\sqrt{\Delta p/\rho} 는 실제로 면적 차원(m2\mathrm{m^2})이고, KvK_v 를 m³/h 로 준 값이라면 Av=Kv/36000A_v = K_v/36000 이다. 즉 CvC_v·KvK_v유효 유동면적을 사람이 읽기 좋은 숫자로 포장한 것이며, 무차원 CdC_d 로 되돌리려면 기준 면적이 따로 필요하다.
  3. 밸브 개도의 함수다. 카탈로그의 CvC_v 는 전개(全開) 값이고, 실제 제어에는 개도-CvC_v 곡선(선형·등비율·퀵오프닝)이 쓰인다. 제어 밸브의 특성 선택이 사실상 이 곡선 선택이다.

또한 압축성 유체와 캐비테이션·초킹 조건에서는 이 식이 그대로 성립하지 않아, IEC 60534 계열이 팽창계수와 초크 판정 인자를 따로 정의한다. 액체가 초크되면(축류 압력이 증기압에 고정) 차압을 더 키워도 유량이 늘지 않으므로 Δp\sqrt{\Delta p} 관계 자체가 깨진다(캐비테이션).

6. 엔진 흡배기 포트 — 기준 면적이 두 개다[편집]

내연기관 쪽은 또 다른 규약을 쓴다. 밸브를 준정상 오리피스로 보고

Aeff=CdArefA_{\text{eff}} = C_d\,A_{\text{ref}}

유효 유동면적을 정의하는 것까지는 같은데, ArefA_{\text{ref}} 로 무엇을 쓰느냐가 두 갈래다.

기준 면적정의리프트 의존관용 명칭
커튼 면적πDvL\pi D_v L있음discharge coefficient CdC_d
밸브 헤드(시트 안쪽) 면적πDv2/4\pi D_v^2/4없음flow coefficient CfC_f

같은 측정에서 나오는 두 숫자가 인자 4L/Dv4L/D_v 만큼 다르다. L/Dv=0.25L/D_v = 0.25 면 네 배 차이다. 커튼 면적 기준 CdC_d 는 작은 리프트에서 0.6~0.7 로 시작해 리프트가 커지면 유로의 목이 시트에서 포트로 옮겨 가며 떨어지고, 헤드 면적 기준 CfC_f 는 리프트에 따라 단조 증가하다 포화한다. 곡선 모양이 정반대라 그래프만 보고도 어느 규약인지 알 수 있을 정도인데, 논문이 규약을 안 밝히면 그 그래프는 재사용이 불가능하다.

측정은 정상류 유동벤치(steady-flow bench)로 한다. 실린더 헤드를 어댑터에 물리고 밸브 리프트를 고정한 뒤, 정해진 차압(업계 관행으로 28 inH₂O ≈ 7 kPa 가 흔하다)에서 유량을 재 리프트별 곡선을 만든다. 여기서 나온 Aeff(L)A_{\text{eff}}(L) 을 크랭크각의 함수로 바꾼 것이 1D 엔진 코드의 밸브 입력이며, 체적효율가변 밸브 타이밍 해석의 실질적 기반이다. 같은 벤치에서 임펠러나 토크미터로 스월·텀블 비도 함께 잰다.

이 접근에는 명백한 근사가 하나 깔려 있다. 정상류로 잰 계수를 비정상 유동에 쓴다는 것이다. 밸브가 실제로는 움직이고 있고 상류에는 압력파가 오가지만, 밸브 통과 시간척도가 파동 시간척도보다 짧다는 준정상 가정으로 정당화한다. 대체로 잘 통하지만 밸브 개폐 직후의 과도 구간과 오버랩 구간에서는 오차가 커진다.

디젤 분사 노즐도 같은 틀이다. 커먼레일 인젝터의 분사 구멍 CdC_d 는 입구 모서리 형상에 강하게 의존해서, 예리한 구멍은 축류와 캐비테이션 때문에 0.6~0.7 대에 머물고, 유체 연마로 입구를 둥글게 만들면 0.85 이상으로 오른다. 인젝터 내부의 제어 오리피스들도 마찬가지로 CdC_d 가 곧 설계 파라미터다.

7. 같은 이름의 다른 무차원수들[편집]

  • 터보기계의 유량계수 ϕ\phi. 축류 기계에서 ϕ=cx/U\phi = c_x/U(축방향 속도 ÷ 익단 속도), 원심 펌프·송풍기에서는 ϕ=Q/(ND3)\phi = Q/(N D^3) 로 정의한다. 부하계수 ψ\psi 와 짝을 이루어 성능 특성곡선의 좌표축 노릇을 하는 상사 파라미터이며, “실제/이상”과는 아무 관계가 없다(터보기계 해석·차원 해석).
  • 웨어·수문의 유량계수. 개수로 유동에서 예연 웨어는 Q=Cd232gbH3/2Q = C_d \frac{2}{3}\sqrt{2g}\,b\,H^{3/2}, 수문(sluice gate)은 게이트 아래 축류를 Cc0.61C_c \approx 0.61 로 놓는 식이 표준이다. 구조는 관로 계측과 같지만 기준 수두를 어디서 재는지가 또 다른 규약이다.
  • 분출계수·유출계수. 안전밸브(PRV)나 노즐 쪽에서 쓰는 이름들로, 대개 같은 CdC_d 를 가리키지만 규격마다 인증된 값이 따로 지정된다.

정리하면 이렇다. “실제/이상”형 계수는 무차원이고 규약 의존적이며, “1 psi 당 유량”형 계수는 차원이 있고 단위 의존적이고, "cx/Uc_x/U"형 계수는 상사 파라미터다. 셋은 같은 문장 안에 들어올 수 있는 종류의 양이 아니다.1

8. CFD로 CdC_d 를 뽑을 때[편집]

CFD 는 CdC_d 를 얻는 유력한 수단이지만, 이 양이 두 수의 비라는 성질 때문에 특유의 함정이 있다.

  • 분모를 실험과 똑같이 계산했는가. 가장 흔하고 가장 조용한 오류다. 실험은 특정 탭 위치의 정압차로 Δp\Delta p 를 정의하는데, CFD 는 전 영역의 압력을 알고 있으므로 아무 데서나 뽑을 수 있다. 입구 정체압에서 출구 정압을 빼면 계측 규약과 다른 Δp\Delta p 가 되고, 결과는 몇 %가 아니라 수십 %가 어긋난다. 탭 위치를 격자에서 재현하고 거기서 샘플링하는 것이 원칙이다. 기준 면적도 마찬가지다.
  • 격자 수렴은 CdC_d 로 본다. 예리연 모서리는 곡률 반경이 0 이라 그 점에서는 어떤 격자로도 수렴하지 않는다. 국소 극한값을 쫓지 말고 적분량인 CdC_d격자 수렴을 확인하는 것이 옳다. 다만 CdC_d 는 유량이라는 적분량이라 격자에 둔감해 보이기 쉬운데, 정작 그것을 결정하는 축류 위치와 박리 전단층은 매우 민감하다. CdC_d 가 안 변하니 수렴했다”가 아니라 축류 위치와 재부착 길이가 안 변하는지를 함께 봐야 한다.
  • 박리 예측이 곧 계수 예측이다. CdC_d 를 정하는 것은 유동박리로 만들어진 축류와 재순환이다. 표준 k-ε 은 재순환 영역을 짧게 예측해 CdC_d 를 낙관적으로 내는 경향이 있고, k-ω SST 나 레이놀즈 응력 모델이 낫다. 저 Re\mathrm{Re} 구간에서는 난류 천이 처리가 결과를 갈라 놓는다(난류 모델링).
  • 정상 해석의 한계. 박리 전단층은 본질적으로 비정상이라 정상 RANS 는 존재하지 않는 정상해를 억지로 만든다. 비정상 해석의 시간평균 CdC_d 와 정상해의 CdC_d 가 몇 % 다른 것은 정상이며, 어느 쪽을 실험과 비교할지 미리 정해야 한다.
  • 캐비테이션과 초킹은 별개 물리다. 액체에서 축류 압력이 증기압에 닿거나 기체에서 유동이 음속에 닿으면 Δp\sqrt{\Delta p} 관계 자체가 무너진다. 단상 비압축 해석으로 뽑은 CdC_d 를 그 영역에 외삽하면 안 된다(다상유동).
  • 검증 케이스로서 훌륭하다. ISO 규격 형상은 CdC_d 가 불확도까지 표로 존재하는 몇 안 되는 유동이라, 새 솔버·격자 전략의 검증 및 확인 문제로 자주 쓰인다.

9. 시스템 해석에서의 위치[편집]

CdC_d 가 결국 하는 일은 3차원 국소 유동을 계수 하나로 축약해 1차원 망 해석에 넘기는 것이다. 배관망·유압회로·엔진 흡배기 1D 모델은 전부 이 구조 위에 서 있다 — 부품 하나를 제곱 법칙 저항 소자로 놓고, 그 계수만 CFD 나 시험으로 뽑아 넣는다.

Δp    m˙2\Delta p \;\propto\; \dot m^2

전기 저항과 달리 제곱 법칙이라는 점이 수치적으로 성가시다. m˙0\dot m \to 0 에서 Δp/m˙0\partial \Delta p / \partial \dot m \to 0 이 되어 뉴턴 반복의 자코비안이 특이해지므로, 원점 부근을 층류식으로 매끄럽게 갈아 끼우는 처리가 사실상 필수다. 밸브가 닫히는 순간 시스템 해석이 통째로 발산하는 사고의 표준 원인이 이것이다. 같은 논리가 압력손실 부품 전반과 밸브·유량계 모델에 그대로 적용된다.23

10. 관련 문서[편집]

11. Footnotes[편집]

  1. 그럼에도 회의 자료에서는 태연히 한 슬라이드에 같이 등장한다. 발표자가 “유동계수가 0.62 에서 320 으로 개선됐습니다”라고 말하는 순간 뒤에서 누군가 커피를 뿜는 것이 이 바닥의 국룰이다. 앞은 무차원 CdC_d, 뒤는 밸브 CvC_v 다.

  2. 처음 유압 회로 코드를 짜면 반드시 한 번은 겪는다. 밸브 개도를 0 으로 보내는 순간 뉴턴이 발산하고, 로그를 열어 보면 자코비안 조건수가 101610^{16} 을 찍고 있다. 물리가 틀린 게 아니라 정의역 끝에서 함수가 미분 불가능해진 것이며, 처방은 언제나 원점 근방 정칙화다. “수렴은 신에게 맡긴다”가 통하지 않는 몇 안 되는 자리.

  3. CdC_d 의 진짜 가치는 정확도가 아니라 이식성이다. 3차원 해석 한 번으로 뽑은 무차원 수 하나가 그 부품이 들어가는 모든 시스템 모델에서 재사용된다. 반대로 말하면 규약을 안 적어 놓은 CdC_d 는 이식성이 0 이라 그냥 버리는 편이 빠르다. 계수를 기록할 때 기준 면적·차압 정의·레이놀즈수 범위를 함께 적는 것이 사실상 그 계수의 사용설명서다.