[논문 리뷰] Periodic homogenization using the Lippmann--Schwinger formalism

Sébastien Brisard, Frédéric Legoll|arXiv (Cornell University)|2014. 11. 02.

Advanced Mathematical Modeling in Engineering참고 문헌 34인용 수 5

한 줄 요약

이 논문은 주기적 균질화를 위한 Lippmann–Schwinger 형식에 대한 수학적 분석을 제시하며, 이를 적절히 일致하는 갈레르킨 이산화된 적분방정식으로 재구성한다. 약간의 가정 하에 수렴성과 사전 오차 추정을 증명하며, 3차원 선형 탄성 응용 사례에서 메쉬 크기 h에 비례하는 오차를 보이며, 효율적인 병렬 FFT 기반 솔버를 사용한다.

ABSTRACT

When homogenizing elliptic partial differential equations, the so-called corrector problem is pivotal to compute the macroscale effective coefficients from the microscale information. To solve this corrector problem in the periodic setting, Moulinec and Suquet introduced in the mid-nineties a numerical strategy based on the reformulation of that problem as an integral equation (known as the Lippmann--Schwinger equation), which is then suitably discretized. This results in an iterative, matrix-free method, which is of particular interest for complex microstructures. Since the seminal work of Moulinec and Suquet, several variants of their scheme have been proposed. The aim of this contribution is twofold. First, we provide an overview of these methods, recast in the language of the applied mathematics community. These methods are presented as asymptotically consistent Galerkin discretizations of the Lippmann--Schwinger equation. The bilinear form arising in the weak form of this integral equation is indeed the sum of a local and a non-local term. We show that most of the variants proposed in the literature correspond to alternative approximations of this non-local term. Second, we propose a mathematical analysis of the discretized problem. In particular, we prove under mild hypotheses the convergence of these numerical schemes with respect to the grid-size. We also provide a priori error estimates on the solution. The article closes on a three-dimensional numerical application within the framework of linear elasticity.

연구 동기 및 목표

• 주기적 균질화에서 Lippmann–Schwinger 방법에 대한 통합된 수학적 프레임워크를 제공하여 기존의 수치적 방법들을 응용수학의 언어로 재구성한다.
• Moulinec–Suquet 방법의 다양한 변종의 수렴성을 분석하기 위해 Lippmann–Schwinger 적분방정식의 갈레르킨 이산화로 해석한다.
• 계수 필드에 대한 미약한 정규성 가정 하에 이산 해에 대한 사전 오차 추정을 수립한다.
• FFT 기반 반복적 솔버를 사용한 3차원 수치적 응용을 통해 이론적 결과를 검증한다.
• 병렬 구현과 성능 분석을 통해 방법의 효율성과 확장성을 입증한다.

제안 방법

• 기준 매질의 그린 함수를 사용하여 주기적 보정 문제를 Lippmann–Schwinger 적분방정식으로 재구성한다.
• 조각상수 유한요소 기저 함수를 사용한 갈레르킨 방법으로 적분방정식을 이산화하여, 매트릭스-프리 구조를 갖는 선형계로 이끌어낸다.
• 약한 형태에서의 비국소 항을 다양한 적분 구조로 근사함으로써, 다양한 방법의 변종에 대응한다.
• 스펙트럼 영역에서 발생하는 컘벌루션 연산을 효율적으로 계산하기 위해 빠른 푸리에 변환(FFT)을 사용한다.
• FFT를 통해 행렬-벡터 곱을 계산하는 켈레르 기반의 공액 기울기 방법으로 최종 선형계를 푸름으로써, 매트릭스-프리 연산을 가능하게 한다.
• 영역 분할과 MPI를 사용한 병렬화를 통해 여러 프로세서에 걸쳐 계산을 확장하며, 통신 최적화된 데이터 분포를 구현한다.

실험 결과

연구 질문

• RQ1 기존의 주기적 균질화를 위한 Lippmann–Schwinger 방법들은 어떻게 적분방정식의 갈레르킨 이산화로 체계적으로 해석될 수 있는가?
• RQ2 계수 필드에 대한 최소한의 가정 하에 이러한 방법들의 수렴 성질은 어떠한가?
• RQ3 메쉬 크기 h에 대한 이산 해에 대한 사전 오차 추정을 유도할 수 있는가?
• RQ4 수치 실험에서 관측된 수렴 속도가 이론적 예측과 일치하는가, 특히 계수 필드가 조각상수일 경우에?
• RQ5 이 방법은 얼마나 효율적으로 병렬화될 수 있으며, 대규모 프로세서 수에 확장할 때의 성능 저하 요인은 무엇인가?

주요 결과

• 제안된 방법은 미약한 가정 하에 메쉬 크기 h에 대해 수렴하며, 갈레르킨 공식화에 대해 엄밀히 수렴성을 증명하였다.
• 사전 오차 추정이 수립되었으며, 오차가 Ch^δ (δ < 1) 에 의해 유계임을 보였고, 수치 결과는 실질적으로 h에 가까운 수렴 속도를 보였다.
• 조각상수 물성의 3차원 선형 탄성 문제에서, 효과적 전단 모odulus A^⋆_xyxy는 메쉬 정밀도 향상에 따라 증가하며, 512³와 1024³ 메쉬 간에 약 2%의 차이를 보였다.
• 효과 계수의 수치 오차는 h에 비례함을 확인하여, 이론적 경계보다 더 빠른 수렴 속도를 보였으며, 이는 계수 필드의 더 높은 정규성 때문일 가능성이 있다.
• 한 번의 공액 기울기 반복의 계산 비용은 N³ log N 비례하며, 총 시간은 N³ log N × 반복 수에 비례한다. 큰 N에 대해 T / (N³ log N × N_iter) 비율이 일정한 값을 향해 수렴함을 관찰하였다.
• 4개의 프로세스까지는 효율성이 약 0.9로 높았으나, 8개 프로세스에서는 효율성이 약 0.6으로 떨어져, 통신 오버헤드가 중간 수준의 프로세서 수를 초과한 확장성 제한 요인임을 시사하였다.