유전파괴 모형

계산물리 전자기해석 시뮬레이션 마지막 수정: 2026-08-15 04:51:07

1. 개요[편집]

유전파괴 모형 Dielectric Breakdown Model
약칭	DBM
제안	1984년, Niemeyer · Pietronero · Wiesmann
지배 방정식	∇²φ = 0 (매 성장 단계마다 재계산)
성장 규칙	p ∝ 국소장 크기의 η 제곱
η = 1	DLA와 같은 보편성 (2D에서 D ≈ 1.71)
대응 현상	리히텐베르크 그림, 전기 트리, 수지상 방전

유전파괴 모형(dielectric breakdown model, DBM)은 절연체가 전기적으로 파괴될 때 생기는 나뭇가지 모양 방전 경로를, 매 단계마다 라플라스 방정식을 실제로 풀고 그 해에서 나온 국소 전기장의 $\eta$ 제곱에 비례하는 확률로 클러스터를 한 칸씩 키우는 확률적 프랙탈 성장 모형이다. 1984년 니마이어·피에트로네로·비스만이 제안했다.¹

확산 제한 응집(DLA)과 형제 사이지만 발상의 출발점이 반대다. DLA는 입자를 무작위로 걷게 해서 확률장을 간접적으로 만들고, DBM은 장 방정식을 직접 풀어 확률을 만든다. 그래서 이 문서는 프랙탈 클러스터가 왜 생기는지, 차폐가 무엇인지 같은 이야기를 반복하지 않는다 — 그건 DLA 문서에 있다. 여기서는 “장을 풀어 확률을 만든다”가 무엇을 가능하게 하는가에 집중한다. 요약하면 세 가지다. 지수 $\eta$ 라는 손잡이가 생기고, 임의의 전극 형상과 경계조건을 그대로 넣을 수 있고, 실제 방전의 물리(문턱 전계, 채널 내부 전위 강하)를 모형에 얹을 수 있다.

2. 알고리즘[편집]

격자 위에서 돌린다고 하면 절차는 다음과 같다.

씨앗(음극/침 전극)을 클러스터로 두고 그 위의 전위를 $\phi = 0$ , 바깥 경계(양극)를 $\phi = 1$ 로 고정한다.
클러스터 밖 영역에서 $\nabla^2\phi = 0$ 을 푼다.
클러스터에 인접한 아직 안 자란 격자점(성장 후보) $j$ 를 모으고, 각 후보에 확률을 부여한다.

p_j = \frac{\lvert \phi_j \rvert^{\eta}}{\displaystyle\sum_{k \in \text{후보}} \lvert \phi_k \rvert^{\eta}}

이 확률로 후보 하나를 뽑아 클러스터에 편입시키고, 그 점의 전위를 0으로 고정한다.
2번으로 돌아간다.

클러스터 위에서 $\phi = 0$ 이고 격자 간격이 $a$ 이므로, 후보점의 전위는 $\phi_j \approx a\,\lvert\nabla\phi\rvert_j$ 다. 즉 위 식은 성장확률이 국소 전기장 크기의 $\eta$ 제곱에 비례한다는 규칙과 같다. 이게 모형의 전부다.

p \;\propto\; \lvert \mathbf{E} \rvert^{\eta} = \lvert \nabla\phi \rvert^{\eta}

무작위성이 들어가는 곳이 딱 한 군데뿐이라는 점을 눈여겨보자. 장 계산은 완전히 결정론적이고, 주사위는 4번에서 단 한 번 굴린다. DLA에서 입자 하나가 클러스터에 도달하기까지 수만 번 난수를 뽑는 것과 대조적이다. 그래서 DBM은 **“확률적 규칙을 가진 이동경계 문제”**라는 성격이 훨씬 노골적으로 드러난다.

$\eta$ 의 물리적 의미는 국소 파괴 기준의 비선형성이다. $\eta = 1$ 은 “성장 속도 = 전류 밀도”에 해당하고, $\eta > 1$ 은 장이 조금만 세도 훨씬 빨리 뚫린다는 뜻이라 실제 기체 방전의 강한 문턱 특성을 흉내 낸다. 원논문이 $\eta$ 를 도입한 것도 실험에서 관측된 프랙탈 차원을 맞추기 위한 현상론적 손잡이로서였다.

3. η 가 만드는 형태의 스펙트럼[편집]

$\eta$ 를 돌리면 같은 알고리즘이 전혀 다른 모양을 낸다.

$\eta$	성장 규칙의 성격	결과	2D 프랙탈 차원
0	후보 전부 동일 확률	에덴 모형형 조밀 덩어리	2
0.5	약한 팁 선호	표면이 거친 조밀 클러스터	2 에 가까움
1	장에 정비례	DLA와 동일	$\approx 1.71$
2~3	강한 팁 선호	가지가 성글어지고 굵은 줄기 소수	1.7 미만, 단조 감소
$\gtrsim 4$	사실상 최대 장 지점만 성장	거의 갈라지지 않는 준선형 침	$\to 1$

