그린 함수

수치해석 물리 마지막 수정: 2026-07-14 04:26:48

1. 개요[편집]

그린 함수 Green's Function
분야	편미분방정식 · 수리물리 · 적분방정식
고안	조지 그린 (1828)
정의	선형 미분연산자의 임펄스 응답
입력	점 소스 (디랙 델타)
주 사용처	경계요소법, 전자기해석, 구조해석

방정식에게 “점 하나 콕 찍으면 어떻게 반응하니?”라고 물어보면, 나머지는 적분이 다 한다.

그린 함수(Green’s function)는 선형 미분연산자에 점 형태의 소스(임펄스)를 넣었을 때 나오는 응답 함수다. 신호처리를 아는 사람이라면 “아, 임펄스 응답(impulse response)이구나” 하면 정확하다. 선형 시불변 시스템의 출력을 입력과 임펄스 응답의 합성곱으로 구하듯, 선형 편미분방정식의 해도 소스 항과 그린 함수의 합성곱(적분)으로 구할 수 있다.

수학적으로, 선형 미분연산자 $\mathcal{L}$ 에 대한 그린 함수 $G(\mathbf{x}, \mathbf{x}')$ 는 다음을 만족한다.

\mathcal{L}\, G(\mathbf{x}, \mathbf{x}') = \delta(\mathbf{x} - \mathbf{x}')

여기서 $\delta$ 는 디랙 델타 함수, 즉 $\mathbf{x}'$ 에 몰빵된 무한히 뾰족한 점 소스다. 이 하나만 알면 임의의 소스 $f(\mathbf{x})$ 에 대한 방정식 $\mathcal{L}u = f$ 의 해가 단번에 적분으로 튀어나온다. 조지 그린이라는, 정규 교육이라곤 1년밖에 못 받은 방앗간 집 아들이¹ 1828년에 자비 출판한 소논문에서 이 개념을 처음 던졌다는 걸 생각하면 더 경이롭다.

2. 임펄스 응답으로서의 그린 함수[편집]

왜 점 소스 하나의 응답만 알면 모든 게 풀릴까? 핵심은 선형성과 중첩 원리다. 임의의 소스 분포 $f(\mathbf{x})$ 는 무수히 많은 점 소스의 가중합으로 볼 수 있다.

f(\mathbf{x}) = \int f(\mathbf{x}')\, \delta(\mathbf{x} - \mathbf{x}')\, d\mathbf{x}'

각 점 소스 $\delta(\mathbf{x} - \mathbf{x}')$ 의 응답이 $G(\mathbf{x}, \mathbf{x}')$ 이고, 연산자가 선형이니 전체 응답은 그냥 다 더하면(적분하면) 된다.

u(\mathbf{x}) = \int G(\mathbf{x}, \mathbf{x}')\, f(\mathbf{x}')\, d\mathbf{x}'

미분방정식을 푸는 문제가 적분 계산으로 바뀌었다. 이 관점은 스펙트럴 방법이나 변분법과는 결이 다른, “응답을 미리 알아 두고 갖다 붙이는” 접근이다. 물론 공짜는 아니다. 그린 함수 자체를 구하는 게 원래 방정식을 푸는 것만큼 어려울 수 있고, 경계가 복잡하면 닫힌 형태로는 거의 안 나온다.

3. 자유공간 그린 함수[편집]

경계가 무한히 멀리 있는(자유공간, free-space) 이상적인 경우에는 그린 함수가 깔끔한 닫힌 형태로 존재한다. 물리에서 가장 유명한 예가 라플라스 방정식·포아송 방정식을 지배하는 연산자 $\mathcal{L} = -\nabla^2$ 의 그린 함수다. 3차원에서는

G(\mathbf{x}, \mathbf{x}') = \frac{1}{4\pi\, |\mathbf{x} - \mathbf{x}'|}

이것을 보고 “어? 쿨롱 퍼텐셜이네?” 했다면 정확하다. 점전하가 만드는 $1/r$ 퍼텐셜이 바로 라플라시안의 자유공간 그린 함수인 것이다. 즉 물리 시간에 외웠던 그 퍼텐셜 공식은, 알고 보면 그린 함수였다.² 2차원에서는 로그 형태 $-\frac{1}{2\pi}\ln|\mathbf{x}-\mathbf{x}'|$ 가 되고, 시간이 들어간 파동방정식이면 지연 시간이 반영된 지연 그린 함수(retarded Green’s function)가 된다. 헬름홀츠 방정식(진동수 영역의 파동)은 $\frac{e^{ikr}}{4\pi r}$ 꼴이라, 파동이 소스에서 퍼져 나가는 모습이 지수항에 그대로 담겨 있다.

4. 경계요소법에서의 역할[편집]

그린 함수가 순수 이론을 넘어 밥벌이가 되는 대표 무대가 경계요소법(Boundary Element Method, BEM)이다. 발상은 이렇다. 지배 방정식의 그린 함수를 이미 알고 있다면, 그린 항등식(Green’s identity)을 써서 도메인 내부 전체의 방정식을 경계 위의 적분만으로 바꿔 쓸 수 있다.

이게 왜 대박이냐면, 유한요소법이나 유한차분법은 3차원 부피 전체를 격자로 채워야 하는데, BEM은 그 표면(2차원)만 이산화하면 되기 때문이다. 미지수 차원이 하나 줄어든다. 무한히 뻗은 영역(음향 방사, 전자기 산란)을 다룰 때 특히 강력한데, 자유공간 그린 함수가 무한대에서의 방사 조건을 이미 만족하고 있어서 인공 경계를 세우고 완전정합층을 붙이는 삽질을 안 해도 되기 때문이다.³ 대신 그린 함수를 쓴 대가로 행렬이 조밀(dense)해져서, 큰 문제에서는 고속 다중극법(FMM) 같은 가속 기법이 필요해진다. 세상에 공짜 점심은 없다.

5. 전자기·구조해석에서의 응용[편집]

그린 함수는 분야를 가리지 않고 선형 문제가 있는 곳이면 어디든 나타난다.

전자기해석: 맥스웰 방정식에서 유도되는 헬름홀츠 방정식의 그린 함수는 안테나 해석과 산란 문제의 뼈대다. 방법의 모멘트(Method of Moments, MoM)라 불리는 적분방정식 해법이 통째로 그린 함수 위에 서 있으며, 유한요소 전자기와 함께 고주파 해석의 양대 축을 이룬다.
구조해석: 탄성체의 점 하중에 대한 응답을 주는 켈빈 해(Kelvin solution)가 탄성역학의 그린 함수다. 이걸 이용한 경계요소법으로 균열 선단의 응력집중 같은 특이점 문제를 정밀하게 다룬다.
양자·통계물리: 다체 문제와 슈뢰딩거 방정식 기반 계산에서 그린 함수는 입자의 전파(propagator)를 나타내는 중심 도구로, 파인만 다이어그램과 섭동 이론의 언어 자체가 되어 버렸다.

정리하면 그린 함수는 “선형 시스템의 유전자”다. 점 하나에 대한 반응만 해독해 두면, 나머지 모든 입력에 대한 반응이 적분 한 번으로 복원된다. 조지 그린이 방앗간에서 굴린 이 아이디어가 200년 뒤 안테나부터 균열까지 시뮬레이션하고 있다.

6. 관련 문서[편집]

조지 그린(George Green, 1793~1841)은 40세에야 케임브리지에 학부생으로 입학한 늦깎이였다. 그전까지는 아버지의 풍차 방앗간을 운영하며 독학으로 수학을 팠다. 정규 교육 1년 남짓으로 근대 수리물리의 언어를 만든 셈이니, 학벌 타령은 여기서 접자. ↩
그래서 물리 전공자는 그린 함수를 처음 배울 때 “이거 우리가 이미 다 하던 거잖아?”라는 기묘한 기시감을 느낀다. 맞다. 정전기학의 절반은 라플라시안 그린 함수 응용이었고, 우리는 그걸 모르고 쓰고 있었을 뿐이다. ↩
유한요소·유한차분법으로 무한 영역을 다루려면 도메인을 억지로 잘라 인공 경계를 세우고 반사를 막는 흡수 경계나 완전정합층(PML)을 발라야 한다. 그린 함수 기반 BEM은 이 성가신 과정을 애초에 건너뛴다. 대신 조밀 행렬이라는 다른 청구서가 날아온다. ↩