[논문 리뷰] On Mixed Isogeometric Analysis of Poroelasticity

Yared W. Bekele, Eivind Fonn외 3명|arXiv (Cornell University)|2017. 06. 05.

Advanced Numerical Analysis Techniques인용 수 5

한 줄 요약

이 논문은 소성계수의 다중 다항식 차수를 사용하여 이격된 이소지오메트릭 형식을 제안하여, 작은 시간 단계에서 발생하는 허위 압력 진동을 억제한다. 이 방법은 동일한 다항식 차수를 사용하는 유한요소법에 비해 진동을 크게 감소시키며, 특히 고차 다항식 기반의 형상 함수와 물질 경계에서의 C0 연속성 조건에서 안정성과 수렴성이 향상됨을 보여준다. 이는 Terzaghi의 일차원 문제와 저투수성 토양 다층 모델 모두에서 입증되었다.

ABSTRACT

Pressure oscillations at small time steps have been known to be an issue in poroelasticity simulations. A review of proposed approaches to overcome this problem is presented. Critical time steps are specified to alleviate this in finite element analyses. We present a mixed isogeometric formulation here with a view to assessing the results at very small time steps. Numerical studies are performed on Terzaghi's problem and consolidation of a layered porous medium with a very low permeability layer for varying polynomial degrees, continuities across knot spans and spatial discretizations. Comparisons are made with equal order simulations.

연구 동기 및 목표

• 소성계수의 유한요소법에서 작은 시간 단계에서 발생하는 압력 진동 문제를 해결하기 위해.
• 이sovolumetric 형상 함수의 다항식 차수를 다르게 적용한 혼합 이소지오메트릭 형식의 성능 평가를 위해.
• 다항식 차수, 키스팬 간 연속성, 메쉬 정밀도가 수치 안정성에 미치는 영향을 조사하기 위해.
• 기본 테스트 문제에서 혼합 다항식 차수 형식과 동일 다항식 차수 형식을 비교하기 위해.
• 저투수성 층 근처에서 메쉬의 점진적 조정이 진동 감소에 기여하는 정도를 평가하기 위해.

제안 방법

• 이sovolumetric 형상 함수를 사용한 혼합 이소지오메트릭 형식을 적용하여, 변위와 다공성 압력에 대해 서로 다른 다항식 차수를 사용한다.
• 일致된 적분을 사용한 갈레르킨 방법을 통해 비트의 응집 방정식의 약한 형태를 구현한다.
• 소용 없는 임계값 이하의 시간 단계를 사용하여 후진 오일러 방법을 적용함으로써 시간 이산화를 수행한다.
• 물질 경계에서 Cpp−1 또는 C0 연속성 조건을 적용하여 물질 성질의 불연속성을 모델링한다.
• 균일한 메쉬와 점진적 메쉬를 사용하며, 저투수성 층 근처에서의 정밀도를 높여 급격한 기울기를 포착한다.
• Terzaghi의 1차원 응집 문제와 저투수성 중간 층을 가진 2차원 다층 토양 모델에 대해 수치 실험을 수행한다.

실험 결과

연구 질문

• RQ1 혼합 다항식 차수 형식이 동일 다항식 차수 형식에 비해 작은 시간 단계에서 다공성 압력 진동을 얼마나 효과적으로 억제하는가?
• RQ2 다항식 차수를 증가시킬 경우 다공성 압력 해의 안정성과 정확도에 어떤 영향을 미치는가?
• RQ3 물질 경계에서 Cpp−1 대비 C0 연속성 조건이 다공성 응력 시뮬레이션의 진동 행동에 어떤 영향을 미치는가?
• RQ4 물질 경계 근처에서 메쉬의 점진적 조정이 저투수성 영역의 진동 감소에 얼마나 기여하는가?
• RQ5 혼합 이소지오메트릭 방법이 복잡한 다층 다공성 탄성 매체에서 진동을 효과적으로 억제할 수 있는가?

주요 결과

• 혼합 다항식 차수 이소지오메트릭 시뮬레이션은 동일 다항식 차수 형식에 비해 다공성 압력 진동을 크게 감소시켰으며, 특히 작은 시간 단계에서 두드러진다.
• 압력의 다항식 차수를 증가시킬수록 진동이 감소하였고, 다공성 압력 오차의 상대적 L2 노름은 고차 근사에서 수렴하였다.
• 물질 경계에서 C0 연속성 조건이 Cpp−1 연속성 조건보다 더 우수한 결과를 보였으며, 특히 저투수성 층에서 두드러졌다.
• 저투수성 층 근처에서 요소를 정밀하게 조정한 점진적 메쉬는 동일 다항식 차수와 혼합 다항식 차수 모두에서 진동을 극적으로 감소시켰다.
• Haga 문제의 경우, pp = 4 이며 C0 연속성 조건을 적용한 혼합 다항식 시뮬레이션은 거의 진동 없이 안정적인 결과를 보였다.
• 결과는 물질 불연속성 기반의 적응형 메쉬 정밀도 전략이 다공성 응력 시뮬레이션의 수치 성능을 크게 향상시킬 수 있음을 시사한다.