Time-Domain Two-Way Fluid–Structure Interaction Analysis of the HIRENASD Transonic Wing: Static Aeroelastic Response and Flutter Prediction

기술 요약: HIRENASD 천음속 익형의 시간 영역 양방향 유체-구조 상호작용 분석

문제 정의
수송용 항공기 익형의 구조적 무결성은 복잡한 에어로엘라스틱 결합, 특히 파괴적인 실패를 초축할 수 있는 플러터(flutter) 현상에 의해 결정된다. Doublet Lattice Method(DLM)와 같은 선형 주파수 영역 방법은 이동 충격파, 충격 유도 경계층 박리, 강한 유체-구조 상호작용을 포함하여 천음속 유동에 내재된 중요한 비선형 현상을 포착하는 데 정밀도가 부족한 경우가 많다. High Reynolds Number Aero-Structural Dynamics (HIRENASD) 프로젝트는 초임계 익형(BAC3-11 에어포일)에 대한 포괄적인 실험 데이터베이스를 제공하지만, 상용 범용 CFD-CSD 솔버를 사용하여 이러한 구성을 검증하는 것과 관련하여 공개 문헌에는 여전히 상당한 공백이 존재한다. 또한, 여러 마하 수에 걸친 체계적인 시간 영역 플러터 경계 식별과 자유류 밀도 및 마하 수가 플러터 동압에 공동으로 미치는 영향에 대한 정량적 분석은 이 벤치마크에 대해 충분히 다뤄지지 않았다.

방법론
본 연구는 ANSYS Workbench 내에 구현된 시간 영역, 느슨하게 결합된(loosely-coupled) 양방향 유체-구조 상호작용(FSI) 프레임워크를 채택한다. 이 접근 방식은 다음과 같이 솔루션을 분할한다:

• CFD: SST $k-\omega$ 난류 모델을 사용하는 밀도 기반 솔버인 ANSYS Fluent를 활용한 RANS(Reynolds-Averaged Navier–Stokes) 방정식. 유체 영역은 $y^+ \approx 1$을 달기 위해 프리즘 레이어를 포함한 약 60만 개의 셀로 구성된 비정렬 사면체 격도를 특징으로 한다.
• CSD: HIRENASD 익형 구조(니켈 마라징강)의 육면체 유한 요소 모델(41,923개 요소)을 사용하는 ANSYS Mechanical. 모드 해석에는 Block-Lanczos 알고리즘을 사용하며, 과도 응답 적분에는 Newmark-$\beta$ 방법을 사용한다.
• 결합(Coupling): ANSYS System Coupling은 힘의 보존을 위한 GGI(General Grid Interface)와 변위 매핑을 위한 Smart Bucket 보간을 통해 데이터 교환을 관리한다. 격자 변형은 국부적 재격자를 포함한 스프링 유사(spring-analogy) 평활화 방식을 통해 처리된다.

본 연구는 NASA FUN3D 공력 벤치마크 및 HIRENASD 실험 구조 주파수에 대해 프레임워크를 검증한다. 정적 에어로엘라스틱 해석은 마하 수($M$) 0.7 및 0.8에서 받음각 $-1.5^\circ$부터 $+4.5^\circ$까지 수행된다. 동적 플러터 해석은 두 가지 뚜렷한 비행 조건인 $M=0.8$ ($\rho = 0.2525$ kg/m³) 및 $M=0.7$ ($\rho = 0.50$ kg/m³)에서 자유류 속도를 체계적으로 변화시켜 감쇠(sub-critical), 임계(critical), 발산(super-critical) 응답 영역을 특성화함으로써 플러터 경계를 식별한다.

주요 기여

