[논문 리뷰] On the existence of self-similar converging shocks for arbitrary equation of state

Zachary M. Boyd, Scott D. Ramsey외 1명|arXiv (Cornell University)|2017. 07. 12.

Computational Fluid Dynamics and Aerodynamics참고 문헌 29인용 수 7

한 줄 요약

이 논문은 임의의 상태방정식에서 자가유사 수축 충격 해가 존재하기 위한 필수 열역학적 조건—구체적으로, 단열 부피모odulus $ K_S = p f(\rho) $ 형태여야 한다—를 규명한다. 충격파가 자가유사해지기 위한 조건을 분석하고, Guderley 문제의 해가 비선형 고유값 문제로 축소됨을 보이며, 일정한 충격 속도($ \lambda = 0 $)는 유계해를 유도하지 못해 비이상적 물질에 대한 표준 Guderley 프레임워크를 무효화함을 밝힌다.

ABSTRACT

We extend Guderley's problem of finding a self-similar scaling solution for a converging cylindrical or spherical shock wave from the ideal gas case to the case of flows with an arbitrary equation of state closure model, giving necessary conditions for the existence of a solution. The necessary condition is a thermodynamic one, namely that the adiabatic bulk modulus, $K_S$, of the fluid be of the form $pf(ρ)$ where $p$ is pressure, $ρ$ is mass density, and $f$ is an arbitrary function. Although this condition has appeared in the literature before, here we give a more rigorous and extensive treatment. Of particular interest is our novel analysis of the governing ordinary differential equations (ODEs), which shows that, in general, the Guderley problem is always an eigenvalue problem. The need for an eigenvalue arises from basic shock stability principles -- an interesting connection to the existing literature on the relationship between self-similarity of the second kind and stability. We also investigate a special case, usually neglected by previous authors, where assuming constant shock velocity yields a reduction to ODEs for every material, but those ODEs never have a bounded, differentiable solution. This theoretical work is motivated by the need for more realistic test problems in the verification of inviscid compressible flow codes that simulate flows in a variety of non-ideal gas materials.

연구 동기 및 목표

• 이deal 기체에서의 Guderley 자가유사 수축 충격 문제를 임의의 상태방정식으로 확장하여, 비이상적 물질에 대한 보다 현실적인 코드 검증을 가능하게 한다.
• 압축성 유동에서 유계이고 미분 가능한 자가유사 해가 존재하기 위한 필수 열역학적 조건—특히 $ K_S = p f(\rho) $—를 규명한다.
• 유도된 상미분방정식(OED) 시스템을 비선형 고유값 문제로 엄밀히 분석하여, 자가유사 해에서 안정성과 척도의 역할을 명확히 한다.
• 일정한 충격 속도($ \lambda = 0 $)의 특수한 경우를 조사하여, 비록 상미분방정식이 단순화되지만, 유계해를 도출하지 못함을 보인다.
• 비이상적 압축성 유동에서 수치해석 코드 검증을 위한 정확한 해석적 시험 문제를 구성하는 기초를 제공한다.

제안 방법

• 임의의 상태방정식에서의 압축성 유동에 대한 유동의 에이러 방정식을 자가유사화하여, 반경 및 시간 척도 기반의 유사성 변수를 사용한다.
• 자기유사 해 존재를 위해 핵심 열역학량으로 단열 부피모odulus $ K_S $ 를 규명하고, $ K_S = p f(\rho) $ 가 필수적임을 보인다.
• 지배하는 편미분방정식(PDE)을 유사성 변수 $ \xi = r / |t|^\lambda $ 에 대한 상미분방정식(OED) 시스템으로 축소하며, $ \lambda $ 는 척도 행동을 결정한다.
• 유계이고 미분 가능한 해를 위해 OED 시스템을 분석하여, 충격 안정성 제약 조건으로 인해 문제 자체가 본질적으로 비선형 고유값 문제임을 밝힌다.
• 충격파 전면에서 OED의 특이성 분석을 수행하여, 고유값 문제로 해결하지 않으면 유계해가 발생하지 않음을 보였다.
• 일정한 충격 속도($ \lambda = 0 $)의 경우를 고려하여, 비록 OED가 단순화되지만 분모가 0이 되고 분자가 0이 아니므로, 어떤 유계이고 미분 가능한 해도 존재하지 않음을 입증한다.

실험 결과

연구 질문

• RQ1 임의의 유체에서 자가유사 수축 충격 해가 존재하기 위해 상태방정식에 필요한 열역학적 조건은 무엇인가?
• RQ2 압축성 유동의 에이러 방정식이 일반 상태방정식에서 자가유사화될 때의 행동은 어떻게 되며, 충격 안정성으로부터 어떤 제약 조건이 유도되는가?
• RQ3 비이상적 물질에 대한 Guderley 문제의 해가 비선형 고유값 문제로 축소되는 이유는 무엇이며, 이에 대한 물리적 근거는 무엇인가?
• RQ4 Guderley 프레임워크에서 $ \lambda = 0 $ (일정한 충격 속도)의 경우가 유계이고 미분 가능한 해를 도출할 수 있는가? 만약 불가능하면 그 이유는 무엇인가?
• RQ5 일정한 충격 속도의 경우가 실패할 경우, 비이상적 물질의 수치해석 코드 검증을 위한 해석적 시험 문제 구성에 어떤 영향을 미치는가?

주요 결과

• 임의의 상태방정식에서 자가유사 수축 충격 해가 존재하기 위한 필수 및 충분한 열역학적 조건은 $ K_S = p f(\rho) $ 이며, 여기서 $ f $ 는 밀도의 임의의 함수이다.
• 임의의 상태방정식에 대해 Guderley 문제의 해는 항상 충격 안정성과 해의 유계성 간의 상호작용로 인해 비선형 고유값 문제로 축소된다.
• 표준 자가유사 척도($ \lambda \neq 0 $)는 고유값 조건이 충족될 경우에만 잘 정의된 상미분방정식 시스템이 되며, 이는 유계이고 미분 가능한 해를 보장한다.
• $ \lambda = 0 $ 경우(일정한 충격 속도)는 상미분방정식 시스템이 단순화되지만, 분모가 0이 되고 분자가 0이 아니므로, 어떤 유계이고 미분 가능한 해도 존재하지 않는다.
• $ \lambda = 0 $ 경우는 충격 속도가 일정할 수 없고, 이는 유계성 조건을 위반하므로 유효한 Guderley 해를 도출하지 못한다. 비록 상미분방정식이 단순화되지만 이는 해의 존재를 보장하지 못한다.
• 이 분석은 비이성적 물질, 특히 복잡하거나 비이성적 상태방정식을 가진 물질에 대해 수치해석 코드 검증을 위한 정확한 해석적 시험 문제를 구성하는 데 엄밀한 기초를 제공한다.