계산화학 시뮬레이션 물리 마지막 수정: 2026-08-26 04:48:09

1. 개요[편집]

노킹
Engine Knock
발생 위치미연 말단가스(end gas)
원인화염 도달 전 자착화
소리의 정체실린더 기체 공진 — 대개 5~10 kHz
연료 지표RON · MON (PRF 기준)
극단 형태초노킹(super-knock) · LSPI
표준 예측리븐굿-우 적분
표준 대응노크센서 기반 점화 지각

엔진이 “따그닥” 소리를 낸다면 그건 부품이 헐거워서가 아니라, 연료가 제멋대로 먼저 터지고 있다는 뜻이다.

노킹(knocking, engine knock)은 스파크 점화 엔진에서 점화 플러그가 만든 화염면이 도달하기 전에 미연 혼합기 — 말단가스(end gas) — 가 스스로 자착화하여, 급격한 열방출과 압력파를 만드는 이상연소 현상이다. 금속성 타격음이 나기 때문에 붙은 이름이고, 심하면 피스톤 링랜드가 깨지고 헤드 개스킷이 날아간다.

원리 자체는 자착화의 응용 문제 하나에 불과하다. 하지만 그 응용 문제가 내연기관의 압축비 상한을 정하고, 따라서 열효율의 상한을 정하고, 결국 주유소 간판의 숫자를 정한다. 100년 넘게 연구되고도 여전히 엔진 캘리브레이션 공수의 상당 부분을 잡아먹는 주제다.

화학은 자착화냉염에, 그 화학을 재는 장비는 급속압축기에 정리돼 있다. 여기서는 엔진 안에서 무슨 일이 벌어지고, 그것을 어떻게 듣고 어떻게 예측하고 어떻게 피하는가만 다룬다.

2. 정상 연소와 무엇이 다른가[편집]

정상적인 SI 엔진 연소는 점화 플러그에서 시작한 화염핵이 난류 화염면으로 자라 연소실 전체를 훑는 전파 과정이다. 화염면이 지나간 뒤 기연가스가 팽창하면서, 아직 안 탄 말단가스를 압축한다. 즉 말단가스는 피스톤에 한 번, 화염 팽창에 또 한 번 압축당한다.

문제는 이 압축이 말단가스를 데운다는 것. 화염면이 도착하기 전에 그 온도·압력 이력의 누적이 임계를 넘으면, 말단가스가 자기 스스로 착화한다. 화염 전파에 수 밀리초가 걸리는 반면 자착화는 거의 동시에 넓은 영역에서 일어나므로, 그만큼의 열이 한꺼번에 풀린다. 급격한 압력 상승 → 압력파 → 벽면 반사 → 공진. 여기까지가 노킹이다.

정리하면 노킹의 조건은 경쟁이다.

τig(말단가스)  <  t화염 도달\tau_{ig}(\text{말단가스}) \;<\; t_{\text{화염 도달}}

화염이 빨리 오거나(tt 감소), 말단가스가 잘 안 터지면(τig\tau_{ig} 증가) 노킹이 안 난다. 그래서 연소 기간을 줄이는 텀블·스월 강화, 화염 전파 거리를 줄이는 중앙 점화, 말단가스를 식히는 냉각수로 배치, EGR로 온도·반응성 낮추기가 전부 같은 부등식을 겨냥한 설계다.

3. 노킹음의 정체[편집]

노킹 소리는 “폭발음”이 아니다. 자착화 자체는 순간이고, 우리가 듣는 것은 그 뒤 실린더 안 기체가 울리는 음향 공진이다.

1938년 드레이퍼(C. S. Draper)가 원통형 경계조건에서 파동방정식을 풀어 이 공진 모드를 정리했다. 지름 BB인 원통 공동의 모드 주파수는 대략

fm,n=ρm,ncπBf_{m,n} = \frac{\rho_{m,n}\, c}{\pi B}

이고, ρm,n\rho_{m,n}은 베셀 함수에서 오는 모드 상수, cc는 기연가스 음속이다. 1차 원주 모드(m=1,n=0m{=}1, n{=}0)의 ρ1.841\rho \approx 1.841이 가장 낮고 가장 잘 여기된다. 승용차 보어(80 mm 남짓)와 기연가스 음속(800 m/s 안팎)을 넣으면 5~7 kHz 정도가 나온다. 고차 모드까지 세면 실측 스펙트럼은 대략 6~20 kHz에 걸쳐 봉우리를 여럿 만든다.

