[논문 리뷰] Numerical analysis of the LDG method for large deformations of prestrained plates

Andrea Bonito, Diane Guignard외 2명|arXiv (Cornell University)|2021. 06. 25.

Advanced Numerical Methods in Computational Mathematics참고 문헌 22인용 수 4

한 줄 요약

이 논문은 비선형이고 비볼록인 거리 제약 조건에 따라 조절되는 예비변형이 가해진 판에서의 큰 변형을 시뮬레이션하기 위해 국소 불연속 갈레르킨(LDG) 방법의 수치 분석을 제시한다. 이는 이산 에너지의 $\gamma$-수렴을 확립하고, 이산 경사 하강법에서 에너지 감소를 증명하며, 거리 제약 위반을 통제함으로써 복잡한 판 형상에 대한 안정적이고 정확한 최소화자 계산을 가능하게 한다.

ABSTRACT

A local discontinuous Galerkin (LDG) method for approximating large deformations of prestrained plates is introduced and tested on several insightful numerical examples in our previous computational work. This paper presents a numerical analysis of this LDG method, focusing on the free boundary case. The problem consists of minimizing a fourth order bending energy subject to a nonlinear and nonconvex metric constraint. The energy is discretized using LDG and a discrete gradient flow is used for computing discrete minimizers. We first show $Γ$-convergence of the discrete energy to the continuous one. Then we prove that the discrete gradient flow decreases the energy at each step and computes discrete minimizers with control of the metric constraint defect. We also present a numerical scheme for initialization of the gradient flow, and discuss the conditional stability of it.

연구 동기 및 목표

• 복잡한 비이sov메트릭 기하형상을 가진 예비변형이 가해진 판에서 큰 변형을 시뮬레이션하기 위한 강력한 수치 방법을 개발하는 것.
• 4차 비틀림 에너지를 비선형이고 비볼록인 거리 제약 조건 하에서 최소화하는 데 도전하는 것.
• 메esh 크기가 0으로 갈수록 이산 에너지가 연속 에너지로 $\gamma$-수렴함을 보장하는 것.
• 이산 경사 하강법이 각 단계에서 에너지를 단조 감소시키며 거리 제약 위반을 통제함으로써 최소화자 계산을 보장하는 것.
• 이중 라플라스 방정식 기반의 사전처리 단계를 이용한 경사 하강법에 대한 안정적인 초기화 방법 제공.

제안 방법

• 재구성된 헤시안과 안정화 항을 포함한 국소 불연속 갈레르킨(LDG) 방법을 사용하여 에너지를 이산화한다.
• 최소화자를 계산하기 위해 이산 경사 하강법을 적용하여 각 반복 단계에서 에너지 감소를 보장한다.
• Nitsche의 방법을 통해 거리 제약을 약하게 강제하며, 안정성을 확보하기 위해 안정화 파rameter를 선택한다.
• 이중 라플라스 방정식 문제를 푸는 사전처리 단계를 통해 초기 변형을 생성하여 예비변형 결함을 최소화한다.
• 자유 경계 조건과 딜레르흐 경계 조건을 모두 처리하며, 경계 조건은 Nitsche의 형식을 통해 약하게 적용한다.
• 시간 스텝 제약 조건 하에서 거리 제약 결함의 오차 통제를 포함한 분석을 수행한다.

실험 결과

연구 질문

• RQ1 메쉬 크기가 0으로 갈수록 이산 LDG 에너지가 연속 에너지로 $\gamma$-수렴하는가?
• RQ2 이산 경사 하강법이 에너지를 단조 감소시키며 거리 제약 위반을 통제하는 방식으로 최소화자를 계산할 수 있는가?
• RQ3 경사 하강법 시작 전에 예비변형 결함을 최소화하기 위해 안정적이고 정확한 초기 변형을 어떻게 구성할 수 있는가?
• RQ4 시간 스텝에서 이산 경사 하강법 스킴의 안정성을 보장하는 조건은 무엇인가?
• RQ5 경계 조건이 이산 스킴의 수렴성과 안정성에 미치는 영향은 어떠한가?

주요 결과

• LDG 방법을 통해 구성된 이산 에너지는 메쉬 크기가 0으로 갈수록 연속 에너지로 $\gamma$-수렴한다.
• 이산 경사 하강법은 각 단계에서 에너지를 엄격히 감소시켜 이산 최소화자로의 수렴을 보장한다.
• 거리 제약 결함은 $D_{h}(\mathbf{y}_{h}^{n}) \leq \varepsilon_{0} + c\tau(E_{h}(\mathbf{y}_{h}^{0}) + \widetilde{c})$로 통제되며, $c$는 도메인과 경계에만 의존한다.
• 이중 라플라스 방정식을 이용한 사전처리 단계는 경계 조건을 약간 근사하고 예비변형 결함을 감소시키는 초기 변형 $\widehat{\mathbf{y}}_{h}$를 생성한다.
• 디릴레흐 경계 조건의 경우, 에너지 및 제약 결함의 경계는 수정된 안정성 상수를 통해 유지되어 강인함을 보장한다.
• 메쉬 크기가 충분히 작을 때 $E^{b}_{h}(\widetilde{\mathbf{y}}_{h}^{n_{h}}) \lesssim 1$ 이며 $D_{h}(\widetilde{\mathbf{y}}_{h}^{n_{h}}) \lesssim \sigma_{h}^{1/2}$ 가 성립하며, $E^{p}_{h}(\widetilde{\mathbf{y}}_{h}^{n_{h}}) \lesssim \sigma_{h}$ 라면 거리 제약 통제가 효과적으로 이루어짐을 확인한다.