[논문 리뷰] A numerical stability investigation of strong ZND detonations for Majda's model

Blake Barker, Kevin Zumbrun|arXiv (Cornell University)|2010. 11. 06.

Combustion and Detonation Processes참고 문헌 26인용 수 7

한 줄 요약

이 논문은 다양한 활성화 에너지에서 아레니우스형 ignitions 함수를 가진 Majda의 모델에서 강한 ZND 폭발의 수치적 안정성을 조사한다. 고급 고유값 해소기법을 사용하여 고활성화 에너지 한계에서도 진동 불안정성이 관찰되지 않았으며, Humpherys–Zumbrun 알고리즘이 고전적 Lee–Stewart 방법보다 100~1,000배 빠르며, 적응형 메쉬가 가장 큰 성능 향상을 제공함을 입증한다.

ABSTRACT

We carry out a systematic numerical stability analysis of ZND detonations of Majda's model with Arrhenius-type ignition function, a simplified model for reacting flow, as heat release and activation energy are varied. Our purpose is, first, to answer a question of Majda whether oscillatory instabilities can occur for high activation energies as in the full reacting Euler equations, and, second, to test the efficiency of various versions of a numerical eigenvalue-finding scheme suggested by Humpherys and Zumbrun against the standard method of Lee and Stewart. Our results suggest that instabilities do not occur for Majda's model with Arrhenius-type ignition function, nor with a modified Arrhenius-type ignition function suggested by Lyng--Zumbrun, even in the high-activation energy limit. We find that the algorithm of Humpherys--Zumbrun is in the context of Majda's model $100$-$1,000$ times faster than the one described in the classical work of Lee and Stewart and $1$-$10$ times faster than an optimized version of the Lee--Stewart algorithm using an adaptive-mesh ODE solver

연구 동기 및 목표

• 고활성화 에너지 한계에서 아레니우스형 점화 함수를 가진 Majda의 모델에서 허프비프 분기와 같은 진동 불안정성이 발생하는지 여부를 규명하는 것.
• 폭발 안정성 분석의 맥락에서 Humpherys–Zumbrun 고유값 찾기 기법이 고전적 Lee–Stewart 방법보다 효율성이 높은지 테스트하는 것.
• 메쉬 조절, 좌표 매개변수화, 기법 설정과 같은 수치적 구현 선택 사항이 계산 성능 및 수렴성에 미치는 영향을 평가하는 것.
• 불안정한 고유값에 대한 엄밀한 경계를 설정하고, 정사각형파 근사와 같은 특이한 극한 상황을 다룰 수 있는 하이브리드 기법을 개발하는 것.
• Majda의 모델을 단순화된 시험기반으로 삼아 더 복잡한 반응성 Euler 시스템에 적용 가능한 수치적 기법의 기준이 되는 것.

제안 방법

• 아레니우스 및 수정 아레니우스 점화 함수를 가진 Majda의 모델에서 ZND 폭발 프로파일에 대한 체계적인 수치적 고유값 분석을 수행하였다.
• Humpherys–Zumbrun Evans 함수 유형 기법을 재구성하여 x=0에서 x=-∞로의 후진 사 Shoot 방법으로 구현하였으며, 스펙트럼 갭의 이점을 활용하였다.
• 공간(x) 및 주파수(λ) 변수에서 적응형 단계 ODE 해석기를 적용하여 수렴성 향상과 계산 비용 절감을 도모하였다.
• z좌표 매개변수화(프로파일이 z에서 명시적으로 해석 가능)와 x좌표 매개변수화(기존 설정)를 적용하여 성능 영향을 비교하였다.
• 고차수 고주파수 점근 해석과 하이브리드 기법을 개발하여 특이한 극한 상황, 특히 정사각형파 근사를 다룰 수 있도록 하였다.
• 다양한 기법 간 광범위한 수치적 비교를 수행하였으며, 이에는 μ(x), 수반 μ(x), 하이브리드, Lee–Stewart(적응형 및 고정 메쉬), 극좌표 수반, 극좌표 수반 반경 기법 등이 포함되었다.

실험 결과

연구 질문

• RQ1 Majda의 모델에서 아레니우스형 점화 함수를 가진 경우, 전체 반응성 Euler 방정식에서 관찰된 것처럼 고활성화 에너지 한계에서 진동 불안정성이 나타나는가?
• RQ2 Humpherys–Zumbrun 고유값 찾기 기법이 ZND 폭발 안정성 분석의 계산 효율성 측면에서 고전적 Lee–Stewart 방법을 능가할 수 있는가?
• RQ3 적응형 메쉬 대비 고정 메쉬, x좌표 대비 z좌표 매개변수화, 기법 설정 방식 등의 다양한 수치적 선택 사항이 계산 속도와 수렴성에 어떤 영향을 미치는가?
• RQ4 Majda의 모델에서 불안정한 고유값의 최대 크기는 얼마이며, 이를 위한 엄밀한 경계를 도출할 수 있는가?
• RQ5 하이브리드 또는 극한 기법이 점화 함수의 특이한 정사각형파 근사 상황을 안정성 계산에서 효과적으로 다룰 수 있는가?

주요 결과

• 고활성화 에너지 한계에서도 아레니우스형 또는 수정 아레니우스형 점화 함수를 가진 Majda의 모델에서 진동 불안정성 또는 기타 불안정성은 관찰되지 않았다.
• Humpherys–Zumbrun 기법은 표준 Lee–Stewart 방법보다 100배에서 1,000배 빠르며, 최적화된 Lee–Stewart 변종 대비 1~10배의 성능 향상을 보였다.
• x진화에서 적응형 단계 ODE 해석기를 사용함으로써 계산 시간이 고정 메쉬 해석기 대비 100배에서 1,000배 감소하였으며, 이는 가장 큰 성능 향상 요소였다.
• z좌표 매개변수화 방법은 x좌표 매개변수화 대비 일관된 요인 6의 성능 향상을 보였으며, 주로 프로파일 평가 비용 감소로 인한 것이다.
• 최적의 수치적 선택을 통해 안정성 계산은 각 파rameter 세트당 몇 초 내로 완료되었으며, 이는 단순화된 모델의 양호한 수치 조건성을 시사한다.
• 고활성화 에너지 한계에서 성능이 저하되었지만, 성장/감쇠 행동의 분해를 통해 이전에 효과적이지 못했던 기법이 실패한 곳에서도 수렴이 가능해졌다.