[논문 리뷰] Exact solution of a differential problem in analytical fluid dynamics: use of Airy's functions

Gianluca Argentini|ArXiv.org|2006. 06. 28.

Fluid Dynamics and Turbulent Flows참고 문헌 7인용 수 4

한 줄 요약

이 논문은 유선을 따라 유체역학적 방정식을 리만-타입의 상미분방정식으로 단순화함으로써 2차원 정 steady 압축성 나비에-스토크스 경계값 문제에 대한 정확한 해석적 해를 제시한다. 공기역학적 함수를 활용한 변환을 통해 해는 Ai(t)와 Bi(t) 및 그 도함수들의 유리적 조합으로 표현되며, 일정한 운동점성계수와 압력-밀도 상관관계에 대한 특정 구성관계를 가정할 경우 닫힌 형태의 속도 분포를 제공한다.

ABSTRACT

Treating a boundary value problem in analytical fluid dynamics, translation of 2D steady Navier-Stokes equations to ordinary differential form leads to a second order equation of Riccati type. In the case of a compressible fluid with constant kinematic viscosity along streamlines, it is possible to find an exact solution of the differential problem by rational combination of Airy's functions and their derivatives.

연구 동기 및 목표

• 유선을 따라 2차원 정 steady 압축성 나비에-스토크스 경계값 문제에 대한 정확한 해석적 해를 유도하기.
• 유선 매개변수화와 속도 분해를 통해 나비에-스토크스 방정식을 상미분방정식으로 단순화하기.
• 일정한 운동점성계수와 압력-밀도 상관관계에 대한 특정 구성관계를 가정할 경우 리카티 유형의 방정식을 에어리 함수를 활용해 정확히 해결하기.
• 초기 및 최종 속도 조건을 통해 결정되는 세 개의 적분상수를 활용한 경계값 해 설정하기.

제안 방법

• 유선을 스칼라 매개변수 s로 매개변수화하고, v₁ ≠ 0 이며 φ₁(s) = s 라는 단순화 조건을 가정한다.
• 사슬법칙과 속도 분해 u₁ = v₁(φ(s))를 활용해 나비에-스토크스 방정식을 유선을 따라 상미분방정식으로 변환한다.
• 비선형 2차 상미분방정식을 선형 2차 상미분방정식(에어리 유형)으로 변환하기 위해 리카티 변환 u₁ = -2ν (ż/z)를 적용한다.
• t가 점성계수 및 압력-밀도 매개변수를 포함하는 s의 선형 함수로 정의되는 조건에서, 변환된 방정식을 에어리 함수 Ai(t)와 Bi(t)를 사용해 해석한다.
• 물리적 타당성이 확보되는 조건 a < 0 하에서 속도 해를 u₁ = 2ν (-a)^{1/3} (c₁Ai’(t) + c₂Bi’(t)) / (c₁Ai(t) + c₂Bi(t)) 형태로 재구성한다.
• 경계 조건 u₁(0) = u₁₀, ẋ₁(0) = ẋ₁₀, u₁(L) = u₁L 를 통해 상수 c₁, c₂, c 를 결정한다.

실험 결과

연구 질문

• RQ1 특수함수를 활용해 유선을 따라 2차원 정 steady 압축성 나비에-스토크스 방정식에 대한 정확한 해석적 해를 도출할 수 있는가?
• RQ2 운동점성계수와 압력 기울기 대비 밀도 비율이 유선을 따라 일정할 경우 속도 분포의 형태는 어떻게 되는가?
• RQ3 나비에-스토크스 축소 과정에서 유도된 리카티 유형의 상미분방정식을 어떻게 에어리 함수를 활용해 정확히 해결할 수 있는가?
• RQ4 어떤 물리적 조건에서 유도된 해가 실수이면서 물리적으로 의미 있는가?
• RQ5 세 점 경계값 문제를 만족시키기 위해 적분상수는 어떻게 결정되는가?

주요 결과

• 유선을 따라 나비에-스토크스 방정식의 정확한 해는 에어리 함수 및 그 도함수들의 유리적 조합 형태로 표현된다.
• 매개변수 a = (1/(2ν²))((q̇/ρ) - f₁) 가 음수일 경우 해가 유효하며, 이는 실제적인 가정 조건인 압력 감소(q̇ < 0)와 양의 외력(f₁ > 0) 조건에서 성립한다.
• 속도 성분 u₁은 유선을 따라 u₂ 성분에 의해 결정되며, 이는 dφ₂/ds 의 도함수를 통해 u₁을 스케일링한다.
• 세 개의 조건(u₁(0), ẋ₁(0), u₁(L))을 갖는 경계값 문제는 대수방정식의 연립방정식을 통해 세 상수 c, c₁, c₂ 를 유일하게 결정한다.
• 유선이 직선이고 압축성이 없을 경우 해는 일정한 속도장으로 축소되며, 고전 유체역학과 일치한다.
• 이 방법은 일정한 운동점성계수와 특정한 압력-밀도 구성관계를 갖는 일군의 정 steady, 2차원 압축성 유동에 대해 닫힌 형태의 해석적 해를 제공한다.