이 구조 덕분에 노킹은 주파수로 잡을 수 있다. 실차 노크센서는 블록에 볼트로 붙인 압전 가속도계이고, ECU는 이렇게 처리한다.

  • 연소 구간에 해당하는 크랭크각 창(window)만 잘라낸다 — 흡배기 밸브 착좌음이나 인젝터 소음을 배제하기 위해.
  • 그 구간을 공진 대역으로 대역통과 필터링하고 에너지를 적분한다. 고속 푸리에 변환이나 필터뱅크로 처리하며, 과도 신호라 웨이블릿 변환 기반 검출도 쓰인다.
  • 그 값을 최근 사이클들의 배경 소음 수준과 비교해 임계 초과 여부를 판정한다.

임계를 넘으면 즉시 점화시기를 몇 도 지각(knock retard)시키고, 노킹이 사라지면 천천히 다시 진각한다. **“빠르게 물러서고 느리게 돌아온다”**가 이 폐루프의 국룰이며, 이 비대칭이 없으면 노킹 경계 근처에서 진동이 생긴다. 여기서 주파수가 온도에 의존한다는 점도 성가시다 — 음속 cTc \propto \sqrt{T}라 공진 주파수가 운전 조건과 사이클마다 흔들린다. 그래서 고정 대역 필터 대신 공진 주파수를 사이클마다 추정하는 방식도 쓰인다.1

4. 세기의 스펙트럼 — 초노킹과 LSPI[편집]

노킹은 있고 없고의 이분법이 아니라 세기의 연속체다.

  • 경노킹(light knock) — 압력 진동 진폭이 작다. 사실 연비에는 유리한 영역이라, 실무 캘리브레이션은 노킹 없는 조건이 아니라 “허용 가능한 노킹” 경계 바로 아래를 노린다.
  • 중노킹 — 지속되면 열부하가 올라 피스톤 상단·링랜드가 손상된다.
  • 초노킹(super-knock) — 통상 노킹과 자릿수가 다르다. 실린더 압력 첨두가 정상의 몇 배까지 튀고, 단 몇 사이클로 엔진을 파괴한다.

초노킹은 다운사이징 과급 직분사(TGDI) 엔진에서 문제가 됐다. 대개 저속조기점화(Low-Speed Pre-Ignition, LSPI)에서 출발한다 — 저회전·고부하(대략 BMEP 15 bar 이상, 1500~2500 rpm 영역)에서 스파크가 튀기도 전에 무언가가 혼합기를 점화시킨다. 유력한 개시원은 피스톤 크레비스에서 빠져나온 엔진오일 액적이나, 연소실 퇴적물이 떨어져 나온 조각이다. 오일 액적은 국소적으로 옥탄가를 떨어뜨리는 셈이라 그 자리에서 자착화한다.

일단 조기점화가 나면 화염이 훨씬 이른 시점부터 퍼지므로 말단가스가 더 심하게 압축되고, 그 결과가 극단적 노킹 → 폭굉으로의 발전이다. 그래서 LSPI 대책은 연소 제어가 아니라 엔진오일 첨가제 조성(칼슘계 청정제 저감 등)에서 나왔다. 화학이 문제를 만들었으니 화학으로 막는다.2

5. 폭연-폭굉 천이와 브래들리 다이어그램[편집]

말단가스 자착화가 왜 어떤 때는 얌전하고 어떤 때는 엔진을 부수는가. 답은 자착화가 어떻게 전파하는가에 있다.

자착화 문서에 나온 젤도비치 구배 메커니즘이 골격이다. 반응성이 공간적으로 분포하면 착화가 가장 빠른 지점(핫스팟)에서 먼저 터지고, 그 사건이 이웃으로 번지는 속도는 화염 속도가 아니라

usp=τigx1u_{sp} = \left|\frac{\partial \tau_{ig}}{\partial x}\right|^{-1}

이다. 이 속도가 음속과 비슷해지면 반응이 압력파와 공명해 서로를 증폭시키고, 폭연-폭굉 천이(DDT)가 일어난다.

브래들리 그룹은 이 판정을 무차원 두 개로 압축했다.

