세포망 파손 모델의 눈덩이 크기와 대기 시간 분포

초록

본 논문은 대나무 조직의 건조 과정과 미세구조를 모사한 Cell Network Model을 이용해 파손 현상을 연구한다. 파손 임계값을 균등분포와 다양한 Weibull 분포로 설정했을 때, 눈덩이(avalanches) 크기 분포가 지수 -3.0의 파워‑러프를 보이며, 파손 사이의 습도 감소(대기 시간)는 모두 지수분포를 따른다. 남은 결절 비율은 시간에 따라 지수적으로 감소하고, 부분 파손과 누적 손상을 도입하면 눈덩이 크기 히스토그램에서 두 개의 파워‑러프가 교차하는 전이 현상이 나타난다. 이러한 결과는 Cell Network Model이 Random Fuse Model과 동일한 보편성 클래스를 공유한다는 가설을 뒷받침한다.

상세 분석

Cell Network Model(CNM)은 대나무 Guadua angustifolia의 속 조직(parenchymatous tissue)에서 수분이 증발하면서 발생하는 미세 균열을 네트워크 형태로 구현한 모델이다. 각 결절(junction)은 일정한 전기 저항을 갖고, 파손 임계값은 재료의 강도 변동성을 반영해 확률적으로 부여된다. 본 연구에서는 기존에 사용된 균등분포 외에도 Weibull 분포(형상 매개변수 k와 규모 매개변수 λ)를 다양한 조합으로 적용하였다. 시뮬레이션은 습도 H를 단계적으로 감소시키며, H가 임계값 이하가 되면 해당 결절이 파손하고, 인접 결절에 전류가 재분배되면서 연쇄 파손, 즉 avalanche가 발생한다. 각 avalanche는 파손된 결절 수 S로 정의되며, 두 연속 avalanche 사이의 습도 감소 ΔH를 waiting time에 대응시켰다.

시뮬레이션 결과는 다음과 같은 핵심 특징을 보인다. 첫째, S에 대한 확률분포 P(S)는 P(S)∝S^−3 형태의 파워‑러프를 나타냈으며, 이는 Weibull 분포의 형태 매개변수 k가 0.5에서 5까지 변해도 변동이 거의 없었다. 즉, 임계값 분포의 꼬리 두께가 달라져도 지수 −3이라는 스케일링 지수는 강인하게 유지된다. 이는 CNM이 임계 현상에 대한 보편적 거동을 나타낸다는 증거이다.

둘째, waiting time ΔH의 분포는 모든 경우에 지수분포 Q(ΔH)=τ^−1 exp(−ΔH/τ) 형태를 보였다. τ는 평균 습도 감소량으로, Weibull 매개변수에 따라 약간의 변동은 있었지만 전반적으로 일정한 값으로 수렴했다. 이는 파손 이벤트가 메모리 없는 포아송 과정에 가까운 동역학을 갖는다는 의미이며, Random Fuse Model(RFM)에서 보고된 결과와 일치한다.

셋째, 남은 결절 비율 R(t) (시간 t는 총 습도 감소량에 비례)는 R(t)=R0 exp(−γt) 형태의 지수 감쇠를 나타냈다. γ는 파손 속도 상수이며, Weibull 분포의 스케일 λ이 클수록 γ가 작아져 파손이 느려지는 경향을 보였다. 이는 파손 임계값이 전반적으로 높을수록 네트워크가 더 오래 유지된다는 물리적 직관과 부합한다.

넷째, 부분 파손(partial breakings)과 누적 손상(cumulative damage) 메커니즘을 도입하면 눈덩이 크기 히스토그램에 두 개의 파워‑러프 구간이 나타난다. 작은 S 구간에서는 기존 −3 지수와 동일하지만, 일정 크기 S_c를 초과하면 지수가 −2.5 정도로 완만해지는 전이 현상이 관찰된다. 이는 손상이 누적될수록 큰 규모의 연쇄 파손이 더 쉽게 발생함을 시사한다. 이러한 전이는 RFM에서 비선형 저항 모델을 적용했을 때 보고된 전이와 유사하며, CNM이 손상 메커니즘에 따라 보편성 클래스가 변할 수 있음을 암시한다.

전체적으로 본 연구는 임계값 분포의 다양성, 대기 시간의 통계적 특성, 그리고 손상 메커니즘이 CNM의 거시적 스케일링 법칙에 미치는 영향을 체계적으로 검증하였다. 결과는 CNM이 Random Fuse Model과 동일한 보편성 클래스를 공유한다는 가설을 강력히 뒷받침한다.