[논문 리뷰] Numerical Methods for Solving Convection-Diffusion Problems

Alexander Churbanov, Petr N. Vabishchevich|arXiv (Cornell University)|2012. 08. 28.

Advanced Numerical Methods in Computational Mathematics참고 문헌 36인용 수 79

한 줄 요약

이 논문은 비압축성, 비압축성, 그리고 비대칭 형태의 대류-확산 방정식에 대해 비조건적 안정성과 단조성을 보장하는 유한차분 및 유한요소 스킴을 제시한다. 힐버트 공간과 바나흐 공간 모두에서 안정성과 수렴성을 확립하였으며, 새로운 국소 1차원 및 덧셈 평균화 분할 스킴을 통해 시간 간격에 관계없이 단조성과 최적 오차 한계를 확보한다. 이는 확산 및 속도 장에 적절한 조건이 있을 경우에 성립한다.

ABSTRACT

Convection-diffusion equations provide the basis for describing heat and mass transfer phenomena as well as processes of continuum mechanics. To handle flows in porous media, the fundamental issue is to model correctly the convective transport of individual phases. Moreover, for compressible media, the pressure equation itself is just a time-dependent convection-diffusion equation. For different problems, a convection-diffusion equation may be be written in various forms. The most popular formulation of convective transport employs the divergent (conservative) form. In some cases, the nondivergent (characteristic) form seems to be preferable. The so-called skew-symmetric form of convective transport operators that is the half-sum of the operators in the divergent and nondivergent forms is of great interest in some applications. Here we discuss the basic classes of discretization in space: finite difference schemes on rectangular grids, approximations on general polyhedra (the finite volume method), and finite element procedures. The key properties of discrete operators are studied for convective and diffusive transport. We emphasize the problems of constructing approximations for convection and diffusion operators that satisfy the maximum principle at the discrete level --- they are called monotone approximations. Two- and three-level schemes are investigated for transient problems. Unconditionally stable explicit-implicit schemes are developed for convection-diffusion problems. Stability conditions are obtained both in finite-dimensional Hilbert spaces and in Banach spaces depending on the form in which the convection-diffusion equation is written.

연구 동기 및 목표

• 유체역학 및 열/질량 이동 현상에서 발생하는 대류-확산 방정식에 대해 강건하고 비조건적 안정성 수치 해법을 개발한다.
• 대류 및 확산 연산자의 단조적 근사로 이산 최대원칙을 만족시키는 것을 확보한다.
• 최적의 안정성과 수렴 성질을 갖는 2단계 및 3단계 시간 이산화 스킴을 구축한다.
• 다양한 형태의 대류-확산 방정식에 대해 유한차원 힐버트 공간과 바나흐 공간에서의 안정성과 단조성 분석을 수행한다.
• 실제 구현에 적합한 효율적이고 병렬 처리 가능한 분할 스킴—국소 1차원 및 덧셈 평균화 스킴—을 설계한다.

제안 방법

• 직각 격자에서의 유한차분 스킴, 일반 다면체 메esh에서의 유한체적 방법, 그리고 일치하는 유한요소 방법을 사용하여 공간 이산화를 수행한다.
• 대류 운반 연산자의 보존형(발산형), 비보존형(특성형), 비대칭형을 사용하여 다양한 물리적 거동을 모델링한다.
• 파arameter σ를 사용한 가중치 시간 스텝 스킴(예: 크랭크-니콜슨 유사)을 적용하여 2차 정밀도와 비조건적 안정성을 확보한다.
• 진화 연산자의 곱셈 분할을 통한 국소 1차원(LOD) 스킴을 도입하여, 시간 간격 제약 조건 하에서 단조성을 보장한다.
• 덧셈 분할을 사용한 덧셈 평균화 LOD 스킴을 제안하여 구성 요소를 별도로 계산할 수 있도록 하여 병렬 처리 성능을 향상시킨다.
• 행렬 분석을 통해 안정성 및 단조성 조건을 유도하며, 반복 행렬이 M-행렬이 되고 그 노름이 1 이하가 되도록 요구한다.

실험 결과

연구 질문

• RQ1 보존형 및 비보존형 모두에서 대류-확산 방정식에 대해 비조건적 안정성과 단조성을 보장하는 수치 해법을 구성할 수 있는가?
• RQ2 대류 운반 연산자의 다양한 형태(발산형, 비발산형, 비대칭형)가 수치 해의 안정성과 정밀도에 어떤 영향을 미치는가?
• RQ3 유한차원 힐버트 공간과 바나흐 공간에서 시간 분할 스킴의 단조성과 안정성에 필요한 충분한 조건은 무엇인가?
• RQ4 국소 1차원 또는 덧셈 평균화 스킴이 임의의 시간 간격에서 단조성과 최적 오차 한계를 유지할 수 있는가?
• RQ5 기존 방법과 비교해 병렬 처리 가능성과 계산 효율성 측면에서 제안된 스킴의 성능은 어떠한가?

주요 결과

• 곱셈 분할을 통한 국소 1차원 스킴은 σ = 1일 경우 임의의 τ > 0에서 L∞ 또는 L1 공간에서 사전 추정치를 만족하며 단조성을 보장한다. σ < 1일 경우 시간 간격 제약 조건 하에서 동일한 성질이 유지된다.
• 덧셈 평균화 LOD 스킴은 유사한 조건 하에서 단조성을 보장하며 동일한 사전 추정치를 확보한다. 구성 요소를 별도로 계산할 수 있어 병렬 처리 성능이 향상된다.
• 반복 행렬 S₁ 및 S₂가 M-행렬이면서 노름이 1 이하임을 확인함으로써 안정성과 단조성을 증명하였으며, 이는 이산 최대원칙의 만족을 보장한다.
• σ = 1일 경우 시간에 대해 1차 정밀도를 확보하고 비조건적 안정성을 가지며, σ < 1일 경우 최대 확산 및 속도 항의 가중합에 대한 시간 간격 제약 조건 하에서 안정성이 유지된다.
• 이 스킴은 열 방정식, 농도 이동, 나비에-스토크스 유사 운동량 방정식 등 다양한 형태의 대류-확산 방정식에 적용 가능하다.
• 이론적 분석은 L∞ 및 L1 노름에서의 사전 추정치를 통해 검증되었으며, 문제 데이터에 대한 일반적인 가정 하에서 수렴성과 강건성을 확인하였다.