핵심은 두 극한 사이의 연속적 다리라는 점이다. $\eta \to 0$ 이면 전기장 정보를 완전히 무시하므로 표면 어디나 균등하게 자라고, 이는 곧 조밀한 에덴 모형이다. $\eta = 1$ 은 성장확률이 조화측도와 정확히 같아져 DLA와 같은 보편성 부류에 들어간다 — 두 모형이 서로 다른 방식으로 정의됐는데 같은 차원을 준다는 사실은 라플라스 성장이 하나의 통일된 대상임을 보여 주는 대표적 증거다.² $\eta$ 를 계속 키우면 차폐가 극단으로 치달아 가장 앞선 팁 하나가 모든 성장을 독식하고, 2차원에서는 대략 $\eta \approx 4$ 부근에서 차원이 1에 도달해 클러스터가 거의 직선 침이 된다. 그 위로는 형태가 더 변하지 않는다.

$D(\eta)$ 곡선을 수치적으로 재는 일 자체는 고약하다. $\eta$ 가 커질수록 확률 분포가 소수의 후보에 몰려 유효 표본이 급감하고, 유한 크기 보정이 커진다. 그래서 문헌마다 $\eta_c$ 값이 조금씩 다르고, 격자 없는 등각사상 기반 계산이 이 논쟁을 정리하는 데 큰 역할을 했다.

4. 실제 방전과 얼마나 닮았나[편집]

원논문의 설득력은 실험과의 직접 비교에서 나왔다. 유전체 판 표면에 침 전극으로 고전압을 걸어 만든 표면 방전 무늬 — 흔히 **리히텐베르크 그림**이라 불리는 그 패턴 — 의 프랙탈 차원을 재면 1.7 언저리가 나오고, 이는 $\eta = 1$ 인 DBM 클러스터와 잘 맞는다. 상자 세기로 잰 숫자 하나가 맞는 것을 두고 물리를 다 이해했다고 할 수는 없지만, 순수하게 기하학적인 규칙 두 줄이 실제 절연 파괴 무늬의 통계를 재현한다는 사실은 충분히 인상적이었다.

그 뒤로 DBM 계열이 실제 공학에 쓰이는 곳은 대체로 다음과 같다.

절연 파괴 경로 예측. 고전압 케이블·부싱·에폭시 몰드에서 시간이 지나며 자라는 전기 트리(electrical treeing)의 형태와 수명 예측. 실제 절연물의 열화는 수년 걸리는데, 그 결과물의 기하는 DBM으로 몇 초 만에 나온다.
뇌방전·스트리머 모형. 리더 채널의 가지치기 통계를 재현하는 데 쓴다. 다만 이쪽은 순정 DBM으로는 부족하고 보정이 필요하다(아래).
전착·수지상 성장. 리튬 전지의 덴드라이트나 전기도금 결함도 전해질 안의 전위장이 지배하므로 같은 틀에 들어간다.

보정 없이는 안 되는 지점이 명확하다는 것도 짚어야 한다.

문턱 전계. 순정 DBM은 아무리 약한 장에서도 확률만 작을 뿐 성장한다. 실제 기체·고체 절연물에는 임계 전계 $E_c$ 가 있어 그 아래에서는 아예 안 뚫린다. 그래서 $\lvert\nabla\phi\rvert < E_c$ 인 후보를 후보 목록에서 제외하는 변형이 표준으로 쓰이고, 이걸 넣으면 무작정 퍼지지 않고 뚜렷한 줄기와 가지가 생겨 실제 스트리머에 훨씬 가까워진다.
채널 내부 전위 강하. 클러스터 전체를 $\phi = 0$ 인 완전 도체로 두는 것은 방전 채널의 저항을 0으로 보는 것이다. 실제 채널은 유한 저항이라 뿌리에서 팁까지 전위가 떨어지고, 그 때문에 긴 가지의 팁이 상대적으로 불리해져 가지치기 패턴이 달라진다. 채널을 따라 $\phi$ 를 선형으로 떨어뜨리는 것만으로도 형태가 크게 개선된다.
시간이 없다. DBM은 한 단계에 격자 한 칸이라는 순서만 줄 뿐 물리적 시간 척도가 없다. 전파 속도나 시간 의존 현상을 보려면 유체형 스트리머 모형(전자 밀도 이송 + 포아송 방정식)으로 넘어가야 한다.

5. 구현에서 실제로 어려운 것[편집]

