On the impact of the turbulent grazing flow development on the acoustic response of an acoustic liner

본 연구는 격자 볼츠만 매우 큰 와류 시뮬레이션(Lattice-Boltzmann Very-Large-Eddy simulations)을 활용하여, 음향 라이너 상부의 난류 흐름의 공간적 발달이 경계층 역학 및 오리피스 흐름 거동을 크게 변화시키며, 이로 인해 위치에 따른 음향 에너지 소산과 현재의 방식으로는 완전히 포착할 수 없는 임피던스 측정값의 차이가 발생함을 입증한다.

핵심

개요: 제트 엔진의 소음을 줄이는 법

제트 엔진을 거대하고 시끄러운 진공청소기라고 상상해 보세요. 엔진의 비명 소리를 멈추기 위해 엔지니어들은 엔진 내부의 공기 통로에 **라이너(liner)**라고 불리는 특수한 "음향 스펀지"를 덧댑니다. 이 라이너는 작은 구멍들이 작은 포켓(공동)으로 이어지는 벌집 구조와 같습니다. 음파가 라이너에 부딪히면 공기가 이 구멍들 사이로 드나들며 마찰과 작은 소용돌이를 만들어내고, 이 과정에서 소리 에너지를 열로 변환하여 엔진의 소음을 효과적으로 줄여줍니다.

하지만 실제 엔진 내부의 공기는 가만히 멈춰 있는 것이 아니라, 매우 빠른 속도로 라이너를 스쳐 지나갑니다(마치 플루트 위로 강한 바람이 불고 있는 것과 같습니다). 이 논문은 음파, 난류(불규칙한 바람), 그리고 이 음향 스펀지가 결합되었을 때 어떤 현상이 발생하는지를 조사합니다.

실험: 디지털 풍동

연구진은 물리적인 엔진을 직접 제작하지 않았습니다. 대신, 대학 연구실의 환경을 재현하기 위해 초정밀 컴퓨터 시뮬레이션("디지털 풍동")을 사용했습니다. 그들은 11줄의 벌집형 공동(cavity)이 있는 덕트 구간을 모델링하고, 바람이 부는 가운데 음파를 쏘아 실험했습니다.

그들은 다음과 같은 다양한 시나리오를 테스트했습니다:

• 풍속: 공기가 얼마나 빠르게 움직이는가.
• 음량: 소리가 얼마나 큰가 (속삭임부터 제트 엔진의 굉음까지).
• 소리의 방향: 소리가 바람을 타고 흐르는가, 아니면 바람을 거슬러 가는가?

주요 결과: "움직이는 카펫" 효과

1. 바람이 공기를 밀어낸다

라이너 표면 바로 옆의 공기를 얇고 끈적한 카펫이라고 생각해 보세요. 바람이 라이너 위를 지나갈 때, 공기는 단순히 매끄럽게 미끄러지는 것이 아니라 라이너의 구멍들이 마치 작은 팬처럼 작동합니다. 이 구멍들이 공기를 표면에서 약간 멀어지게 밀어냅니다.

• 비유: 길가에 서 있는 한 줄의 사람들(구멍)을 상상해 보세요. 강한 바람이 불면 사람들은 뒤로 몸을 기울일 수 있습니다. 만약 사람들이 (소리 때문에) 위아래로 뛰기 시작한다면, 그들은 바람을 더욱 멀리 밀어낼 것입니다.
• 결과: 이로 인해 바람이 흘러가야 할 "더 두꺼운" 공기 층이 형성됩니다. 바람이 구멍들을 따라 이동함에 따라, 이 "공기 카펫"은 점점 더 두꺼워집니다.

2. 하류로 갈수록 바람은 "게을러진다"

공기 카펫이 구멍을 따라 이동할수록 점점 더 두꺼워지기 때문에, 구멍 바로 옆의 풍속은 느려집니다.

