[논문 리뷰] The role of the patch test in 2D atomistic-to-continuum coupling methods

Christoph Ortner|arXiv (Cornell University)|2011. 01. 27.

Microstructure and mechanical properties참고 문헌 37인용 수 8

한 줄 요약

이 논문은 2차원 원자역학-연속체(a/c) 상호작용 방법에서 페치 테스트 일致성—유령력 제거—이 에너지 일치성, 국소적 인터페이스 보정, 보정의 체적 스케일링, 원자역학 영역의 연결성이라는 네 가지 핵심 기술 조건을 만족할 경우, 일阶 정확도를 보장함을 입증한다. 분석은 W^{1,p} 노름에서 새로운 응력 텐서 구성법과 변분 범죄 기법을 통해 일치성과 수렴 속도를 연결함으로써 에너지 기반 a/c 방법에 대한 엄밀한 이론적 기반을 제공한다.

ABSTRACT

For a general class of atomistic-to-continuum coupling methods, coupling multi-body interatomic potentials with a P1-finite element discretisation of Cauchy--Born nonlinear elasticity, this paper adresses the question whether patch test consistency (or, absence of ghost forces) implies a first-order error estimate. In two dimensions it is shown that this is indeed true under the following additional technical assumptions: (i) an energy consistency condition, (ii) locality of the interface correction, (iii) volumetric scaling of the interface correction, and (iv) connectedness of the atomistic region. The extent to which these assumptions are necessary is discussed in detail.

연구 동기 및 목표

• 2차원 a/c 방법에서 페치 테스트 일치성이 일阶 오차 추정치를 유도하는지 여부를 규명하는 것.
• 일치성이 수렴과 연결되지 않았던 에너지 기반 a/c 방법의 이론적 이해에 있어 임계적 간극을 메우는 것.
• 유령력이 없는 방법이 최적의 수렴 속도를 달성할 수 있는 충분한 조건을 설정하는 것.
• W^{1,p} 노름과 변분 범죄 기법을 사용하여 a/c 방법의 오차 추정에 대한 일반적인 분석 프레임워크를 개발하는 것.

제안 방법

• 다체 원자간 퍼텐셜과 Cauchy–Born 탄성에 대한 P1-유한요소를 결합하는 일반적인 a/c 방법의 클래스를 도입한다.
• 2차원 및 3차원 a/c 결합을 위해 특화된 새로운 기호 체계와 분석 프레임워크를 개발하며, 국소적 스무스함을 측정하기 위해 이산 진동 연산자를 포함한다.
• Shapeev의 결합 밀도 보조정리를 활용하여 결합 적분을 체적 적분과 연결함으로써 응력 텐서 비교를 가능하게 한다.
• 비틀림 응력과 유사한 원자역학 및 a/c 응력 함수를 구성함으로써 응력 일치성에 기반한 오차 분해를 가능하게 한다.
• 수정된 a/c 응력 함수와 보정 함수를 사용하여 인터페이스 보정이 존재하는 상황에서도 모델링 오차를 제어한다.
• 보간 오차 추정과 쌍대 노름 기법을 적용하여 W^{-1,p} 노름에서 모델링 오차를 경계하며, 핵심 안정성 과제를 분리한다.

실험 결과

연구 질문

• RQ1 2차원 a/c 방법에서 페치 테스트 일치성이 에너지 기반 설정에서 일阶 수렴을 유도하는가?
• RQ2 유령력 제거가 최적의 오차 추정치를 보장하기 위해 추가로 필요한 기술적 가정은 무엇인가?
• RQ3 원자역학 및 a/c 모델의 응력 텐서를 어떻게 구성하여 엄밀한 오차 분석을 가능하게 할 수 있는가?
• RQ4 에너지 일치성, 국소적 보정, 체적 스케일링, 원자역학 영역의 연결성이라는 가정들이 일阶 수렴을 위해 얼마나 필수적인가?
• RQ5 메쉬 굵기가 조절되더라도 모델링 오차를 독립적으로 경계할 수 있는가? 어떤 조건에서 가능한가?

주요 결과

• 페치 테스트 일치성만으로는 일阶 수렴이 불충분하며, 네 가지 추가 기술적 가정이 필요하다.
• 메서드가 에너지 일치성, 국소적 인터페이스 보정, 보정의 체적 스케일링, 원자역학 영역의 연결성을 만족할 경우 2차원에서 일阶 일致성이 달성된다.
• 원자역학 및 a/c 응력 함수의 구성은 모델링 오차를 응력 일치성 항과 보정 항으로 정밀하게 분해할 수 있게 한다.
• 수정된 a/c 응력 함수의 사용은 인터페이스 보정이 존재하는 상황에서도 오차를 제어할 수 있으며, 보정 함수는 인터페이스에서의 일치성을 보장한다.
• 분석 결과, 안정성 가정이 만족되면 모델링 오차는 메쉬 굵기와 무관하게 경계될 수 있음을 보여준다.
• 주요 결과인 정리 6은 제시된 가정 하에 a/c 방법이 W^{1,p} 노름에서 일阶 일치성을 달성함을 증명한다.