와류격자법

전산유체역학 유체역학 마지막 수정: 2026-08-26 04:47:05

1. 개요[편집]

와류격자법 Vortex Lattice Method (VLM)
모형	두께 없는 캠버면 + 말굽와 격자
미지수	패널마다의 순환 Γ_ij
배치	결합와 패널 1/4 시위선 · 제어점 3/4 시위선
영향계수	비오-사바르 법칙의 유한 선분 해
선형계	조밀·비대칭 N×N (N ~ 수백)
대표 코드	AVL · Tornado · XFLR5 · VSPAERO

날개를 종이처럼 얇게 편 다음, 그 위에 말굽 모양 와류를 격자로 깔고 연립방정식 한 번 푼다. 노트북에서 0.1초.

와류격자법(vortex lattice method, VLM)은 두께를 무시한 양력면(캠버면)을 격자로 나누고 각 패널에 말굽와 하나씩을 얹은 뒤, 모든 제어점에서 유동이 면에 접한다는 조건을 부과해 패널 순환 $\Gamma_{ij}$ 를 미지수로 하는 선형 연립방정식을 푸는 3차원 저차 공력 수치기법이다. 계보로 보면 양력선 이론을 시위 방향으로도 이산화해 후퇴각·저종횡비·복수 양력면까지 다룰 수 있게 확장한 것이고, 패널법에서 두께를 덜어 낸 축소판이기도 하다.

역할 분담을 먼저 정해 두자. 쿠타 조건은 순환을 고정하는 물리 규칙, 얇은 익형 이론은 2차원 해석해, 와류격자법은 그 둘을 3차원에서 수치로 실행하는 기계다. 이 문서는 그 기계의 구조와 한계를 다룬다.

미지수가 수백 개면 충분해서 형상 하나를 푸는 데 걸리는 시간이 밀리초~초 단위다. 그래서 2020년대에도 개념설계실의 기본기로 남아 있다. 최적화 루프를 만 번 돌려야 하는 자리에서 RANS는 애초에 후보가 아니기 때문이다.

2. 모형 — 캠버면과 말굽와[편집]

날개를 시위 방향 $M$ 개, 스팬 방향 $N$ 개로 잘라 $M\times N$ 개의 사각 패널을 만든다. 각 패널에 말굽와(horseshoe vortex)를 하나씩 놓는데, 구조는 이렇다.

결합와(bound vortex): 패널의 1/4 시위선에 놓인 유한 길이 선분. 이 조각이 양력을 담당한다.
자유와(trailing legs): 결합와 양 끝에서 하류 무한대로 뻗는 두 개의 반무한 선분. 자유류 방향으로 흘려보내는 것이 기본(rigid wake)이다.

헬름홀츠 와정리가 “와선은 유체 속에서 끝날 수 없다”고 요구하므로 결합와는 반드시 어딘가로 이어져야 하고, 말굽 모양은 그 요구를 만족시키는 가장 싼 구조다. 고리를 저 멀리 하류에서 닫아 주는 것은 이륙 순간 남긴 시작와류이며, 그 뒷사정은 쿠타 조건 문서에 있다.

경계조건은 각 패널의 3/4 시위선 중앙에 잡은 제어점(control point)에서만 부과한다. 물체 표면이 아니라 캠버면 위에서, 그것도 유한개의 점에서만 거는 콜로케이션 방식이라 패널법보다 훨씬 성기다. 실용적으로는 시위 방향 5~~10개, 반스팬 방향 10~~20개면 공력계수가 1% 안으로 수렴하고, 앞전과 날개끝에는 코사인 간격을 써서 특이성이 강한 곳을 촘촘히 잡는 것이 국룰이다.

3. 왜 하필 1/4 과 3/4 인가[편집]

VLM 을 처음 배우는 사람이 가장 먼저 걸려 넘어지는 지점. 결합와와 제어점의 위치는 임의로 고른 것이 아니라 2차원 정답을 정확히 재현하도록 역산된 배치다. 평판 하나로 확인해 보자.

