[논문 리뷰] Error estimates for semi-discrete finite element approximations for a moving boundary problem capturing the penetration of diffusants into rubber

Surendra Nepal, Yosief Wondmagegne외 1명|arXiv (Cornell University)|2021. 07. 02.

Advanced Numerical Methods in Computational Mathematics인용 수 5

한 줄 요약

이 논문은 고무로의 용매 확산을 모델링하는 일차원 이동 경계 문제에 대해, 반연속 유한요소법에 대한 사전 및 사후 오차 추정을 제시한다. 문제를 고정 도메인으로 변환하고 에너지 및 적분 추정, 보간-경계 부등식을 사용하여 농도와 이동 경계 위치에 대해 각각 𝒪(1)의 수렴 속도를 확립하였으며, 실험 데이터와 일치하는 수치적 검증을 통해 검증되었다.

ABSTRACT

We consider a moving boundary problem with kinetic condition that describes the diffusion of solvent into rubber and study semi-discrete finite element approximations of the corresponding weak solutions. We report on both a priori and a posteriori error estimates for the mass concentration of the diffusants, and respectively, for the a priori unknown position of the moving boundary. Our working techniques include integral and energy-based estimates for a nonlinear parabolic problem posed in a transformed fixed domain combined with a suitable use of the interpolation-trace inequality to handle the interface terms. Numerical illustrations of our FEM approximations are within the experimental range and show good agreement with our theoretical investigation. This work is a preliminary investigation necessary before extending the current moving boundary modeling to account explicitly for the mechanics of hyperelastic rods to capture a directional swelling of the underlying elastomer.

연구 동기 및 목표

• 고무 내 용매 침투를 모델링하는 이동 경계 문제에 대한 반연속 유한요소 근사의 정확도를 분석하기 위해.
• 확산종 농도와 이동 경계 위치에 대한 사전 및 사후 오차 추정을 유도하기 위해.
• 이론적 오차 한계를 수치 실험을 통해 검증하기 위해. 실험 데이터와 일치하는 결과를 얻기 위해.
• 초기 탄성체에서의 방향성 팽창을 고려한 기계-결합 모델로의 향후 확장을 위한 기초를 마련하기 위해.

제안 방법

• 이동 경계 문제를 수치적 처리를 단순화하기 위해 좌표 변환을 사용해 고정 도메인으로 변환하기 위해.
• 변환된 도메인에서 약한 해를 정의하고, [17, 18]의 사전 분석을 통해 존재성과 균일 유계성을 확립하기 위해.
• 고정 도메인 내 비선형 포아송 방정식계에 대해 에너지 및 적분 추정을 적용하여 오차 한계를 도출하기 위해.
• 이동 경계에서 발생하는 경계 항을 오차 분석에서 제어하기 위해 보간-경계 부등식을 사용하기 위해.
• 균일한 시간 이산화를 사용하여 1차원에서 유한요소법을 구현하고, 오픈소스 파이썬 도구를 활용하기 위해.
• 미세 메쉬 기반 참값과 이산 ℓ² 노름을 사용하여 수렴 속도를 계산하고 이론적 추정치를 검증하기 위해.

실험 결과

연구 질문

• RQ1 고무 내 이동 경계 확산 모델에서 확산종 농도를 근사하기 위한 반연속 유한요소법의 수렴 속도는 무엇인가?
• RQ2 이동 경계의 위치와 속도의 수치 근사는 얼마나 정확한가? 그리고 엄밀하게 도출된 오차 한계는 무엇인가?
• RQ3 수치 결과는 조밀하고 거품 고무에 대한 용매 침투 실험 데이터와 어느 정도 일치하는가?
• RQ4 실제로 사전 및 사후 오차 추정은 어떻게 비교되며, FEM 스킴의 신뢰성에 대해 무엇을 드러내는가?
• RQ5 알 수 없는 이동 경계로 인한 비선형 결합은 최적 수렴 속도에 어떤 영향을 미치는가?

주요 결과

• 확산종 농도와 이동 경계 위치에 대한 사전 오차 추정은 비선형 결합으로 인해 최적성에서 벗어나지만, 수렴 속도가 𝒪(1)임을 보여준다.
• 사후 오차 추정식 (59)는 이론적 기대와 일치하는 계산 가능한 오차 한계를 제공하며, 적응형 메쉬 개선 전략을 뒷받침한다.
• 수치 실험은 조밀하고 거품 고무의 경우 모두 시뮬레이션된 농도 프로파일과 실험 데이터 사이에 양호한 일치를 보였다.
• 메쉬 정밀도를 높이는 동안 이산 오차의 로그-로그 플롯을 통해 계산된 수렴 차수는 이론적 비율과 일치하여 오차 분석의 타당성을 확인한다.
• 유한요소 스킴은 실험 범위 내에서 안정적이고 정확한 결과를 도출하였으며, 이론적 예측에 매우 가까운 수렴 속도를 보였다.
• 미세 메쉬 기반 참값을 사용함으로써 오차 노름과 수렴 차수를 신뢰성 있게 계산할 수 있었고, 수치 구현의 타당성이 검증되었다.