• 비유: 바위가 깔린 강물을 상상해 보세요. 처음에는 물살이 빠르고 요동칩니다. 하지만 물이 바위들을 지나면서 아래로 내려갈수록, 바닥 근처의 물은 느려지고 에너지가 줄어듭니다.
• 결과: 라이너의 끝부분에서는 시작 부분에 비해 "전단(shear, 빠른 바람과 느린 공기 사이의 마찰)"이 약해집니다.

3. 소리의 방향에 따라 다르게 반응한다

이 부분이 가장 놀라운 발견입니다. 연구진은 소리가 바람에 대해 어느 방향으로 이동하는지가 매우 중요하다는 것을 발견했습니다.

• 바람을 거슬러 갈 때: 소리가 바람을 거슬러 이동하면, 공기 카펫이 두껍고 바람이 느린 라이너의 "게으른" 끝부분에 먼저 도달합니다. 그 후 "빠른" 끝부분을 향해 이동합니다.
• 바로를 타고 갈 때: 소리가 바람을 타고 이동하면, "빠른" 끝부분에 먼저 도달한 뒤 "게으른" 끝부분을 향해 이동합니다.
• 결과: 바람의 조건이 라이너를 따라 변하기 때문에, 소리 파동은 방향에 따라 서로 다른 "지형"을 경험하게 됩니다. 연구진은 이 두 가지 시나리오에서 라이너의 소리 흡수율이 다르다는 것을 밝혀냈습니다. 이는 마치 언덕을 올라갈 때와 내려갈 때의 차이와 같습니다. 언덕 자체는 같더라도, 당신의 노력과 경험은 달라지는 것과 같습니다.

4. "두 개의 서로 다른 자" 문제

엔지니어들은 보통 라이너가 얼마나 잘 작동하는지 측정할 때 "임피던스(impedance)"라고 불리는 단일 수치(소리에 대한 저항 정도)를 계산합니다.

• 문제점: 이 논문은 만약 라이너의 시작점에서 이 수치를 측정하는 것과 끝점에서 측정하는 것이 서로 다른 결과를 낳는다는 것을 보여줍니다.
• 비유: 방의 "평균 온도"를 측정하려고 하는데, 한쪽은 얼음처럼 차갑고 다른 한쪽은 끓는 것처럼 뜨겁다고 상상해 보세요. 만약 방 전체가 균일하다고 가정하는 자를 사용한다면, 잘못된 답을 얻게 될 것입니다.
• 발견: 컴퓨터 시뮬레이션 결과, 임피던스는 라이너 전체에 걸쳐 고정된 하나의 숫자가 아니었습니다. 바람과 공기 층이 변함에 따라 표면을 따라 이동하며 계속해서 변화합니다.

이 연구가 중요한 이유 (논문에 따르면)

이 논문은 이러한 라이너를 테스트하고 설계하는 현재의 방식들이 종종 바람이 균일하고 공기 층이 얇으며 변하지 않는다고 가정한다는 점을 지적합니다. 이 연구는 그 가정이 틀렸음을 증명합니다.

• 바람이 핵심이다: 바람이 어떻게 발달하는지(공기 층이 두꺼워지고 느려지는 과정)가 소리 흡수 방식에 영향을 미칩니다.
• 방향이 중요하다: 소리가 이동하는 방향이 바람과의 상호작용을 결정합니다.
• 결론: 더 조용하고 성능 좋은 엔진을 설계하기 위해서, 엔지니어들은 라이너를 정적인 물체로 취급하는 것을 멈춰야 합니다. 대신 바람과 공기 층이 표면을 가로질러 이동함에 따라 끊임없이 변화한다는 사실을 반드시 고려해야 합니다.

요약하자면, 음향 라이너는 단순히 정적인 스펀지가 아닙니다. 그것은 바람, 소리, 그리고 공기 층이 함께 어우러져 춤을 추는 역동적인 시스템이며, 그 춤의 방향에 따라 음악이 달라집니다.

기술 요약

기술 요약: 난류 유동 발달이 음향 라이너의 음향 응답에 미치는 영향에 대하여

