패널법

전산유체역학 유체역학 마지막 수정: 2026-08-24 04:23:07

1. 개요[편집]

패널법 Panel method
분류	경계요소법의 포텐셜 유동 특화 버전
지배 방정식	라플라스 방정식 ∇²φ = 0
미지수	표면 패널마다의 소스·더블릿·와류 세기
경계조건	무침투(∂φ/∂n = 0) + 쿠타 조건
선형계	조밀·비대칭 N×N
원조	헤스 & 스미스 (Douglas Aircraft, 1962)

유동장을 안 푼다. 표면만 푼다. 그것도 노트북에서 0.01초에.

패널법(panel method)은 포텐셜 유동의 지배식인 라플라스 방정식을 경계적분방정식으로 바꾼 뒤, 물체 표면을 유한 개의 패널로 나누고 각 패널에 소스·더블릿·와류 같은 특이점을 분포시켜 그 세기를 미지수로 푸는 수치기법이다. 즉 경계요소법을 비점성·비회전 유동에 특화시킨 물건이며, 실제로 항공·조선 분야에서 BEM이 산업 코드로 자리 잡은 것은 거의 전부 이 이름으로였다.

핵심 이득은 차원 하나가 사라진다는 것. 3차원 항공기 주위 공기를 격자로 채우는 대신 표면만 삼각형·사각형으로 덮으면 되고, 원거리 경계조건은 특이점 해가 무한대에서 자동으로 균일류에 수렴하므로 신경 쓸 필요가 없다. 미지수가 수천 개면 충분해서, RANS 한 케이스 돌릴 시간에 패널법은 수만 케이스를 돌린다.

물론 공짜는 아니다. 점성이 없으니 마찰항력도 박리도 실속도 없다. 그럼에도 2020년대의 초기 설계실에서 패널법이 여전히 켜져 있는 이유는, 형상 최적화 루프처럼 “정확도보다 회수(回數)가 중요한” 작업에서 대안이 없기 때문이다.¹

2. 왜 표면만으로 충분한가[편집]

$\nabla^2\phi=0$ 의 그린 함수(3차원 기본해 $-1/4\pi r$ )를 이용해 그린 제3항등식을 쓰면, 영역 내부의 퍼텐셜이 경계 위의 값과 법선 도함수만으로 재구성된다.

\phi(\mathbf{x}) = \phi_\infty(\mathbf{x}) + \frac{1}{4\pi}\oint_S \left[ \sigma\,\frac{1}{r} - \mu\,\frac{\partial}{\partial n}\frac{1}{r} \right] dS

여기서 $\sigma$ 는 소스 세기(source strength, 법선 속도 점프에 대응), $\mu$ 는 더블릿 세기(doublet strength, 퍼텐셜 점프에 대응)다. 이 식이 패널법의 전부다. 물체를 어떻게 생겼든 관계없이, 표면에 이 두 종류의 특이점을 적당한 세기로 발라 놓기만 하면 라플라스 방정식은 이미 자동으로 만족되고, 남은 일은 경계조건을 맞추는 것뿐이다.

소스만 쓰면 물체를 “뚫고 나오는 유량”을 조절할 수 있다. 폐곡면에 발라 두께를 표현한다.
더블릿은 퍼텐셜 점프, 즉 와도를 표면에 실은 것과 동치다. 2차원에서 세기 $\mu$ 의 더블릿 분포는 세기 $\gamma=-d\mu/ds$ 의 와류 분포와 정확히 같다. 순환을 만들어야 하는 양력 문제에는 이쪽이 필수다.

경계조건은 두 가지 방식으로 건다. 노이만 형식은 표면에서 법선 속도를 0으로 두고( $\mathbf{V}\cdot\mathbf{n}=0$ , 무침투), 디리클레 형식은 물체 내부의 퍼텐셜을 상수로 고정한다(내부에 유동이 없다는 뜻이라 등가다). 실무 3차원 코드는 수치적으로 더 얌전한 디리클레 쪽을 자주 택한다.

3. 이산화 — 콜로케이션과 영향계수[편집]

절차는 경계요소법의 그것과 같지만, 이 바닥에서는 용어가 좀 다르다.

패널 분할. 표면을 $N$ 개의 평면 조각으로 근사한다. 2차원이면 익형 윤곽선을 꺾은선으로, 3차원이면 사각형 패널로.
세기 분포 가정. 패널마다 세기를 상수(constant-strength)로 두는 것이 가장 흔하고, 정확도를 원하면 선형(linear-strength)·이차로 올린다.
콜로케이션. 각 패널의 중앙점(control point)을 하나씩 잡아 거기서 경계조건이 정확히 만족되도록 요구한다. 요구조건 $N$ 개, 미지수 $N$ 개.
영향계수 조립. 패널 $j$ 가 세기 1일 때 콜로케이션점 $i$ 에 유도하는 속도(또는 퍼텐셜) $a_{ij}$ 를 계산한다.
선형계 풀이. $\mathbf{A}\boldsymbol{\sigma} = \mathbf{b}$ , 우변은 균일류의 기여.

