안정적인 계면 균열 성장 모델

초록

제한된 기하학적 조건에서 안정적인 계면 균열의 성장 메커니즘을 모사한 모델을 제시한다. 실험에서 관찰되는 두 가지 서로 다른 스케일의 자기친화적 거동을 하나의 시스템 안에서 재현하고, 각각의 스케일에 대응하는 물리적 원인을 분석한다.

상세 분석

본 논문은 계면 균열이 제한된 두께의 재료 내에서 진행될 때 나타나는 복합적인 거동을 이론적으로 설명하기 위해, 비선형 탄성 변형과 랜덤 결함 분포를 결합한 격자 기반 모델을 구축하였다. 모델은 2차원 격자에 각 사이트마다 파괴 강도 임계값을 부여하고, 외부 하중이 서서히 증가함에 따라 응력 집중이 발생하는 방식을 구현한다. 특히, 균열 전면에 존재하는 잔류 응력이 장거리 상호작용을 통해 균열 전진을 억제하거나 촉진하는 메커니즘을 정량화하였다. 이러한 상호작용은 전통적인 라프라시안(∇²) 형태의 탄성 커플링이 아니라, 파워 로우 형태의 비국소 커플링으로 기술되어, 실험에서 보고된 장거리 상관성을 재현한다. 모델 내에서 두 개의 주요 파라미터, 즉 결함 강도 분포의 폭(σ)과 비국소 상호작용 지수(α)가 스케일 전환을 결정한다. 작은 스케일에서는 결함 강도 변동이 주도하여 전형적인 1+1 차원 KPZ(Kardar‑Parisi‑Zhang) 지수와 유사한 거동을 보이며, 큰 스케일에서는 비국소 탄성 상호작용이 지배하여 새로운 자기친화적 지수(ζ≈0.8)를 만든다. 수치 시뮬레이션 결과는 구조 함수와 파워 스펙트럼 분석을 통해 두 개의 뚜렷한 스케일 구간을 확인했으며, 각각의 구간에서 추정된 거칠기 지수는 실험값과 정량적으로 일치한다. 따라서 모델은 두 물리적 메커니즘—결함 주도와 비국소 탄성 주도—을 하나의 프레임워크 안에 통합함으로써, 기존 실험 해석에 필요한 별도 모델링 없이도 복합 스케일링을 설명한다는 점에서 혁신적이다. 또한, 파라미터 공간을 탐색함으로써 전이점이 어떻게 변하는지, 그리고 특정 재료 설계에서 어느 메커니즘이 우세하게 작용할지를 예측할 수 있는 실용적 가이드를 제공한다.