[README 학술 검토] 수식 오류, 논리 모순, 근거 부족 항목 정리

[chestnutbread]

🔴 학술적 오류 (명백한 불일치)

1. 압축기 출구 온도 계산식과 결과값 불일치

README에 명시된 공식:

T₂ = T₁ × (1 + (1/η) × (PR^((γ-1)/γ) - 1))

실제 계산:

8^(0.4/1.4) = 8^0.2857 ≈ 1.811
T₂ = 288.15 × (1 + (1/0.85) × (1.811 - 1))
   = 288.15 × (1 + 1.176 × 0.811)
   = 288.15 × 1.954
   ≈ 563 K  (← README는 ≈613 K라고 주장)

시뮬레이션 값 613.75 K와 공식 결과 563 K 사이에 약 50 K 오차 존재.
시뮬레이션이 613.75 K를 산출하려면 η=0.85 기준에서 PR이 8이 아니라 약 10.6 이어야 함.
→ 코드 내부 실제 파라미터와 README 명시값이 다를 가능성 있음.

---

2. 터빈 출구 온도 판정 기준 모순

기준	시뮬레이션 값	판정
T₄ < 923 K	924.87 K	✅ (오차 0.2%)

924.87 K > 923 K → 기준 미달임에도 ✅ 처리.
"오차 0.2%"는 수치적 근접도를 나타낼 뿐, 기준 통과의 근거가 아님.
→ 기준값이 924 K였다면 말이 되지만, < 923 K라고 명시한 이상 내부 모순.

---

3. 노즐 면적 역산 논리 모순 (섹션 4.2)

• "p_e/p_in = 1.0이면 유량함수 = 0 → 유량이 0이 됨" 이라고 스스로 인정
• 바로 다음 줄에서 "mdot ∝ area (선형 비례 관계 성립)"이라고 결론

유량이 0인 조건에서 비례 관계가 어떻게 성립하는지 설명 없음.
또한 4.4절에서 초기 area = 0.0637 m²에서 mflow = 18.92 kg/s가 역산됐다고 했는데, 이는 4.2절의 "유량 = 0" 결론과 충돌.

---

4. 터빈 동력 계산에서 질량유량 불일치

• 터빈 출구 온도 공식: T₄ = T₃ + ΔT_comp × η_m × (1 + FAR) → FAR로 질량 증가 반영
• 터빈 동력 계산: W_turb = ṁ × cp × (T₃ - T₄) → ṁ = 3.493 kg/s (FAR 미반영)

실제 터빈 통과 질량유량: ṁ_gas = 3.493 × 1.019 = 3.559 kg/s 여야 함.
온도식에서는 (1+FAR)을 반영하면서 동력 계산에서는 미반영 → 내부 일관성 없음.

---

🟡 근거 부족 / 질문 예상 항목

5. 시뮬레이션 구간 200s의 근거 없음

• 엔진 시상수(time constant) 계산 없음 (τ = J/b 형태)
• "정상상태 수렴" 판정 기준치 없음 (예: 신호 변동 < 0.01% for 10s 이상 등)
• UAV 순항 임무에서 해당 조건이 실제로 200s 유지된다는 운용 맥락 없음

예상 질문: "200s가 충분한 이유는? 300s에서도 동일한 결과라는 보장이 있나요?"

---

6. FAR = 0.019 이론값과 불일치

연소기 에너지 평형 기반 이론 FAR:

FAR = cp_air × (T₃ - T₂) / (η_comb × LHV_fuel)
    = 1005 × (1250 - 613.75) / (0.96 × 43,200,000)
    ≈ 0.0154

코드 사용값 0.019 vs 이론값 0.0154 → 약 23% 차이.
LHV 가정이 다른지, 코드가 별도 방식으로 FAR를 계산하는지 설명 없음.

예상 질문: "FAR 0.019는 어디서 나온 값인가요?"

---

7. 핵심 설계 파라미터 출처 전무

