[논문 리뷰] Recent Advances of Isogeometric Analysis in Computational Electromagnetics

Zeger Bontinck, Jacopo Corno외 9명|arXiv (Cornell University)|2017. 09. 18.

Advanced Numerical Analysis Techniques인용 수 9

한 줄 요약

이 논문은 이소지오메트릭 해석(IGA)을 전자기계계산에서 고정밀도, 기하학적 정확도를 갖춘 방법으로 제시하며, B-spline 및 NURBS 기저 함수를 사용하여 CAD 기하학과 유한요소 이산화를 통합한다. 형태 최적화 및 불확실성 정량화 워크플로우에서 뛰어난 수렴성과 강건성을 입증하였으며, 주요 결과로 2D 및 3D 기준시험에서 평균 자기장 기울기 17.5% 향상 및 비균일성 요소 74.4% 감소를 확인하였다.

ABSTRACT

In this communication the advantages and drawbacks of the isogeometric analysis (IGA) are reviewed in the context of electromagnetic simulations. IGA extends the set of polynomial basis functions, commonly employed by the classical Finite Element Method (FEM). While identical to FEM with N\\'ed\\'elec's basis functions in the lowest order case, it is based on B-spline and Non-Uniform Rational B-spline basis functions. The main benefit of this is the exact representation of the geometry in the language of computer aided design (CAD) tools. This simplifies the meshing as the computational mesh is implicitly created by the engineer using the CAD tool. The curl- and div-conforming spline function spaces are recapitulated and the available software is discussed. Finally, several non-academic benchmark examples in two and three dimensions are shown which are used in optimization and uncertainty quantification workflows.

연구 동기 및 목표

• 특히 산업 적용에서 분석 시간의 80%가 메쉬 생성에 소요되는 전자기계계산에서 높은 계산 비용과 메쉬 생성의 복잡성을 해결한다.
• CAD 네이티브 기저 함수를 활용하여 고전적 FEM이 복잡한 기하학과 고차 이산화를 다루는 데 한계를 극복한다.
• 제어점을 통한 기하학과 계산 메쉬의 동시에 변형을 가능하게 하여 형태 최적화 및 불확실성 정량화 워크플로우의 원활한 통합을 가능하게 한다.
• 문제 차원을 낮추고 체적 메쉬 생성이 필요 없도록 하는 경계요소법과의 조합을 통해 이소지오메트릭 해석의 실현 가능성과 이점을 입증한다.
• 맥스웰 방정식에 적용된 컬 및 디버전스-준수 스퍼링 공간의 포괄적 리뷰를 제공하고, 스펙트럼적으로 정확한 이산화를 보장한다.

제안 방법

• CAD 도구에 내장된 B-spline 및 NURBS 기저 함수를 사용하여 곡선과 표면, 특히 원추곡선까지 정확한 매개변수화를 가능하게 하여 기하학과 장 변수를 동시에 표현한다.
• 맥스웰 방정식의 미분형식과 이산적으로 호환되는 컬-준수 및 디버전스-준수 스퍼링 함수 공간을 도입하여 중요한 물리적 성질을 유지한다.
• 이소지오메트릭 프레임워크 내에서 갈레르킨 변분 공식을 적용하여 기존 FEM 코드베이스의 직접 재사용이 가능하며, 기저 함수 루틴만 수정하면 된다.
• 경계에만 기반한 이산화를 사용하는 이소지오메트릭 경계요소법(IGABEM)을 적용하여 문제의 차원을 낮추고 삼변량 체적 매개변수화의 필요성을 제거한다.
• 글로벌 C1-연속성 스퍼링 공간을 통해 고차수이고 부드러운 해 장을 구성함으로써, 자유도 당 수렴 속도와 정밀도를 향상시킨다.
• 경계요소 공식에서 발생하는 밀집 선형 시스템을 효율적으로 압축하여 표준 하드웨어에서도 효율적인 해법을 가능하게 한다.

실험 결과

연구 질문

• RQ1 이소지오메트릭 해석은 CAD 기하학을 직접 사용함으로써 전자기계계산에서 메쉬 생성 노력의 현저한 감소를 이룰 수 있는가?
• RQ2 NURBS 기반 기저 함수의 사용이 고전적 다항식 FEM 대비 전자기 문제에서 정밀도와 수렴성 향상에 어떻게 기여하는가?
• RQ3 IGA는 기하학과 메쉬의 동시에 변형을 가능하게 함으로써 형태 최적화 및 불확실성 정량화 워크플로우를 어느 정도 지원할 수 있는가?
• RQ4 전자기 문제에서 이소지오메트릭 해석과 경계요소법을 조합할 때의 이점과 과제는 무엇인가?
• RQ5 이소지오메트릭 접근법은 맥스웰 방정식의 미분형식 구조와 이산적으로 호환되도록 어떻게 보장하는가?

주요 결과

• B-spline 기저 함수의 높은 전역 연속성 덕분에 고전적 FEM 대비 자유도 당 더 빠른 수렴을 달성한다.
• 원추곡선 포함한 복잡한 기하학을 기하학적 근사 없이도 정확히 표현할 수 있다.
• 3D 형태 최적화 기준시험에서 IGA를 사용한 후 평균 자기장 기울기 12.2% 향상 및 비균일성 요소 60.0% 감소를 기록하였다.
• 2D 기준시험에서 평균 자기장 기울기 17.5% 향상 및 비균일성 요소 74.4% 감소를 기록하여 최적화 워크플로우에서 뛰어난 성능을 입증하였다.
• 이소지오메트릭 경계요소법(IGABEM)은 작고 밀집된 선형 시스템을 생성하여 표준 데스크톱 컴퓨터에서도 해를 구할 수 있으며, 확장성을 확보하기 위한 효율적 압축 기법도 제공된다.
• 제어점을 통한 동시에 기하학과 메쉬의 변형이 가능하여 메쉬 재생이 필요 없이 견고한 형태 감도 분석 및 최적화를 가능하게 한다.