Algorithmic Synthesis and Parametric Modeling of Dual Bio-Inspired Airfoils Featuring Leading-Edge Tubercles and Trailing-Edge Serrations

기술 요약: 이중 생체 모방형 에어포일의 알고리즘 합성 및 파라메트릭 모델링

문제 정의
본 논문은 두 가지 뚜렷한 생체 모방형 변형, 즉 전연 결절(humpback whale의 지느러미에서 영감을 얻음)과 후연 톱니 구조(nocturnal owl의 깃털에서 영감을 얻음)를 포함하는 복잡한 3차원 에어포일 기하학을 효율적으로 생성하는 과제를 다룬다. 표준 CAD 도구는 개별 표면 패치를 모델링할 수는 있지만, 해석적 프로파일 방정식과 스팬 방향(spanwise)의 사인파 스윕(sinusoidal sweep) 및 자동 절단 연산을 동적으로 연결하는 고유한 능력은 부족하다. 이러한 한계는 실속 특성 개선 및 공력 소음 저감을 목표로 하는 전산유체역학(CFD) 연구에 필요한 신속한 반복 설계를 저해한다.

방법론
저자는 Grasshopper(Rhinoceros 3D의 비주얼 프로그래밍 플러그인) 내에서 구현된 완전한 파라메트릭 모델링 프레임워크를 제안한다. 이 워크플로우는 에어포일 프로파일을 생성하기 위해 임베디드된 Python 3 스크립트를 통합하고, 공간 조작 및 불리언(Boolean) 연산을 위해 Grasshopper 컴포넌트를 활용한다. 프로세스는 다음 단계를 따른다:

1. 수학적 기초: 핵심 기하학은 NACA 4-digit 에어포일 시리즈를 기반으로 한다. 두께 포락선(thickness envelope)은 조정 가능한 4차 계수($a_4$)를 통해 수정되어, 고전적인 날카로운 후연에서 벗어나 사용자가 지정한 후연 두께 비율을 강제한다. 평균 캠버 라인(mean camber line)과 표면 좌표는 표준 NACA 방정식을 사용하여 계산되며, 곡률이 가장 높은 곳에서의 해상도를 보장하기 위해 코사인 간격(cosine spacing)을 사용하여 전연 근처에 점들을 클러스터링한다.
2. 전연 결절(Leading-Edge Tubercles): 사인 함수가 블레이드 스팬($z$)을 따라 국부 코드 길이($c$)를 변조한다. 이러한 변화는 특징적인 "혹등고래" 형태의 물결 모양을 생성한다. 진폭($A_{tub}$), 파장($\lambda$), 위상 변화($\delta$)는 사용자 정의 슬라이더에 의해 제어된다.
3. 후연 톱니 구조(Trailing-Edge Serrations): 루트 평면(root plane)에서 생성된 톱니 모양의 폴리라인이 스팬 방향으로 돌출(extrusion)되어 커터(cutter)를 형성한다. 불리언 차집합 연산을 통해 이 부피를 메인 블레이드 본체에서 제거함으로써 톱니형 노치(notch)를 생성한다.
4. 자동화: Python 스크립트는 변조된 코드 길이에 따라 각 스팬 방향 스테이션에 대한 폐쇄형 B-스플라인 곡선을 생성한다. 이 곡선들은 로프트(loft)되어 수밀(watertight) 솔리드 볼륨을 생성하며, 이후 불리언 차집합에 의해 수정된다.

주요 기여

• 재현 가능한 파라메트릭 파이프라인: 본 논문은 사용자가 슬라이더를 조정하는 것만으로 기하학적으로 정확하고 수밀한 블레이드를 재생성할 수 있는 완전한 엔드 투 엔드(end-to-end) 워크플로우를 제공하며, 이를 통해 수동 재구축의 필요성을 제거한다.
• 오픈 소스 구현: 에어포일 생성기를 위한 주석이 달린 전체 Python 소스 코드와 함께 상세한 Grasshopper 컴포넌트 레이아웃을 제공한다.
• 분리된 파라미터 제어: 본 연구는 결절(진폭 및 파장)과 톱니(깊이)를 위한 기하학적 파라미터가 독립적으로 작동함을 입증한다. 이를 통해 설계자는 의도하지 않은 기하학적 간섭 없이 유동 제어 기능(결절)과 소음 저감 기능(톱니)을 각각 별도로 튜닝할 수 있다.
• CFD 준비성: 본 방법론은 표면 연속성($C^2$), 결절 골(trough)에서의 경계층 격자 요구 사항, 호환 가능한 내보내기 형식(STEP, Parasolid)을 강조하며 CFD를 위한 기하학 준비에 관한 구체적인 노트를 포함한다.

결과
주요 생체 모방 파라미터에 대한 블레이드의 기하학적 반응을 분석하기 위해 단일 변수 민감도 연구가 수행되었다:

• 결절 깊이 ($A_{tub}$): 진폭이 증가하면 파형 주기(wave cycles)의 수는 유지하면서 코드 포락선($\Delta c = 2A_{tub}$)이 선형적으로 넓어진다. 그러나 진폭이 매우 높으면 로프트 알고리즘에 어려움을 줄 수 있는 극도로 얇은 단면이 발생할 수 있다.
• 결절 파장 ($\lambda$): 파장을 줄이면 코드 진폭을 유지하면서 스팬 방향을 따르는 주기 수($N_{cycles}$)가 증가한다. 짧은 파장은 표면 패싯(faceting) 현상을 방지하기 위해 더 높은 점 개수($n_{pts}$)를 요구한다.
• 톱니 깊이 ($A_{ser}$): 이 파라미터는 불리언 커터에만 영향을 미친다. 이는 코드 분포나 전연 기하학에 영향을 주지 않는다. 연구에서는 노치 깊이가 프로파일 내부를 관통하지 않도록 베이스 코드의 50%를 초과해서는 안 된다는 실질적인 상한선을 언급하였다.

의의 및 주장
본 논문은 연구자와 엔지니어가 이중 생체 모방형 에어포일의 설계 공간을 탐색할 수 있는 실용적이고 접근 가능한 도구를 제공한다고 주장한다. 복잡한 기하학적 형상을 생성하는 과정을 자동화함으로써, 본 워크플로우는 미래의 CFD 시뮬레이션 및 실험적 검증에 필요한 신속한 프로토타이핑을 용이하게 한다. 저자는 현재의 작업이 기하학적 생성에 국한되어 있으며, 공력 예측, 격자 생성 또는 솔버 구성은 포함하지 않는다고 명시적으로 밝힌다. 본 연구의 의의는 향후 스팬 방향 비틀림 보간(spanwise twist interpolation) 및 사용자 정의 에어포일 좌표 지원을 포함할 계획인, 향-미래의 공력 및 음향 연구를 구축할 수 있는 견고한 제일 원리(first-principles) 기반의 토대를 확립하는 데 있다.