문제 정의
음향 라이너(일반적으로 단일 자유도(SDOF) 구성)는 항공기 엔진의 소음 저감에 매우 중요하다. 순수 음향 여기(excitation) 또는 흐르는 유동(grazing flow) 하에서의 라이너 물리 현상은 부분적으로 이해되어 있으나, 다중 공동(multi-cavity) 라이너 상에서 발달하는 난류 유동과 음파 사이의 상호작용에 대해서는 정량적 공백이 존재한다. 기존의 임피던스 추출 기법은 종종 공간적으로 균질한 임피던스와 단순화된 경계층 프로파일(예: 무한히 얇은 와류 시트)을 가정한다. 최근 문헌들은 임피던스가 유동 프로파일과 음파 전파 방향에 민감하다는 점을 시사하고 있다. 본 연구는 난류 경계층의 횡방향(streamwise) 발달이 라이너에 의해 수정되고 음향 강제력에 의해 영향을 받는 방식이 국부 및 전역 음향 응답에 어떻게 영향을 미치는지에 대한 정량적 분석의 부재를 다룬다.

연구 방법론
저자들은 상용 솔버인 PowerFLOW를 이용한 고충실도 격자 볼츠만 매우 큰 와류 시뮬레이션(LB-VLES)을 채택하였다. 계산 도메인은 브라질 산타 카타리나 연방 대학교(UFSC)에서 수행된 유동 통과 임피던스 튜브(GFIT) 실험을 재현하였으며, 여기에는 8개의 구멍(orifice)을 포함한 11개의 횡방향 정렬 공동을 가진 라이너가 포함되었다.

• 유동 솔버: 시뮬레이션은 D3Q19 격자 스킴과 VLES 접근 방식을 사용하였으며, 하위 격자 규모는 Renormalisation Group (RNG) $k-\epsilon$ 모델을 통해 모델링하고 대규모 난류 스케일은 직접 해상하였다. 압력 구배를 고려하기 위해 확장된 난류 벽 모델(PGE-WM)을 적용하였다.
• 테스트 케이스: 다음을 포함하는 20개의 시뮬레이션을 수행하였다:
  • 매끄러운 벽 vs. 라이너가 있는 벽 구성.
  • 중심선 마하수 0.0 (유동 없음) 및 0.32 (유동 존재).
  • 음향 주파수 800, 1400, 2000 Hz.
  • 음압 레벨(SPL) 130 dB 및 145 dB.
  • 라이너의 상류(co-flow) 및 하류(counter-flow)에 위치한 음향원.
• 임피던스 추출: 두 가지 방법이 비교되었다:
  1. 모드 매칭(Mode Matching, MM): 평균 임피던스를 결정하기 위해 전역 음향장을 피팅하는 역방법(inverse method).
  2. 인시투(In-situ, Dean의 방법): 페이스 시트(face sheet)와 공동 배면(backplate)의 압력 측정에 기반한 지점별 계산.
• 유동 분석: 공동 내부 및 주변의 유동 역학을 분석하기 위해 평균(mean), 코히어런트(음향 유도), 그리고 확률적(난류) 속도 성분을 분리하는 삼중 분해(triple decomposition) 기법을 사용하였다.

주요 결과

1. 임피던스 가변성 및 방향성:

   • 공간적 불균질성: 인시투 방법은 라이너 표면 전체에서 강한 임피던스 공간적 변화를 드러냈다. 저항(resistance) 값은 구멍과의 상대적 위치에 따라 샘플링 위치에 따라 최대 3배까지 차이가 났다.
   • 추출 방법 간의 불일치: 인시투 방법을 통해 얻은 평균 저항 값은 상류 및 하류 파동 전파 사이에 최소한의 차이를 보였다. 반면, 모드 매칭 방법은 상류로 전파되는 파동에 대해 저항이 더 낮다는 점에서 상당한 차이를 나타냈다. 이는 선택한 추출 기법이 인지된 방향성에 중요한 영향을 미친다는 점을 강조한다.
   • 유동의 영향: 유동은 임피 임피던스의 저항 성분을 증가시켰다. 리액턴스(reactance)는 상대적으로 작은 변화를 보였으며, 주로 공동 깊이와 끝단 보정(end correction)의 영향을 받았다.