파라미터	값	문제
공기 유량 ṁ	3.493 kg/s	"TPE331-10 설계점"이라 하나 출처(스펙시트, 논문 등) 없음
압력비 PR	8	TPE331-10 실제 PR은 문헌상 910 수준, 근거 없음
η_comp	0.85	출처 없음
η_t	0.87	출처 없음
TIT	1250 K	실제 TPE331-10 TIT와 비교 없음

---

8. shaft inertia = 1.0 단위 불명확

단위 없이 "inertia=1.0"으로만 명시. 관성 모멘트 단위는 kg·m²이며, 실제 터보프롭 샤프트는 0.05~2 kg·m² 범위.
→ 값의 물리적 근거 없으면 placeholder로 의심받을 수 있음.

---

9. omega₀ = 4370 rad/s = 41,730 RPM의 정체 불명확

• 터보프롭 프로펠러 축: 통상 1,000~2,500 RPM → 41,730 RPM은 비현실적
• 가스발생기(gas generator) 축: 35,000~55,000 RPM → 41,730 RPM은 합리적

어느 축의 속도인지 명시 없음. 감속기어(reduction gearbox) 모델링도 없음.

예상 질문: "41,730 RPM이면 프로펠러 축이 아니라 gas generator 축 아닌가요?"

---

10. 연소가스 γ = 1.31 온도 기준 없음

γ는 온도에 따라 연속 변화 (1000K ≈ 1.32, 1500K ≈ 1.28). 터빈 입구(1250K)~출구(924K) 구간 어느 온도 기준인지, 어떤 문헌을 참조했는지 설명 없음.

---

11. turbine₁.t_in 모니터링 "구조적 제약" 주장 검증 불가

"VHDL-AMS port in quantity는 Twin Builder Results에서 0으로 표시되는 구조적 제약" → Ansys 공식 문서 인용 없음.
역산으로 1250 K를 증명하는 논리는 잘 구성되어 있으나, 제약 자체의 공식 근거가 없으면 심사에서 의심받을 수 있음.

---

요약

분류	항목	심각도
수식 오류	압축기 공식 결과 불일치 (563K vs 613K)	🔴 높음
논리 모순	터빈 T₄ < 923K 기준에 924.87K 합격 처리	🔴 높음
논리 모순	노즐 역산: 유량=0 주장 후 비례관계 성립 주장	🔴 높음
내적 불일치	터빈 동력 계산에서 (1+FAR) 미반영	🟡 중간
근거 부족	200s 시뮬레이션 구간 정당화 없음	🟡 중간
근거 부족	FAR=0.019 이론값과 23% 불일치	🟡 중간
미인용	ṁ, PR, η, TIT 파라미터 출처 없음	🟡 중간
불명확	inertia=1.0 단위 없음	🟡 중간
불명확	41,730 RPM이 어떤 축인지 불명	🟡 중간
불명확	γ=1.31 온도 기준 없음	🟡 중간
검증 불가	port in quantity 제약 공식 문서 없음	🟡 중간

[chestnutbread]

🔧 수정 가능 항목 및 수정 방향

Honeywell TPE331-10 공식 브로셔 (N61-1491-000-000) 및 추가 조사 결과를 바탕으로 수정 가능한 항목을 정리합니다.

---