시위 $c$ 인 평판에 순환 $\Gamma$ 인 와류 하나를 $c/4$ 에 놓고, 제어점을 $3c/4$ 에 잡는다. 둘 사이 거리는 $c/2$ 이므로 제어점에서의 하향 유도속도는

w = \frac{\Gamma}{2\pi \cdot (c/2)} = \frac{\Gamma}{\pi c}

이고, 작은 받음각에서 접선 조건 $w = V_\infty \alpha$ 를 걸면 $\Gamma = \pi c V_\infty \alpha$ 다. 쿠타-주코프스키 정리를 통과시키면

c_l = \frac{\rho V_\infty \Gamma}{\tfrac12 \rho V_\infty^2 c} = 2\pi\alpha

얇은 익형 이론의 결과가 와류 한 개로 정확히 나온다. 게다가 결합와가 $c/4$ 에 있으니 그 점 둘레의 모멘트는 0, 즉 공력중심도 $c/4$ 로 맞는다. 양력 기울기와 공력중심이라는 2차원 익형 이론의 두 결과가 동시에 재현되는 배치는 이것뿐이며, 이 사실을 처음 지적한 사람의 이름을 따 피스톨레지 정리라 부른다.¹

이 배치의 부수 효과가 하나 더 있다. 제어점이 결합와보다 뒤에 있으므로 뒷전에서 와도가 따로 남지 않는 해가 자동으로 선택된다 — 즉 쿠타 조건이 기하학에 이미 구워져 있다. VLM 코드에 쿠타 조건 관련 코드가 한 줄도 없는 이유이자, 처음 읽는 사람이 그 부재에 당황하는 이유다.

패널을 시위 방향으로 여러 개 쓸 때도 규칙은 그대로 유지된다. 각 패널 안에서 1/4-3/4 비율을 지키면 캠버선의 시위 방향 하중 분포까지 근사되고, $M$ 을 늘릴수록 얇은 익형 이론의 $\gamma(x)$ 분포로 수렴한다.

4. 비오-사바르와 영향계수 행렬[편집]

말굽와 $j$ 가 세기 1일 때 제어점 $i$ 에 만드는 속도를 계산해야 한다. 도구는 비오-사바르 법칙의 유한 직선 와선 해다. 선분의 양 끝점에서 제어점까지의 벡터를 $\mathbf{r}_1, \mathbf{r}_2$ , 선분 벡터를 $\mathbf{r}_0$ 라 하면

\mathbf{V} = \frac{\Gamma}{4\pi}\, \frac{\mathbf{r}_1\times\mathbf{r}_2}{|\mathbf{r}_1\times\mathbf{r}_2|^2}\; \mathbf{r}_0\cdot\left(\frac{\mathbf{r}_1}{r_1} - \frac{\mathbf{r}_2}{r_2}\right)

반무한 선분은 한쪽 끝을 무한대로 보낸 극한을 쓴다. 닫힌 형태라 수치적분이 전혀 없다 — 패널법과 마찬가지로 행렬 조립이 대단히 싸다는 것이 저차 기법의 공통된 이점이다. 다만 제어점이 와선에 너무 가까우면 $|\mathbf{r}_1\times\mathbf{r}_2|^2$ 이 0으로 가며 유도속도가 발산하므로, 실무 코드는 반경 $r_c$ 이내에서 속도를 0으로 자르거나 랭킨 와류 코어를 씌운다.

말굽와 $j$ 가 제어점 $i$ 의 법선 방향에 만드는 성분을 $a_{ij}$ 라 두면 접선 조건은 곧바로 선형계가 된다.

\sum_{j=1}^{N} a_{ij}\,\Gamma_j = -\,\mathbf{V}_\infty\cdot\mathbf{n}_i

