경계요소법

1. 개요

경계요소법 Boundary Element Method
약칭	BEM
분야	수치해석 × 적분방정식 × 포텐셜 이론
핵심 아이디어	경계에만 요소를 놓아 차원을 하나 낮춤 (3D→2D)
기반 도구	그린 함수, 기본해, 경계적분방정식
강점 영역	무한·반무한 영역 (음향, 전자기, 포텐셜 유동, 응력)

경계요소법(Boundary Element Method, BEM)은 지배 편미분방정식을 경계적분방정식(boundary integral equation)으로 변환한 뒤 영역 내부가 아니라 경계에만 요소를 배치하여 푸는 수치해석 기법이다. 유한요소법이나 유한차분법이 3차원 물체 내부를 통째로 격자로 채워야 하는 것과 달리, BEM은 3차원 문제를 그 표면(2차원)에서만, 2차원 문제를 그 윤곽선(1차원)에서만 이산화한다. 한마디로 차원을 하나 깎아내는 마법이다.¹

이 마법이 공짜는 아니다. 대가로 얻는 것은 조밀(dense)하고 비대칭(non-symmetric)인 행렬. 하지만 무한히 뻗어나가는 영역을 다뤄야 하는 음향 방사, 전자기 산란, 포텐셜 유동 같은 문제에서 BEM은 유한요소법이 흉내 낼 수 없는 우아함을 보여준다.

2. 왜 경계만 풀 수 있나

핵심은 그린 함수(Green’s function)와 기본해(fundamental solution)다. 선형 미분연산자 $\mathcal{L}$에 대해, 점 소스에 대한 응답인 기본해 $u^*$는 다음을 만족한다.

$$\mathcal{L}\, u^*(\mathbf{x}, \boldsymbol{\xi}) = -\delta(\mathbf{x} - \boldsymbol{\xi})$$

여기서 $\delta$는 디랙 델타, $\boldsymbol{\xi}$는 소스점이다. 예를 들어 라플라스 방정식 $\nabla^2 u = 0$의 3차원 기본해는 $u^* = 1/(4\pi r)$로, 점전하가 만드는 퍼텐셜 그 자체다. 이 기본해가 지배방정식을 이미 만족하기 때문에, 그린의 제2항등식(Green’s second identity)과 가중잔차법을 적용하면 영역 적분이 마법처럼 사라지고 경계적분방정식만 남는다.

$$c(\boldsymbol{\xi})\, u(\boldsymbol{\xi}) + \int_{\Gamma} u\, q^*\, d\Gamma = \int_{\Gamma} q\, u^*\, d\Gamma$$

여기서 $u$는 경계에서의 미지 함수(예: 퍼텐셜), $q = \partial u / \partial n$은 그 법선 도함수(예: 유속), $c$는 경계의 기하학적 각도에 따른 자유항 계수다. 미지수가 오로지 경계 위에만 있다는 점이 이 식의 전부다. 내부는 계산이 끝난 뒤 필요한 지점에서만 후처리로 뽑아낸다.²

3. 이산화 절차

경계적분방정식을 실제로 푸는 과정은 유한요소법과 형식은 닮았지만 무대가 경계라는 점이 다르다.

• 경계 이산화: 표면을 삼각형·사각형 패널(3D) 또는 선분(2D)으로 나눈다. 각 요소 위에서 $u$와 $q$를 기저함수로 보간한다(상수/선형/이차 요소).
• 선점 배치(collocation): 각 노드를 소스점 $\boldsymbol{\xi}$로 삼아 경계적분방정식을 세운다. 노드 수만큼 방정식이 생긴다.
• 적분 계산: 각 요소에 대해 $\int q^*\, d\Gamma$와 $\int u^*\, d\Gamma$를 수치적분한다. 소스점이 요소 위에 있을 때는 적분핵이 특이점(singularity)이 되어 특수한 처리가 필요하다.³
• 선형계 조립·해석: 결과는 $\mathbf{H}\mathbf{u} = \mathbf{G}\mathbf{q}$ 꼴의 선형계다. 경계조건(주어진 $u$ 또는 $q$)을 대입해 미지수만 좌변으로 옮기면 $\mathbf{A}\mathbf{x} = \mathbf{b}$가 완성된다.

4. 유한요소법과의 비교

BEM과 유한요소법은 라이벌이라기보다 상호보완적이다. 각자 잘하는 영역이 뚜렷하게 갈린다.