✅ Fix 1: PR = 8 → 10.55 수정 (이슈 #1 해소)

공식 스펙: Pressure ratio: 10.55 (Honeywell 브로셔)

현재 README:

| Compressor | 등엔트로피 압축 | PR=8, η_comp=0.85, γ=1.4 |

수정:

| Compressor | 등엔트로피 압축 | PR=10.55, η_comp=0.85, γ=1.4 |

근거:
PR=10.55 대입 시 압축기 출구 온도 재계산:

10.55^(0.4/1.4) = 10.55^0.2857 ≈ 1.960
T₂ = 288.15 × (1 + (1/0.85) × (1.960 - 1))
   = 288.15 × 2.130 ≈ 613.8 K

→ 시뮬레이션 값 613.75 K와 일치. 공식은 맞고 PR 표기만 잘못된 것이었음.

---

✅ Fix 2: ṁ = 3.493 kg/s 출처 추가 (이슈 #7 부분 해소)

공식 스펙: Airflow: 7.7 lb/sec → 7.7 × 0.4536 = 3.493 kg/s (정확히 일치)

현재 출처 없이 "TPE331-10 설계점"이라고만 되어 있는 부분에 출처 추가:

출처: Honeywell TPE331-10 Turboprop Engine Brochure (N61-1491-000-000, 2016)

---

✅ Fix 3: omega₀ = 4370 rad/s 설명 보완 (이슈 #9 해소)

공식 스펙: RPM: gas gen 41,730 / shaft output 2,000 (cw) / 1,591 (ccw)

현재 README:

| Shaft | 터빈-압축기 기계적 연결 | inertia=1.0, omega0=4370 rad/s |

수정: omega₀가 가스발생기(gas generator) 축 속도임을 명시하고, 프로펠러 출력축과의 관계 추가:

omega₀ = 4370 rad/s = 41,730 RPM (가스발생기 축, Honeywell 공식 스펙과 일치)
프로펠러 출력축: 2,000 RPM (CW) / 감속비 ≈ 20.9:1

---

✅ Fix 4: T₄ 수렴 기준값 수정 (이슈 #2 해소)

현재: turbine1.t_out < 923 K 기준, 결과 924.87 K → ✅ (모순)

수정 방향 (두 가지 중 선택):

Option A — 기준값을 시뮬레이션 결과에 맞게 수정:

| 터빈 출구 온도 | ~925 K | 924.87 K | ✅ |

이론 근거 추가:

T₄ = T₃ + ΔT_comp × η_m × (1 + FAR)
   = 1250 + (-325.6) × 0.98 × 1.019 ≈ 924.85 K

→ 코드 수식에서 도출된 이론 예측값 ~925 K을 기준으로 설정

Option B — < 923 K를 유지하되 오차 허용 범위 명시:

기준: 925 K 이하 (이론값 대비 ±2 K 허용)
결과: 924.87 K → ✅

---

✅ Fix 5: 노즐 역산 논리 수정 (이슈 #3 해소)

현재 4.2절 문제: "p_e = p_in → 유량 = 0" 인정 후 "mdot ∝ area" 결론 → 자기 모순

수정: p_e ≠ p_in인 실제 상황을 명시하는 방식으로 재작성:

※ 해면고도 정상 작동 시 터빈 출구압력은 대기압보다 높음.
실제 p_in (터빈 출구) > p_back (= p_amb = 101,325 Pa) 이므로 압력차 존재.
이 조건에서 아음속 유량식 적용 시 mdot ∝ area 성립.
p_e = p_in = 101,325 Pa는 이상적 한계 조건이며, 실 모델에서는 터빈 출구압이 더 높아 유량이 정상 발생함.

---

✅ Fix 6: 터빈 동력 계산 질량유량 일관성 수정 (이슈 #4 해소)

현재 불일치:

• 온도식: (1 + FAR) 반영
• 동력식: ṁ = 3.493 kg/s (FAR 미반영)

수정: 동력식에 연료 질량 반영:

ṁ_gas = ṁ_air × (1 + FAR) = 3.493 × 1.019 = 3.560 kg/s
W_turb = ṁ_gas × cp_gas × (T₃ - T₄)
       = 3.560 × 1212.7 × 325.13
       ≈ 1,404,800 W ≈ 1.40 MW

단, 코드 내부에서 ṁ_air 기준으로 동력 계산 시 그 이유(연료 질량 분리 처리 등)를 명시할 것.

---

⚠️ 추가 출처 확보 필요 (조사 중)

항목	상태
η_comp = 0.85, η_t = 0.87	교과서(Mattingly, Cumpsty 등) 인용 필요
TIT = 1250 K	파일럿 매뉴얼 또는 논문 참조 필요
shaft inertia = 1.0 kg·m²	단위 명시 + 물리적 근거 필요
port in quantity 제약	Ansys 포럼/공식 문서 인용 필요
200s 시뮬레이션 근거	시상수 계산 추가 필요

[chestnutbread]

[파라미터 수정] TPE331-10_Research_Parameters_v2 기준 비교 결과

✅ Python code/normal.py 수정 완료 (commit: 2025-06-30)

연구 파라미터 기준으로 normal.py 업데이트됨.

파라미터	구 버전	신 버전 (연구기준)	출처
PR	10.55 (VHDL 역산)	9.6	Garrett TPE331-10; Saravanamuttoo 2017 부록A
η_c	0.85	0.80	Saravanamuttoo 2017 Ch.5
η_b	0.96	0.990	JSSG-2007A §3.2.2.6
FAR	0.019	0.018	Mattingly 2006 Ch.9
TIT	1250 K (고정)	열역학 계산 (~1364 K)	연소기 에너지 평형
η_t	0.87	0.88	Saravanamuttoo 2017 Ch.7
η_pt	—	0.87	Park et al. 2023, Energy 305

---

🔴 VHDL Twin Builder README 수정 필요 항목

VHDL 코드는 수정하지 않고 README에 아래 불일치를 명시해야 함.

1. PR (압력비) — 명백한 기재 오류

• 현재 README: PR = 8
• 연구 기준: PR = 9.6 (Garrett TPE331-10; Saravanamuttoo 2017 부록A)
• VHDL 코드 실제 동작: T₂ = 613.75 K 역산 시 PR ≈ 10.55 (Honeywell FAA TCDS E4WE)
• 권고: README를 PR = 9.6 (문헌 기준) 또는 10.55 (Honeywell 공식)로 수정

2. η_c (압축기 효율)

• 현재 README: η_c = 0.85
• 연구 기준: η_c = 0.80 (Saravanamuttoo 2017 Ch.5, 2단 원심 0.78~0.82)
• 권고: 0.80으로 수정하고 문헌 인용 명시

3. η_b (연소 효율)

• 현재 README: η_b = 0.96
• 연구 기준: η_b = 0.990 (JSSG-2007A §3.2.2.6 신규 엔진 기준)
• 권고: 0.990으로 수정

4. FAR (연료-공기 비)

• 현재 README: FAR = 0.019
• 연구 기준: FAR = 0.018 (Mattingly 2006 Ch.9)
• 권고: 0.018로 수정

5. TIT (터빈 입구 온도) — 구조적 불일치

• 현재 VHDL: TIT = 1250 K 고정 설계값
• 연구 기준 계산: FAR=0.018, η_b=0.990 → T₃ ≈ 1364 K
• TIT=1250 K에 대응하는 이론 FAR: ≈ 0.0157 (VHDL FAR=0.019와 불일치)
• 권고: README에 "TIT = 1250 K는 설계점 고정값이며, 열역학적 평형 TIT (~1364 K)와 상이함" 명시

6. η_t (가스터빈 효율)

• 현재 README: η_t = 0.87
• 연구 기준: η_t = 0.88 (Saravanamuttoo 2017 Ch.7 대표값)
• 권고: 0.88로 수정 (0.87은 η_pt 값과 혼동 가능)

---

📊 파라미터 변경에 따른 출력값 영향

연구 기준 파라미터 적용 시 Python 계산 결과:

항목	연구기준 (v2)	VHDL 실측	차이
T₂ (압축기 출구)	~615 K	613.75 K	+1 K
T₃ (TIT)	~1364 K	1250 K	+114 K
T₄ (터빈 출구)	~1038 K	924.87 K	+113 K
W_comp	~1168 kW	1150 kW	+18 kW

→ TIT 차이(+114 K)가 주요 불일치 원인. VHDL이 TIT를 고정하는 구조적 특성에서 기인.