입문 — off-design 계산이란 무엇인가 입문

1

무슨 문제인가 — 기계는 안 바뀐다, 조건만 바뀐다

이 기계는 이미 만들어져 있습니다. 터빈·열교환기·펌프·배관의 형상은 운전 중에 바뀌지 않습니다 — 터빈의 목(가장 좁은 통로) 넓이와 열교환기의 전열면적도 제작할 때 정해집니다.

바뀌는 건 바깥 조건입니다. 예를 들어 모델에서 냉각수 입구를 25 °C에서 30 °C로 바꾸거나, 열원이 약해지거나, 운전자가 펌프 회전수를 바꾸는 경우입니다. 여기서 25→30 °C는 실제 시험 기록이 아니라 조건 변화의 예입니다.

설계할 때 (순방향) "50 kW 를 내고 싶다"
→ 필요한 기기 형상과 성능을 정한다
→ 만든다
원하는 출력 → 기기 사양
off-design (역방향) "기기 사양은 정해져 있다. 날이 더워졌다"
그래서 얼마나 나오나?
기기 사양 + 바깥 조건 → 나오는 출력
off-design(탈설계) — 설계할 때 정한 그 조건이 아닌 상태에서 운전하는 것입니다. 실제 설비는 1년 내내 설계 조건으로 돌지 않으므로, 이 계산은 정해 둔 모델과 입력 안에서 그때의 상태와 출력을 추정합니다. 실제 성능과 같다고 말하려면 별도의 시험 자료로 검증해야 합니다.

이 강의와 현재 시나리오 코드가 다루는 바깥 조건 변화는 세 갈래입니다 — 증발기/열원 시나리오 · 냉각수 시나리오 · 펌프 회전수 시나리오. 마지막 절에서 이 이름으로 결과를 얼마나 믿을지 가립니다.

2

냉매가 도는 사이클과 12개 지점

냉매(R1336mzz(Z))가 사이클을 한 바퀴 돕니다. 경로 위에서 압력·온도 같은 상태를 읽는 12개 위치에 ①~⑫를 붙였습니다. 기기 번호가 아니라 냉매 상태점 번호입니다.

고압측 — 데우고 끓인다 (약 14.8 ~ 15.8 bar) 저압측 — 식히고 다시 밀어 올린다 (약 1.65 ~ 1.85 bar) 재생기 기기는 하나 찬 쪽 ③ → ④ 뜨거운 쪽 ⑩ ← ⑨ 증발기⑤ → ⑥ 터빈 ★ 응축기⑫ ← ⑪ 펌프 펌프가 밀어 올린 액체가 재생기로 터빈에서 나온 증기가 재생기로 배가스 (열원) 냉각수
재생기가 지점 번호를 네 개 갖는 이유 — 기기는 하나인데 유체가 두 줄기 지나가기 때문입니다. 터빈에서 나온 뜨거운 증기(⑨→⑩)로 펌프에서 나온 찬 액체(③→④)를 미리 데웁니다 — 버릴 열을 다시 쓰는 것이라 재생입니다.

사이클은 펌프와 터빈을 경계로 둘로 나뉩니다. 펌프가 압력을 올리고 터빈이 내립니다. 그 사이 고압측(파란 선)에서 데우고 끓이고, 저압측(빨간 선)에서 식힙니다.

그림에서 번호가 건너뛰는 구간은 빠진 상태점이 아닙니다. 예를 들어 ②→③은 펌프 출구에서 재생기 찬 쪽 입구까지, ⑧→⑨는 터빈 출구에서 재생기 뜨거운 쪽 입구까지의 연결 구간입니다. 화살표를 따라 ①부터 ⑫까지 읽으면 냉매가 다시 ①로 돌아옵니다.

배가스(열원 — 연소를 마친 뜨거운 가스)는 증발기에서 냉매에 열을 건네고, 냉각수는 응축기에서 냉매가 버리는 열을 받아 나갑니다.

