본 연구는 액체로켓엔진에서 나타나는 비선형 연소불안정을 해석하기 위한 수치해석 기법을 개하는데 목적이 있다. 액체로켓엔진에서 연소불안정의 비신형적 거동은 한계사이클의 존재와 비선형 촉발 불안정으로 특징지어진다. 예비해석 단계로서 로켓노즐의 형상을 고려한 준1차원 방정식에 대한 해석기법 개발연구가 진행되었고 대류항에 대하여 근사리만해법이 적용되었다. 공진관내에서 일어나는 비선형 음향거동의 결과는 연소볼안정에서 흔히 수반되는 N-파형 충격파의 발생메커니즘을 잘 묘사할 수 있음을 보여주었다. 우주왕복선의 주엔진에 대한 연소불안정이 n-τ 연소응답모텔에 의해서 조시되었다. 초기에 주어진 압력교란은 천이구간을 거쳐 스스로 주어진 자동조건에 맞는 연소불안정 모드로 변화되는 것을 발견하였으며, 초기 교란의 특성(크기, 모드)에 상관없이 동일한 한계사이클에 도달하는 것을 알 수 있었다. 이를 통해 연소파라미터에 따른 선형안정한계를 구할 수 있었으나, 선형화된 n-τ 연소응답모델의 한계로 인해 비선형안정한계는 예측할 수 없었다.