이 $\mathbf{A}$ 를 공력영향계수 행렬(AIC, aerodynamic influence coefficient matrix)이라 부른다. 조밀·비대칭이고, $N$ 이 수백 수준이라 LU 분해로 한 번에 푼다. 우변의 $\mathbf{n}_i$ 에 캠버 기울기·비틀림·플랩 편향·받음각·롤 각속도까지 전부 실을 수 있다는 점이 실무에서 결정적이다. 행렬은 형상에만 의존하므로 한 번 분해해 두면 비행조건 수천 개를 뒤에서 전진·후진 대입만으로 훑는다. 조종면 미계수 표를 초 단위로 뽑는 힘이 여기서 나온다.

5. 힘 — 쿠타-주코프스키와 트레프츠 평면[편집]

$\Gamma$ 를 얻었으면 힘은 두 갈래로 계산한다.

근접장(near-field). 각 결합와 조각에 쿠타-주코프스키 정리를 적용한다. 시위 방향 $i$ 번째, 스팬 방향 $j$ 번째 패널의 양력은

\Delta L_{ij} = \rho\, V_\infty \left(\Gamma_{ij} - \Gamma_{i-1,j}\right)\Delta y_j

시위 방향으로 앞 패널의 순환을 빼는 이유는, 겹쳐 놓인 말굽와들의 결합와 조각이 상쇄되고 남는 정미 와도만이 그 자리의 하중이기 때문이다. 이 값을 스팬 방향으로 모으면 $\Gamma(y)$ 분포와 국소 양력 분포가 그대로 나온다. 구조 하중 해석에 넘겨야 하는 것이 정확히 이 분포다.

원격장(far-field). 유도항력은 근접장에서 계산하면 이산화 오차에 민감하다. 대신 하류 무한대의 트레프츠 평면에서 후류 시트가 남긴 운동에너지를 적분한다.

D_i = \frac{\rho}{2}\int_{-b/2}^{b/2} \Gamma(y)\, w_T(y)\, dy

여기서 $w_T$ 는 트레프츠 평면에서의 하강기류이고, 날개 위치에서의 하강기류의 정확히 두 배다(반무한 와선이 자기 시작점 평면에서는 절반만 유도하기 때문). 이 적분은 후류의 스팬 방향 분포만 보므로 시위 방향 격자에 둔감하고, 그래서 근접장 합보다 훨씬 빨리 수렴한다. 계산 결과를 관례대로 정리하면

C_{D,i} = \frac{C_L^2}{\pi\, e\, AR}

이고, 스팬효율 $e$ 가 VLM 이 내놓는 가장 값진 숫자 중 하나다. 자세한 유도와 $e$ 의 의미는 양력선 이론 문서를 참고.

6. 스팬 순환 분포와 종횡비[편집]

VLM 결과를 볼 때 사람이 실제로 들여다보는 그림은 계수 하나가 아니라 $\Gamma(y)$ 곡선이다. 이 곡선 하나에 설계 정보가 다 들어 있다.

곡선이 타원에 가까울수록 $e \to 1$ , 유도항력이 최소. 사다리꼴 평면형에 워시아웃을 주는 것은 타원 평면형을 제작하지 않고도 이 곡선을 타원에 붙이려는 시도다.
곡선을 국소 시위로 나눈 $c_l(y) = 2\Gamma(y)/(V_\infty c(y))$ 가 어디서 먼저 실속하는가를 알려 준다. 이 값이 날개끝에서 먼저 최대에 닿으면 끝단 실속 → 에일러론이 먼저 죽는 최악의 시나리오다. 워시아웃 각도는 사실상 이 곡선을 보고 정한다.
곡선 아래 면적이 $C_L$ , 곡선의 1차 모멘트가 뿌리 굽힘 모멘트. 공력과 구조가 만나는 지점이 이 그림이다.

종횡비 $AR = b^2/S$ 의 효과도 여기서 나온다. 후류가 만드는 하강기류가 유효 받음각을 깎으므로 3차원 양력 기울기는 2차원 값보다 항상 작다. 타원 하중일 때

