[논문 리뷰] A nonvariational finite element method for fully nonlinear elliptic problems

Omar Lakkis, Tristan Pryer|arXiv (Cornell University)|2011. 03. 15.

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

한 줄 요약

이 논문은 Lakkis & Pryer (2011)에서 제안한 이산화 기법을 활용하여 완전히 비선형 타원형 PDE를 해결하기 위한 비변분 유한요소법을 제시한다. 이 기법은 선형 PDE의 강한 형태를 직접 해석할 수 있도록 해주며, 자연스럽게 유한요소 헤시안을 부산물로 생성한다. 두 가지 다른 접근 방식이 두 번째 차수의 완전히 비선형 문제에 대해 개발되었으며, 수치 기준을 통해 몽주-암페르 및 푸치 방정식에 대한 수렴성이 확인되었다.

ABSTRACT

We present a continuous finite element method for some examples of fully nonlinear elliptic equation. A key tool is the discretisation proposed in Lakkis & Pryer (2011, SISC) allowing us to work directly on the strong form of a linear PDE. An added benefit to making use of this discretisation method is that a recovered (finite element) Hessian is a biproduct of the solution process. We build on the linear basis and ultimately construct two different methodologies for the solution of second order fully nonlinear PDEs. Benchmark numerical results illustrate the convergence properties of the scheme for some test problems including the Monge-Ampere equation and Pucci's equation.

연구 동기 및 목표

• 완전히 비선형 타원형 PDE를 위한 변분 공식을 필요로 하지 않는 유한요소법을 개발하기 위해.
• 약한 공식을 필요로 하지 않고도 직접 PDE를 해석할 수 있도록 강한 형태의 이산화 기법을 활용하기 위해.
• 해결 과정에서 유도되는 유한요소 헤시안을 자연스럽게 부산물로 활용하기 위해.
• 두 번째 차수의 완전히 비선형 PDE를 해결하기 위한 두 가지 서로 다른 방법론을 구축하기 위해.
• 기본 문제에 대한 제안된 방법의 수렴성과 안정성을 수치적으로 검증하기 위해.

제안 방법

• 방법은 Lakkis & Pryer (2011)에서 소개한 선형 PDE의 강한 형태 기반 연속 유한요소 이산화를 사용한다.
• 해결을 근사하고, 요소별 복원을 통해 헤시안 행렬을 재구성하기 위해 노드 기반 유한요소 기저를 사용한다.
• 해결 과정에서 변분 공식을 필요로 하지 않고, PDE의 강한 형태를 직접 다룸으로써 이를 회피한다.
• 완전히 비선형 PDE를 다루기 위해, 복원된 헤시안을 활용한 두 가지 서로 다른 해결 전략이 개발되었다.
• 이 방법은 몽주-암페르 및 푸치 방정식을 포함한 시험 문제에 적용되었다.
• 구조적 메esh에서의 수치 기준을 통해 수렴성을 평가하였다.

실험 결과

연구 질문

• RQ1 비변분 유한요소법이 완전히 비선형 타원형 PDE에 효과적으로 적용될 수 있는가?
• RQ2 선형 PDE의 강한 형태를 유한요소 프레임워크 내에서 어떻게 이산화할 수 있으며, 이를 통해 비선형 문제를 직접 해석할 수 있는가?
• RQ3 복원된 헤시안이 완전히 비선형 PDE의 해결 과정에서 얼마나 기여하는가?
• RQ4 몽주-암페르 및 푸치 방정식과 같은 기본 비선형 PDE에 대해 제안된 방법의 수렴 성질은 어떠한가?
• RQ5 두 가지 제안된 방법론은 정확도와 안정성 측면에서 어떻게 비교되는가?

주요 결과

• 이 방법은 변분 공식을 필요로 하지 않고도 완전히 비선형 타원형 PDE를 성공적으로 해결한다.
• 해결 과정에서 자연스럽게 유한요소 헤시안이 부산물로 생성되며, 이는 비선형 솔버에서 직접 활용 가능하다.
• 수치 결과는 몽주-암페르 및 푸치 방정식에 대해 제안된 방법에 의한 수렴성을 입증한다.
• 이산화 접근 방식은 표준 유한요소 메쉬에서 두 번째 차수의 완전히 비선형 PDE를 안정적이고 정확하게 해석할 수 있도록 한다.
• 제안된 두 가지 방법론은 기준 시험 문제 전반에서 강력한 수렴 행동을 보인다.
• 이 방법은 유한요소 맥락에서 완전히 비선형 문제에 대해 직접 강한 형태 이산화의 실현 가능성을 보여준다.