여기서 패널법이 유한요소·유한체적과 결정적으로 다른 점 하나. 영향계수를 수치적분하지 않는다. 평면 패널 위에 상수 세기를 두면 적분이 손으로 풀린다. 2차원 상수 소스 패널의 유도 퍼텐셜은 패널 국소좌표에서

\phi = \frac{\sigma}{2\pi}\int_0^{\ell} \ln\sqrt{(x-s)^2+y^2}\; ds

인데, 이 적분이 $t\ln r$ 꼴의 로그항과 $y\arctan(t/y)$ 꼴의 아크탄젠트항의 조합으로 손으로 떨어진다. 3차원 상수 소스·더블릿 사각패널도 헤스가 1962년에 정리한 해석식이 있고 구조는 같다(로그항 + 아크탄젠트항). 그래서 행렬 조립이 대단히 빠르고, 경계요소법 일반에서 대학원생을 괴롭히는 특이적분 문제도 패널 자기 자신에 대해서는 해석적 극한값( $\sigma/2$ 같은 자유항)으로 깔끔하게 정리된다.

대가는 행렬 구조다. 모든 패널이 모든 패널에 영향을 주므로 $\mathbf{A}$ 는 꽉 찬 비대칭 행렬이고, LU 분해로 정공법을 쓰면 $O(N^3)$ · $O(N^2)$ 메모리다. $N\sim3000$ 까지는 그냥 풀고, 그 이상은 GMRES + 고속 다중극자법으로 행렬-벡터 곱을 $O(N\log N)$ 에 처리한다. 패널이 아주 얇게 붙거나 종횡비가 극단적이면 조건수가 폭발하므로, 패널 분할 품질이 곧 해의 품질이라는 점은 격자 기반 CFD와 다르지 않다.

4. 소스 패널법 — 양력 없는 물체[편집]

비양력 물체(원기둥, 동체, 잠수함 선체, 자동차 형상)는 소스만으로 끝난다. 무침투 조건을 콜로케이션점 $i$ 에서 쓰면

\frac{\sigma_i}{2} + \sum_{j\ne i} \sigma_j\, a_{ij} + V_\infty \cos\beta_i = 0

이고, 좌변 첫 항 $\sigma_i/2$ 가 자기 패널의 자유항이다. 이걸 풀면 표면 접선속도 $V_t$ 가 나오고, 베르누이로 압력계수를 얻는다.

C_p = 1 - \left(\frac{V_t}{V_\infty}\right)^2

원기둥에 적용하면 정확해 $C_p = 1-4\sin^2\theta$ 를 패널 수십 개만으로 소수점 셋째 자리까지 재현한다. 그리고 앞뒤가 완벽히 대칭인 $C_p$ 분포가 나오므로 압력을 적분한 항력은 정확히 0 — 달랑베르 역설이 수치 결과로 눈앞에 찍힌다. 자세한 사연은 포텐셜 유동 문서 참고.

여기까지가 1962년 헤스와 스미스가 더글러스 항공에서 만든 원조 코드다. 임의 형상 3차원 물체를 표면 이산화만으로 풀 수 있다는 것이 당시로서는 혁명이었고, 이후 20년간 항공산업의 공력 해석은 사실상 이 계열이 담당했다.

5. 와류 패널법과 쿠타 조건 — 양력이 생기는 자리[편집]

익형에 소스만 바르면 양력이 0이다. 순환 $\Gamma$ 가 없으면 쿠타-주코프스키 정리에 의해 $L=\rho V_\infty \Gamma = 0$ 이기 때문. 그런데 순환을 얼마로 잡아야 하는지는 라플라스 방정식이 알려주지 않는다 — 비점성 문제의 해가 순환 값마다 하나씩, 무한히 많다.

이 부정성을 끊는 물리적 처방이 쿠타 조건이다. “실제 유동은 뾰족한 뒷전을 돌아 나가지 못한다. 따라서 뒷전에서 속도가 유한해야 한다.” 이산화 형태로는 뒷전 위·아래 패널의 와류 세기 합을 0으로 두거나( $\gamma_{\text{upper}} + \gamma_{\text{lower}} = 0$ ), 뒷전 위·아래의 접선속도 크기를 같게 놓는다.

고전적인 헤스-스미스 익형 패널법의 구성이 우아하다.

