[논문 리뷰] High order gradient, curl and divergence conforming spaces, with an application to NURBS-based IsoGeometric Analysis

René R. Hiemstra, R.H.M. Huijsmans외 1명|arXiv (Cornell University)|2012. 09. 09.

Advanced Numerical Analysis Techniques참고 문헌 22인용 수 6

한 줄 요약

이 논문은 NURBS 기반의 고차 미분, 커플, 발산을 고려한 유한요소 공간을 제안하여 이산 수준에서 보존 법칙(가우스의 정리 및 스토크스 정리 등)을 정확히 만족시킨다. 기저 함수를 단위 분할에서 구성함으로써, 곡선 메esh에서도 스토크스 유동에서 정확히 발산이 없는 속도장을 보장하며, 다항수의 순서에 관계없이 최적 수렴성과 뛰어난 안정성을 확보한다.

ABSTRACT

Conservation laws, in for example, electromagnetism, solid and fluid mechanics, allow an exact discrete representation in terms of line, surface and volume integrals. We develop high order interpolants, from any basis that is a partition of unity, that satisfy these integral relations exactly, at cell level. The resulting gradient, curl and divergence conforming spaces have the property that the conservation laws become completely independent of the basis functions. This means that the conservation laws are exactly satisfied even on curved meshes. As an example, we develop high order gradient, curl and divergence conforming spaces from NURBS - non uniform rational B-splines - and thereby generalize the compatible spaces of B-splines developed by Buffa et al.[1]. We give several examples of 2D Stokes flow calculations which result, amongst others, in a point wise divergence free velocity field.

연구 동기 및 목표

• 곡선 메쉬에서 기초적 적분 법칙(가우스 정리, 스토크스 정리, 미적분학의 기본정리)을 정확히 만족하는 고차 유한요소 공간을 개발하는 것.
• 이소지오메트릭 분석(IsoGeometric Analysis, IGA)에 활용하기 위해 호환성 있는 B-스플라인 공간을 NURBS로 일반화하는 것.
• 기저 함수에 관계없이 이산 보존 성질을 보장함으로써 안정성과 물리적 일致성을 향상시키는 것.
• 압축성 없는 스토크스 및 레이놀즈 평형 문제를 정확히 발산이 없는 속도장으로 해결하는 데에 이 방법의 효과성을 입증하는 것.
• 테일러-쿠에트 및 리드드라이브 캐비티 유동과 같은 기준 문제를 통해 접근법을 검증하여 최적 수렴성과 뛰어난 안정성을 보여주는 것.

제안 방법

• 모든 단위 분할 기반 기저 함수로부터 고차 보간을 구성함으로써 기본 적분 정리(Gauss, Stokes, 미적분학의 기본정리)를 정확히 이행한다.
• 기하학적 요소(노드, 모서리, 면, 체적 요소)와 관련된 자유도를 정의하여 이산 설정에서 보존 법칙의 물리적 구조를 유지한다.
• NURBS에 이 프레임워크를 적용하여, 정확한 기하 구조 표현을 위한 비균일 유리 B-스플라인으로 호환성 있는 B-스플라인 공간을 확장한다.
• 속도, 압력, 유량 변수를 위한 혼합 유한요소 공간을 사용하는 혼합 갈레르킨 형식을 적용하여, 기초적으로 인프-서프 스테이블리티를 확보한다.
• 비균일 가중치를 가진 NURBS 패치를 사용하여 곡선 도메인을 정확히 모델링하며, 고리형 영역과 같은 복잡한 기하 구조도 포함한다.
• 다양한 유동 구성에서 속도와 압력의 최적 수렴률을 검증하기 위해 h-선택 연구를 수행한다.

실험 결과

연구 질문

• RQ1 임의의 단위 분할 기반 기저 함수로부터 유한요소 공간을 구성하여 셀 수준에서 정확히 적분 보존 법칙을 만족시킬 수 있는가?
• RQ2 NURBS 기반 이소지오메트릭 분석은 어떻게 압축성 없는 흐름에서 정확히 발산이 없는 속도장을 지원할 수 있는가?
• RQ3 비균일 곡선 메쉬에서 제안된 방법이 최적 수렴률과 이산 보존 성질을 유지하는가?
• RQ4 특수 처리 없이도 흐름 문제의 특이점(예: 리드드라이브 캐비티)을 처리할 수 있으며 물리적 일치성을 유지하는가?
• RQ5 이 방법의 보존 성질이 복잡한 유체역학 시뮬레이션에서 수치적 안정성과 정확도를 어느 정도 향상시키는가?

주요 결과

• 2D 스토크스 유동에서 정확히 발산이 없는 속도장이 달성되었으며, 전체 도메인에서 점별로 발산이 정확히 0이었다.
• 테일러-쿠에트 흐름 시험에서 다항식 차수 p = 1에서 p = 4까지 속도의 L2 노름에서 최적의 h-수렴률 p+1이 관측되었다.
• 모든 메쉬 세분화 수준과 다항식 차수에서 압력이 일정하며 기계 정밀도 수준에서 정확하였다.
• 낮은 차수의 근사와 비교적 굵은 9×9 메쉬를 사용한 리드드라이브 캐비티 흐름 시뮬레이션에서도 빌트인 보존 성질 덕분에 기준 결과와 뛰어난 일치를 보였다.
• 비균일 가중치(예: 모서리에서 0.5√2)를 가진 NURBS 패치에서도 최적 수렴률과 정확한 기하 구조 표현이 유지되었다.
• 물리적 구조를 이산 공간에 통합함으로써, 기저 함수 선택과 무관하게 자동으로 인프-서프 안정성이 보장된다.