3

왜 바로 못 푸나 — 사이클이 닫혀 있다

계산을 시작하려면 어딘가 아는 값이 있어야 합니다. 그런데 여기는 출발점이 없습니다.

유량을 알아야압력강하를 안다 압력을 알아야터빈을 푼다 터빈을 풀어야유량이 나온다 제자리로 돌아온다 — 어디서도 시작할 수 없다
서로가 서로를 필요로 해서 순서대로 풀 수가 없습니다 — 이런 것을 연립방정식(따로따로가 아니라 한꺼번에 풀어야 하는 식들)이라고 합니다.

그래서 방법을 바꿉니다. 모르는 값을 일단 추정해서 넣어 봅니다. 그 값이 맞다고 치고 끝까지 계산해 본 뒤, 결과가 처음 추정한 값과 맞는지 확인합니다.

12지점 × (압력·온도) = 최대 24개를 다 찾아야 할 것 같지만, 아는 법칙들이 하나씩 지워줍니다:

법칙지워주는 것
질량보존 — 사이클에 갈림길이 없다유량은 어디서나 같다 → 유량 미지수 1개
구간 압력강하 = 그 구간의 설계 압력강하 × 배율 ⚠️(압력강하는 마찰·굽힘·밸브·기기 내부 등 구간을 지나며 잃는 압력입니다. 배율은 지금 유량 ÷ 설계 유량에서 나옵니다 — 1단계에서 구합니다)한 지점 압력만 알면 이웃 압력이 차례로 나온다
포화관계 — 끓는 중이면 압력↔온도 1:1⑫ 압력만 알면 ⑫ 온도가 나온다(+ 과냉도 고정 ⚠️ — 과냉도란 응축기 출구에서 냉매가 끓는 온도보다 몇 도 더 식었는지입니다)
물성표 — 냉매의 성질을 담은 표(REFPROP. 냉매 조성과 유효한 독립 상태량 두 개를 주면 다른 상태량을 조회합니다)현재 모델이 쓰는 유효한 입력쌍으로 나머지 상태량을 계산

남는 미지수는 4개입니다. 바깥에서 되풀이해 맞추는 L1·L2·L3와, 각 시험값을 평가할 때 8단계 안에서 다시 푸는 L0입니다.

미지수왜 모르나채점 기준
L1냉매 유량 조건이 바뀌면 터빈을 지나는 양이 달라짐 현재 시험 유량 = 터빈 모델이 계산한 통과 유량 (질량유량이 스스로 맞나)
L2터빈 입구 ⑦ 온도 증발기가 데워 내보내는 온도는 유량·열원·증발기 모델에 따라 달라짐 총괄 UA(시나리오 기본): 필요 UA = 모델 UA. 도면(옵션): 1차원이 준 ⑥ 엔탈피로 ⑦ 온도를 대입
L3응축기 출구 ⑫ 압력 냉각수가 뜨거워지면 ⑫ 압력이 올라감 현재 상태에 필요한 응축기 UA = 모델 응축기 UA (열전달 능력에 맞나)
L0재생기 열량 열량↔출구온도가 서로 물려 있음 현재 열량에 필요한 재생기 UA = 모델 재생기 UA 가 되도록 8단계 안에서 그 자리에서 푼다. 바깥 L1·L2·L3 반복에는 남기지 않는다
왜 현재 계산은 ⑫ 압력을 시험하는가 — ⑫는 응축기 출구입니다. 이 압력에서 냉매의 응축 포화온도를 찾고 고정 과냉도를 빼면 ⑫ 온도와 액체 상태가 정해집니다. 또한 ⑫ 압력에서 저압측 압력강하를 따라 ①·⑧ 압력까지 계산할 수 있어, 응축기 UA 잔차를 움직이는 L3 변수로 사용합니다. 응축기는 이 solved 경로입니다. 증발기는 총괄 UA 또는 도면입니다.
응축기 L3는 늘 필요 UA = 모델 UA로 채점합니다. 증발기는 경로가 둘입니다. 시나리오 기본인 총괄 UA 경로는 증발기 모듈이 같은 채점을 안에서 되풀이해 출구를 찾고 ⑦ 온도를 대입하며, 옵션인 도면 경로는 1차원이 ⑦ 온도를 대입합니다. 1차원 ΣUoA와 LMTD 필요 UA는 정의가 달라 맞추지 않습니다.
⑥ 엔탈피를 시험한다증발기 모듈 안 옮길 열량입구·출구 엔탈피 차 ÷ 양쪽 온도차배가스 온도는 아는 값 = 필요한 UA열량 ÷ 평균온도차 모델 증발기 UA이 운전점에 배정된 값 비교 두 UA가 다르면 → 출구 엔탈피를 옮겨 다시 (모듈 안 이분법)

