[논문 리뷰] A nonlocal operator method for solving partial differential equations

Huilong Ren, Xiaoying Zhuang외 1명|arXiv (Cornell University)|2018. 10. 04.

Numerical methods in engineering인용 수 10

한 줄 요약

이 논문은 편미분방정식(PDE)을 해결하기 위해 전통적인 미분 연산자를 테일러 급수 전개로부터 유도된 적분 기반 비국소 연산자로 대체하는 비국소 연산자 방법을 제안한다. 이 방법은 변분 원리와 가중 잔차 방법과의 일致성을 보장하며, 행렬 곱셈을 통한 잔차와 탄성 스텝 행렬의 효율적 계산을 가능하게 하며, 노드 적분 기법에서의 영에너지 모드를 억제하기 위해 시계계단 에너지 함수를 도입한다.

ABSTRACT

We propose a nonlocal operator method for solving partial differential equations (PDEs). The nonlocal operator is derived from the Taylor series expansion of the unknown field, and can be regarded as the integral form "equivalent" to the differential form in the sense of nonlocal interaction. The variation of a nonlocal operator is similar to the derivative of shape function in meshless and finite element methods, thus circumvents difficulty in the calculation of shape function and its derivatives. {The nonlocal operator method is consistent with the variational principle and the weighted residual method, based on which the residual and the tangent stiffness matrix can be obtained with ease.} The nonlocal operator method is equipped with an hourglass energy functional to satisfy the linear consistency of the field. Higher order nonlocal operators and higher order hourglass energy functional are generalized. The functional based on the nonlocal operator converts the construction of residual and stiffness matrix into a series of matrix multiplications on the nonlocal operators. The nonlocal strong forms of different functionals can be obtained easily via support and dual-support, two basic concepts introduced in the paper. Several numerical examples are presented to validate the method.

연구 동기 및 목표

• 형상 함수 구성과 도함수 평가의 복잡성을 피하기 위해 PDE를 해결하기 위한 메쉬프리 유사 방법을 개발하기 위해.
• 노드 적분에서 발생하는 수치적 불안정성, 예를 들어 시계계단 모드를 해결하기 위해 시계계단 에너지 함수를 도입하기 위해.
• 고차 비국소 연산자와 고차 일치 조건으로의 일반화를 위해.
• 신뢰할 수 있는 잔차와 스텝 행렬 계산을 위해 변분 원리와 가중 잔차 방법과의 일치성을 보장하기 위해.
• 불연속성이나 강한 기울기를 갖는 문제, 예를 들어 균열 또는 손상 역학에서의 PDE에 대해 안정적이고 일치적이며 계산적으로 효율적인 프레임워크를 제공하기 위해.

제안 방법

• 비국소 연산자는 테일러 급수 전개로부터 유도되며, 지지 도메인을 통해 정의되어, 미분 연산자의 적분 기반 근사화를 가능하게 한다.
• 이 방법은 '이중지지(dual-support)' 개념을 사용하여 비균일 이산화를 가능하게 하고, 변분 원리를 통한 비국소 강형 형태의 유도를 촉진한다.
• 잔차와 탄성 스텝 행렬은 비국소 연산자와 물리적 구성 법칙 항목을 포함하는 행렬 곱셈을 통해 계산되며, 형상 함수 도함수의 명시적 평가를 피한다.
• 노드 적분의 안정성을 높이기 위해 영에너지 모드를 억제하는 시계계단 에너지 함수를 도입하여 선형 일치성을 확보한다.
• 정확도와 일치성을 향상시키기 위해 고차 비국소 연산자와 그에 해당하는 고차 시계계단 에너지 함수를 일반화한다.
• 지지 도메인과 이중지지 개념을 사용하여 비국소 강형 형태를 구성함으로써, 약형 형태를 대수적 시스템으로 직접 이산화할 수 있다.

실험 결과

연구 질문

• RQ1 PDE의 미분 연산자를 형상 함수 구성 없이 적분 기반 비국소 연산자로 어떻게 대체할 수 있는가?
• RQ2 비국소 연산자 방법은 변분 원리와 가중 잔차 방법과의 일치성을 유지할 수 있는가?
• RQ3 비국소 프레임워크 내에서 노드 적분에서의 시계계단 모드를 효과적으로 억제할 수 있는가?
• RQ4 이중지지의 역할은 비균일 이산화와 안정된 공식화를 가능하게 하는 데서 무엇인가?
• RQ5 고차 비국소 연산자와 에너지 함수는 방법의 정확도와 일치성을 어떻게 향상시키는가?

주요 결과

• 비국소 연산자 방법은 국소 미분 연산자를 적분 기반 형태로 성공적으로 대체하여, 명시적 형상 함수 도함수 없이도 안정적이고 정확한 PDE 해법을 가능하게 한다.
• 이 방법은 변분 원리와 가중 잔차 방법과의 일치성을 확보하여, 행렬 곱셈을 통한 잔차와 탄성 스텝 행렬의 직접 계산이 가능하다.
• 시계계단 에너지 함수는 노드 적분에서의 영에너지 모드를 효과적으로 억제하여, 명시적 및 암시적 분석 모두에서 수치적 안정성을 확보한다.
• 비국소 연산자는 무한소 지지의 극한에서 국소 미분 연산자로 수렴하여 고전적 PDE와의 일致성을 보장한다.
• 단순비, 포아송 방정식, 고유값 문제를 포함한 수치 예제들은 불연속성이 있는 문제에서의 정확도와 강건성을 입증한다.
• 이중지지의 사용은 민첩하고 비균일 이산화를 가능하게 하여, 복잡한 기하학 및 비정규 메쉬에 대한 방법의 적용 가능성을 높인다.