Supersonic jet impingement on concave surfaces

기술 요약: 오목한 표면에 대한 초음속 제트 충돌

문제 정의
비이상적으로 팽창된 초음속 제트가 표면에 충돌할 때 발생하는 에어로어쿠스틱 공진(스크리치, screech)은 잘 알려진 현상이지만, 벽면 곡률이 이러한 공진에 미치는 영향은 여전히 명확히 이해되지 않고 있다. 자유 제트(free jet)와 평판 충돌(flat-wall impingement)에 대한 피드백 루프 메커니즘은 하류로 이동하는 불안정성 파동과 상류로 전파되는 음향 또는 유도파를 포함하는 것으로 확립되어 있으나, 오목한 표면이 피드백 경로, 톤 진폭 및 모드 선택을 수정하는 구체적인 역할은 불분명하다. 오목한 벽면은 음향 거울 역할을 하여 상류 파동을 집중시키거나, 유효 루프 길이를 변화시키거나, 소스 강도를 변화시킬 가능성이 있다. 본 연구는 가우시안 오목 함몰(Gaussian concave indentation)의 확산 정도가 흐름의 불균일성, 표면 하중 및 스크리치 특성에 미치는 영향을 조사하며, 전역 루프 길이와 상류 파동 분산의 역할을 구분한다.

연구 방법론
본 연구는 이상적 팽창 마하수 $M_j = 1.56$ 및 레이놀즈 수 $6 \times 10^4$ 조건에서 6가지 충돌 형상을 조사하기 위해 압축성 대규모 와류 시뮬레이션(LES)을 채택한다. 기하학적 구조는 노즐-벽 간 거리 $L/D = 2.08$ 및 $2.58$인 두 개의 평판과, 고정된 깊이를 가지되 확산 매개변수($\sigma \in {0.4, 0.8, 1.6, 4.0}$)가 변하는 가우시안 프로파일을 가진 네 개의 오목한 표면을 포함한다.

수치 해석에는 구조화된 멀티 블록 격(약 1억 2,700만 셀) 상에서 3차원 비정상 압축성 나비에-스토크스 방정식을 풀기 위한 고차 저분산 유한 차분 솔버가 사용된다. 이 시뮬레이션은 전단층 발달, 충격 셀 구조, 마하 디스크 역학을 포함한 전체 피드백 루프를 포착한다.

주요 분석 도구는 다음과 같다:

• 고유 직교 분해 (POD): 지배적인 결합 구조 및 진동 모드를 식별하기 위함.
• Powell-Tam 소스-전달 예산 (Source-Transfer Budget): 곡률에 의한 증폭을 마하 디스크 소스 강도($\eta_g$)와 상류 전달/수용 효율($\eta_u \eta_r$)의 기여도로 분할하기 위함.
• 와류 시트 모델링 (Vortex-Sheet Modeling): 유도 제트 모드(GJM) 분산 관계를 예측하고 이를 LES 데이터와 비교하기 위함.
• 파수-주파수 분석 (Wavenumber-Frequency Analysis): 내부(플룸 제한) 및 외부(주변부) 상류 전파 파동을 구분하고 이들의 분산 특성을 결정하기 위함.

주요 결과

1. 곡률 의존적 톤 증폭:
   함몰 확산(indentation spread, $\sigma$)이 좁아질수록(즉, 오목한 벽이 제트 축 주위에 더 국부적으로 위치할수록), 일차 스크리치 톤 진폭이 크게 증가한다. 가장 좁은 함몰($\sigma=0.4$)은 $L/D=2.6$인 평판 기준 대비 최대 23 dB의 일차 톤 증폭을 보인다. 이 증폭은 두 가지 상보적인 메커니즘에 기인한다:

   • 소스 증폭: 제한된 재순환 영역 내에서의 강력한 충격-와류 상호작용으로 인한 마하 디스크의 진폭 증가.
   • 상류 전달 효율: 상류 파동의 노즐로의 회귀 강화. Powell-Tam 예산에 따르면, 가장 좁은 함몰의 경우 소스가 전체 증폭의 약 62%, 상류 전달이 약 38%를 기여한다.
   • 집중 메커니즘: 상류로 전파되는 파동은 주변부 음향 파동보다는 기둥 내에 제한된 성분(제트 플룸 내부)에 의해 지배되는 것으로 나타났다. 함몰이 제트 기둥 아래에 국한될 때, 주변 채널보다 기둥 채널이 훨씬 더 큰 진폭(최대 15.8 dB)을 얻는 것은 오목한 벽이 제트 기둥을 통해 회귀하는 파동에 주로 작용함을 시사한다.