1. 프레임워크 검증: 본 연구는 독점적인 인하우스 코드나 전용 솔버 없이도 상용 도구(ANSYS Fluent, Mechanical, System Coupling)만을 사용하여 천음속 에어로엘라스틱성을 위한 재현 가능한 산업 접근 방식의 FSI 워크플로우를 성공적으로 구축하였다.
2. 정적 에어로엘라스틱 벤치마킹: 본 프레임워크는 NASA FUN3D 참조 데이터 대비 1.87%에서 3.52% 이내의 오차로 정상 상태 공력 계수($C_L, C_D, C_M$)를 재현하며, 실험적 지면 진동 시험 결과와 1.94% 이내의 오차로 기본 굽힘 주파수(26.504 Hz)를 예측한다.
3. 동적 플러터 특성화: 본 연구는 시간 영역에서의 체계적인 플러터 경계 결정을 제공하며, 감쇠 진동에서 발산 진동으로의 전이를 명시적으로 특성화한다. 또한, 테스트된 두 마하 영역 모두에서 플러터가 단일 모드 굽힘 메커니즘에 의해 지배됨을 확인하였다.
4. 스케일링에 대한 물리적 통찰: 본 연구는 마하 수, 자유류 밀도 및 플러터 동압 사이의 관계를 정량화하며, 플러터 속도는 낮은 밀도(고고도)에서 더 높을 수 있지만, 스케일링 법칙 $q_f \propto \rho V_f^2$에 따라 플러터 동압은 높은 밀도에서 현저히 높다는 것을 입증한다.

결과

• 정적 응답: $M=0.8$에서 양력 곡선 기울기는 약 $0.107/^\circ$이다. 최대 익단 변위는 받음각 $4.5^\circ$에서 12 mm($\delta/b = 0.0093$)에 도달한다. $M=0.8$에서의 항력 폴라(drag polar)는 급격한 충격 유도 항력 상승을 보이는 반면, $M=0.7$ 폴라는 더 포물선 형태이며 유도 항력에 의해 지배된다.
• 모드 해석: 제1차 수직 굽힘 모드는 26.504 Hz로 예측되었다(실험 대비 1.94% 오차). 고차 모드들은 유한 요소 모델의 구조적 단순화(예: 체결부의 순응성 결여)로 인해 9~13%의 체계적인 과대 예측을 보였다.
• 플러터 경계:
  • $M=0.8$ ($\rho = 0.2525$ kg/m³)에서: 플러터 속도($V_f$)는 266.56 m/s이며 플러터 동압($q_f$)은 8.97 kPa이다. 플러터 주파수는 약 22 Hz로, 이는 천음속 공력 강성 변화로 인해 진공 상태의 구조 주파수로부터 16.9% 감소한 수치이다.
  • $M=0.7$ ($\rho = 0.50$ kg/m³)에서: 플러터 속도는 230.46 m/s이며 플러터 동압은 13.27 kPa이다. 플러터 주파수는 약 25.5 Hz로, 구조 주파수로부터 3.8%의 더 작은 감소를 보였으며, 이는 천음속 비선형성이 더 약함을 나타낸다.
• 스케일링 법칙: $M=0.7$에서 플러터 속도가 더 낮음에도 불구하고, 98% 더 높은 자유류 밀도로 인해 플러터 동압은 $M=0.8$보다 47.9% 더 높다. 이는 차원적 스케일링 $q_f \propto \rho V_f^2$를 확인시켜 준다.

의의 및 주장
본 논문은 검증된 ANSYS Workbench 프레임워크가 전용 인하우스 코드나 독점 솔버 없이도 수송 범주 항공기 익형의 예비 플러터 경계 평가를 위한 실용적이고 산업적으로 유의미한 경로를 제공한다고 주장한다. 본 연구는 상용 도구만을 사용하는 것이 산업 설계 환경으로의 즉각적인 전이 가능성을 보장한다는 점을 강조한다.

저자들은 밀도 의존적 플러터 동압 스케일링에 대한 정량적 확인이, 특히 고도에 따른 플러터 위험과 관련하여 윈드 터널 모델 결과를 실제 비행 조건으로 외삽하기 위한 과학적으로 엄격한 기초를 제공한다고 단언한다. 본 연구는 동일한 마하 수에서도 고도에 따른 효과적인 공력 강성 감소로 인해, 저고도(고밀도)에서는 플러터가 발생하지 않더라도 순항 고도(저밀도)에서는 플러터 경계에 도달할 수 있음을 강조한다. 본 연구는 식별된 플러터 속도가 AGARD 445.6 익형에 대해 검증된 유사 프레임워크와의 비교를 바탕으로 5~8%의 추정 방법론적 불확실성을 가진다고 결 결론지으며, 관찰된 단일 모드 굽림 메커니즘은 이 특정 구성에서 굽힘과 비틀림 모드 사이의 큰 주파수 분리와 일치한다.