위 그림은 총괄 UA 경로(시나리오 기본)입니다. 두 UA 가 맞는 출구가 나오면 ⑦ 온도를 대입합니다. 옵션인 도면 경로는 1차원 모델이 ⑥ 엔탈피를 직접 줍니다.

총괄 UA 경로 — 시나리오 기본 증발기 모듈이 “필요 UA = 모델 UA”인 출구 ⑥ 엔탈피를 안에서 찾아 ⑦ 온도를 대입한다. 바깥에서 ⑦ 온도를 따로 맞추지 않는다.
도면 경로 — 옵션 도면 1차원이 ⑥ 엔탈피와 증발기 Δp를 준다. ⑦ 온도는 대입한다. LMTD UA로 맞추지 않는다.
⚠️ 과냉도 고정 · 응축기 모델 UA · 냉각수 조건은 가정 🔶 냉각수·열원의 펌프 고정 2800 rpm은 모델이 닫히는 값 — 논문 설계속도가 아님
4

추정하고, 한 바퀴 돌고, 채점하고, 다시 추정한다

앞 절에서 추정하기로 한 세 값(L1·L2·L3)을 넣고 한 바퀴 돌리면, 그 값이 맞았는지 채점할 재료가 나옵니다.

세 값을추정한다 1∼9단계 한 번 평가한다 채점 세 값이 안 움직이나? 안 맞으면 값을 고쳐서 다시 — 셋 다 안 움직일 때까지

여기서 한 번의 평가는 9단계만 실행한다는 뜻이 아닙니다. 1단계부터 9단계까지 전부 다시 계산하고, 8단계에서는 재생기 열량 L0도 현재 조건으로 다시 풉니다. L1·L2·L3을 하나씩 고칠 때마다 같은 평가가 반복됩니다.

5

한 바퀴에 무엇을 하나 — 9단계

한 바퀴는 9단계입니다. 큰 줄기만 보면 이렇습니다.

단계무엇을 하나
1지금 추정해 둔 유량이 설계의 몇 배인지 구한다 — 이 비율로 압력강하와 재생기 UA 를 보정한다
2저압측 압력을 다 구한다. ⑫ 하나를 기준으로 양쪽으로 퍼뜨린다
3⑫ 압력에서 응축 포화온도를 찾고, 고정 과냉도를 빼서 ⑫의 액체 상태를 만든다
4펌프의 현재 유량·회전수와 QH 곡선으로 압력 상승과 필요한 동력을 계산한다
5고압측 압력을 다 구한다
6총괄 UA면 L2가 둔 ⑦ 시험 온도로 ⑦ 상태를 만든다. 도면 경로면 1차원이 준 ⑥ 엔탈피로 ⑦ 온도를 대입한다
7터빈 — 냉매 1 kg당 실제 에너지 낙차와 ⑧ 상태를 구하고, ★ 모델 통과 유량을 계산한다
8재생기 열량 L0를 안쪽 반복으로 풀어 ④와 ⑩의 상태를 정한다
9유량과 엔탈피 차로 열량·터빈·펌프 유체동력·순유체동력을 계산한다. 응축기 UA 잔차는 L3 채점. 증발기 UA 잔차는 총괄 UA 경로에서만 L2 채점에 쓴다
7단계가 핵심입니다. 여기서 터빈 모델이 계산한 유량이 나옵니다. 터빈 노즐에는 가장 좁은 통로인 이 있으며, 입구 조건과 유효 목면적이 정해지면 현재 노즐 모델이 통과 유량을 계산합니다. 이 값이 처음 추정한 유량과 맞는지가 L1 채점입니다.
터빈의 총동력은 7단계가 아니라 9단계에서 계산합니다. 7단계가 만드는 것은 ⑧ 상태와 냉매 1 kg당 에너지 낙차, 그리고 터빈 모델 유량입니다. 9단계가 최종 유량을 곱해 kW 단위의 터빈동력을 만듭니다.
6