ξ=ausp,ε=r0aτe\xi = \frac{a}{u_{sp}}, \qquad \varepsilon = \frac{r_0}{a\,\tau_e}
  • ξ\xi — 음속과 자착화 전파 속도의 비. 너무 크면(전파가 느리면) 그냥 순차 착화로 끝나고, 너무 작으면(전파가 음속보다 훨씬 빠르면) 사실상 동시 폭발이라 압력파가 조직되지 않는다.
  • ε\varepsilon — 핫스팟 크기 r0r_0를 음파가 가로지르는 시간과, 화학에너지의 대부분이 풀리는 여기 시간(excitation time) τe\tau_e의 비. 이 값이 커야 음파가 핫스팟 안에 머무는 동안 충분한 열이 공급된다.

ξ\xiε\varepsilon 평면에 실험·계산 결과를 찍으면 폭굉이 발생하는 영역이 반도(peninsula) 모양으로 닫힌 곳에 몰린다 — 폭굉 반도(detonation peninsula)다. 실무적 함의가 크다. τe\tau_eτig\tau_{ig}와 다른 양이며, 착화지연이 같아도 여기 시간이 짧은 연료가 훨씬 위험하다. RON/MON만으로는 초노킹 경향을 못 잡는다는 지적이 여기서 나온다.3

6. 옥탄가와 압축비의 흥정[편집]

연료의 내노킹성을 재는 표준 척도가 **옥탄가**다. 기준 연료는 딱 두 개다.

  • iso-옥테인(2,2,4-트라이메틸펜테인) = 100
  • n-헵테인 = 0

이 둘을 섞은 것이 1차 기준연료(Primary Reference Fuel, PRF)이고, PRF 몇 대 몇과 같은 노킹 거동을 보이는지가 그 연료의 옥탄가다. 가지가 많은 iso-옥테인이 저온 반응성이 낮고 직쇄인 n-헵테인이 높다는, 냉염 문서의 6원환 이성질화 이야기가 그대로 척도의 물리적 근거다.

시험은 가변 압축비 단기통 CFR 엔진에서 하며, 조건이 둘이라 값도 둘이다.

지표규격회전수흡기/혼합기 온도
RON (리서치 옥탄가)ASTM D2699600 rpm약 52 °C
MON (모터 옥탄가)ASTM D2700900 rpm약 149 °C

MON 쪽이 훨씬 가혹하고, 둘의 차이 S=RONMONS = \mathrm{RON} - \mathrm{MON}감도(sensitivity)라 부른다. 감도가 크다는 것은 온도 이력에 따라 반응성이 크게 달라진다는 뜻이고, 곧 그 연료의 NTC 거동이 PRF와 다르다는 신호다. 실제 엔진 조건이 RON·MON 시험 조건 어디에도 정확히 맞지 않기 때문에, 운전 조건을 하나의 계수 KK로 요약해

OI=RONKS\mathrm{OI} = \mathrm{RON} - K\,S

로 유효 옥탄가를 잡는 옥탄 지수(Octane Index) 접근이 널리 쓰인다. 현대 과급 엔진에서는 KK가 음수가 되는 조건도 흔한데, 그러면 감도가 큰 연료(에탄올 블렌드 등)가 오히려 유리해진다.

왜 이걸 신경 쓰는가. 이상 오토 사이클의 열효율은 압축비 rr에 대해

η=11rγ1\eta = 1 - \frac{1}{r^{\gamma-1}}

이다. 압축비를 올리면 효율이 오른다. 그런데 압축비를 올리면 말단가스의 온도·압력 이력도 같이 올라가 τig\tau_{ig}가 짧아진다. 효율과 노킹은 같은 손잡이의 양쪽 끝이고, 그래서 엔진 설계의 압축비는 “쓸 수 있는 연료의 옥탄가”가 정한다. 밀러/앳킨슨 사이클, 물분사, 쿨드 EGR, 가변 압축비가 전부 이 흥정에서 조금이라도 유리한 자리를 사려는 시도다.

7. 어떻게 계산하는가[편집]

노킹 예측의 표준은 놀랍게도 70년 된 한 줄이다. 1955년 리븐굿과 우가 제안한 리븐굿-우 적분.

0tknockdtτig(T(t),p(t))=1\int_0^{t_{\text{knock}}} \frac{dt}{\tau_{ig}\big(T(t),\,p(t)\big)} = 1