패널마다 독립적인 상수 소스 세기 $\sigma_j$ ( $N$ 개) — 두께를 담당.
전 패널에 공통인 상수 와류 세기 $\gamma$ (1개) — 순환을 담당.
방정식: 콜로케이션 $N$ 개 + 쿠타 조건 1개 = $N+1$ 개.

$(N+1)\times(N+1)$ 계를 한 번 풀면 끝이다. 익형 하나에 패널 100개면 밀리초 단위로 $C_p$ 분포가 나온다.

NACA 0012 표면을 N개 직선 패널로 쪼개, 패널마다 소스강도 하나와 익형 전체에 공통인 와류강도 하나를 미지수로 놓고 무침투 N식 + 쿠타 조건 1식의 (N+1)차 조밀 선형계를 부분 피벗 가우스 소거로 푼다. 영향계수는 직선 패널 위 균일 소스·와류의 해석 유도속도식이고, α=5°·N=80에서 c_l = 0.607이 나온다 — 얇은 익형 이론 2πα = 0.548보다 11% 큰 만큼이 두께 효과다. 점성이 없어 항력은 0이고(화면의 c_d = 5.6e−4는 이산화 잔차) 실속도 나타나지 않는다.

얇은 익형 이론과의 대조가 이 지점에서 재미있다. 얇은 익형 이론은 두께를 0으로, 받음각을 작다고 가정해 캠버선에만 와류를 발라 손으로 푸는 이론이고, 결과는 그 유명한

C_l = 2\pi\left(\alpha - \alpha_{L=0}\right)

이다. 패널법은 같은 물리를 두께까지 포함해 수치로 푸는 것이므로, 두 결과를 비교하면 두께가 하는 일이 정량적으로 드러난다. 실제로 두께비가 커질수록 양력곡선 기울기가 $2\pi$ 보다 조금 커지고(익형 표면이 유효 캠버선보다 길어 국소 유속이 빨라지는 효과), 반대로 실제 실험값은 경계층 배제두께가 뒷전을 두껍게 만드는 탓에 $2\pi$ 보다 작게 나온다. 비점성 패널법이 실험보다 양력을 과대예측하는 고전적 패턴이 여기서 나온다.

6. 3차원 — 양력면법과 와류격자법[편집]

3차원에서는 후류(wake)가 새로운 골칫거리로 등장한다. 날개 뒷전에서 흘러나간 와류면은 순환을 하류로 실어 나르며, 그 위치 자체가 해의 일부다.

와류격자법(Vortex Lattice Method, VLM). 두께를 무시하고 캠버면만 격자로 나눈 뒤, 각 패널의 1/4 시위선에 말굽와류를 놓고 3/4 시위선 중앙에 콜로케이션점을 잡는다. 이 1/4-3/4 배치가 2차원 얇은 익형 이론의 결과를 정확히 재현한다는 것이 피스톨레지(Pistolesi)의 정리다. 미지수가 수백 개면 충분해서 항공기 전기체 공력미계수를 초 단위로 뽑는다. AVL이 대표적.
전면 패널법(full panel method). 두께를 포함한 폐곡면에 소스+더블릿을 바른다. PANAIR(A502), VSAERO, PMARC 계열. 동체-날개 간섭이나 흡입구처럼 두께가 중요한 형상에 쓴다.
후류 이완(wake relaxation). 처음에는 후류를 자유류 방향 평면으로 고정해 놓고(rigid wake) 푼 뒤, 계산된 속도장을 따라 후류 격자점을 옮기고 다시 푸는 것을 반복한다. 델타익의 앞전 와류나 헬리콥터 로터처럼 후류가 자기 자신을 감아 올리는 문제에서는 이 비선형 반복이 필수다.

항력에 대한 오해를 여기서 정리해야 한다. “패널법은 항력을 못 준다”는 말은 반만 맞다. 3차원 양력 유동에서는 후류가 하류 무한대에 남기는 운동에너지가 있고, 이를 트레프츠 평면(Trefftz plane)에서 적분하면 유도항력이 정확히 나온다. VLM으로 계산한 $C_{D_i} = C_L^2/(\pi e AR)$ 은 신뢰할 만하며, 실제로 날개 평면형 최적화의 표준 도구다. 패널법이 원리적으로 못 주는 것은 점성 마찰항력과 박리에 의한 형상항력이다. 2차원 비양력 물체에서 항력이 정확히 0이 되는 달랑베르 역설은 그중 후자의 이야기다.

7. 점성과의 결합 — 그래서 XFOIL[편집]

비점성 해가 실험과 벌어지는 원인이 경계층 하나라면, 경계층만 따로 붙이면 되지 않을까? 이 발상이 점성-비점성 상호작용(viscous-inviscid interaction)이다.