DBM의 계산 비용은 전적으로 매 스텝 라플라스 방정식을 다시 푸는 데 있다. 클러스터에 $N$ 개의 점을 붙이려면 원칙적으로 경계값 문제를 $N$ 번 풀어야 하고, 격자점이 $M$ 개면 순진하게 가우스-자이델이나 SOR로 매번 수렴시킬 때 총 비용이 감당이 안 된다. 실무의 요령은 다음과 같다.

직전 해를 초기 추정치로 쓴다. 격자점 하나가 추가됐을 뿐이므로 해는 국소적으로만 바뀐다. 반복법을 0에서 시작하지 않고 이전 해에서 이어 돌리면 몇 회 만에 다시 수렴한다. 이 한 가지만으로 비용이 자릿수 단위로 줄어든다.
**다중격자법**으로 각 스텝의 라플라스 풀이를 $O(M)$ 으로 만든다.
그린 함수 갱신. 새 점이 하나 추가될 때 전위장의 변화를 그린 함수로 표현해 갱신하면 반복법 자체를 피할 수 있다. 경계요소적 관점과 통한다.
등각사상 기반 오프-래티스 계산. 2차원에서는 성장 단계 하나를 복소평면의 등각사상 한 번으로 표현할 수 있다(헤이스팅스-레비토프 방식). 격자를 아예 쓰지 않으므로 격자 이방성 문제가 사라지고, $\eta$ 를 연속적으로 바꿔 가며 $D(\eta)$ 를 재는 데 특히 유리하다.

격자 위에서 굴릴 때의 함정은 DLA와 동일하다. 정사각 격자는 축 방향을 미세하게 편애하고, 성장의 양의 되먹임이 그 편향을 증폭해서 큰 클러스터가 십자 모양으로 기운다. 차원을 소수 둘째 자리까지 주장하려면 오프-래티스로 가거나, 최소한 씨앗을 바꿔 수십 회 이상 평균 내고 초기·말기 구간을 회귀에서 빼야 한다.

6. 라플라스 성장 가족 안에서의 위치[편집]

DBM의 진짜 값어치는 특정 현상 하나를 잘 맞히는 데 있지 않고, 하나의 손잡이로 여러 성장 모형을 잇는다는 데 있다. 경계에서 $\phi$ 를 0으로 고정하고 법선 속도를 $\lvert\nabla\phi\rvert^\eta$ 로 주는 틀 안에서 확산 제한 응집, 에덴 모형, 새프먼-테일러 불안정의 점성 손가락, 전착, 절연 파괴가 전부 특수한 경우로 들어온다(대응 관계 표는 DLA 문서에 있다).

여기서 자연스럽게 나오는 질문 하나. 표면장력이 없는 순수 라플라스 성장은 유한 시간에 경계에 첨점이 생겨 해가 깨지는 것으로 알려져 있는데, DBM은 왜 안 깨지나? 답은 DLA와 같다 — 격자 간격이라는 최소 길이와 성장 규칙의 무작위성이 정칙화 역할을 한다. 잡음이 특이점을 미리 잘라 주는 셈인데, 뒤집어 말하면 DBM 클러스터의 모양에는 실제 물리(표면장력, 결정 이방성)가 아니라 격자와 난수가 정한 몫이 섞여 있다는 뜻이기도 하다. 그래서 실제 재료의 수지상 성장을 정량적으로 재현하려는 목적이라면 위상장 모형처럼 계면 에너지를 명시적으로 담는 틀로 넘어가는 것이 정공법이다.³

7. 관련 문서[편집]

Niemeyer, L., Pietronero, L. & Wiesmann, H. J. (1984). “Fractal Dimension of Dielectric Breakdown”. Phys. Rev. Lett. 52, 1033. 저자 중 둘이 스위스 ABB 계열 고전압 연구소 소속이었다는 점이 이 모형의 성격을 잘 말해 준다. 통계물리 논문처럼 생겼지만 출발점은 “케이블이 언제 어디로 뚫리는가”라는 지극히 산업적인 질문이었다. ↩
더 강한 결과도 있다. $\eta = 1$ 인 DBM은 이른바 라플라스 무작위걷기와 정확히 동등하다는 것이 보여졌고, 이것이 DBM( $\eta=1$ )과 DLA가 같은 차원을 주는 이유를 규칙 수준에서 설명한다. 물론 “왜 하필 1.71인가”는 여전히 아무도 해석적으로 못 푼 상태다. 두 모형이 같은 미지수로 수렴한다는 것만 아는 셈. ↩
그래서 “번개는 프랙탈이다”라는 대중적 서술은 DBM이 준 통찰의 절반만 담고 있다. 나머지 절반은 어떤 프랙탈인지가 국소 파괴 법칙의 비선형성( $\eta$ )과 채널의 전기적 성질에 달려 있다는 것이고, 실무에서 돈이 걸린 부분은 정확히 이 나머지 절반이다. 차원 숫자 하나로는 케이블이 언제 뚫릴지 알 수 없다. ↩