말단가스가 겪는 온도·압력 이력을 따라가면서, 매 순간의 조건에 해당하는 τig\tau_{ig}의 역수를 누적한다. 누적량이 1이 되면 착화. **“유도 반응이 쌓이는 양”**이라는 직관을 그대로 수식으로 옮긴 것이고, 근거는 유도기간 동안 라디칼 풀이 조건에 무관하게 비슷한 임계까지 자란다는 가정이다.

실무 파이프라인은 이렇게 굴러간다.

  1. 상세 화학으로 0차원 균질 반응기를 돌려 τig(T,p,ϕ,EGR)\tau_{ig}(T,p,\phi,\mathrm{EGR}) 테이블 또는 상관식을 만든다(자착화).
  2. 엔진 사이클 시뮬레이션 또는 CFD에서 말단가스의 T(t),p(t)T(t), p(t) 이력을 뽑는다.
  3. 위 적분을 누적해 1에 도달하는 크랭크각을 노킹 개시로 본다.

CFD에서는 적분값 자체를 수송되는 스칼라로 취급하는 transported Livengood–Wu 형태로 확장해, 말단가스 덩어리가 이동·혼합하는 효과까지 담는다. 대와류 모사로 사이클 간 변동(cycle-to-cycle variation)까지 보면 “몇 번째 사이클에서 노킹하는가”라는 통계적 질문에도 답할 수 있다 — 노킹은 본질적으로 꼬리 사건이라 평균 사이클 해석으로는 안 잡힌다.

한계도 분명하다. 리븐굿-우는 원래 단일 단계 착화를 염두에 둔 근사라, 냉염이 있는 2단 착화 연료에서는 예측력이 떨어진다. 1단 발열이 이미 말단가스를 데워 놓기 때문에 τig\tau_{ig} 테이블의 입력 조건 자체가 흔들리기 때문. 그래서 저온 화학이 중요한 조건에서는 적분식 대신 말단가스에서 축소 메커니즘을 직접 적분하는 쪽으로 간다. 비싸지만 정직하다.

그리고 리븐굿-우가 답하는 것은 언제 터지는가뿐이다. 얼마나 세게 터지는가는 앞 절의 ξ\xiε\varepsilon, 즉 반응성 구배와 여기 시간의 문제라 별도 모형이 필요하다. 노킹 개시 예측은 꽤 성숙했지만 노킹 강도 예측은 여전히 연구 주제다.4

8. 관련 문서[편집]

9. Footnotes[편집]

  1. 노크센서 캘리브레이션은 엔진 개발에서 가장 지루하면서 가장 안 끝나는 작업 중 하나로 꼽힌다. 임계를 낮추면 멀쩡한 소음에 반응해 점화를 물려 출력이 깎이고, 높이면 진짜 노킹을 놓쳐 엔진이 상한다. 게다가 판정 창·필터 대역·임계가 회전수와 부하마다 다 다르다. 결국 맵이 수천 점짜리 테이블이 되고, 그 테이블을 채우는 것이 사람의 현생이다.

  2. LSPI가 발견됐을 때 업계 반응이 재미있었다. 연소 연구자들은 “화염 전파 문제”로, 윤활 연구자들은 “오일 소모 문제”로 접근했는데, 결과적으로 실효적인 해법의 상당 부분이 오일 처방 쪽에서 나왔다. 엔진 파괴의 원인이 연료가 아니라 오일이었다는 결말은 지금도 회자된다.

  3. ξ\xi가 너무 작아도 폭굉이 안 난다는 것이 처음엔 반직관적이다. 하지만 자착화 전파가 음속보다 압도적으로 빠르면 압력파가 반응면과 보조를 맞출 시간이 없다 — 그냥 전체가 동시에 정적 폭발한다. 압력은 오르되 조직된 폭굉파는 안 생긴다. 폭굉은 화학과 음향이 박자를 맞출 때만 성립하는 현상이라, 반도 모양의 좁은 영역이 나오는 것이 오히려 자연스럽다.

  4. 노킹 예측 논문에서 “실험과 잘 맞는다”의 의미를 확인할 필요가 있다. 대개는 노킹 개시 크랭크각을 몇 도 안에서 맞췄다는 뜻이고, 압력 진동 진폭까지 맞췄다는 뜻인 경우는 드물다. 진폭은 핫스팟의 크기·개수·분포에 지배되는데 그것들은 사이클마다 난수에 가깝다. “언제”는 결정론, “얼마나”는 통계 — 이 분리를 인정하는 것이 이 분야의 성숙도 지표다.