[논문 리뷰] On a finite element approximation of the Stokes problem under leak or slip boundary conditions of friction type

Takahito Kashiwabara|arXiv (Cornell University)|2010. 12. 22.

Contact Mechanics and Variational Inequalities참고 문헌 18인용 수 5

한 줄 요약

이 논문은 마찰 유형의 슬립 또는 누출 경계 조건을 갖는 스토크스 방정식에 대한 유한요소 방법을 제시하며, 존재성과 유일성을 보장하기 위해 변분부등식 공식을 사용한다. 정규성 조건 하에서 수렴 속도가 $O(h^2)$까지 확립되며, 수치적 해를 위한 우자와 유형의 반복 알고리즘을 제안하고, 최적 수렴을 보여주는 수치 예제로 검증한다.

ABSTRACT

A finite element approximation of the Stokes equations under a certain nonlinear boundary condition, namely, the slip or leak boundary condition of friction type, is considered. We propose an approximate problem formulated by a variational inequality, prove an existence and uniqueness result, present an error estimate, and discuss a numerical realization using an iterative Uzawa-type method. Several numerical examples are provided to support our theoretical results.

연구 동기 및 목표

• 비선형 마찰 유형 경계 조건 하에서 스토크스 문제에 대한 종합적인 유한요소 이론을 개발하기 위해.
• 슬립 또는 누출 경계 조건이 마찰 유형인 스토크스 방정식의 유한요소 근사에 대해 존재성, 유일성 및 오차 추정을 확립하기 위해.
• 이산 변분부등식 문제의 수치적 해를 위한 우자 유형의 반복 알고리즘을 설계하고 분석하기 위해.
• 다양한 매개변수 설정 하에서 이론적 수렴 속도를 수치 실험을 통해 검증하기 위해.
• 누출 경계 조건에서 압력의 비유일성 문제와 그 오차 추정에 미치는 영향을 다루기 위해.

제안 방법

• 슬립 및 누출 경계 조건이 마찰 유형인 연속 문제를 변분부등식을 통해 공식화하기 위해.
• 이차 다각형 영역에서 P2/P1 유한요소 쌍(Taylor-Hood 요소)을 사용하여 문제를 이산화하기 위해.
• 변분부등식 이론을 통해 SBCF 및 LBCF 모두에 대해 이산 해의 존재성과 유일성을 증명하기 위해.
• 오차 추정을 유도: $H^{1+\epsilon}$ 정규성 조건 하에서 $O(h^{\text{min}{\epsilon,1/4\u007d})$, 추가적인 접선 속도 부호 조건이 있으면 $O(h^{\text{min}\u007b\epsilon,1\u007d})$로 향상된다.
• 이산 변분부등식을 해결하기 위한 우자 유형의 반복 알고리즘을 제안하고, 반복 수렴이 해로 수렴함을 증명하기 위해.
• LBCF에서 압력의 비유일성을 처리하기 위해 기준 점을 통해 압력 상수를 고정하여 오차 비교를 가능하게 하기 위해.

실험 결과

연구 질문

• RQ1 마찰 유형의 슬립 또는 누출 경계 조건 하에서 스토크스 방정식에 대해 엄밀하게 유한요소 방법을 공식화하고 분석할 수 있는가?
• RQ2 이러한 비선형 경계 조건 하에서 유한요소 근사의 수렴 속도는 무엇인가?
• RQ3 이산 문제를 효율적이고 신뢰성 있게 반복 알고리즘을 통해 어떻게 해결할 수 있는가?
• RQ4 누출 경계 조건의 경우 압력 필드에 덧셈 상수가 존재할 때 오차 분석에 어떤 영향을 미치는가?
• RQ5 적절한 정규성 조건과 접선 속도의 부호 조건 하에서 최적 수렴 속도를 달성할 수 있는가?

주요 결과

• 적절한 유한요소 이산화 및 변분부등식 공식화 하에서 슬립 및 누출 경계 조건 모두에 대해 이산 문제는 유일한 해를 갖는다.
• 슬립 경계 조건의 경우 오차 추정은 $O(h^{\text{min}\u007b\epsilon,1\u007d})$이며, 접선 속도에 유리한 부호 조건이 있으면 $O(h^{\epsilon})$로 향상되며, $\epsilon=2$일 때 최적의 $O(h^2)$ 수렴을 달성한다.
• 누출 경계 조건의 경우 압력 상수의 불확실성으로 인해 오차 추정은 $O(h^{\text{min}\u007b\epsilon/2,1/4\u007d})$이지만, 압력 상수를 적절히 처리하면 $O(h^2)$로 향상된다.
• 수치 실험은 이론적 수렴 속도를 확인하며, $N=120$일 때 기준 해와 비교해 SBCF 및 LBCF 모두에서 $O(h^2)$ 수렴을 보여준다.
• 우자 유형의 반복 알고리즘은 이산 해로 수렴하며, 반복 횟수를 통해 문제 매개변수와 초기 추측 값에 따라 수렴 속도가 달라짐을 수치 실험에서 확인하였다.
• 경계에서의 라그랑주 승수 $\lambda_h$는 $u_h(M)=0$이면 $|\lambda_h(M)| \leq 1$을 만족하고, 그 외에는 $\lambda_h(M) = \pm 1$이 되며, 이는 이론적 한계와 일치한다.