2. 경계층 발달 및 변위 두께($\delta^*$):

   • 라이너와 그 구멍의 존재는 유동을 페이스 시트로부터 멀어지게 하여, 횡방향을 따라 경계층 변위 두께($\delta^*$)를 크게 증가시켰다(매끄러운 벽 대비 최대 30% 증가).
   • $\delta^*$는 각 구멍의 하류에서 국부적인 "험프(hump)"를 나타냈다.
   • $\delta^*$의 성장은 음향 파라미터에 민리하게 반응하였다: 높은 SPL, 공진 근처의 주파수에서 증가하였으며, 음파의 전파 방향에 의존적이었다.

3. 전단층 약화 및 유동 역학:

   • $\delta^*$의 횡방향 성장은 상류 구멍에 비해 하류 구멍의 전단층을 약화시키는 결과를 초래했다.
   • 이러한 약화는 공동 내부의 유동 역학을 변화시켰다. 구체적으로, 하류 공동에서의 감소된 전단 강도는 음향 유도 속도장이 공동 내부로 더 많이 침투할 수 있게 하였다.
   • 결과적으로, 음파가 평균 유동에 대항하여 진행할 때(하류 음원), 파동이 더 약한 전단층과 더 큰 유효 구멍 면적을 만나게 되어 공동을 통하는 음향 유도 질량 유량이 더 높아졌다.

4. 음향 유도 유동 비대칭성:

   • 유동장은 평균 유동에 대한 음파의 전파 방향에 따라 비대칭적이다. 상류로 이동하는 파동은 라이너 위에서 아직 완전히 발달하지 않은 경계층(더 얇은 $\delta^*$, 더 강한 전단)을 마주하는 반면, 하류로 이동하는 파동은 완전히 발달하여 더 두꺼워진 경계층(더 약한 전단)을 마주한다.
   • 이러한 비대칭성은 음향 에너지 소산이 라이너를 따라 균일하지 않으며, 진화하는 유동 토폴로지에 따라 공간적으로 변한다는 것을 시사한다.

의의 및 결론
본 연구는 유동이 존재하는 환경의 음향 라이너를 특성화할 때 공간적으로 균질한 임피던스와 정적인 경계층을 가정하는 것은 불충분하다고 결론짓는다. 난류 경계층의 발달은 국부 전단층 강도를 크게 변화시키며, 이는 다시 음향 유도 질량 유량과 에너지 소산을 조절한다.

저자들은 다음과 같이 주장한다:

• 경계층 변위 두께($\delta^*$)는 라이너를 따라 크게 변하는 핵심 파라미터이며, 임피던스 모델에서 반드시 고려되어야 한다.
• 균일한 임피던스나 단순화된 경계 조건을 가정하는 현재의 임피던스 추출 방법들은 공간적으로 진화하는 유동-음향 상호작용을 포착하는 데 실패할 수 있다.
• 평균 유동에 대한 음파의 전파 방향은 서로 다른 유동 환경(서로 다른 $\delta^*$ 및 전단 프로파일)을 생성하며, 이는 특히 모드 매칭과 같은 전역 피팅 방법을 사용할 때 측정된 임피던스에 영향을 미친다.

본 논문은 이러한 고충실도 시뮬레이션에 대한 포괄적인 오픈 액세스 데이터베이스를 제공하여, 음향-유동 결합 및 더 견고한 반경험적 임피던스 추정 모델 개발을 위한 벤치마크 역할을 한다. 연구 결과는 향-후 라이너 특성화 및 추출 기법이 균일한 상태를 가정하기보다 공간적으로 진화하는 난류 유동을 고려해야 함을 시사한다.