[논문 리뷰] $C^0$ Discontinuous Galerkin Methods for a Kirchhoff Plate Contact Problem

Fei Wang, Tianyi Zhang외 1명|arXiv (Cornell University)|2014. 10. 25.

Advanced Numerical Methods in Computational Mathematics참고 문헌 17인용 수 6

한 줄 요약

이 논문은 마찰이 있는 키르히호프 플레이트 접촉 문제를 모델링하는 4계 미분형 타원형 변분부등식에 대해 $C^0$ 비연속 갈레르킨(DG) 방법을 개발하고 분석한다. 저자들은 기존의 다섯 가지 DG 스킴—C0 IP, NIPG, Wells-Dung, Baker-DG, LCDG—를 확장하여 엄밀한 일致성 및 안정성 분석을 수행하고, 에너지 노름과 $H^1$ 노름에서 수렴률이 약 1에 가까운 최적 순서 오차 추정치를 유도하였다. 이는 수치적으로 확인되었다.

ABSTRACT

Discontinuous Galerkin (DG) methods are considered for solving a plate contact problem, which is a 4th-order elliptic variational inequality of second kind. Numerous $C^0$ DG schemes for the Kirchhoff plate bending problem are extended to the variational inequality. Properties of the DG methods, such as consistency and stability, are studied, and optimal order error estimates are derived. A numerical example is presented to show the performance of the DG methods; the numerical convergence orders confirm the theoretical prediction.

연구 동기 및 목표

• 4계 미분형 타원형 변분부등식이 발생하는 키르히호프 플레이트 마찰 접촉 문제를 해결하기 위해 기존의 $C^0$ 비연속 갈레르킨 방법을 확장한다.
• 이러한 유형의 변분부등식에 대해 제안된 DG 스킴의 일치성 및 안정성 성질을 확립한다.
• 약한 정규성 조건 하에서 $C^0$ DG 방법에 대한 최적 순서 오차 추정치를 도출한다.
• 메쉬 크기 $h$를 줄이는 방식의 수렴 분석을 통해 이론적 수렴률이 수치적으로 검증된다.

제안 방법

• 수치적 경계값과 페널티 항을 포함하여 4계 미분형 변분부등식에 대해 기존의 다섯 가지 $C^0$ DG 형식—C0 IP, NIPG, Wells-Dung, Baker-DG, LCDG—를 적응시킨다.
• 이산 문제를 다음과 같은 볼록 최소화 문제로 공식화한다: $\min_{u^h \in V^h} \frac{1}{2}a(u^h,u^h) + j_h(u^h) - (f,u^h)$, 여기서 $j_h$ 는 접촉 잠재에너지이다.
• 이산 최적화 문제를 해결하기 위해 원-이중 고정점 알고리즘을 사용하며, 이중 단계에서 $\ell_1$-노름 항에 대해 프록시멀 연산자를 적용한다.
• 오차 수렴을 측정하기 위해 에너지 노름 $\interleave \cdot \interleave$ 과 $H^1(\Omega)$ 반노름을 정의한다.
• 안정성 조건이 명확하게 기술된 일관되고 안정적인 $C^0$ DG 형식(LCDG)을 적용하여, 큰 페널티 계수에 의존하지 않도록 한다.
• 가장 미세한 메쉬($h = 1/64$)의 해를 기준으로 더 넓은 메쉬의 오차를 계산하여 수렴 순서 평가를 수행한다.

실험 결과

연구 질문

• RQ1 기존의 4계 미분형 방정식에 대한 $C^0$ DG 방법은 플레이트 접촉 문제에서 발생하는 4계 미분형 변분부등식을 다룰 수 있는가?
• RQ2 확장된 $C^0$ DG 스킴의 이론적 일치성 및 안정성 성질은 무엇인가?
• RQ3 제안된 DG 방법은 에너지 노름과 $H^1$ 노름에서 최적 수렴률을 달성하는가?
• RQ4 C0 IP, LCDG, Wells-Dung 등의 다양한 DG 형식 간 정확도 및 수렴 행동은 어떻게 비교되는가?
• RQ5 원-이중 고정점 알고리즘이 $\ell_1$-정규화가 포함된 이산 볼록 최적화 문제를 효율적으로 해결할 수 있는가?

주요 결과

• 제안된 $C^0$ DG 방법은 에너지 노름에서 최적 수렴 순서를 달성하며, 모든 다섯 가지 스킴에서 수치적 수렴률이 약 1에 가까운 것으로 확인되었다.
• 수치 결과는 $H^1(\Omega)$ 반노름에서의 수렴 순서도 이론적 예측과 유사하게 약 1에 가까운 것으로 확인되었다.
• LCDG 방법은 명확한 페널티 계수 조건을 제공하여 안정성을 확보하면서도, 안정성을 확보하기 위해 매우 큰 페널티 값을 사용할 필요가 없었다.
• Wells-Dung 및 Baker-DG 형식은 C0 IP 및 NIPG와 유사한 성능을 보였으며, 다양한 페널티 계수 $\eta$ 에서도 일관된 수렴 행동을 보였다.
• 원-이중 고정점 알고리즘은 이산 문제를 효율적으로 해결하였으며, 프록시멀 연산자가 $\ell_1$-정규화 항에 대해 명시적 갱신을 가능하게 하였다.
• 모든 방법에서 $h$ 를 반으로 줄일 때 에너지 노름의 오차가 약 2배 감소함으로써, 일阶 수렴이 확인되었다.