C_{L_\alpha} = \frac{a_0}{1 + a_0/(\pi\,AR)} = \frac{2\pi\,AR}{AR + 2} \qquad (a_0 = 2\pi)

$AR = 6$ 이면 2차원 값의 75%, $AR=2$ 면 50%다. 다만 이 식은 높은 종횡비 근사이므로 $AR \lesssim 4$ 에서는 헬름볼트 보정

C_{L_\alpha} = \frac{2\pi\,AR}{2 + \sqrt{AR^2+4}}

쪽이 실측에 가깝다. VLM 의 미덕은 이런 보정식을 외울 필요 없이 격자만 깔면 저종횡비·후퇴각·복합 평면형에서도 같은 답을 직접 계산해 준다는 것이다. 후퇴각이 붙는 순간 양력선 이론의 전제가 깨지는 것과 대비된다.

얇은 날개 근사와 평면 강체 후류로 말굽와류 96개(스팬 24 × 시위 4)를 깔고, 96×96 영향계수 행렬을 갱신마다 부분피벗 가우스 소거로 풀어 스팬방향 순환 Γ(y) 를 얻는다. 재생하면 종횡비를 2↔12 로 훑으며 매 프레임 계를 새로 푼다 — 직사각 날개 AR = 8 에서 CL_α = 4.6802/rad(헬름볼트 근사 대비 −4.6 %), 스팬효율 e = 0.9494 이고, 테이퍼를 λ ≈ 0.34 로 주면 같은 AR 에서 e = 0.9881 까지 올라간다. 타원 평면형에서는 e 가 0.99937 로 이론값 1 에서 1e−3 안이다.

7. 무엇을 못 하는가[편집]

가정을 되짚으면 한계가 그대로 나온다.

두께 없음. 두께가 만드는 표면 압력 분포, 동체 간섭, 흡입구 유동은 표현 불가. 그런 것이 필요하면 패널법으로 올라가야 한다.
점성 없음. 마찰항력·박리·실속이 전부 없다. $\alpha$ 를 40°로 넣어도 양력이 선형으로 자란다. 실무에서는 2차원 익형 데이터를 스트립별로 덧씌우는 준정상 보정(스트립 이론)을 얹어 실속을 흉내 내지만, 이는 VLM 이 아니라 그 위에 바른 반경험 모형이다.
선형화된 경계조건. 접선 조건을 평면(또는 얕은 캠버면)에 부과하므로 받음각이 크면 무너진다. 실용 범위는 대체로 부착유동 영역, $\alpha \lesssim 10^\circ$ .
강체 평면 후류. 후류를 자유류 방향 직선으로 고정한 기본형은 델타익의 앞전 와류나 헬리콥터 로터처럼 후류가 스스로 말리는 문제를 못 푼다. 후류 격자점을 계산된 속도장을 따라 옮기며 반복하는 후류 이완(free-wake)을 붙이면 다룰 수 있지만, 그 순간 문제가 비선형이 되고 계산량이 자릿수로 뛴다.
비압축. 아음속은 프란틀-글라우어트 보정으로 좌표를 늘려 얹지만 천음속은 불가.
앞전 흡입력 가정. 유도항력 계산은 앞전에서 흡입력이 온전히 발생한다고 전제한다. 앞전이 날카롭거나 큰 받음각에서 앞전 박리가 일어나면 그 흡입력이 실현되지 않아 항력이 과소평가된다. 델타익처럼 앞전 와류가 오히려 양력을 만드는 형상에서는 폴하무스(Polhamus)의 앞전 흡입 유추로 그 간극을 반경험적으로 메운다.²

8. 그래서 언제 쓰나[편집]

상황	VLM	패널법	RANS CFD
개념설계 형상 훑기	최적	가능	불가(느림)
안정·조종 미계수 표	최적	가능	과잉
스팬 하중 → 구조 입력	최적	가능	가능
정적 공탄성(비틀림 발산·조종면 반전)	최적	가능	과잉
두께·동체 간섭	불가	가능	가능
박리·실속·충격파	불가	불가	가능

