[논문 리뷰] Highly Accurate Solutions of the Blasius and Falkner-Skan Boundary Layer Equations via Convergence Acceleration

B. D. Ganapol|arXiv (Cornell University)|2010. 06. 19.

Iterative Methods for Nonlinear Equations참고 문헌 15인용 수 14

한 줄 요약

이 논문은 페르크너-스칸 및 블라시우스 경계층 방정식을 해결하기 위한 매우 정확하고 강력한 수치 알고리즘을 제시한다. 매클로린 급수 전개와 재귀적으로 유도된 계수, 윈-에프스실 수렴 가속화, 연속적인 해석적 계속을 사용하며, 고정밀도 산술이나 기호 계산을 사용하지 않고 이중 정밀도 포트란에서 10자리 이상의 정확도를 달성한다. 이는 피부 마찰 계수 측면에서 새로운 기준을 설정한다.

ABSTRACT

A new highly accurate algorithm for the solution of the Falkner-Skan equation of boundary layer theory is presented. The algorithm, based on a Maclaurin series representation, finds its coefficients from recurrence. In addition, Wynn-epsilon convergence acceleration and continuous analytical continuation enable an accurate evaluation. The most accurate skin friction coefficients (shooting angle) to date are presented along with comparisons to past and present values found in the literature. The algorithm, coded in FORTRAN, uses neither enhanced precision arithmetic beyond quadruple precision nor computer algebra to achieve results in a timely fashion. Key Words: Falkner-Skan flow; Blasius flow; Wynn-epsilon acceleration; Romberg acceleration; Continuous analytical continuation

연구 동기 및 목표

• 페르크너-스칸 및 블라시우스 경계층 방정식을 해결하기 위한 강력하고 매우 정확한 수치 알고리즘을 개발하는 것.
• 세부 조정이 필요한 초기 추정치와 수치적 불안정성 문제를 야기하는 기존 수치 방법의 한계를 극복하는 것.
• 사용자 정의 정밀도 산술이나 컴퓨터 대수 시스템에 의존하지 않고도 높은 정확도(10자리 이상)를 달성하는 것.
• 기존 문헌 값과의 비교를 통해 피부 마찰 계수를 기준화하고 이전 결과를 초월하는 정확도를 확보하는 것.
• 전진 및 역류 영역 모두에서의 효과성을 입증하며, 특히 음수 β 값에 대한 적용도 포함하는 것.

제안 방법

• 페르크너-스칸 비선형 미분방정식에서 유도된 재귀 관계를 통해 계수를 생성하는 매클로린 급수로 해를 표현한다.
• 유한 수렴 반경을 초월하는 급수의 수렴을 향상시키기 위해 윈-에프스실 알고리즘을 사용한 수렴 가속화를 적용한다.
• 속도 프로파일이 1에 수렴하는 渐近 영역으로 해를 연장하기 위해 연속적인 해석적 계속 기법을 사용한다.
• 무한대에서의 경계 조건을 만족시키기 위해 이분법 기반의 射撃 방법을 사용하여 초기 기울기(사격 각도)를 반복적으로 조정한다.
• 뉴턴-라프슨 반복을 피하고 물리적 통찰력과 이분법을 기반으로 하여 강력성과 수렴성을 확보한다.
• 모든 계산은 표준 이중 정밀도 포트란에서 수행되며, 정확도를 확인하기 위해 일부 부분은 사중 정밀도 검증을 실시한다.

실험 결과

연구 질문

• RQ1 수렴 가속화와 해석적 계속 기법을 적용한 매클로린 급수 기반 접근법이, 고정밀도 산술을 사용하지 않고도 페르크너-스칸 방정식에 대해 매우 정확한 해를 도출할 수 있는가?
• RQ2 기존의 룬게-쿠타 해법에 의존하고 초기 추정치를 세밀하게 조정이 필요한 수치 방법과 비교해 볼 때, 제안된 알고리즘의 정확도와 강력성은 어떠한가?
• RQ3 특히 역류를 나타내는 음수 β 값에 대해, 블라시우스 및 페르크너-스칸 방정식의 가장 정확한 피부 마찰 계수(사격 각도)는 무엇인가?
• RQ4 이 방법은 전진 및 역류 경계층 유동을 신뢰성 있게 처리할 수 있는가? 특히 β ≈ -0.198837735와 같은 도전적인 경우도 포함하여.
• RQ5 수렴 가속화와 해석적 계속 기법이 매클로린 급수의 자연적 수렴 반경을 초월하여 유효한 해 평가 영역을 얼마나 넓힐 수 있는가?

주요 결과

• 이 알고리즘은 고정밀도 산술을 사용하지 않고도 페르크너-스칸 방정식에 대해 10자리 이상의 유의 숫자 정확도를 달성한다.
• 블라시우스의 경우(β = 0) 가장 정확한 피부 마찰 계수는 α = -0.469599988361로 보고되었으며, 새로운 기준 값으로 -0.46959998836100000000으로 확인되었다.
• β = -0.198837735일 경우 피부 마찰 계수는 α = -0.19883773500000000000로 계산되어 유동 반전의 임계값 존재를 확인한다.
• 음수 β 값에 대해, 즉 β = -0.16 및 β = -0.198837735에 대해도 역류 해를 일관되고 정확하게 계산하였다.
• 알고리즘은 나쁜 초기 추정치(예: β = 1000일 때 α = 200)에서도 수렴을 보이며, 높은 신뢰도로 α = 36.5171968을 도출한다.
• 본 연구는 블라시우스의 원래 매클로린 급수 접근법이 현대 수렴 기법과 결합될 경우 여전히 경계층 문제 해결에 강력하고 정확한 방법임을 확인한다.