답을 어디까지 믿나

계산이 수렴했다고(= 값이 안 움직인다고) 그 답이 물리적으로 가능한 것은 아닙니다. 수학적으로는 풀리는데 기계가 부서지는 답일 수 있습니다.

그래서 수렴한 뒤에 물리 검사를 합니다 — 터빈에 액체가 들어가지는 않는지, 열이 차가운 쪽에서 뜨거운 쪽으로 거꾸로 흐르지는 않는지 같은 것들입니다.

수렴·물리 가능성·실기 신뢰도는 서로 다릅니다. 계산이 수렴하고 물리 검사를 통과해도, 넣은 자료와 모델이 검증되지 않았다면 실제 설비의 정확한 성능이라고 할 수 없습니다.

아래 등급은 계산 성공 여부가 아니라 현재 자료로 결과를 해석할 수 있는 범위입니다. 「경향만」은 조건을 바꿔 가며 본 추세까지만, 「방향만」은 결과가 오르내리는 방향까지만 사용한다는 뜻입니다.

등급어떤 결과가
모델 경향만증발기/열원 시나리오배가스 열용량률은 그림 판독값. 총괄 UA 경로가 필요 UA = 모델 UA로 ⑦ 온도를 맞추는데, 그 UA 자체가 설계점 역산값
모델 방향만냉각수 시나리오 · 응축기 관련 값냉각수 조건이 논문에 없어 가정. 펌프는 2800 rpm으로 고정 — 논문 설계속도가 아니라 모델이 닫히게 고른 값
절대값 금지펌프 회전수 시나리오제조사 3539 rpm 원곡선은 3204 rpm 시험 양정을 31.3% 높게 예측. 3539·3204는 설계 UA와 만나면 근이 없고, 시나리오는 2800 근처만 훑는다
🔴 순유체동력 54.312 kW(터빈 57.002 − 펌프 2.690)는 총괄 UA solved 기준점의 유체 계산값입니다. 펌프까지 곡선으로 푼 시나리오 기준점 값도, 전기출력 kWe도 아닙니다. 논문 전기 피크 46.7 kW 와 직접 비교하지 않습니다 (이유는 표준 편).
🔴 터빈 84 %: 달성 여부가 확인되지 않은 모델값 🔴 냉각수 조건: 논문에 없어 가정 🔵 펌프 QH: 제조사 3539 rpm이 곡선 r=1 · 3600은 논문 명목값 · 3204는 시험점 · 시나리오 고정은 2800 🔵 배가스: 335 °C만 논문 인쇄값, 열용량률은 그림 판독값
7

다음은 어디로

표준으로 넘어가기 전 확인 — ①∼⑫는 냉매 상태점, 1∼9는 한 번의 계산 순서, L1∼L3은 바깥 시험값, L0는 8단계 안에서 푸는 재생기 열량입니다. 이 네 종류를 구분할 수 있으면 다음 문서를 읽을 준비가 된 것입니다.