관 유동 교차검산 재확인: 2026-10-22
계산으로 답할 질문을 먼저 정하기
지름 10 mm인 원형 관에 평균 속도 0.2 m/s로 유체가 흐른다고 가정해 봅니다. 여기서 알고 싶은 것은 관성과 점성의 상대적인 크기, 그리고 열이 흐름에 실려 가는 속도와 전도로 퍼지는 속도의 비율입니다. 두 질문은 서로 다르므로 레이놀즈 수 하나만 계산해서 열전달까지 결론을 내리면 안 됩니다. 이 노트에서는 같은 입력으로 세 계산기를 연결해 결과가 서로 맞는지 확인하고, 이어서 실제 물성표의 물 값으로 두 온도에서 계산을 반복합니다.
예제 1은 손으로 따라 계산하기 쉽도록 반올림한 상수(ρ = 1000 kg/m³, μ = 0.001 Pa·s, cₚ = 4200 J/(kg·K), k = 0.6 W/(m·K))를 씁니다. 이 값은 특정 온도의 물을 정확히 재현하는 자료가 아닙니다. 예제 2와 3은 IAEA 포화수 물성표(표 3.3)의 20 °C, 60 °C 행에서 ρ, μ, cₚ, k를 가져오며, 포화액 값을 저압 액체 물의 근삿값으로 사용합니다. 관 지름과 유속은 세 예제 모두 가정한 운전 조건입니다. 길이 L은 처음부터 끝까지 관 안지름으로 고정합니다.
1. 입력표를 SI 단위로 정리하기
| 물리량 | 입력 | 선택 이유 |
|---|---|---|
| 밀도 ρ | 1000 kg/m³ | 점성과 열확산의 환산에 같은 밀도 사용 |
| 평균 속도 v | 0.2 m/s | 관 단면의 평균값 |
| 안지름 L | 0.01 m | 10 mm ÷ 1000 |
| 점성 계수 μ | 0.001 Pa·s | 운동 점성 계수와 구분 |
| 정압 비열 cₚ | 4200 J/(kg·K) | 열확산 계산에 사용 |
| 열전도율 k | 0.6 W/(m·K) | 같은 유체의 열전도율 |
지름을 10으로 입력하면 계산기는 10 m로 해석합니다. Re와 Pe가 모두 1000배 커져도 두 결과 사이의 관계는 그대로 유지됩니다. 따라서 교차검산이 맞는다고 해서 입력 단위까지 맞았다고 볼 수는 없으며, 단위는 별도로 확인해야 합니다. 비열도 마찬가지입니다. 4200 J 대신 4.2 kJ의 숫자만 넣으면 Re는 그대로인데 Pr과 Pe가 1000분의 1로 작아집니다.
같은 10 mm 관으로 계산한 세 가지 사례
예제 1. 반올림 상수로 Re, Pr, Pe를 두 경로로 구하기
레이놀즈 수
운동 점성 계수로도 확인할 수 있습니다. ν = μ/ρ = 10⁻⁶ m²/s이므로 Re = vL/ν = 0.002/10⁻⁶ = 2000입니다. Pa·s 단위의 μ를 넣는 칸에 m²/s 단위의 ν를 넣으면 두 계산은 서로 다른 뜻이 됩니다.
2000은 매끈한 원형 관의 완전히 발달한 유동을 생각할 때 층류 범위에 드는 안내값입니다. 입구 길이, 교란, 거칠기는 별도로 검토해야 합니다. 계산기의 구간 표시는 그 조건들을 측정하거나 판정한 결과가 아닙니다.
Re 계산기에 입력하기프란틀 수와 열확산 계수
Pr = μcₚ/k = 0.001 × 4200 / 0.6 = 7
α = k/(ρcₚ) = 0.6/(1000 × 4200) = 1.428571… × 10⁻⁷ m²/s
독립 확인: ν/α = 10⁻⁶ / (1.428571… × 10⁻⁷) = 7
Pr = 7은 이 입력에서 운동량 확산율이 열확산율의 7배라는 뜻입니다. 물성값을 고정한 채 속도만 두 배로 바꾸면 Pr은 7로 유지됩니다. 속도만 바꿨는데 Pr이 달라졌다면 다른 물성이나 단위가 함께 바뀌었는지 확인하세요.
Pr 계산기에 입력하기페클레 수를 두 경로로 확인하기
다른 경로는 Re × Pr = 2000 × 7 = 14000입니다. 곱을 전개하면 (ρvL/μ)(μcₚ/k) = ρcₚvL/k = vL/α가 됩니다. 길이·속도·물성 조건을 일치시키면 이 등식은 반올림 오차를 제외하고 성립합니다. α를 먼저 1.43 × 10⁻⁷로 반올림하면 직접 계산한 Pe는 약 13986이 됩니다. 표시값이 조금 다르다면 계산기 오류를 의심하기 전에 입력 정밀도부터 확인하세요.
시간 척도로도 읽을 수 있습니다. t이송 = L/v = 0.05 s, t확산 = L²/α = 700 s이고 그 비가 14000입니다. 이 비가 크다는 것은 유체가 지름만큼 흘러가는 동안 열은 지름 방향으로 거의 퍼지지 못한다는 뜻입니다. 두 시간은 선택한 L에서 만든 척도일 뿐이며, 실제 출구 온도나 관 전체를 가열하는 데 걸리는 시간을 예측한 값은 아닙니다.
Pe 계산기에 입력하기예제 2. 같은 관에 IAEA 물성표의 20 °C 물을 넣기
주어진 값. IAEA 표 3.3(포화수, 20 °C)에서 ρ = 998.2 kg/m³, μ = 1.0000 × 10⁻³ Pa·s, cₚ = 4.185 kJ/(kg·K) = 4185 J/(kg·K), k = 0.5995 W/(m·K)를 가져옵니다. 같은 행에는 ν = 1.0020 × 10⁻⁶ m²/s, α = 1.435 × 10⁻⁷ m²/s, Pr = 6.981도 실려 있습니다. 표의 열 머리글에는 배율이 붙어 있으므로 읽을 때 주의해야 합니다. 점성 열의 10.000은 10⁴μ 값이어서 μ = 1.0000 × 10⁻³ Pa·s가 되고, 열전도율 열의 599.5는 10³k 값이어서 k = 0.5995 W/(m·K)가 됩니다. 비열은 kJ/(kg·K)로 실려 있으므로 1000을 곱해 J 단위로 바꿉니다. L = 0.01 m와 v = 0.2 m/s는 가정한 조건입니다.
Re = 998.2 × 0.2 × 0.01 / 0.001 = 1996.4
Pr = 0.001 × 4185 / 0.5995 = 6.9808 (표: 6.981)
α = 0.5995/(998.2 × 4185) = 1.4351 × 10⁻⁷ m²/s (표: 1.435 × 10⁻⁷)
Pe = vL/α = 0.002/(1.4351 × 10⁻⁷) ≈ 13936.5, Re × Pr = 1996.4 × 6.9808 ≈ 13936.5
t확산 = L²/α ≈ 696.8 s
해석. 예제 1의 반올림 상수는 표의 20 °C 값과 약 0.4% 이내로 차이가 나므로 Re, Pr, Pe도 거의 변하지 않습니다. 눈여겨볼 점이 하나 있습니다. 표에 실린 ν로 계산하면 Re = 0.002/(1.0020 × 10⁻⁶) = 1996.0이 되어 1996.4와 다릅니다. μ/ρ를 직접 계산하면 1.0018 × 10⁻⁶인데 표는 ν를 따로 반올림해서 싣기 때문입니다. 두 경로의 결과가 다섯째 유효숫자에서 어긋난다면 공식이 틀린 것이 아니라 물성표의 반올림 흔적으로 보는 것이 맞습니다. 등식을 깔끔하게 검산하려면 ν, α, Pr을 모두 한 벌의 ρ, μ, cₚ, k에서 계산하세요.
예제 3. 같은 속도에서 물을 60 °C로 데우기
주어진 값. IAEA 표 3.3의 60 °C 행에서 ρ = 983.2 kg/m³, μ = 4.677 × 10⁻⁴ Pa·s, cₚ = 4183 J/(kg·K), k = 0.6510 W/(m·K)를 가져옵니다. 표에 실린 Pr은 3.005, α는 1.583 × 10⁻⁷ m²/s입니다. L = 0.01 m, v = 0.2 m/s는 예제 2와 같게 가정하고 모든 물성을 60 °C에서 평가합니다. 단위 환산은 예제 2와 같으며, 점성 열의 4.677은 10⁴μ 값이므로 4.677 × 10⁻⁴ Pa·s로 읽습니다.
Re = 983.2 × 0.2 × 0.01 / (4.677 × 10⁻⁴) = 4204.4
Pr = 4.677 × 10⁻⁴ × 4183 / 0.6510 = 3.0052
α = 0.6510/(983.2 × 4183) = 1.5829 × 10⁻⁷ m²/s
Pe = 0.002/(1.5829 × 10⁻⁷) ≈ 12635 = Re × Pr
t확산 = L²/α ≈ 631.8 s
해석. 표를 보면 20 °C에서 60 °C로 올라가면 점성 계수가 절반 아래로 떨어지므로 같은 속도에서도 Re가 두 배 넘게 커집니다. 반면 Pe는 오히려 약 9% 줄어듭니다. 열전도율이 0.5995에서 0.6510 W/(m·K)로 약 9% 커지고 ρcₚ는 조금 줄어서, 열확산 계수가 약 10% 늘었기 때문입니다. Re와 Pe는 서로 다른 질문에 답하는 수이며, 이 예제에서는 변하는 방향까지 반대입니다. 또한 Re = 4204는 원형 관 천이의 대표값인 약 2300을 넘습니다. 따라서 Nu = 3.66이나 f = 64/Re 같은 층류 결과를 실제 유동 상태를 확인하지 않고 그대로 가져다 쓰면 안 됩니다. 가열되는 관에서는 길이를 따라 유체 온도가 변하므로, 60 °C 사례는 특정 가열기를 묘사한 결과가 아니라 범위를 가늠하는 경계 계산으로 읽어야 합니다.
세 사례의 결과 비교
| 사례 | Re | Pr | α (m²/s) | Pe |
|---|---|---|---|---|
| 예제 1, 반올림 상수 | 2000 | 7 | 1.4286 × 10⁻⁷ | 14000 |
| 예제 2, IAEA 물 20 °C | 1996.4 | 6.9808 | 1.4351 × 10⁻⁷ | 13936.5 |
| 예제 3, IAEA 물 60 °C | 4204.4 | 3.0052 | 1.5829 × 10⁻⁷ | 12635 |
예제 1에서 v만 0.4 m/s로 바꾸면 Re = 4000, Pr = 7, Pe = 28000입니다. 표에서 α 열은 세 사례 사이의 차이가 10% 안팎에 그치지만 Re 열은 두 배 이상 벌어집니다. 어떤 수가 온도에 민감한지 비교하려면 이처럼 같은 형식의 표로 정리해 두는 편이 좋습니다. v만 바꾼 변화는 물성을 고정했을 때 공식이 보이는 민감도일 뿐이며, 물성이 변하는 실제 가열 유동의 해답을 대신하지 않습니다. 온도가 물성을 바꾸고 그 결과 Re, Pr, Pe가 함께 움직이는 더 현실적인 경우는 예제 3이 보여 줍니다.
층류이고 완전히 발달한 유동인지 확인하기
2300은 단면이 일정한 원형 관에서 쓰는 천이의 대표값이며 날카로운 경계가 아닙니다. MIT 생리유체역학 강의 자료는 이런 관의 유동이 대개 Re 약 2300을 넘으면 난류가 된다고 설명하고, 좁은 오리피스에서 넓은 관으로 나오는 유동은 Re가 10 정도로 낮아도 난류가 될 수 있다고 적고 있습니다. 경계를 다르게 잡는 자료도 있습니다. TU Delft의 슬러리 수송 교재는 2320을 기준으로 64/Re와 Colebrook–White 마찰계수를 나누고, WIPP 성능 평가 모델은 Re 2100 미만을 층류로 취급합니다. 예제 2의 Re ≈ 1996은 세 기준보다 모두 낮고, 예제 3의 Re ≈ 4204는 모두 높습니다.
완전 발달 여부는 따로 따져야 하는 조건입니다. 같은 MIT 자료는 층류 입구 길이를 Le/D = 0.065 Re로 제시합니다. 예제 1에 적용하면 0.065 × 2000 = 130으로 지름의 130배이고, 10 mm 관에서는 약 1.3 m입니다. 길이 0.5 m인 시험 구간이라면 전체가 발달 구간에 들어가며, 그곳에서는 벽에서 자라는 경계층이 아직 얇아서 국소 마찰과 국소 열전달이 완전 발달 값보다 큽니다. 따라서 짧은 구간에서 측정한 평균 마찰계수나 Nu를 완전 발달 상관식과 비교하면 측정값이 체계적으로 높게 나옵니다. 완전 발달 결과를 인용하기 전에 확보한 직관 길이부터 확인하세요.
결과가 맞지 않을 때 확인할 순서
- L을 한 계산에서는 지름, 다른 계산에서는 관 길이로 썼는지 확인합니다.
- μ와 ν의 단위를 구분하고 비열이 kJ가 아닌 J 단위인지 확인합니다.
- α를 같은 ρ·cₚ·k에서 계산했는지, 또는 모든 물성을 물성표의 같은 행에서 가져왔는지 확인합니다. 서로 다른 자료의 값을 섞으면 등식이 넷째 자리에서부터 어긋날 수 있습니다.
- 열 Pe와 물질 Pe를 구분합니다. 물질 확산 계수 D로 계산했다면 대응 관계는 Re × Sc입니다.
- 같은 유효숫자 자릿수로 비교합니다. 다섯째 자리의 차이는 예제 2처럼 대개 물성표의 반올림에서 생깁니다.
이 결과로 알 수 없는 것
이 세 비율만으로 압력 손실, 대류 열전달 계수 h, 벽 온도 또는 전달 열량을 얻을 수는 없습니다. Nu = hL/k를 계산하려면 먼저 형상과 경계 조건에 맞는 측정값 또는 검증된 상관식으로 h를 정해야 합니다. 계산값과 함께 물성의 출처, 평가 온도, L의 정의를 기록해 두면 다른 계산과 같은 기준으로 비교할 수 있습니다.
Re·Pr·Pe 교차검산에서 자주 묻는 질문
μcₚ/k로 계산한 Pr이 물성표 값과 정확히 같지 않은 이유는 무엇인가요?
물성표는 열마다 따로 반올림해서 싣습니다. IAEA 20 °C 행으로 계산하면 μcₚ/k = 6.9808이고, 이를 반올림하면 표에 실린 6.981이 됩니다. 넷째나 다섯째 유효숫자의 차이는 정상입니다. 10% 이상 차이가 난다면 kJ와 J를 섞었거나 μ 칸에 ν를 넣은 단위 실수일 가능성이 큽니다.
L에는 관 지름과 관 길이 중 무엇을 넣어야 하나요?
관 유동의 Re와 Pe에는 안지름을 넣습니다. Pr에는 길이가 아예 들어가지 않습니다. 관 길이는 길이/지름 비나 입구 길이 Le/D 같은 별도의 비로 반영합니다. 예를 들어 길이 2 m, 지름 10 mm인 관이면 길이/지름 비는 200입니다. 한 수에는 지름을, 다른 수에는 관 길이를 쓰면 Pe = Re × Pr 등식이 깨집니다.
Re = 2000이면 항상 층류인가요?
그렇지 않습니다. 매끈한 원형 관에 대해서도 자료마다 기준값이 2100, 약 2300, 2320으로 다르고, 입구 교란이나 진동, 상류의 오리피스가 있으면 훨씬 낮은 Re에서도 난류가 생길 수 있습니다. Re = 2000은 입구가 잔잔하고 관이 길고 곧을 때 층류일 가능성이 높다는 뜻으로 읽고, 층류를 보장하는 값으로 보지 않는 편이 안전합니다.
물성은 어느 온도에서 평가해야 하나요?
가열이나 냉각이 있는 관에서는 입구와 출구 유체 온도의 평균을 흔히 사용합니다. 이 선택이 왜 중요한지는 예제 3에서 드러납니다. 20 °C에서 60 °C로 바꾸기만 해도 Re는 두 배가 되고 Pr은 절반 아래로 줄어듭니다. 벽 온도가 평균 유체 온도와 크게 다르면 벽 근처의 점성 변화가 속도 분포까지 바꾸므로 평균 온도 하나로 대표하기 어려워집니다. 결과마다 물성을 평가한 온도를 함께 적어 두세요.
같은 교차검산을 물질 전달에도 쓸 수 있나요?
쓸 수 있습니다. α 자리에 화학종의 확산 계수 D를 넣으면 물질 전달 페클레 수는 Re × Sc가 되고, 여기서 Sc = μ/(ρD)입니다. 액체에서는 열 Pe와 물질 Pe가 수십에서 수천 배까지 차이 날 수 있으므로 두 값을 다른 열에 나누어 기록하세요.
참고 자료
- NASA Glenn Research Center: Reynolds Number. Re를 관성력과 점성력의 비로 정의합니다.
- Penn State CE 574: 6.6 Dimensionless Numbers. Pe를 이송과 확산의 시간 척도로 비교하는 설명입니다.
- Michigan Technological University, CM3120: Module 2 Lecture IV. 열 Pe = Re Pr, 물질 Pe = Re Sc 관계를 확인했습니다.
- IAEA (2008): Thermophysical Properties of Materials for Nuclear Engineering. 표 3.3 포화수의 20 °C, 60 °C 행(ρ, μ, ν, cₚ, k, α, Pr)을 사용했습니다.
- MIT OpenCourseWare, HST.542J: Physiological Fluid Mechanics. 균일한 관의 천이 Re 약 2300, 오리피스 유동의 낮은 천이 Re, 층류 입구 길이 Le/D = 0.065 Re를 확인했습니다.
- Engineering LibreTexts (S. Miedema, TU Delft): 3.2 The Darcy-Weisbach Friction Factor. Re 2320 미만은 64/Re, 그 이상은 Colebrook–White 식을 씁니다.
- U.S. DOE WIPP: CRA-2009 Appendix PA. 층류·난류 경계를 Re 2100으로 두는 모델입니다.