패널법으로 표면 속도 분포 $u_e(s)$ 를 얻는다.
그 $u_e$ 를 입력으로 적분형 경계층 방정식(운동량 적분식 + 형상계수 방정식)을 풀어 배제두께 $\delta^*(s)$ 를 얻는다.
물체 표면을 $\delta^*$ 만큼 부풀리거나(변위 두께 개념), 표면에 등가 소스 $\sigma = d(u_e\delta^*)/ds$ 를 추가로 분포시킨다.
다시 패널법을 푼다. 수렴할 때까지 반복.

이 아이디어를 가장 완성도 높게 구현한 것이 마크 드렐라의 XFOIL(1986)이다. 선형 와류 패널법 + 2방정식 적분 경계층 + $e^N$ 천이 예측을 하나의 뉴턴 반복으로 완전 결합해 푸는 구조라, 약한 박리 버블까지 통과해 수렴한다. 순차 반복 방식이 박리 근처에서 발산하는(골드슈타인 특이점) 고질병을 정면 돌파한 설계다. 천이 판정의 $N_{\text{crit}}$ 기본값 9는 저난류 풍동에 대응하는 값으로, 자유대기나 거친 표면에서는 낮춰 잡아야 한다.²

덕분에 XFOIL은 저레이놀즈수 익형 설계(모형항공기, 소형 UAV, 풍력 블레이드 뿌리 단면)에서 40년째 사실상의 표준이다. 실속 근처처럼 박리가 커지면 당연히 무너지지만, 부착 유동 영역에서의 $C_l$ · $C_d$ 예측 정확도는 난류 모델링을 얹은 RANS와 비교해도 밀리지 않으면서 계산 시간이 만 배 싸다.

8. 아직 살아 있는 이유[편집]

전산유체역학이 나비에-스토크스 방정식을 통째로 푸는 시대에 왜 60년 된 비점성 코드가 남아 있는가. 답은 비용 구조다.

항목	패널법	RANS
격자 준비	표면만, 분 단위	체적 격자, 시간~일 단위
1케이스 계산	밀리초~초	분~시간
형상 변경 대응	표면 재분할만	격자 재생성
유도항력	정확	정확
마찰·형상항력	불가 (결합 필요)	가능
충격파·박리	불가	가능

수천 번의 형상 평가가 필요한 최적설계 루프, 수만 개 비행조건을 훑어야 하는 공력 데이터베이스 구축, 실시간 비행 시뮬레이터의 공력 모델, 플러터 해석의 더블릿-격자법(doublet-lattice, 비정상 패널법의 주파수 영역 버전) — 전부 “빠르고 미분 가능한 근사”가 정확한 한 점보다 값진 곳들이다. 선박 쪽에서도 조파저항 계산과 프로펠러 성능 예측에 패널법이 여전히 현역이다.

요약하면 이렇다. 패널법은 최종 검증용 도구가 아니다. 설계 공간을 훑는 도구이고, 그 역할에서는 아직 대체재가 없다. 최근에는 패널법으로 만든 대량 데이터가 서로게이트 모델의 학습 데이터로 재활용되면서 오히려 수요가 늘었다는 이야기도 나온다.³

9. 관련 문서[편집]

패널법이 등장하기 전 익형 해석의 정석은 등각사상이었다. 원기둥의 해를 알고 있으니 익형을 원기둥으로 사상하면 된다는 발상인데, 문제는 “임의의 형상”을 사상하는 함수를 매번 찾아야 한다는 것. 조코프스키 익형이 실용성보다 수학적 편의 때문에 유명해진 이유이기도 하다. 패널법은 이 제약을 통째로 없앴다 — 형상이 뭐든 그냥 꺾은선으로 근사하면 되니까. ↩
$e^N$ 방법의 $N$ 은 톨미엔-슐리히팅 파의 증폭 인자 $\ln(A/A_0)$ 이고, $N=9$ 는 $e^9\approx8100$ 배 증폭에서 천이가 일어난다는 경험칙이다. 물리적 유도가 아니라 풍동 데이터에 맞춘 값이라, “ $N_{\text{crit}}$ 을 몇으로 놨느냐”가 곧 “어느 풍동을 믿느냐”와 같은 말이 된다. 논문에서 XFOIL 결과를 볼 때 이 숫자를 안 적어 놨으면 일단 의심하는 게 안전하다. ↩
물론 그 서로게이트가 비점성 데이터로 배웠다는 사실은 어디 안 간다. 패널법이 못 보는 박리를 신경망이 대신 봐 줄 리는 없다. “데이터가 많으면 물리가 따라온다”는 믿음이 가장 조용히 배신하는 지점이 이런 곳이다. ↩

패널법