[논문 리뷰] A numerical methodology for enforcing maximum principles and the non-negative constraint for transient diffusion equations

K. B. Nakshatrala, Harsha Nagarajan외 1명|arXiv (Cornell University)|2012. 06. 04.

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

한 줄 요약

이 논문은 시간에 따라 변화하는 비등방성 확산 방정식에서 최대 원칙과 비음성 제약 조건을 만족시키기 위해 수평선 방법(Rothe 방법)에 기반한 새로운 수치적 방법론을 제안한다. 먼저 시간을 이산화하고, 각 시간 단계에서 비등방성 확산-감쇠 문제를 볼록 2차 프로그래밍을 통해 풀어, 일반적인 비정렬 메쉬에서 시간 스텝 제약 없이 물리적으로 일관되고 음이 아닌 해를 보장한다. 이는 표준 유한요소 형식이 이러한 기본 수학적 성질을 위반하는 것과 대비된다.

ABSTRACT

Transient diffusion equations arise in many branches of engineering and applied sciences (e.g., heat transfer and mass transfer), and are parabolic partial differential equations. It is well-known that, under certain assumptions on the input data, these equations satisfy important mathematical properties like maximum principles and the non-negative constraint, which have implications in mathematical modeling. However, existing numerical formulations for these types of equations do not, in general, satisfy maximum principles and the non-negative constraint. In this paper, we present a methodology for enforcing maximum principles and the non-negative constraint for transient anisotropic diffusion equation. The method of horizontal lines (also known as the Rothe method) is applied in which the time is discretized first. This results in solving steady anisotropic diffusion equation with decay equation at every discrete time level. The proposed methodology for transient anisotropic diffusion equation will satisfy maximum principles and the non-negative constraint on general computational grids, and with no additional restrictions on the time step. We illustrate the performance and accuracy of the proposed formulation using representative numerical examples. We also perform numerical convergence of the proposed methodology. For comparison, we also present the results from the standard single-field semi-discrete formulation and the results from a popular software package, which all will violate maximum principles and the non-negative constraint.

연구 동기 및 목표

• 일시적 확산 방정식에 대한 표준 수치적 방법에서 최대 원칙과 비음성 제약 조건을 만족하지 못하는 문제를 해결하기 위해.
• 일시적 비등방성 확산 문제에 대해 물리적으로 일관되고 음이 아닌 해를 보장하는 방법론을 개발하기 위해.
• 비음성 조건을 유지하기 위해 일반적으로 필요로 하는 엄격한 시간 스텝 제약 조건을 제거하기 위해.
• 일般적인 비정렬 메쉬에서 작동하는 견고하고 메쉬에 의존하지 않는 설정을 제공하기 위해.
• 제안된 방법이 표준 반디스cret화 형식과 상용 소프트웨어에 비해 수학적 및 물리적 제약 조건을 더 잘 유지하는지 보여주기 위해.

제안 방법

• 시간을 먼저 이산화하여 포물형 편미분방정식을 각 시간 단계에서의 정 steady 비등방성 확산-감쇠 문제의 연속으로 변환함으로써 수평선 방법(Rothe 방법)을 적용한다.
• 각 시간 단계에서 유도된 정상 상태 문제는 비음성 제약 조건과 최대 원칙을 강제하는 볼록 2차 프로그래밍 설정을 통해 해결된다.
• 공간 이산화가 예각 삼각형과 같은 메쉬 품질 조건을 요구하지 않는 일반적인 비정렬 계산 메쉬에 적용된다.
• 가중치가 부여된 시간 수준(t_{n+γ})에서 시간 이산화를 수행함으로써, 보간을 통한 정수 시간 단계에서의 해 재구성에 안정적이고 정확한 방법을 확보한다.
• 가중치가 부여된 설정의 비자기성 문제를 해결하기 위해 첫 번째 시간 단계에서는 후진 오일러 스킴(γ=1)을 사용하는 자가 시작 전략을 채택한다.
• 이 방법은 이산화 시스템의 우변이 음이 아니게 유지됨을 보장하여, 해에서 비음성 제약 조건을 유지하는 데 필수적인 조건을 충족시킨다.

실험 결과

연구 질문

• RQ1 일반적인 비정렬 메쉬에서 일시적 비등방성 확산 방정식에 대해 비음성 제약 조건과 최대 원칙을 강제할 수 있는 수치적 방법을 개발할 수 있는가?
• RQ2 제안된 방법이 비음성 조건을 유지하기 위해 일반적으로 요구되는 엄격한 시간 스텝 제약 조건이 필요 없도록 할 수 있는가?
• RQ3 최적화 기반의 제안된 방법은 표준 단일장 반디스크리트화 형식과 상용 소프트웨어에 비해 물리적 제약 조건을 얼마나 잘 유지하는가?
• RQ4 예각 삼각형 요소와 같은 메쉬 품질 조건이 필요 없이도 안정성과 정확성을 유지할 수 있는가?
• RQ5 이 방법은 복잡한 기하 구조와 이질적, 비등방성 확산 계수를 가진 문제에 적용 가능한가?

주요 결과

• 제안된 방법은 표준 유한요소 형식과는 달리 일반적인 비정렬 메쉬에서 일시적 비등방성 확산 방정식에 대해 비음성 제약 조건과 최대 원칙을 성공적으로 강제한다.
• 비음성 조건을 유지하기 위해 시간 스텝 제약 조건이 필요 없어, 물리적 제약 조건을 위반하지 않으면서도 더 큰 시간 스텝을 허용한다.
• 수치적 예제들은 제안된 설정이 비등방성 및 이질적 확산 계수를 가진 다양한 시험 케이스에서 정확성과 안정성을 유지함을 보여준다.
• 제안된 방법은 표준 반디스크리트 유한요소 형식과 인기 있는 상용 소프트웨어 패키지보다 뛰어나며, 모두 비음성 제약 조건과 최대 원칙을 위반함을 보여준다.
• 볼록 2차 프로그래밍 접근법은 우변이 음일 경우에도 이산화 시스템이 항상 잘 정의되고 물리적으로 의미 있는 상태를 유지함을 보장한다.
• 가중치가 부여된 시간 수준에 기반한 보간 기법은 정수 시간 단계에서의 해 재구성에 정확성을 제공하며, 첫 번째 단계에서 후진 오일러 스킴을 함께 사용함으로써 방법이 자가 시작 가능해진다.