항목	경계요소법 (BEM)	유한요소법 (FEM)
이산화 대상	경계만 (차원 −1)	영역 전체
미지수 개수	적음	많음
행렬 구조	조밀·비대칭	희소·대칭
무한 영역	자연스럽게 처리	인위적 절단 경계 필요
비선형·비균질	불리 (영역 적분 부활)	유리

BEM의 행렬은 꽉 찬 $N \times N$이라 정공법으로 풀면 $O(N^3)$ 비용에 $O(N^2)$ 메모리가 든다. 노드 수가 수만을 넘어가면 이게 발목을 잡는다. 반면 FEM 행렬은 희소행렬이라 반복법으로 값싸게 풀린다. 그래서 “적은 미지수 vs 값싼 행렬”의 트레이드오프가 성립한다.⁴

5. 어디에 쓰나

BEM이 진가를 발휘하는 곳은 무한·반무한 영역 문제다. 절단 경계를 어디에 둘지 고민할 필요가 없다는 게 결정적이다.

• 포텐셜 유동: 비점성·비회전 유동의 속도 퍼텐셜을 패널법(panel method)으로 푼다. 항공기 초기 설계에서 날개 주위 유동을 빠르게 추정하는 고전적 도구.
• 음향·소음 방사: 헬름홀츠 방정식으로 기술되는 외부 음장. 스피커·엔진 소음이 무한 공간으로 퍼져나가는 문제에 이상적이다.
• 안테나 해석·전자기 산란: 도체 표면의 전류 분포를 미지수로 두는 모멘트법(Method of Moments)이 사실상 BEM의 전자기 버전이다.
• 탄성·응력집중: 균열 선단이나 표면 응력을 경계에서 직접 다뤄 파괴역학에서 유용하다.
• 역문제: 미지 형상·소스 추정처럼 경계 정보만으로 내부를 역추정하는 문제와 궁합이 좋다.

6. 현실에서의 BEM

BEM은 “만능 대체재”는 아니다. 오히려 특정 문제군의 스페셜리스트에 가깝다. 조밀 행렬 문제를 완화하기 위해 고속 다중극법(Fast Multipole Method, FMM)이나 계층 행렬(H-matrix) 같은 가속 기법이 개발되어, 이제는 수백만 미지수 규모도 다룰 수 있다.⁵ 그럼에도 비선형·비균질·시간의존 문제 앞에서는 영역 적분이 부활하며 “차원 하나 깎는 마법”이 풀려버린다. 결국 도구는 문제에 맞춰 고르는 것이고, 무한 영역의 선형 문제라면 BEM만한 우아함이 없다는 게 이 바닥의 국룰이다.

7. 관련 문서

• 유한요소법 · 유한차분법 · 유한체적법
• 포텐셜 유동
• 안테나 해석 · 모멘트법
• 라플라스 방정식 · 포아송 방정식
• 기저함수 · 수치적분
• 희소행렬 · 반복법
• 파괴역학 · 응력집중

8. Footnotes

1. 물론 “하나 낮춘다”는 표현은 이산화 차원 이야기다. 물리 문제 자체가 3차원이면 답도 3차원이다. BEM은 그 3차원 답을 경계 정보만으로 재구성할 수 있다는 뜻이지, 세상이 2차원이 되는 건 아니다. 착각하면 곤란하다. ↩

2. 그래서 BEM은 “내부의 특정 몇 점에서만 값이 필요한” 문제에 특히 경제적이다. 영역 전체를 채운 뒤 필요 없는 절점까지 다 계산하는 FEM과 대조적이다. ↩

3. 특이적분(singular integral) 처리는 BEM 구현의 난관 중 하나다. 코시 주치(Cauchy principal value), 정칙화(regularization), 해석적 적분 공식 등 온갖 기법이 동원된다. “경계만 풀면 되니 쉽겠지” 하고 덤볐다가 이 특이점에서 무릎 꿇는 대학원생이 매년 나온다. ↩

4. 정리하면 이렇다. BEM은 미지수는 적지만 행렬이 꽉 차 있고, FEM은 미지수는 많지만 행렬이 텅 비어 있다. 세상에 공짜 점심은 없다. ↩

5. FMM은 그린가드·로클린(Greengard & Rokhlin, 1987)이 다체 문제 가속을 위해 개발한 알고리즘으로, 20세기 10대 알고리즘에 꼽힌다. BEM의 $O(N^2)$ 병목을 $O(N \log N)$ 수준으로 끌어내려 BEM을 실무 규모로 되살렸다. ↩