정적 공탄성이 VLM 의 숨은 주력 무대다. 날개가 공력하중으로 비틀리면 받음각이 변하고, 변한 받음각이 다시 하중을 바꾸는 되먹임 문제인데, AIC 행렬 $\mathbf{A}$ 와 구조 유연도 행렬을 곱해 고유값 하나만 보면 비틀림 발산 동압이 나온다. 조종면 반전 속도도 같은 틀에서 떨어진다. 이 계산에 CFD 를 쓰는 것은 낭비이고, 실제로 항공탄성학 실무는 이 계열의 저차 모형 위에 서 있다. 주파수 영역으로 확장한 사촌이 플러터 해석의 산업 표준인 더블릿-격자법(DLM)이다.

대표 구현으로는 마크 드렐라와 하럴드 영그렌의 AVL 이 사실상 표준이다. 텍스트 파일 하나로 전기체 형상과 조종면을 정의하면 6자유도 공력미계수를 통째로 뱉어 주며, 소형 무인기·모형항공기·대학 설계 경진대회 판을 30년째 장악하고 있다. 오픈소스 계열로 Tornado(MATLAB), XFLR5(XFOIL 2차원 데이터와 VLM 을 결합), OpenVSP 의 VSPAERO 등이 있다.

정리하면 이렇다. 와류격자법은 답을 확정하는 도구가 아니라 설계 공간을 훑는 도구다. 스팬과 비틀림과 꼬리날개 크기를 정하는 것은 여전히 이쪽 몫이고, CFD 는 그렇게 정해진 형상 하나를 나중에 검증한다. 값싼 저차 모형이 없으면 최적화 루프는 시작조차 못 한다는 사실이 60년 넘게 변하지 않았다.³⁴

9. 관련 문서[편집]

이탈리아의 엔리코 피스톨레지가 1930년대에 정리한 결과라 그 이름이 붙었지만, 영어권에서는 그냥 “1/4-3/4 rule” 이라고 부르는 쪽이 압도적으로 흔하다. 정리 이름을 아는 것보다 중요한 건 이게 우연이 아니라는 감각이다. 격자를 대충 깔고 결합와를 패널 중앙에 놓으면 $c_l/\alpha$ 가 $2\pi$ 에서 눈에 띄게 벗어나며, 그때 대부분의 사람은 격자를 더 쪼개는 잘못된 처방으로 시간을 태운다. ↩
앞전 흡입 유추(leading-edge suction analogy)는 “앞전에서 실현되지 못한 흡입력의 크기만큼이 와류 양력으로 방향만 바뀌어 나타난다”고 놓는 발상이다. 유도라기보다 에너지 장부를 억지로 맞춘 것에 가까운데, 델타익 실험 데이터와 놀랄 만큼 잘 맞아서 60년째 쓰이고 있다. 이 바닥에서 “왜 맞는지 모르지만 맞는” 모형이 살아남는 전형적인 사례. ↩
VLM 이 초기 설계에서 죽지 않는 진짜 이유는 속도만이 아니라 미분 가능성이다. AIC 행렬이 형상의 매끄러운 함수라 인접법(adjoint)이나 유한차분으로 설계 민감도를 안정적으로 뽑을 수 있다. 격자를 다시 만들 때마다 목적함수가 미세하게 튀는 CFD 최적화에 데어 본 사람이라면 이 성질의 가치를 안다. ↩
물론 값싼 도구의 고전적 함정도 그대로다 — 결과가 예쁘게 나온다는 것. 실속도 박리도 없으니 곡선이 항상 매끈하게 선형이고, 어디서도 경고를 뱉지 않는다. 학부 설계 프로젝트에서 받음각 25°짜리 VLM 결과를 자랑스럽게 발표하는 장면이 매년 반복되는 이유다. ↩