2. 모드 선택 및 대칭성:
   본 연구는 함몰 확산에 따른 진동 모드의 뚜렷한 분기(bifurcation)를 식-별한다:

   • 축대칭 모드 (Mode A): 가장 좁은 함몰($\sigma=0.4$)과 $L/D=2.08$인 평판에서 관찰된다. 이 경우들은 마하 디스크의 강력하고 주기적인 축 방향 맥동과 토로이드형 제트 기둥 진동을 보인다.
   • 나선형 모드 (Mode C): 더 넓은 함몰($\sigma \ge 0.8$)과 $L/D=2.58$인 평판에서 관찰된다. 이 경우들은 약하고 광대역적인 마하 디스크 운동과 회전하는 나선형 제트 기둥 진동을 보인다.
   • 축대칭 모드에서 나선형 모드로의 전이는 함몰이 제트 기둥 반경 너머로 확산되는 시점과 일치한다.

3. 표면 하중 특성:
   스크리치 모드의 대칭성은 비정상 표면 하중의 성격을 결정한다:

   • 축대칭 스크리치: 톤 성분의 법선(축) 방향 힘을 생성하며, 횡방향 힘이나 굽힘 모멘트는 무시할 만한 수준이다.
   • 나선형 스크리치: 광대역 법선 힘을 생성하지만, 톤 에너지는 회전하는 면내(in-plane) 굽히 모멘트에 집중된다. 이 모멘트는 제트 축을 중심으로 타원 궤도를 따라 세차 운동을 하며, 이는 축 방향 하중을 갖는 축대칭 모드와 구별되는 특정 피로 위험을 야기한다. 최대 굽힘 모멘트는 반드시 가장 강한 곡률에서 발생하는 것이 아니라, 중간 정도의 함몰 확산에서 발생한다.

4. 주파수 선택 메커니즘:
   본 연구는 스크리치 주파수를 결정하는 두 가지 뚜렷한 메커니즘을 밝혀냈다:

   • 나선형 모드: 나선형 사례(오목한 경우 및 평판 모두)의 주파수는 와류 시트 모델에 의해 예측된 $H_1$ 유도 제트 모드(GJM) 분지(branch)의 저주파 한계와 밀접하게 일치한다. 이는 와류 시트 모델의 전단층 프로파일에 의해 결정되며, 벽면 기하학이나 루프 길이에 거의 영향을 받지 않는다.
   • 축대칭 모드: 축대칭 사례의 주파수는 Powell의 고전적인 루프 길이 기준($N/f = \ell/u_c + \ell/a_\infty$)을 따르며, 와류 시트 모델의 중립 $A_n$ 분지와 일치하지 않는다. 이는 축대칭 스크리치의 경우 루프 길이가 주요 선택 기준이며, 유도 모드 제약 조건이 적용되지 않거나 충족되지 않음을 시사한다.

의의 및 주장
본 논문은 벽면 곡률이 스크리치 진폭과 표면 하중을 제어하는 효과적인 메커니즘을 제공하며, 가장 큰 증폭은 오목한 함몰이 제트 기둥에 반경 방향으로 국한될 때 발생함을 입증한다고 주장한다.

결정적으로, 본 연구는 충돌하는 제트의 나선형 스크리치와 축대칭 스크리치 사이의 물리적 메커니즘을 구분한다. 나선형 톤은 유도 제트 모드(특히 $H_1$ 분지)의 분산 특성에 의해 지배되어 벽면 기하학 변화에 견고한 반면, 축대칭 톤은 전역 피드백 루프 길이에 의해 지배됨을 확립하였다. 또한, 상류 피드백 경로가 이론적으로 음향 거울 역할을 할 수 있는 오목한 벽이 존재함에도 불구하고, 주로 제트 플룸 내부(기둥 제한)를 통해 전달된다는 점을 명확히 하였다.

이 연구 결과는 충돌하는 제트의 피드백 루프에 대한 정교한 이해를 제공하며, 소스 강도(마하 디스크 역학)와 상류 전달 효율을 분리하고, 선택된 진동 모드에 기초하여 피로 관련 하중 패턴(축 방향 vs 회전 굽힘)을 예측할 수 있는 근거를 제공한다.