제1부 비행체 기술 기반
1.5 공력 해석·풍동·CFD 설계 도구
1.5.1 이 장에서 배우는 것 / 선수 지식 / 특히 중요한 독자
앞 장들의 계산은 양항비, 효율, 최대 양력 계수 같은 숫자를 입력으로 가정했다. 이 숫자를 어디서 얻고, 얼마나 믿을 수 있는가가 이 장의 질문이다. 설계자는 손으로 푸는 식에서 수백만 격자의 CFD(전산유체역학)까지, 풍동 시험에서 비행시험까지 여러 도구를 단계마다 바꿔 사용한다. 도구마다 무엇을 예측하고 무엇을 못 하는지를 아는 것이 도구를 사용하는 능력의 절반이다.
배우는 것:
- 해석 도구의 계층 — 양력선·패널법·와류격자법(VLM)에서 레이놀즈 평균 나비에-스토크스(RANS)·큰에디모사(LES)까지 정확도와 비용의 사다리.
- 소형 무인기의 지배 문제 — 낮은 레이놀즈수의 층류 박리 기포와 천이 예측의 어려움.
- 로터·프로펠러 CFD — 회전 기준틀(MRF)·슬라이딩 메시·액추에이터 모델의 선택.
- 풍동 시험 — 회로 형식, 막힘·벽 간섭 보정, 로터 시험대, 불확도, 비행시험과의 연결.
- 검증과 확인(V&V) — 격자 수렴, 난류 모델 불확도, 오차와 불확도의 구분.
- 도구 선택과 대리 모델 — 오픈소스 대 상용, 머신러닝 보조 설계의 용도과 위험.
선수 지식: 1.1(양력·항력·경계층의 기본), 1.2·1.3(로터·프로펠러 공기역학). 식은 대수 수준이다. 비엔지니어는 각 절의 결론 문장과 표를 이어 읽으면 어떤 수치를 얼마나 믿어도 되는지 감을 잡을 수 있다.
특히 중요한 독자: 엔지니어(도구와 한계), 인증 담당자(해석·시험 자료를 입증 근거로 사용하는 조건), 사업·정책 담당자(CFD 결과를 보고서에서 읽는 법).
기준 시점과 정직성. 이 장은 2026년 10월을 기준으로 하고, 집필에 사용한 지식은 2026년 중반까지다. 이 장의 수치는 교육용 코드(와류격자법·경계층 방정식)와 교과서 경험식으로 직접 계산한 값이며, 실제 CFD 상용 코드의 결과나 시험 자료가 아니다. 계산은 파이썬(
math만)으로 했다. 도구의 소속·기능은 확인한 범위에서만 적고 나머지는(확인 필요)로 둔다.
1.5.2 개념·직관 — 모든 해석은 무언가를 버리고 얻는다
공기의 흐름은 나비에·스토크스 방정식이 정한다. 이 방정식을 그대로 풀려면(직접수치모사, DNS) 가장 작은 소용돌이까지 해상해야 해서 기체 하나에 슈퍼컴퓨터로도 현실적이지 않다. 그래서 해석 도구는 모두 무언가를 버린다. 점성을 버리거나, 난류를 평균으로 대신하거나, 3차원을 2차원으로 줄인다. 버린 것이 결과에 영향을 주지 않는 영역에서는 싸고 정확하고, 영향을 주는 영역에서는 빠르게 틀린다.
그래서 설계에서는 계층적으로 사용한다. 초기에는 값싼 도구로 수백 개 형상을 훑고, 후보가 좁혀지면 정교한 도구로 몇 개만 확인하고, 마지막에 풍동 시험과 비행시험으로 닫는다. 위쪽 도구로 정해진 결론을 아래쪽 도구가 뒤집는 지점이 어디인지를 아는 것이 핵심이다.
세 축은 같은 도구를 다른 조건에서 사용한다.
- 군용 드론. 소형 정찰기는 낮은 레이놀즈수, 고고도 장기 체공기는 낮은 밀도로 인한 낮은 레이놀즈수, 고속 기체는 압축성이 문제다. 시험 시설과 해석 코드의 접근이 수출통제 대상일 수 있다(19.16).
- 민간 UAS. 값싼 풍동·오픈소스 도구와 비행시험으로 반복 개량한다. 작은 기체는 실제 크기로 풍동에 넣을 수 있다는 이점이 있다.
- AAM. 인증이 해석과 시험을 입증 수단으로 요구하므로 도구의 검증 이력과 불확도가 논점이 된다(20.6).
1.5.3 해석 도구의 계층 — 정확도와 비용
| 계층 | 가정 | 예측하는 것 | 못 하는 것 |
|---|---|---|---|
| 양력선 이론(Prandtl) | 날개를 하나의 양력선으로, 단면은 2차원 자료 사용 | 큰 종횡비 직선 날개의 양력 분포·유도 항력 | 작은 종횡비·큰 후퇴각·실속 |
| 패널법·와류격자법(VLM) | 비점성·비압축 퍼텐셜 흐름. 표면이나 캠버면에 와류·특이점 배치 | 임의 형상의 양력·모멘트·양력 분포·유도 항력 | 점성·분리·실속 |
| 경계층 결합 패널법(XFOIL 등) | 퍼텐셜 흐름에 적분 경계층과 $e^N$ 천이를 결합(2차원) | 낮은 $Re$ 익형 극선, 작은 층류 박리 기포 | 큰 분리, 3차원 |
| RANS | 난류를 평균하고 난류 모델(Spalart–Allmaras 모델 SA, $k$-$\omega$ 전단응력수송 모델 SST 등)로 닫음. 선택적으로 천이 모델 | 3차원 전기체·로터의 점성 흐름 | 큰 분리·실속에서 모델 의존 |
| URANS·DES | 비정상 계산(URANS는 비정상 RANS). DES(분리와류모사)는 벽 근처는 RANS, 분리 영역은 LES | 비정상 분리·후류·간섭 | 격자·시간 비용 |
| LES·DNS | 큰 소용돌이(LES)나 모든 척도(DNS)를 직접 풂 | 난류 구조·소음원 | 실기체에는 비용이 비현실적 |
결론부터. 모든 도구가 같은 문제를 같게 푸는 것은 아니다. 가장 아래쪽 비용을 사용한다고 해서 항상 더 정확한 것도 아니다. 값싼 도구의 오차가 문제의 어느 영역에서 커지는지 보이는 간단한 예가 있다. 직사각 날개의 양력 기울기를 세 방법으로 계산했다(평판, 단면 기울기 $a_0=2\pi$).
$$ a_{\mathrm{LLT}}=\frac{a_0}{1+a_0/(\pi AR)},\qquad a_{\mathrm{Helmbold}}=\frac{a_0}{\sqrt{1+\left(a_0/(\pi AR)\right)^2}+a_0/(\pi AR)} $$
| 종횡비 $AR$ | 양력선(타원 분포) | Helmbold 근사 | VLM(교육용 코드, 스팬 40패널) | 세 값의 퍼짐 |
|---|---|---|---|---|
| 4 | 4.19 | 3.88 | 3.63 | 15% |
| 6 | 4.71 | 4.53 | 4.25 | 11% |
| 10 | 5.24 | 5.15 | 4.89 | 7% |
| 20 | 5.71 | 5.69 | 5.49 | 4% |
(단위는 1/rad. VLM은 이 책의 교육용 코드로, 시위 방향에 패널 하나를 두었다.) 세 값의 퍼짐은 (최대-최소)/최소다. 종횡비가 클수록 세 방법이 가까워지고, 작을수록 벌어진다. 양력선 이론이 큰 종횡비를 가정했기 때문이다. 이 퍼짐 자체가 도구의 모델 형식 불확도의 한 척도다. 어느 값이 맞는지는 시험이 판정한다. 소형 무인기의 날개는 종횡비 6~10이라 이 퍼짐이 7~11%에 이르고, 양력 기울기 하나가 안정성 해석(1.1)·제어 설계(2.2)에 곧바로 들어간다.
1.5.4 소형 무인기의 지배 문제 — 낮은 레이놀즈수
소형 무인기의 날개와 프로펠러는 레이놀즈수 $Re=Vc/\nu$가 $10^4\sim10^6$대다. 이 영역의 흐름은 층류가 쉽게 떨어진다. 층류 경계층은 압력이 오르는 구간(날개 윗면 뒤쪽)에서 부착을 유지하지 못하고 표면에서 떨어지고, 떨어진 얇은 전단층이 난류로 바뀌며 다시 붙는다. 이 안쪽이 층류 박리 기포다. 기포는 항력을 키우고 양력 곡선을 일그러뜨리며, 받음각을 올릴 때와 내릴 때 값이 달라지는 이력(히스테리시스)을 만든다.
| 비행체(가정 입력) | 속도 [m/s] | 시위 [m] | $Re$ |
|---|---|---|---|
| 마이크로 비행체 | 5 | 0.05 | 약 1.7만 |
| 소형 UAS | 15 | 0.15 | 약 15만 |
| 중형 UAS | 30 | 0.4 | 약 82만 |
| 대형 중고도 장기 체공기 | 60 | 1.0 | 약 410만 |
| AAM 날개 | 70 | 1.5 | 약 720만 |
| 고고도 장기 체공기(18 km, $\rho\approx0.12$ kg/m³) | 100 | 1.0 | 약 85만 |
(공기의 동점성 계수 $\nu=1.46\times10^{-5}$ m²/s로 계산, 마지막 줄은 표준대기 18 km에서 $\nu\approx1.2\times10^{-4}$ m²/s. 입력은 예시이며 제품 값이 아니다.) 마지막 줄이 눈에 띈다. 큰 기체라도 높이 날면 소형 기체의 레이놀즈수로 날아 같은 문제가 생긴다(1.15).
그림은 평판 마찰만으로 본 최대 양항비의 상한이다. 층류 유지(파랑)는 높고, 난류(주황)는 낮다. 천이가 시작되는 $Re\approx5\times10^5$ 부근에서 검은 곡선이 꺾이는 것은 층류 이점을 잃는 대가다. 낮은 $Re$일수록 마찰이 크므로 상한 자체가 낮다. 실제 익형은 층류 박리 기포와 압력 항력 때문에 이 상한보다 낮고, 그 간격이 낮은 $Re$에서 가장 크다. 측정값은 UIUC의 낮은 속도 익형 시험 자료집 등을 직접 읽어야 한다(1.5.14).
천이를 예측하기 어렵다. 천이 위치는 자유류 난류 강도, 표면 거칠기, 소음·진동, 받음각에 민감하다. 해석에서는 $e^N$ 법(경계층 안 교란이 $N$ 배 증폭되면 천이로 본다. XFOIL의 기본값은 $N_{crit}=9$로 알려져 있다)이나 천이 모델(Langtry와 Menter의 $\gamma$-$Re_\theta$ 모델 등)을 사용한다. 값 하나를 바꾸면 항력이 크게 바뀌므로 같은 코드 안에서도 입력(난류 강도)이 결과를 결정한다. 풍동 시험에서는 천이를 일부러 일으키는 트립(트립 와이어·거친 띠)을 사용하기도 하고, 이는 비행 조건과의 대응을 따로 따져야 한다(1.5.6).
1.5.5 로터·프로펠러 CFD
회전하는 깃은 CFD에서 가장 까다로운 문제 중 하나다. 깃이 움직이며 후류와 만나고 기체와 간섭하기 때문이다. 방법은 비용과 보는 것에 따라 달라진다.
| 방법 | 원리 | 보는 것 | 못 보는 것·비용 |
|---|---|---|---|
| 액추에이터 디스크·라인 | 깃을 체적력으로 대체. 하중은 BEMT(1.3)에서 | 로터 후류, 기체와 로터의 간섭, 지면·건물 근처 흐름 | 깃 주위 흐름·깃끝 와류의 세부. 비용이 가장 작다 |
| 회전 기준틀(MRF) | 로터 둘레 영역을 회전 좌표계에서 정상으로 풂 | 호버 성능, 로터 한 개의 정상 상태 | 로터와 기체의 상대 위치가 고정되어 비정상 간섭은 못 본다 |
| 슬라이딩 메시 | 로터 영역 격자가 실제로 회전하며 비정상 계산 | 로터-기체 간섭, 소음원의 시간 변화 | 시간 간격이 작고 비용이 크다 |
| 오버셋(중첩) 격자 | 로터·기체 격자를 중첩해 상대 운동 | 여러 로터·전환 비행 | 격자 사이 보간 오차, 구현 복잡 |
호버 성능 예측의 정확도는 낮다. AIAA 호버 예측 워크숍(2014년부터 SciTech 포럼에서 열림)을 정리한 NASA·육군 자료는 최상의 CFD 호버 성능 예측 오차가 대체로 2~5%이고 0.5%가 목표라고 말한다 (2차 출처, 확인 필요). 성능지수 5%의 오차는 호버 동력 5%의 오차이고, 로터 지름을 5% 바꾸는 효과와 같은 크기다(1.2.3). 낮은 $Re$ 소형 로터는 이 정확도도 보장되지 않는다.
실무 선택. 로터 한 개의 호버 성능은 BEMT(1.3)와 MRF로, 기체와 로터의 간섭이나 건물 근처의 흐름은 액추에이터 모델이나 슬라이딩 메시로 본다. 소음은 CFD와 별도의 음향 해석이 필요하다(1.31).
1.5.6 풍동 시험과 비행시험에서의 식별
풍동은 공기를 흐르게 해 모형을 고정하고 힘을 재는 장치다. 회로 형식은 둘이다. 개방 회로는 공기를 밖에서 빨아들여 시험부를 지나 내보내고, 값싸고 단순하지만 외기와 소음의 영향을 받는다. 폐쇄 회로는 공기가 고리를 돈다. 흐름의 질이 좋고 전력이 적게 들지만 시설이 크다. 시험부도 벽이 있는 폐쇄형과 터진 개방형이 있다. 폐쇄형은 벽이 흐름을 좁혀 모형 둘레의 속도가 빨라지고, 개방형은 반대 방향의 오차를 만든다. 둘 다 보정이 필요하다.
막힘과 벽 간섭 보정. 모형이 시험부를 가로막으면 모형 둘레의 속도가 올라간다. 교과서(Barlow·Rae·Pope)의 폐쇄형 보정은 고체 막힘 $\varepsilon_{sb}$(모형 부피), 후류 막힘 $\varepsilon_{wb}=\dfrac{S}{4C}C_{D,u}$(후류 때문)를 더해 총 막힘 $\varepsilon_t$로 속도와 동압을 고친다. 양력이 있으면 벽이 받음각과 항력도 바꾼다.
$$ V_c=V_u(1+\varepsilon_t),\quad q_c=q_u(1+\varepsilon_t)^2,\quad C_{L,c}=\frac{C_{L,u}}{(1+\varepsilon_t)^2},\quad \Delta\alpha=\delta\,\frac{S}{C}\,C_{L,u}\ (\text{rad}),\quad \Delta C_D=\delta\,\frac{S}{C}\,C_{L,u}^2 $$
여기서 $S$는 모형 날개 면적, $C$는 시험부 단면적이다. $\delta$는 시험부 형상과 모형 스팬비로 정해지는 계수다. 소형 무인기 날개를 가정한 계산이다.
| 입력(가정) | 값 |
|---|---|
| 모형 | 날개 면적 $S=0.09$ m²(스팬 0.75 m·시위 0.12 m), 시험부 $1.2\times0.9$ m($C=1.08$ m²) |
| 시험 조건 | $V=20$ m/s, $C_{L,u}=0.8$, $C_{D,u}=0.05$, $Re\approx1.6\times10^5$ |
| 보정 계수 | $\varepsilon_{sb}=0.005$(가정), $\delta=0.1$(가정, 확인 필요) |
계산하면 $S/C=0.083$, 후류 막힘 $\varepsilon_{wb}=0.0010$, 총 $\varepsilon_t=0.0060$이다. 속도가 0.6%, 동압이 1.2% 오르고 보정한 양력 계수는 0.790이다. 벽 간섭으로 받음각은 $0.38^\circ$ 보정되고 항력 계수는 0.0053(10.7%) 보정된다. 작은 보정으로 보이는 막힘보다 벽 간섭의 항력 보정이 크다. 시험부에 비해 모형이 클수록 두 보정이 모두 커진다.
소형 무인기의 이점과 제약. 시위 0.12 m의 실물 날개를 20 m/s로 시험하면 $Re=1.6\times10^5$로 비행과 같다. 대형 기체는 축소 모형으로 비행 $Re$를 못 맞춘다. 대신 소형 날개의 저 $Re$ 시험은 풍동의 난류 강도와 표면 상태에 매우 민감하다. 시험 결과를 비행에 옮길 때는 풍동의 난류 강도를 확인한다 (경험, 확인 필요).
로터 시험대. 호버 시험은 추력·토크·회전수를 재며, 막힘이 아니라 재순환(후류가 벽과 바닥에 닿아 돌아옴)과 지면 효과(1.2.7)가 오차 원인이다. 시험대의 구성과 오차는 1.3.9에 정리했다. 전진 비행의 성능은 풍동에서 얻는다.
불확도. 측정에는 반드시 불확도를 붙인다. $C_L=L/(qS)$를 예로, 양력 1.0%, 속도 0.25%(동압 0.5%), 공기 밀도 0.3%, 면적 0.2%의 상대 표준 불확도가 서로 독립이면 합성 표준 불확도는 $\sqrt{1.0^2+0.58^2+0.2^2}=1.17\%$, 포함인자 2의 확장 불확도는 2.3%다(계산 입력은 가정). 방법은 GUM(JCGM 100)과 풍동 시험에 맞춘 AIAA S-071A가 정한다.
비행시험에서 식별. 풍동·CFD는 형상 안의 공력 모델을 주지만 실제 기체와는 어긋난다. 비행시험에서 입력을 주고 응답을 재어 공력 도함수를 추정하는 것이 시스템 식별이다(2.1·8.20). 풍동·CFD 값은 사전 모델로 사용하고, 비행시험 자료로 보정한다. 순서는 해석→풍동→비행이고, 앞 단계의 불확도가 클수록 비행시험 범위를 넓게 잡는다.
1.5.7 검증과 확인(V&V) — 이 숫자를 얼마나 믿을 수 있는가
두 낱말을 구분한다. 검증(verification)은 수학 모델을 올바르게 풀었는가(코드 오류·격자·수렴), 확인(validation)은 올바른 수학 모델을 골랐는가(시험과 비교)이다. 둘을 섞으면 "시험과 맞으니 격자는 충분하다" 같은 틀린 결론에 이른다.
격자 수렴. 격자를 세밀하게 해 값이 어떻게 변하는지로 이산화 오차를 가늠한다. 세 격자(정밀 $f_1$, 중간 $f_2$, 조대 $f_3$, 세분 비율 $r$)로 관측 차수 $p$, 외삽값, 격자 수렴 지수(GCI, Roache)를 구한다.
$$ p=\frac{\ln\left|\dfrac{f_3-f_2}{f_2-f_1}\right|}{\ln r},\qquad f_{ext}=f_1+\frac{f_1-f_2}{r^{p}-1},\qquad \mathrm{GCI}=\frac{1.25\,\left|(f_2-f_1)/f_1\right|}{r^{p}-1} $$
절차가 맞는지는 정답을 아는 문제에서 먼저 확인한다. 평판 위 층류 경계층(Blasius 방정식 $f'''+\tfrac12 ff''=0$)의 벽 기울기 $f''(0)$는 정답이 0.332057이다. 2차 유한차분으로 격자 간격 $h=0.4$, 0.2, 0.1에서 풀면 다음과 같다.
| 격자 간격 $h$ | $f''(0)$ | 정답과의 오차 |
|---|---|---|
| 0.4 | 0.333342 | 0.387% |
| 0.2 | 0.332268 | 0.063% |
| 0.1 | 0.332096 | 0.0117% |
관측 차수는 $p=2.64$(이론 2 부근에서 점근 영역 앞이라 조금 크다), 외삽값은 0.332063이다. 정밀 격자의 $\mathrm{GCI}=0.0123\%$는 실제 오차 0.0117%를 상회한다. GCI가 정직한 상한 구실을 한다는 확인이다.
정답을 모르는 문제에서는 같은 절차가 불확도를 붙여 준다. 앞 절의 직사각 날개($AR=6$)의 VLM 양력 기울기를 스팬 패널 수 10·20·40·80개로 풀면 4.439·4.315·4.250·4.216 /rad다. 가장 고운 세 격자($r=2$)로 $p=0.96$(VLM의 스팬 수렴은 1차에 가깝다), 외삽값 4.18, $\mathrm{GCI}=1.07\%$를 얻는다. 결과는 "4.22 ±1.1%(이산화 불확도)"로 적고, 이 불확도는 다른 방법과의 차이($AR=6$에서 6~10%, 1.5.3)에 비해 작다. 격자를 더 촘촘히 해도 그 차이는 줄지 않는다. 그 차이는 모델의 가정에서 오기 때문이다.
난류 모델 불확도는 격자로 안 줄어든다. 같은 격자에서 난류 모델(SA·SST 등)을 바꾸었을 때의 결과 퍼짐은 모델 형식 불확도의 단서다. 분리·실속 영역에서는 이 퍼짐이 격자 불확도보다 훨씬 크다. 그래서 설계 자료로 CFD 값을 사용할 때는 격자 수렴 지수와 모델 퍼짐을 둘 다 사용한다.
오차와 불확도. 계산 $S$와 시험 $D$의 차이 $E=S-D$를 ASME V&V 20의 틀로 해석한다. 확인 불확도(validation uncertainty)는 세 가지의 합성이다.
$$ u_{val}=\sqrt{u_{num}^2+u_{input}^2+u_D^2},\qquad \text{모델 형식 오차의 추정}\;\;\delta_{model}\approx E\pm u_{val} $$
$\lvert E\rvert$가 $u_{val}$보다 크면 모델 형식의 오차가 신호로 드러난 것이고, 작으면 이 시험의 정밀도로는 모델 오차를 알 수 없다. 이 판단은 시험한 조건에서만 유효하다. 다른 받음각·$Re$로 외삽하는 것은 공학적 판단이다.
보고서를 읽을 때 묻는 다섯 가지. 해석 결과를 근거로 사용하는 정책·인증·사업 담당자도 이 질문으로 신뢰도를 가늠할 수 있다.
- 격자 수렴을 보였는가(격자 세 개 이상, 수렴 지수)?
- 난류·천이 모델은 무엇이고 선택 이유가 있는가? 다른 모델로 바꾸면 결과가 얼마나 변하는가?
- 입력 조건(난류 강도·레이놀즈수·경계 조건)이 실제 비행과 맞는가?
- 시험과 비교했는가? 어떤 조건에서, 시험 불확도는 얼마인가?
- 확인한 범위 밖(다른 받음각·속도·형상)으로 외삽했는가?
1.5.8 도구 선택 — 오픈소스 대 상용, 대리 모델과 머신러닝
| 도구 | 성격 | 용도 | 확인 수준 |
|---|---|---|---|
| XFOIL | M. Drela(MIT). 패널법과 경계층·천이 결합 | 낮은 $Re$ 익형 | 문헌으로 확인 |
| AVL | M. Drela·H. Youngren(MIT). 양력선·VLM | 비행역학용 공력 도함수 | 문헌으로 확인 |
| OpenVSP·VSPAERO | NASA가 2012년 공개한 형상 도구와 그 공력 해석기(VLM·패널법) | 형상에서 해석까지 | 1.4.8 참조 |
| XFLR5 | XFOIL 기반의 오픈소스 해석 도구 | 소형 날개·익형 | (확인 필요) |
| OpenFOAM | 오픈소스 CFD. 재단법인 계열과 ESI(OpenCFD) 계열이 있음 | RANS·LES·로터 MRF | (계열별 확인 필요) |
| SU2 | 스탠퍼드 기원의 오픈소스 CFD와 설계 최적화(Palacios 외, AIAA 2013-0287 (번호 확인 필요)) | 공력 최적화 | 문헌으로 확인 |
| 상용 CFD(ANSYS Fluent·CFX, Siemens STAR-CCM+ 등) | 상업 코드. 지원·검증 자료·전처리 통합 | 산업 전반 | 이름은 예시. 기능·라이선스는 (확인 필요) |
| NASA 코드(FUN3D·OVERFLOW 등) | 정부 개발 비정형·중첩 격자 RANS | 로터·전기체 연구 | 배포가 제한될 수 있다 (확인 필요, 19.16) |
오픈소스는 투명성과 재현성이 강점이다. 코드를 읽고 같은 결과를 다시 만들 수 있으므로 검증 이력과 논문이 쌓인다. 상용은 전처리·지원·사용자 교육이 강점이고, 인증 자료에서 사용된 이력이 협상에서 도움이 될 수 있다. 어느 쪽이든 코드가 아니라 사용 절차(격자·모델·수렴 기준)가 결과를 결정한다. 같은 코드를 다른 사람이 사용하면 다른 답이 나온다는 것이 CFD 워크숍의 반복된 교훈이다.
대리 모델과 머신러닝. 수천 번의 해석이 필요한 최적화와 불확실성 분석에서는 싼 대리 모델이 사용된다. 고충실도 해석 자료로 가우스 과정(크리깅)·신경망·고유직교분해(POD)를 학습해 수 초에 값을 낸다.
- 용도. 설계 공간 탐색, 최적화(1.4.8), 실시간 성능 예측, 해석 결과의 보간.
- 위험 1: 외삽. 학습 영역 밖에서 조용히 틀린다.
- 위험 2: 자료의 대표성. 학습 자료의 오차와 편향이 그대로 들어간다.
- 위험 3: 물리 위반. 보존 법칙을 어기는 값을 낼 수 있다.
- 위험 4: 불확도의 부재. 값은 주지만 얼마나 믿을지를 모른다.
견해: 머신러닝은 CFD를 대체하지 않고 CFD의 비용을 줄이는 수단으로 사용하는 것이 합리적이다. 근거는 학습 자료가 결국 CFD·시험이어서 그 불확도를 넘지 못한다는 점이다. 대리 모델의 결론은 최종 후보에서 원래 도구로 재확인한다.
1.5.9 설계 트레이드오프 — 설계 단계에 어떤 충실도를 쓸까
① 저충실도 중심(양력선·VLM·XFOIL) + 시험 보정.
- 근거: 한 번에 몇 초라 수백 개 형상을 훑고 민감도 분석이 가능하다. 소형 무인기는 시제품이 싸서 시험으로 보정하기 쉽다.
- 대가: 분리·실속·간섭이 지배하는 영역에서 오차가 크다. 보정 시험이 늘어난다.
- 선택하는 쪽: 소규모 민간 UAS 개발, 초기 개념 설계.
② 중충실도 RANS 중심.
- 근거: 전기체·로터 간섭 같은 3차원 점성 효과를 보고, 시험 전에 문제를 찾는다.
- 대가: 비용과 전문성이 필요하고 난류 모델 의존이 남는다. 격자 의존성 검증이 필수다.
- 선택하는 쪽: 소음·간섭이 중요한 AAM, 고속 군용 기체.
③ 풍동·비행시험 중심.
- 근거: 모델 형식 불확도가 없다. 인증이 시험 자료를 가장 높이 평가한다.
- 대가: 시설·시간·모형 제작 비용이 크고, 한 번 시험하면 형상을 바꾸기 어렵다.
- 선택하는 쪽: 인증이 가까운 후기 개발, 해석으로 못 예측하는 현상이 있는 기체.
견해: 어느 한 계층만 사용하는 것보다 계층 사이의 불일치를 설계 위험으로 보는 것이 합리적이다. 근거는 1.5.3의 예다. 도구 사이의 퍼짐이 7~11%이면 그만큼의 마진을 가산하거나 시험으로 좁혀야 한다. 시험 예산은 불일치가 가장 큰 곳(분리·간섭)에 우선 사용한다.
1.5.10 세 축 비교
| 항목 | 군용 드론 | 민간 UAS | AAM |
|---|---|---|---|
| 지배 흐름 문제 | 저 $Re$(소형), 고고도 저밀도, 고속 압축성 | 저 $Re$, 로터 간섭 | 로터·날개 간섭, 천이 비행, 소음 |
| 주 도구 | 내부 코드·수출통제 대상 도구 포함 | 오픈소스·상용·자체 풍동 | 다충실도 CFD, 대형 풍동 |
| 시험 시설 | 국방 시험 시설(접근 제한) | 대학·소형 풍동·비행시험 | 대형 풍동·시험 비행장 |
| V&V 요구 | 임무 요구 대비 충분성 | 제조사 판단 | 인증 입증 수단으로 요구 |
| 보고 | 비공개가 많음 | 공개 논문·오픈 자료 | 인증기관과 협의 |
1.5.11 관련 규격·프로토콜
CFD와 풍동 시험의 절차 규격은 소수이고, 항공기 형식증명이 이를 입증 방법으로 인용한다.
| 문서 | 발행 | 무엇을 정하나 | 상세 |
|---|---|---|---|
| AIAA G-077-1998 | AIAA | CFD 검증·확인 지침(용어·방법) | 20.6 |
| ASME V&V 20-2009 | ASME | CFD·열전달의 검증·확인 표준(불확도 틀) (개정 확인 필요) | 20.6 |
| ANSI/AIAA S-071A-1999 | AIAA | 풍동 시험의 실험 불확도 평가 | 20.6 |
| JCGM 100 (GUM) | BIPM 등 | 측정 불확도 표현 지침 | 이 장 1.5.6 |
| NASA-STD-7009A | NASA | 모델·시뮬레이션 표준(신뢰도 평가) (판 확인 필요) | 8.19 |
| MIL-STD-3022 | 미 국방부 | 모델·시뮬레이션의 검증·확인·인정(VV&A) 문서화 (판 확인 필요) | 8.19 |
1.5.12 관련 규제·인증 요구와의 연결
- 해석·시험을 입증 수단으로. 인증은 시험·해석·검사를 적합성 입증 방법으로 사용한다. 해석을 사용하려면 검증·확인 이력이 필요하다(20.6).
- 모델 신뢰도의 설명. 안전 사례(Safety Case)가 해석의 불확도를 논증에 연결한다(8.7).
- 시험 시설·자료의 통제. 군용 해석 코드와 시험 자료는 수출통제 대상일 수 있다(19.16).
- 소음. 소음 인증은 측정이 기준이므로 CFD는 설계 예측에만 사용된다(1.31·20.11).
- 비행시험과 모델. 비행시험 승인은 시험 전 모델의 불확도를 근거로 범위를 정한다(8.20).
1.5.13 함정과 한계
| 오해 | 실제 | 더 볼 곳 |
|---|---|---|
| "CFD 결과가 시험과 맞으니 틀림없다" | 격자 의존성과 모델 형식을 따로 확인해야 한다 | 1.5.7 |
| "격자를 촘촘히 하면 정확해진다" | 난류 모델·천이 모델의 오차는 격자로 줄지 않는다 | 1.5.7 |
| "정교한 도구일수록 항상 낫다" | 비용 대비 의미가 없는 영역이 있고, 값싼 도구가 더 안정적일 때도 있다 | 1.5.3 |
| "소형 무인기도 대형 기체 방법으로 해석하면 된다" | 낮은 $Re$는 층류 박리 기포와 천이가 지배한다 | 1.5.4 |
| "풍동 데이터는 그대로 비행값이다" | 막힘·벽 간섭·난류 강도·$Re$ 보정이 필요하다 | 1.5.6 |
| "머신러닝이 CFD를 대체한다" | 학습 자료의 불확도를 못 넘고 외삽에서 조용히 틀린다 | 1.5.8 |
| "이 장의 수치가 상용 CFD의 결과다" | 교육용 코드와 경험식의 값이다 | 1.5.1 |
미해결 문제. 낮은 $Re$ 천이의 신뢰할 만한 예측, 로터-기체 간섭의 비정상 해석, 난류 모델 불확도의 정량화, 머신러닝 대리 모델의 불확도 부여, 해석을 입증 수단으로 인정받는 기준의 정립은 연구 중이다 (확인 필요).
1.5.14 체크포인트·연구 질문
체크포인트
- 양력선·VLM·RANS·LES가 각각 무엇을 버리고 무엇을 예측하는지 표로 설명하라.
- 소형 UAS와 고고도 장기 체공기의 레이놀즈수가 비슷해지는 이유를 설명하라.
- 층류 박리 기포가 항력과 양력 곡선에 미치는 영향을 말하라.
- 풍동 시험의 후류 막힘 $\varepsilon_{wb}=SC_D/(4C)$를 이 장의 입력으로 계산하라.
- 세 격자의 값 4.3154·4.2495·4.2156($r=2$)에서 관측 차수 $p$를 구하라.
- 오차 $E$가 $u_{val}$보다 작을 때 무엇을 알 수 있고 무엇을 알 수 없는가?
연구 질문 — 이 장에서 파생되는 열린 주제(25.7에 모은다)
- 저 $Re$ 천이의 불확도 — 같은 익형·같은 코드에서 난류 강도와 거칠기 입력이 항력 예측에 미치는 영향을 공개 시험 자료로 정량화할 수 있는가?
- 로터-기체 간섭의 설계용 모델 — 액추에이터 모델과 슬라이딩 메시의 차이가 설계 결론을 바꾸는 영역은 어디인가?
- 대리 모델의 불확도 부여 — 머신러닝 대리 모델에 신뢰 구간을 붙여 인증의 해석 입증에 사용할 수 있는가?
- 풍동·CFD·비행시험의 불확도 연쇄 — 앞 단계의 불확도가 비행시험 범위 설정에 주는 영향을 정량적 규칙으로 사용할 수 있는가?
1차 출처(이 장)
- J. Katz, A. Plotkin, Low-Speed Aerodynamics, Cambridge University Press (2판, 2001)
- M. Drela, Flight Vehicle Aerodynamics, MIT Press (2014)
- J. B. Barlow, W. H. Rae, A. Pope, Low-Speed Wind Tunnel Testing, Wiley (3판, 1999)
- H. Schlichting, K. Gersten, Boundary-Layer Theory, Springer
- P. J. Roache, 1994, Perspective: A Method for Uniform Reporting of Grid Refinement Studies, J. Fluids Engineering 116(3)
- W. L. Oberkampf, C. J. Roy, Verification and Validation in Scientific Computing, Cambridge University Press (2010)
- T. J. Mueller (편), Fixed and Flapping Wing Aerodynamics for Micro Air Vehicle Applications, AIAA Progress in Astronautics and Aeronautics 195 (2001); T. J. Mueller, J. D. DeLaurier, 2003, Aerodynamics of Small Vehicles, Annual Review of Fluid Mechanics 35
- P. B. S. Lissaman, 1983, Low-Reynolds-Number Airfoils, Annual Review of Fluid Mechanics 15
- UIUC Low-Speed Airfoil Tests(M. S. Selig 연구진 시험 자료집)
- AIAA G-077-1998, ASME V&V 20-2009, ANSI/AIAA S-071A-1999, JCGM 100(GUM)
다음 장: 1.6에서 제1부의 지도를 펼친다. 고정익·회전익·하이브리드 VTOL을 어떤 축으로 분류하고, 호버 효율과 순항 효율의 트레이드가 형상을 어떻게 구분하는지, 각 형상 장(1.7~1.15)이 이 지도 위 어디에 놓이는지를 본다.